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

Occurrence of fenbendazole resistance in Parascaris spp. on breeding farms in Sweden

Anthelmintic resistance is an increasing problem in many gastrointestinal parasites of grazing animals. Among these, the equine roundworm, Parascaris spp., has developed wide-spread resistance to macrocyclic lactones over the past decades. Additionally, there are recent observations of emerging treatment failure of both tetrahydropyrimidine and fenbendazole. Therefore, the aims of this study were to further investigate the occurrence of fenbendazole resistance on breeding farms and to explore potential management-related risk factors associated with resistance in Parascaris spp. in Sweden. Eleven farms with 92 foals positive for Parascaris spp. were included in a faecal egg count reduction test during the years 2021-2023. According to the clinical protocol of the guidelines of the World Association for the Advancement of Veterinary Parasitology, fenbendazole resistance was present on four farms with efficacies varying from 45 % to 96 %. Having previously reported reduced efficacy on one of these farms, we can now confirm that fenbendazole resistance in Parascaris spp. has established. Farms with more than 40 yearly born foals had a significantly higher probability of having resistant Parascaris spp. Populations compared with smaller farms, (generalized linear model (GLM), t = 70.39, p < 0.001). In addition, there was a correlation between the number of foals on the farm and the frequency of yearly treatments showing that farms with < 20 foals were notably inclined to administer treatments twice during the first year (GLM, t=2.76, p < 0.05) in contrast to larger farms with > 40 foals that were using more frequent treatment intervals. In conclusion, this study confirms the establishment of fenbendazole resistance in Parascaris spp. populations on Swedish stud farms with the number of foals on the farm identified as a risk factor for development of anthelmintic resistance.

New Insights on Tools for Detecting β-Tubulin Polymorphisms in Trichuris trichiura Using rhAmpTM SNP Genotyping

Soil-transmitted helminth (STH) infections, commonly treated with benzimidazoles, are linked to resistance through single nucleotide polymorphisms (SNPs) at position 167, 198, or 200 in the β-tubulin isotype 1 gene. The aim of this study was to establish a novel genotyping assay characterized by its rapidity and specificity. This assay was designed to detect the presence of SNPs within the partial β-tubulin gene of Trichuris trichiura. This was achieved through the biallelic discrimination at codons 167, 198, and 200 by employing the competitive binding of two allele-specific forward primers. The specificity and reliability of this assay were subsequently confirmed using Trichuris samples isolated from captive primates. Furthermore, a molecular study was conducted to substantiate the utility of the β-tubulin gene as a molecular marker. The assays showed high sensitivity and specificity when applied to field samples. Nevertheless, none of the SNPs within the β-tubulin gene were detected in any of the adult worms or eggs from the analyzed populations. All specimens consistently displayed an SS genotype. The examination of the β-tubulin gene further validated the established close relationships between the T. trichiura clade and Trichuris suis clade. This reaffirms its utility as a marker for phylogenetic analysis.

Worldwide absence of canonical benzimidazole resistance-associated mutations within β-tubulin genes from Ascaris

BACKGROUND

The giant roundworm Ascaris is an intestinal nematode, causing ascariasis by infecting humans and pigs worldwide. Recent estimates suggest that Ascaris infects over half a billion people, with chronic infections leading to reduced growth and cognitive ability. Ascariasis affects innumerable pigs worldwide and is known to reduce production yields via decreased growth and condemnation of livers. The predominant anthelminthic drugs used to treat ascariasis are the benzimidazoles. Benzimidazoles interact with β-tubulins and block their function, and several benzimidazole resistance-associated mutations have been described in the β-tubulins of ruminant nematodes. Recent research on ascarids has shown that these canonical benzimidazole resistance-associated mutations are likely not present in the β-tubulins of Ascaris, Ascaridia or Parascaris, even in phenotypically resistant populations.

METHODS

To further determine the putative absence of key β-tubulin polymorphisms, we screened two β-tubulin isotypes of Ascaris, highly expressed in adult worms. Using adult and egg samples of Ascaris obtained from pigs and humans worldwide, we performed deep amplicon sequencing to look for canonical resistance-associated mutations in Ascaris β-tubulins. Subsequently, we examined these data in closer detail to study the population dynamics of Ascaris and genetic diversity within the two isotypes and tested whether genotypes appeared to partition across human and pig hosts.

RESULTS

In the 187 isolates, 69 genotypes were found, made up of eight haplotypes of β-tubulin isotype A and 20 haplotypes of isotype B. Single nucleotide polymorphisms were seen at 14 and 37 positions for β-tubulin isotype A and isotype B, respectively. No evidence of any canonical benzimidazole resistance-associated mutations was found in either human- or pig-derived Ascaris isolates. There was, however, a difference in the genetic diversity of each isotype and distribution of β-tubulin genotypes between human- and pig-derived Ascaris. Statistical tests of population differentiation show significant differences (p < 0.001) between pig- and human-derived worms; however, more diversity was seen between worms from different populations than worms from different hosts.

CONCLUSIONS

Our work suggests an absence of canonical β-tubulin mutations within Ascaris, but alternative modes of anthelminthic resistance may emerge necessitating continued genetic scrutiny alongside monitoring of drug efficacy.

Understanding anthelmintic resistance in livestock using "omics" approaches

Widespread and improper use of various anthelmintics, genetic, and epidemiological factors has resulted in anthelmintic-resistant (AR) helminth populations in livestock. This is currently quite common globally in different livestock animals including sheep, goats, and cattle to gastrointestinal nematode (GIN) infections. Therefore, the mechanisms underlying AR in parasitic worm species have been the subject of ample research to tackle this challenge. Current and emerging technologies in the disciplines of genomics, transcriptomics, metabolomics, and proteomics in livestock species have advanced the understanding of the intricate molecular AR mechanisms in many major parasites. The technologies have improved the identification of possible biomarkers of resistant parasites, the ability to find actual causative genes, regulatory networks, and pathways of parasites governing the AR development including the dynamics of helminth infection and host-parasite infections. In this review, various "omics"-driven technologies including genome scan, candidate gene, quantitative trait loci, transcriptomic, proteomic, and metabolomic approaches have been described to understand AR of parasites of veterinary importance. Also, challenges and future prospects of these "omics" approaches are also discussed.

Population genomics of helminth parasites

Next generation sequencing technologies have facilitated a shift from a few targeted loci in population genetic studies to whole genome approaches. Here, we review the types of questions and inferences regarding the population biology and evolution of parasitic helminths being addressed within the field of population genomics. Topics include parabiome, hybridization, population structure, loci under selection and linkage mapping. We highlight various advances, and note the current trends in the field, particularly a focus on human-related parasites despite the inherent biodiversity of helminth species. We conclude by advocating for a broader application of population genomics to reflect the taxonomic and life history breadth displayed by helminth parasites. As such, our basic knowledge about helminth population biology and evolution would be enhanced while the diversity of helminths in itself would facilitate population genomic comparative studies to address broader ecological and evolutionary concepts.

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.

Exploring the β-tubulin gene family in a benzimidazole-resistant Parascaris univalens population

Benzimidazole (BZ) drugs are frequently used to treat infections with the equine ascarid Parascaris univalens due to increasing resistance to macrocyclic lactones and pyrantel. Benzimidazole resistance is rare in ascarids in contrast to strongyle parasites where this resistance is widespread. In strongyles, single nucleotide polymorphisms (SNPs) at codons 167, 198 and 200 in a β-tubulin gene have been correlated to BZ resistance, but little is known about the β-tubulin genes and their possible involvement in BZ resistance in P. univalens and other ascarids. Previously two β-tubulin genes have been identified in P. univalens. In this study, we present five additional β-tubulin genes as well as the phylogenetic relationship of all seven genes to β-tubulins of other clade III and V nematodes. In addition, the efficacy of fenbendazole for treatment of P. univalens on a Swedish stud farm was studied in 2019 and 2020 using faecal egg count reduction test. Reductions varied from 73% to 88%, indicating the presence of a resistant P. univalens population on the farm. The emergence of BZ resistance emphasizes the need for development of molecular markers for rapid and more sensitive detection of resistant populations. We therefore investigated whether possible SNPs at positions 167, 198 or 200 in any of the β-tubulin genes could be used to distinguish between resistant and susceptible P. univalens populations. Amplicon sequencing covering the mutation sites 167, 198 and 200 in all seven β-tubulin genes revealed an absence of SNPs in both resistant and susceptible populations, suggesting that the mechanism behind BZ resistance in ascarids is different from that in strongyle nematodes and the search for a molecular marker for BZ resistance in P. univalens needs to continue.

•

First case of fenbendazole resistance in Parascaris univalens in Europe.

•

The P. univalens β-tubulin family contains seven genes.

•

P. univalens β-tubulin genes cluster with β-tubulins from other clade V nematodes.

•

No resistance associated SNPs were identified in P. univalens β-tubulin genes.

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.

Transcriptional responses in Parascaris univalens after in vitro exposure to ivermectin, pyrantel citrate and thiabendazole

Background

Parascaris univalens is a pathogenic parasite of foals and yearlings worldwide. In recent years, Parascaris spp. worms have developed resistance to several of the commonly used anthelmintics, though currently the mechanisms behind this development are unknown. The aim of this study was to investigate the transcriptional responses in adult P. univalens worms after in vitro exposure to different concentrations of three anthelmintic drugs, focusing on drug targets and drug metabolising pathways.

Methods

Adult worms were collected from the intestines of two foals at slaughter. The foals were naturally infected and had never been treated with anthelmintics. Worms were incubated in cell culture media containing different concentrations of either ivermectin (10⁻⁹ M, 10⁻¹¹ M, 10⁻¹³ M), pyrantel citrate (10⁻⁶ M, 10⁻⁸ M, 10⁻¹⁰ M), thiabendazole (10⁻⁵ M, 10⁻⁷ M, 10⁻⁹ M) or without anthelmintics (control) at 37 °C for 24 h. After incubation, the viability of the worms was assessed and RNA extracted from the anterior region of 36 worms and sequenced on an Illumina NovaSeq 6000 system.

Results

All worms were alive at the end of the incubation but showed varying degrees of viability depending on the drug and concentration used. Differential expression (Padj < 0.05 and log2 fold change ≥ 1 or ≤ − 1) analysis showed similarities and differences in the transcriptional response after exposure to the different drug classes. Candidate genes upregulated or downregulated in drug exposed worms include members of the phase I metabolic pathway short-chain dehydrogenase/reductase superfamily (SDR), flavin containing monooxygenase superfamily (FMO) and cytochrome P450-family (CYP), as well as members of the membrane transporters major facilitator superfamily (MFS) and solute carrier superfamily (SLC). Generally, different targets of the anthelmintics used were found to be upregulated and downregulated in an unspecific pattern after drug exposure, apart from the GABA receptor subunit lgc-37, which was upregulated only in worms exposed to 10⁻⁹ M of ivermectin.

Conclusions

To our knowledge, this is the first time the expression of lgc-37 and members of the FMO, SDR, MFS and SLC superfamilies have been described in P. univalens and future work should be focused on characterising these candidate genes to further explore their potential involvement in drug metabolism and anthelmintic resistance.

Optimizing and Evaluating the Antihelminthic Activity of the Biocompatible Zinc Oxide Nanoparticles Against the Ascaridid Nematode, Parascaris equorum In Vitro

PURPOSE

In the present study, the effect of different biocompatible concentrations from ZnO nanoparticles (ZnO NPs) on the physiological state and surface topography of the nematode P. equorum was determined in vitro.

METHODS

Different concentrations of ZnO NPs (100, 200, 300 and 400 mg/l) synthesized using the egg white were prepared followed by the incubation of parasitic worms with these concentrations in vitro. The physiological state of treated worms such as oxidative stress markers, enzymatic activities and biochemical parameters in addition to the surface topography was determined and compared with control untreated worms.

RESULTS

In comparison to control worms, it was observed that at high concentrations of ZnO NPs, most of the treated worms showed an increase in the levels of ALT, AST and ALP (worm muscle damage, and gonad injury); enhancement of the total protein content (worm cellular dysfunction); significant increase in MDA level (free radical-mediated worm cell membrane damage); depletion in GST and GSH activities (reduced ability to clear toxic compounds like lipid peroxides); CAT depletion (superoxide dismutase and hydrogen peroxide toxicity) and NO increase (detoxification activity and stressful conditions on worms). SEM showed that there was a modified morphological appearance in the surface of treated worms; lips were wrinkled with irregularly arranged denticles, weathering of cuticle, bursts of cuticle layers, disruption of surface annulations and erosion of surface papillae of male around the cloacal opening.

CONCLUSION

ZnO NPs at environmentally relevant concentrations achieved a significant antihelminthic activity against P. equorum which represents a successful model used in parasite control experiments.

First identification of the benzimidazole resistance-associated F200Y SNP in the beta-tubulin gene in Ascaris lumbricoides

The main control strategy for Ascaris lumbricoides is mass drug administration (especially with benzimidazoles), which can select strains of parasites resistant to treatment. Mutations in the beta-tubulin isotype-1 gene at codons 167, 198 and 200 have been linked to benzimidazole resistance in several nematodes. The mutation in codon 200 is the most frequent in different species of parasites, as previously observed in Necator americanus and Trichuris trichiura; however, this mutation has never been found in populations of A. lumbricoides. This study aimed to screen for single nucleotide polymorphisms (SNPs) in the beta-tubulin isotype-1 gene at codon 200 in A. lumbricoides. We developed a technique based on an amplification refractory mutation system (ARMS-PCR) for the analysis of 854 single A. lumbricoides eggs collected from 68 human stool samples from seven Brazilian states. We detected the mutation in codon 200 at a frequency of 0.5% (4/854). This is the first report of this mutation in A. lumbricoides. Although the observed frequency is low, its presence indicates that these parasite populations have the potential to develop high levels of resistance in the future. The methodology proposed here provides a powerful tool to screen for the emergence of anthelmintic resistance mutations in parasitic nematode populations.

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.