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

Distribution of large lungworms (Nematoda: Dictyocaulidae) in free-roaming populations of red deer Cervus elaphus (L.) with the description of Dictyocaulus skrjabini n. sp

Lungworms of the genus Dictyocaulus are causative agents of parasitic bronchitis in domestic and wild ungulates. This study investigates the distribution, morphology and genetic diversity of D. cervi and a new lungworm species, Dictyocaulus skrjabini n. sp. infecting red deer Cervus elaphus, fallow deer Dama dama and moose Alces alces in Poland and Sweden. The study was conducted on 167 red deer from Poland and on the DNA of lungworms derived from 7 fallow deer, 4 red deer and 2 moose collected in Sweden. The prevalence of D. cervi and D. skrjabini n. sp. in dissected red deer in Poland was 31.1% and 7.2%, respectively. Moreover, D. skrjabini n. sp. was confirmed molecularly in 7 isolates of fallow deer lungworms and 1 isolate of red deer lungworms from Sweden. Dictyocaulus skrjabini n. sp. was established based on combination of their distinct molecular and morphological features; these included the length of cephalic vesicle, buccal capsule (BC), buccal capsule wall (BCW), distance from anterior extremity to the nerve ring, the width of head, oesophagus, cephalic vesicle, BC and BCW, as well as the dimensions of reproductive organs of male and female. Additionally, molecular analyses revealed 0.9% nucleotide sequence divergence for 1,605 bp SSU rDNA, and 16.5–17.3% nucleotide sequence divergence for 642 bp mitochondrial cytB between D. skrjabini n. sp. and D. cervi, respectively, and 18.7–19% between D. skrjabini n. sp. and D. eckerti, which translates into 18.2–18.7% amino acid sequence divergence between D. skrjabini n. sp. and both lungworms.

Large lungworms (Nematoda: Dictyocaulidae) recovered from the European bison may represent a new nematode subspecies

Although the Dictyocaulus lungworm, the agent of dictyocaulosis, is one of parasitological threats to European bison, its systematic position remains unclear. The aim of the present study was to evaluate the morphological features of the lungworm and the pathological lesions it induces, and to analyse mitochondrial (mt) genetic markers for systematic and molecular epidemiological studies. The morphological findings indicate that Dictyocaulus lungworms of European bison can be distinguished from those of cattle on the basis of differences in buccal capsule wall length, total body length, and spicules length in males, all of which were significantly longer in those of European bison. Nucleotide diversity calculated from pairwise sequence alignments of partial cytochrome c oxidase subunit 1 (cox1), cytochrome B (cytB) and NADH dehydrogenase subunit 5 (nad5) of specimens from cattle and European bison varied from 1.7% for nad5, 2.1% for cytB, to 3.7% for cox1 gene. Thus, among the lungworms of European bison and cattle, nad5 and cytB were the most conserved proteins, whereas cox1 was the most diverse. The mt cytB marker gene may be a suitable candidate for distinguishing between the two genotypes, as nad5 demonstrated the greatest within-genus sequence variation. The lung tissue of infected European bison manifests signs of verminous pneumonia characterized by interstitial pneumonia, bronchitis and bronchiolitis. Therefore, it appears that European bison and cattle are infected with slightly diverged, morphologically-different, genotypes of D. viviparus, indicating they belong to two separate worm populations. We propose, therefore, that the lungworm of European bison should be classified as D. viviparus subsp. bisontis.

•

European bison harbour a unique morphotype and genotype of Dictyocaulus viviparus.

•

Mt cytB is an efficient genetic marker for studying large lungworms in bovids.

•

European bison lungworm can be classified as D. viviparus subsp. bisontis.

•

Pathologies induced by a newly described nematode subspecies mirrored those of lungworm in cattle.

An assessment of the use of cox1 and cox3 mitochondrial genetic markers for the identification of Dictyocaulus spp. (Nematoda: Trichostrongyloidea) in wild ruminants

Lungworms of the genus Dictyocaulus Railliet and Henry, 1907 (Nematoda: Trichostrongyloidea) are the causative agents of parasitic bronchitis (dictyocaulosis, husk) of various ungulate hosts, including domestic and wild ruminants. Correct diagnosis of lungworm species and a better understanding of the transmission patterns of Dictyocaulus spp. are crucial in minimising the risk of its cross transmission between wildlife and livestock, and for the control of dictyocaulosis. The study was conducted on large lungworms collected from European bison, roe deer and red deer. The study resulted in 14 sequences of the partial cox1 region of Dictyocaulus spp. and 10 novel DNA sequences of partial cox3 region, including the first available mt cox3 sequence, of the roe deer lungworm (D. capreolus). The European bison was infected with bison genotype of D. viviparus, whereas red deer and roe deer were infected with D. cervi and D. capreolus respectively. The current study revealed that the cox3 nucleotide sequences of D. capreolus and D. viviparus were 100% homologous to each other. Our findings indicate that the mt cox3 gene does not serve as an efficient mt marker for systematic, population genetic or molecular epidemiological studies of Dictyocaulus lungworms.

Interrelationships of Dictyocaulus spp. in Wild Ruminants with Morphological Description of Dictyocaulus cervi n. sp. (Nematoda: Trichostrongyloidea) from Red Deer, Cervus elaphus

Lungworms from the genus Dictyocaulus cause parasitic bronchitis (dictyocaulosis) characterized by coughing and severe lung pathology in both domestic and wild ruminants. In this study we investigated the interrelationships of Dictyocaulus spp. from European bison (Bison bonasus L.), roe deer (Capreolus capreolus), and red deer (Cervus elaphus) by nucleotide sequence analysis spanning the 18S RNA gene (small subunit [SSU]) and internal transcribed spacer 2 (ITS2) regions of the ribosomal gene array as well as the mitochondrial cytochrome c oxidase subunit 1 (cox1). Molecular analyses of sequence data obtained partly with novel primers from between 10 and 50 specimens from each host were carried out. Bayesian inference analysis revealed that each host species was infected with different genotypes. Analysis of cox1 sequence data showed a diverse genetic background and high evolutionary potential of Dictyocaulus taxa. Data from lungworms of European bison revealed a distinct genotype of Dictyocaulus viviparus, whereas Dictyocaulus capreolus was only found in roe deer. In contrast, red deer were infected with a taxon with unique SSU, ITS2, and cox1 sequences. These results indicate the occurrence of a novel genotype from red deer, which differs significantly from the National Center for Biotechnology Information reference sequence of Dictyocaulus eckerti. The molecular evidence was consistent with a morphological study with description and imaging of Dictyocaulus cervi n. sp. recovered from red deer. Dictyocaulus cervi n. sp. can be distinguished from D. eckerti on the basis of the absence of cervical papillae, the occurrence of a single ring of 4 symmetrical submedian cephalic papillae, length of the tail in females, morphometry of the female reproductive system, and measurements of gubernacula in males. In conclusion, our findings further strengthen the idea that the genetic complexity and diversity among Dictyocaulus lungworms infecting wildlife ruminants is larger than previously believed and warrants further investigation.

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.

Absence of genetic structure in Baylisascaris schroederi populations, a giant panda parasite, determined by mitochondrial sequencing

Background

Infection with the parasitic nematode, Baylisascaris schroederi (Ascaridida: Nematoda), is one of the most important causes of death in giant pandas, and was responsible for half of deaths between 2001 and 2005. Mitochondrial (mt) DNA sequences of parasites can unveil their genetic diversity and depict their likely dynamic evolution and therefore may provide insights into parasite survival and responses to host changes, as well as parasite control.

Methods

Based on previous studies, the present study further annotated the genetic variability and structure of B. schroederi populations by combining two different mtDNA markers, ATPase subunit 6 (atp6) and cytochrome c oxidase subunit I (cox1). Both sequences were completely amplified and genetically analyzed among 57 B. schroederi isolates, which were individually collected from ten geographical regions located in three important giant panda habitats in China (Minshan, Qionglai and Qinling mountain ranges).

Results

For the DNA dataset, we identified 20 haplotypes of atp6, 24 haplotypes of cox1, and 39 haplotypes of atp6 + cox1. Further haplotype network and phylogenetic analyses demonstrated that B. schroederi populations were predominantly driven by three common haplotypes, atp6 A1, cox1 C10, and atp6 + cox1 H11. However, due to low rates of gene differentiation between the three populations, both the atp6 and cox1 genes appeared not to be significantly associated with geographical divisions. In addition, high gene flow was detected among the B. schroederi populations, consistent with previous studies, suggesting that this parasite may be essentially homogenous across endemic areas. Finally, neutrality tests and mismatch analysis indicated that B. schroederi had undergone earlier demographic expansion.

Conclusions

These results confirmed that B. schroederi populations do not follow a pattern of isolation by distance, further revealing the possible existence of physical connections before geographic separation. This study should also contribute to an improved understanding of the population genetics and evolutionary biology of B. schroederi and assist in the control of baylisascariasis in giant pandas.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-014-0606-3) contains supplementary material, which is available to authorized users.

Unexpected absence of genetic separation of a highly diverse population of hookworms from geographically isolated hosts

The high natal site fidelity of endangered Australian sea lions (Neophoca cinerea) along the southern Australian coast suggests that their maternally transmitted parasitic species, such as hookworms, will have restricted potential for dispersal. If this is the case, we would expect to find a hookworm haplotype structure corresponding to that of the host mtDNA haplotype structure; that is, restricted among geographically separated colonies. In this study, we used a fragment of the cytochrome c oxidase I mitochondrial DNA (mtDNA) gene to investigate the diversity of hookworms (Uncinaria sanguinis) in N. cinerea to assess the importance of host distribution and ecology on the evolutionary history of the parasite. High haplotype (h=0.986) and nucleotide diversity (π=0.013) were seen, with 45 unique hookworm mtDNA haplotypes across N. cinerea colonies; with most of the variation (78%) arising from variability within hookworms from individual colonies. This is supported by the low genetic differentiation co-efficient (GST=0.007) and a high gene flow (Nm=35.25) indicating a high migration rate between the populations of hookworms. The haplotype network demonstrated no clear distribution and delineation of haplotypes according to geographical location. Our data rejects the vicariance hypothesis; that female host natal site fidelity and the transmammary route of infection restrict hookworm gene flow between N. cinerea populations and highlights the value of studies of parasite diversity and dispersal to challenge our understanding of parasite and host ecology.

The mitochondrial genome of Protostrongylus rufescens - implications for population and systematic studies

Background

Protostrongylus rufescens is a metastrongyloid nematode of small ruminants, such as sheep and goats, causing protostrongylosis. In spite of its importance, the ecology and epidemiology of this parasite are not entirely understood. In addition, genetic data are scant for P. rufescens and related metastrongyloids.

Methods

The mt genome was amplified from a single adult worm of P. rufescens (from sheep) by long-PCR, sequenced using 454-technology and annotated using bioinformatic tools. Amino acid sequences inferred from individual genes of the mt genomes were concatenated and subjected to phylogenetic analysis using Bayesian inference.

Results

The circular mitochondrial genome was 13,619 bp in length and contained two ribosomal RNA, 12 protein-coding and 22 transfer RNA genes, consistent with nematodes of the order Strongylida for which mt genomes have been determined. Phylogenetic analysis of the concatenated amino acid sequence data for the 12 mt proteins showed that P. rufescens was closely related to Aelurostrongylus abstrusus, Angiostrongylus vasorum, Angiostrongylus cantonensis and Angiostrongylus costaricensis.

Conclusions

The mt genome determined herein provides a source of markers for future investigations of P. rufescens. Molecular tools, employing such mt markers, are likely to find applicability in studies of the population biology of this parasite and the systematics of lungworms.

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.

The mitochondrial genome of Aelurostrongylus abstrusus-diagnostic, epidemiological and systematic implications

Aelurostrongylus abstrusus (Railliet, 1898) is a metastrongylid nematode of major clinical relevance in felids, causing aelurostrongylosis. In spite of its clinical importance in cats, the genetics, epidemiology and biology of this parasite are not entirely understood. mt DNA can provide markers for studies of these areas, but genetic data are scant for A. abstrusus and related lungworms. Here, the mt genome was amplified by long-range polymerase chain reaction (long-PCR) from a single male adult of A. abstrusus, sequenced using 454 technology and annotated using an established bioinformatic pipeline. This circular mt genome is 13,913 bp and contains two ribosomal RNA, 12 protein-coding and 22 transfer RNA genes, consistent with most other chromadorean nematodes. This genome should provide a source of markers for future investigations of the epidemiology and ecology of A. abstrusus. Molecular tools, employing such mt markers, are likely to find utility for explorations into the epidemiology, biology and systematics of this parasite, and the diagnosis of feline aelurostrongylosis.

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.

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.

Using data-display networks for exploratory data analysis in phylogenetic studies

Exploratory data analysis (EDA) is a frequently undervalued part of data analysis in biology. It involves evaluating the characteristics of the data "before" proceeding to the definitive analysis in relation to the scientific question at hand. For phylogenetic analyses, a useful tool for EDA is a data-display network. This type of network is designed to display any character (or tree) conflict in a data set, without prior assumptions about the causes of those conflicts. The conflicts might be caused by 1) methodological issues in data collection or analysis, 2) homoplasy, or 3) horizontal gene flow of some sort. Here, I explore 13 published data sets using splits networks, as examples of using data-display networks for EDA. In each case, I performed an original EDA on the data provided, to highlight the aspects of the resulting network that will be important for an interpretation of the phylogeny. In each case, there is at least one important point (possibly missed by the original authors) that might affect the phylogenetic analysis. I conclude that EDA should play a greater role in phylogenetic analyses than it has done.

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.

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.