Infant Survival Among Free-Living Bonnet Macaques (Macaca radiata) in South India

Female reproductive success depends to a large extent on infants’ ability to survive to maturity. While most studies of female reproductive success have focused on the effects of individuals’ sociodemographic factors (e.g., age/parity, dominance rank) on offspring survival among wild primates living in less disturbed habitats, little research has focused on offspring survival in urban or periurban animals. Here we investigated sociodemographic and anthropogenic determinants of infant survival (up to 1 yr of age) in free-ranging bonnet macaques (Macaca radiata) living in a periurban environment in Southern India. We conducted the study from November 2016 to May 2018, on two groups of bonnet macaques at the Thenmala tourist site in the state of Kerala. Fifty infants were born across two birth seasons. Of these infants, 29.2% died or disappeared in 2017 and 26.9% died or disappeared in 2018. We found that infant survival was strongly influenced by the mother’s parity: infants of experienced mothers had a better chance of survival than those of first-time mothers. We also found that male infants were more likely to die than female infants. However, we found no effects of mothers’ dominance rank, or of frequency of mothers’ interactions with humans and time spent foraging on anthropogenic food, on infant survival. Our results, consistent with findings from other wild primate species, show that even in challenging human-impacted environments, experienced bonnet macaque mothers have greater success than inexperienced ones.


Introduction
Infant survival is a major contributor to animal fitness, which makes the causal factors that affect infant survival in wild animals such as primates of interest to both evolutionary and conservation biologists. Broadly, the likelihood of primates surviving beyond infancy is determined by differential maternal investment and the quality of maternal care (van Noordwijk 2012). Maternal behavior may, in turn, depend on several factors, primarily access to food resources (e.g., Lee 1987), but also the mother's experience in raising offspring, her age, and her social status (Abello and Colell 2006;Harcourt 1987;Pusey et al. 1997;Robbins et al. 2007;Wasser et al. 2004). Thus, infant survival may depend on many factors that are attributable to individual animals (e.g., mothers' age, parity, sex of the infant), their social environment (e.g., mothers' dominance status), and their ecological environment (e.g., access to food resources, predation, human impact on the natural environment).
Female parity is often correlated with age, as first-time (primiparous) females are usually younger animals than those with a good deal of maternal experience. As primiparous females are often still growing, they must allocate their energy to both their offspring's growth and their own growth (Bercovitch et al. 1998;Setchell et al. 2002). Consequently, younger mothers may not be able to invest sufficiently in their offspring to keep them alive. Moreover, as the reproduction rate of older multiparous females is lower, they are more likely than younger mothers to invest care and protection in their offspring (e.g., Arlet et al. 2014Arlet et al. , 2019, which may in turn lead to higher infant survival. Indeed, a literature review found that younger mothers' infants were more likely to die than those of older mothers in 15 of 21 primate species (Pusey 2012). In female rhesus macaques (Macaca mulatta), age-related change in milk quality (or yield) may be responsible for differences in infant size (Hinde and Capitanio 2010). As females become older and stop growing, the energy available in their milk increases (Hinde and Capitanio 2010). Consistent with this finding is that first-born yellow baboon (Papio cynocephalus) infants have lower body mass than later-born infants, which may suggest low maternal investment by mothers that are still growing (Altmann and Alberts 2005).
In many primates, male infants require more maternal investment than female infants (e.g., chimpanzees, Pan troglodytes: Nishida et al. 2003; mountain gorillas, Gorilla beringei: Robbins et al. 2007; Hanuman langurs, Semnopithecus entellus: Ostner et al. 2005; rhesus macaques: Bercovitch et al. 2000;Dettmer et al. 2016). Greater maternal investment is suggested by longer interbirth intervals following the birth of sons (Bercovitch and Berard 1993;Bercovitch et al. 2000;Ostner et al. 2005;Robbins et al. 2007), longer lactational periods with sons (Ostner et al. 2005), greater body mass in sons (Bercovitch et al. 2000;Hinde 2009), and higher quality milk for sons (Hinde 2007), compared to daughters. Therefore, greater male mortality in infancy has been related to male infants needing more maternal investment than female infants, but not receiving enough maternal investment to thrive (Drevenstedt et al. 2008;Hinde 2007). In rhesus macaques, for instance, the overall energy available in young mothers' milk is lower than in multiparous females' milk, and therefore primiparous females are often unable to sustain that investment in sons (Hinde 2007).
Prior experience of infant care can also have a positive impact on a mother's own infant survival (chimpanzee, Pan troglodytes: Pusey 2012; vervets, Chlorocebus pygerythrus : Fairbanks 1990; lion tamarins, Leontopithecus spp. : French et al. 1996; Japanese macaques, Macaca fuscata: Schino and Troisi 2005; great apes: Abello and Colell 2006). Females that survive long enough to have multiple offspring gain experience with care of their own infants, which may improve infant survival and thus reproductive success (e.g., Arlet et al. 2019). Another way to gain experience with infant care is alloparenting, which may enhance individual animals' own parenting skills (Silk 1999). Young females usually handle infants more clumsily and become more skilled with age and experience (Hrdy 1976;Meaney et al. 1990;Silk 1990). Moreover, help from others enhances the likelihood that a young female will raise her first offspring successfully (Fairbanks 1990). However, rough handling by group conspecifics Isbell et al. 2009), including subadult and adult females with no infants, may lead to infant death. This may explain why, among bonnet macaques, for example, mothers are reluctant to allow their offspring to be handled, although females of all ages are strongly attracted to infants (Silk 1990;Simonds 1974).
Female dominance rank may impact a female's energetic balance and behavior, and thereby influence the survival of her offspring. Higher-ranking females may be able to provide greater maternal investment than lower-ranking females, achieving higher infant survival (van Noordwijk andvan Schaik 1987, 1999). As mothers vary in their quality, they invest differential amounts of energy in infants (Garcia et al. 2009). Infants of high-ranking female olive baboons (Papio anubis) and chimpanzees at Gombe (Tanzania) tend to have higher survival than infants of low-ranking females (Packer et al. 1995;Pusey et al. 1997). Dominance rank also influences aggression directed toward the mother and her infant. In Japanese macaques, for instance, lowranking mothers receive twice as much aggression as high-ranking mothers and are therefore more prone to abandoning their infants (Schino and Troisi 2005) and that, in consequence, leads to the infants' death.
Aside from demographic attributes and social factors, ecological factors such as females' access to resources and exposure to predation may affect infant survival (Sterck et al. 1997). While most studies in this regard have focused on wild, relatively unperturbed primate populations, we know little about how urban ecology may impact infant survival among primates living in human-impacted or anthropogenic environments. Anthropogenic habitats can impose challenging, stressful living conditions on wild primates (Fuentes 2012;McKinney 2015), causing individuals to be more vigilant in order to avoid antagonistic, stressful interactions with humans (Balasubramaniam et al. 2020a;Kaburu et al. 2019a, b;Marechal et al. 2011;Marty et al. 2019a), which may negatively impact reproductive success. Furthermore, primates in anthropogenic habitats, particularly in urban or periurban environments, may change their foraging strategies to rely more on human-provisioned food, which is generally more predictable and of higher calorific value than natural food McLennan and Ganzhorn 2017;Oro et al. 2013;Riley et al. 2013 ). Access to human-provisioned or anthropogenic food may also have consequences for reproductive success and has been previously linked to shorter interbirth intervals and lower infant mortality in primates (Jaman and Huffman 2013;Kurita et al. 2008;Warren et al. 2011). However, anthropogenic food has also been associated with health-related costs (e.g., Japanese macaques, Macaca fuscata: Sugiyama et al. 2014; olive baboons: Sapolsky and Share 2004), and therefore may negatively (rather than positively) impact infant survival over longer periods.
Bonnet macaques (Macaca radiata) are well suited for studying the effects of maternal demographic and social attributes on infant survival, in the relative absence of alloparenting (Rahaman and Parthasarathy 1969;Silk 1990;Simonds 1974). They predominantly live in multimale, multifemale groups, in which females are philopatric and males tend to disperse when they reach sexual maturity (Singh et al. 1984). Females form stable and matrilineally heritable dominance ranks (Silk et al. 1981a). They generally live in seasonal environments (dry and wet season) in Southern India, and most females within a group resume cycling synchronously and give birth during the dry season (typically January-April), just before the monsoon (Rahaman and Parthasarathy 1969;Singha 2001). After the infant is born, the mother-infant relationship continues for ca. 6-10mo and is strongest during the first few weeks of life (Nathan and Kaufman 1972;Rahaman and Parthasarathy 1969;Simonds 1965;Singh et al. 1980). Bonnet macaque mothers are often protective of their infants during this period, and allomothering is rare (Rahaman and Parthasarathy 1969;Simonds 1974).
Bonnet macaques are also a suitable model species to assess the impact of urban ecology, specifically anthropogenic factors, on infant survival. Many populations of this species inhabit roadsides, temples, villages, agricultural fields, and small towns of Southern India, where they live close to humans; a relatively smaller portion of these macaque populations inhabit forests and protected areas (Kumara et al. 2010;Sinha 2013;Sinha and Mukhopadhyay 2013). Over the years, the habitat of bonnet macaques in India has reduced substantially due to agricultural expansion, urbanization, and road development (Sinha 2013;Sinha and Mukhopadhyay 2013). At the same time, the probability of mortality due to factors such as collision with vehicles and electrocution has increased (Erinjery et al. 2017). Such anthropogenic factors have been implicated in a drastic population decline of bonnet macaques in recent years (Erinjery et al. 2017;Kumar et al. 2011;Singh and Rao 2004;Singh et al. 2011).
Here we assessed the social, demographic, and anthropogenic determinants of infant survival in wild bonnet macaques living in periurban environments in Southern India. We predicted that if prior experience of maternal care positively impacts infant survival among bonnet macaques, then infant survival will be greater among older and/or multiparous mothers compared to among younger, primiparous mothers. If high dominance rank is associated with increased reproductive success, then dominance rank will show a positive effect on infant survival. If mothers tend to be more protective of their daughters than their sons, which are also likely to be more mobile or exploratory in their behavior, then female infants would be more likely to survive compared to male infants. Finally, if exposure to anthropogenic environments is stressful and imposes costs on primates, then we predicted that infant survival will be negatively associated with the mothers' frequencies of interactions with humans. In contrast, if access to anthropogenic foods that are generally higher in calorific value than natural foods offset the costs of interacting with humans, then females' tendencies to forage on anthropogenic foods will positively (rather than negatively) impact infant survival.

Study Area and Subjects
We conducted the study between November 2016 and May 2018 at Thenmala dam and ecotourism recreational area, in Kerala, Southern India (8.9°N, 77°E). Aside from a few small forest patches, the area covered three villages, forest offices, and an ecotourism center. The area is home to several groups of bonnet macaques. We focused on two groups that were different in size, whose ranges partially overlapped with each other and covered forest patches and anthropogenic landscapes. During the study period, the home range of the large group (LG; Table I) was ca. 0.32 km 2 . The home range of the small group (SG) was ca. 0.16 km 2 (Balasubramaniam et al. 2020a).
The subjects of this study were all females with infants: in 2017, we observed 26 mother-infant pairs, and in 2018 we observed 24 mother-infant pairs. To determine age classes, we adapted the approach originally proposed by Strum and Western (1982) for olive baboons. Specifically, we categorized females into the following three categories based on reproductive status: young female (5-8 yr), middle-aged female (9-13 yr), older female (14-18 yr), and old female (>19 yr). We estimated these age classes based on relative size, nipple length, and general marks of aging, such as wrinkled skin on the face and the neck, saggy skin around the face, and thinning hair.
When we noticed a new infant in a group, we counted it as a birth. We defined infants as individuals <12 mo old. We designated an infant as "dead" when we first noticed that it was either dead or no longer present in a group (since infants cannot survive by themselves).

Data Collection
We collected demographic data from November 2016 to May 2018, by recording information on group composition and membership daily. All adult macaques were well habituated to human presence and were individually recognizable using natural markings such as face characteristics (color, scars, wrinkles), relative body size, and nipple color and size. We collected behavioral data on macaque-macaque agonistic interactions, human-macaque interactions, and macaque foraging on natural and anthropogenic food for 11 mo between July 2017 and May 2018. MA, RS, and four other field assistants collected 10-min focal animal sampling on all adult macaques (Altmann LG SG Young adult females 5-8 6 6 Middle-aged females 9-13 9 7 Older females 14-18 5 1 1974). The observers reached interobserver reliability in focal sampling prior to commencing data collection, with MA as the standard (range of Cohen's κ: 0.83-0.96). On each day, we followed a predetermined, randomized sequence to find and sample the focal animals, aiming to obtain two 10-min focal samples of each individual each week, one in the morning and one in the afternoon. During focal sessions, we recorded dyadic agonistic interactions (aggressive interactions and submissive signaling; Boccia et al. 1988;Silk et al. 1981a) continuously. We also collected human-macaque interactions continuously, i.e., behaviors initiated by a focal animal toward one or more humans (e.g., approach, contact and noncontact aggression, begging for food), or by one or more human(s) toward a focal animal (e.g., approach, aggression, provisioning with food, taunting) that elicited a response from the recipient. Finally, they also recorded the focal animals' main activity (resting/moving/socializing/foraging on natural or anthropogenic food) once every 2 min in the form of point-time scans within the focal sessions (Altmann 1974). We provide further details and definitions of protocols, behaviors, and definitions/criteria used to classify anthropogenic and natural foods in our previous publications (Balasubramaniam et al. 2020a;Kaburu et al. 2019a, b;Marty et al. 2019a, b).

Dominance Rank of Mothers
From all dyadic agonistic interactions with a clear winner and loser, we calculated the dominance ranks of individuals of each group, using the Perc package (Fujii et al. 2015) in R (R Core Team 2020). To control for group size, we standardized the ordinal ranks for each group and sex to create a rank index, hereafter "percentile rank," that ranged between 0 (lowest ranked individual) and 1 (highest ranked individual) (as in Balasubramaniam et al. 2020a).

Anthropogenic Factors: Human-Macaque Interactions and Foraging on Anthropogenic Food
From the continuously collected focal data, we calculated the overall frequency of interactions with humans relative to the total observation time for each female. From the point time samples of macaques' main activity, we calculated the use of anthropogenic food for each individual as the number of point-time samples in which they foraged on anthropogenic food divided by the total number of point-time samples they foraged on either natural or anthropogenic food (Balasubramaniam et al. 2020a;Marty et al. 2019b).

Statistical Analysis
We first ran diagnostic tests to confirm that the age and dominance rank of the females correlated with parity. To determine the relationship between macaque age class and parity, we ran a chi-squared test. To determine the effect of rank on parity, we ran a generalized linear mixed model (GLMM) with a binomial function, using the lme4 package in R (Bates et al. 2015; R Core Team 2020). The outcome variable was the parity of the female coded as "0" for primiparous, and "1" for multiparious. As predictors, we included dominance rank percentile as a main effect and group ID as a random effect.
To test our predictions, we used a Cox proportional hazards model to model survival to 1 yr. We allowed different baselines for the 2 yr. This approach allowed the year effects to vary in the model while making no attempt to statistically quantify this variation (as 2 yr is not enough to do this reliably). We used sex of the infant, anthropogenic factors (i.e., human-macaque interactions and preference for anthropogenic food), and demographic and social characteristics of the mother (choices based on the results of the above diagnostic tests) in the model. We also included group ID as a factor in the model to account for potential intergroup differences in our outcome. We used likelihood ratio tests (LRTs) to compare the full model with one without each predictor variable in turn, to assess the effects of each predictor variable while accounting for the effects of all other predictor variables.
We performed all analyses using the statistical software R (Bates et al. 2015; R Core Team 2020). We fitted the Cox proportional hazards model using survival R package (Therneau 2011) and calculated the variance inflation factors (VIFs) using the rms package (Harrell 2020).

Ethical Note
All data were collected from habituated, free-living animals without interfering with their activities. We conducted the study outside protected areas and complied with Indian laws and those of the Kerala Forest Department. The protocols were approved by the Forest Department in Thenmala at the time the field research was conducted and reviewed and approved by the University of California Davis Institutional Animal Care and Use Committee. The authors declare that they have no conflict of interest.

Data Availability
The datasets analyzed during the current study are available from the corresponding author on reasonable request.

General Patterns of Infant Mortality
During two birth seasons (2017-2018), 34 females gave birth to 50 infants (N = 24 females and N = 26 males; Table II). Ten infants were born to primiparous females and 40 infants were born to multiparous females. There were 31 births in LG group and 19 in SG group. Fourteen infants (41%) died or disappeared during the study (Table II).

Diagnostic Tests
As we predicted, females' age was strongly associated with parity (χ 2 (2) = 29.2, P < 0.01). This informed our decision not to include parity and age in the same Cox hazards model (due to the strong overlap of the variable categories), and only use parity. The GLMM showed that dominance rank percentile also strongly predicted parity: high- ranking females were more likely to be multiparous than low-ranking females (β = 4.90, z = 2.62, P = 0.01). Nevertheless, we still included parity and rank in the same Cox hazards model; VIF values for these variables (1.9 for both rank and parity) suggested that their collinearity was not an issue. We could not calculate the VIF for age in a model including parity because parameter estimates for the model did not converge, presumably because of the strong association between age and parity.

Infant Survival Relative to Mother's Parity, Rank, Infant Sex, and Anthropogenic Factors
Our proportional hazards model (LRT: 15.9, df = 6, P = 0.01; Wald: 13.3, df = 6, P = 0.04; score: 17.1, df = 6, P = 0.01; Table III) showed that parity significantly predicted infant survival: infants of primiparous mothers had higher mortality than those of multiparous mothers (LRT: 9.01, df = 1, P = 0.01; Figs. 1 and 2). We also found that male infants had lower survival than female infants (LRT: 4.37, df = 1, P = 0.04; Figs. 1 and 2). In contrast, rank was not a significant predictor of mortality. We found no significant effect of human-macaque interactions or foraging on anthropogenic food on infant survival. Finally, group ID had no significant effect on survival.

Discussion
Consistent with findings in other primate species, we found that parity strongly influenced infant survival, and that infants of multiparous mothers were three times more likely to survive than those of primiparous mothers. This implies a difference in the quality of maternal care between primiparous and multiparous mothers, which may be related to their experience. Moreover, younger mothers are more limited in the energy they can allocate to reproduction, because they may not have stopped growing (Bercovitch et al. 1998;Clutton-Brock 1984;Descamps et al. 2008;Hinde 2007;Hinde and Capitanio 2010;Setchell et al. 2002). The relationship between parity and infant survival may be linked to different energetic needs in young and old females. Younger females may also be less able to carry their offspring everywhere because the females are smaller than older females. As alloparental interactions are rare in bonnet macaques (Silk 1990;Simonds 1974), first-time mothers may not receive much help with care for their infants. With subsequent births, however, the quality of their maternal care may increase (Maestripieri and Carroll 1998).
We also found that male infants were more likely to die than female infants. Greater male mortality in infancy has also been reported in captive rhesus macaques and may  be related to male infants needing more maternal investment than female infants but not receiving sufficient maternal investment to thrive (Hinde 2007). A similar pattern has also been documented in humans (Drevenstedt et al. 2008). As a proximate explanation, higher mortality of infant males may be a consequence of infant temperament. For instance, in squirrel monkeys (Saimiri sciureus), rhesus macaques, and patas monkeys (Erythrocebus patas), male infants are more exploratory, and engage in play behavior more than female infants (Biben 1986;Brown and Dixson 2000;Rowell and Chism 1986). Such behaviors may bear particular risks in anthropogenic environments, where exploring infants may be more exposed to roads, vehicles, electric wires, and interactions with humans and feral dogs. Our results differ from earlier work on captive bonnet macaques, which did not show such sex differences in infant mortality (Silk et al. 1981b). Here, such a result may be a consequence of living in wild, and particularly human-impacted environments, where macaques face more variable, sometimes unpredictable or risky socioecological conditions, like differential access to natural and anthropogenic food, male migration, encounters with humans and other urban animal predators, as well as exposure to traffic, sharp fences and substrates, and electric wires. We found no evidence that having a higher rank affording greater infant survival. Studies that do show a correlation between rank and infant survival suggest that this effect may be more pronounced under certain conditions, such as low or no predation  pressure and/or rank-related effects on access to food (e.g., olive baboons: Packer et al. 1995;chimpanzees: Pusey et al. 1997;vervet monkeys: Cheney and Seyfarth 1987; long-tailed macaques Macaca fascicularis: van Noordwijk and van Schaik 1987). At our study site, however, bonnet macaques coexist with feral dogs, and occasionally (in parts of their home range) with leopards or Indian rock pythons, which are natural predators of the macaques (Coss et al. 2002). They also experience potentially stressful interactions with humans and anthropogenic environments and coping with this may require greater intragroup cohesion and social tolerance (Balasubramaniam et al. 2020a;Marechal et al. 2011). In many primate species (including bonnet macaques), dominance rank offers better access to preferred food items (e.g., Barton 1993;Chancellor and Isbell 2009;Koenig 2000;Marty et al. 2019b;Whitten 1983). Although dominance rank may generally impact macaques' access to anthropogenic food , in our population rank was not correlated to access to anthropogenic food resources (McCowan, unpubl. data). This may have been because both natural and anthropogenic food was abundant both spatially and temporally throughout the study. Under such ecological conditions of moderate predation pressure, high human impact, and greater abundance of food, dominance rank does not have a strong impact on female fitness and/or infant survival.
We found that many infant deaths were side effects of living in an anthropogenic environment, and multiple infants (N = 5) died from being hit by cars or electrocuted. Anthropogenically modified habitats can present challenging, stressful environments for primates (McKinney 2015;McKinney and Dore 2018). Mothers that are relatively naïve to infant behavior might find it more difficult to combine infant care and vigilance while walking on the roads, passing bridges, or climbing poles with electrical wires than more experienced mothers. Indeed, previous studies on bonnet macaque populations in anthropogenic habitats showed that the survival of immature individuals was markedly lower than those in forest habitats (Singh and Rao 2004). In this light, a comparative study of the predictors of female reproductive success and infant survival among bonnet macaques living in human-impacted environments (like our study population) to macaques in relatively unperturbed forest habitats would be a logical next step.
The anthropogenic causes of infant deaths contrast with our finding that the anthropogenic exposure of individual females, specifically the frequencies of interactions with humans and time spent foraging on anthropogenic foods, had no effect on infant survival. One reason for this discrepancy may be that, in the short term, the costs females incurred when interacting with humans (e.g., injuries, stress, illness) may have been offset by the benefits of access to high-energy, abundantly available anthropogenic food . In other words, these two anthropogenic factorsinteractions with humans and foraging on anthropogenic food-may have negated each other's impact on female fitness and hence infant survival. A second reason could be that frequencies of direct interactions with humans, as measured and evaluated here, were uniformly and relatively (to other macaque populations) low in our study population, possibly too low to impact interindividual differences in females' fitness (Balasubramaniam et al. 2020b). Finally, it is conceivable that the impact of interactions with humans and foraging on human foods on offspring survival may be more discernible 1) over longer time frames that affect chronic health and fitness outcomes, or through 2) comparative studies of other populations facing varying degrees of anthropogenic environmental impact. The short time frame of this study, which covered just two birth seasons, was one of its limitations.
Our study has implications for the conservation of bonnet macaques. Generally, understanding socioecological factors that affect reproduction and survival is important in increasing our ability to effectively manage populations threatened with extinction, even when models are built in contexts where the study species is not threatened. The fact that naïve mothers in particular may be more vulnerable to losing offspring in unpredictable anthropogenic environments suggests that the Thenmala bonnet macaque population remains vulnerable. As a species, bonnet macaques are restricted to peninsular Southern India, where an expanding human population is associated with a decreasing population trend (Sinha 2013;Sinha and Mukhopadhyay 2013). It will therefore be important to design and implement cost-effective tools for monitoring population health; gain more insights into the mortality of bonnet macaques; and more generally better distinguish the effects of social, ecological, and anthropological stress on reproduction and infant survival.