Maternal stress effects on infant development in wild Verreaux's sifaka (Propithecus verreauxi)

Maternal effects mediated by nutrients or specific endocrine states of the mother can affect infant development. Specifically, pre- and postnatal maternal stress associated with elevated glucocorticoid (GC) output is known to influence the phenotype of the offspring, including their physical and behavioral development. These developmental processes, however, remain relatively poorly studied in wild vertebrates, including primates with their relatively slow life histories. Here, we investigated the effects of maternal stress, assessed by fecal glucocorticoid output, on infant development in wild Verreaux’s sifakas (Propithecus verreauxi), a group-living Malagasy primate. In a first step, we investigated factors predicting maternal fecal glucocorticoid metabolite (fGCM) concentrations, how they impact infants’ physical and behavioral development during the first 6 months of postnatal life as well as early survival during the first 1.5 years of postnatal life. We collected fecal samples of mothers for hormone assays and behavioral data of 12 infants from two birth cohorts, for which we also assessed growth rates. Maternal fGCM concentrations were higher during the late prenatal but lower during the postnatal period compared to the early/mid prenatal period and were higher during periods of low rainfall. Infants of mothers with higher prenatal fGCM concentrations exhibited faster growth rates and were more explorative in terms of independent foraging and play. Infants of mothers with high pre- and postnatal fGCM concentrations were carried less and spent more time in nipple contact. Time mothers spent carrying infants predicted infant survival: infants that were more carried had lower survival, suggesting that they were likely in poorer condition and had to be cared for longer. Thus, the physical and behavioral development of these young primates were impacted by variation in maternal fGCM concentrations during the first 6 months of their lives, presumably as an adaptive response to living in a highly seasonal, but unpredictable environment. The early development of infants can be impacted by variation in maternal condition. These maternal effects can be mediated by maternal stress (glucocorticoid hormones) and are known to have downstream consequences for behavior, physiology, survival, and reproductive success well into adulthood. However, the direction of the effects of maternal physiological GC output on offspring development is highly variable, even within the same species. We contribute comparative data on maternal stress effects on infant development in a Critically Endangered primate from Madagascar. We describe variation in maternal glucocorticoid output as a function of ecological and reproductive factors and show that patterns of infant growth, behavioral development, and early survival are predicted by maternal glucocorticoids. Our study demonstrates how mothers can influence offspring fitness in response to challenging environmental conditions.


Introduction
The survival and reproductive success of animals are determined by adequate responses to challenges and opportunities in their physical and social environments (Alberts 2019). These responses are modulated by phenotypic differences Communicated by M. A van Noordwijk.
1 3 among individuals, which in turn are shaped during development by additive genetic and maternal effects (Moore et al. 2019). While these developmental processes remain relatively poorly studied in wild vertebrate populations, it has been established that maternal effects are mediated by nutrients or specific endocrine states of the mother that affect offspring morphology, but also life history traits and behavior, and that they tend to influence juvenile traits more than adult traits (Bernardo 1996;Moore et al. 2019;Mouton and Duckworth 2021). In addition, it is now widely recognized that events experienced by individuals during early development and adolescence not only have instantaneous fitness effects, but that they can also have downstream consequences for behavior, physiology, survival, and reproductive success well into adulthood, in taxa ranging from insects to humans (Patin et al. 2002;Lea et al. 2015;Tung et al. 2016;Morimoto et al. 2017;Langenhof and Komdeur 2018;Sachser et al. 2018;Schülke et al. 2019;Rosenbaum et al. 2020;Snyder-Mackler et al. 2020; Thompson and Cords 2020;Weibel et al. 2020). In fact, the environment individuals experience during early development can even impact the phenotypes and fitness of their own later offspring (Zipple et al. 2019(Zipple et al. , 2021. Ultimate explanations for maternal effects focus on phenotypic adaptations to environmental factors, either instantly or in the future, but the effects themselves may also result from intrinsic factors, such as maternal age (Monaghan 2008). These effects are based on the fact that, as a consequence of anisogamy, female fitness is primarily limited by mothers' ability to provide pre-and postnatal maternal investment, which, in the end, is linked to maternal condition and access to resources. If a mother is exposed to favorable environmental conditions (e.g., in terms of climate, food, predators, or parasite exposure), and if these conditions are stable and predictable, her offspring should enjoy relatively high fitness under similarly good conditions. However, if the environmental conditions are more variable and less predictable, there is a risk of a mismatch between the maternal environment and that of her then adult young (Nettle and Bateson 2015). If offspring that developed under good conditions achieve the highest fitness regardless of future environmental conditions, they enjoy a so-called silver-spoon effect (Grafen 1988). If maternal effects adapt offspring to the expected future environmental conditions, they exhibit an external predictive adaptive response (PAR). When infants develop instead under developmental constraints that affect their somatic state negatively, and they adjust to this unfavorable internal state (e.g., by accelerating reproduction), regardless of future environmental conditions, this represents an internal PAR (Gluckman et al. 2005;Sheriff and Love 2013;Hanson and Gluckman 2014;Schülke et al. 2019). Because the predictability of later-life environments decreases with increasing longevity, external PARs are not expected in long-lived species (Wells 2007a), such as most primates (van Noordwijk 2012;van Noordwijk et al. 2013;Berghänel et al. 2016;Snyder-Mackler et al. 2020). Here, an alternative model has been proposed, postulating that earlylife effects promote immediate survival, independent of any costs this may incur later in life (Lea et al. 2015(Lea et al. , 2018, though these costs may be partly compensated for by internal PARs (Nettle et al. 2013).
In addition to ecological crises, infants may experience stressful social events, like weaning-related rejection, maternal death, or birth of a sibling (Tung et al. 2016;Maestripieri 2018), or they are impacted by maternal senescence (Ivimey-Cook and Moorad 2020), all of which tend to reduce the amount of maternal investment they receive. There is empirical evidence from studies of both humans and other vertebrates for rippling effects of early life adversity on offspring development and fitness (Lummaa and Clutton-Brock 2002;Weinstock 2008;Giesing et al. 2011;Berghänel et al. 2017;Crockford et al. 2020;McGhee et al. 2021), but the proximate mechanisms underlying this phenomenon remain poorly known, especially in wild animals (Edes and Crews 2017;Langenhof and Komdeur 2018;Lea et al. 2018;Hawkley and Capitanio 2020). However, one factor that has been identified as important mediator of the relationship between individuals' early experiences and health outcomes in several species and taxa is the exposure to maternal stress hormones, i.e., glucocorticoids (GC; Love et al. 2013), both pre-and postnatally.
In birds and other oviparous species, mothers can adjust offspring´s exposure to prenatal GC levels only during egg production (Love and Williams 2008). In mammals, however, maternal GCs are transferred through the placenta and later on through mother's milk and can therefore presumably be fine-tuned over longer periods of early offspring development (Hinde 2013). Maternal GC levels in free-living birds and mammals are predictably increased in response to unpredictable changes in food availability or quality (Love et al. 2005;Welcker et al. 2009), but also to other potential stressors, like reproduction or predation risk, that make an increase in mobilized energy beneficial (Beehner and Bergman 2017). Pre-, and in the case of mammals, also postnatal maternal stress is known to influence several aspects of offspring's phenotype, including their physical development (Berghänel et al. 2017), behavioral development (Maestripieri 2018), skill acquisition (Berghänel et al. 2015), or immune function (Veru et al. 2014). However, the direction of the effects of maternal physiological GC output on offspring development is highly variable, even within the same species (Hauser et al. 2007;Berghänel et al. 2017;Schülke et al. 2019; for intra-specific plasticity, see Dantzer et al. 2013), leading to opposite conclusions about their ultimate function. Accordingly, the observed offspring phenotypic responses to maternal GCs have either been interpreted as unavoidable negative outcomes, such as reduced offspring size or growth , or as evidence for adaptive maternal programming that prepares offspring for their predictable future environment (Sheriff et al. 2017), or as adjustments of their internal somatic state (Nettle et al. 2013). In addition, there is also evidence for sex-specific intraspecific variation in the response to maternal stress. For example, infant male and female rhesus monkeys exhibited different temperament responses to variation in maternal milk cortisol. Infant males of mothers with higher milk GC concentrations showed a more confident temperament (more active, playful, and curious), whereas female infants did not (Sullivan et al. 2011). Similarly, in humans, higher milk cortisol levels were found to induce "negative affectivity" (fear, sadness, frustration) in daughters but not in sons (Grey et al. 2013). Consequently, there is a need for studies contributing comparative information from additional taxa to document patterns and drivers of maternal stress effects on infant development.
Here, we contribute new data on the effects of maternal stress as well as environmental and social adversaries on infant development in wild Verreaux's sifakas (Propithecus verreauxi), a group-living Critically Endangered (Louis et al. 2020) primate from the seasonal dry forests of southwestern Madagascar. Specifically, we first examine the effects of variation in local climate, food availability, and previous reproduction on maternal GCs. In a second step, we explore how maternal GC output is linked to the physical and behavioral development as well as to the survival of two cohorts of infants.
Moreover, social variables were also found to induce changes in GC levels, albeit not consistently so. For example, members of larger groups exhibited elevated GC level in ring-tailed lemurs (Pride 2005;Fardi et al. 2018;Gabriel et al. 2018), but not in Verreaux's sifakas (Rudolph et al. 2019). Moreover, increased GC levels characterized topranking females in ring-tailed lemurs (Cavigelli 1999;Cavigelli et al. 2003) and dominant males in Verreaux's sifakas (Fichtel et al. 2007;Rudolph et al. 2020), but rank did not predict GC levels in male red-fronted lemurs (Ostner et al. 2008) and male ring-tailed lemurs (Gould et al. 2005). Further, fluctuations in group aggression rates associated with seasonal reproduction did not influence GC levels in male ring-tailed lemurs (Gould et al. 2005), but increases in male GC levels during the mating and/or birth season were reported for Verreaux's sifakas (Fichtel et al. 2007;Brockman et al. 2009;Rudolph et al. 2020), red-fronted lemurs (Ostner et al. 2008), and collared lemurs (Balestri et al. 2014). Finally, GC levels increased during late gestation and the birth season in female sifakas (Tecot et al. 2019;Ross 2020;Rudolph et al. 2020), collared (Balestri et al. 2014), and ring-tailed lemurs (Cavigelli 1999). Thus, lemurs apparently respond to challenging ecological and social factors with changes in GC levels, and some of these responses appear to be variable among species and between the sexes. Previous lemur studies have not linked these stress responses to infant development, however, even though this variability in GC responsiveness offers a promising opportunity for examining whether pre-and postnatal maternal GC output affects infant development and if so, in which way.
Here, we predicted that relatively harsh climatic conditions (i.e., low temperatures and rainfall) with low food availability will increase maternal average GC levels, both prenatally and postnatally. Further, focusing on developmental constraints, we predicted that elevated pre-and postnatal maternal GC concentrations will be associated with reduced maternal investment (e.g., in terms of infant carrying) and retarded infant physical and behavioral development (but see Dantzer et al. 2013). Specifically, we predicted that, as in some other species (Sullivan et al. 2011;Grey et al. 2013), infants from mothers with higher GCs will grow and develop more slowly and exhibit reduced levels of playing, foraging, locomotion, and grooming compared to mothers exhibiting lower GC concentrations. Finally, given the detrimental effects of prolonged exposure to high GC levels in adult primates , we predicted that high maternal GC levels are also associated with an increased mortality risk for infants.

Study site and subjects
Data on wild Verreaux's sifaka were collected during two successive breeding seasons (2017 and 2018) at the field station of the German Primate Center at Kirindy Forest/ CNFEREF, located in central western Madagascar (Kappeler and Fichtel 2012). The area is dominated by dry deciduous forest and a climate characterized by a hot wet season from November to March and a cool dry season from April to October. A local population of Verreaux's sifakas has been studied in Kirindy Forest since 1995. These lemurs are habituated to human observers and are being individually marked with unique neck collars when they are about 8 months old or when they immigrate into the study population.
Verreaux's sifakas are 3 kg arboreal, foli-frugivorous primates that can live up to 30 years (Richard et al. 2002). They live in groups of about 6 individuals, and there is no sexual size dimorphism (Lawler 2009;Kappeler and Fichtel 2012). Females become mature at 5 years and give birth to a single infant at the height of the austral winter (July to August) after a gestation period of about 5 months. Weaning takes place at the height of the wet season (January to February), when food is most abundant. Infants are entirely dependent on their mother during their first 3 months of life (Jolly 1966), and mothers constitute their primary caregiver and invest more in transport than the other group members (Patel 2007;Malalaharivony et al. 2021). Infant mortality is very high, with about 62% of infants dying during their first year of life (Kappeler and Fichtel 2012).

Behavioral and morphometric data
A total of 12 infants, whose exact birth dates were known from daily group censuses, and their mothers from 7 study groups served as focal animals. All of the mothers included in this study were multiparous, and four of them had an infant in both study years. Infants were born between July and August, with most being born in July. We observed infants until February for the first cohort (N = 6 infants, 2017) and until December for the second cohort (N = 6 infants, 2018). An effort was made to record the behavior of mothers and infants simultaneously for 5 daily observation hours, balanced between 3 h in the morning (8:00-11:00 h) and 2 h in the afternoon (14:00-16:00 h), using continuous focal animal sampling (Altmann 1974). Focal subjects were chosen in a randomized but counter-balanced manner throughout morning and afternoon observations. As our study involved focal animal observations, it was not possible to record data blind.
Operational definitions of all recorded behavioral elements are provided in SEM Table 1 (see also Malalaharivony et al. 2021). To quantify patterns of infant development, we specifically recorded the occurrence and duration of all infant activities, such as feeding on solid food and the time spent in nipple contact, which provides a proxy for the opportunity to obtain access to milk, solitary locomotor, and object play as well as independent locomotion. We also recorded all interactions, including social play and grooming with other group members to characterize the social development of infants. Finally, to assess maternal investment, we recorded the proportion of time infants were carried by their mothers.
We used two independent methods to assess infant physical growth rates. One non-invasive estimate was based on changes in infants' lower arm length, which was estimated once a month via photogrammetry (Breuer et al. 2007). A second estimate was based on records of body mass obtained during infants' first capture when they were about 9 months old and related to this species' estimated mass at birth (Kappeler and Fichtel 2012). Two infants died postnatally before either estimate could be obtained. Both measures were positively correlated with each other and described and discussed in detail elsewhere (Malalaharivony et al. 2021).

Phenology and climatic variables
Food availability, estimated as the relative abundance of leaves, flowers, and fruits (i.e., the main constituents of sifakas' diet, Guevara et al. 2021), was assessed during biweekly phenology transects of 784 trees in the study area as detailed elsewhere (Koch et al. 2017). Briefly, we applied a semi-quantitative method in which the availability for each food item was scored on an ordinal scale ranging from 0 to 4, where 0 reflects the complete absence of the item and 4 represents its maximum abundance. Average monthly scores for leaves, fruits, and flowers were summed up to obtain one score for food availability per month. A fully automatic weather station (Lamprecht, Göttingen, Germany) recorded data on ambient temperature and rainfall at the field station.

Collection of fecal samples and GC analyses
We attempted to collect weekly fecal samples for glucocorticoid metabolite (fGCM) analyses from mothers (N = 377) throughout gestation and lactation (SEM Table 2). Samples uncontaminated with urine were collected immediately after defecation in the morning between 07:00 and 11:00 h and transferred into a 15-ml tube containing 5 ml of 80% ethanol in water. FGCM extractions took place in our field laboratory using a validated field extraction method (Shutt et al., 2012) applied successfully in several other studies on wild primates (Rimbach et al. 2013;Hämäläinen et al. 2015;Kalbitzer et al. 2015), including Verreaux´s sifakas (Rudolph et al., 2019). In brief, samples were manually homogenized in their original ethanolic solvent, subsequently vortexed, and then centrifuged manually (Rudolph et al., 2019). About 1.5 ml of the supernatant were finally decanted into 2-ml polypropylene safe-lock tubes (Eppendorf®, Hamburg, Germany) for storage at ambient temperatures in the dark. Within 1-6 months following sample collection, fecal extracts were shipped to the German Primate Center and stored at − 20 °C until fGCM analysis. We assessed fGCM concentrations by using a validated enzyme immunoassay (EIA) for the measurement of immunoreactive 11ß-hydroxyetiocholanolone, which is a group-specific measure of 5ß-reduced cortisol metabolites that has been shown to reliably track adrenocortical activity from fecal samples in several primate and other mammal species (Braga Goncalves et al. 2016;Heistermann et al. 2006), including Verreaux's sifakas (Fichtel et al. 2007;Rudolph et al. 2019Rudolph et al. , 2020. For details of assay methodology, see Heistermann et al. (2006). Inter-assay coefficients of variation (CVs) of high-and low-value quality controls were 11.2% (high) and 9.5% (low), while both intra-assay CVs were < 10%. All hormone concentrations are expressed as ng/g wet weight of feces.

Statistical analyses
In general, we fitted statistical models (LMM, GlmmTMB) in R (version 4.0.3; R Core Team 2018). All covariates in the models were z-transformed in order to obtain models that are easier to interpret (Schielzeth 2010) and to facilitate model convergence. For LMMs, we checked the assumptions of normality distributions and homogeneity by visual inspection of a QQ-plot (Field 2009) of residuals and residuals plotted against fitted values (Queen et al. 2002). For Glm-mTMBs, we checked for overdispersion. For both types of models, we assessed model stability through the level of estimated coefficients and standard deviations (Nieuwenhuis et al. 2012). Furthermore, we checked collinearity issues by deriving variance inflation factors (VIF; Field 2009) of the standard linear model lacking the random effects. To test the significance of the predictors as a whole, we compared the fit of the full model with that of the null model comprising the intercept and the random effects (Forstmeier and Schielzeth 2011). Statistics for these comparisons are given in the summary tables of the respective models.
According to when fecal samples were collected, we categorized log-transformed fGCM concentrations into a preor postnatal period. For analyses examining the effects of maternal fGCM levels on infant development, the 6 months of the prenatal period were further divided into three periods between the estimated date of conception (calculated back from the birth date of the infant) and the date of birth because maternal GCs are known to increase during late gestation (Rudolph et al. 2019), presumably also because of increasing fetal GC production (see below). We therefore distinguished between the early-mid prenatal period, which lasted from February to May, and the late prenatal period, which included data from June and July. The postnatal period was defined as the period between the birth date and the end of data collection, spanning the months from July/ August to December for most infants of both birth cohorts.

Predictors of variation in fGCM concentrations and influence of FGCMs on maternal care
To assess variation in pre-and postnatal fGCM concentrations, we fitted a Linear Mixed Model using the package lme4 (Bates et al. 2014). Since one infant died on day 1 postnatally, we excluded postnatal samples collected from this female, resulting in a sample size of N = 359. We fitted Table 2 Factors predicting variation in time mothers spent carrying their infants (GlmmTMB; N = 317; N ID = 11; estimates, together with standard errors). Likelihood ratio test comparing the full with the null model: χ 2 = 167.35, df = 6, p < 0.001 a z-transformed, mean ± SD of the original values are: infant age, 10.14 ± 6.2 weeks; postnatal fGCm: 0.21 ± 0.1 (ng/g); early/mid prenatal fGCm, 0.24 ± 0.1 (ng/g); fruit availability, 0.11 ± 0.06; mother's age, 12.63 ± 4.63 years b Manually dummy-coded with year 2017 as reference category c Not shown as has no meaningful interpretation log-fGCM concentrations as response, period (early/mid, late prenatal and postnatal), monthly cumulative rainfall, cumulative rainfall during the preceding wet season until tentative conception (November to January), and mother's age (in years) as fixed factors. Mothers' identity was included as a random factor, and all mothers were either the only or the dominant adult female of their group. We included rainfall as a proxy for food availability in this model because fruit availability, which is associated with variation in fGCM concentrations (Rudolph et al. 2020), was colinear with increases in fGCM concentrations during early pregnancy. To avoid overconfident model estimates and to keep type I error rate at the nominal level of 0.05, we included random slopes (Schielzeth and Forstmeier 2009;Barr et al. 2013) of mother's age, rainfall before conception, and rainfall within mothers' identity. Originally, we also included parameters for the correlations between random slopes and intercepts, but as the model did not converge, we had to exclude them again.
To examine whether maternal prenatal fGCM concentrations predict postnatal fGCM concentrations, we fitted another LMM. As response we fitted postnatal fGCM concentrations, average prenatal fGCM concentrations as fixed factor and mother's identity as random factor. To control for autocorrelation, we fitted two additional models, using the package "nlme" (Pinheiro and Bates 2000), by including sample date as correlation. Samples that were collected on the same day (N = 18) were excluded in these models. Since the results did not deviate between these and the original models without controlling for autocorrelation, we present the original models and depict the results of the new models in the SEM (Table S3, S4).
To quantify inter-individual variation in maternal care, we calculated the proportion of time (minutes per hour) for each behavioral activity based on observational data collected between July and December of each study year. To estimate variation in the proportion of time mothers spent carrying infants, we fitted a Generalized Linear Mixed Model (GLMM; Baayen 2008) with a beta error distribution and logit link function (package glmmTMB; McCullagh and Nelder 1989;Bolker 2008). The proportion of time mothers spent carrying infants was fitted as response, and infants' age in weeks, infant's sex, average early/mid prenatal and monthly postnatal fGCM concentrations, study year, fruit availability, and mother's age as fixed factors. Individual and maternal identities were fitted as random factors, including monthly postnatal fGCM concentrations and mothers' age within maternal identity and fruit availability as well as monthly postnatal fGCM concentrations within individual identity as random slopes. Originally, we also included parameters for the correlations between random slopes and intercepts, but as the model did not converge, we had to exclude them again.

Influence of maternal fGCM concentrations on infant physical and behavioral development
To estimate the influence of maternal early/mid prenatal and postnatal fGCM concentrations on infants' body mass and arm length growth rates, we calculated the mean maternal pre-and postnatal fGCM concentrations and correlated them with individual growth rates derived from a generalized least square regression (Malalaharivony et al. 2021) by using Spearman rank correlations.
To estimate variation in infant behavioral development (proportion of time spent in nipple contact, foraging, and playing), we fitted three Generalized Linear Models (GLMM; Baayen et al. 2008) with a beta error distribution and logit link function (package glmmTMB; McCullagh and Nelder 1989;Bolker 2008). The proportion of time spent in nipple contact, foraging, or playing per observation hour was fitted as response, and infants´ age in weeks, average early/mid prenatal and monthly postnatal fGCM concentrations, study year, fruit availability, and mother's age were included as fixed factors. For the model investigating variation in time spent playing, we included another fixed factor differentiating between play types (solitary locomotor, solitary object, and social play). In the models fitting the time spent in nipple contact and foraging, individual and maternal identities were fitted as random factors including mother's age and fruit availability within mother's identity and postnatal fGCM concentrations and fruit availability as random slopes within individual identity without correlations between random slopes and intercepts. In the model fitting time spent playing, we had to exclude the random slope of mother's age with mother's identity because the model did not converge.
The sample sizes for the proportion of time infants spent grooming or locomoting independently were relatively small because these behaviors were more regularly shown towards the end of the study period when infants were older. We therefore calculated the mean of the time spent grooming or locomoting independently and correlated them with mean average early/mid prenatal and postnatal fGCM concentrations by using Spearman rank correlations.
To assess factors predicting infants' survival during the first 1.5 years, we calculated Cox proportional-hazards model (function "Coxph," "survival" package; Therneau et al. 2020) with the function Surv (age at 30.01.2020 | survived 1.5 years) as proportional hazard. One infant died a day after birth, which was excluded in these analyses. All other infants (n = 5) died after the end observation period at a mean age of 11.8 ± SD 4.4 months. Since some mothers gave birth in both years, we included maternal identity as a cluster to control for repeated measures. As predictors, we included average early/mid prenatal and postnatal fGCM concentrations, time being carried by the mother, body mass and arm length growth rates, time spent in nipple contact and foraging, and year of the birth cohort.
Because of the small sample size of 12 infants, we fitted a separate model for each predictor variable.

Predictors of variation in prenatal fGCM concentrations and influence of fGCM levels on maternal care
FGCM concentrations differed across periods, with higher levels during the late prenatal period in comparison to the early/mid prenatal period. (Table 1, S3, Fig. 1a). Average postnatal fGCM concentrations were lower than early/ mid concentrations, though absolute values appeared to be in between early/mid fGCM concentrations. Rainfall was negatively correlated with fGCM concentrations, with females having higher fGCM concentrations during periods with little rainfall (SEM, Fig. S1). Mothers that had given birth to an infant the year before experienced higher average fGCM concentrations (Fig. 1b). Rainfall during the wet season before conception and mother's age did not co-vary with fGCM concentrations (Table 1, S3). Since rainfall is a proxy for food availability (Rudolph et al. 2020), food shortages might have resulted in reduced maternal condition and were associated with higher stress. Prenatal fGCM concentrations predicted postnatal fGCM concentrations (LMM: N = 194, likelihood ratio test comparing the full with the null model: χ 2 = 7.61, df = 1, p = 0.006; estimate = 0.28, SE = 0.09, p = 0.016, Table S4).
The time spent carrying infants decreased with infants' age (Table 2). Both, early/mid and postnatal fGCM concentrations, predicted carrying time, with mothers having lower fGCM concentrations spending more time carrying their infants (Figs. 2a, b and 3). Fruit availability, birth cohort, and mother's age did not influence the time spent carrying infants ( Table 2).
Postnatal fGCM concentrations and fruit availability predicted the time spent in nipple contact (Table 3). Infants of mothers with higher postnatal fGCM concentrations spent more time in nipple contact. In addition, infants spent more time in nipple contact when more fruits were available. Neither infant's nor mother's age, sex, birth cohort, or early/mid prenatal fGCM concentrations covaried with time spent in nipple contact (Table 3). Infants spent more time foraging with increasing age. Early to mid fGCM concentrations predicted the proportion of time infants spent foraging; infants of mothers with higher fGCM concentrations spent more time foraging (Fig. 4, Table 4). Postnatal fGCM concentrations, birth cohort, mother's age, and fruit availability did not co-vary with time spent foraging (Table 4).
The time infants spent playing increased with age, and they spent more time with solitary locomotor than with solitary object or social play (Fig. 5, Table 5). Early/mid prenatal fGCM concentrations predicted time infants spent  playing, with infants of mothers exhibiting higher prenatal fGCM concentrations spending more time playing (Fig. 5). Infants of the birth cohort 2018 spent more time playing than infants of the birth cohort 2017, but postnatal fGCM concentrations, fruit availability, and mother's age did not predict the time infants spent playing ( Table 5)    did not correlate with either maternal pre-or postnatal fGCM concentrations (Spearman rank correlation, prenatal fGCM, N = 11, r = 0.36, p = 0.313; postnatal fGCM, N = 10, r = 0.22, p = 0.537). Infant survival was neither predicted by mid-early prenatal nor postnatal fGCM concentrations, but by the proportion of time being carried by the mother and arm length growth rate. Infants that were carried more often died earlier ( Fig. 6; Table 6). In addition, infants with higher arm length growth rates lived longer, but body mass growth rate did not predict survival. Mother's age, year of the birth cohort, and time spent foraging by the infants or being in nipple contact did not predict their survival either (Table 6).

Discussion
In this study, we investigated factors predicting maternal fGCM concentrations and their impact on infant physical and behavioral development in a wild population of Verreaux's sifakas. We found that maternal fGCM concentrations were higher in the late prenatal period, compared to the early/mid prenatal period. Average postnatal fGCM concentrations were lower than early/mid concentrations, but absolute postnatal fGCM values appear to be in between early/ mid concentrations. This dynamic might be due to the steep rise in fGCM concentrations from early to mid-gestation, which might have also been influenced by changes in rainfall and food availability at the beginning of the dry season in April/May. fGCM concentrations were higher during months with little rainfall. Since rainfall is a proxy for food availability (Koch et al. 2017;Rudolph et al. 2019), mothers might have been in poorer condition during periods of lower food availability, a state typically associated with elevated fGCM concentrations. Mothers that gave birth to an infant the previous year had higher fGCM concentrations, indicating that these mothers had to mobilize more energy or higher HPA activity and that generally higher GC levels in the mother (potentially due to her own early life experience) lead to a faster pace of life.
Most importantly, infants' physical and behavioral developments were influenced by mothers' prenatal fGCM concentrations. Specifically, infants of mothers with higher prenatal fGCM concentrations exhibited higher growth rates. Infants of mothers with higher early/mid prenatal fGCM concentrations were more active by spending more time foraging independently and being more active in all categories of play. Even though solitary locomotor and object play emerged earlier than social play (Malalaharivony et al. 2021), fGCM concentrations did not affect play types differently. In contrast, infants of mothers with low early/mid and postnatal fGCM concentrations were carried more, and those of mothers with higher postnatal GC levels spent more Kaplan-Meier survival curve. Survival probability for infants surviving during the first 1.5 years depending on the time mother's spent carrying infants (N = 6 infants surviving and N = 6 infants that died). For illustrative purposes, we divided the time spent carrying into much carrying (black line; above median time spent carrying and low (red line); below median time spent carrying) time in nipple contact. Time spent carrying infants predicted infants' survival, suggesting that infants that were potentially in poorer condition had to be carried more. Hence, infants that were carried less often and were more independent were able to successfully navigate their first ecological bottleneck, i.e., the first dry season after weaning, a period in which about 62% of infants die (Kappeler and Fichtel 2012). Thus, infants' physical and behavioral developments were presumably impacted by maternal fGCM concentrations, but it is equally possible that the demands of caring for an infant in relatively poor condition caused an increase in maternal fGCM. While our study was not designed to formally test predictions of ultimate hypotheses about maternal effects, the observed patterns of growth and behavioral development during early infancy are possibly indicative of an internal PAR response. However, "the PAR and developmental constraints models can only be distinguished by comparing fitness-related traits in individuals from high-and low-quality early environments, when each of these sets of individuals experience both high-and low-quality adult conditions" (Lea et al. 2015), which was not possible in this study. Additional future research will be required to determine how prenatally stressed offspring manage to grow faster and play more, and whether the extra time spent foraging is sufficient to provide the necessary energetic support.

Variation in maternal fGCM concentrations
Maternal fGCM concentrations of Verreaux's sifaka can be examined with respect to the costs of reproduction or seasonality, which are not independent in sifakas, and therefore challenging to disentangle (see also Charpentier et al. 2018). First, maternal GCs increased during the prenatal period and were higher than during the postnatal period, confirming results of previous studies in the same population (Rudolph et al. 2019(Rudolph et al. , 2020. Even though lactation is the energetically most challenging reproductive activity for mammals (Hinde and Milligan 2011), several other studies also reported higher maternal GC levels during gestation (particularly late gestation) than during lactation (e.g., Rimbach et al. 2013;Berghänel et al. 2016; but see Starling et al. 2010). In fact, GC levels increase across gestation in most mammals (Edwards and Boonstra 2018; but see Scott et al. 2013 for artiodactyls).
These pregnancy-related elevations in GC output are likely caused by interactions of increased metabolic demands during pregnancy and activation of the HPA axis with concomitant stimulation of corticosteroid-binding globulin due to elevated levels of estrogens as well as placental production of cortisol (Edwards and Boonstra 2018). High levels of GCs may also prepare the body for the stressful period surrounding parturition (Edwards and Boonstra 2018) and probably trigger developmental changes that prepare the late fetus for the extra-uterine environment, particularly with respect to organ maturation (Fowden et al. 1998) and may therefore not reflect direct adaptations to females' energetic costs of reproduction. Our finding of a positive link between pre-and postnatal maternal GCs, in contrast, points towards intrinsic factors related to individual energetic investments into reproduction. Our finding that mothers that had an infant the previous year had marginally higher GC levels than those that did not also indicates that some long-term costs of reproduction contributed to variation in maternal GC levels, suggesting that the physiological costs of reproduction go beyond direct and indirect costs previously characterized for mammals (Wells 2007b;Speakman 2008).
Second, maternal fGCM concentrations were higher during periods with low rainfall, which serves as a proxy for food availability in this seasonal habitat (Koch et al. 2017;Rudolph et al. 2019), but rainfall in the months before conception did not co-vary with fGCM concentrations. This may reflect the position of Verreaux's sifakas along the continuum from capital to income breeders (Lewis and Kappeler 2005). A reduction in the quantity or quality of food has been shown to increase GC levels in several other species, including lemurs (Behie andPavelka 2013, Chapman et al. 2015;Berghänel et al. 2016;Tecot et al. 2019;Rudolph et al. 2020), and may therefore reflect a response to compensate for reduced energy intake. Reduced food availability is particularly challenging for gregarious species because of feeding competition, which can be stressful for reproducing lowranking individuals (Creel et al. 2013); subordinate reproduction has been regularly observed in our study population in previous years (Kappeler and Fichtel 2012). Because these sifakas inhabit forests that are highly seasonal, but in which rainfall and food availability is highly unpredictable from year to year (Dewar and Richard 2007;Lawler et al. 2009; but see Federman et al. 2017), it is likely that physiological adaptations are fine-tuned and superimposed upon predictable seasonal cycles (Reeder and Kramer 2005;Creel et al. 2013).

Effects of maternal GC levels on maternal investment and infant development
Our study revealed that variation in maternal GC levels was associated with several aspects of maternal care and infant development. First, variation in maternal stress levels was associated with patterns of maternal care. High maternal GC levels during the prenatal phase correlated with more explorative infants and low carry effort because infants of these mothers spent more time foraging independently and more time playing. High maternal GC during the postnatal phase correlated with infants who have more nipple contact and forage more while resting and being carried less. While mothers with lower fGCM levels carried their infants more, maternal GC levels had no effect on independent infant locomotion, presumably because infants also moved when their mothers did not. It would therefore be interesting to estimate maternal carrying effort as the proportion of total maternal locomotion to assess different maternal strategies: infant carrying vs. nursing plus independent infant locomotion (Kramer 1998;Harel et al. 2021). Further, time spent carrying infants predicted infant survival in the first 1.5 years, suggesting that infants that were potentially in poorer condition had to be carried for longer, but more detailed analyses of mothering styles are required to substantiate this idea. Moreover, long-term effects of early maternal adversity on maternal effort should also be considered in this context (Patterson et al. 2021).
Second, infants exhibited faster growth rates when mothers had higher prenatal GC levels, echoing findings of previous studies reporting a positive correlation between postnatal infant growth rates and maternal GC levels (Swolin-Eide et al. 2002;Hauser et al. 2007;Dantzer et al. 2013;Hinde et al. 2015;Berghänel et al. 2016). However, other studies reported the opposite effect, prompting Berghänel et al. (2017) to propose a model that reconciled these conflicting findings by suggesting that the effects of maternal GC output on infant growth depend on the infant's specific developmental period. Because we have no data on fetal growth and our estimate of postnatal body mass growth is based on a single datum obtained after weaning, and arm length growth rates were measured only until an age of 6 months, we cannot assess the potential specific dynamics of infant growth patterns in this species. While it is likely that mothers elevate their GC levels during gestation as a result of reduced energetic intake and/or increased expenditure, this response has no detrimental effects on their current reproduction, i.e., their infants' growth rates. It has been argued that mothers, through increasing prenatal GCs, adjust their infant's life trajectory by accelerating infant's life history pace altogether (Nettle and Bateson 2015), but we lack data on the long-term consequences on other life history traits, such as sexual maturation and longevity to evaluate this possibility. Because the infants of mothers with high prenatal GC levels developed faster and began foraging independently earlier, these stressed mothers benefit from an overall reduction of maternal investment, thereby making the best of a bad start of the current reproductive cycle (Jones 2005).
Finally, we found that arm length growth and time spent carrying predicted infant survival during the first 1.5 years of life. In Verreaux's sifaka, survival during the first year of life has been found to represent the most challenging life stage, with more than half of infants dying (Richard et al. 2002;Lawler et al. 2009;Kappeler and Fichtel 2012). Thus, this relationship is important for the fitness of both mothers and infants, but it is difficult to evaluate because in our study, infant survival was not consistently linked to the two different growth estimates. Moreover, survival was not predicted by maternal GC levels or our estimates of infant's energy input, i.e., time spent foraging or in nipple contact and therefore appears to be unrelated to how much energy infants obtained. Finally, mother's age also did not predict early survival, although mothers in these two birth cohorts were rather old, with a mean age of 12.6 ± 4.6 years, including an 18-, 19-, and 20-year-old female. Female sifakas have an average life expectancy of 5 years (Colchero et al. 2021), and the oldest female in our long-term records was 21 years old. Hence, maternal senescence as described for many other mammals (Ivimey-Cook and Moorad 2020) did not influence infant pre-adult survival. However, since our sample size with 12 infants is small, limiting our analyses to one factor per model, and we had not many young mothers during the study period, the relative importance of each effect remains difficult to assess without additional data.

Conclusions
This study of a wild primate population in a seasonal habitat revealed that ecological and intrinsic effects on maternal GC concentrations are difficult to disentangle because responses to current climatic conditions, reproductive history, and current reproductive activities are superimposed. Maternal GC levels were profoundly linked to infant behavior and development, however. Maternal effects mediated by hormones are therefore prevalent in this species and appear to reveal a trade-off between maternal benefits and infant costs. As such, our study provides mixed evidence for beneficial effects of extended maternal GC levels and provides one of the first studies simultaneously examining the connections between environmental factors, maternal stress, and fitness proxies (Beehner and Bergman 2017;Campos et al. 2021).