Genetic diversity assessed by amplified fragment length polymorphism analysis of the parasitic nematode Dictyocaulus viviparus the lungworm of cattle

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.

A review of methods for nematode identification

Nematodes are non-segmented roundworms found in soil, aquatic environment, plants, or animals. Either useful or pathogenic, they greatly influence environmental equilibrium, human and animal health, as well as plant production. Knowledge on their taxonomy and biology are key issues to answer the different challenges associated to these organisms. Nowadays, most of the nematode taxonomy remains unknown or unclear. Several approaches are available for parasite identification, from the traditional morphology-based techniques to the sophisticated high-throughput sequencing technologies. All these techniques have advantages or drawbacks depending on the sample origin and the number of nematodes to be processed. This review proposes an overview of all newly available methods available to identify known and/or unknown nematodes with a specific focus on emerging high-throughput molecular techniques.

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.

Dictyocaulus viviparus genome, variome and transcriptome elucidate lungworm biology and support future intervention

The bovine lungworm, Dictyocaulus viviparus (order Strongylida), is an important parasite of livestock that causes substantial economic and production losses worldwide. Here we report the draft genome, variome, and developmental transcriptome of D. viviparus. The genome (161 Mb) is smaller than those of related bursate nematodes and encodes fewer proteins (14,171 total). In the first genome-wide assessment of genomic variation in any parasitic nematode, we found a high degree of sequence variability in proteins predicted to be involved host-parasite interactions. Next, we used extensive RNA sequence data to track gene transcription across the life cycle of D. viviparus, and identified genes that might be important in nematode development and parasitism. Finally, we predicted genes that could be vital in host-parasite interactions, genes that could serve as drug targets, and putative RNAi effectors with a view to developing functional genomic tools. This extensive, well-curated dataset should provide a basis for developing new anthelmintics, vaccines, and improved diagnostic tests and serve as a platform for future investigations of drug resistance and epidemiology of the bovine lungworm and related nematodes.

PGP expression in Cooperia oncophora before and after ivermectin selection

The aim of this study was to investigate genetic selection and P-glycoprotein (PGP) expression in three different isolates of Cooperia oncophora before treatment and after ivermectin (IVM) injection. Adult parasites were recovered from nine calves experimentally infected with the isolates represented by one IVM susceptible laboratory isolate, and two field isolates showing signs of phenotypic macrocyclic lactone resilience according to the faecal egg count reduction test. Five males and five females per isolate were examined both pre- and post-IVM treatment giving a total of 60 worms. A sequence from C. oncophora (Con-pgp) was identified, showing 83 % similarity to Pgp-9 of Caenorhabditis elegans. Primers specific to putative Con-pgp-9 mRNA were designed, generating a 153-bp PCR product. Total RNA was prepared from all worms, and Con-pgp-9 expression was measured by quantitative real-time reverse transcription PCR. Our results showed that mean PGP concentrations were four to five times higher in female as compared to male worms. No significant differences in gene expression between experimental groups pre- and post-IVM selection were detected. However, PGP gene expression tended to be increased by IVM treatment in male worms (p = 0.091), with 70 % higher mean expression in treated than in untreated male worms. Amplified fragment length polymorphism analysis did not demonstrate any bottleneck effect within the different isolates post-treatment. The total mean gene diversity values were 0.158 and 0.153 before and after treatment, respectively. Inbreeding coefficient in subpopulations compared to total population F_ST was 0.0112, suggesting no genetic differentiation between or within the investigated isolates in relation to treatment. In conclusion, comparison of Con-pgp-9 expression showed no significant difference before and after treatment, but some tendency towards increasing expression in male worms.

Population genetic structure of Ascaridia galli re-emerging in non-caged laying hens

Background

The poultry roundworm Ascaridia galli has reappeared in hens kept for egg production in Sweden after having been almost absent a decade ago. Today this is a frequent intestinal nematode parasite in non-caged laying hens. The aim of this study was to investigate the genetic diversity (F_st) in A. galli collected from different poultry production sites in southern Sweden, to identify possible common routes of colonization.

Methods

Adult parasites (n = 153) from 10 farms, including both broiler breeder parents and laying hens, were investigated by amplified restriction fragment length polymorphism analysis (AFLP). Worms from a Danish laying hen farm were also included for comparison. Most of the farms were represented by worms from a single host, but on two farms multiple samples from different hosts were assessed in order to study flock variation.

Results

A total of 97 fragments (loci) were amplified among which 81% were variable alleles. The average genetic diversity was 0.13 (range = 0.09-0.38), which is comparable to other AFLP studies on nematodes of human and veterinary importance. Within-farm variation showed that worms harboured by a single hen in a flock covered most of the A. galli genetic variation within the same flock (F_st = 0.01 and 0.03 for two farms). Between-farm analysis showed a moderate population genetic structure (F_st = 0.13), along with a low mutational rate but high gene flow between different farms, and absence of strong genetic selection. Network analysis showed repeated genetic patterns among the farms, with most worms on each farm clustering together as supported by high re-allocation rates.

Conclusions

The investigated A. galli populations were not strongly differentiated, indicating that they have undergone a genetic bottlenecking and subsequent drift. This supports the view that the investigated farms have been recently colonized, and that new flocks are reinfected upon arrival with a stationary infection.

Control of nematode parasites with agents acting on neuro-musculature systems: lessons for neuropeptide ligand discovery

There are a number of reasons why the development of novel anthelmintics is very necessary. In domestic animals, parasites cause serious loss of production and are a welfare concern. The control of these parasites requires changes in management practices to reduce the spread of infection and the use of therapeutic agents to treat affected animals. The development of vaccines against parasites is desirable but their development so far has been very limited. One notable exception is the vaccination of calves against infection by Dictyocaulus viviparous (lungworm) which has proved to be very effective. In domestic animals, the total market for anti-parasitic agents (both ecto- and endo-parasites) is in excess of a billion U.S. dollars. In humans there are serious problems ofmorbidity and mortality associated with parasite infections. 1.6 billion People throughout the world are infected with ascariasis (Fig. 1A) and/or hookworm. Approximately one-third of the world's population is suffering from the effects of intestinal nematode parasites, causing low growth-rates in infants, ill-thrift, diarrhea and in 2% of cases, loss of life. Despite the huge number of affected individuals, the market for anti-parasitic drugs for humans is not big enough to foster the development of anthelmintics because most infestations that occur are in undeveloped countries that lack the ability to pay for the development of these drugs. The major economic motivator then, is for the development of animal anthelmintics. In both domestic animals and now in humans, there is now a level of resistance to the available anthelmintic compounds. The resistance is either: constitutive, where a given species of parasite has never been sensitive to the compound; or acquired, where the resistance has developed through Darwinian selection fostered by the continued exposure to the anti-parasitic drugs. The continued use of all anthelmintics has and will, continue to increase the level of resistance. Cure rates are now often less than 100% and resistance of parasites to agents acting on the neuromuscular systems is present in a wide range of parasites of animals and humans hosts. In the face of this resistance the development of novel and effective agents is an urgent and imperative need. New drugs which act on the neuromuscular system have an advantage for medication for animals and humans because they have a rapid therapeutic effect within 3 hours of administration. The effects on the neuromuscular system include: spastic paralysis with drugs like levamisole and pyrantel; flaccid paralysis as with piperazine; or disruption of other vital muscular activity as with ivermectin. Figure 1 B and C, illustrates an example ofa spastic effect oflevamisole on infectious L3 larvae of Ostertagia ostertagiae, a parasite of pigs. The effect was produced within minutes of the in vitro application oflevamisole. In this chapter we comment on the properties of existing agents that have been used to control nematode parasites and that have an action on neuromuscular systems. We then draw attention to resistance that has developed to these compounds and comment on their toxicity and spectra of actions. We hope that some of the lessons that the use of these compounds has taught us may to be applied to any novel neuropeptide ligand that may be introduced. Our aim is then is to provide some warning signs for recognized but dangerous obstacles.

Limited sequence variation in the major sperm protein 1 (MSP) gene within populations and species of the genus Dictyocaulus (Nematoda)

Populations of the bovine lungworm, Dictyocaulus viviparus, are genetically structured based on variation in mtDNA and AFLP data. Our aim was to investigate if this genetic variability also is reflected in a protein recognized by the host immune system. We focused on the major sperm protein (MSP), a small and abundant protein used in diagnostic immunoassays, which has been shown to be variable in some nematodes but not others. MSP was sequenced using worm DNA from eight adult worms from each of nine populations whose genetic structure previously had been quantified. For comparison, we also analyzed MSP sequences of the closely related Dictyocaulus eckerti and Dictyocaulus capreolus and from nematodes with sequences deposited in GenBank. In contrast to previous results, this study shows that the MSP ofD. viviparus is similar to that of other nematodes. Almost no sequence variation, and thus no antigenic diversity, was detected in MSP between worms from different sub-populations or in the other Dictyocaulus species investigated. A functional test of a recombinant variant of the MSP showed that the expressed protein was recognized by antibodies in sera from infected cattle. This has practical implications for the development of species-specific markers, recombinant vaccines, and immunodiagnostics.

AFLP fingerprinting for the analysis of genetic diversity within Necator americanus

In the present study, we utilised the method of AFLP to screen for genetic variation within and among individuals of the blood-feeding human hookworm Necator americanus (Nematoda) from Africa, Asia and South America. A total of 45 adult worms (i.e. 20 from Ghana, 16 from Colombia and 9 from Nepal) were subjected to analysis using the restriction enzyme/primer combination HindIII+AG/BglII+AC. Cluster analysis divided N. americanus into multiple, genetically distinct groups, consistent with previous findings using ribosomal and mitochondrial DNA data sets. The results demonstrated the usefulness of AFLP fingerprinting for establishing genetic variation within N. americanus and reinforce its applicability to other parasitic helminths of human and/or veterinary health importance.

Global patterns reveal strong population structure in Haemonchus contortus, a nematode parasite of domesticated ruminants

We have examined the global population genetic structure of Haemonchus contortus. The genetic variability was studied using both amplified fragment length polymorphism (AFLP) and nad4 sequences of the mitochondrial genome. To examine the performance and information content of the two different marker systems, comparative assessment of population genetic diversity was undertaken in 19 isolates of H. contortus, a parasitic nematode of small ruminants. A total of 150 individual adult worms representing 14 countries from all inhabited continents were analysed. Altogether 1,429 informative AFLP markers were generated using four different primer combinations. Also, the genetic variation was high, which agrees with results from previous AFLP studies of nematode parasites of livestock. The genetic structure was high, indicating limited gene flow between the different isolates and populations from each continent mostly formed monophyletic groups in the phylogenetic analysis. However, for isolates representing Australia, Greece and one laboratory strain that originated from South Africa (WRS), there was no clear genetic relationship between the isolates and the distance between their geographical origins. Basically the same pattern was observed for the mitochondrial marker, although the phylogenetic analysis was less resolved than for AFLP. In contrast with previous findings on the population genetic structure of H. contortus, the calculation of population structure gave high values (Nst=0.59). The strong structure was present also for the four Swedish isolates (Nst=0.16) representing a small geographical area.

Seasonal dynamics of population genetic structure in cryptic taxa of the Pellioditis marina complex (Nematoda: Rhabditida)

The distribution patterns and genetic structure of the Pellioditis marina species complex in Belgium and The Netherlands were compared between four consecutive seasons. Different types of habitats (coast, estuary, semi-estuary and lake) with different degrees of connectivity were sampled. In addition, each habitat type was characterised by either temporal or permanent algal deposits. We screened 426 bp of the mitochondrial cytochrome oxidase c (COI) gene with the single-strand conformation polymorphism (SSCP) method in 1615 individuals of Pellioditis marina. The 51 haplotypes were divided into four (sympatric) lineages, with divergences ranging from 0.25 to 10.6%. Our results show that the lineages have different temporal dynamics, which may be linked to abiotic factors. Analysis of Molecular Variance (AMOVA) indicated a significant structuring in the PmI lineage, which correlated with habitat characteristics and which changed over time (Mantel, r = 0.51; p = 0.126). Intrapopulational diversity was similar in all locations, and temporal changes in haplotype frequencies were not higher in temporary than in permanent algal deposits. Instead, the results of the temporal survey indicated that (some) P. marina populations are characterised by a metapopulation structure. It is emphasized that a complete and correct interpretation of processes causing genetic structuring within species and of the genetic structure itself can only be done when analyses are performed at several time points.

Targeted selective treatment of lungworm infection in an organic dairy herd in Sweden

The effect of targeted selective anthelmintic treatment on the seroprevalence of the lungworm Dictyocaulus viviparus in cattle was investigated. The study was commenced on an organic dairy enterprise in Sweden in November 1998 after the observation of an outbreak dictyocaulosis in the herd, and then continued for almost 3 years. The first year sampling was conducted on a monthly basis and then biannually with the exception of between August and November 2000 when sampling was performed monthly following a second outbreak of dictyocaulosis. Throughout the study, blood samples were examined for specific IgG(1) levels from all animals in the herd that had been grazing for more than 3 months. At the first sampling occasion, 13% out of the 90 blood samples were seropositive. One month later, after targeted selective treatment with eprinomectin (Eprinex), Merial), the whole herd was seronegative. Seroprevalence then gradually increased and 1 year later it returned to levels similar to those observed at the start of the study. At turnout in April 2000, seroprevalence was 1.3% but it then rapidly increased to 28% and 30% in August and September, respectively. This increase was mainly due to an increase in FSG animals of which many were coughing. Consequently, all seropositive animals were injected with ivermectin (Ivomec), Merial) at 0.05 mg/kg body weight in late August 2000. Although all animals recovered, seroprevalence was only reduced to 12% 1 month later. The differences in seroprevalence after both of these anthelmintic treatments were probably attributed to the timing. The first deworming with eprinomectin was conducted in November when the infection already was transient, whereas ivermectin in connection with the second outbreak was injected in a more acute phase of the infection cycle. Infection levels in 2001 were low with seroprevalences of 2.3% and 5.6% in May and September, respectively. These results show that dictyocaulosis in Sweden can be effectively controlled by the use of macrocyclic lactones. However, the infection was not eradicated from the herd despite close monitoring of the seroprevalence and targeted selective treatment of every seropositive animal on two occasions.

Population structure in Ascaris suum (Nematoda) among domestic swine in Denmark as measured by whole genome DNA fingerprinting

We here analyze the population structure in the pig roundworm, Ascaris suum, among domestic pigs in Denmark using a whole-genome DNA fingerprinting technique, "amplified fragment length polymorphism" (AFLP) analysis. With these data, we can extract absolute gene frequency variance components and G-statistics for 135 independent nucleotide polymorphisms. The average proportion of total variance partitioned between Jutland and Zealand is less than 3% of the total variance, implying no restriction in gene flow between worms from different regions in Denmark. The average gene frequency difference between two farms widely separated in Jutland represents 5% of the total genetic variance of these two farms combined. Conversely, worms from different hosts within these two farms are more subdivided, with an average of 12% of the total variance in gene frequencies within farms being distributed between hosts. This result implies substantial single generation inbreeding due to founder effects in the establishment of adult worms in single hosts. Absolute variance components extracted from the gene diversities also showed significant differences, with the among-host variance being greater that the between-farm and between-region values. This little geographical variation is discussed in relation to the hierarchic structure of the Danish swine production system. Comparison of our results with other studies on parasitic roundworms, suggests that patterns of host dispersal effectively control patterns of worm gene flow. Furthermore, the potential spread of anthelminth resistance among A. suum may thus be rapid, due to the flow of infected hosts within the domestic swine stocks in Denmark.

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.

Molecular tools--advances, opportunities and prospects

Modern molecular technologies are having a substantial impact in many fundamental and applied areas of parasitology. In particular, polymerase chain reaction (PCR)-coupled approaches have found broad applicability because their sensitivity permits the enzymatic amplification of gene fragments from minute quantities of nucleic acids from tiny amounts of parasite material. Also, high-resolution electrophoretic and genomic methods are finding increased utility. This paper briefly discusses some developments and applications of DNA methods to parasites and highlights their usefulness or potential for those of veterinary importance. Selected examples of applications with implications in fundamental (systematics, population genetics, epidemiology and ecology) and applied (diagnosis, prevention and control) areas are presented. The focus is mainly on tools for the accurate identification of parasitic nematodes and protozoa of socio-economic importance, the diagnosis of infections and the detection of genetic variability using PCR-coupled mutation scanning technology.

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.

Effect of moxidectin selection on the genetic variation within Cylicocyclus nassatus based on amplified fragment length polymorphism (AFLP)

Cyathostomins are among the most important intestinal nematodes of horses, yet, the literature on the molecular genetics of these worms is scarce. In this study, the technique of amplified fragment length polymorphism (AFLP) was applied to study the genetic diversity as well as to determine the effect of moxidectin selection on the population genetic diversity for Cylicocyclus nassatus, one of the most common cyathostomin species. Genomic DNAs from 30 individual male worms were used from each of two populations: an avermectin-milbemycin (AM)-naive population (Population-S) and a population derived from Population-S following 21 treatments with moxidectin (Population-Mox). Three selective primer pairs were used for each worm, yielding a total of 229 AFLP markers. Calculation of average pair wise Jaccard indices revealed a high degree of genetic variation within both populations using all three primer combinations. In addition, selection by moxidectin during a 3-year period caused a significant decrease in the level of genetic diversity as evidenced by analysis of AFLP markers for two primer combinations but not for the third. A dendrogram of relationships among individuals based on AFLP markers did not show a clear classification of individuals in separate groups. It was concluded that a high degree of genetic intrapopulation variation exists in C. nassatus and that moxidectin selection has a significant effect on the genetic composition of C. nassatus.

High resolution DNA fingerprinting by AFLP to study the genetic variation among Oesophagostomum bifurcum (Nematoda) from human and non-human primates from Ghana

An AFLP approach was established to investigate genetic diversity within Oesophagostomum bifurcum (order Strongylida) from human and non-human primates. Evaluation of different combinations of restriction enzymes (n = 8) and primers (n = 29) demonstrated that the use of HindIII/BglII digested templates and primers with the selective nucleotides + AG/ +AC, respectively, was the most effective for the analysis of O. bifurcum DNA. A total of 63 O. bifurcum adults from human, Patas monkey, Mona monkey and Olive baboon hosts from different geographical regions in Ghana were subjected to analysis using this method. Cluster analysis revealed 4 genetically distinct groups, namely O. bifurcum from the Patas monkey (I), from the Mona monkey (II), from humans (III) and from the Olive baboon (IV). These findings were concordant with those achieved previously using RAPD analysis and supports population genetic substructuring within O. bifurcum according to host species. The results demonstrated the effectiveness of the present AFLP method for establishing genetic variation within O. bifurcum, and indicates its applicability to other parasitic nematodes of human and/or veterinary health importance.

Networks in phylogenetic analysis: new tools for population biology

Phylogenetic analysis has changed greatly in the past decade, including the more widespread appreciation of the idea that evolutionary histories are not always tree-like, and may, thus, be best represented as reticulated networks rather than as strictly dichotomous trees. Reconstructing such histories in the absence of a bifurcating speciation process is even more difficult than the usual procedure, and a range of alternative strategies have been developed. There seem to be two basic uses for a network model of evolution: the display of real but unobservable evolutionary events (i.e. a hypothesis of the true phylogenetic history), and the display of character conflict within the data itself (i.e. a summary of the data). These two general approaches are briefly reviewed here, and the strengths and weaknesses of the different implementations are compared and contrasted. Each network methodology seems to have limitations in terms of how it responds to increasing complexity (e.g. conflict) in the data, and therefore each is likely to be more appropriate for one of the two uses than for the other. Several examples using parasitological data sets illustrate the uses of networks within the context of population biology.