Life cycle variation and the genetic structure of nematode populations

Genomic Epidemiology in Filarial Nematodes: Transforming the Basis for Elimination Program Decisions

Onchocerciasis and lymphatic filariasis are targeted for elimination, primarily using mass drug administration at the country and community levels. Elimination of transmission is the onchocerciasis target and global elimination as a public health problem is the end point for lymphatic filariasis. Where program duration, treatment coverage, and compliance are sufficiently high, elimination is achievable for both parasites within defined geographic areas. However, transmission has re-emerged after apparent elimination in some areas, and in others has continued despite years of mass drug treatment. A critical question is whether this re-emergence and/or persistence of transmission is due to persistence of local parasites-i.e., the result of insufficient duration or drug coverage, poor parasite response to the drugs, or inadequate methods of assessment and/or criteria for determining when to stop treatment-or due to re-introduction of parasites via human or vector movement from another endemic area. We review recent genetics-based research exploring these questions in Onchocerca volvulus, the filarial nematode that causes onchocerciasis, and Wuchereria bancrofti, the major pathogen for lymphatic filariasis. We focus in particular on the combination of genomic epidemiology and genome-wide associations to delineate transmission zones and distinguish between local and introduced parasites as the source of resurgence or continuing transmission, and to identify genetic markers associated with parasite response to chemotherapy. Our ultimate goal is to assist elimination efforts by developing easy-to-use tools that incorporate genetic information about transmission and drug response for more effective mass drug distribution, surveillance strategies, and decisions on when to stop interventions to improve sustainability of elimination.

The population genetics of parasitic nematodes of wild animals

Parasitic nematodes are highly diverse and common, infecting virtually all animal species, and the importance of their roles in natural ecosystems is increasingly becoming apparent. How genes flow within and among populations of these parasites - their population genetics - has profound implications for the epidemiology of host infection and disease, and for the response of parasite populations to selection pressures. The population genetics of nematode parasites of wild animals may have consequences for host conservation, or influence the risk of zoonotic disease. Host movement has long been recognised as an important determinant of parasitic nematode population genetic structure, and recent research has also highlighted the importance of nematode life histories, environmental conditions, and other aspects of host ecology. Commonly, factors influencing parasitic nematode population genetics have been studied in isolation, such that an integrated view of the drivers of population genetic structure of parasitic nematodes is still lacking. Here, we seek to provide a comprehensive, broad, and integrative picture of these factors in parasitic nematodes of wild animals that will be a useful resource for investigators studying non-model parasitic nematodes in natural ecosystems. Increasingly, new methods of analysing the population genetics of nematodes are becoming available, and we consider the opportunities that these afford in resolving hitherto inaccessible questions of the population genetics of these important animals.

Population structure analysis of the neglected parasite Thelazia callipaeda revealed high genetic diversity in Eastern Asia isolates

Background

Thelazia callipaeda is the causative agent of thelaziasis in canids, felids and humans. However, the population genetic structure regarding this parasite remains unclear.

Methodology/principal findings

In this study, we first explored the genetic variation of 32 T. callipaeda clinical isolates using the following multi-molecular markers: cox1, cytb, 12S rDNA, ITS1 and 18S rDNA. The isolates were collected from 13 patients from 11 geographical locations in China. Next, the population structure of T. callipaeda from Europe and other Asian countries was analyzed using the cox1 sequences collected during this study and from the GenBank database. In general, the Chinese clinical isolates of T. callipaeda expressed high genetic diversity. Based on the cox1 gene, a total of 21 haplotypes were identified. One only circulated in European countries (Hap1), while the other 20 haplotypes were dispersed in Korea, Japan and China. There were five nucleotide positions in the cox1 sequences that were confirmed as invariable among individuals from Europe and Asia, but the sequences were distinct between these two regions. Population differences between Europe and Asian countries were greater than those among China, Korea and Japan. The T. callipaeda populations from Europe and Asia should be divided into two separate sub-populations. These two groups started to diverge during the middle Pleistocene. Neutrality tests, mismatch distribution and Bayesian skyline plot (BSP) analysis all rejected possible population expansion of T. callipaeda.

Conclusions

The Asian population of T. callipaeda has a high level of genetic diversity, but further studies should be performed to explore the biology, ecology and epidemiology of T. callipaeda.

Thelazia callipaeda is the causative agent of thelaziasis canids, felids and humans. Despite the existing threat of thelaziosis in China, the genetic diversity of T. callipaeda has not been investigated across its wide geographical distribution in China, yet such information may provide insight into the disease epidemiology and the development of specific control measures. In this study, the genetic variation of 32 T. callipaeda clinical isolates collected from 13 patients from 11 geographical locations in China were explored using the following multi-molecular markers: cox1, cytb, 12S rDNA, ITS1 and 18S rDNA. In addition, the population structure of T. callipaeda from Europe and other Asian countries was analyzed using the cox1 sequences collected during this study and from the GenBank database. In general, the Chinese clinical isolates of T. callipaeda demonstrated high genetic diversity. Based on the cox1 gene, a total of 21 haplotypes were identified, one circulated in European countries (Hap1), while the other 20 haplotypes were dispersed in Korea, Japan and China. There were five nucleotide positions in the cox1 that were confirmed as invariable among individuals from Europe and Asia, but the sequences were distinct between these two regions. Population differences between Europe and Asian countries were greater than those among China, Korea and Japan, such that the T. callipaeda population from Europe and Asia should be divided into two separate sub-populations. These two groups started to diverge during the middle Pleistocene.

Belowground Plant-Herbivore Interactions Vary among Climate-Driven Range-Expanding Plant Species with Different Degrees of Novel Chemistry

An increasing number of studies report plant range expansions to higher latitudes and altitudes in response to global warming. However, consequences for interactions with other species in the novel ranges are poorly understood. Here, we examine how range-expanding plant species interact with root-feeding nematodes from the new range. Root-feeding nematodes are ubiquitous belowground herbivores that may impact the structure and composition of natural vegetation. Because of their ecological novelty, we hypothesized that range-expanding plant species will be less suitable hosts for root-feeding nematodes than native congeneric plant species. In greenhouse and lab trials we compared nematode preference and performance of two root-feeding nematode species between range-expanding plant species and their congeneric natives. In order to understand differences in nematode preferences, we compared root volatile profiles of all range-expanders and congeneric natives. Nematode preferences and performances differed substantially among the pairs of range-expanders and natives. The range-expander that had the most unique volatile profile compared to its related native was unattractive and a poor host for nematodes. Other range-expanding plant species that differed less in root chemistry from native congeners, also differed less in nematode attraction and performance. We conclude that the three climate-driven range-expanding plant species studied varied considerably in their chemical novelty compared to their congeneric natives, and therefore affected native root-feeding nematodes in species-specific ways. Our data suggest that through variation in chemical novelty, range-expanding plant species may vary in their impacts on belowground herbivores in the new range.

Genetic variation in mitochondrial cox2 of Heterakis gallinarum from poultry in Sichuan, China

Heterakis gallinarum is one of the common parasitic nematodes found in the caecum of poultry. To investigate the genetic diversity and genetic structure of the H. gallinarum population in Sichuan, we amplified and sequenced the complete mitochondrial (mt) cytochrome c oxidase subunit II (cox2) gene of 59 H. gallinarum isolates from seven different geographical regions, then analyzed their genetic polymorphisms. All cox2 genes of the 59 H. gallinarum isolates were 696 bp in length, with an average A + T content of 67.1%. Fifty-nine sequences contained 34 variable sites, and were classified into 23 haplotypes (HS1-HS23). The values of haplotype diversity (Hd) and nucleotide diversity (π) were 0.688 and 0.00288, respectively. Based on values of F_ST and Nm (F_ST = 0.01929, Nm = 12.71), there was a frequent gene flow but no significant genetic differentiation observed among the populations. The network map showed that the most prominent haplotype was HS1, and the other haplotypes (HS2-HS23) were centered on HS1 with a star-like topology, indicating that H. gallinarum had previously experienced a population expansion. To our knowledge, this is the first research on the population genetics of H. gallinarum based on mitochondrial cox2.

Temporal sampling helps unravel the genetic structure of naturally occurring populations of a phytoparasitic nematode. 2. Separating the relative effects of gene flow and genetic drift

Studying wild pathogen populations in natural ecosystems offers the opportunity to better understand the evolutionary dynamics of biotic diseases in crops and to enhance pest control strategies. We used simulations and genetic markers to investigate the spatial and temporal population genetic structure of wild populations of the beet cyst nematode Heterodera schachtii on a wild host plant species, the sea beet (Beta vulgaris spp. maritima), the wild ancestor of cultivated beets. Our analysis of the variation of eight microsatellite loci across four study sites showed that (i) wild H. schachtii populations displayed fine-scaled genetic structure with no evidence of substantial levels of gene flow beyond the scale of the host plant, and comparisons with simulations indicated that (ii) genetic drift substantially affected the residual signals of isolation-by-distance processes, leading to departures from migration-drift equilibrium. In contrast to what can be suspected for (crop) field populations, this showed that wild cyst nematodes have very low dispersal capabilities and are strongly disconnected from each other. Our results provide some key elements for designing pest control strategies, such as decreasing passive dispersal events to limit the spread of virulence among field nematode populations.

Genetic diversity and population genetics of large lungworms (Dictyocaulus, Nematoda) in wild deer in Hungary

Dictyocaulus nematode worms live as parasites in the lower airways of ungulates and can cause significant disease in both wild and farmed hosts. This study represents the first population genetic analysis of large lungworms in wildlife. Specifically, we quantify genetic variation in Dictyocaulus lungworms from wild deer (red deer, fallow deer and roe deer) in Hungary, based on mitochondrial cytochrome c oxidase subunit 1 (cox1) sequence data, using population genetic and phylogenetic analyses. The studied Dictyocaulus taxa display considerable genetic diversity. At least one cryptic species and a new parasite–host relationship are revealed by our molecular study. Population genetic analyses for Dictyocaulus eckerti revealed high gene flow amongst weakly structured spatial populations that utilise the three host deer species considered here. Our results suggest that D. eckerti is a widespread generalist parasite in ungulates, with a diverse genetic backround and high evolutionary potential. In contrast, evidence of cryptic genetic structure at regional geographic scales was observed for Dictyocaulus capreolus, which infects just one host species, suggesting it is a specialist within the studied area. D. capreolus displayed lower genetic diversity overall, with only moderate gene flow compared to the closely related D. eckerti. We suggest that the differing vagility and dispersal behaviour of hosts are important contributing factors to the population structure of lungworms, and possibly other nematode parasites with single-host life cycles. Our findings are of relevance for the management of lungworms in deer farms and wild deer populations.

Species mtDNA genetic diversity explained by infrapopulation size in a host-symbiont system

Understanding what shapes variation in genetic diversity among species remains a major challenge in evolutionary ecology, and it has been seldom studied in parasites and other host‐symbiont systems. Here, we studied mtDNA variation in a host‐symbiont non‐model system: 418 individual feather mites from 17 feather mite species living on 17 different passerine bird species. We explored how a surrogate of census size, the median infrapopulation size (i.e., the median number of individual parasites per infected host individual), explains mtDNA genetic diversity. Feather mite species genetic diversity was positively correlated with mean infrapopulation size, explaining 34% of the variation. As expected from the biology of feather mites, we found bottleneck signatures for most of the species studied but, in particular, three species presented extremely low mtDNA diversity values given their infrapopulation size. Their star‐like haplotype networks (in contrast with more reticulated networks for the other species) suggested that their low genetic diversity was the consequence of severe bottlenecks or selective sweeps. Our study shows for the first time that mtDNA diversity can be explained by infrapopulation sizes, and suggests that departures from this relationship could be informative of underlying ecological and evolutionary processes.

Insight into species diversity of the Trichostrongylidae Leiper, 1912 (Nematoda: Strongylida) in ruminants

This paper focuses on the species diversity among the Trichostrongylidae Leiper, 1912 (Nematoda: Strongylida), and complexity of the family systematics. Polymorphism (subfamilies: Ostertagiinae, Cooperiinae and Haemonchinae), the presence of cryptic species (genus: Teladorsagia) and hybridization (genera: Cooperia, Haemonchus and Ostertagia) are presented and discussed, considering both morphological and molecular evidence. Some of these phenomena are common, nevertheless not sufficiently understood, which indicates the need for expanding the current state of knowledge thereof. Within the Trichostrongylidae, species distinction supported merely by morphological features is difficult, and requires confirmation by means of molecular methods. The parasitic nematode taxonomy is complicated mainly by the genus Teladorsagia, but complexity may also be expected among other Ostertagiinae (e.g. in the genera Ostertagia and Marshallagia). The data presented here show that the members of the Trichostrongylidae can significantly complicate unambiguous species identification. Hence, it is essential to consider the phenomena mentioned, to gather valid and comparable data on the biodiversity of this family.

Higher relatedness mitigates mortality in a nematode with lethal male fighting

According to kin selection theory, individuals show less aggression towards their relatives. Limited dispersal promotes interactions among relatives but also increases competition among them. The evolution of cooperation in viscous populations has been subject of mainly theoretical exploration. We investigated the influence of relatedness on aggression in males of entomopathogenic nematode Steinernema longicaudum that engage in lethal fighting. In a series of in vitro experiments, we found that both competitor male group size and relatedness influence male mortality rates. Higher relatedness led to progressively lower rates of male mortality. In experimentally infected insects, wherein large numbers of males and females interact, the proportion of dead and paralysed (= terminally injured) males was higher when infection was established by infective juveniles originating from a mixture of three lines than in those infected by a single line. The results collectively show that Steinernema longicaudum males recognize their kin and consequently male mortality rates are lower in groups consisting of more related males. Furthermore, this monotonic negative relationship between aggression and relatedness suggests that kin selection benefits are still substantial even under extreme competition. Our experiments also suggest that kin recognition in entomopathogenic nematodes has a genetic basis rather than being strictly based on environmental cues. We discuss our findings within the theoretical context of the evolution of altruistic/cooperative behaviour in structured populations.

Species discrimination in the subfamily Ostertagiinae of Northern China: assessment of DNA barcode in a taxonomically challenging group

Gastrointestinal nematodes within the subfamily Ostertagiinae (Teladorsagia, Ostertagia, and Marshallagia et al.) are among the most common infections of domesticated livestock. These parasites are of particular interest, as many of the species within this group are of economic importance worldwide. Traditionally, nematode species designations have been based on morphological criteria. However, this group possesses poorly defined species. There is an urgent need to develop a reliable technique that can distinguish species of Ostertagiinae. DNA barcoding has been proved to be a powerful tool to identify species of birds, mammals, and arthropods, but this technique has not yet been examined for identifying species of Ostertagiinae. In this study, a total of 138 mitochondrial DNA (mtDNA) cytochrome c oxidase subunit I (COI) sequences from individuals representing 11 species of Ostertagiinae were acquired by PCR for the first time. The specimens were collected from pastoral area of northern China. Genetic divergence analyses showed that mean interspecific Kimura two-parameter distances of COI (13.61 %) were about four times higher than the mean value of the intraspecific divergence (3.69 %). Then, the performance of the COI to identify species of Ostertagiinae was evaluated by identification success rates using nearest neighbor (NN) and BLASTn. The results indicated that the rates of correct sequence identification for COI were high (>80 %) when using the NN and BLASTn methods. Besides, the deep lineage divergences are detected in Teladorsagia circumcincta. Meanwhile, the analyses also detected no genetic differentiation between some species such as Ostertagia hahurica and Ostertagia buriatica. These results indicate that the traditional status of species within Ostertagiinae should be closely examined based on the molecular data.

Coevolution between Contracaecum (Nematoda, Anisakidae) and Austrolebias (Cyprinodontiformes, Rivulidae) host-parasite complex from SW Atlantic coastal basins

In recent years, molecular studies in host-parasite interactions in terms of coevolution have become important. Larvae (L3) of two species of Contracaecum were found parasitizing species of Rivulidae in the Atlantic coastal basins from Uruguay. The aim of this study is to determine the patterns of differentiation of this host-parasite complex in order to clarify possible coevolutionary events in such interaction throughout phylogeographic approach using both nuclear and mitochondrial molecular markers (internal transcribed spacers (ITS) and cytochrome oxidase subunit 1 (cox-1)). Based on both markers, intraspecific variation in Contracaecum species was lower than 2 %, while interspecific variation was greater than 10 %. Both species of Contracaecum constitute monophyletic groups. Contracaecum resulted in a paraphyletic genus when incorporating other Contracaecum species and closely related nematode sequences from GenBank. ITS regions showed that Contracaecum sp. 1 is more closely related to other species of the same genus than with their counterparts from Atlantic coastal basins in Uruguay. Haplotype network for both markers corroborate the existence of two distinct taxa. While ITS pairwise FST comparisons and the indirect estimate of gene flow confirm the existence of two distinct Contracaecum species, mitochondrial gene detected low levels of migrants between some of the populations from both species. Our results suggest that coevolution in this host-parasite complex species is plausible. Parasite cladogenetic events occur almost simultaneously with the separation of the hypothetical ancestors of each species complex of Austrolebias during Pliocene. Additionally, the two lineages of Contracaecum colonize differently the species within each of the Austrolebias complexes.

Comparative Estimation of Genetic Diversity in Population Studies using Molecular Sampling and Traditional Sampling Methods

Entomopathogenic nematodes (EPN) are efficient biological pest control agents. Population genetics studies on EPN are seldom known. Therefore, it is of interest to evaluate the significance of molecular sampling method (MSM) for accuracy, time needed, and cost effectiveness over traditional sampling method (TSM). The study was conducted at the Mohican Hills golf course at the state of Ohio where the EPN H. bacteriophora has been monitored for 18 years. The nematode population occupies an area of approximately 3700 m(2) with density range from 0.25-2 per gram soil. Genetic diversity of EPN was studied by molecular sampling method (MSM) and traditional sampling method (TSM) using the mitochondrial gene pcox1. The MSM picked 88% in compared to TSM with only 30% of sequenced cox 1 gene. All studied genetic polymorphism measures (sequence and haplotype) showed high levels of genetic diversity of MSM over TSM. MSM minimizes the chance of mitochondrial genes amplification from non target organisms (insect or other contaminating microorganisms). Moreover, it allows the sampling of more individuals with a reliable and credible representative sample size. Thus, we show that MSM supersedes TSM in labour intensity, time consumption and requirement of no special experience and efficiency.

Population genetic structure of the acanthocephalan Acanthosentis cheni in anadromous, freshwater, and landlocked stocks of its fish host, Coilia nasus

The acanthocephalan Acanthosentis cheni was found in anadromous, freshwater, and landlocked stocks of its fish host, Coilia nasus. To examine the genetic variations of the acanthocephalan among the 3 populations with the adaptation of the host to the freshwater, the genetic structure of the helminth was investigated in anadromous (Zhoushan and Chongming islands, and Anqing), freshwater (Anqing, Ezhou, and Poyang Lake), and landlocked (Tian'ezhou Reserve) populations by sequencing intergenic transcribed spacers (ITS) of the ribosomal RNA coding genes. Low Fst values and high gene flow were found among the 7 populations (Fst = 0.0135, P = 0.2723; Nm = 36.48) and the 3 ecotypes of Acanthosentis cheni (Fst = 0.0178, P = 0.1044; Nm = 27.67). On the other hand, significant genetic differentiation of the C. nasus host populations was detected between the upstream and downstream areas of Xiaogu Mountain (Fst = 0.1961, P = 0.0030; Nm = 2.05), which is the farthest location of spawning migration for C. nasus . However, the migration break of the fish host appeared not to cause significant genetic differentiation of A. cheni populations between the upper and lower reaches of Xiaogu Mountain. Other factors might promote genetic exchange of A. cheni populations such as dispersal of the intermediate host by flooding or other fish species serving as the definitive or paratenic hosts. In Anqing, nucleotide diversity of the acanthocephalan was highest in the freshwater population (0.0038) and lower in the anadromous population (0.0026). This suggested that new mutations may have occurred in the freshwater A. cheni population in Anqing when adapting to a freshwater environment.

Reductions in genetic diversity of Schistosoma mansoni populations under chemotherapeutic pressure: the effect of sampling approach and parasite population definition

Detecting potential changes in genetic diversity in schistosome populations following chemotherapy with praziquantel (PZQ) is crucial if we are to fully understand the impact of such chemotherapy with respect to the potential emergence of resistance and/or other evolutionary outcomes of interventions. Doing so by implementing effective, and cost-efficient sampling protocols will help to optimise time and financial resources, particularly relevant to a disease such as schistosomiasis currently reliant on a single available drug. Here we explore the effect on measures of parasite genetic diversity of applying various field sampling approaches, both in terms of the number of (human) hosts sampled and the number of transmission stages (miracidia) sampled per host for a Schistosoma mansoni population in Tanzania pre- and post-treatment with PZQ. In addition, we explore population structuring within and between hosts by comparing the estimates of genetic diversity obtained assuming a 'component population' approach with those using an 'infrapopulation' approach. We found that increasing the number of hosts sampled, rather than the number of miracidia per host, gives more robust estimates of genetic diversity. We also found statistically significant population structuring (using Wright's F-statistics) and significant differences in the measures of genetic diversity depending on the parasite population definition. The relative advantages, disadvantages and, hence, subsequent reliability of these metrics for parasites with complex life-cycles are discussed, both for the specific epidemiological and ecological scenario under study here and for their future application to other areas and schistosome species.

Elucidating the temporal and spatial dynamics of Biomphalaria glabrata genetic diversity in three Brazilian villages

OBJECTIVE

The freshwater snail Biomphalaria glabrata is the principal intermediate host for the parasite Schistosoma mansoni within Brazil. We assessed the potential effects of snail population dynamics on parasite transmission dynamics via population genetics.

METHODS

We sampled snail populations located within the confines of three schistosome-endemic villages in the state of Minas Gerais, Brazil. Snails were collected from individual microhabitats following seasonal periods of flood and drought over the span of 1 year. Snail spatio-temporal genetic diversity and population differentiation of 598 snails from 12 sites were assessed at seven microsatellite loci.

RESULTS

Average genetic diversity was relatively low, ranging from 4.29 to 9.43 alleles per locus, and overall, subpopulations tended to exhibit heterozygote deficits. Genetic diversity was highly spatially partitioned among subpopulations, while virtually, no partitioning was observed across temporal sampling. Comparison with previously published parasite genetic diversity data indicated that S. mansoni populations are significantly more variable and less subdivided than those of the B. glabrata intermediate hosts.

DISCUSSION

Within individual Brazilian villages, observed distributions of snail genetic diversity indicate temporal stability and very restricted gene flow. This is contrary to observations of schistosome genetic diversity over the same spatial scale, corroborating the expectation that parasite gene flow at the level of individual villages is likely driven by vertebrate host movement.

Assessment of the genetic relationship between Dictyocaulus species from Bos taurus and Cervus elaphus using complete mitochondrial genomic datasets

BACKGROUND

Dictyocaulus species are strongylid nematodes of major veterinary significance in ruminants, such as cattle and cervids, and cause serious bronchitis or pneumonia (dictyocaulosis or "husk"). There has been ongoing controversy surrounding the validity of some Dictyocaulus species and their host specificity. Here, we sequenced and characterized the mitochondrial (mt) genomes of Dictyocaulus viviparus (from Bos taurus) with Dictyocaulus sp. cf. eckerti from red deer (Cervus elaphus), used mt datasets to assess the genetic relationship between these and related parasites, and predicted markers for future population genetic or molecular epidemiological studies.

METHODS

The mt genomes were amplified from single adult males of D. viviparus and Dictyocaulus sp. cf. eckerti (from red deer) by long-PCR, sequenced using 454-technology and annotated using bioinformatic tools. Amino acid sequences inferred from individual genes of each of the two mt genomes were compared, concatenated and subjected to phylogenetic analysis using Bayesian inference (BI), also employing data for other strongylids for comparative purposes.

RESULTS

The circular mt genomes were 13,310 bp (D. viviparus) and 13,296 bp (Dictyocaulus sp. cf. eckerti) in size, and each contained 12 protein-encoding, 22 transfer RNA and 2 ribosomal RNA genes, consistent with other strongylid nematodes sequenced to date. Sliding window analysis identified genes with high or low levels of nucleotide diversity between the mt genomes. At the predicted mt proteomic level, there was an overall sequence difference of 34.5% between D. viviparus and Dictyocaulus sp. cf. eckerti, and amino acid sequence variation within each species was usually much lower than differences between species. Phylogenetic analysis of the concatenated amino acid sequence data for all 12 mt proteins showed that both D. viviparus and Dictyocaulus sp. cf. eckerti were closely related, and grouped to the exclusion of selected members of the superfamilies Metastrongyloidea, Trichostrongyloidea, Ancylostomatoidea and Strongyloidea.

CONCLUSIONS

Consistent with previous findings for nuclear ribosomal DNA sequence data, the present analyses indicate that Dictyocaulus sp. cf. eckerti (red deer) and D. viviparus are separate species. Barcodes in the two mt genomes and proteomes should serve as markers for future studies of the population genetics and/or epidemiology of these and related species of Dictyocaulus.

The molecular epidemiology of parasite infections: tools and applications

Molecular epidemiology, broadly defined, is the application of molecular genetic techniques to the dynamics of disease in a population. In this review, we briefly describe molecular and analytical tools available for molecular epidemiological studies and then provide an overview of how they can be applied to better understand parasitic disease. A range of new molecular tools have been developed in recent years, allowing for the direct examination of parasites from clinical or environmental samples, and providing access to relatively cheap, rapid, high throughput molecular assays. At the same time, new analytical approaches, in particular those derived from coalescent theory, have been developed to provide more robust estimates of evolutionary processes and demographic parameters from multilocus, genotypic data. To date, the primary application of molecular epidemiology has been to provide specific and sensitive identification of parasites and to resolve taxonomic issues, particularly at the species level and below. Population genetic studies have also been used to determine the extent of genetic diversity among populations of parasites and the degree to which this diversity is associated with different host cycles or epidemiologically important phenotypes. Many of these studies have also shed new light on transmission cycles of parasites, particularly the extent to which zoonotic transmission occurs, and on the prevalence and importance of mixed infections with different parasite species or intraspecific variants (polyparasitism). A major challenge, and one which is now being addressed by an increasing number of studies, is to find and utilize genetic markers for complex traits of epidemiological significance, such as drug resistance, zoonotic potential and virulence.

Genetic structure of Schistosoma mansoni in western Kenya: The effects of geography and host sharing

We examined the spatial structure of Schistosoma mansoni, a parasite of humans, from natural infections at two levels: across the Lake Victoria basin of Kenya and among snail hosts. Using 20 microsatellite markers we examined geographic patterns of relatedness and population structure of cercariae and found weak, but significant structure detected by some, but not all analyses. We hypothesise structure created by aggregations of clonal individuals or adherence of hosts to local transmission sites is eroded by high amounts of gene flow in the region. This finding also supports previous hypotheses concerning the evolution of drug resistance in the region. Intrasnail dynamics were investigated in the context of aggregation and kin selection theory to determine how relatedness and also sex influence host sharing and host exploitation. Cercarial production did not differ significantly between snails with one or two genotypes suggesting that mixed infections resulted in decreased individual fitness and provides a framework for reproductive competition. Coinfection patterns in snails were independent of parasite relatedness indicating that schistosomes were not aggregated according to their relatedness and that kin selection was not influencing host sharing. Additionally, host exploitation in coinfections (measured by cercarial production) was not negatively correlated with relatedness, as predicted by classical models due to increased competition and thus exploitation when parasites are unrelated. Because of the low levels of relatedness within the population, schistosomes may rarely encounter close relatives and kin selection mechanisms that influence the distribution of individuals within snails or the virulence mode of the parasites may simply have not evolved.

Genetics of mating and sex determination in the parasitic nematode Haemonchus contortus

Genetic analysis of parasitic nematodes has been a neglected area of research and the basic genetics of this important group of pathogens are poorly understood. Haemonchus contortus is one of the most economically significant livestock parasites worldwide and is a key experimental model for the strongylid nematode group that includes many important human and animal pathogens. We have undertaken a study of the genetics and the mode of mating of this parasite using microsatellite markers. Inheritance studies with autosomal markers demonstrated obligate dioecious sexual reproduction and polyandrous mating that are reported here for the first time in a parasitic helminth and provide the parasite with a mechanism of increasing genetic diversity. The karyotype of the H. contortus, MHco3(ISE) isolate was determined as 2n = 11 or 12. We have developed a panel of microsatellite markers that are tightly linked on the X chromosome and have used them to determine the sex chromosomal karyotype as XO male and XX female. Haplotype analysis using the X-chromosomal markers also demonstrated polyandry, independent of the autosomal marker analysis, and enabled a more direct estimate of the number of male parental genotypes contributing to each brood. This work provides a basis for future forward genetic analysis on H. contortus and related parasitic nematodes.

Genetic variation of fish parasite populations in historically connected habitats: undetected habitat fragmentation effect on populations of the nematode Procamallanus fulvidraconis in the catfish Pelteobagrus fulvidraco

Habitat fragmentation may have some significant effects on population genetic structure because geographic distance and physical barriers may impede gene flow between populations. In this study, we investigated whether recent habitat fragmentation affected genetic structure and diversity of populations of the nematode Procamallanus fulvidraconis in the yellowhead catfish, Pelteobagrus fulvidraco. The nematode was collected from 12 localities in 7 floodplain lakes of the Yangtze River. Using 11 intersimple sequence repeat markers, analysis of molecular variance showed that genetic diversity occurred mainly within populations (70.26%). Expected heterozygosity (He) of P. fulvidraconis was barely different between connected (0.2105) and unconnected lakes (0.2083). Population subdivision (Fst) between connected lakes (0.2177) was higher than in unconnected lakes (0.1676). However, the connected and unconnected lakes did not cluster into 2 clades. A Mantel test revealed significant positive correlation between genetic and geographic distances (R = 0.5335, P < 0.01). These results suggest that habitat fragmentation did not cause genetic differentiation among populations or a reduction of diversity in isolated populations of P. fulvidraconis. At least 2 factors may increase the dispersal range of the nematode, i.e., flash flooding in summer and other species of fish that may serve as the definitive hosts. Moreover, lake fragmentation is probably a recent process; population size of the nematode in these lakes is large enough to maintain population structure.

Genetic diversity and population structuring of Schistosoma mansoni in a Brazilian village

The digenean trematode Schistosoma mansoni is responsible for chronic schistosomiasis worldwide, and in Brazil alone an estimated 35 million people are at risk. To evaluate epidemiological patterns among human definitive hosts, we assessed genetic diversity and population subdivision of S. mansoni infrapopulations in human hosts from the highly endemic village of Virgem das Graças in the state of Minas Gerais, Brazil. We believe this is the largest such survey to date. Genetic diversity of parasites, measured over eight polymorphic microsatellite loci, was relatively high and standard measures of inbreeding indicated that the population was panmictic. Furthermore, there was no significant isolation-by-distance of parasite infrapopulations, and measures of population subdivision indicated significant but low to moderate levels of population differentiation. We conclude that patients within this village sample from a broad range of schistosome genetic diversity and effectively act as "genetic mixing bowls" for the parasites. These results contrast with those previously observed in the Brazilian village of Melqui ades and thus provide the opportunity for comparisons of environmental and epidemiological differences that are likely to influence host-parasite coevolution and parasite transmission.

Macrogeographic population structure in a parasitic nematode with avian hosts

Much remains to be discovered about the population genetic structure of parasites, despite the importance of such knowledge to understanding the processes involved in the spread of drug resistance through populations. Here we present a study of population genetic diversity in Trichostrongylus tenuis, an avian parasitic nematode infecting both poultry and game birds, where anthelmintic use is common. We examined diversity of nicotinamide dehydrogenase subunit 4 (nad4) mtDNA sequences within and between seven locations: five in the UK (red grouse hosts), one in Iceland (domestic goose) and one in Norway (willow grouse). Within-UK comparisons showed high nucleotide diversity (pi=0.015, n=23) but no structure between locations (Phi(ST)=0.022, P=0.27), with over 97% of variation distributed within-hosts. The highest diversity was found in Iceland (pi=0.043, n=4), and the lowest in Norway (pi=0.003, n=4). Differentiation between countries was considerable (Phi(CT)=0.44, P<0.05), in spite of the potential mixing effects of gene flow via migrating wild hosts and the poultry trade. However, significant pairwise F(ST) values were found only between Norway and the other locations. Phylogenetic analysis provided statistical support for a separate clade for Norwegian samples only, with unresolved diversity leading to a star-shaped relationship between Icelandic and UK haplotypes. These results suggest that Norwegian T. tenuis are isolated, but that there is some connectivity between UK and Icelandic populations. Although anthelmintic resistance has not yet been reported for T. tenuis, the population structure is such that emerging resistance has the potential to spread by gene flow over a large geographic scale.

Parasites: proxies for host genealogy and ecology?

Genetic information is used extensively to reconstruct the evolutionary and demographic history of organisms. Recently, it has been suggested that genetic information from some parasites can complement genetic data from their hosts. This approach relies upon the hypothesis that such parasites share a common history with their host. In some cases, parasites provide an additional source of information because parasite data can better reconstruct the common history. Here, we discuss which parasite traits are important in determining their usefulness for analysing host history. The key is the matching of the traits of the parasite (e.g. effective population size, generation time, mutation rate and level of host specificity) with the timescales (phylogenetic, phylogeographic and demographic) that are relevant to the issues of concern in host history.

Evidence for multiple mitochondrial lineages of Fasciola hepatica (liver fluke) within infrapopulations from cattle and sheep

The economic, veterinary, and medical impact of the parasite Fasciola hepatica, liver fluke, is difficult to alleviate due to increasing incidences of resistance to the principal anthelmintic drugs. These have occurred in widely separated regions. The rate of response to selection imposed by such drugs will be dependent on the genetic variation present in the F. hepatica gene pool, but this is at present unknown. We have assessed the genetic diversity of mitochondrial haplotypes found in the infrapopulation of flukes recovered from a calf of known provenance and from six other cattle and sheep hosts located in Ireland and four from elsewhere. Our results revealed that at least ten different mitochondrial composite PCR-restriction fragment length polymorphism haplotypes had been acquired by a single animal in 1 year, and there was comparable diversity in six other definitive hosts carrying field-acquired infections. The extent of divergence between these fluke lineages suggests that they predate the last ice age and, thus, cannot have developed in Northern Europe. A consequence of this high level of diversity is that there will be frequent selection for anthelmintic resistance and rapid responses to climatic changes.

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.

Patterns of nucleotide polymorphism distinguish temperate and tropical wild isolates of Caenorhabditis briggsae

Caenorhabditis briggsae provides a natural comparison species for the model nematode C. elegans, given their similar morphology, life history, and hermaphroditic mode of reproduction. Despite C. briggsae boasting a published genome sequence and establishing Caenorhabditis as a model genus for genetics and development, little is known about genetic variation across the geographic range of this species. In this study, we greatly expand the collection of natural isolates and characterize patterns of nucleotide variation for six loci in 63 strains from three continents. The pattern of polymorphisms reveals differentiation between C. briggsae strains found in temperate localities in the northern hemisphere from those sampled near the Tropic of Cancer, with diversity within the tropical region comparable to what is found for C. elegans in Europe. As in C. elegans, linkage disequilibrium is pervasive, although recombination is evident among some variant sites, indicating that outcrossing has occurred at a low rate in the history of the sample. In contrast to C. elegans, temperate regions harbor extremely little variation, perhaps reflecting colonization and recent expansion of C. briggsae into northern latitudes. We discuss these findings in relation to their implications for selection, demographic history, and the persistence of self-fertilization.

Spiculopteragia mathevossianiRuchliadev, 1948 is the minor morph ofSpiculopteragia spiculoptera(Gushanskaya, 1931): molecular evidence

Although Spiculopteragia spiculoptera is primarily a parasite of cervids, it can also but less often contaminate domestic livestock. Little is known about its epidemiology and its pathogenicity in domestic ruminants and other unusual cervid species. Its taxonomic status remains unclear since the hypothesis of morphological polymorphism among males has been proposed. However, accurate taxonomy is fundamental in the identification and survey of potentially pathogenic species of parasites. The second internal transcribed spacer of rDNA (ITS-2) and the mitochondrial (mt) DNA-derived ND4 gene were used to study the polymorphism hypothesis for S. spiculoptera. ND4 evolves more quickly than ITS-2 and is considered to be more discriminant in the characterization of closely related species. DNA sequences of ITS-2 and ND4 were studied in 18 individual males of morphological type spiculoptera and in 3 of morphological type of mathevossiani from Red deer (Cervus elaphus), Roe deer (Capreolus capreolus) and Chamois (Rupicapra rupicapra). Intraindividual ITS-2 variations were detected within and between each morphotype of Spiculopteragia but these differences did not separate the two morphs mathevossiani and spiculoptera. Similarly, although ND4 showed a high level of nucleotide substitution, the morphotypes S. mathevossiani and S. spiculoptera were clustered together. Our genetic data support the dimorphic male hypothesis for the species S. spiculoptera.

GENETIC VARIATION AND RELATIONSHIPS BETWEEN ISOLATES AND SPECIES OF THE ENTOMOPATHOGENIC NEMATODE GENUS HETERORHABDITIS DECIPHERED THROUGH ISOZYME PROFILES

We studied variation in isozyme patterns of 8 metabolic enzymes in 5 species of Heterorhabditis (H. bacteriophora, H. indica, H. marelata, H. megidis, and H. zealandica) comprising 18 isolates. Isozyme banding patterns of all the 8 enzymes were species specific; however, 3 enzymes, i.e., arginine kinase, fumarate hydratase, and malate dehydrogenase, displayed distinct patterns among all the 18 isolates. Phylogenetic analysis of the isozyme patterns produced dendrograms depicting a high degree of genetic variation among Heterorhabditis species, with the average pairwise distance of 0.2000. Trees constructed using different phylogenetic methods showed a relatively close genetic relationship between H. megidis and H. zealandica and between H. bacteriophora and H. indica. Also, H. bacteriophora HP88 was the most distant species from H. megidis (UK isolate), H. marelatus (Oregon isolate), and H. zealandica (X1 isolate) with pairwise distance of 0.1957, 0.2228, and 0.2120, respectively. Phylogenetic analysis also revealed genetic variation among H. bacteriophora isolates with the average pairwise distance of 0.1507. GPS2 and GPS3 were the most closely related isolates with the average distance of only 0.0870, followed by GPS1 and GPS2 with average distance of 0.1087. In contrast, KMD19 and HP88, OH25, and HP88, and OH25 and Acows isolates were the most divergent populations with a pairwise distance of 0.2011 and 37 character differences. Pairwise distance analysis also revealed that genetic divergence among populations of H. bacteriophora is relatively independent of geographic distance. Overall, these results demonstrate strong subspecies structuring in H. bacteriophora.

Population genetics of the bovine/cattle lungworm (Dictyocaulus viviparus) based on mtDNA and AFLP marker techniques

Mitochondrial DNA (mtDNA) sequence data and amplified fragment length polymorphism (AFLP) patterns were compared for the lungworm Dictyocaulus viviparus, a nematode parasite of cattle. Eight individual D. viviparus samples from each of 8 herds in Sweden and 1 laboratory isolate were analysed, with the aim of describing the diversity and genetic structure in populations using different genetic markers on exactly the same DNA samples. There was qualitative agreement between the whole-genome AFLP data and the mtDNA sequence data, both indicating relatively strong genetic differentiation among the Swedish farms. However, the AFLP data detected much more genetic variation than did the mtDNA data, even after allowing for the different inheritance patterns of the markers, and indicated that there was much less differentiation among the populations. The mtDNA data therefore seemed to be more informative about the most recent history of the parasite populations, as the general patterns were less obscured by detailed inter-relationships among individual worms. The 4 mtDNA genes sequenced (1542 bp) showed consistent patterns, although there was more genetic variation in the protein-coding genes than in the structural RNA genes. Furthermore, there appeared to be at least 3 distinct genetic groups of D. viviparus infecting Swedish cattle, 1 of which was predominant and showed considerable differentiation between farms, but not necessarily within farms. Second, the 2 smaller genetic groups occurred on farms where the predominant group also occurred, suggesting that these farms have had multiple introductions of D. viviparus.

The complete mitochondrial genome of the entomopathogenic nematode Steinernema carpocapsae: insights into nematode mitochondrial DNA evolution and phylogeny

We determined the complete sequence of the mitochondrial DNA of the entomopathogenic nematode Steinernema carpocapsae and analyzed its structure and composition as well as the secondary structures predicted for its tRNAs and rRNAs. Almost the complete genome has been amplified in one fragment with long PCR and sequenced using a shotgun strategy. The 13,925-bp genome contains genes for 2 rRNAs, 22 tRNAs, and 12 proteins and lacks an ORF encoding ATPase subunit 8. Four initiation codons were inferred, TTT, TTA, ATA, and ATT, most of the genes ended with TAA or TAG, and only two had a T as an incomplete stop codon. All predicted tRNAs showed the nonconventional secondary structure typical of Secernentea. Although we were able to fold the sequences of trnN, trnD, and trnC into more conventional cloverleaf structures after adding adjacent nucleotides, northern blot experiments showed that the nonstandard tRNAs are actually expressed. Phylogenetic and comparative analyses showed that the mitochondrial genome of S. carpocapsae is more closely related to the genomes of A. suum and C. elegans than to that of Strongyloides stercoralis. This finding does not support the phylogeny based on nuclear small subunit ribosomal DNA sequences previously published. This discrepancy may result from differential reproductive strategies and/or differential selective pressure acting on nuclear and mitochondrial genes. The distinctive characteristics observed among mitochondrial genomes of Secernentea may have arisen to counteract the deleterious effects of Muller's ratchet, which is probably enhanced by the reproductive strategies and selective pressures referred to above.

Mitochondrial genomes of parasitic nematodes – progress and perspectives

Mitochondria are subcellular organelles in which oxidative phosphorylation and other important biochemical functions take place within the cell. Within these organelles is a mitochondrial (mt) genome, which is distinct from, but cooperates with, the nuclear genome of the cell. Studying mt genomes has implications for various fundamental areas, including mt biochemistry, physiology and molecular biology. Importantly, the mt genome is a rich source of markers for population genetic and systematic studies. To date, more than 696 mt genomes have been sequenced for a range of metazoan organisms. However, few of these are from parasitic nematodes, despite their socioeconomic importance and the need for fundamental investigations into areas such as nematode genetics, systematics and ecology. In this article, we review knowledge and recent progress in mt genomics of parasitic nematodes, summarize applications of mt gene markers to the study of population genetics, systematics, epidemiology and evolution of key nematodes, and highlight some prospects and opportunities for future research.

Analysis of genetic variability within Thelazia callipaeda (Nematoda: Thelazioidea) from Europe and Asia by sequencing and mutation scanning of the mitochondrial cytochrome c oxidase subunit 1 gene

This study investigated genetic variability within the 'eyeworm'Thelazia callipaeda (Nematoda: Thelazioidea) from Europe and Asia by polymerase chain reaction (PCR)-coupled sequencing and mutation scanning of the mitochondrial cytochrome c oxidase subunit 1 gene (cox 1). Eight different sequence variants of cox 1 (haplotypes) were determined for the 50 individual adult specimens of T. callipaeda (from dogs, foxes or cats from Italy, Germany and the Netherlands and from dogs from China and Korea). Nucleotide variation (0.3--2%) was detected at 23 of 649 positions in the cox 1. Six of these positions were invariable among all 37 individuals from Europe and among the 13 individuals from Asia (irrespective of host origin) but differed (five G<-->A and one C<-->T changes) between Europe and Asia. PCR-based single-strand conformation polymorphism (SSCP) analysis of the most variable portion (v-cox 1) of the cox 1 was validated (for a subset of samples) as a tool to rapidly screen for genetic (haplotypic) variability. The results for the SSCP analysis and sequencing were concordant, indicating that the mutation scanning approach provides a useful tool for investigating the population genetics and molecular ecology of T. callipaeda.

Speciation in parasites: a population genetics approach

Parasite speciation and host-parasite coevolution should be studied at both macroevolutionary and microevolutionary levels. Studies on a macroevolutionary scale provide an essential framework for understanding the origins of parasite lineages and the patterns of diversification. However, because coevolutionary interactions can be highly divergent across time and space, it is important to quantify and compare the phylogeographic variation in both the host and the parasite throughout their geographical range. Furthermore, to evaluate demographic parameters that are relevant to population genetics structure, such as effective population size and parasite transmission, parasite populations must be studied using neutral genetic markers. Previous emphasis on larger-scale studies means that the connection between microevolutionary and macroevolutionary events is poorly explored. In this article, we focus on the spatial fragmentation of parasites and the population genetics processes behind their diversification in an effort to bridge the micro- and macro-scales.

Testing the hypothesis of recent population expansions in nematode parasites of human-associated hosts

It has been predicted that parasites of human-associated organisms (eg humans, domestic pets, farm animals, agricultural and silvicultural plants) are more likely to show rapid recent population expansions than are parasites of other hosts. Here, we directly test the generality of this demographic prediction for species of parasitic nematodes that currently have mitochondrial sequence data available in the literature or the public-access genetic databases. Of the 23 host/parasite combinations analysed, there are seven human-associated parasite species with expanding populations and three without, and there are three non-human-associated parasite species with expanding populations and 10 without. This statistically significant pattern confirms the prediction. However, it is likely that the situation is more complicated than the simple hypothesis test suggests, and those species that do not fit the predicted general pattern provide interesting insights into other evolutionary processes that influence the historical population genetics of host-parasite relationships. These processes include the effects of postglacial migrations, evolutionary relationships and possibly life-history characteristics. Furthermore, the analysis highlights the limitations of this form of bioinformatic data-mining, in comparison to controlled experimental hypothesis tests.

Phylogeography of a nematode (Heligmosomoides polygyrus) in the western Palearctic region: persistence of northern cryptic populations during ice ages?

This study establishes the continental phylogeographical pattern of a European nematode, Heligmosomoides polygyrus (Dujardin, 1845; Heligmosomoidea). We sequenced 687 base pairs of the mitochondrial DNA (mtDNA) cyt b gene for 136 individuals collected in 22 localities. The results revealed that H. polygyrus populations are separated into five major units corresponding to the Italian, northern European (Denmark and Ireland), Iberian, western European, and Balkan populations. Different subclades were also observed within the first two groups. Based on the rate of molecular evolution of H. polygyrus cyt b gene-estimated to 3.5%-3.7% divergence per million years (Myr) in a previous study--the isolation time of the five clades was estimated between 2.5 +/- 0.24 and 1.5 +/- 0.23 million years ago. Moreover, H. polygyrus presents a higher genetic variability in the Mediterranean peninsulas as compared to northwestern Europe, highlighting the role of these regions as refuge areas. Like its specific host, the wood mouse Apodemus sylvaticus, H. polygyrus' pattern of postglacial recolonization of northwestern Europe was initiated from Iberian populations, while Italian and Balkan populations did not expand to the north. The results also suggest the existence of forested and temperate refuges in the southern British Isles during the Quaternary. Finally, the genetic diversity as well as the level of genetic divergence between the lineages of H. polygyrus are compared to those observed in other vertebrate and invertebrate phylogeographical studies: the existence of highly differentiated lineages in H. polygyrus (5%-10% of genetic divergence) highlights that the effects of Pleistocene climate changes on free-living organisms are also reflected in their obligate parasites.

Host isolation and patterns of genetic variability in three populations of Teladorsagia from sheep

We have used a mitochondrial marker to explore the population genetics of an economically important parasite of sheep, Teladorsagia. We examined diversity within and between parasites from three very different host populations, as well as within and between individual hosts. One of our study populations, the Soay sheep on Hirta, St Kilda, is unusually isolated with no sheep having been introduced to the island since 1932. Worm haplotypes from Hirta were compared with those from two other host populations. Remarkably, despite its historical isolation the Hirta population shows similar levels of within-population diversity to the other study populations. No divergence between the three Teladorsagia populations was found, consistent with gene flow between the populations. The high diversity within Teladorsagia populations provides compelling evidence that this variability is a general feature of parasitic nematode populations. Such diversity may be caused by high effective population size, coupled with an increased mutation rate for mtDNA, which has important implications for the spread of anthelmintic resistance in nematode populations.