Population genetics of Parascaris equorum based on DNA fingerprinting

The equine ascarids: resuscitating historic model organisms for modern purposes

The equine ascarids, Parascaris spp., are important nematode parasites of juvenile horses and were historically model organisms in the field of cell biology, leading to many important discoveries, and are used for the study of chromatin diminution. In veterinary parasitology, Parascaris spp. are important not only because they can cause clinical disease in young horses but also because they are the only ascarid parasites to have developed widespread anthelmintic resistance. Despite this, much of the general biology and mechanisms of anthelmintic resistance are poorly understood. This review condenses known basic biological information and knowledge on the mechanisms of anthelmintic resistance in Parascaris spp., highlighting the importance of foundational research programs. Although two variants of this parasite were recognized based on the number of chromosomes in the 1870s and suggested to be two species in 1890, one of these, P. univalens, appears to have been largely forgotten in the veterinary scientific literature over the past 100 years. We describe how this omission has had a century-long effect on nomenclature and data analysis in the field, highlighting the importance of proper specimen identification in public repositories. A summary of important basic biology, including life cycle, in vitro maintenance, and immunology, is given, and areas of future research for the improvement of knowledge and development of new systems are given. Finally, the limited knowledge regarding anthelmintic resistance in Parascaris spp. is summarized, along with caution regarding assumptions that resistance mechanisms can be applied across clades.

Equus roundworms (Parascaris univalens) are undergoing rapid divergence while genes involved in metabolic as well as anthelminic resistance are under positive selection

Background

The evolution of parasites is often directly affected by the host's environment. Studies on the evolution of the same parasites in different hosts are of great interest and are highly relevant to our understanding of divergence.

Methods

Here we performed whole-genome sequencing of Parascaris univalens from different Equus hosts (horses, zebras and donkeys). Phylogenetic and selection analyses were performed to study the divergence and adaptability of P. univalens.

Results

At the genetic level, multiple lines of evidence indicate that P. univalens is mainly separated into two clades (horse-derived and zebra & donkey-derived). This divergence began 300–1000 years ago, and we found that most of the key enzymes related to glycolysis were under strong positive selection in zebra & donkey-derived roundworms, whereas the lipid-related metabolic system was under positive selection in horse-derived roundworms, indicating that the adaptive evolution of metabolism has occurred over the past few centuries. In addition, we found that some drug-related genes showed a significantly higher degree of selection in diverse populations.

Conclusions

This work reports the adaptive evolution and divergence trend of P. univalens in different hosts for the first time. Its results indicate that the divergence of P. univalens is a continuous, dynamic process. Furthermore, the continuous monitoring of the effects of differences in nutritional and drug histories on the rapid evolution of roundworms is conducive to further understanding host-parasite interactions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08702-6.

Very low intraspecific sequence variation in selected nuclear and mitochondrial Parascaris univalens genes

Equines were over decades considered to be infected by two morphologically virtually indistinguishable ascarid species, Parascaris univalens and Parascaris equorum. Reliable species discrimination is only possible using enzyme isoelectric focussing and karyotyping with P. univalens having one and P. equorum two chromosome pairs. However, presumably the complexity of both methods prevented their routine use in nearly all previous studies about prevalence and drug resistance of Parascaris spp. These have barely been performed on the species level although most studies stated presence of one or the other species. Recently, only P. univalens has been identified by karyotyping and the last published study identifying P. equorum dates back to 1989. In order to improve species-specific detection, molecular markers are required. Here, partial 12S rRNA, cytochrome oxidase I (COI) and complete internal transcribed spacer (ITS)-1 and - 2 sequences were obtained from 24 karyotyped Parascaris specimens from Poland and 6 German specimens (not karyotyped) and used in phylogenetic analyses with orthologous sequences from GenBank. All karyotyped specimens were identified as P. univalens. In the phylogenetic analysis, they formed very homogenous clusters for all target genes and in a multi-locus analysis. Within this cluster, almost all sequences from GenBank were also included, no matter if they had been assigned to P. univalens or P. equorum. However, a small number of P. univalens ITS and COI sequences originating from donkeys from a single farm in China formed a highly supported sister cluster suggesting that they might represent another Parascaris genotype or species. Our data also strongly suggest that nearly all ITS and COI sequences previously deposited in GenBank and assigned to P. equorum actually represent P. univalens. The fact that significantly different sequences can be found in Parascaris spp. suggests that PCR-based species diagnosis will be possible once molecular markers have been identified for P. equorum from karyotyped specimens.

First Report of Resistance to Ivermectin in Parascaris univalens in Iceland

Horses in Iceland have been isolated for more than 1,000 yr but still harbor a similar range of gastrointestinal parasites as do horses across the world. The long isolation of the horses and their parasites presumably means that no resistance genes have been introduced into the Parascaris spp. population. It is therefore of particular interest to investigate the efficacy of ivermectin on Parascaris spp. infecting Icelandic foals. Potential treatment failure of ivermectin in Iceland will add substantial new information on how resistance can arise independently. This study aimed to determine the efficacy of subcutaneous injection of ivermectin for the treatment of Parascaris spp. infection in foals and to identify the Parascaris species present in the west and north of Iceland. A fecal egg count reduction (FECR) test (FECRT) was performed on 50 foals from 8 farms, including an untreated control group of 6 foals, from September to November 2019. The foals were between 3 and 5 mo of age at the start of the study and had not previously been treated with anthelmintic drugs. Each foal was treated subcutaneously with off-label use of Ivomec® injection 10 mg/ml or Noromectin® 1% at a dose of 0.2 mg/kg. The FECR for each farm was calculated in 2 ways, by the eggCounts package in R and by the Presidente formula (FECRT). Both calculation methods resulted in efficacy levels between 0% and 80.78%, indicating ivermectin resistance on all farms. We also confirmed, by karyotyping, that the species of equine ascarid present in the west and north of Iceland is Parascaris univalens. This study provides evidence for treatment failure of ivermectin against P. univalens infection in foals. Since Icelandic horses have been isolated on the island for more than 1,000 yr, this implies that resistance alleles have developed independently in the Icelandic Parascaris population. The actual clinical impact of ivermectin resistance is unknown but another drug of choice should be considered to treat Parascaris infection in foals in Iceland.

A survey of ivermectin resistance in Parascaris species infected foals in south-eastern Poland

Parascaris spp. are major gastro-intestinal nematodes that infect foals and can lead to respiratory symptoms, poor growth, and in some cases obstruction of the small intestine and death. Ivermectin resistance has been reported for Parascaris spp. in many countries. In Poland, the knowledge of the level of resistance against ivermectin in Parascaris spp. is limited. The aim of this study was to examine the efficacy of ivermectin against Parascaris spp. in foals from south-eastern Poland. Foals (n = 225 = reared in 7 stud farms) were treated orally with ivermectin paste. Faecal samples were collected from the rectum of each foal or from the environment straight after defaecation on 1 day prior and 2 weeks after deworming. A faecal egg count (FEC) was performed using the McMaster method with a minimum detection limit of 50 eggs/g. FEC reduction (FECR) was calculated using the Faecal Egg Count Reduction Test. The statistical analysis was limited to foals excreting more than 150 eggs/g before treatment and to stud farms with at least 6 foals excreting at or above this level. Confidence intervals were determined by 1000 bootstraps at farm level and the contribution of sex and age to FECR was quantified using a generalized equation estimation procedure. Parascaris spp. eggs were found in 40% of the foals. Following ivermectin treatment, Parascaris spp. eggs were identified in 28.4% of the foals. The mean estimated FECR ranged from 44% to 97% and average efficacy was 49.3%. FECR was more pronounced in older foals (P-values = 0. 003). The FECR was more pronounced in males than in females (P value = 0.028). This study is the first to indicate a reduced efficacy of ivermectin against Parascaris spp. in foals in Poland.

Molecular characterization of β-tubulin gene associated with benzimidazole resistance in larvae of field isolates of Parascaris (Nematoda: Ascarididae)

Since the past 2 decades, an increasing number of resistance to the Benzimidazoles (BZs) have been reported in nematode parasites of livestock. More recently, detection of single-nucleotide polymorphisms (SNPs) at codons of 167, 198 or 200 of the β-tubulin gene has been attributed to the occurrence of resistance. In the present study, we investigated the presence of those SNPs in the β-tubulin isotype-1 gene in different isolates of Parascaris in horse. Also, the mitochondrial (mt) and ribosomal genes were sequenced for species confirmation of the isolates. The analysis of sequences inferred from COII gene confirmed that those isolates were P. equorum. The distance between mt genes obtained here and several ascarid species in equids and other hosts suggests the need for the combination of more genetic data with morphologic and other diagnostic measures. The analysis on β-tubulin isotype-1 gene revealed no resistance-related SNPs or substitutions at the expected codon positions and selection pressure with BZs has not occurred for Parascaris worms. Although the molecular data showed the susceptibility of Parascaris isolates against BZs, other mechanisms of resistance should be also investigated to confirm the validity of molecular results.

An annotated catalogue of the ascaridoid nematode parasites of Chinese vertebrates

A catalogue, based on both examined specimens and the published literature, of all the ascaridoid nematodes recorded in China is presented. A total of 95 recognised species, representing 26 genera in five families, are reported. Detailed information on the type-host, type-locality, original reference, synonyms, annotated subsequent references of taxonomic importance, other host records, site of infection, location of type-specimens and distribution are listed for each recognised species. Additional comments on the taxonomic status of some species are also given. Moreover, some nomenclatural changes are proposed: (i) Toxascaris selenarctis Wang, 1965 and T. ailuri Wu, He & Hu, 1987 are placed in synonymy with Baylisascaris transfuga (Rudolphi, 1819); (ii) Raphidascaris lophii Wang & Wu, 1991 is a secondary homonym of R. lophii (Wu, 1949) and a replacement name, R. wangi nom. nov., is proposed for the former species; (iii) Aliascaris aetoplatea Luo, 2001 is transferred to Terranova Leiper & Atkinson, 1914, as T. aetoplatea (Luo, 2001) n. comb., and should be considered a species inquirenda; (iv) Ophidascaris orientalis (Wang, 1965) is resurrected as a valid species; (v) Phocascaris longispiculum Wang & Wu, 1991 and Ophidascaris agkistrodontis Wang, 1979 are treated as incertae sedis; and (vi) Hysterothylacium sauridae Li, Xu & Zhang, 2008 is listed as a nomen nudum.

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.

Deep amplicon sequencing of preselected isolates of Parascaris equorum in β-tubulin codons associated with benzimidazole resistance in other nematodes

Background

The development of anthelmintic resistance (AR) to macrocyclic lactones in the equine roundworm Parascaris equorum has resulted in benzimidazoles now being the most widely used substance to control Parascaris infections. However, over-reliance on one drug class is a risk factor for the development of AR. Consequently, benzimidazole resistance is widespread in several veterinary parasites, where it is associated with single nucleotide polymorphisms (SNPs) in drug targets encoded by the β-tubulin genes. The importance of these SNPs varies between different parasitic nematodes, but it has been hypothesised that they occur, at low allele frequencies, even in unselected populations. This study investigated whether these SNPs exist in the P. equorum population and tested the hypothesis that BZ resistance can develop from pre-existing SNPs in codons 167, 198 and 200 of the β-tubulin isotype 1 and 2 genes, reported to be associated with AR in strongylids. The efficacy of the oral paste formula fenbendazole on 11 farms in Sweden was also assessed.

Methods

Two isotype-specific primer pairs were designed, one on either side of the codon 167 and one on either side of codons 198 and 200. A pool of 100 000 larvae was sequenced using deep amplicon sequencing by Illumina HiSeq. Faecal egg count reduction test was used to assess the efficacy of fenbendazole.

Results

No SNPs were observed in codons 167, 198 or 200 of the β-tubulin isotype 1 or 2 genes of P. equorum, even though 100 000 larvae were sequenced. Faecal egg count reduction testing of fenbendazole showed that this anthelmintic was still 100% effective, meaning that the likelihood of finding high allele frequency of SNPs associated with benzimidazoles resistance in P. equorum was low. Unexpectedly, the allele frequencies observed in single worms were comparable to those in pooled samples.

Conclusions

We concluded that fenbendazole does not exert selection pressure on the β-tubulin genes of isotypes 1 and 2 in P. equorum. The fact that no pre-existing SNPs were found in codons 167, 198 and 200 in P. equorum also illustrates the difficulties in generalising about AR mechanisms between different taxonomic groups of nematodes.

Electronic supplementary material

The online version of this article (doi:10.1186/1756-3305-7-410) contains supplementary material, which is available to authorized users.