In vitro selection of Haemonchus contortus for benzimidazole resistance reveals a mutation at amino acid 198 of β-tubulin

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.

Transcriptomic analysis of albendazole resistance in human diarrheal parasite Giardia duodenalis

Benzimidazole-2-carbamates (BZ, e.g., albendazole; ALB), which bind β-tubulin to disrupt microtubule polymerization, are one of two primary compound classes used to treat giardiasis. In most parasitic nematodes and fungi, BZ-resistance is caused by β-tubulin mutations and its molecular mode of action (MOA) is well studied. In contrast, in Giardia duodenalis BZ MOA or resistance is less well understood, may involve target-specific and broader impacts including cellular damage and oxidative stress, and its underlying cause is not clearly determined. Previously, we identified acquisition of a single nucleotide polymorphism, E198K, in β-tubulin in ALB-resistant (ALB-R) G. duodenalis WB-1B relative to ALB-sensitive (ALB-S) parental controls. E198K is linked to BZ-resistance in fungi and its allelic frequency correlated with the magnitude of BZ-resistance in G. duodenalis WB-1B. Here, we undertook detailed transcriptomic comparisons of these ALB-S and ALB-R G. duodenalis WB-1B cultures. The primary transcriptional changes with ALB-R in G. duodenalis WB-1B indicated increased protein degradation and turnover, and up-regulation of tubulin, and related genes, associated with the adhesive disc and basal bodies. These findings are consistent with previous observations noting focused disintegration of the disc and associated structures in Giardia duodenalis upon ALB exposure. We also saw transcriptional changes with ALB-R in G. duodenalis WB-1B consistent with prior observations of a shift from glycolysis to arginine metabolism for ATP production and possible changes to aspects of the vesicular trafficking system that require further investigation. Finally, we saw mixed transcriptional changes associated with DNA repair and oxidative stress responses in the G. duodenalis WB-1B line. These changes may be indicative of a role for H₂O₂ degradation in ALB-R, as has been observed in other G. duodenalis cell cultures. However, they were below the transcriptional fold-change threshold (log₂FC > 1) typically employed in transcriptomic analyses and appear to be contradicted in ALB-R G. duodenalis WB-1B by down-regulation of the NAD scavenging and conversion pathways required to support these stress pathways and up-regulation of many highly oxidation sensitive iron-sulphur (FeS) cluster based metabolic enzymes.

In Silico Docking of Nematode β-Tubulins With Benzimidazoles Points to Gene Expression and Orthologue Variation as Factors in Anthelmintic Resistance

The efficacy of benzimidazole anthelmintics can vary depending on the target parasite, with Ascaris nematodes being highly responsive, and whipworms being less responsive. Anthelmintic resistance has become widespread, particularly in strongyle nematodes such as Haemonchus contortus in ruminants, and resistance has recently been detected in hookworms of humans and dogs. Past work has shown that there are multiple β-tubulin isotypes in helminths, yet only a few of these contribute to benzimidazole interactions and resistance. The β-tubulin isotypes of ascarids and soil-transmitted helminths were identified by mining available genome data, and phylogenetic analysis showed that the ascarids share a similar repertoire of seven β-tubulin isotypes. Strongyles also have a consistent pattern of four β-tubulin isotypes. In contrast, the whipworms only have two isotypes, with one of these clustering more basally and distinct from any other group. Key β-tubulin isotypes selected based on previous studies were the focus of in silico molecular docking simulations to look at the interactions with benzimidazoles. These showed that all β-tubulins had similar interactions with benzimidazoles and maintained the key bond with residue E198 in all species, indicating similar mechanisms of action. However, the interaction was stronger and more consistent in the strongyles and whipworms than it was in the ascarids. Alteration of β-tubulin isotypes with the common resistance-associated mutations originally identified in H. contortus resulted in similar interaction modeling for all species. In conclusion, ascarids, strongyles, and whipworms all have their own unique repertoire of β-tubulins, which could explain why benzimidazole resistance and susceptibility varies between these groups of parasites. These data complement recent work that has highlighted the roles of essential residues in benzimidazole drug binding and shows that there is a separation between strongyle parasites that frequently develop resistance and ascarid parasites, which have been much less prone to developing resistance.

Characterization of the β-tubulin gene family in Ascaris lumbricoides and Ascaris suum and its implication for the molecular detection of benzimidazole resistance

Background

The treatment coverage of control programs providing benzimidazole (BZ) drugs to eliminate the morbidity caused by soil-transmitted helminths (STHs) is unprecedently high. This high drug pressure may result in the development of BZ resistance in STHs and so there is an urgent need for surveillance systems detecting molecular markers associated with BZ resistance. A critical prerequisite to develop such systems is an understanding of the gene family encoding β-tubulin proteins, the principal targets of BZ drugs.

Methodology and principal findings

First, the β-tubulin gene families of Ascaris lumbricoides and Ascaris suum were characterized through the analysis of published genomes. Second, RNA-seq and RT-PCR analyses on cDNA were applied to determine the transcription profiles of the different gene family members. The results revealed that Ascaris species have at least seven different β-tubulin genes of which two are highly expressed during the entire lifecycle. Third, deep amplicon sequencing was performed on these two genes in more than 200 adult A. lumbricoides (Ethiopia and Tanzania) and A. suum (Belgium) worms, to investigate the intra- and inter-species genetic diversity and the presence of single nucleotide polymorphisms (SNPs) that are associated with BZ resistance in other helminth species; F167Y (TTC>TAC or TTT>TAT), E198A (GAA>GCA or GAG>GCG), E198L (GAA>TTA) and F200Y (TTC>TAC or TTT>TAT). These particular SNPs were absent in the two investigated genes in all three Ascaris populations.

Significance

This study demonstrated the presence of at least seven β-tubulin genes in Ascaris worms. A new nomenclature was proposed and prioritization of genes for future BZ resistance research was discussed. This is the first comprehensive description of the β-tubulin gene family in Ascaris and provides a framework to investigate the prevalence and potential role of β-tubulin sequence polymorphisms in BZ resistance in a more systematic manner than previously possible.

Benzimidazole (BZ) drugs remain the standard of treatment in large-scale deworming programs that aim to control the morbidity caused by intestinal worms. As these deworming programs are expanding world-wide, there is an increasing risk of worms becoming resistant to BZ drugs, highlighting the necessity for tools to detect gene mutations associated with drug resistance. However, the development of such tools is impeded by a lack of insights into the genes that are coding for β-tubulin proteins, which are the principal targets of BZ drugs. The aim of this study was to comprehensively characterize these genes in the worm species Ascaris lumbricoides and Ascaris suum. The findings highlight that these species have at least seven β-tubulin genes. Only two genes are highly expressed throughout the different life stages of the worm, and hence are more likely to be involved in the development of BZ resistance. No mutations that have previously been associated with BZ resistance in other intestinal worms were found. This study provides a baseline towards more efficient and accurate monitoring of drug resistance in large-scale deworming programs.

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.

In vitro selection of Giardia duodenalis for Albendazole resistance identifies a β-tubulin mutation at amino acid E198K

Benzimidazole-2-carbamate (BZ) compounds, including Albendazole (Alb), are one of just two drug classes approved to treat the gastrointestinal protist Giardia duodenalis. Benzimidazoles bind to the tubulin dimer interface overlapping the colchicine binding site (CBS) of β-tubulin, thereby inhibiting microtubule polymerisation and disrupting microtubule networks. These BZ compounds are widely used as anthelmintic, anti-fungal and anti-giardial drugs. However, in helminths and fungi, BZ-resistance is widespread and caused by specific point mutations primarily occurring at F167, E198 and F200 in β-tubulin isoform 1. BZ-resistance in Giardia is reported clinically and readily generated in vitro, with significant implications for Giardia control. In Giardia, BZ mode of action (MOA) and resistance mechanisms are presumed but not proven, and no mutations in β-tubulin have been reported in association with Alb resistance (AlbR). Herein, we undertook detailed in vitro drug-susceptibility screens of 13 BZ compounds and 7 Alb structural analogues in isogenic G. duodenalis isolates selected for AlbR and podophyllotoxin, another β-tubulin inhibitor, as well as explored cross-resistance to structurally unrelated, metronidazole (Mtz). AlbR lines exhibited co-resistance to many structural variants in the BZ-pharmacophore, and cross-resistance to podophyllotoxin. AlbR lines were not cross-resistant to Mtz, but MtzR lines had enhanced survival in Alb. Lastly, Alb analogues with longer thioether substituents had decreased potency against our AlbR lines. In silico modelling indicated the Alb-β-tubulin interaction in Giardia partially overlaps the CBS and corresponds to residues associated with BZ-resistance in helminths and fungi (F167, E198, F200). Sequencing of Giardia β-tubulin identified a single nucleotide polymorphism resulting in a mutation from glutamic acid to lysine at amino acid 198 (E198K). To our knowledge, this is the first β-tubulin mutation reported for protistan BZ-resistance. This study provides insight into BZ mode of action and resistance in Giardia, and presents a potential avenue for a genetic test for clinically resistance isolates.

•

In vitro albendazole-resistant Giardia were broadly resistant to benzimidazole-2-carbamates.

•

Cross-resistance to structurally unrelated microtubule inhibitors was observed.

•

The predicted Giardia benzimidazole binding overlaps the colchicine binding site.

•

The E198K β-tubulin mutation was identified in the albendazole-resistant line.

•

Metronidazole-resistance may enhance acquisition of albendazole-resistance.

Haemonchosis: A Challenging Parasitic Infection of Sheep and Goats

The paper reviews the challenges about haemonchosis-a significant and common parasitic infection of small ruminants. Haemonchus contortus is a highly pathogenic parasite that localises in the abomasum of affected animals and exerts its pathogenicity by blood-sucking activity, adversely affecting the health and productivity of animals. The first challenge is the uneven distribution of the infection globally, this being more prevalent in tropical and subtropical and warm temperate and summer rainfall regions than in cool and cold temperate and arid regions; hence, this leads in differences in the approaches required for its control. Another challenge is the widespread presence of Haemonchus strains resistant to the various anthelmintics available: Benzimidazoles, imidazothiazoles, macrocyclic lactones, closantel and monepantel, which makes the control of the infection difficult. The third challenge refers to the difficulty of diagnosing the disease, given that field evidence can provide suspicion about the infection, which needs to be subsequently confirmed by laboratory tests through parasitological or molecular techniques. The final challenge relates to the difficulties in the control of the infection and the necessity to use pharmaceutical products cautiously and with a planned approach, to avoid further development of anthelmintic resistance, also given that use of a recently licenced vaccine is not widespread. In conclusion, at the moment, we should be concerned, but not worried, about this infection, and apply correctly the appropriate health management plans.

Contrasting patterns of isotype-1 β-tubulin allelic diversity in Haemonchus contortus and Haemonchus placei in the southern USA are consistent with a model of localised emergence of benzimidazole resistance

The benzimidazoles are one of the most important broad-spectrum anthelmintic drug classes for parasitic nematode control in domestic animals and humans. They have been widely used in livestock, particularly in small ruminants for over 40 years. This has resulted in widespread resistance in small ruminant gastrointestinal nematode parasite species, especially Haemonchus contortus. Benzimidazole resistance mutations have also been reported in Haemonchus placei, but only at low frequencies, suggesting resistance is at a much earlier stage of emergence than is the case for H. contortus. Here, we investigate the haplotype diversity of isotype-1 β-tubulin benzimidazole resistance mutations and the population genetic structure of H. contortus and H. placei populations from sheep and cattle from the southern USA. Microsatellite genotyping revealed a low level of genetic differentiation in six H.placei and seven H. contortus populations examined. This is consistent with several previous studies from other regions, mainly in H. contortus, supporting a model of high gene flow between parasite populations. There was a single F200Y(TAC) haplotype present in all six H. placei populations across Georgia, Florida and Arkansas. In contrast, there were at least two different F200Y(TAC) haplotypes (up to four) and two different F167Y(TAC) haplotypes across the seven H. contortus populations studied. These results provide further evidence to support a model for benzimidazole resistance in Haemonchus spp, in which resistance mutations arise from a single, or the small number of locations, in a region during the early phases of emergence, and subsequently spread due to animal movement.

Nicotinic acetylcholine receptors: Ex-vivo expression of functional, non-hybrid, heteropentameric receptors from a marine arthropod, Lepeophtheirus salmonis

Nicotinic acetylcholine receptors (nAChRs) are ligand-gated ion channels mostly located in the post-synaptic membrane of cholinergic synapses. The natural neurotransmitter is acetylcholine, but they are also the direct targets for neonicotinoids, chemicals widely used against ectoparasites, arthropod vectors and agricultural pests. There are significant concerns regarding adverse effects of neonicotinoids on beneficial insects. In arthropods, functional nAChRs made of α subunits have been expressed from Drosophila genes, and hybrid receptors (sometimes also referred to as chimeric receptors) using species-specific α subunits and vertebrate β subunits have been expressed ex-vivo. Arthropod-specific nAChRs made of both α and β subunits from the target species have not been expressed ex-vivo. The aim of the current study was to express such receptors in Xenopus oocytes using only genes from Lepeophtheirus salmonis, to characterize them and study their modulation. Genes encoding α and β subunits of the nAChRs and three ancillary proteins, RIC-3, UNC-50 and UNC-74 were identified in the L. salmonis genome, subjected to RACE-PCR, cloned into an expression vector and the cRNA produced was then injected into Xenopus laevis oocytes. Co-expression of the ancillary proteins was essential for the successful expression of the L. salmonis nAChRs with both α and β subunits. Two functional nAChRs were identified: Lsa-nAChR1 consisting of α1, α2, β1 and β2 subunits, reconstituted to one distinct receptor, while Lsa-nAChR2, consisting of α3, β1 and β2 subunits reconstitutes receptors with two distinct characteristics. Out of seven neonicotinoids tested, six worked as partial agonist of Lsa-nAChR1 while only three did so for Lsa-nAChR2. Four non-neonicotinoid compounds tested had no effect on either of the nAChRs. The study demonstrated that fully functional, non-hybrid nAChRs containing both α and β subunits from an arthropod can be reconstituted ex-vivo by co-expression of essential ancillary proteins. Such models would be valuable for in-depth studies of effects by neonicotinoids and other compounds on target pests, as well as for studies of adverse effects on non-target arthropods.

Nicotinic acetylcholine receptors, nAChRs, respond to the neurotransmitter acetylcholine or drugs like nicotine. These receptors are targets for neonicotinoids, the most commonly used compounds against ectoparasites and agricultural pests. In-depth studies of the function of these channels in arthropods are sparse, as no groups managed to reconstitute functional nAChRs made of both α and β subunits using genes only from the target arthropod in an ex-vivo system. We report the successful assembly of non-hybrid, fully functional nAChRs containing both α and β subunits from a marine arthropod, assembled and expressed in Xenopus laevis oocytes. We identified two possible combinations of α and β subunits producing functional receptors. We found ancillary proteins to be essential for successful expression and assembly of both α and β subunits into a functional receptor. The findings of the present study provide a basis for studying native nAChRs from arthropods, with a switch from hybrid nAChRs to species-specific native nAChRs.

Frequency of Resistance to Benzimidazoles of Haemonchus contortus Helminths from Dairy Sheep, Goats, Cattle and Buffaloes in Greece

The study investigated the presence of resistance to benzimidazoles in Haemonchus contortus helminths from ruminant species in Greece through the detection of the Phe/Tyr polymorphism in the amino acid at position 200 of the β-tubulin protein. In total, 288 adult female H. contortus helminths collected from the abomasum of various ruminant animals in Greece were tested. Of these, 96 were collected from sheep, 96 from goats, 48 from cattle, and 48 from buffaloes. The frequencies of the homozygous and heterozygous resistant genotypes at the position 200 of the β-tubulin gene of helminths recovered from sheep were 96.9% and 3.1%, respectively. The frequencies of the homozygous and heterozygous resistant genotypes, respectively, were 100.0% and 0.0% in helminths from goats, 25.0% and 75.0% in helminths from cattle and 8.3% and 91.7% in helminths from buffaloes. In all parasitic populations, no homozygous susceptible genotypes were detected. The present study highlighted, for the first time, the emergence of benzimidazole-resistant H. contortus in goats, cattle, and buffaloes in Greece, using an allele-specific PCR. It is postulated that benzimidazole-resistant alleles were transferred from sheep or goats to cattle and buffaloes at the commonly grazing pastures in Greece.

A perspective on the discovery of selected compounds with anthelmintic activity against the barber's pole worm-Where to from here?

Parasitic roundworms (nematodes) cause substantial morbidity and mortality in animals worldwide. Anthelmintic treatment is central to controlling these worms, but widespread resistance to most of the commercially available anthelmintics for veterinary and agricultural use is compromising control, such that there is an urgency to discover new and effective drugs. The purpose of this article is to review information on parasitic nematodes, the treatment and control of parasitic nematode infections and aspects of discovering new anthelmintics in the context of anthelmintic resistance problems, and then to discuss some progress that our group has made in identifying selected compounds with activity against nematodes. The focus of our recent work has been on discovering new chemical entities and known drugs with anthelmintic activities against Haemonchus contortus as well as other socioeconomically important parasitic nematodes for subsequent development. Using whole worm-based phenotypic assays, we have been screening compound collections obtained via product-development-partnerships and/or collaborators, and active compounds have been assessed for their potential as anthelmintic candidates. Following the screening of 15,333 chemicals from five distinct compound collections against H. contortus, we have discovered one new chemical entity (designated SN00797439), two human kinase inhibitors (SNS-032 and AG-1295), 14 tetrahydroquinoxaline analogues, one insecticide (tolfenpyrad) and two tolfenpyrad (pyrazole-5-carboxamide) derivatives (a-15 and a-17) with anthelmintic activity in vitro. Some of these 20 'hit' compounds have selectivity against H. contortus in vitro when compared to particular human cell lines. In our opinion, some of these compounds could represent starting points for 'lead' development. Accordingly, the next research steps to be pursued include: (i) chemical optimisation of representative chemicals via structure-activity relationship (SAR) evaluations; (ii) assessment of the breadth of spectrum of anthelmintic activity on a range of other parasitic nematodes, such as strongyloids, ascaridoids, enoplids and filarioids; (iii) detailed investigations of the absorption, distribution, metabolism, excretion and toxicity (ADMET) of optimised chemicals with broad nematocidal or nematostatic activity; and (iv) establishment of the modes of action of lead candidates.

Novel and selective inactivators of Triosephosphate isomerase with anti-trematode activity

Trematode infections such as schistosomiasis and fascioliasis cause significant morbidity in an estimated 250 million people worldwide and the associated agricultural losses are estimated at more than US$ 6 billion per year. Current chemotherapy is limited. Triosephosphate isomerase (TIM), an enzyme of the glycolytic pathway, has emerged as a useful drug target in many parasites, including Fasciola hepatica TIM (FhTIM). We identified 21 novel compounds that selectively inhibit this enzyme. Using microscale thermophoresis we explored the interaction between target and compounds and identified a potent interaction between the sulfonyl-1,2,4-thiadiazole (compound 187) and FhTIM, which showed an IC₅₀ of 5 µM and a K_d of 66 nM. In only 4 hours, this compound killed the juvenile form of F. hepatica with an IC₅₀ of 3 µM, better than the reference drug triclabendazole (TCZ). Interestingly, we discovered in vitro inhibition of FhTIM by TCZ, with an IC₅₀ of 7 µM suggesting a previously uncharacterized role of FhTIM in the mechanism of action of this drug. Compound 187 was also active against various developmental stages of Schistosoma mansoni. The low toxicity in vitro in different cell types and lack of acute toxicity in mice was demonstrated for this compound, as was demonstrated the efficacy of 187in vivo in F. hepatica infected mice. Finally, we obtained the first crystal structure of FhTIM at 1.9 Å resolution which allows us using docking to suggest a mechanism of interaction between compound 187 and TIM. In conclusion, we describe a promising drug candidate to control neglected trematode infections in human and animal health.

Development and evaluation of a loop-mediated isothermal amplification (LAMP) assay for the detection of the E198A SNP in the isotype-1 β-tubulin gene of Haemonchus contortus populations in China

Haemonchus contortus is one of the most important gastrointestinal nematodes of small ruminants around the world, seriously hampering the healthy development of the sheep industry. The control of this parasite mainly depends on anthelmintics, however, drug resistance of H. contortus has become a serious problems worldwide. Previous studies demonstrated that the E198A (GAA to GCA), a single nucleotide polymorphism (SNP) in the isotype-1 β-tubulin gene is associated with benzimidazole resistance in H. contortus. However, only PCR-RFLP and ARMS-PCR methods have been previously used for the detection of the E198A mutation. In the present study, a loop-mediated isothermal amplification (LAMP) assay was established for rapid detection of the E198A SNP in H. contortus. The results showed that optimization of LAMP reaction reagents and conditions could achieve this. The resulting amplicons were visualized by adding hydroxynaphthol blue dye (HNB) prior to amplification. The color of LAMP products amplified without DNA or from DNA from worms with the E198A homozygous susceptible genotype was still violet, but the products with DNA from worms with the E198A heterozygous genotype or the E198A resistant homozygous genotype changed to sky blue. The specificity of this method was further verified by sequencing, which confirmed the successful LAMP detection of the E198A mutation with high specificity. In conclusion, the developed LAMP method has high specificity and good reproducibility for screening the E198A SNP of isotype-1 β-tubulin gene of H. contortus of field samples without using sophisticated equipment, providing useful technique for the rapid detection and thus prevention and control of benzimidazole resistant H. contortus infections.

The uptake, effects and biotransformation of monepantel in meadow plants used as a livestock feed

Drugs are potentially dangerous environmental contaminants, as they are designed to have biological effects at low concentrations. Monepantel (MOP), an amino-acetonitrile derivative, is frequently used veterinary anthelmintics, but information about MOP environmental circulation and impact is almost non-existent. We studied the phytotoxicity, uptake and biotransformation of MOP in two fodder plants, Plantago lanceolata and Medicago sativa. The seeds and whole plant regenerants were cultivated with MOP. The plant roots and the leaves were collected after 1, 2, 3, 4, 5 and 6 weeks of cultivation. The lengths of roots and proline concentrations in the roots and leaves were measured to evaluate MOP phytotoxicity. The UHPLC-MS/MS technique with a Q-TOF mass analyser was used for the identification and semi-quantification of MOP and its metabolites. Our results showed no phytotoxicity of MOP. However, both plants were able to uptake, transport and metabolize MOP. Comparing both plants, the uptake of MOP was much more extensive in Medicago sativa (almost 10-times) than in Plantago lanceolate. Moreover, 9 various metabolites of MOP were detected in Medicago sativa, while only 7 MOP metabolites were found in Plantago lanceolata. Based on metabolites structures, scheme of the metabolic pathways of MOP in both plants was proposed. MOP and its main metabolite (MOP sulfone), both anthelmintically active, were present not only in roots but also in leaves that can be consumed by animals. This indicates the potential for undesirable circulation of MOP in the environment, which could lead to many pharmacological and toxicological consequences.

Development of amplicon sequencing for the analysis of benzimidazole resistance allele frequencies in field populations of gastrointestinal nematodes

Anthelmintic resistant gastrointestinal helminths have become a major cause of poor health in sheep and goats. Sensitive and specific molecular markers are needed to monitor the genotypic frequency of resistance in field parasite populations. Gastrointestinal nematode resistance to benzimidazole is caused by a mutation in one of three positions within the isotype 1 β-tubulin gene. In the absence of markers for resistance to the other broad spectrum anthelmintic classes, these provide a relevant study example. Determination of the prevalence of these single nucleotide polymorphisms in field nematode populations can be impractical using conventional molecular methods to examine individual parasites; which can be laborious and lack sensitivity in determining low levels of resistance in parasite populations. Here, we report the development of a novel method based on an Illumina MiSeq deep amplicon sequencing platform to sequence the isotype 1 β-tubulin locus of the small ruminant gastrointestinal nematode, Teladorsagia circumcincta, and determine the frequency of the benzimidazole resistance mutations. We validated the method by assessing sequence representation bias, comparing the results of Illumina MiSeq and pyrosequencing, and applying the method to populations containing known proportions of resistant and susceptible larvae. We applied the method to field samples collected from ewes and lambs on over a period of one year on three farms, each highlighting different aspects of sheep management and approaches to parasite control. The results show opportunities to build hypotheses with reference to selection pressures leading to differences in resistance allele frequencies between sampling dates, farms and ewes or lambs, and to consider the impact of their genetic fixation or otherwise. This study provides proof of concept of a practical, accurate, sensitive and scalable method to determine frequency of anthelmintic resistance mutations in gastrointestinal nematodes in field studies and as a management tool for livestock farmers.

In vivo selection for Haemonchus contortus resistance to monepantel

Gastrointestinal nematodes significantly affect the ovine industry, and Haemonchus contortus is considered the most pathogenic parasite in tropical regions. This situation is aggravated when the main strategy to control worms fails because of the genetic resistance that parasites acquire against anthelmintics. Aiming to anticipate the events involved in anthelmintic resistance, we induced monepantel resistance in H. contortus by in vivo subdosing of sheep hosts. Four successive passages of a monepantel-susceptible H. contortus isolate in Santa Ines or Ile de France sheep hosts resulted in three monepantel-resistant (efficacy varying from 0 to 58.5%) H. contortus isolates. Sheep hosts were treated from 0.075 mg/kg to the therapeutic dose of 2.5 mg/kg of monepantel in 19-26 rounds of selection for 112-133 weeks. Success in inducing H. contortus resistance to monepantel may have been affected by worm burden and by host-parasite interactions, including a possible effect of the breed of sheep hosts. We conclude that subdosing of sheep, although time-consuming, is an efficient in vivo strategy for the induction of monepantel resistance in H. contortus. The resistant parasites can be used in further studies to elucidate the genetic and biochemical events involved in the acquisition of anthelmintic resistance.

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.

Isolation and characterization of a naturally occurring multidrug-resistant strain of the canine hookworm, Ancylostoma caninum

Soil-transmitted nematodes infect over a billion people and place several billion more at risk of infection. Hookworm disease is the most significant of these soil-transmitted nematodes, with over 500 million people infected. Hookworm infection can result in debilitating and sometimes fatal iron-deficiency anemia, which is particularly devastating in children and pregnant women. Currently, hookworms and other soil-transmitted nematodes are controlled by administration of a single dose of a benzimidazole to targeted populations in endemic areas. While effective, people are quickly re-infected, necessitating frequent treatment. Widespread exposure to anthelmintic drugs can place significant selective pressure on parasitic nematodes to generate resistance, which has severely compromised benzimidazole anthelmintics for control of livestock nematodes in many areas of the world. Here we report, to our knowledge, the first naturally occurring multidrug-resistant strain of the canine hookworm Ancylostoma caninum. We reveal that this isolate is resistant to fenbendazole at the clinical dosage of 50 mg/kg for 3 days. Our data shows that this strain harbors a fixed, single base pair mutation at amino acid 167 of the β-tubulin isotype 1 gene, and by using CRISPR/Cas9 we demonstrate that introduction of this mutation into the corresponding amino acid in the orthologous β-tubulin gene of Caenorhabditis elegans confers a similar level of resistance to thiabendazole. We also show that the isolate is resistant to the macrocyclic lactone anthelmintic ivermectin. Understanding the mechanism of anthelmintic resistance is important for rational design of control strategies to maintain the usefulness of current drugs, and to monitor the emergence of resistance. The isolate we describe represents the first multidrug-resistant strain of A. caninum reported, and our data reveal a resistance marker that can emerge naturally in response to heavy anthelminthic treatment.

Molecular detection of benzimidazole resistance levels associated with F167Y and F200Y polymorphisms in Haemonchus contortus of goats from Mozambique

Benzimidazole (BZ) resistance of Haemonchus contortus has been associated with single nucleotide polymorphisms (SNPs) in codons 200 (F200Y) and 167 (F167Y) and, to a lesser extent, in codon E198A, of the β-tubulin isotype 1 gene. The present study was undertaken to survey the status of BZ resistance in naturally infected goats in smallholder farms in southern Mozambique by real-time PCR (qPCR) using TaqMan® assays. H. contortus-infective larvae (L3; n = 432) from 12 populations were individually genotyped for F200Y and F167Y SNPs to detect BZ resistance. For the F200Y SNP, the results revealed an overall mean percentages of 18.8% homozygous resistant (RR), 47.8% homozygous susceptible (SS) and 33.4% heterozygous (RS) H. contortus. For the F167Y SNP, the overall mean percentages were 1.6% RR, 94.9% SS and 3.5% RS. The percentage of resistant alleles (%R) for the F200Y and F167Y SNPs was 35.7 and 3.4%, respectively. Genotype combinations of the two mutations indicate resistant percentages ranging from 0.0 to 52.9%. From the four herds with high RR individuals, three farms dewormed the animals monthly, while the fourth farm dewormed the animals every 3 months. In farms where animals were dewormed every 6 months, low percentages of RR individuals were found, whereas no RR individuals were discovered in herds where animals were dewormed annually. These results suggest that the F200Y SNP is more significant in BZ resistance development of the surveyed population compared with the F167Y SNP.

Emergence and the spread of the F200Y benzimidazole resistance mutation in Haemonchus contortus and Haemonchus placei from buffalo and cattle

Benzimidazoles have been intensively (for over 40 years) used in the livestock sector, particularly in small ruminants. This has been led to the widespread emergence of resistance in a number of small ruminant parasite species, especially Haemonchus contortus. In many countries benzimidazole resistance has severely compromised the control of H. contortus in small ruminants; but there is a little information on benzimidazole resistance in H. contortus infecting buffalo and cattle. Resistance to benzimidazoles have also been reported in the large ruminant parasite, Haemonchus placei, but again there is relatively little information on its prevalence. Hence it is very important to understand how resistance-conferring mutations emerge and spread in both parasites in buffalo and cattle hosts in order to develop approaches for the recognition of the problem at an early stage of its development. The present study suggests that the F200Y (TAC) mutation is common in H. contortus, being detected in 5/7 populations at frequencies between 7 and 57%. Furthermore, 6/10 H. placei populations contained the F200Y (TAC) mutation, albeit at low frequencies of between 0.4 and 5%. The phylogenetic analysis suggests that the F200Y (TAC) mutation in H. contortus has emerged on multiple occasions in the region, with at least three independent emergences across the populations. In contrast, the F200Y (TAC) resistance-conferring mutation in H. placei is only seen on a single haplotype. A high level frequency of the resistance haplotypes in the region, suggests that the unique resistance conferring-mutation has spread from a single emergence; likely by anthropogenic animal movement. Overall, these results provide the first clear genetic evidence for the spread of benzimidazole resistance-conferring mutations to multiple different locations from a single emergence in H. placei; while being consistent with previous small ruminant-based observations of multiple emergence of resistance mutations in H. contortus.

Molecular analysis of benzimidazole-resistance associated SNPs in Haemonchus contortus populations of Uruguay

Haemonchus contortus is one of the most important parasite nematodes of small ruminants around the world and causes great economic losses in livestock production. Control of gastrointestinal nematode infections, like haemonchosis, relies mainly on anthelmintic drugs, but its excessive and inappropriate use has caused serious drug resistance issues in many countries, including Uruguay, where sheep production occupies an important place in the country's economy. Benzimidazole (BZ) anthelmintics have been used for decades to treat sheep against H. contortus infection and resistance to this anthelmintic group has been widely described. Molecularly, BZ resistance in H. contortus has been correlated with single nucleotide polymorphisms (SNPs) of the β-tubulin isotype 1 gene at codon 200 and 167 (both TTC to TAC, F167Y and F200Y) and at codon 198 (GAA to GCA, E198A).The aim of this work was to explore the presence of these tubulin SNPs in H. contortus adult worms recovered from sheep abomasa from a slaughterhouse in Uruguay. The mean resistant allelic frequencies at positions F167Y and F200Y were 20.25 and 47.45%, respectively, for worms recovered from naturally infected sheep slaughtered in 2013, while those that were slaughtered in 2014 presented only F200Y SNP with a frequency of 86.89%. Also H. contortus Kirby adult worms (anthelmintic- susceptible McMaster isolate), recovered from artificially infected sheep, were analyzed as reference for comparative purposes This analysis showed susceptible genotype at 167 and 198 position, and a low level of the resistance allele at the 200 position (3.66%). This is the first study for the presence of SNPs in the isotype-1 β-tubulin gene of H. contortus populations in Uruguay, which is consistent with the previous epidemiological studies carried out through the method of fecal egg count reduction test (FECRT), thus confirming the serious resistance levels to BZ anthelmintics also in this country.

Development of a tetra-primer ARMS-PCR for detecting the E198A SNP in the isotype-1 β-tubulin gene of Haemonchus contortus populations in China

The tetra-primer ARMS-PCR is a rapid, simple and low cost method for single nucleotide polymorphism (SNP) genotyping and has been used to detect SNPs associated with diseases and drug resistance. E198A in the isotype-1 β-tubulin gene is one of the three SNPs associated with benzimidazole resistance in parasitic nematode Haemonchus contortus. However, up to now, only PCR-RFLP method was used to test E198A in H. contortus. In the present study, we developed a tetra-primer ARMS-PCR to detect E198A in H. contortus and the accuracy of the results was compared with that of PCR-coupled sequencing. The results showed that optimization of PCR reaction system, especially the proportion of the amount of inner and outer primers, could achieve desirable amplification effect. Three different profiles displaying three distinct genotypes could be identified clearly and intuitively on the agarose gel where the samples with amplified PCR products containing two bands of 433 bp and 200 bp in size indicated susceptible homozygous (SS), those with PCR products containing two bands of 433 bp and 284 bp in length indicated resistant homozygous (RR) and the samples with amplified PCR products containing three bands of 433 bp, 284 bp and 200 bp in size indicated heterozygous (RS). The results showed that the established method can be successfully applied to the detection of E198A in H. contortus, which has high accuracy and is easy to perform.

RNA-Seq de novo assembly and differential transcriptome analysis of the nematode Ascaridia galli in relation to in vivo exposure to flubendazole

The nematode Ascaridia galli (order Ascaridida) is an economically important intestinal parasite responsible for increased food consumption, reduced performance and elevated mortality in commercial poultry production. This roundworm is an emerging problem in several European countries on farms with laying hens, as a consequence of the recent European Union (EU) ban on conventional battery cages. As infection is associated with slow development of low levels of acquired protective immunity, parasite control relies on repeated use of dewormers (anthelmintics). Benzimidazoles (BZ) are currently the only anthelmintic registered in the EU for use in controlling A. galli and there is an obvious risk of overuse of one drug class, selecting for resistance. Thus we developed a reference transcriptome of A. galli to investigate the response in gene expression before and after exposure to the BZ drug flubendazole (FLBZ). Transcriptional variations between treated and untreated A. galli showed that transcripts annotated as mitochondrial glutamate dehydrogenase and cytochrome P450 were significantly down-regulated in treated worms, whereas transcripts homologous to heat shock proteins (HSP), catalase, phosphofructokinase, and a multidrug resistance P-glycoprotein (PGP1) were significantly up-regulated in treated worms. Investigation of candidate transcripts responsible for anthelmintic resistance in livestock nematodes led to identification of several tubulins, including six new isoforms of beta-tubulin, and several ligand-gated ionotropic receptors and ABC-transporters. We discovered several transcripts associated with drug binding and processing genes, but further characterisation using a larger set of worms exposed to BZs in functional assays is required to determine how these are involved in drug binding and metabolism.

Recent Research Progress in China on Haemonchus contortus

Haemonchus contortus is one of the most important parasites of ruminants with worldwide distribution that can bring huge economic losses to the breeding industry of cattle, sheep, and goats. In recent 20 years, studies on H. contortus in China mainly focused on the epidemiology, population genetics, anthelmintic resistance, structural and functional studies of important genes regulating the development of this parasite, interaction between parasite molecules and host cells and vaccine development against haemonchosis, and achieved good progress. However, there is no systematic review about the studies by Chinese researchers on H. contortus in China. The purpose of this review is to bring together the findings from the studies on H. contortus in China in order to obtain the knowledge gained from the recent studies in China and provide foundation for identifying future research directions to establish novel diagnostic methods, discover new drug targets and vaccine candidates for use in preventing and controlling H. contortus in China.

Genotypic and phenotypic evaluation for benzimidazole resistance or susceptibility in Haemonchus contortus isolates

Haemonchus contortus isolates were evaluated for benzimidazole (BZ) resistance or susceptibility by allele-specific PCR based on β-tubulin isotype 1 gene polymorphisms at the F167Y, E198A, and F200Y sites. Two isolates, one presumed susceptible from wild pronghorn antelope (PH) and one known to be resistant from goats (VM), were also assayed phenotypically for BZ resistance or susceptibility in the larval development assay (Drenchrite®). The BZ EC50 was 0.198 μM (intermediate between susceptible and weak resistant) for PH with critical well 5 (intermediate between susceptible and weak resistant) and 1.456 μM (intermediate weak resistant and resistant) for VM with critical well 8.5 (resistant). Genotypically, DNA extracted from pooled VM L3 larvae in the Drenchrite® wells with the highest BZ concentration was homozygous susceptible (SS) at the F167Y and E198A sites and homozygous resistant (RR) at the F200Y site by PCR, and sequence analysis bore this out. PH L3 larvae DNA from a control well (no BZ) was SS at all three sites by PCR, confirmed by sequence analysis. All single adult worm samples (N = 21) from PH, VM, Egypt goat (EG), and a Texas llama were SS at F167Y and E198A by PCR; however, only 3 PH worms and 1 EG worm were SS at F200Y. Three additional PH worms were RS and upon cloning two clones were identified as resistant by sequencing and two as susceptible. Clones from single adult worms VM, llama, and EG samples that were RR by PCR at F200Y were sequence verified as resistant. In this study, F200Y was the most frequently found genotypic marker for BZ resistance or susceptibility in the different Haemonchus isolates.

Two benzimidazole resistance-associated SNPs in the isotype-1 β-tubulin gene predominate in Haemonchus contortus populations from eight regions in China

Haemonchus contortus is one of the most important parasitic nematodes of small ruminants around the world, particularly in tropical and subtropical regions. The control of haemonchosis relies mainly on anthelmintics, but the excessive and prolonged use of anthelmintics is causing serious drug resistance issues in many countries. As benzimidazole (BZ) anthelmintics have been broadly used in China, we hypothesized that resistance is widespread. Given the link between three known single nucleotide polymorphisms (SNPs, designated F167Y, E198A and F200Y) in the isotype-1 β-tubulin gene and BZ resistance, our goal here was to explore the presence of these mutations in H. contortus from small ruminants (sheep and goats) from eight provinces in China using PCR-coupled sequencing. In addition, the genetic diversity and genetic relationship of isotype-1 β-tubulin sequence haplotypes were also investigated. Among 192 H. contortus adult individuals representing the eight populations, we identified six distinct sequence types, five of which had SNP E198A (GCA) and/or F200Y (TAC). Sequence analysis showed that the frequencies of SNPs E198A and F200Y were 0-70% and 0-31%, respectively. SNP F167Y (TAC) was not detected in any population. In addition, high haplotype diversities (0.455-0.939) and nucleotide diversities (0.018-0.039) were calculated. A network analysis of the isotype-1 β-tubulin gene sequences showed that SNPs E198A and F200Y occurred in multiple distinct groupings, suggesting multiple independent origins of these SNPs. The findings of this first study of SNPs in the isotype-1 β-tubulin gene of H. contortus populations suggest that BZ resistance is prevalent in some regions of China, and that any control strategy might focus on monitoring BZ resistance in this country.

Benzimidazole resistance in helminths: From problem to diagnosis

Helminth parasites cause significant morbidity and mortality in endemic countries. Given the severity of symptoms that helminths may elicit in the host, intervention with prophylactic and therapeutic measures is imperative. Treatment with benzimidazoles is the most widely used means of combatting these parasites. However, widespread use of these drugs can select for drug-resistant parasite strains. In this review, we approach the problem of benzimidazole resistance in helminths in both humans and animals, focusing on the properties of the drug, the molecular mechanisms of drug resistance and how resistance is diagnosed.

Genotyping of benzimidazole resistant and susceptible isolates of Haemonchus contortus from sheep by allele specific PCR

Extensive and indiscriminate use of the benzimidazole class of drugs has led to the onset of anthelmintic resistance. In tropical countries like India, Haemonchus contortus is the most pathogenic parasite infecting sheep and goats. The widespread presence of resistant helminths (especially H. contortus) threatens the livestock farming. The use of various drugs has led to single nucleotide polymorphism that causes specific amino acid substitutions in β-tubulin protein of H. contortus to confer resistance. This emphasizes the need for a survey on the present status of resistance in India. In this study, allele specific PCR was employed to screen the presence of a SNP, a thymine-to-adenine transversion which leads to substitution of amino acid in codon 200 of β-tubulin gene that is correlated specifically with BZ resistance. Third stage larvae (L3) from pooled faecal cultures of four organized sheep farms served as a source of genomic DNA for identification of H. contortus and further genotype analysis. A total of 1000 larvae was screened, out of which 673 larvae were identified as H. contortus. Among 673 H. contortus larvae, 539 larvae (80 %) were genotyped as homozygous resistant (rr) and remaining 134 (20 %) were heterozygous susceptible (Sr) by allele specific PCR. The concluded resistance status reasons out the failure of anthelmintic drug in treating ruminants. Immediate steps are needed to avoid further aggravation of the problem. Target selective treatment by reviewing the resistance status of individual drugs, appropriate use of anthelmintic drugs and other control strategies will provide a pragmatic option for delaying the further spread of anthelmintic resistance.

Sympatric species distribution, genetic diversity and population structure of Haemonchus isolates from domestic ruminants in Pakistan

Haemonchus species are major gastro-intestinal parasites affecting ruminants across the world. The present study aimed to assess the sympatric species distribution, genetic diversity, population structure and frequency of β-tubulin isotype 1 alleles associated with benzimidazole resistance. Internal transcribed spacer 2 (ITS2) sequences revealed three sympatric species of Haemonchus, H. contortus, H. placei and H. longistipes with 12 distinct genotypes circulating among ruminant hosts in Pakistan. High genetic variability was observed in Pakistani Haemonchus isolates at nicotine amide dehydrogenase subunit 4 (ND4) and cytochrome oxidase subunit 1 (COI) gene loci. Intra-population diversity parameters were higher in H. contortus isolates than H. placei. Phylogenetic analysis of ND4 and COI sequences did not reveal clustering of haplotypes originating from a particular host indicating high rate of gene flow among Haemonchus parasites infecting sheep, goat and cattle in Pakistan. ND4 and COI haplotypes from Pakistan were compared to sequences of Haemonchus isolates from 11 countries to elucidate the population structure. Multidimensional scaling (MDS) plot of pairwise FST derived from 531 ND4 haplotypes revealed clustering together of H. contortus from Pakistan, China, Malaysia and Italy while the isolates from Yemen and United States were found to be genetically distinct. With respect to H. placei, isolates from Pakistan were found to be genetically differentiated from isolates of other countries. The tests for selective neutrality revealed negative D statistics and did not reveal significant deviations in Pakistani Haemonchus populations while significant deviation (P < 0.05) was observed in Brazilian and Chinese H. contortus populations. Median Joining (MJ) network of ND4 haplotypes revealed Yemenese H. contortus being closer to H. placei cluster. β-tubulin isotype 1 genotyping revealed 7.86% frequency of Y allele associated with benzimidazole resistance at F200Y locus in Pakistani Haemonchus isolates.

Understanding Haemonchus contortus Better Through Genomics and Transcriptomics

Parasitic roundworms (nematodes) cause substantial mortality and morbidity in animals globally. The barber's pole worm, Haemonchus contortus, is one of the most economically significant parasitic nematodes of small ruminants worldwide. Although this and related nematodes can be controlled relatively well using anthelmintics, resistance against most drugs in common use has become a major problem. Until recently, almost nothing was known about the molecular biology of H. contortus on a global scale. This chapter gives a brief background on H. contortus and haemonchosis, immune responses, vaccine research, chemotherapeutics and current problems associated with drug resistance. It also describes progress in transcriptomics before the availability of H. contortus genomes and the challenges associated with such work. It then reviews major progress on the two draft genomes and developmental transcriptomes of H. contortus, and summarizes their implications for the molecular biology of this worm in both the free-living and the parasitic stages of its life cycle. The chapter concludes by considering how genomics and transcriptomics can accelerate research on Haemonchus and related parasites, and can enable the development of new interventions against haemonchosis.

The emergence of resistance to the benzimidazole anthlemintics in parasitic nematodes of livestock is characterised by multiple independent hard and soft selective sweeps

Anthelmintic resistance is a major problem for the control of parasitic nematodes of livestock and of growing concern for human parasite control. However, there is little understanding of how resistance arises and spreads or of the "genetic signature" of selection for this group of important pathogens. We have investigated these questions in the system for which anthelmintic resistance is most advanced; benzimidazole resistance in the sheep parasites Haemonchus contortus and Teladorsagia circumcincta. Population genetic analysis with neutral microsatellite markers reveals that T. circumcincta has higher genetic diversity but lower genetic differentiation between farms than H. contortus in the UK. We propose that this is due to epidemiological differences between the two parasites resulting in greater seasonal bottlenecking of H. contortus. There is a remarkably high level of resistance haplotype diversity in both parasites compared with drug resistance studies in other eukaryotic systems. Our analysis suggests a minimum of four independent origins of resistance mutations on just seven farms for H. contortus, and even more for T. circumincta. Both hard and soft selective sweeps have occurred with striking differences between individual farms. The sweeps are generally softer for T. circumcincta than H. contortus, consistent with its higher level of genetic diversity and consequent greater availability of new mutations. We propose a model in which multiple independent resistance mutations recurrently arise and spread by migration to explain the widespread occurrence of resistance in these parasites. Finally, in spite of the complex haplotypic diversity, we show that selection can be detected at the target locus using simple measures of genetic diversity and departures from neutrality. This work has important implications for the application of genome-wide approaches to identify new anthelmintic resistance loci and the likelihood of anthelmintic resistance emerging as selection pressure is increased in human soil-transmitted nematodes by community wide treatment programs.

Recent advances in candidate-gene and whole-genome approaches to the discovery of anthelmintic resistance markers and the description of drug/receptor interactions

Anthelmintic resistance has a great impact on livestock production systems worldwide, is an emerging concern in companion animal medicine, and represents a threat to our ongoing ability to control human soil-transmitted helminths. The Consortium for Anthelmintic Resistance and Susceptibility (CARS) provides a forum for scientists to meet and discuss the latest developments in the search for molecular markers of anthelmintic resistance. Such markers are important for detecting drug resistant worm populations, and indicating the likely impact of the resistance on drug efficacy. The molecular basis of resistance is also important for understanding how anthelmintics work, and how drug resistant populations arise. Changes to target receptors, drug efflux and other biological processes can be involved. This paper reports on the CARS group meeting held in August 2013 in Perth, Australia. The latest knowledge on the development of molecular markers for resistance to each of the principal classes of anthelmintics is reviewed. The molecular basis of resistance is best understood for the benzimidazole group of compounds, and we examine recent work to translate this knowledge into useful diagnostics for field use. We examine recent candidate-gene and whole-genome approaches to understanding anthelmintic resistance and identify markers. We also look at drug transporters in terms of providing both useful markers for resistance, as well as opportunities to overcome resistance through the targeting of the transporters themselves with inhibitors. Finally, we describe the tools available for the application of the newest high-throughput sequencing technologies to the study of anthelmintic resistance.

Resistência anti-helmíntica em nematoides gastrintestinais de pequenos ruminantes: avanços e limitações para seu diagnóstico

A seleção e a crescente disseminação de nematoides resistentes aos anti-helmínticos mais comumente utilizados, benzimidazóis (BZs), imidazotiazóis e lactonas macrocíclicas (LMs), constituem um sério entrave na produção de pequenos ruminantes em todo o mundo. O uso de métodos eficientes e sensíveis para a detecção e o monitoramento da resistência anti-helmíntica no campo torna-se urgente, especialmente para os grupos de BZs e LMs, devido aos constantes relatos de resistência. A obtenção de um diagnóstico preciso e precoce da resistência é extremamente importante para auxiliar a tomada de decisão em programas de controle parasitário, com o objetivo de preservar a vida útil dos produtos e limitar o desenvolvimento da resistência nas populações de nematoides. Os testes in vivo e, mais recentemente, os testes in vitro têm sido desenvolvidos para a detecção de nematoides resistentes aos principais grupos de anti-helmínticos. No entanto, a disponibilidade de testes in vitro validados e o seu uso prático ainda são muito limitados. Embora o teste de redução na contagem de ovos nas fezes (TRCOF, in vivo - indireto) seja o principal método de escolha para a detecção de resistência no campo, vem recebendo críticas quanto à validade dos resultados, e passa por significativas modificações. Além disso, o desenvolvimento de técnicas moleculares a partir de alterações genômicas gerou avanços consideráveis nessa área de investigação, com o uso de mutações nos códons 167, 198 e 200 do gene da β-tubulina como principais SNPs (polimorfismos de nucleotídeo único; do inglês Single Nucleotide Polymorphisms) associados à resistência aos BZs. A presente revisão tem o objetivo de discutir os métodos de diagnóstico disponíveis para a detecção de resistência anti-helmíntica em nematoides de pequenos ruminantes, destacando progressos e obstáculos para seu uso na rotina laboratorial e no campo.

Molecular basis for benzimidazole resistance from a novel β-tubulin binding site model

Benzimidazole-2-carbamate derivatives (BzCs) are the most commonly used antiparasitic drugs for the treatment of protozoan and helminthic infections. BzCs inhibit the microtubule polymerization mechanism through binding selectively to the β-tubulin subunit in which mutations have been identified that lead to drug resistance. Currently, the lack of crystallographic structures of β-tubulin in parasites has limited the study of the binding site and the analysis of the resistance to BzCs. Recently, our research group has proposed a model to explain the interaction between the BzCs and a binding site in the β-tubulin. Herein, we report the homology models of two susceptible (Haemonchus contortus and Giardia intestinalis) parasites and one unsusceptible (Entamoeba histolytica) generated using the structure of the mammal Ovis aries, considered as a low susceptible organism, as a template. Additionally, the mechanism by which the principal single point mutations Phe167Tyr, Glu198Ala and Phe200Tyr could lead to resistance to BzCs is analyzed. Molecular docking and molecular dynamics studies were carried out in order to evaluate the models and the ligand-protein complexes' behaviors. This study represents a first attempt towards understanding, at the molecular level, the structural composition of β-tubulin in all organisms, also suggesting possible resistance mechanisms. Furthermore, these results support the importance of benzimidazole derivative optimization in order to design more potent and selective (less toxic) molecules for the treatment of parasitic diseases.

Characterization and comparative analysis of the complete Haemonchus contortus β-tubulin gene family and implications for benzimidazole resistance in strongylid nematodes

Parasitic nematode β-tubulin genes are of particular interest because they are the targets of benzimidazole drugs. However, in spite of this, the full β-tubulin gene family has not been characterized for any parasitic nematode to date. Haemonchus contortus is the parasite species for which we understand benzimidazole resistance the best and its close phylogenetic relationship with Caenorhabditis elegans potentially allows inferences of gene function by comparative analysis. Consequently, we have characterized the full β-tubulin gene family in H. contortus. Further to the previously identified Hco-tbb-iso-1 and Hco-tbb-iso-2 genes, we have characterized two additional family members designated Hco-tbb-iso-3 and Hco-tbb-iso-4. We show that Hco-tbb-iso-1 is not a one-to-one orthologue with Cel-ben-1, the only β-tubulin gene in C. elegans that is a benzimidazole drug target. Instead, both Hco-tbb-iso-1 and Hco-tbb-iso-2 have a complex evolutionary relationship with three C. elegans β-tubulin genes: Cel-ben-1, Cel-tbb-1 and Cel-tbb-2. Furthermore, we show that both Hco-tbb-iso-1 and Hco-tbb-iso-2 are highly expressed in adult worms; in contrast, Hco-tbb-iso-3 and Hco-tbb-iso-4 are expressed only at very low levels and are orthologous to the Cel-mec-7 and Cel-tbb-4 genes, respectively, suggesting that they have specialized functional roles. Indeed, we have found that the expression pattern of Hco-tbb-iso-3 in H. contortus is identical to that of Cel-mec-7 in C. elegans, being expressed in just six "touch receptor" mechano-sensory neurons. These results suggest that further investigation is warranted into the potential involvement of strongylid isotype-2 β-tubulin genes in mechanisms of benzimidazole resistance.

Relative level of thiabendazole resistance associated with the E198A and F200Y SNPs in larvae of a multi-drug resistant isolate of Haemonchus contortus

While the F200Y SNP in the beta-tubulin gene is most commonly associated with benzimidazole resistance in trichostrongylid nematodes, other SNPs as well as drug efflux pathways have been implicated in the resistance. The relative contributions of all these mechanisms are not understood sufficiently to allow expected drug efficacy to be inferred from molecular data. As a component of developing better means to interpret molecular resistance tests, the present study utilised a drug resistant Haemonchus contortus isolate which possesses two of the principal benzimidazole resistance SNPs (E198A and F200Y) in order to assess the relative degree of resistance conferred by the two SNPs. We exposed larvae to a range of thiabendazole concentrations in in vitro development assays, and collected the surviving L3 larvae at each drug concentration to establish sub-populations showing increasing levels of resistance. We then sequenced the isotype 1 beta-tubulin gene in pooled larval samples, and measured allele frequencies at the two SNP positions. The frequency of the resistance allele at the 198 position increased as the thiabendazole concentration increased, while the frequency of the resistance allele at the 200 position decreased. Genotyping of individual larvae showed that the highest drug concentration was associated with the removal of all genotypes except for homozygous resistance at the 198 position alongside homozygous susceptible at the 200 position. This indicates that, at least for larval life stages, the E198A SNP is able to confer higher levels of resistance to benzimidazole drugs than the F200Y SNP, and that the homozygosity at 198 in the highly resistant individuals is mutually exclusive with heterozygosity or resistant homozygosity at the 200 position. This study illustrates the need to understand the relative contributions of different resistance mechanisms in order to maximise the degree to which molecular tests are able to inform on drug resistance phenotype.

Effects of infection intensity with Strongyloides papillosus and albendazole treatment on development of oxidative/nitrosative stress in sheep

The objective of this study was to estimate and evaluate oxidative/nitrosative stress parameters in sheep infected with Strongyloides papillosus and after antihelminthic treatment with albendazole (ABZ). This parasite, especially during development stages can seriously damage parenchaematous organs during migration within the host. The presence of parasites leads to increased productions of reactive oxygen species (ROS) and reactive nitrogen species (RNS). It is also well known that certain drugs can be very harmful for the delicate oxidant-antioxidant equilibrium, provoking oxidative stress during their biotransformation. ABZ is a broad spectrum antihelminthic drug, frequently used in veterinary medicine for therapy of parasitic infections. The current research was performed on female Württemberg sheep (n=48). The distribution of parasites in sheep was evaluated using the native smear coprological technique, by sedimentation and flotation methods, revealing the presence of S. papillosus. The degree of infection intensity per sheep was quantitatively established by the method of McMaster, the animals having been divided into three groups according to the intensity of infection; mild, moderate and high. The control group consisted of sheep negative to the parasites. After determining the type of parasite infection, the sheep were treated with ABZ, per orally, in single doses of 5mg/kg per body weight. Sampling of feces for parasitological and blood for biochemical assaying was performed on the 0 and 21st day after treatment with ABZ. The oxidative stress parameters were measured for catalase activity (CAT), the red cell membrane damage by level of malondialdehyde (MDA), while carbonyl and thiol plasma protein group concentrations were used as indicators of the degree of protein oxidative modification. The activity of Cu,Zn-superoxide dismutase (SOD) and relative distribution of lactate dehydrogenase (LDH(1)-LDH(5)) activity were determined electrophoretically. The distribution of LDH isoenzymes in sheep moderately and highly infected with S. papillosus revealed that the parasite induced damage to the myocardial (LDH(2)), lung (LDH(3)) and liver cells (LDH(5)) in infected animals, while ABZ treatment only damaged liver cells (LDH(5)). The MDA concentration revealed that lipid peroxidation increased both in the presence of parasites and the antihelminthic formulation tested (p<0.001) when compared to the control sheep, while the increase of carbonyl concentration (p<0.001), as well as the observed decrease of thiol concentration (p<0.001) indicated significant oxidative damage of plasma proteins in experimental sheep, when compared to the control animals. Our results indicate that S. papillosus induces oxidative/nitrosative stress in sheep. The antihelminthic treatment with ABZ further promotes the disbalance of oxidative-antioxidative equilibrium in all tested sheep.

In vitro effect of seven antiparasitics on Acolpenteron ureteroecetes (Dactylogyridae) from largemouth bass Micropterus salmoides (Centrarchidae)

Few drugs are approved by the United States Food and Drug Administration for treating parasite infections in minor species such as fish, due in part to the high cost of developing such drugs and to a relatively small market share for drug sponsors. Because in vivo effectiveness trials for antiparasitic drugs are costly, time consuming, and use many animals, a systematic in vitro screening approach to describe parasite motility could help find promising drug candidates. We evaluated the effects of 7 antiparasitics on the activity and survival of the endoparasitic monogenean Acolpenteron ureteroecetes (Dactylogyridae) collected from the posterior kidneys of juvenile largemouth bass Micropterus salmoides (Lacepede, 1802) (Centrarchidae) held in the laboratory. Tests were conducted in 12 well tissue culture plates; each well had 3 parasites, and we tested 3 concentrations and 1 control for each of the 7 antiparasitics. The parasites were observed immediately after adding the drug, at 1 to 3 h, and 17 to 26 h, and video recordings were made. Drug effects were recorded by documenting morbidity (reduced movement, tremors, contracted body, abnormal morphology) and mortality. A. ureteroecetes was strongly affected by the quinoline praziquantel, the imidazothiazide levamisole, and the organophosphates dichlorvos and trichlorfon. The parasites were moderately affected by the macrocyclic lactones ivermectin and emamectin, and generally unaffected by the benzimidazole mebendazole. Our study demonstrates the utility of characterizing in vitro responses with video microscopy to document responses of fish parasites for initial screens of drug effects on a fish monogenean.

A novel high throughput assay for anthelmintic drug screening and resistance diagnosis by real-time monitoring of parasite motility

Background

Helminth parasites cause untold morbidity and mortality to billions of people and livestock. Anthelmintic drugs are available but resistance is a problem in livestock parasites, and is a looming threat for human helminths. Testing the efficacy of available anthelmintic drugs and development of new drugs is hindered by the lack of objective high-throughput screening methods. Currently, drug effect is assessed by observing motility or development of parasites using laborious, subjective, low-throughput methods.

Methodology/Principal Findings

Here we describe a novel application for a real-time cell monitoring device (xCELLigence) that can simply and objectively assess anthelmintic effects by measuring parasite motility in real time in a fully automated high-throughput fashion. We quantitatively assessed motility and determined real time IC₅₀ values of different anthelmintic drugs against several developmental stages of major helminth pathogens of humans and livestock, including larval Haemonchus contortus and Strongyloides ratti, and adult hookworms and blood flukes. The assay enabled quantification of the onset of egg hatching in real time, and the impact of drugs on hatch rate, as well as discriminating between the effects of drugs on motility of drug-susceptible and –resistant isolates of H. contortus.

Conclusions/Significance

Our findings indicate that this technique will be suitable for discovery and development of new anthelmintic drugs as well as for detection of phenotypic resistance to existing drugs for the majority of helminths and other pathogens where motility is a measure of pathogen viability. The method is also amenable to use for other purposes where motility is assessed, such as gene silencing or antibody-mediated killing.

Parasitic worms cause untold morbidity and mortality on billions of people and livestock. Drugs are available but resistance is problematic in livestock parasites and is a looming threat for human helminths. Currently, new drug discovery and resistance monitoring is hindered as drug efficacy is assessed by observing motility or development of parasites using laborious, subjective, low-throughput methods evaluated by eye using microscopy. Here we describe a novel application for a cell monitoring device (xCELLigence) that can simply and objectively assess real time anti-parasite efficacy of drugs on eggs, larvae and adults in a fully automated, label-free, high-throughput fashion. This technique overcomes the current low-throughput bottleneck in anthelmintic drug development and resistance detection pipelines. The widespread use of this device to screen for new therapeutics or emerging drug resistance will be an invaluable asset in the fight against human, animal and plant parasitic helminths and other pathogens that plague our planet.

Anthelmintic resistance: markers for resistance, or susceptibility?

The Consortium for Anthelmintic Resistance and Susceptibility (CARS) brings together researchers worldwide, with a focus of advancing knowledge of resistance and providing information on detection methods and treatment strategies. Advances in this field suggest mechanisms and features of resistance that are shared among different classes of anthelmintic. Benzimidazole resistance is characterized by specific amino acid substitutions in beta-tubulin. If present, these substitutions increase in frequency upon drug treatment and lead to treatment failure. In the laboratory, sequence substitutions in ion-channels can contribute to macrocyclic lactone resistance, but there is little evidence that they are significant in the field. Changes in gene expression are associated with resistance to several different classes of anthelmintic. Increased P-glycoprotein expression may prevent drug access to its site of action. Decreased expression of ion-channel subunits and the loss of specific receptors may remove the drug target. Tools for the identification and genetic analysis of parasitic nematodes and a new online database will help to coordinate research efforts in this area. Resistance may result from a loss of sensitivity as well as the appearance of resistance. A focus on the presence of anthelmintic susceptibility may be as important as the detection of resistance.