Metabarcoding in two isolated populations of wild roe deer (Capreolus capreolus) reveals variation in gastrointestinal nematode community composition between regions and among age classes
Parasites & Vectors volume 14, Article number: 594 (2021)
Gastrointestinal nematodes are ubiquitous for both domestic and wild ungulates and have varying consequences for health and fitness. They exist as complex communities of multiple co-infecting species, and we have a limited understanding of how these communities vary in different hosts, regions and circumstances or of how this affects their impacts.
We have undertaken ITS2 rDNA nemabiome metabarcoding with next-generation sequencing on populations of nematode larvae isolated from 149 fecal samples of roe deer of different sex and age classes in the two isolated populations of Chizé and Trois Fontaines in France not co-grazing with any domestic ungulate species.
We identified 100 amplified sequence variants (ASVs) that were assigned to 14 gastrointestinal nematode taxa overall at either genus (29%) or species (71%) level. These taxa were dominated by parasites classically found in cervids—e.g. Ostertagia leptospicularis, Spiculopteragia spp. Higher parasite species diversity was present in the Trois Fontaines population than in the Chizé population including the presence of species more typically seen in domestic livestock (Haemonchus contortus, Bunostomum sp., Cooperia punctata, Teladorsagia circumcincta). No differences in parasite species diversity or community composition were seen in the samples collected from three zones of differing habitat quality within the Chizé study area. Young roe deer hosted the highest diversity of gastrointestinal nematodes, with more pronounced effects of age apparent in Trois Fontaines. The effect of host age differed between gastrointestinal nematode species, e.g. there was little effect on O. leptospicularis but a large effect on Trichostrongylus spp. No effect of host sex was detected in either site.
The presence of some livestock parasite species in the Trois Fontaines roe deer population was unexpected given the isolation of this population away from grazing domestic livestock since decades. Overall, our results illustrate the influence of host traits and the local environment on roe deer nemabiome and demonstrate the power of the nemabiome metabarcoding approach to elucidate the composition of gastrointestinal nematode communities in wildlife.
Gastrointestinal nematode parasites of wildlife and livestock have a worldwide distribution. Multiple species usually co-infect a single host, and they have major ecological, economic and animal welfare impacts [1,2,3]. Gastrointestinal nematodes of wildlife can have negative impacts on the health and individual fitness of hosts  as well as on their population dynamics , which, in extreme cases, can lead to population extinction . Consequently, understanding host-parasite relationships and determining parasite diversity and abundance are crucial for sustainable wildlife management [7, 8].
The immunological response can help the host to control parasite infestations but it incurs a significant cost to its resources and must be balanced with other energy demanding functions [4, 9], especially when the quality of habitat is poor . Therefore, the physiological and behavioral attributes of individual hosts influence infection level, parasite community composition and their effects [11, 12]. The age, sex and reproductive status of hosts are important factors in this context and determine the parasite community structure, i.e. the abundance and diversity of parasites.
As their immune system is still immature, the highest abundance and diversity of parasites are generally observed in young individuals (e.g. roe deer, Capreolus capreolus ; Dall’s sheep, Ovis dalli dalli ; feral sheep, Ovis aries [3, 15]). Similarly, because of the deterioration of their immune system, senescent individuals usually host a higher abundance and diversity of parasites than younger adults [16, 17]. Variation in parasite abundance and diversity is also observed between sexes in various species, including ungulates (e.g. roe deer , African buffalo, Syncerus caffer , mouflon, Ovis gmelini musimon × Ovis sp. , feral sheep ), and is generally higher in males . The higher susceptibility of males to infection has been ascribed to the immunosuppressive effects of sex steroid hormones in males , sexual size dimorphism  and also increased likelihood of pathogen exposure . In addition to the influences of life history traits, host populations inhabiting different areas can be infected by various parasites because of the direct effect of environment on parasite communities (e.g. temperature, precipitations)  or via the effect of environment on hosts .
No relation (e.g. moose, Alces alces)  or negative relationship (e.g. chamois, Dall’s sheep) [14, 27] on the health and individual fitness of host with increased parasite richness were reported. This might be explained by the nematode community composition with different levels of pathogenicity among parasites, even among closely phylogenetically related species, and may have synergistic or antagonistic interactions [18, 28,29,30]. Thus, refining the identification of the parasite communities would help to improve our knowledge on host-parasite relationships. However, investigating parasite diversity in wildlife can be challenging . Recently, next-generation sequencing metabarcoding approaches have been developed to investigate gastrointestinal nematode communities [31, 32], although rarely applied to wild ungulates so far [33, 34]. Metabarcoding has a number of advantages including the study of endoparasite communities, without prior assumption of which species are likely to be present. It also allows large numbers of samples to be processed simultaneously and can be highly sensitive to detect rare species . Metabarcoding does not quantify the absolute number of parasites in a sample, but provides information on the proportions of each parasite species.
The aim of this study was to investigate the gastrointestinal nematode communities present in two contrasting natural populations of roe deer (Capreolus capreolus), a common and widely distributed ungulate species in Europe. The gastrointestinal nematodes of roe deer have been previously investigated in wild populations [36,37,38,39] and have been succinctly described in one of the populations studied here . The two populations of roe deer examined in this study have been spatially isolated in fenced areas for more than 60 years, preventing ungulate movement. Therefore, we expected the diversity of nematodes to be dominated by roe deer-specific parasites (H1). We also expected to observe a higher parasitic diversity in young individuals of a few months of age than in adults (H2) and in males than in females (H3). Because differences in immune responses and differences of body conditions were reported between Chizé and Trois Fontaines , we predicted differences of nemabiome with a higher diversity of gastrointestinal nematodes in roe deer living in poor (i.e. Chizé) compared with good environmental conditions (i.e. Trois Fontaines). We predicted similar results at local scale, among the different sectors of Chizé , with a higher diversity in the sector with the poorest environmental conditions (H4).
Materials and methods
The two study sites, Chizé and Trois Fontaines, are located in western (46°05’ N, 0°25’ W) and northeastern France (48°43’ N, 2°61’ W), respectively. Both sites are fenced forests of 2614 ha in Chizé (since 1978) and 1360 ha in Trois Fontaines (since 1976), preventing roe deer and other ungulates to move into or out of the study sites. No ungulates other than roe deer are present within the two sites, except wild boar. Although some rare hare Lepus europaeus inhabit the site of Chizé, there are no lagomorphs in Trois Fontaines (Office Français de la Biodiversité, personal communication). The two sites have different climate, soil and forest attributes. The oceanic climate in Chizé is characterized by mild winters and dry summers in contrast to Trois Fontaines where the continental climate results in rainy warm summers and cold winters. In both sites, forests are dominated by oak (Quercus sp.), beech (Fagus sylvatica) and hornbeam (Carpinus betulus), but with a more heterogeneous spatial distribution of species and coppice abundance and quality in Chizé  compared to Trois Fontaines. In Chizé, the variation of coppice type is categorized into three sectors characterized by a rich, medium or poor environment for roe deer (Fig. 1). With its low productive forest, Chizé is still considered a poor habitat for roe deer, in contrast to the highly productive forest of Trois Fontaines . In fact, roe deer in Trois Fontaines have access to high quantity and quality of food in the whole study area, while roe deer in Chizé experience a trade-off between high quantity or quality of forage . Consequently, in Trois Fontaines, the population of roe deer is increasing, whereas the population growth rate is close to 0 in Chizé . Both populations are kept stable by yearly removals (Office Français de la Biodiversité, personal communication).
Roe deer data
Roe deer of known age are captured annually as part of a long-term monitoring program. Briefly, roe deer were caught using long nets surrounding forest patches and beaters and maintained in wooden boxes until the end of the drive. Details of the capture procedure are available in Gaillard et al. . Fecal samples were collected rectally from immobilized roe deer captured between January and March in 2018 and 2019 under a protocol approved by the Director of Food, Agriculture and Forest (Prefectoral order 2009–14 from Paris) and under authority of the Office Français de la Biodiversité. For each individual, we recorded the sex and age using tooth eruption pattern  as well as the approximate location of capture. We also verified that individuals of different age and sex are equally distributed throughout the capture period and spatial localization (data not shown).
Parasitological analyses and culture
Feces collected from roe deer per rectum were placed in plastic bags and air removed before sealing. In the field, bags were stored in an insulated box to avoid exposure to low temperatures, which could influence survival and development of some parasite species . Samples were then either directly transported or mailed to the parasitology laboratory of VetAgro Sup (Lyon, France) arriving within 24–48 h of collection. On receipt, we set up coprocultures with 5.4 ± 3.7 g (minimum = 1 g, maximum = 23 g) of feces mixed with vermiculite. After an incubation period of 10–12 days at 24 °C with regular mixing and moisture, the larvae present were collected with a Baermann apparatus after 24 h of sedimentation in water. After determining the total number of L3 collected by counting an aliquot, we stored them in 70% ethanol until DNA extraction.
We also counted eggs of gastrointestinal nematodes following a modified McMaster protocol , with a solution of ZnSO4 (density = 1.36) and a McMaster slide, providing a theoretical sensitivity of 15 eggs per gram (epg) of feces. We also prepared a “control slide” to improve sensitivity of the technique by fully filling a 14-ml tube with the remaining solution and covering it with a coverslip. After centrifugation (5 min at 1200 rpm), the coverslip was transferred on a microscope slide before microscopical observation. We attributed the value 7.5 epg for parasite species detected on the control slide but not on the McMaster slide. Coproscopic data were reported only for samples containing enough L3 collected from coproculture and analyzed by metabarcoding.
DNA extraction and metabarcoding
DNA extraction was performed using the Qiagen DNeasy® powerSoil kit (Qiagen, Hilden, Germany) following the manufacturer’s instructions with an elution volume of 50–150 µl. ITS2 rDNA metabarcoding was performed using the protocol described by Avramenko et al. . Briefly, the rDNA ITS2 target was PCR-amplified from a 1:10 dilution of genomic DNA template. Negative experimental controls were included consisting of water added to the PCR reaction instead of template DNA. These controls were handled in exactly the same way as the actual samples. All samples (n = 149) and negative controls (n = 22) were tagged with unique barcode identifiers to allow the pooling (or normalization) into a single 100 ng amplicon library. The final concentration of the pooled library was assessed with the KAPA qPCR Library Quantification Kit (KAPA Biosystems, USA) following the manufacturer's recommended protocol. The prepared pooled library was run on an Illumina MiSeq Desktop Sequencer using a 600-cycle pair-end reagent kit (MiSeq Reagent Kits v3) at a concentration of 15 nM with the addition 25% PhiX Control v3 (Illumina). A standard demultiplexing protocol was employed that generated FASTQ files for all barcoded samples in the pooled library.
All analyses were carried out using R 3.6 (R Core team, 2020). We curated DNA sequence data in two steps (see Additional file 1: Figure S1). First, the sequence reads were passed through a DADA2-based pipeline . Specific details are available at www.nemabiome.ca. Primers were removed with Cutadapt  after the removal of ambiguous bases. Next, reads were discarded based on the number of errors (> 2 and > 5 in forward and reverse reads, respectively), length (< 50 bp) and quality (truncG = 2 and phiX genome). Forward and reverse reads were merged only if they overlapped by > 12 bp with low mismatches (i.e. maximum 3 for > 100 bp, maximum 2 for > 50 bp and 1 for < 50 bp reads). Putative chimeras were automatically removed but we verified the effect of this filtering on the data set post-analysis, as these programs are known to produce many false positives . A conservative approach was used for taxonomic assignment of amplicon sequence variants (ASV) in which three different assignment methods were used (Additional file 1: Table S1): (1) IDTaxa using the default conservative threshold of 60, which corresponds to the confidence at which to truncate the output taxonomic classification ; (2) assignTaxonomy [46, 50] and (3) BLASTn . The nematode ITS2 rDNA database 1.1.0  containing 9811 complete and non-redundant sequences, with a few additional manual corrections of species assignation corresponding to morphotypes, was used for both IDTaxa and assignTaxonomy. For BLASTn, the three best hits were considered only if they reached a threshold of 90% similarity against NCBI database sequences (2021/04/13). All species-level assignments were ascribed a ‘confidence level' (Additional file 1: Table S1): high confidence indicates that all three methods assigned to the same species and moderate confidence indicates that two out of the three methods agree at the species level and also the third method was in agreement but could only be assigned to the genus level. Taxonomic assignment was limited to the genus level when there was disagreement across any of the three methods at the species level (but consistency at the genus level) or if there was no similar sequence on GenBank with > 98% of similarity and a query cover > 85% with the corresponding ASV.
The raw ASV data set was curated following a procedure developed by Taberlet et al. (2018). ASVs observed only once in the whole data set were considered spurious and were removed . To delete potential contaminants , we discarded any ASVs more abundant in the negative controls (maximum 55 reads) than in other samples and reads with < 0.13% (maximizing the removal of sequences with low abundance in controls) of the total ASV abundance in the entire data set. We only kept ASVs identified as gastrointestinal nematodes at least to the genus level. Finally, any samples with < 1000 reads were removed.
A phylogenetic analysis of ASV was conducted using a maximum likelihood tree with the HKY + G  model according to “ModelTest” function [56, 57]. We used FigTree 1.4.4  to visualize ASV phylogeny. Although the generation of ASVs enables genetic diversity to be characterized to the highest resolution, the biological relevance of much of this diversity remains unknown. A limitation of using ASVs to characterize genetic diversity is that large sample sizes are required to detect significant statistical differences among samples (sensitive to type II error). To overcome this hurdle, we also clustered ASVs into “taxa” and OTUs (operational taxonomic units) using the “IdCluster” function of the package DECIPHER  and a posteriori cutoff of 0.03.
Statistical analyses on measures of parasitic community
We removed the lungworm Varestrongylus from the data set because the focus of the present study was on gastrointestinal nematodes. To control the representativeness of the data, we drew the accumulation curves for the ASV and taxa with the “exact” method of the function “specaccum” of the package vegan . Instead of the exact read count per sample and items, we used the relative frequencies of reads. The read relative frequencies in samples were plotted with ggtree . After controlling for the year, the capture date, the epg and the number of L3 used for DNA extraction, we assessed the effects of age (i.e. young [< 1 year] or adults) and sex (i.e. males or females) on two proxies of diversity of gastrointestinal nematodes community: alpha (measure of community diversity within individual samples) and beta (measure of comparison of dissimilarity between each pair of samples). We measured the beta diversity by using the Bray-Curtis dissimilarity matrices and the alpha diversity by computing three indices: richness (the number of different items), Shannon-Weaver index  and Simpson index (1-D)  indices.
To test for differences in nemabiome among individual hosts of different age, sex and site/sector of Chizé, we tested for the influence of these factors, and their interactions, on the alpha and beta diversity. We ensured that all samples come from a unique individual to avoid individual effect in statistical analyses. Doing so, we removed 19 samples (17 from Chizé and 2 from Trois Fontaines; Additional file 1: Figure S2), selecting for each repeated measure the samples with the higher mass of feces used for the coprocultures. Age and sex effects were first assessed at broad scale with the two sites (Chizé and Trois Fontaines) and, second, at the local scale, comparing the three sectors of Chizé (rich coppices, medium coppices and poor beech stand). We used generalized linear models in a model selection approach to test the effects of site, sex and age of individuals on the alpha diversity, with a Poisson family for richness, a Gaussian family for Shannon and a Gaussian (with logit link because Simpson index ranges from 0 to 1) family for Simpson. For the beta diversity, we used a custom function to compute Akaike’s information criterion corrected for small sample size (AICc) based on residual sums of squares . Following the recommendations of Burnham and Anderson , we identified models that are biologically meaningful and considering limits due to the sample size. It resulted in 18 competitive models. We then ranked the models for each set of candidate models using the AICc and calculated ΔAICc and AICc weights. We selected the model with the lowest AICc value. Models with ΔAICc ≤ 2 were considered equivalent , and in this case, we considered the most parsimonious one, i.e. the model with the lowest degree of freedom.
Differences in nemabiome among samples (beta diversity) were visualized on two dimensions with a non-metric multidimensional scaling (NMDS). We used specifically the NMDS because Bray-Curtis dissimilarities are non-Euclidean embeddable distances . We verified the consistency of patterns between samples resulting from the three taxonomic levels with Procrustes analyses.
ITS2 rDNA nemabiome analyses were performed on 149 samples (n = 36 and 113 samples in Trois Fontaines and Chizé, respectively; Table 1, Fig. 1). The year of sampling had no influence on the nemabiome (Adonis test and Kruskal-Wallis rank sum test; see Additional file 1: Table S2), and this factor was therefore not considered in the following analyses.
The median number eggs per gram of feces was low and heterogeneous (7.5 eggs/g [0; 114]95%IQR) as well as the number of L3 used for DNA extraction (80 [12; 993]95%IQR; Table 1, Fig. 2). Both had no significant effects on alpha and beta diversity (Pearson correlation test and Adonis test; see Additional file 1: Table S2).
Taxonomic and community diversity
Considering the whole data set, we detected 100 ASVs resolved across 14 taxa of gastrointestinal nematodes and 2 taxa of lungworms (Table 1 and Additional file 1: Table S1). ASVs were identified either at the genus (43%) or species (57%) level (Additional file 1: Table S1). Ten genera, including ten species, were identified, and their phylogeny was consistent with previous knowledge on nematode taxonomy (Additional file 1: Figure S3). Sampling captured most of the taxa as the taxa accumulation curves rapidly reached a plateau (Additional file 1: Figure S4). Forty-seven ASVs were found only in one sample for which we were not able to disentangle whether their origin was intra-specific genetic diversity or PCR/sequencing error, but following Callahan et al. , we considered those ASVs in analyses whatever their frequencies in the data set were. Procrustes analyses performed on ASV with taxa have tight correlations (P = 0.001). Thus, we only show results at the ASV taxonomic level, but results from taxa and OTUs are available in Additional file 1: Table S3–S4, Figure S5 and Additional file 2, respectively.
According to the best models on the data subset without repeated measures (Table 2; Additional file 1: Table S3), the richness and Simpson diversity index were higher in Trois Fontaines than in Chizé (estimate ± standard error, richness: β = 0.20 ± 0.06, P = 0.002; Simpson: β = 0.24 ± 0.12, P = 0.041). The Shannon index, which considers the relative abundances of ASV, included significant interactions between age and site, indicating a higher alpha diversity in younger roe deer in Trois Fontaines (β = 0.32 ± 0.13, P = 0.017) (Fig. 3; Table 2). At the local scale of Chizé, no differences among the three sectors were detected for all alpha diversity indexes (Fig. 2).
According to the most parsimonious model selected for the Bray-Curtis dissimilarity, the differences of nemabiome between samples were mostly explained by the site (F1,127 = 35.92, P = 0.001) and age (F1,127 = 12.68, P = 0.001) variables. Locally in Chizé, the selected model for Bray-Curtis dissimilarity depended only on age (F1,94 = 6.36, P = 0.002). Similarly, we observed distinct groups in the NMDS, with partially separated clouds of points for Chizé and Trois Fontaines and a marked separation between young and adult roe deer in Trois Fontaines (Fig. 4).
The most common species or genera in the entire data set were Ostertagia leptospicularis, Trichostrongylus spp. and Spiculopteragia spiculoptera present in 139, 121 and 120 samples of roe deer out of 149 (93%, 81% and 81%, respectively) with a mean relative abundance among those samples of 44%, 29% and 17%, respectively (Fig. 2). In Chizé, O. leptospicularis, Trichostrongylus spp. and S. spiculoptera were present in 98%, 71% and 79% of samples in the sector 1; 93%, 80% and 90% in the sector 2; and 100%, 86% and 100% in the sector 3. The mean relative abundances of species among those samples were 53%, 22% and 17% (sector 1), 46%, 28% and 23% (sector 2) and 45%, 33% and 19% (sector 3). In Trois Fontaines, O. leptospicularis, Trichostrongylus spp. and S. spiculoptera were present in 81%, 97% and 64% of samples with a mean relative abundance among those samples of 24%, 42% and 9%. The ovine species Chabertia ovina, Haemonchus contortus and Bunostomum sp. were only detected in Trois Fontaines (Figs. 2, 5), explaining the significantly higher alpha diversity in Trois Fontaines (Table 2) and the significative difference of nemabiome between the two sites (Table 3). We also observed variation in the prevalence and relative abundance of five Trichostrongylus spp. ASVs between the two sites and between the host age groups, with greater differences of abundance between age in Trois Fontaines compared to Chizé (Fig. 5, A11), as indicated by the significant interaction between site and age on Shannon diversity index.
ITS2 rDNA nemabiome metabarcoding provided fine resolution data on the gastrointestinal nematode parasite community composition of two populations of roe deer in France. The gastrointestinal nematode species detected in Trois Fontaines and Chizé correspond largely to species previously detected on necropsies of roe deer in the two study sites [13, 67]. We also detected species not previously described in these two populations [13, 67], but previously found in other roe deer populations (e.g. in Spain  and in the Netherlands ). This study provided new insights into site- and age-specific gastrointestinal nematode diversity because we showed that young roe deer have different nemabiomes compared to adults, with more pronounced effects of age in one site (i.e. Trois Fontaines).
Diversity of gastrointestinal nematodes in roe deer
In accordance with previous studies on helminths of roe deer (e.g. [37, 38]), the main species detected were O. leptospicularis, Trichostrongylus spp. and Spiculopteragia spp. These gastrointestinal nematodes are commonly found in roe deer, confirming our hypothesis H1 that the dominant nematode species were roe deer-specific parasites given that these populations are enclosed in fenced forest areas. However, we found some differences between the two study sites. We have detected generalist species (i.e. gastrointestinal nematodes species found in different host species) at a relatively high frequency in Trois Fontaines but not in Chizé. In addition to C. ovina previously detected in Trois Fontaines , we have detected for the first time the presence of the highly pathogenic species H. contortus and B. trigonocephalum in this population. In Chizé and Trois Fontaines, we noted the presence of T. axei and O. venulosum, at low relative abundance, also considered as generalist species, and T. colubriformis (only in Chizé) and T. circumcincta, which generally infect sheep and goats [30, 68]. The identifications of H. contortus and T. axei are at a “high confidence level" because these species were well represented in the nemabiome database and the three different methods of assignment provided the same species assignment. For Bunostomum sp., ASVs were only assigned at the genus level as a result of a relatively low level of identity (92% over 283 bp [query cover of 100%, e-value = 9e−110]) to the best BLAST hit in GenBank (Additional file 1: Table S1). To our knowledge, B. trigonocephalum of sheep and B. phlebotomum of calves are the only Bunostomum species which have been identified in roe deer in Europe (e.g. ). These two species are highly divergent, and considering the low level of intra-specific variation observed within B. trigonocephalum , it is possible that the Bunostomum sp. identified in Trois Fontaines is a previously unidentified, or cryptic, species closely related to B. trigonocephalum.
For some ASVs, the taxonomic identification did not reach the species level, e.g. Trichostrongylus spp. (67% among this genus), Spiculopteragia spp. (6%) or Ostertagia spp. (67%). This may be largely due to the lack of reference sequence available for these species in the nemabiome database and ultimately in GenBank. This issue is accentuated when applying nemabiome metabarcoding to wildlife hosts compared with domestic livestock where parasite communities are better defined. For example, there is just a single ITS2 reference sequence for Trichostrongylus capricola compared to 299 ITS2 reference sequences for Haemonchus contortus in the nemabiome database (see also Aivelo and Medlar 2018 for the SILVA database). The potential presence of cryptic species and/ or hybrids of closely related species  may further complicate accurate species assignation. Consequently, the generation of reference sequences from morphologically validated parasites should continue to be a research priority for the scientific community as a whole.
Ashworthius sidemi was not detected in Trois Fontaines, although it has previously been isolated on roe deer carcasses during winter in this study area [13, 71]. This may be due to the hypobiosis in this species  or to the fact that during winter females contained only non-embryonated eggs . Further investigations on the seasonal nemabiome will be required to determine the complete gastrointestinal nematode species richness and variations during the year and species.
Effects of site, age and sex class on the nemabiome
In Trois Fontaines, the age of individuals had a major influence on the gastrointestinal nematode species community found in roe deer, with a higher alpha diversity (β = 0.32 ± 0.13, P = 0.017) observed in young compared to adult animals in accordance with our prediction (H2). This age effect on beta diversity is mostly driven by some specific variants of Trichostrongylus spp. in both sites and of H. contortus in Trois Fontaines that are more prevalent, and of higher relative abundance, in young roe deer compared to adults. The limited influence of age on prevalence and relative abundance observed for some of the other parasite species (see also e.g. ) may result from survival mechanisms of these parasite species, such as the dysregulation of the immune response of the host , from an already developed and efficient immune response of the host when young roe deer were captured against these parasite species (i.e. at 8–10 months of age in our study) or from an immune response that is inefficient.
The higher prevalence and relative abundance of some gastrointestinal nematode species found in adults in Chizé but not in adults in Trois Fontaines suggest a higher ability to resist to certain parasite species (e.g. some Trichostrongylus variants) in adult roe deer in Trois Fontaines. There is probably a synergistic interaction among the body condition of roe deer, their immune response and the richness of parasites to which they have been exposed . Individuals in lower body condition are more susceptible to parasite infection or/and high levels of parasitism are suspected to negatively impact body condition . Such negative correlation between physiological performance and parasitic burden has previously been observed to be more pronounced in Chizé than in Trois Fontaines . In addition, in Trois Fontaines, the higher body condition of roe deer  and strongyle species diversity may have allowed roe deer to invest in innate immunity , which could prevent the establishment or facilitate the expulsion of strongyles . In contrast, roe deer in Chizé, exposed to a lower species richness of strongyles and being in lower body condition, invest more in adaptive response, which is considered less costly than innate response .
Contrary to our expectation (H3), the alpha diversity of the parasite communities did not differ between males and females in this study. Indeed, the probability of exposure to gastrointestinal species is equal between males and females as they do not spatially segregate , and few sexual differences in their feeding behavior or activity have been observed between males and females [81, 82]. In addition, no sex-specific physiological status occurred during the study period, such as nursing young  or rutting , that could have diverted or limited the allocation of resources to parasite defenses. Sampling during such sex-specific energetically costly periods could be pertinent to consider the mechanisms underlying sex and age differences in gastrointestinal nematodes communities.
Despite the differences in habitat quality, the nemabiome of roe deer was not different among the three sectors of Chizé (H4). This suggests that local ecological variation and the resulting differences in body conditions among roe deer inhabiting these different habitats  do not affect the nemabiome. This may also be explained by the movement of roe deer among sectors that can spread the parasites among the feeding areas [40, 84].
Detection of parasite species generally associated with domestic ungulates
The occurrence in roe deer of Trois Fontaines of parasites commonly detected in small domestic ruminants, including highly pathogenic species (e.g. H. contortus, B. trigonocephalum), confirmed previous studies showing that wildlife can share some parasites found in livestock (e.g. [33, 38]). To our knowledge, only four to six domestic sheep have been bred for several years on the other side of a road surrounding the fenced study site of Trois Fontaines (C. Warnant, personal communication). Even if the fence of the study site is regularly checked for damages, the fence does not totally prevent a few roe deer from leaving and/or entering the reserve (e.g. after the Lothar hurricane in 1999). Such movements may have facilitated the introduction of generalist parasites in the population of Trois Fontaines. The consequences of these pathogenic species on roe deer fitness remain however to be determined.
In the present study we investigated the community of gastrointestinal nematode species in two isolated populations of roe deer using ITS2 rDNA nemabiome metabarcoding. Most gastrointestinal nematodes detected in the two sites were parasites classically found in roe deer. Even though the two populations inhabit fenced areas, we also observed gastrointestinal nematodes usually found in domestic ungulates in one of the populations (Trois Fontaines). Despite a fence restraining movements of animals from and to the study site, migrations of a few individuals may have contributed to the movements of parasites among populations, including generalist parasites from domestic ungulates. This could have important implications for the management of wild populations and the introduction of non-native species to wildlife communities. Therefore, it appears highly relevant to study roe deer and domestic ungulates living in sympatry in anthropogenic landscapes dominated by farmland to better understand the dynamic of gastrointestinal nematodes among wild and domestic populations and the risk of cross-transmission. This work also illustrated the power and limitations of ITS2 nemabiome metabarcoding for wildlife studies. It allows the screening of large numbers of samples for a large number of parasites, with identification to the species level in many cases. However, some parasites could still only be identified to the genus level emphasizing the need to continue to enrich the DNA sequence databases for more comprehensive identification.
Availability of data and materials
The roe deer samples information and the sequence data generated and analyzed during the current study are available in the NCBI Short Read Archive under project accession number PRJNA765944 (SUB9936230).
Internal transcribed spacer 2
Amplified sequence variant
Eggs per gram
Hoberg EP, Kocan AA, Rickard LG. Gastrointestinal strongyles in wild ruminants. In: William SM, Pybus MJ, Kocan AA, editors. Parasitic diseases of wild mammals. Ames: Iowa State University Press; 2001. p. 193–227.
Roeber F, Jex AR, Gasser RB. Impact of gastrointestinal parasitic nematodes of sheep, and the role of advanced molecular tools for exploring epidemiology and drug resistance—an Australian perspective. Parasites Vectors. 2013;6:153.
Craig BH, Pilkington JG, Pemberton JM. Gastrointestinal nematode species burdens and host mortality in a feral sheep population. Parasitology. 2006;133:485–96.
Leivesley JA, Bussière LF, Pemberton JM, Pilkington JG, Wilson K, Hayward AD. Survival costs of reproduction are mediated by parasite infection in wild Soay sheep. Ecol Lett. 2019;22:1203–13.
Kelehear C, Brown GP, Shine R. Influence of lung parasites on the growth rates of free-ranging and captive adult cane toads. Oecologia. 2011;165:585–92.
Tompkins DM, Wilson K. Wildlife disease ecology: from theory to policy. Trends Ecol Evol. 1998;13:476–8.
Gunn A, Irvine RJ. Subclinical parasitism and ruminant foraging strategies—a review. Wildl Soc Bull. 2003;31:117–26.
Thompson RCA, Lymbery AJ, Smith A. Parasites, emerging disease and wildlife conservation. Int J Parasitol. 2010;40:1163–70.
Lochmiller RL, Deerenberg C. Trade-offs in evolutionary immunology: just what is the cost of immunity? Oikos. 2000;88:87–98.
Gilot-Fromont E, Jégo M, Bonenfant C, Gibert P, Rannou B, Klein F, et al. Immune phenotype and body condition in roe deer: individuals with high body condition have different, not stronger immunity. PLoS ONE. 2012;7:e45576.
Beldomenico PM, Begon M. Disease spread, susceptibility and infection intensity: vicious circles? Trends Ecol Evol. 2010;25:21–7.
Tompkins DM, Dunn AM, Smith MJ, Telfer S. Wildlife diseases: from individuals to ecosystems. J Anim Ecol. 2011;80:19–38.
Body G, Ferté H, Gaillard J-M, Delorme D, Klein F, Gilot-Fromont E. Population density and phenotypic attributes influence the level of nematode parasitism in roe deer. Oecologia. 2011;167:635–46.
Aleuy OA, Ruckstuhl K, Hoberg EP, Veitch A, Simmons N, Kutz SJ. Diversity of gastrointestinal helminths in Dall’s sheep and the negative association of the abomasal nematode, Marshallagia marshalli, with fitness indicators. PLoS ONE. 2018;13:e0192825.
Sinclair R, Melville L, Sargison F, Kenyon F, Nussey D, Watt K, et al. Gastrointestinal nematode species diversity in Soay sheep kept in a natural environment without active parasite control. Vet Parasitol. 2016;227:1–7.
Gruver AL, Hudson LL, Sempowski GD. Immunosenescence of ageing. J Pathol. 2007;211:144–56.
Peters A, Delhey K, Nakagawa S, Aulsebrook A, Verhulst S. Immunosenescence in wild animals: meta-analysis and outlook. Ecol Lett. 2019;22:1709–22.
Budischak SA, O’Neal D, Jolles AE, Ezenwa VO. Differential host responses to parasitism shape divergent fitness costs of infection. Funct Ecol. 2018;32:324–33.
Bourgoin G, Portanier E, Poirel M-T, Itty C, Duhayer J, Benabed S, et al. Reproductive females and young mouflon (Ovis gmelini musimon × Ovis sp.) in poor body condition are the main spreaders of gastrointestinal parasites. Parasitology. 2021;148:809–18.
Hayward AD. Causes and consequences of intra- and inter-host heterogeneity in defence against nematodes. Parasite Immunol. 2013;35:362–73.
Klein SL. The effects of hormones on sex differences in infection: from genes to behavior. Neurosci Biobehav Rev. 2000;24:627–38.
Moore SL, Wilson K. Parasites as a viability cost of sexual selection in natural populations of mammals. Science. 2002;297:2015–8.
Markle JG, Fish EN. SeXX matters in immunity. Trends Immunol. 2014;35:97–104.
Froeschke G, Harf R, Sommer S, Matthee S. Effects of precipitation on parasite burden along a natural climatic gradient in southern Africa—implications for possible shifts in infestation patterns due to global changes. Oikos. 2010;119:1029–39.
Cheynel L, Lemaître J-F, Gaillard J-M, Rey B, Bourgoin G, Ferté H, et al. Immunosenescence patterns differ between populations but not between sexes in a long-lived mammal. Sci Rep. 2017;7:13700.
Davidson RK, Ličina T, Gorini L, Milner JM. Endoparasites in a Norwegian moose (Alces alces) population—Faunal diversity, abundance and body condition. Int J Parasitol Parasites Wildl. 2015;4:29–36.
Oliver-Guimerá A, Martínez-Carrasco C, Tvarijonaviciute A, Ruiz de Ybáñez MR, Martínez-Guijosa J, López-Olvera JR, et al. The physiological cost of male-biased parasitism in a nearly monomorphic mammal. Parasites Vectors. 2017;10:200.
Dobson A, Lafferty KD, Kuris AM, Hechinger RF, Jetz W. Homage to Linnaeus: how many parasites? How many hosts? Proc Natl Acad Sci USA. 2008;105:11482–9.
Goater TM, Goater CP, Esch GW. Parasitism: the diversity and ecology of animal parasites. Cambridge: Cambridge University Press; 2013.
Taylor MA, Coop RL, Wall RL. Veterinary parasitology. 4th ed. Chichester: Wiley Blackwell; 2015.
Aivelo T, Medlar A. Opportunities and challenges in metabarcoding approaches for helminth community identification in wild mammals. Parasitology. 2018;145:608–21.
Avramenko RW, Redman EM, Lewis R, Yazwinski TA, Wasmuth JD, Gilleard JS. Exploring the gastrointestinal “Nemabiome”: deep amplicon sequencing to quantify the species composition of parasitic nematode communities. PLoS ONE. 2015;10(12):e0143559.
Barone CD, Wit J, Hoberg EP, Gilleard JS, Zarlenga DS. Wild ruminants as reservoirs of domestic livestock gastrointestinal nematodes. Vet Parasitol. 2020;279:109041.
Obanda V, Maingi N, Muchemi G, Ng’ang’a CJ, Angelone S, Archie EA. Infection dynamics of gastrointestinal helminths in sympatric non-human primates, livestock and wild ruminants in Kenya. PLoS ONE. 2019;14:e0217929.
Walker JG, Morgan ER. Generalists at the interface: nematode transmission between wild and domestic ungulates. Int J Parasitol Parasites Wildl. 2014;3:242–50.
Borgsteede FHM, Jansen J, Van Nispen tot Pannerden HPM, Van Der Burg WPJ, Noorman N, Poutsma J, et al. Untersuchungen über die helminthen-fauna beim reh (Capreolus capreolus L.) in den Niederlanden. Z Jagdwiss. 1990;36:104–9.
Pato FJ, Vázquez L, Díez-Baños N, López C, Sánchez-Andrade R, Fernández G, et al. Gastrointestinal nematode infections in roe deer (Capreolus capreolus) from the NW of the Iberian Peninsula: assessment of some risk factors. Vet Parasitol. 2013;196:136–42.
Zaffaroni E, Teresa Manfredi M, Citterio C, Sala M, Piccolo G, Lanfranchi P. Host specificity of abomasal nematodes in free ranging alpine ruminants. Vet Parasitol. 2000;90:221–30.
Segonds-Pichon A, Ferté H, Gaillard J-M, Lamarque F, Duncan P. Nematode infestation and body condition in roe deer (Capreolus capreolus). Game Wildl Sci. 2000;17:241–58.
Pettorelli N, Gaillard J-M, Duncan P, Maillard D, Van Laere G, Delorme D. Age and density modify the effects of habitat quality on survival and movements of roe deer. Ecology. 2003;84:3307–16.
Gaillard J-M, Duncan P, Delorme D, Van Laere G, Pettorelli N, Maillard D, et al. Effects of hurricane Lothar on the population dynamics of European roe deer. J Wildl Manag. 2003;67:767–73.
Gaudry W, Gaillard J-M, Saïd S, Bonenfant C, Mysterud A, Morellet N, et al. Same habitat composition but different use: evidence of context-dependent habitat selection in roe deer females. Sci Rep. 2018;8:5102.
Gaillard J-M, Delorme D, Boutin J-M, Van Laere G, Boisaubert B, Pradel R. Roe deer survival patterns: a comparative analysis of contrasting populations. J Anim Ecol. 1993;62:778–91.
O’Connor LJ, Walkden-Brown SW, Kahn LP. Ecology of the free-living stages of major trichostrongylid parasites of sheep. Vet Parasitol. 2006;142:1–15.
Raynaud J-P, William G, Brunault G. Etude de l’efficacité d’une technique de coproscopie quantitative pour le diagnostic de routine et le contrôle des infestations parasitaires des bovins, ovins, équins et porcins. Ann Parasitol Hum Comp. 1970;45:321–42.
Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJA, Holmes SP. DADA2: High-resolution sample inference from Illumina amplicon data. Nat Methods. 2016;13:581–3.
Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 2011;17:10–2.
Taberlet P, Bonin A, Zinger L, Coissac E. Environmental DNA: for biodiversity research and monitoring. Oxford: Oxford University Press; 2018.
Murali A, Bhargava A, Wright ES. IDTAXA: a novel approach for accurate taxonomic classification of microbiome sequences. Microbiome. 2018;6:140.
Wang Q, Garrity GM, Tiedje JM, Cole JR. Naïve bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol. 2007;73:5261–7.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403–10.
Workentine ML, Chen R, Zhu S, Gavriliuc S, Shaw N, de Rijke J, et al. A database for ITS2 sequences from nematodes. BMC Genet. 2020;21:74.
Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJA, Holmes SP. DADA2: High resolution sample inference from Illumina amplicon data. Nat Methods. 2016;13:581–3.
Schnell IB, Bohmann K, Gilbert MTP. Tag jumps illuminated—reducing sequence-to-sample misidentifications in metabarcoding studies. Mol Ecol Resour. 2015;15:1289–303.
Hasegawa M, Kishino H, Yano T. Dating of the human-ape splitting by a molecular clock of mitochondrial DNA. J Mol Evol. 1985;22:160–74.
Posada D, Crandall KA. modeltest: testing the model of DNA substitution. Bioinformatics. 1998;14:817–8.
Schliep KP. phangorn: phylogenetic analysis in R. Bioinformatics Oxford Academic. 2011;27:592–3.
Rambaut A. (2014). FigTree v1.4.4. http://tree.bio.ed.ac.uk/software/figtree/. Accessed 4 Dec 2020.
Wright ES. Using DECIPHER v2.0 to analyze big biological sequence data in R. R J. 2016;8:352–9.
Oksanen J, Blanchet FG, Friendly M, Roeland K, Legendre P, McGlinn D, et al. (2020). vegan: Community Ecology Package. R package v2.5–7. https://CRAN.R-project.org/package=vegan.
Yu G, Smith DK, Zhu H, Guan Y, Lam TTY. ggtree: an R package for visualization and annotation of phylogenetic trees with their covariates and other associated data. Methods Ecol Evol. 2017;8:28–36.
Shannon CE, Weaver W. A mathematical theory of communication. Urbana: University of Illinois Press; 1948.
Simpson EH. Measurement of diversity. Nature. 1949;163:688.
Dyson K. Custom community ecology helper R scripts [Internet]. 2018. https://github.com/kdyson/R_Scripts
Burnham KP, Anderson DR. Model selection and multimodel inference: a practical information-theoretic approach. 2nd ed. New York: Springer; 2002.
Legendre P, Legendre L. Numerical ecology. 2nd ed. Amsterdam: Elsevier; 1998.
Lienard E, Depaquit J, Ferté H. Spiculopteragia mathevossiani Ruchliadev, 1948 is the minor morph of Spiculopteragia spiculoptera (Gushanskaya, 1931): molecular evidence. Vet Res. 2006;37:683–94.
Wyrobisz-Papiewska A, Kowal J, Nosal P, Chovancová G, Rehbein S. Host specificity and species diversity of the Ostertagiinae Lopez-Neyra, 1947 in ruminants: a European perspective. Parasites Vectors. 2018;11:369.
Wang C, Gao J-F, Chang QC, Zou FC, Zhao Q, Zhu X-Q. Sequence variability in four mitochondrial genes among Bunostomum trigonocephalum isolates from four provinces in China. J Helminthol. 2013;87:416–21.
Wyrobisz A, Kowal J, Nosal P. Insight into species diversity of the Trichostrongylidae Leiper, 1912 (Nematoda: Strongylida) in ruminants. J Helminthol. 2016;90:639–46.
Lehrter V, Jouet D, Liénard E, Decors A, Patrelle C. Ashworthius sidemi Schulz, 1933 and Haemonchus contortus (Rudolphi, 1803) in cervids in France: integrative approach for species identification. Infect Genet Evol. 2016;46:94–101.
Dróżdż J, Demiaszkiewicz A, Lachowicz J. Expansion of the Asiatic parasite Ashworthius sidemi (Nematoda, Trichostrongylidae) in wild ruminants in Polish territory. Parasitol Res. 2003;89:94–7.
Vadlejch J, Kyriánová I, Rylková K, Zikmund M, Langrova I. Health risks associated with wild animal translocation: a case of the European bison and an alien parasite. Biol Invasions. 2017;19:1121–5.
Santín-Durán M, Alunda JM, Hoberg EP, de la Fuente C. Age distribution and seasonal dynamics of abomasal helminths in wild red deer from central Spain. J Parasitol. 2008;94:1031–7.
Else KJ. Have gastrointestinal nematodes outwitted the immune system? Parasite Immunol. 2005;27:407–15.
Jégo M, Ferté H, Gaillard J-M, Klein F, Crespin L, Gilot-Fromont E, et al. A comparison of the physiological status in parasitized roe deer (Capreolus capreolus) from two different populations. Vet Parasitol. 2014;205:717–20.
Bordes F, Morand S. Coevolution between multiple helminth infestations and basal immune investment in mammals: cumulative effects of polyparasitism? Parasitol Res. 2009;106:33–7.
Inclan-Rico JM, Siracusa MC. First responders: innate immunity to helminths. Trends Parasitol. 2018;34:861–80.
Martin LB, Weil ZM, Nelson RJ. Seasonal changes in vertebrate immune activity: mediation by physiological trade-offs. Philos Trans R Soc Lond B Biol Sci. 2008;363:321–39.
Bonenfant C, Gaillard J-M, Dray S, Loison A, Royer M, Chessel D. Testing sexual segregation and aggregation: old ways are best. Ecology. 2007;88:3202–8.
José CS, Lovari S, Ferrari N. Temporal evolution of vigilance in roe deer. Behav Processes. 1996;38:155–9.
Merceron G, Viriot L, Blondel C. Tooth microwear pattern in roe deer (Capreolus capreolus L.) from Chizé (Western France) and relation to food composition. Small Rumin Res. 2004;53:125–32.
Habig B, Doellman MM, Woods K, Olansen J, Archie EA. Social status and parasitism in male and female vertebrates: a meta-analysis. Sci Rep. 2018;8:3629.
Portanier E, Garel M, Devillard S, Duhayer J, Poirel M-T, Henri H, et al. Does host socio-spatial behavior lead to a fine-scale spatial genetic structure in its associated parasites? Parasite. 2019;26:64.
The authors warmly thank all the professionals from the Office Français de la Biodiversité and all the trainees for data collection. This project was founded by Results Driven Agriculture Research (RDAR; grant no. 2017R021R), University of Calgary, the Office Français de la Biodiversité, the NSERC-CREATE (Canada) Host Parasite Interactions (HPI) graduate training program and VetAgro Sup-Pôle d’Expertise Vétérinaire et Agronomique des Animaux Sauvages (EVAAS ; http://evaas.vetagro-sup.fr/; DGAL—VetAgro Sup -INRAE funding). G Bourgoin was supported by the AgreenSkills+ fellowship program (EU program; Marie-Curie FP7 COFUND People Programme; grant agreement n_609398).
See Acknowledgments section.
Ethics approval and consent to participate
Protocol approved by the Director of Food, Agriculture and Forest (Prefectoral order 2009–14 from Paris).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Flowchart summarizing the data curation from Illumina amplicon fastq files to the final curating dataset. Table S1. Taxonomic assignment of ASVs using a combination of three separate methods. Figure S2. Read relative frequencies of GIN at taxa level in samples of 19 recaptured roe deer. Table S2. Results of statistical tests of the relation between diversity index and year of samples, number of L3 and epg. Figure S3. Unrooted maximum likelihood tree of ASVs in samples. Figure S4. Sample accumulation curves. Table S3. Set of generalized linear models and perMANOVA models sorted by AICc value. Table S4. Generalized linear models and perMANOVA models selected for taxa. Figure S5. Non-metric multidimensional scaling (Taxa) of nemabiome.
About this article
Cite this article
Beaumelle, C., Redman, E.M., de Rijke, J. et al. Metabarcoding in two isolated populations of wild roe deer (Capreolus capreolus) reveals variation in gastrointestinal nematode community composition between regions and among age classes. Parasites Vectors 14, 594 (2021). https://doi.org/10.1186/s13071-021-05087-5