Population biology of parasitic nematodes: applications of genetic markers

Haemonchus contortus, an obligatory haematophagus worm infection in small ruminants: Population genetics and genetic diversity

Haemonchus contortus, a stomach worm, is prevalent in ruminants worldwide. They particularly hamper profitable small ruminant production. Here, we estimate the genetic variation of H. contortus collected from slaughtered goats and sheep from various geographic zones of Bangladesh using multiple genes. To perform this, adult parasites were isolated from the abomasum of slaughtered animals (sheep and goats). Among them, 79 male H. contortus were identified by microscopy. Following the extraction of DNA, ITS-2 and cox1 genes were amplified and subsequently considered for sequencing. After alignment and editing, sequences were analyzed to find out sequence variation, diversity pattern of genes, and population genetics of isolates. Among the sequence data, the analyses identified 19 genotypes of ITS-2 and 77 haplotypes of cox1 genes. The diversity of nucleotides was 0.0103 for ITS-2 and 0.029 for cox1 gene. The dendogram constructed by the genotype and haplotype sequences of H. contortus revealed that two populations were circulating in Bangladesh without any demarcation of host and geographic regions. Analysis of population genetics demonstrated a high flow of genes (89.2 %) within the population of the worm in Bangladesh. The Fst value showed very little amount of genetic difference among the worm populations of Bangladesh but marked genetic variation between different continents. The findings are expected to help explain the risks of anthelmintic resistance and the transmission pattern of the parasite, and also provide a control strategy against H. contortus.

Comparison of the first time detected Oesophagostomum asperum with Oesophagostomum columbianum in sheep and goats in Bangladesh based on the trinity: Morphology, morphometry and genetic diversity

Oesophagostomum spp. (Family: Chabertiidae) is keeping a low profile in terms of severity in Bangladesh while maintaining economic loss through disguise within sheep and goats. The study was performed to identify prevalence, confirmation of species through morphology and morphometry followed by phylogeny using ITS2 and COX1 genes. In total 384 slaughterhouse-sourced small and large intestines were pooled from Mymensingh, Kishoreganj, Netrokona, Sherpur and Tangail districts of Mymensingh division. Followed by isolation, O. columbianum and O. asperum were identified following their key morphological features. Notably, O. asperum was first time detected in Bangladesh. The overall prevalence of Oesophagostomum spp. was found 60.93%. The prevalence of O. columbianum (64.95%) was almost double than that of O. asperum (35.04%). Among several characters, only the distance between anus to tail tip showed a significant morphological disparity in female. The Neighbor-joining (NJ) phylogenic trees based on ITS2 and COX1 genes confirmed the study species. The first time identified O. asperum along with morphometry and phylogeny will add value to the fact that nematodes are invisibly present with high prevalence in this country. This study will help to draw specific attention to command a practical control strategy for intervening in economic loss.

Genetic variability of Haemonchus contortus isolates in small ruminants from slaughterhouses in Bangladesh

Haemonchus contortus is a blood-sucking gastrointestinal nematode that infects all ruminants and causes significant economic losses in production. Characterizing the genetic variability of H. contortus populations is crucial for understanding patterns of disease transmission and developing effective control strategies against haemonchosis. This study aimed to identify the genetic variability of H. contortus isolates in small ruminants from slaughterhouses in Bangladesh. During January to December 2015, 400 abomasa samples were collected and 186 were found to be positive for Haemonchus. A 321-bp fragment of the second internal transcribed spacer (ITS-2) of nuclear ribosomal DNA and an 800-bp fragment of the mitochondrial nicotinamide dehydrogenase subunit-4 gene (nad4) were amplified using polymerase chain reaction (PCR) and directly sequenced. The results showed 10 genotypes (ITS-2) and 45 haplotypes (nad4) among the 186 worms. The sequences were 98.5 to 100% identical to reference sequences from the GenBank database. ITS-2 sequence analysis revealed four nucleotide substitutions at positions 30, 41, 42, and 216. There was one transition (C/T) at position 42 and three transversions (C/A at position 30, G/C at position 41, and T/A at position 216). The nad4 gene sequences showed 15 substitutions, all of which were transitions. The pairwise distance of ITS-2 between H. contortus populations ranged from 0.005 to 1.477. The nucleotide diversity (μ) among the populations was 0.009524 using ITS-2 and 0.00394 using nad4. This study indicated low genetic deviation among H. contortus populations in Bangladesh.

What lies behind the curtain: Cryptic diversity in helminth parasites of human and veterinary importance

Parasite cryptic species are morphologically indistinguishable but genetically distinct organisms, leading to taxa with unclear species boundaries. Speciation mechanisms such as cospeciation, host colonization, taxon pulse, and oscillation may lead to the emergence of cryptic species, influencing host-parasite interactions, parasite ecology, distribution, and biodiversity. The study of cryptic species diversity in helminth parasites of human and veterinary importance has gained relevance, since their distribution may affect clinical and epidemiological features such as pathogenicity, virulence, drug resistance and susceptibility, mortality, and morbidity, ultimately affecting patient management, course, and outcome of treatment. At the same time, the need for recognition of cryptic species diversity has implied a transition from morphological to molecular diagnostic methods, which are becoming more available and accessible in parasitology. Here, we discuss the general approaches for cryptic species delineation and summarize some examples found in nematodes, trematodes and cestodes of medical and veterinary importance, along with the clinical implications of their taxonomic status. Lastly, we highlight the need for the correct interpretation of molecular information, and the correct use of definitions when reporting or describing new cryptic species in parasitology, since molecular and morphological data should be integrated whenever possible.

Anthelmintic resistance: is a solution possible?

More than 50 years after anthelmintic resistance was first identified, its prevalence and impact on the animal production industry continues to increase across the world. The term "anthelmintic resistance" (AR) can be briefly defined as the reduction in efficacy of a certain dose of anthelmintic drugs (AH) in eliminating the presence of a parasite population that was previously susceptible. The main aim of this study is to examine anthelmintic resistance in domestic herbivores. There are numerous factors playing a role in the development of AR, but the most important is livestock management. The price of AH and the need to treat a high number of animals mean that farmers face significant costs in this regard, yet, since 1981, little progress has been made in the discovery of new molecules and the time and cost required to bring a new AH to market has increased dramatically in recent decades. Furthermore, resistance has also emerged for new AH, such as monepantel or derquantel. Consequently, ruminant parasitism cannot be controlled solely by using synthetic chemicals. A change in approach is needed, using a range of preventive measures in order to achieve a sustainable control programme. The use of nematophagous fungi or of plant extracts rich in compounds with anthelmintic properties, such as terpenes, condensed tannins, or flavonoids, represent potential alternatives. Nevertheless, although new approaches are showing promising results, there is still much to do. More research focused on the control of AR is needed.

Ventral dermatitis in rowi (Apteryx rowi) caused by cutaneous capillariasis

In 2013 there was an outbreak of crusting ventral dermatitis among a group of juvenile rowi (Apteryx rowi), a species of the endangered New Zealand kiwi, that were being raised on an off-shore island sanctuary. Biopsies taken at the time found nematodes migrating within the epidermis of affected skin but the specific identity and origin of the organisms was not established, and sporadic cases of similar skin disease continue to occur on the island. On examination of additional sections from the original skin biopsies, adult nematodes and eggs were identified, the histomorphology of which was consistent with Capillaria sensu lato. PCR was performed on DNA extracted from archived formalin-fixed, paraffin-embedded tissue blocks of skin from eight affected rowi, using primers targeting the 18S region of nuclear ribosomal DNA and the COI gene of mitochondrial DNA of capillarid nematodes. The 18S sequences from all rowi samples were identical and matched sequences from members of the genus Eucoleus. In contrast, two distinct capillarid COI sequences were obtained, in one case both from the same rowi skin biopsy. While there were no close matches, both COI sequences also aligned nearest to sequences identified as Eucoleus spp. It is considered unlikely that two different nematode species are involved in the rowi skin lesions and the possible amplification of a COI pseudogene or “numt” is discussed. A species-level identification of the capillarid nematodes causing skin disease in rowi was not obtained, however based on histological evaluation the infections include reproductively-active adult nematodes. This finding indicates the possibility of perpetuation of the skin disease in the absence of the original source, as well as raising potential for the transfer of infection from the island when the juvenile rowi are translocated to their new habitats.

•

Dermatitis associated with nematode migration has been diagnosed in juvenile rowi kiwi.

•

Histological and molecular analysis indicates a member of Capillaria sensu lato.

•

The possible amplification of a COI pseudogene is discussed.

Genetic Diversity of the root-knot nematode Meloidogyne enterolobii in Mulberry Based on the Mitochondrial COI Gene

This study explores the genetic diversity and structure of Meloidogyne enterolobii in mulberry in China. The COI mitochondrial gene (mtCOI) in M.enterolobii populations in Guangdong, Guangxi, and Hunan Provinces was PCR-amplified, sequenced, and analyzed for genetic diversity. The total number of variations, haplotypes (Hap), the average number of nucleotide differences (k), haplotype diversity (H), and nucleotide diversity (π) of mtCOI were 25, 11, 4.248, 0.900, and 0.00596, respectively. Insignificant differences in Fst value (0.0169) and a high level of gene flow (7.02) were detected among the 19-mulberry root-knot nematode populations, and high genetic variation within each population and a small genetic distance among populations were observed. Both phylogenetic analyses and network mapping of the 11 haplotypes revealed a dispersed distribution pattern of 19 mulberry root-knot nematode populations and an absence of branches strictly corresponding to the 19 range sampling sites. The neutrality test and mismatch analysis indicated that mulberry root-knot nematode populations experienced a population expansion in the past. The analysis of molecular variance (AMOVA) revealed that the genetic differentiation of M. enterolobii was mainly contributed by the variation within each group. No significant correlation was found between the genetic distance and geographical distance of M. enterolobii populations. The findings of this study provide a profound understanding of the M. enterolobii population and will inform the development of strategies to combat and manage root-knot nematodes in mulberry.

Phylogenetic relationships among Toxocara spp. and Toxascaris sp. from different regions of the world

Toxocara and Toxascaris are parasitic nematodes that infect canids and felids although species of the genus Toxocara also infect humans. This work aimed to establish the phylogenetic and phylogeographic relationship between specimens of T. canis, T. cati, T. malaysiensis and Toxascaris leonina and to evaluate the degree of host specificity. In total, 437 samples (adults and pools of eggs) were collected from canids and felids from eight countries. Parasites were identified by morphology, PCR linked Restriction Fragment Length Polymorphism (PCR-RFLP) and partial sequencing of the mitochondrial gene cox1. Phylogenetic trees were constructed and genetic distance among isolates was estimated. Based on the molecular characterization all worms were identified in agreement with their respective hosts with the exception of three samples; two from cats and one from dogs identified as T. canis and T. cati, respectively. There was no clear geographical clustering of the samples despite this study including parasites from three continents. This is the first study, to our knowledge, to use molecular methods to identify T. canis in cats and T. cati in dogs with host specificity being the most common finding. Our developed PCR-RFLP method was found to be a facile and reliable method for identifying Toxocara species.

First record of Hysterothylacium fabri (Rudolphi 1819) Deardorff and Overstreet 1980 from Scomber Colias of the South Atlantic waters

Among several marine pelagic species of the Brazilian coast, Scomber colias Gmelin, 1789 (Perciformes: Scombridae) stands out for having great economic importance, since it is widely used as a food resource and presents moderate vulnerability. Twenty specimens of S. colias were purchased from October 2015 to October 2016 from the coast of Santa Catarina, Brazil. In the present study, we recorded Hysterothylacium fabri (Rudolphi, 1819) (Deardorff and Overstreet, Proc Biol Soc Wash 93(4):1035-1079 1980) from the S. colias intestine using an integrative taxonomy approach, where morphological data are used in combination with partial sequences of the ITS gene, to validate the taxonomic status of the species and establish their relationships with other members of the genus. This species is being recorded for the first time in the South Atlantic and S. colias. The specimens of H. fabri collected in this study parasitizing S. colias presented morphology similar to the other specimens already registered parasitizing other hosts. The distance matrix generated showed that the partial sequences obtained in this study were more similar to sequences of Hysterothylacium sp. collected in China. In phylogenetic analysis, the two detected haplotypes of this study were grouped with H. fabri haplotypes deposited in GenBank in a monophyletic subclade.

Metabarcoding of soil nematodes: the importance of taxonomic coverage and availability of reference sequences in choosing suitable marker(s)

For many organisms, there is agreement on the specific genomic region used for developing barcode markers. With nematodes, however, it has been found that the COI region designated for most animals lacks the taxonomic coverage (ability to amplify a diverse group of taxa) required of a metabarcoding marker. For that reason, studies on metabarcoding of nematodes thus far have utilized primarily regions within the highly conserved 18S ribosomal DNA. Two popular markers within this region are the ones flanked by the primer pairs NF1-18Sr2b and SSUF04-SSUR22. The NF1-18Sr2b primer pair, especially, has been critiqued as not being specific enough for nematodes leading to suggestions for other candidate markers while the SSUF04-SSUR22 region has hardly been tested on soil nematodes. The current study aimed to evaluate these two markers against other alternative ones within the 28S rDNA and the COI region for their suitability for nematode metabarcoding. The results showed that the NF1-18Sr2b marker could offer wide coverage and good resolution for characterizing soil nematodes. Sufficient availability of reference sequences for this region was found to be a significant factor that resulted in this marker outperforming the other markers, particularly the 18S-based SSUFO4-SSUR22 marker. None of the other tested regions compared with this marker in terms of the proportion of the taxa recovered. The COI-based marker had the lowest number of taxa recovered, and this was due to the poor performance of its primers and the insufficient number of reference sequences in public databases. In summary, this study highlights how dependent the success of metabarcoding is on the availability of a good reference sequence collection for the marker of choice as well as its taxonomic coverage.

A detection of benzimidazole resistance-associated SNPs in the isotype-1 β-tubulin gene in Haemonchus contortus from wild blue sheep (Pseudois nayaur) sympatric with sheep in Helan Mountains, China

Background

Benzimidazole (BZ) resistance is an increasingly serious problem due to the excessive use of this anthelmintic for controlling Haemonchus contortus, which is one of the major gastrointestinal nematodes infecting small ruminants worldwide. Three known single nucleotide polymorphisms (SNPs), F167Y (TAC), E198A (GCA) and F200Y (TAC), in the isotype-1 β-tubulin gene of H. contortus are associated with BZ resistance. Comprehending the spread and origins of BZ resistance-associated SNPs has important implications for the control of this nematode.

Results

Twenty-seven adult H. contortus were harvested from wild blue sheep (Pseudois nayaur), small wild ruminants sympatric with domestic ruminants, inhabiting the Helan Mountains, China, to monitor the status of BZ resistance. In addition, 20 adult H. contortus from domestic sheep sympatric with this wild ruminant and 36 isotype-1 β-tubulin haplotype sequences of H. contortus (two of these haplotypes, E198A3 and E198A4, possessed resistance-associated SNP E198A (GCA) from domestic ruminants in eight other geographical regions of China were used to further define the origins of BZ resistance-associated SNPs within the worms collected from blue sheep. The BZ resistance-associated SNP E198A was detected, whereas SNPs F167Y (TAC) and F200Y (TAC) were not found within the worms collected from blue sheep, and the frequency of homozygous resistant E198A (GCA) was 7.40%. The evolutionary tree and network showed consistent topologies for which there was no obvious boundary among the worms from the wild and domestic hosts, and two haplotypes (E198A1 and E198A2) possessing E198A from the wild blue sheep had two different independent origins. E198A1 had the same origin with E198A3 but E198A2 had a different origin with them. Population genetic analyses revealed a low level of Fst values (ranging from 0 to 0.19749) between all H. contortus worm groups in China.

Conclusions

Results of the current study of the three BZ resistance-associated SNPs of H. contortus from wild blue sheep suggested that only E198A (GCA) was present within the worms collected from the wild ruminants and had multiple independent origins.

Electronic supplementary material

The online version of this article (10.1186/s12917-019-1838-4) contains supplementary material, which is available to authorized users.

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.

Genetic Diversity of Echinococcus granulosus Genotype G1 in Xinjiang, Northwest of China

Cystic echinococcosis (CE) caused by E. granulosus is a serious helminthic zoonosis in humans, livestock and wildlife. Xinjiang is one of high endemic province for CE in China. A total of 55 sheep and cattle livers containing echinococcal cysts were collected from slaughterhouses in Changji and Yining City, northern region of Xinjiang. PCR was employed for cloning 2 gene fragments, 12S rRNA and CO1 for analysis of phylogenetic diversity of E. granulosus. The results showed that all the samples collected were identified as G1 genotype of E. granulosus. Interestingly, YL5 and CJ75 strains were the older branches compared to those strains from France, Argentina, Australia. CO1 gene fragment showed 20 new genotype haploids and 5 new genotype haplogroups (H1-H5) by the analysis of Network 5.0 software, and the YLY17 strain was identified as the most ancestral haplotype. The major haplotypes, such as CJ75 and YL5 strains, showed identical to the isolates from Middle East. The international and domestic trade of livestock might contribute to the dispersal of different haplotypes for E. granulosus evolution.

First morphological and molecular analysis of Eucoleus boehmi like eggs in dogs from Argentina

The canid parasites Eucoleus aerophilus (syn. Capillaria aerophila) and Eucoleus boehmi (syn. Capillaria boehmi) parasitize the lower and the upper respiratory tract, respectively. Reports and descriptions of these nematodes are scarce in Argentina, possibly due to misdiagnosis of morphologically similar trichuroids eggs, and the lack of knowledge about the species of Eucoleus in this geographical area. Scanning electron microscopy is a useful tool for identification of E. boehmi eggs based on the characteristics of the shell structure which differentiate between species. Molecular analysis complements morphological identification. Until now, there are no studies based on the analysis of E. boehmi eggs in Argentina. The aim of the present work was to study by morphological, morphometric, and molecular analysis, eggs attributable to E. boehmi isolated from dogs naturally infected in Mar del Plata city, Argentina. Eggs isolated from two dog fecal samples were analyzed by light and scanning electron microscopy. A fragment of mitochondrial DNA (mtDNA) of the cytochrome c oxidase subunit I gene (cox1) from eggs was sequenced, and phylogenetic analysis was performed in this study. According to morphological results based on the wall surface ultrastructure, the eggs studied were assigned to E. boehmi. Molecular analysis supported the morphological identification. The divergence of 9-12% with the European isolated could suggest a new geographical genetic variation of E. boehmi, but also question the possible existence of cryptic species. This is the first characterization of E. boehmi eggs in dogs from Argentina.

Genetic polymorphism of Baylisascaris procyonis in host infrapopulations and component populations in the Central USA

Baylisascaris procyonis is a nematode of significant concern to public and domestic animal health as well as wildlife management. The population genetics of B. procyonis is poorly understood. To gain insights into patterns of genetic diversity within (infrapopulation level) and among (component population level) raccoon (Procyon lotor) hosts, and specifically to assess the relative importance of indirect and direct transmission of the parasite for explaining observed population structure, we collected 69 B. procyonis from 17 wild raccoons inhabiting five counties in Missouri and Arkansas, USA. Informative regions of mitochondrial (CO1, CO2) and nuclear (28S, ITS2) genes were amplified and the distribution and genetic variability of these genes were assessed within and across raccoons. Concatenation of the CO1 and CO2 mtDNA sequences resulted in 5 unique haplotypes, with haplotype diversity 0.456 ± 0.068. The most common haplotype occurred in 94% of raccoons and 72.5% of B. procyonis. Sequences for 28S rDNA revealed four unique nuclear genotypes, the most common found in 100% of raccoons and 82.6% of B. procyonis. ITS2 genotypes were assessed using fragment analysis, and there was a 1:1 correspondence between 28S and ITS-2 genotypes. Infrapopulation variation in haplotypes and genotypes was high and virtually all hosts infected with multiple sequenced nematodes also harbored multiple haplotypes and genotypes. There was a positive relationship between the size of the analyzed infrapopulation (i.e., the number of nematodes analyzed) and the number of haplotypes identified in an individual. Collectively this work emphasizes the importance of indirect transmission in the lifecycle to this parasite.

Combination deworming for the control of double-resistant cyathostomin parasites - short and long term consequences

Equine cyathostomin are pervasive gastrointestinal parasites with wide-spread resistance to the benzimidazole and tetrahydropyrimidine drug classes worldwide. Combination deworming has been proposed as a more sustainable parasite control strategy. Simulation studies have found combination deworming to be effective in controlling drug resistant ovine trichostrongylid parasites. One equine study demonstrated an additive effect of a combination of oxibendazole and pyrantel pamoate against cyathostomins. However, this is the only equine study evaluating combination therapy, and the effects of repeated combination treatments administered over time remain unknown. The purpose of this study was to observe the efficacy of repeated oxibendazole/pyrantel pamoate combination therapy administered over one year against a cyathostomin population with resistance to benzimidazole and pyrantel products. Fecal egg counts were determined for the entire herd (N = 21) at the day of anthelmintic treatment and at two-week intervals for eight weeks post treatment. Starting efficacies of oxibendazole (OBZ, 10 mg/kg) and pyrantel pamoate (PYR, 6.6 mg base/kg) were 66.7% and 63.3%, respectively. Hereafter, the herd was treated four times with an oxibendazole/pyrantel pamoate combination, eight weeks apart, followed by repeating the single active treatments before concluding the study. While the first combination treatment exhibited an additive effect of the two active ingredients, this efficacy was not sustained over the course of the study. Mean fecal egg count reduction (FECR) was significantly greater for the first combination treatment (76.6%) than the second (42.6%, p = 0.0454), third (41.6%, p = 0.0318), and fourth (40.7%, p = 0.0372) combination treatments. The final single active mean FECRs were 42.3% for oxibendazole, and 42.7% for pyrantel pamoate. These efficacies were not significantly different from the initial single active efficacies (OBZ, p = 0.4421; PYR, p = 0.8361). These results suggest that combination therapy against double resistant equine cyathostomin populations is not sustainable, when using actives with markedly decreased starting efficacies.

Genetic diversity of Haemonchus contortus isolated from sympatric wild blue sheep (Pseudois nayaur) and sheep in Helan Mountains, China

Background

Haemonchus contortus is known among parasitic nematodes as one of the major veterinary pathogens of small ruminants and results in great economic losses worldwide. Human activities, such as the sympatric grazing of wild with domestic animals, may place susceptible wildlife hosts at risk of increased prevalence and infection intensity with this common small ruminant parasite. Studies on phylogenetic factors of H. contortus should assist in defining the amount of the impact of anthropogenic factors on the extent of sharing of agents such as this nematode between domestic animals and wildlife.

Methods

H. contortus specimens (n = 57) were isolated from wild blue sheep (Pseudois nayaur) inhabiting Helan Mountains (HM), China and additional H. contortus specimens (n = 20) were isolated from domestic sheep that were grazed near the natural habitat of the blue sheep. Complete ITS2 (second internal transcribed spacer) sequences and partial sequences of the nad4 (nicotinamide dehydrogenase subunit 4 gene) gene were amplified to determine the sequence variations and population genetic diversities between these two populations. Also, 142 nad4 haplotype sequences of H. contortus from seven other geographical regions of China were retrieved from database to further examine the H. contortus population structure.

Results

Sequence analysis revealed 10 genotypes (ITS2) and 73 haplotypes (nad4) among the 77 specimens, with nucleotide diversities of 0.007 and 0.021, respectively, similar to previous studies in other countries, such as Pakistan, Malaysia and Yemen. Phylogenetic analyses (BI, MP, NJ) of nad4 sequences showed that there were no noticeable boundaries among H. contortus populations from different geographical origin and population genetic analyses revealed that most of the variation (94.21%) occurred within H. contortus populations. All phylogenetic analyses indicated that there was little genetic differentiation but a high degree of gene flow among the H. contortus populations among wild blue sheep and domestic ruminants in China.

Conclusions

The current work is the first genetic characterization of H. contortus isolated from wild blue sheep in the Helan Mountains region. The results revealed a low genetic differentiation and high degree of gene flow between the H. contortus populations from sympatric wild blue sheep and domestic sheep, indicating regular cross-infection between the sympatrically reared ruminants.

Electronic supplementary material

The online version of this article (10.1186/s13071-017-2377-0) contains supplementary material, which is available to authorized users.

Genomic introgression mapping of field-derived multiple-anthelmintic resistance in Teladorsagia circumcincta

Preventive chemotherapy has long been practiced against nematode parasites of livestock, leading to widespread drug resistance, and is increasingly being adopted for eradication of human parasitic nematodes even though it is similarly likely to lead to drug resistance. Given that the genetic architecture of resistance is poorly understood for any nematode, we have analyzed multidrug resistant Teladorsagia circumcincta, a major parasite of sheep, as a model for analysis of resistance selection. We introgressed a field-derived multiresistant genotype into a partially inbred susceptible genetic background (through repeated backcrossing and drug selection) and performed genome-wide scans in the backcross progeny and drug-selected F2 populations to identify the major genes responsible for the multidrug resistance. We identified variation linking candidate resistance genes to each drug class. Putative mechanisms included target site polymorphism, changes in likely regulatory regions and copy number variation in efflux transporters. This work elucidates the genetic architecture of multiple anthelmintic resistance in a parasitic nematode for the first time and establishes a framework for future studies of anthelmintic resistance in nematode parasites of humans.

Teladorsagia circumcincta is an economically significant nematode (roundworm) pathogen affecting sheep and goats in temperate regions of the world. The widespread use of prophylactic treatment has resulted in rapid selection for anthelmintic (anti-worm drug) resistance in this and other species of livestock parasites. The mechanism of resistance is not well understood because most studies have focused on the role of candidate genes using simplistic models of single gene selection, despite evidence that the evolution of resistance is more complex. Here, we report on a comprehensive whole-genome analysis that elucidated resistance-associated genes, which was facilitated by developing a pair of T. circumcincta strains sharing a largely common genetic background but differing markedly in their susceptibility to anthelmintic drugs. The results show that multiple genetic factors contribute to anthelmintic resistance in a variety of ways, including possible reduction/modulation in target site sensitivity, reduced target site expression, and increased drug efflux, to name a few. This suggests that drug resistance in these parasites is a multifactorial quantitative trait rather than a simple discrete Mendelian character. With this study, we established a genomics-based experimental paradigm for investigating anthelmintic resistance, at a time when its medical importance is rapidly increasing.

Population structure of Haemonchus contortus from seven geographical regions in China, determined on the basis of microsatellite markers

Background

Studying genetic variation within and among Haemonchus contortus populations can inform some aspects of this parasite’s population genetics and epidemiology. However, almost nothing is known about such variation in China.

Methods

Adult males of H. contortus (n = 184) representing seven distinct populations in China were collected, and genetic variation within and among these populations was explored using eight distinct microsatellite markers.

Results

Genetic parameters, such as heterozygosity and inbreeding coefficient (F_IS) indicated that all eight microsatellites were highly polymorphic. Various analyses (AMOVA, F_ST, phylogenetic, structure, mantel test and population dynamics) revealed high within-population variation, low population genetic differentiation and high gene flow for H. contortus in China.

Conclusions

This study provides a first snapshot of the genetic substructuring of H. contortus populations in China using polymorphic markers, and might provide a starting point for assessing genetic changes over space and time during or following the implementation of particular treatment or control strategies, or changes as a consequence of environmental, management and climatic factors.

Electronic supplementary material

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

Phylogenetic patterns of Haemonchus contortus and related trichostrongylid nematodes isolated from Egyptian sheep

Haemonchus contortus is a major parasite of small ruminants and its blood-feeding behaviour causes effects ranging from mild anaemia to death. Knowledge of the genetic variation within and among H. contortus populations can provide the foundation for understanding transmission patterns and aid in the control of haemonchosis. Adult male H. contortus were collected from three geographical regions in Egypt. The second internal transcribed spacer (ITS2) of nuclear ribosomal DNA was amplified using the polymerase chain reaction (PCR) and sequenced directly. The population genetic diversity and sequence variations were determined. Nucleotide sequence analyses revealed one genotype (ITS2) in all worms, without genetic differentiation. The similarity in population genetic diversity and genetic patterns observed among the three geographical regions could be attributed to possible movement between the sites. This is the first study of genetic variation in H. contortus in Egypt. The present results could have implications for the rapid characterization of H. contortus and other trichostrongyloid nematodes, and evaluation of the epidemiology of H. contortus in Egypt.

Harnessing the Toxocara Genome to Underpin Toxocariasis Research and New Interventions

Parasitic worms, such as flatworms (platyhelminths) and roundworms (nematodes), cause substantial morbidity and mortality in animals and people globally. The ascaridoid nematode Toxocara canis is a zoonotic parasite of socioeconomic significance worldwide. In humans, this worm causes toxocariasis (disease) mainly in underprivileged communities in both the developed and developing worlds. While reasonably well studied from clinical and epidemiological perspectives, little is understood about the molecular biology of T. canis, its relationship with its hosts and the disease that it causes. However, a recent report of the draft genome and transcriptomes of T. canis should underpin many fundamental and applied research areas in the future. The present article gives a background on Toxocara and toxocariasis, a brief account of diagnostic approaches for specific identification and genetic analysis, and gives a perspective on the impact that the genome of T. canis and advanced molecular technologies could have on our understanding of the parasite and the diseases that it causes as well as the design of new and improved approaches for the diagnosis, treatment and control of toxocariasis.

Fast assembly of the mitochondrial genome of a plant parasitic nematode (Meloidogyne graminicola) using next generation sequencing

Little is known about the variations of nematode mitogenomes (mtDNA). Sequencing a complete mtDNA using a PCR approach remains a challenge due to frequent genome reorganizations and low sequence similarities between divergent nematode lineages. Here, a genome skimming approach based on HiSeq sequencing (shotgun) was used to assemble de novo the first complete mtDNA sequence of a root-knot nematode (Meloidogyne graminicola). An AT-rich genome (84.3%) of 20,030 bp was obtained with a mean sequencing depth superior to 300. Thirty-six genes were identified with a semi-automated approach. A comparison with a gene map of the M. javanica mitochondrial genome indicates that the gene order is conserved within this nematode lineage. However, deep genome rearrangements were observed when comparing with other species of the superfamily Hoplolaimoidea. Repeat elements of 111 bp and 94 bp were found in a long non-coding region of 7.5 kb, as similarly reported in M. javanica and M. hapla. This study points out the power of next generation sequencing to produce complete mitochondrial genomes, even without a reference sequence, and possibly opening new avenues for species/race identification, phylogenetics and population genetics of nematodes.

Morphological and molecular characterization of Mermis nigrescens Dujardin, (Nematoda: Mermithidae) parasitizing the introduced European earwig (Dermaptera: Forficulidae) in New Zealand

Parasitic nematodes of the family Mermithidae were found to be infecting the introduced European earwig Forficula auricularia (Dermaptera: Forficulidae) in Dunedin, South Island, New Zealand. Adult females were later collected from various garden plants while depositing eggs. These mermithid specimens were identified morphologically as Mermis nigrescens Dujardin, 1842. A genetic distance of 0.7% between these specimens and a M. nigrescens isolate from Canada (18S rRNA gene), suggests that they have diverged genetically, but there are currently no available comparable sequences for the European M. nigrescens. Two additional nuclear fragments were also amplified, the 28S rRNA and the ribosomal DNA first internal transcribed spacer (ITS1), providing a basis for future studies. Bearing in mind the morphological similarity with other reported M. nigrescens and the lack of sequence data from other parts of the world, we retain the name M. nigrescens, and suggest that the species may be found to represent a complex of cryptic species when more worldwide data are available. Herein, we present a brief description of the post-parasitic worms and adult females, along with an inferred phylogeny using 18S rRNA gene sequences.

Advances in the diagnosis of key gastrointestinal nematode infections of livestock, with an emphasis on small ruminants

Parasitic nematodes (roundworms) of livestock have major economic impact globally. In spite of the diseases caused by these nematodes and some advances in the design of new therapeutic agents (anthelmintics) and attempts to develop vaccines against some of them, there has been limited progress in the establishment of practical diagnostic techniques. The specific and sensitive diagnosis of gastrointestinal nematode infections of livestock underpins effective disease control, which is highly relevant now that anthelmintic resistance (AR) is a major problem. Traditional diagnostic techniques have major constraints, in terms of sensitivity and specificity. The purpose of this article is to provide a brief background on gastrointestinal nematodes (Strongylida) of livestock and their control; to summarize conventional methods used for the diagnosis and discuss their constraints; to review key molecular-diagnostic methods and recent progress in the development of advanced amplification-based and sequencing technologies, and their implications for epidemiological investigations and the control of parasitic diseases.

Analyses of mitochondrial amino acid sequence datasets support the proposal that specimens of Hypodontus macropi from three species of macropodid hosts represent distinct species

BACKGROUND

Hypodontus macropi is a common intestinal nematode of a range of kangaroos and wallabies (macropodid marsupials). Based on previous multilocus enzyme electrophoresis (MEE) and nuclear ribosomal DNA sequence data sets, H. macropi has been proposed to be complex of species. To test this proposal using independent molecular data, we sequenced the whole mitochondrial (mt) genomes of individuals of H. macropi from three different species of hosts (Macropus robustus robustus, Thylogale billardierii and Macropus [Wallabia] bicolor) as well as that of Macropicola ocydromi (a related nematode), and undertook a comparative analysis of the amino acid sequence datasets derived from these genomes.

RESULTS

The mt genomes sequenced by next-generation (454) technology from H. macropi from the three host species varied from 13,634 bp to 13,699 bp in size. Pairwise comparisons of the amino acid sequences predicted from these three mt genomes revealed differences of 5.8% to 18%. Phylogenetic analysis of the amino acid sequence data sets using Bayesian Inference (BI) showed that H. macropi from the three different host species formed distinct, well-supported clades. In addition, sliding window analysis of the mt genomes defined variable regions for future population genetic studies of H. macropi in different macropodid hosts and geographical regions around Australia.

CONCLUSIONS

The present analyses of inferred mt protein sequence datasets clearly supported the hypothesis that H. macropi from M. robustus robustus, M. bicolor and T. billardierii represent distinct species.

Mitochondrial genomes of Anisakis simplex and Contracaecum osculatum (sensu stricto)--comparisons with selected nematodes

Anisakid nematodes parasitize mainly fish, marine mammals and/or fish-eating birds, and can be transmitted to a range of fish-eating mammals, including humans, where they can cause gastrointestinal disease linked to larval infection or allergic responses. In spite of the animal and human health significance of these parasites, there are still gaps in our understanding of the systematics, biology, epidemiology and ecology of anisakids. Mitochondrial (mt) DNA provides useful genetic markers for investigations in these areas, but complete mt genomic data have been lacking for most anisakids. In the present study, the mt genomes of Anisakis simplex sensu stricto and Contracaecum osculatum sensu stricto were amplified from genomic DNA by long-range polymerase chain reaction and sequenced using 454 technology. The circular mt genomes of these species were 13,926 and 13,823 bp, respectively, and each of them contained 12 protein-coding, 22 transfer RNA, and 2 ribosomal RNA genes consistent for members of the Ascaridida, Oxyurida, Spirurida, Rhabditida and Strongylida. These mt genomes provide a stepping-stone for future comparative analyses of a range of anisakids and a basis for reinvestigating their genetic relationships. In addition, these markers might be used in prospecting for cryptic species and exploring host affiliations.

Genetic variability within and among Haemonchus contortus isolates from goats and sheep in China

BACKGROUND

Haemonchus contortus (order Strongylida) is a common parasitic nematode infecting small ruminants and causing significant economic losses worldwide. Knowledge of genetic variation within and among H. contortus populations can provide a foundation for understanding transmission patterns, the spread of drug resistance alleles and might assist in the control of haemonchosis.

METHODS

152 H. contortus individual adult worms were collected from seven different geographical regions in China. The second internal transcribed spacer (ITS-2) of the nuclear ribosomal DNA and mitochondrial nicotinamide dehydrogenase subunit 4 gene (nad4) were amplified by polymerase chain reaction (PCR) and sequenced directly. The sequence variations and population genetic diversities were determined.

RESULTS

Nucleotide sequence analyses revealed 18 genotypes (ITS-2) and 142 haplotypes (nad4) among the 152 worms, with nucleotide diversities of 2.6% and 0.027, respectively, consistent with previous reports from other countries, including Australia, Brazil, Germany, Italy, Malaysia, Sweden, the USA and Yemen. Population genetic analyses revealed that 92.4% of nucleotide variation was partitioned within populations; there was no genetic differentiation but a high gene flow among Chinese populations; some degree of genetic differentiation was inferred between some specimens from China and those from other countries.

CONCLUSIONS

This is the first study of genetic variation within H. contortus in China. The results revealed high within-population variations, low genetic differentiation and high gene flow among different populations of H. contortus in China. The present results could have implications for studying the epidemiology and ecology of H. contortus in China.

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.

Analysis of the genetic diversity of the nematode parasite Baylisascaris schroederi from wild giant pandas in different mountain ranges in China

Background

Baylisascaris schroederi is one of the most common nematodes of the giant panda, and can cause severe baylisascarosis in both wild and captive giant pandas. Previous studies of the giant pandas indicated that this population is genetically distinct, implying the presence of a new subspecies. Based on the co-evolution between the parasite and the host, the aim of this study was to investigate the genetic differentiation in the B. schroederi population collected from giant pandas inhabiting different mountain ranges, and further to identify whether the evolution of this parasite correlates with the evolution of giant pandas.

Methods

In this study, 48 B. schroederi were collected from 28 wild giant pandas inhabiting the Qinling, Minshan and Qionglai mountain ranges in China. The complete sequence of the mitochondrial cytochrome b (mtCytb) gene was amplified by PCR, and the corresponding population genetic diversity of the three mountain populations was determined. In addition, we discussed the evolutionary relationship between B. schroederi and its host giant panda.

Results

For the DNA dataset, insignificant Fst values and a significant, high level of gene flow were detected among the three mountain populations of B. schroederi, and high genetic variation within populations and a low genetic distance were observed. Both phylogenetic analyses and network mapping of the 16 haplotypes revealed a dispersed pattern and an absence of branches strictly corresponding to the three mountain range sampling sites. Neutrality tests and mismatch analysis indicated that B. schroederi experienced a population expansion in the past.

Conclusions

Taken together, the dispersed haplotype map, extremely high gene flow among the three populations of B. schroederi, low genetic structure and rapid evolutionary rate suggest that the B. schroederi populations did not follow a pattern of isolation by distance, indicating the existence of physical connections before these populations became geographically separated.

Genetic variability of Baylisascaris schroederi from the Qinling subspecies of the giant panda in China revealed by sequences of three mitochondrial genes

The present study examined the variations in three mitochondrial (mt) DNA sequences, namely cytochrome b (cytb), cytochrome c oxidase subunit 3 (cox3) and NADH dehydrogenase subunit 5 (nad5), among Baylisascaris schroederi isolates from the Qinling subspecies of the giant panda in Shaanxi province, northwestern China. No differences in length were detected in the three mt fragments from different isolates. The intra-specific sequence variations within all B. schroederi samples were 0-2.6% for pcytb, 0-1.8% for pcox3 and 0-2.1% for pnad5, while the inter-specific sequence differences among members of the genus Baylisascaris were 8.2-15.2%, 6.2-15.9% and 8.4-16.0% for pcytb, pcox3, pnad5, respectively. A phylogenetic analysis of the combined sequences of pcytb, pcox3 and pnad 5 showed that all B. schroederi samples in the present study were located in two large clusters, with one cluster containing samples from giant pandas in Sichuan province. These findings provide basic information for further study of molecular epidemiology and control of B. schroederi infection in the Qinling subspecies of the giant panda and throughout China.

Trichuris spp. (Nematoda: Trichuridae) from two rodents, Mastomys natalensis and Gerbilliscus vicinus in Tanzania

During a survey of the helminth community of several rodent species in the Morogoro region (Tanzania), Trichuris whipworms (Nematoda: Trichuridae) were found in the ceca of the Natal multimammate mouse, Mastomys natalensis and a gerbil, Gerbilliscus vicinus (both Rodentia: Muridae). The taxonomic literature regarding Trichuris from African native rodents describes 10 species, but includes few metric and morphologic characters that discriminate between some of the pairs. The whipworms we sampled in Tanzanian Natal multimammate mice and gerbils were morphologically identified, respectively, as Trichuris mastomysi Verster, 1960 and Trichuris carlieri Gedoelst, 1916 sensu lato, but with characters that overlap or partially overlap with the cosmopolitan Murinae whipworm, Trichuris muris , already reported from several rodents in Africa. To clarify our identification, we sequenced the ITS-1, 5.8S, and ITS-2 ribosomal DNA region of the worms' nuclear genome. The genetic analyses clearly distinguish the whipworms we found in M. natalensis from those found in the gerbil, and both of these from T. muris whipworm reference sequences. The overlap of morphological characters between rodent whipworms suggests that reports of T. muris from rodent species not closely related to Murinae in other parts of Africa should be treated with caution.

Next-generation molecular-diagnostic tools for gastrointestinal nematodes of livestock, with an emphasis on small ruminants: a turning point?

Parasitic nematodes of livestock have major economic impact worldwide. Despite the diseases caused by these nematodes, some advances towards the development of new therapeutic agents and attempts to develop effective vaccines against some of them, there has been limited progress in the development of practical diagnostic methods. The specific and sensitive diagnosis of parasitic nematode infections of livestock underpins effective disease control, which is now particularly important given the problems associated with anthelmintic resistance in parasite populations. Traditional diagnostic methods have major limitations, in terms of sensitivity and specificity. This chapter provides an account of the significance of parasitic nematodes (order Strongylida), reviews conventional diagnostic techniques that are presently used routinely and describes advances in polymerase chain reaction (PCR)-based methods for the specific diagnosis of nematode infections. A particular emphasis is placed on the recent development of a robotic PCR-based platform for high-throughput diagnosis, and its significance and implications for epidemiological investigations and for use in control programmes.

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.

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.