The physiological cost of male-biased parasitism in a nearly monomorphic mammal

Background Even though male-biased parasitism is common in mammals, little effort has been made to evaluate whether higher parasitic burden in males results in an extra biological cost, and thus a decrease in fitness. Body condition impairment and the augmentation of oxidative stress can be used as indicators of the cost of parasite infections. Here, we examined relationships between gastrointestinal and respiratory helminths, body condition and oxidative stress markers (glutathione peroxidase, paraoxonase-1) in 28 Pyrenean chamois (Rupicapra p. pyrenaica) sampled in autumn. Results Only male chamois showed a reduction in body condition and higher oxidative stress due to parasite infection, likely because of the extremely high parasite burdens observed in males. Conclusions This study made evident a disparity in the physiological cost of multiple parasitism between sexes in a wild mammal, mainly due to parasitic richness. Because of the similar life expectancy in male and female chamois, we suggest that males may have developed natural mechanisms to compensate for higher parasite loads during the rut.


Background
Sex-based differences in exposure, susceptibility [1] and tolerance lead to differences in parasite prevalence [2], intensity [3], and pathology [4] in a broad range of vertebrate species. Males are often more infected than females as observed in birds [5], rodents, bats [6], ungulates [1,3] and humans [4]. This male-biased parasitism has been linked to both hormonal and behavioural differences. Generally, males invest more energy than females in the development of testosterone-mediated traits, such as secondary sexual traits or courtship displays. Although testosterone is necessary for developing of secondary sex characteristics in males, high testosterone levels have been linked to a suppressed immune system [4,7] and increased parasitism. Malebiased infection is common in sexually dimorphic and polygynous species [8], and might reduce male fitness. However, in sexually dimorphic species it might be difficult to separate the physiological costs of parasitism from the costs of the production and maintenance of sexually-selected traits [9]. Monomorphic species, in which males and females are morphologically similar, provide ideal systems for examining the physiological costs of male-biased parasitism.
Here we compare the physiological costs of malebiased parasitism in a nearly monomorphic ungulate (only small differences in body mass can be observed in autumn), the Pyrenean chamois (Rupicapra pyrenaica pyrenaica). Found in the mountains of southern Europe, chamois are medium-sized polygynous ungulates [10]. Males are more heavily parasitized [3], however, life expectancy is similar between the sexes [10].
We assess the physiological costs of infection using both general indicators of nutritional status, such as indicators of body condition (e.g. kidney fat stores, KFs [11]) and insulin-like growth factor 1 (IGF1), and oxidative stress (OS) markers, such as paraoxonase 1 (PON1), total antioxidant capacity (TAC) and glutathione peroxidase (GPX). Oxidative stress originates from an imbalance between the production of reactive oxygen or nitrogen species, and the capacity of the organism to compensate for their detrimental effects or destroy these harmful products [12]. As result, OS biomarkers are widely used to measure infection damage in veterinary sciences [13] and conservation physiology [14] and are considered good animal health indicators of animal health. The proper assessment of oxidative status must consider not only the concentration of antioxidants (e.g. antioxidant vitamins or endogenous enzymes), but also oxidative damage. For example, an elevation in endogenous antioxidants with stable oxidative damage concentrations suggests that individuals are successfully dealing with oxidative stress. In contrast, the depletion of antioxidant substances and stable oxidative conditions suggest oxidative damage.
In this work, we assess the endoparasite (gut and respiratory helminths) and physiological status (oxidative stress and nutritional status) of 28 adult Pyrenean chamois. Using these data, we tested whether, for a given parasite load, males have a higher physiological cost than females. Because male chamois are prone to suffer from malnourishment due to rut-induced hypophagia (mainly territorial individuals, see [15]), we expect higher physiological costs (i.e. lower nutritional status and higher oxidative stress) in males than in females for a given parasite burden. On the other hand, because both female and male chamois have similar life expectancies we expect that males may have developed physiological adaptations to compensate for the damage caused by parasites during the rutting period.

Chamois sampling
The respiratory and gastrointestinal systems of 28 adult Pyrenean chamois, 17 females (mean age 9.2 years, range 3-21) and 11 males (mean age 7.5 years, range 3-12) from the Freser-Setcases National Game Reserve (Catalan Pyrenees, northeast Spain) were collected. Animals were hunter-harvested from October to December in 2012 and 2013, coinciding with the rut and a period of poor food quality [16]. Once the animal was dead, heart puncture was performed to obtain blood samples, which were allowed to clot at room temperature and centrifuged at 1800 G. Serum was transported at 4°C and frozen at −80°C until analysis. Abdominal and thoracic viscera were extracted and kidneys were stored with their fat. The organs were placed in labelled plastic bags and transported in a refrigerated box at 4°C to the laboratory, where the material was stored at −20°C until parasitological examination. Kidney fat reserves were used as a proxy for body condition in sampled chamois [11].

Parasitological data
Adult nematodes were recovered from the lungs, preserved in 70% ethanol and later mounted in lactophenol for identification [17]. Lung larvae were obtained using the Baermann method on 30 g of lung parenchyma. Liquid was collected and centrifuged (800× g). Larvae were identified under a microscope using a Favatti chamber [18].

Biochemical analysis
PON-1 activity was determined using p-nitrophenyl acetate as a substrate [19]. Serum paraoxonase (PON1) is an enzyme able to neutralise lipid-peroxides (oxidised lipids) that cause inflammation among other organic and inorganic substrates. The TAC (ability of antioxidants to clear harmful free radicals) was measured using a method based on 2.2′-azinobis-(3-ethylbenzothiazoline-6-sulfonate) decolorization by antioxidants [20]. TAC considers the cumulative action of all the antioxidants present in body fluids. Thus, a low TAC can be indicative of oxidative stress or increased susceptibility to oxidative damage. GPX was measured according to a previously described method [21]. The main biological role of GPX is to reduce lipid hydroperoxides to their corresponding alcohols and free hydrogen peroxide to water. All three analyses were measured in serum with an automated biochemistry analyser (Olympus AU600; Olympus Europe GmbH, Hamburg, Germany). IGF-1 is a hormone that plays an important role in the growth of juveniles and anabolism in adults and it was analysed with an automated solid-phase, enzyme-labelled chemiluminescent immunometric assay (Immulite System, Siemens Health Diagnostics, Deerfield, IL, USA). All analyses showed inter-and intra-assay coefficients of variation below 15%, and recovery ranged between 80 and 120% in all cases.

Statistical analyses
The lung nematode load obtained in the present study was analysed in combination with the data on gastrointestinal parasite load recorded for the same individuals and previously published by Martínez-Guijosa et al. [3].
The partial least squares (PLS) regression method was used to assess the relationships between parasitism and health in female and male chamois. As a response variable, we defined a "Cost of infection block", which included markers of oxidative stress (GPX, PON1), antioxidant capacity (TAC), and the nutritional status of individuals (IGF-1 and KFs). The endoparasite infection block (set of predictor variables) was represented by the intensity of each parasite species (number of parasite individuals of each species), the total intensity (including all species) and the species richness of lung and gastrointestinal parasites. The PLS regression is suitable when sample size is low relative to the number of variables and when collinearities occur among predictor variables [22]. The loads, weights (i.e. relative contribution of each variable to the derived factors), and cross-correlations (e.g. the correlation between each X variable and the whole Y component, i.e. PLS' scores, and vice versa) were also estimated. Finally, the Stone-Geisser test (Q 2 ) at a threshold of 0.0975 was used to assess the fit of the PLS model. The Q 2 statistic can be viewed as a jackknife version of the R 2 statistic. Analyses were performed with the "plspm" package [23] in R, version 3.3.2 [24].

Results
Similar to gastrointestinal parasites [3], prevalence and intensity of lung nematodes and their larvae were higher in males than in females (Table 1). Regarding the health indicators, however, only the PON1 concentration was significantly different between sexes after a Bonferroni correction (t (24.8) = 5.03, P-value < 0.001, Table 2).
After a preliminary analysis, we decided to exclude the IGF-1 and TAC from the Y block due to their low contribution (<2%) to the health block. Along the same lines, predictor variables with little or no influence on the X's component (e.g. burdens of specific parasite species, such as cestodes or Cooperia oncophora) were also discarded to simplify our PLS model. In female chamois, the health block showed no correlation with parasite infection (Q2 = -0.57). In contrast, the male group analysis revealed a marked correlation (5.7 times greater than females; Q2 = 0.22), where 57.1% of the observed variability of chamois health status was accounted for by parasite infections. As a result, gastrointestinal helminth richness, respiratory nematode richness and total gastrointestinal intensity, in decreasing order of importance, had the greatest negative effect on the Y block (Table 3, Fig. 1). In the case of females, only 1.9% of the health status was related to parasite infections.

Discussion
Chamois with low parasitic intensities and richness showed higher PON1 and GPX concentrations (for a given TAC concentration) and were in better body condition than their  [25], we found that multiple infections have a greater physiological cost than single-parasite infections. However, interpreting sex differences in the physiological cost of parasitism is more challenging for several reasons, the first being the fact that females were by far less parasitized than males. In males, the poor body condition and reduced antioxidant defences in highly parasitized individual chamois suggest an overload in the host's capacity to overcome the infestation. Lower parasite burdens and richness show higher levels of antioxidant biomarkers (PON1 and GPX), resulting from up-regulation of endogenous antioxidants or mobilization of exogenous antioxidants, along with better body condition. Nonetheless, both sexes have similar lifespan, suggesting that male chamois employ compensatory mechanisms to deal with the extracost of parasite infections during the rut season. Females, however, do not seem to need these physiological compensations at this time of the year.
Many parasite infections in wildlife tend to be relatively benign, remaining subclinical until the host experiences additional stress [26]. Thus, many of the pathogenic mechanisms would be driven by the host's own immune system, which is also energetically demanding. In these cases, when the cost of maintaining resistance outweighs the benefits, it has been postulated that tolerance is preferred [27]. This strategy, based on reducing the fitness or health cost of an increasing parasite burden, has been described as naturally selected in ungulates. It is especially advantageous when resources must be allocated in confronting more severe threats [28]. Hence, male chamois  may need to increase their tolerance (through raising PON1 and GPX concentrations) during the rut period, when parasite burdens increase due to the immunosuppressive influence of their reduced feeding coupled with the high testosterone levels. Males could then invest their resources in eliminating parasites or reducing their burdens, and repairing tissue lesions when the rut period is finished, their appetite returns to normal and their immune system is relieved of the aforementioned challenges. The increases in body mass that males show in summer [29] may also contribute to compensate a further cost of parasitism during the rut. The lack of any correlation between health status and parasite infections in females could be explained by their low parasite burden at this time of year (Table 1 and [3]). Some evidence supports this hypothesis. Hypophagia displayed by males during the rutting period, especially in territorial individuals [15], would diminish the probability of introducing new infective (L3) larvae from contaminated grass. Additionally, L3 larvae of most gastrointestinal nematode species in non-naïve Caprinae take a minimum of approximately 3-4 weeks to reach adulthood [30] and, if ingested in late summer/early fall, will most likely be hypobiotic for months unless reawakened by a stressor [31]. Hence, both hypophagia and the long time required by larvae to reach adulthood, may suggest that the higher parasite richness and burden of males at that period is not only due to high testosterone concentrations, but the consequence of previously acquired infections. In fact, only one female displayed a parasite burden (2209 nematodes) close to the mean value of males (2797 nematodes). Thus, this work demonstrates that being a male has an extra physiological cost due to sex-biased parasitism that is most likely exacerbated in the autumn. This situation results from the combination of multiple parasite infections and a debilitated immune system due to malnourishment and high levels of testosterone during rut. We hypothesized that the correlation among health parameters and parasite burden is most likely initiated when tolerance to parasite infection is exceeded.  Table 3. Bubble size indicates total (gut and lung) nematode richness ranging from 0 to 19 nematode species. Bubble colour indicates nematode burden. Chamois silhouettes represent a male chamois with low parasite load and hence in good condition and with a good antioxidant response (top left) and a highly parasitized male chamois in poorer condition and lower antioxidant response (bottom right) Other hypotheses should also be considered. Individual male chamois differ extensively in their energy investment in rutting activities in specific years [32]. Hence, many individuals may not be exposed to some of the detrimental effects, therefore, reducing the overall differences in mortality between males and females. Moreover, females could have their own health depletion period, such as during lactation or pregnancy, differentiated in time from the male depletion periods. This could lead males and females to be under similar stress pressures over a year-long period.

Conclusions
Despite the moderate sample size, this study showed a sex-biased physiological cost of multiple infections by gastrointestinal and respiratory nematode parasites in a wild mammal. Interestingly, this extra-cost is mainly related to the number of different parasite species and, to a lesser extent, with the intensity of infection. The similar lifespan of females and males of this mammal suggest that the latter can deal with this cost of infection, at least during the rutting period. Further research should investigate whether this sex-biased cost of infection persists over the year and whether females suffer from these consequences of parasite infection at other times of year. On the other hand, individual-based information will be required to better understand the mechanisms driving the sex-biased parasitism, its physiological cost, and in turn, survival probabilities of this mountain ungulate.