www.nature.com/scientificreports OPEN Parasite intensity drives fetal development and sex allocation in a wild ungulate O. Alejandro Aleuy1,2*, Emmanuel Serrano3,4, Kathreen E. Ruckstuhl1, Eric P. Hoberg5 & Susan Kutz2 An understanding of the mechanisms influencing prenatal characteristics is fundamental to comprehend the role of ecological and evolutionary processes behind survival and reproductive success in animals. Although the negative influence of parasites on host fitness is undisputable, we know very little about how parasitic infection in reproductive females might influence prenatal factors such as fetal development and sex allocation. Using an archival collection of Dall’s sheep (Ovis dalli dalli), a capital breeder that depends on its body reserves to overcome the arctic winter, we investigated the direct and indirect impacts of the parasite community on fetal development and sex allocation. Using partial least squares modelling, we observed a negative effect of parasite community on fetal development, driven primarily by the nematode Marshallagia marshalli. Principal component analysis demonstrated that mothers with low parasite burden and in good body condition were more likely to have female versus male fetuses. This association was primarily driven by the indirect effect of M. marshalli on ewe body condition. Refining our knowledge of the direct and indirect impact that parasite communities can have on reproduction in mammals is critical for understanding the effects of infectious diseases on wildlife populations. This can be particularly relevant for species living in ecosystems sensitive to the effects of global climate change. Parasites can negatively affect pregnancy ­rate1–3, offspring ­survival4–6, and offspring ­size7. However, very little is known about how prenatal characteristics, such as fetal development and sex allocation, can be influenced by parasitic infections in the mother. Determining the factors that affect prenatal characteristics is key to under- standing the ecological and evolutionary mechanisms behind survival and reproductive success in ­animals8,9. For instance, fetal development has direct consequences for postnatal physiology, metabolism, growth and immune response (e.g. b­ irds9,10, ­people11,12, wild ­mammals13, and various domestic ­species14–19), while sex allocation can influence population dynamics through a variety of sex-biased mechanisms (e.g. sex-biased mortality and sex- biased reproduction rate)20. Parasites can affect fetus development through several pathways including (i) a direct effect of parasite vir- ulence on fetal development or survival causing pregnancy failure (e.g. abortion, fetal mummification, fetal reabsorption), as seen with bacteria such as Brucella abortus and protozoa such as Neospora caninum21,22, (ii) an indirect effect of the energetic cost of parasitic infection on the mother (e.g. parasite nutrition, host immune response, appetite depression), resulting in decreased nutrition to the fetus and a developmental d ­ elay23,24, and (iii) an indirect effect of parasites interrupting or delaying pregnancy due to a negative effect on host condition resulting in delayed ovulation and conception, which can translate to smaller fetus ­size25. The few attempts to quantify the effects of parasites on fetal development have focused on humans, where the infection with gastro- intestinal helminths can cause delayed fetal growth and premature parturition (reviewed i­ n26). There is theoretical and empirical evidence supporting the hypothesis that sex allocation can be directly influenced by extrauterine ­factors27. The Trivers-Willard ­hypothesis28, for instance, suggests that in polygynous 1 Department of Biological Sciences, University of Calgary, 3330 Hospital Drive, Health Science Centre 2559, Calgary, AB, Canada. 2Department of Ecosystem and Public Health, Faculty of Veterinary Medicine, University of Calgary, Calgary, Canada. 3Wildlife Ecology and Health Group (WE&H) and Servei d’Ecopatologia de Fauna Salvatge (SEFaS), Departament de Medicina i Cirurgia Animals, Universitat Autònoma de Barcelona, Bellaterra, Barcelona, Spain. 4Dipartimento di Scienze Veterinarie, Universitá di Torino, Grugliasco, Turin, Italy. 5Museum of Southwestern Biology and Department of Biology, University of New Mexico, Albuquerque, NM, USA. *email: oaleuy@ucalgary.ca Scientific Reports | (2020) 10:15626 | https://doi.org/10.1038/s41598-020-72376-x 1 Vol.:(0123456789) www.nature.com/scientificreports/ mating systems, mothers in good condition will produce more sons than daughters because they can afford to provide the needed maternal care to produce a high-quality son. This strategy yields the greatest fitness return in species where variance in reproductive success is strongly correlated with body size in males but not in females. Conversely, the Local Resource Competition hypothesis suggests that females in poor condition will produce more sons because sons will disperse, and thus are less likely to compete with the mother for resources during ­adulthood29. Many efforts have been made to understand the physiological mechanisms behind these sex-biased trends and their relationship with different factors, including mother condition and stress (e.g.30–34.). Although the role of mother’s parasitic infection on sex allocation has not been documented in wild mammals, studies investigating the influence of protozoan parasite, Toxoplasma gondii, in sex allocation of people and mice have demonstrated an important link between this parasite and sex a­ djustments35–37. We took advantage of an unprecedented historical collection of parasites and associated data in Dall’s sheep (Ovis dalli dalli, see Fig. 2A) collected from the Mackenzie Mountains, Canada in 1971–197238,39, to assess the effect of gastrointestinal parasites on fetal development and sex allocation. Dall’s sheep are highly philopat- ric, females normally remain in the same general area their entire lives, whereas subadult males disperse from their mothers’ group to join bachelor g­ roups40. Also, these free-ranging wild sheep take advantage of summer pastures to fuel their body reserves in anticipation of winter shortages and, as a capital breeder, use the stored energy to support the cost of reproduction during w ­ inter41–43. Under the current conditions of accelerating global warming, seasonal patterns of climate and vegetation growth are being altered in Alpine and Arctic ecosystems ­worldwide44. Concurrently, patterns of parasitism are changing, with amplification and range expansion of some parasitic species already demonstrated at higher l­atitudes45–47. These changes shape body condition, survival, and reproductive patterns of ungulates living in these e­ nvironments48,49. Yet the synergistic effects of parasitism on body reserves, and ultimately on the reproductive outcomes (fetal development and sex allocation), might go unnoticed. We had two main research objectives: (i) to determine the association of fetal development with the mother’s age, body condition, and infection intensity with gastrointestinal parasites, and (ii) to determine the association of fetus sex allocation with the same extrauterine factors. We hypothesized that, due to the high energetic costs of gastrointestinal parasites to their host, they are an important extrauterine factor determining fetal character- istics in Dall’s sheep. We predicted that, due to their negative impact on host condition, higher parasite burdens would be associated with smaller and lighter fetuses and a higher probability of having daughters than sons. To test this, we used a combination of Partial Least Square Path Modeling approach, an innovative regression procedure used to study causal relationships based on observational data, and traditional regression approaches. Materials and methods Sample collection. Dall’s sheep sampling and data collection. Scientific collections of wild Dall’s sheep occurred between February 17–19, 1971, and February 17–22, 1972, in the Mackenzie Mountains, Northwest Territories (NWT), Canada. This work was led by Dr. Norman Simmons with the Canadian Wildlife Service, with the primary objective of doing a demographic study of ­sheep38. The parasite specimens and the associated data collected from these sheep were located in the spring of 2000 at the Canadian Museum of Nature (CMN) in Ottawa (Ontario, Canada) and are permanently archived as “The Simmons Collection”39. The Simmons Collection contained information from 70 pregnant ewes. The information documented from each ewe included: age in years based on tooth cementum annuli; body weight (kg); body length and chest girth, the species and number of adult helminth parasites in the abomasum, small intestine, large intestine and caecum and; the sex, weight (kg) and total body length (cm) of the fetus. Summary statistics of these variables between 1971 and 1972 can be found in Supplementary Table S1. Data analysis. Fetus development. We used Partial Least Square Path Modeling (PLS-PM), also known as “soft modelling” due to the liberal distribution assumptions and sample size requirements, to explore the as- sociation of gastrointestinal helminths and ewe characteristics (i.e. age and body condition) with fetus length and weight (used as a proxy for fetus development and/or fetus age), in Dall’s sheep. The PLS-PM approach is the intersection of Regression Models, Structural Equations Models, and Multiple Table ­Analysis50,51. Briefly, this approach quantifies the network relationship between a set of unobservable latent variables (LV) and a set of manifest variables (MV, i.e. parameters directly measured in the field or in the laboratory). The LVs are concep- tual variables defined by one or several MVs and organized in a network of relationships where the connections among LVs are assumed to represent a cause-effect process. This network of relationships among LVs forms the inner model (also called structural model) while each group of MVs linked to a LV represents the outer, or measurement, model. The links among LV are quantified through path coefficients while the links between LV and MV are quantified through ­weights52. The analyses included ten MVs organized in the following three LVs (Fig. 1A); (i) age of mother in years, (ii) gastrointestinal parasites of the mother, defined by the intensity of each helminth species with prevalences higher than 10% and by the diversity of gastrointestinal parasite species (i.e. number of species infecting each ewe), and iii) ewe body condition, defined by a scale mass index (SMI) calculated for each ewe following Peig and Green (2009)53. The SMI is a versatile index due to its independence from body size and because it can be used to compare individuals from different ­populations54. The morphometric used to calculate the SMI was the ewe’s chest girth, as girth had the highest correlation with body weight (r = 0.69, p < 0.001, n = 67) among all the parameters available. After fitting the first model including all the variables, a model simplification was performed by removing those MVs uncorrelated with their own LVs. The PLS-PM approach does not depend on any distributional assumptions, therefore, a resampling procedure or bootstrap validation was used to get Scientific Reports | (2020) 10:15626 | https://doi.org/10.1038/s41598-020-72376-x 2 Vol:.(1234567890) www.nature.com/scientificreports/ Figure 1.  Initial path models describing phenotypical characteristics of Dall’s sheep fetuses collected in the Mackenzie Mountains, Canada during 1971 and 1972. (A) Path model for intrauterine fetal development. (B) Path model for fetus sex. Figure created using Microsoft PowerPoint, Version 14.1. confidence intervals for evaluating the precision of the PLS parameter estimates (e.g. path coefficients, total effects and fit indices such as ­R2). Fetus sex. The relationship between fetus sex and extrauterine characteristics in the ewe was investigated in three steps. First, a univariate comparison was performed to determine differences in sex distribution of the fetuses (chi-squared test) and differences in fetus morphometrics (i.e. weight and length; t-tests for differences in means, under the assumption of a normal distribution of means). Second, we performed a principal compo- nent analysis (PCA), including not only the burden of the most common parasite species but also other vari- ables potentially linked to the sex of the fetus (e.g. age of the mother, or body condition of the mother). Then, we evaluated whether the values in the linear combination with the maximum variance (first PCA dimension) differed between female and male fetuses. The SMI index was used as a proxy for body condition of females and the log intensity of each parasite species with prevalence > 10% and/or median intensity of > 10 worms (i.e. M. marshalli, Nematodirus spp., and Trichuris sp.), as parasite burden. We also included the addition of the inten- sity of all parasite species with prevalence < 10% and/or median intensity < 10 worms (i.e. intensity of Moniezia sp. + Trichostrongylus sp. + Skrjabinema ovis) in the PCA analysis. Third, to explore indirect effects of parasites on fetus sex, a PLS-PM analysis was performed including body condition, parasite intensity of the mother (i.e. M. marshalli, Nematodirus spp., and Trichuris sp., and parasites with low intensity/prevalence), and fetus sex as LVs (Fig. 1B). After fitting the first model that included all the variables, a model simplification was performed by removing those MVs uncorrelated with their own LVs. The significant effects observed in the final PLS-PM were confirmed using Generalized Linear Models to fit models among the variables included in each significant effect. All the analyses were performed using R (Version 3.5.2) (R Core Team, 2013. R: A language and environ- ment for statistical computing. R foundation for statistical computing Vienna, Austria) and using the “plspm”55 and “FactoMineR”56 packages. Results Descriptive statistics for the variables used in the full PLS-PM to investigate fetus development are in Table 1. The goodness of fit for the final PLS-PM was 0.31 and explained 14.2% (­ R2) of the observed fetal development variability. Parasite intensity of ewes (i.e. M. marshalli intensity) was negatively associated with fetus development, Scientific Reports | (2020) 10:15626 | https://doi.org/10.1038/s41598-020-72376-x 3 Vol.:(0123456789) www.nature.com/scientificreports/ Descriptive statistics x̅ or median, (range) Prevalence (%), confidence Latent variables Manifest Variables n interval (CI) Ewe age Age (years) 70 Median = 5.8 (1.75- 11.75) Ewe parasite intensity Nematodirus spp. 68 Median = 59 (1–1,158) 92.3 (82.9–97.3) Moniezia sp. 68 Median = 1 (1–3) 13.2 (6.6–24.1) T. schumakovitschi 68 Median = 13.5 (1–89) 85.3 (74.1–92.3) S. ovis 68 Median = 5.5 (1–45) 64.7 (52.1–75.6) M. marshalli 50 Median = 207 (42–2,106) 100 (91.1–100) Species diversity 50 Median = 4 (2–5) Ewe condition Scale Mass Index (SMI) 70 x̅ = 51.1 (37.2–60.21) Fetus development Weight (gr) 70 x̅ = 235.0 (11.4–574.9) Total length (cm) 69 x̅ = 24.4 (10.5–33.7) Table 1.  Latent variables and descriptive statistics of the manifest variables used for fitting the causal model for fetus development and sex allocation of Dall’s sheep from the Mackenzie Mountains, Canada. Median, Median parasite intensity calculated considering only parasite counts equal or higher than 1. x̅, mean calculated using all the values. min, minimum value. Not considering null counts. max, maximum value. explaining the variation of this LV to a greater extent (68%) than mother’s age (30.6%) and mother’s body con- dition (0.8%) (Fig. 2B, Table 2). The correlation of each parasite species with its LV in the full model (Fig. 1A) was very low (λ < 0.3) with the exception of M. marshalli (λ = 0.86), therefore, M. marshalli was the only para- site species retained in the final model (Fig. 2B) (Supplementary Table S2). Although ewe body condition was negatively associated with the LV parasite intensity, there was no significant association between the ewe body condition and fetus development. The bootstrap results indicated that all relevant path coefficients, total effects and fit indices (R2) in the final model were significantly non-random at the ≤ 5% level (Table 3). The sex of the fetuses were evenly distributed with 37 females and 33 males (Chi-squared, X2 = 0.36, df = 1, p value = 0.547). The weights for female (mean = 246.31 g, SE = 20.23) and male (mean = 222.29g, SE = 16.69) fetuses did not differ significantly (t-test, weight, t 0.91, 66, p value = 0.363), nor did the total length of fetuses differ among sexes (Females mean = 24.56 cm, SE = 0.79; Males mean = 24.15 cm, SE = 0.67, t 0.39, 65, p value = 0.693). The infection intensity of almost all parasite species was higher in ewes carrying a male fetus than in ewes with a female fetus. This trend was particularly evident for M. marshalli, N. archari, all Nematodirus spp. together, and all gastrointestinal helminths together (Supplementary Figure S1). The first two principal components (PC) explained 47.7% of the PCA variance. The scores of the first component were significantly different between male and female fetuses ­(R2 = 16%, p = 0.004, Fig. 3). The first PC was significantly associated with, in decreasing values of importance (loading), M. marshalli intensity (loading = 0.76), mother body condition (loading = − 0.72), Nematodirus spp. intensity (loading = 0.47), and Trichuris sp. intensity (loading = 0.29). The final PLS-PM model explained 14.2% ­(R2) of the variability observed in fetus sex, showing that ewes in good body condition and low parasite intensity were more likely to have female fetuses. A significant negative effect of parasites on ewe body condition was also observed (Table 4, Fig. 2C and Supplementary Table S3). In the final model, LV parasites only included M. marshalli (Fig. 2C). Discussion Fetal characteristics are influenced by complex biological interactions between intrauterine (e.g. nutritional exchange through the placenta, position in the uterus) and extrauterine factors (e.g. environment, physiologi- cal condition of the mother, genetic potential of the fetus)14. We found that gastrointestinal helminths, and particularly M. marshalli, were both directly and indirectly associated with lighter and smaller fetuses in Dall’s sheep. Additionally, contrary to what is typically predicted by ecological theory for dimorphic species like Dall’s ­sheep28,57, we observed that ewes in good body condition, which typically have low parasite intensity, were more likely to be carrying female than male fetuses (Fig. 2D). Although the association of parasites with smaller fetuses is very intuitive reports in the wildlife literature are rare. The cost of parasites for wild ungulates can be substantial, particularly in species from environments marked by high seasonality or extreme seasonal environmental conditions such as Dall’s sheep, bighorn sheep, and ­others58,59. Dall’s sheep have relatively low pregnancy rates (~ 75%)38, compared to other wild sheep (e.g. bighorn sheep ~ 90%60,61), with pronounced interannual variation in timing and synchrony of p ­ arturition41. Scientific Reports | (2020) 10:15626 | https://doi.org/10.1038/s41598-020-72376-x 4 Vol:.(1234567890) www.nature.com/scientificreports/ Figure 2.  (A) Dall’s sheep ewe and lamb. Final path model describing, (B) the fetal development and (C) sex allocation of Dall’s sheep. Blue and red arrows represent positive and negative associations between latent variables (LV). Solid and dashed arrows indicate significant and non-significant associations among LVs determined by bootstrapping validation. The ­R2 are coefficients of determination for each latent variable and indicate the amount of variance explained by their independent latent variables. β represents the path coefficient between LV (e.g. β1-2 path coefficient between LV 1 “Ewe age” and LV 2 “Ewe parasite intensity (D) summarizes the results from (B) and (C). Ewes in better condition and with low intensity of M. marshalli carry more female fetuses and those fetuses are larger than in ewes in low body condition and with high M. marshalli intensity. Ewes with low body condition and high M. marshalli intensity carry more male fetuses and their fetuses are smaller. Figure created using Microsoft PowerPoint, Version 14.1. Explanatory variables for Fetus development β Correlation Contribution to R2 (%) Ewe parasite intensity − 0.296 − 0.329 68.6 Ewe age − 0.186 − 0.233 30.6 Ewe condition 0.008 0.135 0.8 Table 2.  Regression coefficients and contribution (%) of each latent variable to the explained variability observed in the PLS-PM model describing fetus development of Dall’s sheep. β = path coefficient estimated by bootstrapping. This may suggest that this species is particularly susceptible to costly energetic factors like parasitic infections, especially during periods of energetic stress like breeding season and pregnancy. The gestation period for Dall’s sheep, which occurs under the extreme cold and dry conditions of the Subarctic and Arctic ­winter41,42, coincides with the seasonal peak of M. marshalli abundance in this s­ pecies62–64. Here we suggest two hypotheses that could explain the association of fetus size with the intensity of parasite infection. First, M. marshalli may have a direct effect on the fetus by decreasing the nutrient availability for development. Marshallagia marshalli causes a variety of structural changes in the abomasum as a result of larval migration through the gastric glands during maturation and the mechanical damage of adult parasites attached to the ­mucosa65,66. This results in an increase in abomasal pH, and an increase in both serum pepsinogen and serum Scientific Reports | (2020) 10:15626 | https://doi.org/10.1038/s41598-020-72376-x 5 Vol.:(0123456789) www.nature.com/scientificreports/ Effects Relationship Direct Indirect Total 95% CI Ewe age—> Parasites 0.162 0 0.162 − 0.04–0.35 Ewe age—> Ewe condition 0 − 0.071 − 0.071 − 0.16–0.01 Ewe age—> Fetal development − 0.186 − 0.048 − 0.234 − 0.5–− 0.04* Parasites—> Ewe condition − 0.437 0 − 0.437 − 0.66–− 0.08* Parasites—> Fetal development − 0.296 − 0.004 − 0.3 − 0.54–− 0.06* Ewe condition—> Fetal development 0.009 0 0.009 − 0.36–0.4 Table 3.  Direct, indirect and total effects among latent variables (LV) in the final PLS-PM describing the development of Dall’s sheep fetuses (n = 49). The 95% confidence intervals were determined by bootstrapping validation. Ewe body condition (SMI) 3 55 B 2000 C 2 * M . marsh al l i Scores PC 1 50 ** 1500 intensity 1 1000 0 45 1 500 40 2 0 Female Male Female Male Female Male 1200 T . sch umak ovi t sch i D E N emat od i rus spp. 1000 80 800 60 intensity intensity 600 40 400 200 20 0 0 Female Male Female Male Figure 3.  Comparison between male and female fetuses of Dall’s sheep. (A) Scores of the first component from the PCA analysis to determine the association between extra uterine characteristics and fetus sex in Dall’s sheep ewe. (B)–(E) Variables included in the first component of the same analysis. Significant differences are indicated with different number of * (p < 0.05). Figure created using R (Version 3.5.2) (R Core Team, 2013. R: A language and environment for statistical computing). Effects Relationship Direct Indirect Total 95% CI Parasites—> Ewes condition − 0.421 0 − 0.421 − 0.614–− 0.0003* Parasites—> Fetus Sex 0.036 0.152 0.187 − 0.04–0.385 Ewe condition—> Fetus sex − 0.36 0 − 0.36 − 0.601–− 0.077* Table 4.  Direct, indirect and total effects among latent variables (LV) in the final PLS-PM describing fetus sex in Dall’s sheep (n = 50). The 95% confidence intervals were determined by bootstrapping validation. gastrin concentration, which in turn decreases nutrient metabolism in the abomasum. This effect of parasites on the metabolism of the gestating female is substantial, particularly in terms of decreased protein availability for the growing f­ etus67. These factors, together with the appetite depression typically produced by ostertagiine ­nematodes68 and the environmental constraints during winter within Dall’s sheep range, could have an additive effect on the host, limiting the resources for fetus development, resulting in smaller fetuses. Not mutually exclusive, in our second hypothesis we suggest an indirect effect of M. marshalli on time of conception, as critical reproductive variables, such as the timing of ovulation, depend on the body condition of the ­mother69. Here, we make the assumption that parasite intensity of ewes in the winter (our study period) is Scientific Reports | (2020) 10:15626 | https://doi.org/10.1038/s41598-020-72376-x 6 Vol:.(1234567890) www.nature.com/scientificreports/ correlated to that of the previous summer. Parasites are rarely distributed evenly within a population due to dif- ferences in exposure, immune response, and tolerance that are individually determined (through genetics and/ or behaviour). This results in individuals that consistently have higher (or lower) infection intensities relative to other i­ ndividuals70,71. Perhaps supporting a hypothesis of delayed conception, the reduced fetal size in ewes with high parasite intensity may be a consequence of parasite-induced energetic constraints and consequent later conception dates resulting in shorter development time for the fetuses of these highly parasitized m ­ others25,72,73. In our previous study, non-pregnant females from the same group of animals had significantly higher parasite infection intensities and lower body condition than pregnant ­females39. Assuming similar gestation time among fetuses, delayed time of conception can result in births occurring later in the season and, as a consequence, to a variety of post-birth costs for both the newborn and the mother. For instance, in bighorn sheep, lambs that are born later in the season have higher mortality than early-born ones, probably because of the energetic cost of a shortened access that late-lambs would have to high-quality forage or/and lactation occurring when the quality of the forage declines (e.g. low quality or insufficient milk production)74,75. We did not observe a direct effect of ewe body condition on fetal development. This may be because, in capital breeders like Dall’s sheep, breeding only occurs when body reserves reach a threshold condition. Heavily parasitized ewes may have taken longer during the previous fall to reach this critical threshold, yet were able to maintain this weight through gestation, thus masking any effect of body condition on fetal development. The association of sex and maternal condition remains as a central paradigm to understanding life-history strategies and evolutionary theory. Our analyses revealed that pregnant Dall’s sheep in better condition, with fewer parasites, were more likely to have a female than a male fetus. This is consistent with the Local Resource Competition theory which suggests that, in philopatric species, mothers in poor condition will produce more sons because they often are the dispersing sex and thus less likely to compete for resources in the f­ uture29. This effect can be enhanced when a strong maternal transmission of condition to daughters leads to fitness maxi- mization by producing good-quality daughters instead of average quality ­sons76. For instance, a high-quality territory can influence reproductive value (fraction of a future population that has descended from a female) more than it does immediate reproductive success (number of offspring produced in a lifetime), therefore it may be beneficial for mothers in good condition to produce a high-quality daughter who can maintain the territory and also produce more high-quality daughters in ­time77–79. Our results are contrary to what is predicted by the Trivers and Willard hypothesis, in which case, Dall’s sheep in good condition would be expected to carry sons instead of daughters. An explanation for this may be found in the extension and generalization of the Trivers and Willard hypothesis suggested by Schindler et al.80. They incorporated male-specific demographic parameters commonly observed in polygynous species (e.g. high mortality rate, different sex-specific reproductive schedules, and more risk-prone life-history tactics) into two-sex modelling approaches and demonstrated that changes in these demographic parameters result in different sex allocation tactics. For instance, in a model parameterized with data from Columbian ground squirrels (Urocitellus columbianus) where the identical mortality rate of sexes was assigned, a reversal of the Trivers and Willard effect was predicted, with mothers in good condition producing more daughters than sons. Reversal of the Trivers and Willard effect has been observed also in bighorn s­ heep81 and several other u ­ ngulates76. This context-dependent framework predicts changes to the optimal sex allocation tactics over a relatively small range of demographic values highlighting the importance of the species-specific reproductive value of daughters or sons in order to determine the optimal sex allocation tactics for a mother. Most of the evidence describing physiological mechanisms for facultative adjustment of sex ratios is related to the role of maternal glucose and stress hormones, yet the results are ambiguous and still far from c­ onsensus82. For instance, high glucose concentrations around conception can influence sex ratios by increasing the survival of male blastocysts through its interaction with luteinizing hormone (LH)31,83,84. Similarly, high-glucose-mediated mechanisms have also been associated with an increased likelihood of mothers in good condition having daugh- ters instead of ­sons32. Although our results are aligned with the latter mechanism, as parasites, through their negative influence on host condition, can lower glucose ­concentrations85, the evidence is far from conclusive and more investigation focused on the physiological, behavioural and ecological mechanisms behind facultative adjustment of sex ratios is needed. As a highly seasonal ungulate, Dall’s sheep are particularly vulnerable to climate warming and its effect on the phenology of many ecological ­processes48,49. While changing climatic conditions can certainly favour the development and transmission of parasites at high latitudes, they can also cause the opposite effect by, for instance, changing parasite (or host) phenology leading to a spatial and/or temporal mismatch between host and ­parasites86–88. Predicting the consequences of climate change on parasite dynamics remains a challenge as complex interactions among host, environment and the parasites need to be considered and a variety of outcomes can be ­expected88,89. Regardless of the direction of these changes in parasite dynamics, our results suggest they can have a variety of subtle consequences in Dall’s sheep population dynamics, not only through a direct influ- ence on fetal development but also through an indirect impact on fetus sex allocation. In this research, we have demonstrated subtle and important mechanisms through which parasites can affect host population dynamics. It is well documented that M. marshalli can negatively impact the fitness of wild and domestic ungulate ­species39,90. Our results suggest that this negative effect also extends to fetal development and fetal sex allocation, via a combination of direct and indirect effects of M. marshalli on the fetus and its mother. These results gain even more relevance in the current context of strong environmental disturbance and acceler- ated changes in climate, as we showed that shifts in parasite dynamics, independent of the direction of these changes, can have direct and indirect consequences on host population dynamics through impacts on fetal development and sex allocation. Refining our understanding of the impact that individual parasite species, as well as parasite communities, can have on host dynamics is critical to comprehend complex ecological processes. Scientific Reports | (2020) 10:15626 | https://doi.org/10.1038/s41598-020-72376-x 7 Vol.:(0123456789) www.nature.com/scientificreports/ Received: 17 January 2020; Accepted: 17 August 2020 References 1. Stien, A. et al. The impact of gastrointestinal nematodes on wild reindeer: experimental and cross-sectional studies. J. Anim. Ecol. 71, 937–945 (2002). 2. Budischak, S. A., O’Neal, D., Jolles, A. E. & Ezenwa, V. O. Differential host responses to parasitism shape divergent fitness costs of infection. Funct. Ecol. 32, 324–333 (2018). 3. Albon, S. D. et al. The role of parasites in the dynamics of a reindeer population. Proc. R. Soc. Lond. B 269, 1625–1632 (2002). 4. Festa-Bianchet, M. Numbers of lungworm larvae in feces of bighorn sheep: yearly changes, influence of host sex, and effects on host survival. Can. J. Zool. 69, 547–554 (1991). 5. Richner, H., Oppliger, A. & Christe, P. Effect of an ectoparasite on reproduction in great tits. J. Anim. Ecol. 62, 703–710 (1993). 6. Fitze, P. S., Tschirren, B. & Richner, H. Life history and fitness consequences of ectoparasites. J. Anim. Ecol. 73, 216–226 (2004). 7. Patterson, J. E. H., Neuhaus, P., Kutz, S. J. & Ruckstuhl, K. E. Parasite removal improves reproductive success of female North American red squirrels (Tamiasciurus hudsonicus). PLoS ONE 8, e55779 (2013). 8. Gilbert, S. F. Ecological developmental biology: developmental biology meets the real world. Dev. Biol. 233, 1–12 (2001). 9. Monaghan, P. Early growth conditions, phenotypic development and environmental change. Philos Trans. R. Soc. Lond. B Biol. Sci. 363, 1635–1645 (2008). 10. Bowers, E. K. et al. Neonatal body condition, immune responsiveness, and hematocrit predict longevity in a wild bird population. Ecology 95, 3027–3034 (2014). 11. Gluckman, P. D., Hanson, M. A., Morton, S. M. B. & Pinal, C. S. Life-long echoes–a critical analysis of the developmental origins of adult disease model. Neonatology 87, 127–139 (2005). 12. Gluckman, P. D., Hanson, M. A. & Beedle, A. S. Early life events and their consequences for later disease: A life history and evo- lutionary perspective. Am. J. Hum. Biol. 19, 1–19 (2007). 13. Lindström, J. Early development and fitness in birds and mammals. Trends Ecol. Evol. 14, 343–348 (1999). 14. Wu, G., Bazer, F. W., Wallace, J. M. & Spencer, T. E. Board-invited review: intrauterine growth retardation: implications for the animal sciences. J. Anim. Sci. 84, 2316–2337 (2006). 15. Greenwood, P. L. & Bell, A. W. Prenatal nutritional influences on growth and development of ruminants. Recent Adv. Animal Nutr. Aust. 14, 57 (2003). 16. Alexander, G. & Williams, D. Heat stress and development of the conceptus in domestic sheep. J. Agric. Sci. 76, 53–72 (1971). 17. Holland, M. D. & Odde, K. G. Factors affecting calf birth weight: a review. Theriogenology 38, 769–798 (1992). 18. Reynolds, L. P., Ferrell, C. L., Nienaber, J. A. & Ford, S. P. Effects of chronic environmental heat stress on blood flow and nutrient uptake of the gravid bovine uterus and foetus. J. Agric. Sci. 104, 289–297 (1985). 19. Johnson, J. S. et al. The impact of in utero heat stress and nutrient restriction on progeny body composition. J. Therm. Biol. 53, 143–150 (2015). 20. Lindström, J. & Kokko, H. Sexual reproduction and population dynamics: the role of polygyny and demographic sex differences. Proc. Biol. Sci. 265, 483–488 (1998). 21. de Figueiredo, P., Ficht, T. A., Rice-Ficht, A., Rossetti, C. A. & Adams, L. G. Pathogenesis and Immunobiology of Brucellosis: Review of Brucella-Host Interactions. Am. J. Pathol. 185, 1505–1517 (2015). 22. Donahoe, S. L., Lindsay, S. A., Krockenberger, M., Phalen, D. & Šlapeta, J. A review of neosporosis and pathologic findings of Neospora caninum infection in wildlife. Int. J. Parasitol. Parasites Wildl. 4, 216–238 (2015). 23. Robbins, C. T. & Robbins, B. L. Fetal and Neonatal Growth Patterns and Maternal Reproductive Effort in Ungulates and Subun- gulates. Am. Nat. 114, 101–116 (1979). 24. Martin, R. D. & MacLarnon, A. M. Gestation period, neonatal size and maternal investment in placental mammals.pdf. Nature 313, 220–223 (1985). 25. O’Callaghan, D. & Boland, M. P. Nutritional effects on ovulation, embryo development and the establishment of pregnancy in ruminants. Anim. Sci. 68, 299–314 (1999). 26. Blackwell, A. D. Helminth infection during pregnancy: insights from evolutionary ecology. Int. J. Womens Health 8, 651–661 (2016). 27. Booksmythe, I., Mautz, B., Davis, J., Nakagawa, S. & Jennions, M. D. Facultative adjustment of the offspring sex ratio and male attractiveness: a systematic review and meta-analysis. Biol. Rev. Camb. Philos. Soc. 92, 108–134 (2017). 28. Trivers, R. L. & Willard, D. E. Natural selection of parental ability to vary the sex ratio of offspring. Science 179, 90–92 (1973). 29. Silk, J. B. Local Resource Competition and Facultative Adjustment of Sex Ratios in Relation to Competitive Abilities. Am. Nat. 121, 56–66 (1983). 30. Ryan, C. P., Anderson, W. G., Gardiner, L. E. & Hare, J. F. Stress-induced sex ratios in ground squirrels: support for a mechanistic hypothesis. Behav. Ecol. 23, 160–167 (2012). 31. Cameron, E. Z. Facultative adjustment of mammalian sex ratios in support of the Trivers-Willard hypothesis: evidence for a mechanism. Proc. Biol. Sci. 271, 1723–1728 (2004). 32. Schwanz, L. E. & Robert, K. A. Proximate and ultimate explanations of mammalian sex allocation in a marsupial model. Behav. Ecol. Sociobiol. 68, 1085–1096 (2014). 33. Silk, J. B. & Brown, G. R. Local resource competition and local resource enhancement shape primate birth sex ratios. Proc. Biol. Sci. 275, 1761–1765 (2008). 34. Ruckstuhl, K. E., Colijn, G. P., Amiot, V. & Vinish, E. Mother’s occupation and sex ratio at birth. BMC Public Health 10, 269 (2010). 35. Flegr, J. & Kaňková, Š. The effects of toxoplasmosis on sex ratio at birth. Early Hum. Dev. 141, 104874 (2020). 36. Kanková, S. et al. Women infected with parasite Toxoplasma have more sons. Naturwissenschaften 94, 122–127 (2007). 37. Kanková, S. et al. Influence of latent toxoplasmosis on the secondary sex ratio in mice. Parasitology 134, 1709–1717 (2007). 38. Simmons, N. M., Bayer, M. B. & Sinkey, L. O. Demography of Dall’s Sheep in the Mackenzie Mountains Northwest Territories. J. Wildl. Manage 48, 156–162 (1984). 39. Aleuy, O. A. et al. Diversity of gastrointestinal helminths in Dall’s sheep and the negative association of the abomasal nematode, Marshallagia marshalli, with fitness indicators. PLoS ONE 13, e0192825 (2018). 40. Geist, V. Mountain Sheep: A Study in Behavior and Evolution (University of Chicago Press, Chicago, 1971). 41. Rachlow, J. L. & Bowyer, R. T. Interannual Variation in Timing and Synchrony of Parturition in Dall’s Sheep. J. Mammal. 72, 487–492 (1991). 42. Goodrowe, K. L., Smak, B., Presley, N. & Nlonfort, S. L. Reproductive, behavioral, and endocrine characteristics of the Dall’s Sheep (Ovis dalli). Zoo Biol. 15, 45–54 (1996). 43. Bunnell, F. L. & Nichols, L. Natural history of thinhorn sheep. In Mountain sheep of North America (ed. Valdez, R.) 23–77 (Uni- versity of Arizona Press, Arizona, 1999). 44. Ernakovich, J. G. et al. Predicted responses of arctic and alpine ecosystems to altered seasonality under climate change. Glob. Chang. Biol. 20, 3256–3269 (2014). Scientific Reports | (2020) 10:15626 | https://doi.org/10.1038/s41598-020-72376-x 8 Vol:.(1234567890) www.nature.com/scientificreports/ 45. Kutz, S. J. et al. Invasion, establishment, and range expansion of two parasitic nematodes in the Canadian Arctic. Glob. Chang. Biol. 19, 3254–3262 (2013). 46. Kutz, S. J. et al. The Arctic as a model for anticipating, preventing, and mitigating climate change impacts on host–parasite interac- tions. Vet. Parasitol. 163, 217–228 (2009). 47. Altizer, S., Ostfeld, R. S., Johnson, P. T. J., Kutz, S. & Harvell, C. D. Climate change and infectious diseases: from evidence to a predictive framework. Science 341, 514–519 (2013). 48. Parker, K. L., Barboza, P. S. & Gillingham, M. P. Nutrition integrates environmental responses of ungulates. Funct. Ecol. 23, 57–69 (2009). 49. Pettorelli, N., Pelletier, F. & von Hardenberg, A. Early onset of vegetation growth vs. rapid green-up: impacts on juvenile mountain ungulates. Ecology 88(2), 381–390 (2007). 50. Sanchez, G. PLS Path Modeling with R. (Trowchez Editions, Berkeley, 2013). http://www.gasto​nsanc​hez.com/PLSPa​thMod​eling​ withR​.pdf. 51. Tenenhaus, M., Vinzi, V. E., Chatelin, Y.-M. & Lauro, C. PLS path modeling. Comput. Stat. Data Anal. 48, 159–205 (2005). 52. Hair, J. F., Ringle, C. M. & Sarstedt, M. Partial least squares structural equation modeling: rigorous applications, better results and higher acceptance. Long Range Plann. 46, 1–12 (2013). 53. Peig, J. & Green, A. J. New perspectives for estimating body condition from mass/length data: the scaled mass index as an alterna- tive method. Oikos 118, 1883–1891 (2009). 54. Labocha, M. K., Schutz, H. & Hayes, J. P. Which body condition index is best?. Oikos 123, 111–119 (2014). 55. Sanchez, G., Trinchera, L. & Russolillo, G. plspm: Tools for partial least squares path modeling (PLS-PM). R package version 0.4. https​://doi.org/10.1111/j.1600-0706.2013.00755​.x (2017). 56. Lê, S., Josse, J., Husson, F. Facto. & Mine, R. An R Package for multivariate analysis. J. Stat. Softw. https:​ //doi.org/10.18637/​ jss.v025. i0 (2008). 57. Clutton-Brock, T. H., Albon, S. D. & Guinness, F. E. Maternal dominance, breeding success and birth sex ratios in red deer. Nature 308, 358–360 (1984). 58. De Roos, A. M., Galic, N. & Heesterbeek, H. How resource competition shapes individual life history for nonplastic growth: ungulates in seasonal food environments. Ecology 90, 945–960 (2009). 59. Festa-Bianchet, M. Individual Differences, Parasites, and the Costs of Reproduction for Bighorn Ewes (Ovis canadensis). J. Anim. Ecol. 58, 785–795 (1989). 60. Festa-Bianchet, M., Jorgenson, J. T. & Wuhart, W. D. Early weaning in bighorn sheep, Ovis canadensis affects growth of males but not of females. Behav. Ecol. 5, 21–27 (1994). 61. Singer, F. J., Williams, E., Miller, M. W. & Zeigenfuss, L. C. Population Growth, Fecundity, and Survivorship in Recovering Popula- tions of Bighorn Sheep. Restor. Ecol. 8, 75–84 (2000). 62. Simmons, N. M. Seasonal Ranges of Dall’s Sheep, Mackenzie Mountains Northwest Territories. Arctic 35, 512–518 (1982). 63. Neilsen, C. & Neiland, K. Sheep Disease Report, Project Progress Report, Federal Aid in Wildlife Restoration. (1974). 64. Kutz, S. J. et al. Chapter 2: parasites in ungulates of Arctic North America and Greenland—a view of contemporary diversity, ecology, and impact in a world under change. In Adv Parasit (ed. Rollinson, D.) 99–252 (Academic Press, Cambridge, 2012). 65. Moradpour, N., Borji, H., Razmi, G., Maleki, M. & Kazemi, H. The effect of Marshallagia marshalli on Serum Gastrin concentra- tions in experimentally infected lambs. J. Parasitol. 102, 436–439 (2016). 66. Moradpour, N., Borji, H., Razmi, G., Maleki, M. & Kazemi, H. Pathophysiology of Marshallagia marshalli in experimentally infected lambs. Parasitology 140, 1762–1767 (2013). 67. Simcock, D. C. et al. Hypergastrinaemia, abomasal bacterial population densities and pH in sheep infected with Ostertagia cir- cumcincta. Int. J. Parasitol. 29, 1053–1063 (1999). 68. Jacobs, D., Fox, M., Gibbons, L. & Hermosilla, C. Principles of Veterinary Parasitology (Wiley, Hoboken, 2015). 69. Berger, T. Fertilization in ungulates. Anim. Reprod. Sci. 42, 351–360 (1996). 70. Hayward, A. D. Causes and consequences of intra- and inter-host heterogeneity in defence against nematodes. Parasite Immunol. https​://doi.org/10.1111/pim.12054​ (2013). 71. Hayward, A. D. et al. Natural selection on individual variation in tolerance of gastrointestinal nematode infection. PLoS Biol. 12, e1001917 (2014). 72. Reimers, E. Growth rate and body size differences in Rangifer, a study of causes and effects. Rangifer 3, 3–15 (1983). 73. Sontakke, S. D. Monitoring and controlling ovarian activities in wild ungulates. Theriogenology 109, 31–41 (2018). 74. Festa-Bianchet, M. Birthdate and survival in bighorn lambs (Ovis canadensis). J. Zool. 214, 653–661 (1988). 75. Feder, C., Martin, J. G. A., Festa-Bianchet, M., Bérubé, C. & Jorgenson, J. Never too late? Consequences of late birthdate for mass and survival of bighorn lambs. Oecologia 156, 773–781 (2008). 76. Hewison, A. J. M. & Gaillard, J.-M. Successful sons or advantaged daughters? The Trivers-Willard model and sex-biased maternal investment in ungulates. Trends Ecol. Evol. 14, 229–234 (1999). 77. Leimar, O. Life-history analysis of the Trivers and Willard sex-ratio problem. Behav. Ecol. 7, 316–325 (1996). 78. Sheldon, B. C. & West, S. A. Maternal dominance, maternal condition, and offspring sex ratio in ungulate mammals. Am. Nat. 163, 40–54 (2004). 79. Julliard, R. Sex-specific dispersal in spatially varying environments leads to habitat-dependent evolutionary stable offspring sex ratios. Behav. Ecol. 11, 421–428 (2000). 80. Schindler, S. et al. Sex-specific demography and generalization of the Trivers-Willard theory.PDF. Nature 526, 249–252 (2015). 81. Festa-Bianchet, M. Offspring sex ratio studies of mammals: Does publication depend upon the quality of the research or the direc- tion of the results?. Écoscience 3, 42–44 (1996). 82. Douhard, M. Offspring sex ratio in mammals and the Trivers-Willard hypothesis: In pursuit of unambiguous evidence. Bioessays 39(9), 1700043 (2017). 83. Larson, M. A., Kimura, K., Michael Kubisch, H. & Michael Roberts, R. Sexual dimorphism among bovine embryos in their ability to make the transition to expanded blastocyst and in the expression of the signaling molecule IFN-τ. Proc. Natl. Acad. Sci. U. S. A. 98, 9677–9682 (2001). 84. Cameron, E. Z., Lemons, P. R., Bateman, P. W. & Bennett, N. C. Experimental alteration of litter sex ratios in a mammal. Proc. Biol. Sci. 275, 323–327 (2008). 85. Shea-Donohue, T., Qin, B. & Smith, A. Parasites, nutrition, immune responses and biology of metabolic tissues. Parasite Immunol. 39, e12422 (2017). 86. Lafferty, K. D. The ecology of climate change and infectious diseases. Ecology 90, 888–900 (2009). 87. Kutz, S. J., Hoberg, E. P., Molnár, P. K., Dobson, A. & Verocai, G. G. A walk on the tundra: Host–parasite interactions in an extreme environment. Int. J. Parasitol. Parasites Wildl. 3, 198–208 (2014). 88. Hoar, B. M., Ruckstuhl, K. & Kutz, S. Development and availability of the free-living stages of Ostertagia gruehneri, an abomasal parasite of barrenground caribou (Rangifer tarandus groenlandicus), on the Canadian tundra. Parasitology 139, 1093–1100 (2012). 89. Rose, H., Hoar, B., Kutz, S. J. & Morgan, E. R. Exploiting parallels between livestock and wildlife: Predicting the impact of climate change on gastrointestinal nematodes in ruminants. Int. J. Parasitol. Parasites Wildl. 3, 209–219 (2014). 90. Morgan, E. R. et al. Assessing risks of disease transmission between wildlife and livestock: The Saiga antelope as a case study. Biol. Conserv. 131, 244–254 (2006). Scientific Reports | (2020) 10:15626 | https://doi.org/10.1038/s41598-020-72376-x 9 Vol.:(0123456789) www.nature.com/scientificreports/ Acknowledgements We would like to thank the Killam Program and NSERC CREATE Host-Parasite Interaction Program for fund- ing O. Alejandro Aleuy. This research was supported by NSERC Discovery grant and NSERC Research Tools and Instruments grants to S Kutz. ES was funded by the Spanish Ministerio de Economia y Competitividad (MINECO) through a Ramon y Cajal agreement (RYC-2016-21120). KER acknowledges the financial support of the Natural Sciences and Engineering Research Council of Canada (NSERC), [Funding Reference Number 316189-2012-RGPIN]. Cette recherche a été financée par le Conseil de recherches en sciences naturelles et en génie du Canada (CRSNG), [numéro de référence 316189-2012-RGPIN]. Initial assessments of the Simmons Collections were supported by the Government of the NWT (Alasdair Veitch) and through the Beringian Coevo- lution Project (BCP), funded by the National Science Foundation (DEB 0196095 and 0415668), an exploration of parasite diversity and history across the northern roof of the world, planned and coordinated by Joseph A. Cook (University of New Mexico) and Eric P. Hoberg (former United States National Parasite Collection). Author contributions O.A.A., K.E.R., S.K., E.S. conceptualized and designed the study. S.K. and E.P.H. performed parasite identification and quantification. O.A.A. and E.S. performed statistical analyses and interpretation of results. O.A.A. wrote the first draft of the manuscript. All the authors contributed to writing the final draft. The data generated during and/or analyzed during the current study are available from the corresponding author on reasonable request. Competing interests The authors declare no competing interests. Additional information Supplementary information is available for this paper at https​://doi.org/10.1038/s4159​8-020-72376​-x. Correspondence and requests for materials should be addressed to O.A.A. Reprints and permissions information is available at www.nature.com/reprints. Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creat​iveco​mmons​.org/licen​ses/by/4.0/. © The Author(s) 2020 Scientific Reports | (2020) 10:15626 | https://doi.org/10.1038/s41598-020-72376-x 10 Vol:.(1234567890)