Population genetic analysis of the ovine parasitic nematode Teladorsagia circumcincta and evidence for a cryptic species

An understanding of genetic variation in parasite populations, and how it is partitioned, is required to underpin many areas of basic and applied research. Population genetic studies on parasitic nematode populations are still in their infancy and have been dominated by the use of single locus markers. We have used a panel of five microsatellite markers to undertake a genetic study of a number of field and laboratory populations of Teladorsagia circumcincta. High levels of polymorphism were seen in all the populations examined with the majority of diversity being within rather than between populations. There was no detectable genetic differentiation between the UK populations examined although they included both laboratory passaged and field isolates derived from different geographical regions and host species. This broadly supports previous mtDNA sequence diversity studies of this parasite in the UK and USA. However, some between-population genetic differentiation was apparent when several populations from French goats and a laboratory population from New Zealand were examined. Most notably, a population from a French goat farm, which has previously been suggested to contain a cryptic species, showed very high levels of genetic differentiation from all the other populations. Analysis of multi-locus genotypes suggested the presence of two sympatric parasite populations on this farm with little or no gene flow between them. This supports the hypothesis that parasites currently defined as T. circumcincta by routine morphological criteria comprise more than a single species.     

Structure, biodiversity, and historical biogeography of nematode faunas in holarctic ruminants: morphological and molecular diagnoses for Teladorsagia boreoarcticus n. sp. (Nematoda: Ostertagiinae), a dimorphic cryptic species in muskoxen (Ovibos moschatus)

Discovery of the ostertagiine nematode Teladorsagia boreoarcticus n. sp. in muskoxen, Ovibos moschatus, from the central Canadian Arctic highlights the paucity of knowledge about the genealogical and numerical diversity of nematode faunas characteristic of artiodactyls at high latitudes across the Holarctic. Teladorsagia boreoarcticus is a dimorphic cryptic species distinguished from Teladorsagia circumcincta/Teladorsagia trifurcata in domestic sheep by a 13% divergence in the ND4 region of mitochondrial DNA, constant differences in the synlophe, and significantly longer esophageal valve, spicules, gubernaculum, and bursa. Teladorsagia boreoarcticus represents an archaic component of the North American fauna and may have a Holarctic distribution in muskoxen and caribou. Recognition of T. boreoarcticus in muskoxen, in part, corroborates hypotheses for the existence of a cryptic species complex of Teladorsagia spp. among Caprinae and Cervidae at high latitudes and indicates the importance of climatological determinants during the late Tertiary and Pleistocene on diversification of the fauna. Also reinforced is the concept of the North American fauna as a mosaic of endemic and introduced species. Discovery of a previously unrecognized species of Teladorsagia has additional implications and clearly indicates that (1) our knowledge is incomplete relative to potentially pathogenic nematodes that could be exchanged among domestic and wild caprines; (2) we do not have sufficient knowledge of the fauna to understand the ecological control mechanisms (limitations) on dissemination and host range; and (3) an understanding of historical and geographical influences on the genealogical diversity and distribution of nematode faunas in domestic and wild ruminants is requisite to define the interface between agricultural and natural ecosystems across the Holarctic.     

Host Movement and the Genetic Structure of Populations of Parasitic Nematodes

Mitochondrial DNA (mtDNA) sequence data were used to compare the population genetic structures of five species of parasitic nematodes from three different hosts: Ostertagia ostertagi and Haemonchus placei from cattle, H. contortus and Teladorsagia circumcincta from sheep, and Mazamastrongylus odocoilei from white-tailed deer. The parasites of sheep and cattle showed a pattern consistent with high gene flow among populations. The parasite of deer showed a pattern of substantial population subdivision and isolation by distance. It appears that host movement is an important determinant of population genetic structure in these nematodes. High gene flow in the parasites of livestock also indicates great opportunity for the spread of rare alleles that confer resistance to anthelmintic drugs. All species, including the parasite of deer, had unusually high within-population diversities (averages of 0.019-0.027 substitutions per site between pairs of individuals from the same population). Large effective population sizes (Ne), perhaps in combination with rapid mtDNA evolution, appear to be the most likely explanation for these high within-population diversities.     

Generalist nematodes dominate the nemabiome of roe deer in sympatry with sheep at a regional level

The growth of livestock farming and the recent expansion of wild ungulate populations in Europe favor opportunities for direct and/or indirect cross-transmission of pathogens. Comparatively few studies have investigated the epidemiology of gastro-intestinal nematode parasites, an ubiquitous and important community of parasites of ungulates, at the wildlife/livestock interface. In this study, we aimed to assess the influence of livestock proximity on the gastrointestinal nematode community of roe deer in a rural landscape located in southern France. Using ITS-2 rDNA nemabiome metabarcoding on fecal larvae, we analysed the gastrointestinal nematode communities of roe deer and sheep. In addition, we investigated Haemonchus contortus nad4 mtDNA diversity to specifically test parasite circulation among domestic and wild host populations. The dominant gastrointestinal nematode species found in both the roe deer and sheep were generalist species commonly found in small ruminant livestock (e.g. H. contortus), whereas the more specialised wild cervid nematode species (e.g. Ostertagia leptospicularis) were only present at low frequencies. This is in marked contrast with previous studies that found the nemabiomes of wild cervid populations to be dominated by cervid specialist nematode species. In addition, the lack of genetic structure of the nad4 mtDNA of H. contortus populations between host species suggests circulation of gastrointestinal nematodes between roe deer and sheep. The risk of contact with livestock only has a small influence on the nemabiome of roe deer, suggesting the parasite population of roe deer has been displaced by generalist livestock parasites due to many decades of sheep farming, not only for deer grazing close to pastures, but also at a larger regional scale. We also observed some seasonal variation in the nemabiome composition of roe deer. Overall, our results demonstrate significant exchange of gastrointestinal nematodes between domestic and wild ungulates, with generalist species spilling over from domestic ungulates dominating wild cervid parasite communities.     

"Emerging" parasitic infections in arctic ungulates

Important drivers for emergence of infectious disease in wildlifeinclude changes in the environment, shrinking habitats or concentrationof wildlife, and movement of people, animals, pathogens, orvectors. In this paper we present three case-studies of emergingparasitic infections and diseases in ungulates in the Canadiannorth. First we discuss climate warming as an important driverfor the emergence of disease associated with Umingmakstrongyluspallikuukensis, a nematode lungworm of muskoxen. Then we examinehow Protostrongylus stilesi, the sheep lungworm, emerged (orre-emerged) in muskoxen after re-introduction of this host intoits historical range made it sympatric with Dall's sheep. Finally,we consider Teladorsagia boreoarcticus, a newly described andcommon abomasal nematode of muskoxen that is emerging as a disease-causingparasite and may be an important regulator for muskox populationson Banks Island, Northwest Territories. These and other arctichost-parasite systems are exquisitely tuned and constrainedby a harsh and highly seasonal environment. The dynamics ofthese systems will be impacted by climate change and other ecologicaldisruptions. Baseline knowledge of parasite biodiversity andparasite and host ecology, together with predictive models andlong-term monitoring programs, are essential for anticipatingand detecting altered patterns of host range, geographic distribution,and the emergence of parasitic infections and diseases.     

PARASITE-MEDIATED SELECTION AGAINST INBRED SOAY SHEEP IN A FREE-LIVING ISLAND POPULATON

Parasites are thought to provide a selective force capable of promoting genetic variation in natural populations. One rarely considered pathway for this action is via parasite-mediated selection against inbreeding. If parasites impose a fitness cost on their host and the offspring of close relatives have greater susceptibility to parasites due to the increased homozygosity that results from inbreeding, then parasite-mediated mortality may select against inbred individuals. This hypothesis has not yet been tested within a natural vertebrate population. Here we show that relatively inbred Soay sheep (Ovis aries), as assessed by microsatellite heterozygosity, are more susceptible to parasitism by gastrointestinal nematodes, with interactions indicating greatest susceptibility among adult sheep at high population density. During periods of high overwinter mortality on the island of Hirta, St. Kilda, Scotland, highly parasitised individuals were less likely to survive. More inbred individuals were also less likely to survive, which is due to their increased susceptibility to parasitism, because survival was random with respect to inbreeding among sheep that were experimentally cleared of their gastrointestinal parasite burden by anthelminthic treatment. As a consequence of this selection, average microsatellite heterozygosity increases with age in St. Kildan Soay sheep. We suggest that parasite-mediated selection acts to maintain genetic variation in this small island population by removing less heterozygous individuals.     

Modelling non-additive and nonlinear signals from climatic noise in ecological time series: Soay sheep as an example

Understanding how climate can interact with other factors in determining patterns of species abundance is a persistent challenge in ecology. Recent research has suggested that the dynamics exhibited by some populations may be a non-additive function of climate, with climate affecting population growth more strongly at high density than at low density. However, we lack methodologies to adequately explain patterns in population growth generated as a result of interactions between intrinsic factors and extrinsic climatic variation in non-linear systems. We present a novel method (the Functional Coefficient Threshold Auto-Regressive (FCTAR) method) that can identify interacting influences of climate and density on population dynamics from time-series data. We demonstrate its use on count data on the size of the Soay sheep population, which is known to exhibit dynamics generated by nonlinear and non-additive interactions between density and climate, living on Hirta in the St Kilda archipelago. The FCTAR method suggests that climate fluctuations can drive the Soay sheep population between different dynamical regimes--from stable population size through limit cycles and non-periodic fluctuations.     

Genetic diversity and differentiation of 12 eastern Adriatic and western Dinaric native sheep breeds using microsatellites

Nuclear genetic diversity and differentiation of 341 sheep belonging to 12 sheep breeds from Croatia and Bosnia and Herzegovina were examined. The aim of the study was to provide the understanding of the genetic structure and variability of the analysed pramenka sheep populations, and to give indications for conservation strategies based on the population diversity and structure information. The genetic variation of the sheep populations, examined at the nuclear level using 27 microsatellite loci, revealed considerable levels of genetic diversity, similar to the diversity found in other European indigenous low-production sheep breeds. Population-specific alleles were detected at most loci and in breeds analysed. The observed heterozygosity ranged from 0.643 (in Lika pramenka) to 0.743 (in Vlasic pramenka), and the expected heterozygosity ranged from 0.646 (in Lika pramenka) to 0.756 (in Dalmatian pramenka). Significant inbreeding coefficients were found for half of the populations studied and ranged from 0.040 (Pag island sheep) to 0.091 (Kupres pramenka). Moderate genetic differentiation was found between the studied sheep populations. The total genetic variability observed between different populations was 5.29%, whereas 94.71% of the variation was found within populations. Cres island sheep, Lika pramenka and Istrian sheep were identified as the most distinct populations, which was confirmed by the factorial analysis of correspondence and supported through a bootstrapping adjustment to correct for the difference in the sample sizes. The population structure analysis distinguished 12 clusters for the 12 sheep breeds analysed. However, the cluster differentiation was low for Dalmatian, Vlasic, Stolac and Krk pramenka. This systematic study identified Lika pramenka and Rab island sheep as those with the lowest diversity, whereas Istrian sheep and Pag island sheep had the highest. Conservation actions are proposed for Istrian, Rab and Cres island sheep, Lika and Kupres pramenka because of high estimated coefficients of inbreeding.     

Mitochondrial DNA variation in benzimidazole-resistant and -susceptible populations of the small ruminant parasite Teladorsagia circumcincta

The genetic diversity of the mtDNA ND4 gene in 11 Teladorsagia circumcincta populations from France and Morocco was assessed by sequencing. Some of these nematode populations were resistant to benzimidazole (BZ) anthelmintics, while others were susceptible. The nucleotide diversity in all populations studied was very high, probably due to a high mutation rate in nematodes, but there was no significant difference between them. This suggests that no strong, recurrent bottlenecks occur during the acquisition of BZ resistance by a worm population. The conservation of genetic variations during the acquisition of BZ resistance is probably due to the fact that anthelmintic treatments do not kill all the susceptible adult worms and to the presence of numerous free-living larvae that are not submitted to this anthelmintic pressure. There was no genetic subdivision between worm populations on a small geographical scale (less than 200 km), but significant F(ST)s were found on a larger geographical scale. This kind of subdivision cannot be explained by different genetic flows between populations because all these populations were isolated from each other. This subdivision is probably due to the breeding management practices and the large size of these worm populations, which limit genetic drift.     

Benzimidazole resistance allele haplotype diversity in United Kingdom isolates of Teladorsagia circumcincta supports a hypothesis of multiple origins of resistance by recurrent mutation

Polymorphisms in the isotype I β-tubulin gene are important genetic determinants of benzimidazole (BZ) resistance in a number of parasitic nematode species including Teladorsagia circumcincta, a major gastrointestinal nematode of sheep. This study investigates the genetic diversity at this locus in a BZ-resistant isolate of T. circumcincta (MTci5) derived from a sheep farm in the United Kingdom (UK) that was open to animal, and therefore parasite, migration. Pyrosequencing was used to determine the frequency of single nucleotide polymorphisms (SNPs) known to be associated with BZ resistance. This was followed by a combination of single strand conformation polymorphism (SSCP) analysis and nucleotide sequencing to sample allelic diversity in a 276 bp fragment immediately surrounding the isotype I β-tubulin F200Y mutation. The genetic diversity at this locus was extremely high, with seven different haplotypes found to contain the resistant F200Y polymorphism in this single resistant isolate. Genotyping by SSCP interfaced with pyrosequencing demonstrated that the P200^Y mutation is also present on multiple haplotypes in two other BZ-resistant T. circumcincta isolates from the UK. This contrasts with much lower levels of haplotype diversity in BZ-resistant alleles present in T. circumcincta isolates from French goat farms that are closed to parasite migration. Taken together with our knowledge of T. circumcincta population genetic structure, these results are most consistent with multiple independent origins of resistance and mixing of alleles due to the large amount of livestock movement in the UK.     

Unfolding of population structure in Baltic sheep breeds using microsatellite analysis

Studies of domestic animals are performed on breeds, but a breed does not necessarily equate to a genetically defined population. The division of sheep from three native and four modern Baltic sheep breeds was studied using 21 microsatellite loci and applying a Bayesian clustering method. A traditional breed-wise approach was compared to that relying on the pattern of molecular diversity. In this study, a breed was found to be inconsistent with a distinct genetic population for three reasons: (i) a lack of differentiation between modern Baltic breeds, since the majority of the studied sheep formed a single population; (ii) the presence of individuals of foreign ancestry within the breed; and (iii) an undefined local Saaremaa sheep was referred to as a breed, but was shown to consist of separate populations. In the breed-wise approach, only the clearly distinct Ruhnu sheep demonstrated low within-breed variation, although the newly identified Saaremaa populations also have low variability. Providing adequate management recommendations for the Saaremaa sheep is not possible without further studies, but the potential harmful effects of inbreeding in the Ruhnu sheep could be reduced through the use of two genetically related Saaremaa populations. In other breeds, excessive crossing appears to be a larger concern than inbreeding. Assigning individuals into populations based on the pattern of genetic diversity offers potentially unbiased means of elucidating the genetic population structure of species. Combining these genetic populations with phenotypic and aetiological data will enable formulation of the most informed recommendations for gene resource management.     

Density-dependent selection on horn phenotype in Soay sheep.

In the naturally regulated population of Soay sheep (Ovis aries L.) on Hirta (St Kilda), 12% of males and 24% of females have scurred horns. This trait reduces the competitive ability of males in the mating season but is associated with higher overwinter survival rates in both sexes (Moorcroft et al. 1996). In this paper, we show that scurred females also show higher conception rates and weaning rates than non-scurred ones, and that these differences are associated with heavier maternal body weight as well as higher birth weights of offspring. Selection pressures favouring scurredness in females vary with population density and are generally more pronounced among younger animals than adults. We discuss these results with reference to recent suggestions that temporal fluctuation in selection pressures may help to maintain genetic diversity within populations. We suggest that selection against scurredness operating through male mating success is opposed by selection favouring this trait, operating through the breeding success of females and the survival of both sexes.     

A molecular diagnostic tool to replace larval culture in conventional faecal egg count reduction testing in sheep

The accurate diagnosis of parasitic nematode infections in livestock (including sheep and goats) is central to their effective control and the detection of the anthelmintic resistance. Traditionally, the faecal egg count reduction test (FECRT), combined with the technique of larval culture (LC), has been used widely to assess drug-susceptibility/resistance in strongylid nematodes. However, this approach suffers from a lack of specificity, sensitivity and reliability, and is time-consuming and costly to conduct. Here, we critically assessed a specific PCR assay to support FECRT, in a well-controlled experiment on sheep with naturally acquired strongylid infections known to be resistant to benzimidazoles. We showed that the PCR results were in close agreement with those of total worm count (TWC), but not of LC. Importantly, albendazole resistance detected by PCR-coupled FECRT was unequivocally linked to Teladorsagia circumcincta and, to lesser extent, Trichostrongylus colubriformis, a result that was not achievable by LC. The key findings from this study demonstrate that our PCR-coupled FECRT approach has major merit for supporting anthelmintic resistance in nematode populations. The findings also show clearly that our PCR assay can be used as an alternative to LC, and is more time-efficient and less laborious, which has important practical implications for the effective management and control strongylid nematodes of sheep.     

Functional consequences of genetic diversity in Strongyloides ratti infections

Parasitic nematodes show levels of genetic diversity comparable to other taxa, but the functional consequences of this are not understood. Thus, a large body of theoretical work highlights the potential consequences of parasite genetic diversity for the epidemiology of parasite infections and its possible implications for the evolution of host and parasite populations. However, few relevant empirical data are available from parasites in general and none from parasitic nematodes in particular. Here, we test two hypotheses. First, that different parasitic nematode genotypes vary in life-history traits, such as survivorship and fecundity, which may cause variation in infection dynamics. Second, that different parasitic nematode genotypes interact within the host (either directly or via the host immune system) to increase the mean reproductive output of mixed-genotype infections compared with single-genotype infections. We test these hypotheses in laboratory infections using genetically homogeneous lines of Strongyloides ratti. We find that nematode genotypes do vary in their survivorship and fecundity and, consequently, in their dynamics of infection. However, we find little evidence of interactions between genotypes within hosts under a variety of trickle- and single-infected infection regimes.     

Assessing the relative importance of intrinsic and extrinsic influence on sheep population dynamics on Hirta Island, UK

Identifying the relative importance of intrinsic factors and extrinsic environmental variations on population or ecosystem dynamics is important for ecological conservation research. Here, we use a systematic method proposed by De Menezes and Barabási [2004. Separating internal and external dynamics of complex systems. Physical Review Letters, 93, 068701] to reanalyse the long-term monitoring data of Soay sheep population fluctuations and climate variations on Hirta Island, UK. Our results indicate that the climate conditions have a higher impact than internal factors on the fluctuations of sheep population. The sheep population dynamics are internally self-regulating. The scaling relationships between sheep population and external and internal standard deviations are similar. The threshold of the sheep population on this island as determined by our study is around 1197, which is consistent with previous studies by other methods. Our study indicates that this systematic method may help to understand some of the complicated aspects of population dynamics about which detailed knowledge is limited.     

A 4 year observation of gastrointestinal nematode egg counts,nemabiomes and the benzimidazole resistance genotypes of Teladorsagia circumcincta on a Scottish sheep farm

Anthelmintic resistance threatens the sustainability of sheep production globally. Advice regarding strategies to reduce the development of anthelmintic resistance incorporates the outcomes of modelling exercises. Further understanding of gastrointestinal nematode species diversity, and population dynamics and genetics (which may vary between species) is required to refine these models; and field studies combining faecal egg outputs, species composition and resistance genetics are needed to calibrate them. In this study, faecal samples were taken from ewes and lambs on a commercial farm in south-eastern Scotland at approximately 3 t-4 week intervals between spring and autumn over a period of 4 years. Faecal egg counts were performed on these samples, and L₃ were collected from pooled coprocultures. Deep amplicon sequencing was used to determine both the species composition of these L₃ and the proportions of benzimidazole-resistant single nucleotide polymorphisms in the isotype-1 β-tubulin locus of the predominant species, Teladorsagia circumcincta L₃. Despite consistent management throughout the study, the results show variation in gastrointestinal nematode species composition with time and between age groups, that was potentially associated with weather conditions. The F200Y benzimidazole resistance mutation is close to genetic fixation in the T. circumcincta population on this farm. There was no evidence of variation in isotype-1 β-tubulin single nucleotide polymorphisms frequency between age groups, and no genetic evidence of reversion to benzimidazole susceptibility, despite targeted benzimidazole usage. This study highlights the need to include speciation when investigating gastrointestinal nematode epidemiology and anthelmintic resistance, and serves as an example of how genetic data may be analysed alongside species diversity and faecal egg counts, when markers for other anthelmintic classes are identified.     

Biogeography of Parasitic Nematode Communities in the Galápagos Giant Tortoise: Implications for Conservation Management

The Galápagos giant tortoise is an icon of the unique, endemic biodiversity of Galápagos, but little is known of its parasitic fauna. We assessed the diversity of parasitic nematode communities and their spatial distributions within four wild tortoise populations comprising three species across three Galápagos islands, and consider their implication for Galápagos tortoise conservation programmes. Coprological examinations revealed nematode eggs to be common, with more than 80% of tortoises infected within each wild population. Faecal samples from tortoises within captive breeding centres on Santa Cruz, Isabela and San Cristobal islands also were examined. Five different nematode egg types were identified: oxyuroid, ascarid, trichurid and two types of strongyle. Sequencing of the 18S small-subunit ribosomal RNA gene from adult nematodes passed with faeces identified novel sequences indicative of rhabditid and ascaridid species. In the wild, the composition of nematode communities varied according to tortoise species, which co-varied with island, but nematode diversity and abundance were reduced or altered in captive-reared animals. Evolutionary and ecological factors are likely responsible for the variation in nematode distributions in the wild. This possible species/island-parasite co-evolution has not been considered previously for Galápagos tortoises. We recommend that conservation efforts, such as the current Galápagos tortoise captive breeding/rearing and release programme, be managed with respect to parasite biogeography and host-parasite co-evolutionary processes in addition to the biogeography of the host.     

Population biology and biological control of nematodes.

We studied the population biology of the nematophagous fungus Hirsutella rhossiliensis to understand its potential as a biological control agent. Because the fungus is an infectious and transmissible parasite, we framed our study within an epidemiological context. Field observations, theory, and experiments demonstrated that (i) parasitism of nematodes by H. rhossiliensis is dependent on nematode density, (ii) local populations of the fungus will go extinct unless supplied with some minimum number of nematodes (the host threshold density), and (iii) natural epidemics of this fungus in populations of nematodes develop slowly and only after long periods of high host density. Additional in-depth research on population biology is needed to explain other biological control systems and to guide future research. The most effective research will combine field observation, theory, and experimentation.     

A long-term genetic survey of an ungulate population reveals balancing selection acting on MHC through spatial and temporal fluctuations in selection

We explored a 13-year genetic survey of the major histocompatibility complex (MHC) and neutral loci of the Soay sheep population of St Kilda to test the existence and causes of balancing selection at the MHC. The sheep population experiences demographic fluctuations, partly driven by the nematode Teladorsagia circumcincta. The spatial differentiation detected at the MHC was comparable to that at neutral loci between 1988 and 1996, but significantly lower between 1996 and 2000. The rate of temporal genetic differentiation was higher at the MHC, but within the Eastern heft only. These comparisons of spatial and temporal divergence at MHC and non-MHC loci provide strong evidence of balancing selection at the MHC, acting through spatial and temporal heterogeneity in selection pressure. This heterogeneity could be due to fluctuations in the selection imposed by parasites, either directly, because the prevalence in T. circumcincta varies in space and time, or indirectly, because the fitness costs of infection may vary with resource availability.     

Genetic signatures of a Mediterranean influence in Iberian Peninsula sheep husbandry

Highly adaptable and versatile populations of domestic sheep, the result of millennia of intense husbandry, are found in almost every corner of the world. Here we describe a genetic survey of sheep from the western fringe of its European distribution. We studied the mitochondrial DNA control region sequences from 161 individuals belonging to 7 Portuguese sheep breeds. Our study revealed a high level of genetic diversity, with an average breed haplotype diversity of 0.983, substantially above that observed in central European breeds, as well as the presence of maternal lineages until now only found in the Middle East and Asia. A broad north-south pattern describes the most important trend in the Portuguese sheep population with a southern population clearly distinct from most other breeds. A recurrent influx of new genetic diversity, probably via the Mediterranean Sea, may explain these patterns and appears to corroborate the importance of this maritime route in the history of both mankind and livestock. Zooarchaeological studies of sheep bones from southern Portugal indicate a marked size increase during the Moslem period that may reflect an improvement of this animal--perhaps part of the well known "Arab agricultural revolution" in Andalusia. This could have been a time when the gene pool of Iberian sheep was substantially enriched and may help to explain the history of modern sheep breeds in this peninsula.