Similarities and differences in the biotransformation and transcriptomic responses of Caenorhabditis elegans and Haemonchus contortus to five different benzimidazole drugs

Systems biology of Haemonchus contortus - Advancing biotechnology for parasitic nematode control

Parasitic nematodes represent a substantial global burden, impacting animal health, agriculture and economies worldwide. Of these worms, Haemonchus contortus - a blood-feeding nematode of ruminants - is a major pathogen and a model for molecular and applied parasitology research. This review synthesises some key advances in understanding the molecular biology, genetic diversity and host-parasite interactions of H. contortus, highlighting its value for comparative studies with the free-living nematode Caenorhabditis elegans. Key themes include recent developments in genomic, transcriptomic and proteomic technologies and resources, which are illuminating critical molecular pathways, including the ubiquitination pathway, protease/protease inhibitor systems and the secretome of H. contortus. Some of these insights are providing a foundation for identifying essential genes and exploring their potential as targets for novel anthelmintics or vaccines, particularly in the face of widespread anthelmintic resistance. Advanced bioinformatic tools, such as machine learning (ML) algorithms and artificial intelligence (AI)-driven protein structure prediction, are enhancing annotation capabilities, facilitating and accelerating analyses of gene functions, and biological pathways and processes. This review also discusses the integration of these tools with cutting-edge single-cell sequencing and spatial transcriptomics to dissect host-parasite interactions at the cellular level. The discussion emphasises the importance of curated databases, improved culture systems and functional genomics platforms to translate molecular discoveries into practical outcomes, such as novel interventions. New research findings and resources not only advance research on H. contortus and related nematodes but may also pave the way for innovative solutions to the global challenges with anthelmintic resistance.

Biotransformation of anthelmintics in nematodes in relation to drug resistance

In all organisms, the biotransformation of xenobiotics to less toxic and more hydrophilic compounds represents an effective defense strategy. In pathogens, the biotransformation of drugs (used for their elimination from the host) may provide undesirable protective effects that could potentially compromise the drug's efficacy. Accordingly, increased drug deactivation via accelerated biotransformation is now considered as one of the mechanisms of drug resistance. The present study summarizes the current knowledge regarding the biotransformation of anthelmintics, specifically drugs used to treat mainly nematodes, a group of parasites that are a significant health concern for humans and animals. The main biotransformation enzymes are introduced and their roles in anthelmintics metabolism in nematodes are discussed with a particular focus on their potential participation in drug resistance. Similarly, the inducibility of biotransformation enzymes with sublethal doses of anthelmintics is presented in view of its potential contribution to drug resistance development. In the conclusion, the main tasks awaiting scientists in this area are outlined.

Differences in constitutive gene expression of cytochrome P450 enzymes and ATP-binding cassette transporter gene expression between a susceptible and a highly macrocyclic lactone-resistant Haemonchus contortus isolate in the absence of drug-inducible expression

BACKGROUND

Anthelmintic resistance in ruminants is a widespread problem that has a severe impact on productivity and animal welfare. The helminth Haemonchus contortus is generally considered the most important parasite in small ruminants due to its high pathogenicity and the widespread occurrence of anthelmintic resistance in it. Although the molecular mechanisms associated with resistance against the anthelmintics benzimidazoles (BZs) and levamisole are relatively well understood, the resistance mechanisms against the widely used anthelmintic macrocyclic lactones (MLs) ivermectin (IVM) and moxidectin (MOX) remain poorly understood. Detoxifying enzymes and xenobiotic transporters have been frequently proposed to play a role in ML resistance in multiple organisms, including nematodes.

METHODS

The reference genome of H. contortus was screened for cytochrome P450 genes (cyp genes) by using the Basic Local Alignment Search Tool, and maximum-likelihood phylogenetic analysis was used to assign the sequences to gene families. Fourth-stage larvae of the susceptible (McMaster) and the ML-resistant (Berlin-selected) H. contortus isolates were generated in vitro and compared regarding basal expression levels of cyp genes and ATP-binding cassette (ABC) transporters by using RNA sequencing. The resistant isolate was further incubated with 100 nM IVM or MOX for 3, 6 and 12 h, and the effects of incubation time and drugs were evaluated.

RESULTS

Twenty-five cyp genes were identified in the H. contortus genome and assigned to 13 different families. The ML-resistant isolate showed significantly higher and lower constitutive expression of 13 and four cyp genes, respectively. Out of the 50 ABC transporter genes, only six showed significantly higher expression in the ML-resistant isolate, while 12 showed lower expression. The fold changes were in general low (range 0.44-5.16). Only pgp-13 showed significant downregulation in response to IVM (0.77 fold change at 6 h, 0.96 fold change at 12 h) and MOX (0.84 fold change at 12 h). In contrast, mrp-5 was significantly, albeit minimally, upregulated in the presence of IVM, but not MOX, after 12 h (1.02 fold change).

CONCLUSIONS

Despite little observable ML-inducible gene expression in the isolate examined here, some of the changes in the baseline expression levels might well contribute to ML resistance in the context of additional changes in a multigenic resistance model. However, neither cyp genes nor the ABC transporters appear to be the main drivers that can explain the high levels of resistance observed in the resistant isolate examined here.

Integrating metabolomics into the diagnosis and investigation of anthelmintic resistance

Anthelmintic resistance (AR) in parasitic nematodes poses a global health problem in livestock and domestic animals and is an emerging problem in humans. Consequently, we must understand the mechanisms of AR, including target-site resistance (TSR), in which mutations affect drug binding, and non-target site resistance (NTSR), which involves alterations in drug metabolism and detoxification processes. Because much of the focus has been on TSR, NTSR has received less attention. Here, we describe how metabolomics approaches using Caenorhabditis elegans offer the ability to disentangle nematode drug metabolism, identify metabolic changes associated with resistance, uncover novel biomarkers, and enhance diagnostic methods.

CYP14 family in Caenorhabditis elegans: Mitochondrial function, detoxification, and lifespan

CYP-14 members of the Caenorhabditis elegans (C. elegans) Cytochrome P450 (CYP) enzyme family, plays important roles in mitochondrial dysfunction, detoxification, lipid metabolism, defense and lifespan regulation. The review identifies CYP-14 members: cyp-14A1, cyp-14A2, cyp-14A3, cyp-14A4, cyp-14A5 and their homology with human CYP families. Despite limited studies on C. elegans cyp-14 members, the findings unraveled their complex crosstalk between mitochondrial stress, detoxification mechanisms, and lifespan regulation, emphasizing the complexity of these interconnected pathways as well as how their regulation depends on environmental cues changes including pH, nutrients, ROS and chemical stressors. The review underscores the translational relevance to human health, shedding light on potential human homologues and their implications in age-related, metabolic and respiratory diseases. Among other genes, cyp-14A2 and cyp-14A4 predominate the mitochondrial function, heat resistance, lipid metabolism, detoxification and lifespan pathways. In conclusion, these insights pave the way for future research, offering promising avenues for therapeutic interventions targeting CYP-14 activity to address age-related diseases and promote healthy aging.

Multiple UDP glycosyltransferases modulate benzimidazole drug sensitivity in the nematode Caenorhabditis elegans in an additive manner

Xenobiotic biotransformation is an important modulator of anthelmintic drug potency and a potential mechanism of anthelmintic resistance. Both the free-living nematode Caenorhabditis elegans and the ruminant parasite Haemonchus contortus biotransform benzimidazole drugs by glucose conjugation, likely catalysed by UDP-glycosyltransferase (UGT) enzymes. To identify C. elegans genes involved in benzimidazole drug detoxification, we first used a comparative phylogenetic analysis of UGTs from humans, C. elegans and H. contortus, combined with available RNAseq datasets to identify which of the 63 C. elegans ugt genes are most likely to be involved in benzimidazole drug biotransformation. RNA interference knockdown of 15 prioritized C. elegans genes identified those that sensitized animals to the benzimidazole derivative albendazole (ABZ). Genetic mutations subsequently revealed that loss of ugt-9 and ugt-11 had the strongest effects. The "ugt-9 cluster" includes these genes, together with six other closely related ugts. A CRISPR-Cas-9 deletion that removed seven of the eight ugt-9 cluster genes had greater ABZ sensitivity than the single largest-effect mutation. Furthermore, a double mutant of ugt-22 (which is not a member of the ugt-9 cluster) with the ugt-9 cluster deletion further increased ABZ sensitivity. This additivity of mutant phenotypes suggest that ugt genes act in parallel, which could have several, not mutually exclusive, explanations. ugt mutations have different effects with different benzimidazole derivatives, suggesting that enzymes with different specificities could together more efficiently detoxify drugs. Expression patterns of ugt-9, ugt-11 and ugt-22 gfp reporters differ and so likely act in different tissues which may, at least in part, explain their additive effects on drug potency. Overexpression of ugt-9 alone was sufficient to confer partial ABZ resistance, indicating increasing total UGT activity protects animals. In summary, our results suggest that the multiple UGT enzymes have overlapping but not completely redundant functions in benzimidazole drug detoxification and may represent "druggable" targets to improve benzimidazole drug potency.

Role of UDP-Glycosyltransferase (ugt) Genes in Detoxification and Glycosylation of 1-Hydroxyphenazine (1-HP) in Caenorhabditis elegans

Caenorhabditis elegans is a useful model organism to study the xenobiotic detoxification pathways of various natural and synthetic toxins, but the mechanisms of phase II detoxification are understudied. 1-Hydroxyphenazine (1-HP), a toxin produced by the bacterium Pseudomonas aeruginosa, kills C. elegans. We previously showed that C. elegans detoxifies 1-HP by adding one, two, or three glucose molecules in N2 worms. Our current study evaluates the roles that some UDP-glycosyltransferase (ugt) genes play in 1-HP detoxification. We show that ugt-23 and ugt-49 knockout mutants are more sensitive to 1-HP than reference strains N2 or PD1074. Our data also show that ugt-23 knockout mutants produce reduced amounts of the trisaccharide sugars, while the ugt-49 knockout mutants produce reduced amounts of all 1-HP derivatives except for the glucopyranosyl product compared to the reference strains. We characterized the structure of the trisaccharide sugar phenazines made by C. elegans and showed that one of the sugar modifications contains an N-acetylglucosamine (GlcNAc) in place of glucose. This implies broad specificity regarding UGT function and the role of genes other than ogt-1 in adding GlcNAc, at least in small-molecule detoxification.

Exploring the role of detoxification genes in the resistance of Bursaphelenchus xylophilus to different exogenous nematicidal substances using transcriptomic analyses

Bursaphelenchus xylophilus (Pine wood nematode, PWN) has become a worldwide forest disease due to its rapid infection ability, high lethality and difficulty in control. The main means of countering B. xylophilus is currently chemical control, but nematicides can present problems such as environmental pollution and drug resistance. The development of novel environmentally-friendly nematicides has thus become a focus of recent research. In this study, BxUGT3 and BxUGT34, which might be related to detoxification, were investigated by comparing transcriptomic and WGCNA approaches. Three other genes with a similar expression pattern, BxUGT13, BxUGT14, and BxUGT16, were found by gene family analysis. Further bioassays and qPCR assays confirmed that these five genes showed significant changes in transcript levels upon exposure to α-pinene and carvone, demonstrating that they respond to exogenous nematicidal substances. Finally, RNAi and bioassays showed that B. xylophilus with silenced BxUGT16 had increased mortality in the face of α-pinene and carvone stress, suggesting that BxUGT16 plays an important role in detoxification. Taken together, this study used novel molecular research methods, explored the detoxification mechanism of B. xylophilus at a transcriptomic level, and revealed a molecular target for the development of novel biopesticides.

Repurposing degradation pathways for modular metabolite biosynthesis in nematodes

Recent studies have revealed that Caenorhabditis elegans and other nematodes repurpose products from biochemical degradation pathways for the combinatorial assembly of complex modular structures that serve diverse signaling functions. Building blocks from neurotransmitter, amino acid, nucleoside and fatty acid metabolism are attached to scaffolds based on the dideoxyhexose ascarylose or glucose, resulting in hundreds of modular ascarosides and glucosides. Genome-wide association studies have identified carboxylesterases as the key enzymes mediating modular assembly, enabling rapid compound discovery via untargeted metabolomics and suggesting that modular metabolite biosynthesis originates from the 'hijacking' of conserved detoxification mechanisms. Modular metabolites thus represent a distinct biosynthetic strategy for generating structural and functional diversity in nematodes, complementing the primarily polyketide synthase- and nonribosomal peptide synthetase-derived universe of microbial natural products. Although many aspects of modular metabolite biosynthesis and function remain to be elucidated, their identification demonstrates how phenotype-driven compound discovery, untargeted metabolomics and genomic approaches can synergize to facilitate the annotation of metabolic dark matter.

Assembly and Analysis of Haemonchus contortus Transcriptome as a Tool for the Knowledge of Ivermectin Resistance Mechanisms

Haemonchus contortus (Hc) is an important parasitic nematode of small ruminants. In this study we assembled the transcriptome of Hc as a model to contribute to the knowledge about the profile of the differential gene expression between two Mexican Hc strains under different anthelmintic resistance statuses, one susceptible and the other resistant to ivermectin (IVMs and IVMr, respectively), in order to improve and/or to have new strategies of control and diagnosis. The transcript sequence reads were assembled and annotated. Overall, ~127 Mbp were assembled and distributed into 77,422 transcript sequences, and 4394 transcripts of the de novo transcriptome were matched base on at least one of the following criteria: (1) Phylum Nemathelminthes and Platyhelminthes, important for animal health care, and (2) ≥55% of sequence identity with other organisms. The gene ontology (GO) enrichment analysis (GOEA) was performed to study the level of gene regulation to IVMr and IVMs strains using Log Fold Change (LFC) filtering values ≥ 1 and ≥ 2. The upregulated-displayed genes obtained via GOEA were: 1993 (for LFC ≥ 1) and 1241 (for LFC ≥ 2) in IVMr and 1929 (for LFC ≥ 1) and 835 (for LFC ≥ 2) in IVMs. The enriched GO terms upregulated per category identified the intracellular structure, intracellular membrane-bounded organelle and integral component of the cell membrane as some principal cellular components. Meanwhile, efflux transmembrane transporter activity, ABC-type xenobiotic transporter activity and ATPase-coupled transmembrane transporter activity were associated with molecular function. Responses to nematicide activity, pharyngeal pumping and positive regulation of synaptic assembly were classified as biological processes that might be involved in events related to the anthelmintic resistance (AR) and nematode biology. The filtering analysis of both LFC values showed similar genes related to AR. This study deepens our knowledge about the mechanisms behind the processes of H. contortus in order to help in tool production and to facilitate the reduction of AR and promote the development of other control strategies, such as anthelmintic drug targets and vaccines.

Transcriptomics and Proteomics of Haemonchus contortus in Response to Ivermectin Treatment

A major problem faced by the agricultural industry is the resistance of Haemonchus contortus to anthelmintic drugs. For a better understanding of the response of H. contortus to IVM and for the screening of drug-resistance-related genes, we used RNA sequencing and isobaric tags for relative and absolute quantification (iTRAQ) technology to detect the transcriptomic and proteomic changes in H. contortus after ivermectin treatment. An integrated analysis of the two omics showed that the differentially expressed genes and proteins were significantly enriched in the pathways of amino acid degradation, the metabolism of xenobiotics by cytochrome P450, the biosynthesis of amino acids, and the tricarboxylic acid cycle. We found that the upregulated UDP-glycosyltransferases (UGT), glutathione S-transferase (GST), cytochrome P450 (CYP), and p-glycoprotein (Pgp) genes play important roles in drug resistance in H. contortus. Our work will help in the understanding of the transcriptome and proteome changes in H. contortus after IVM and will facilitate the discovery of genes related to drug resistance. This information can be further applied to increase the understanding of the response of IVM in relation to H. contortus.

Occurrence and fate of contaminants of emerging concern and their transformation products after uptake by pak choi (Brassica rapa subsp. chinensis)

Recovery and reuse of nutrients is a major challenge in agriculture. A new process contributing to a circular economy is the anaerobic digestion of food waste, which is a sustainable way of recycling nutrients as the digestate can be used as fertiliser in agriculture and horticulture. However, the digestate may be polluted with contaminants of emerging concern (CECs) that can be circulated back into the food chain, posing a risk to the environment and human health. In this work, the nutrient solution was spiked with 18 selected CECs frequently detected in food waste biogas facilities, and subsequent uptake and fate of these CECs were evaluated in pak choi grown in two different nutrient solutions (mineral and organic). All spiked compounds except two (propylparaben, fenbendazole) were taken up by pak choi plants, with perfluorobutanoic acid (PFBA) and sertraline displaying the highest concentrations (270 and 190 μg/kg fresh weight, respectively). There were no statistically significant differences in uptake between mineral and organic nutrient solutions. Uptake of per- and polyfluoroalkyl substances (PFAS) was negatively correlated with perfluorocarbon chain length and dependent on the functional group (r = -0.73). Sixteen transformation products (TPs) were tentatively identified using suspect screening, most of which were Phase II or even Phase III metabolites. Six of these TPs were identified for the first time in plant metabolism and their metabolic pathways were considered.

Molecular Defense Response of Bursaphelenchus xylophilus to the Nematophagous Fungus Arthrobotrys robusta

Bursaphelenchus xylophilus causes pine wilt disease, which poses a serious threat to forestry ecology around the world. Microorganisms are environmentally friendly alternatives to the use of chemical nematicides to control B. xylophilus in a sustainable way. In this study, we isolated a nematophagous fungus-Arthrobotrys robusta-from the xylem of diseased Pinus massoniana. The nematophagous activity of A. robusta against the PWNs was observed after just 6 h. We found that B. xylophilus entered the trap of A. robusta at 24 h, and the nervous system and immunological response of B. xylophilus were stimulated by metabolites that A. robusta produced. At 30 h of exposure to A. robusta, B. xylophilus exhibited significant constriction, and we were able to identify xenobiotics. Bursaphelenchus xylophilus activated xenobiotic metabolism, which expelled the xenobiotics from their bodies, by providing energy through lipid metabolism. When PWNs were exposed to A. robusta for 36 h, lysosomal and autophagy-related genes were activated, and the bodies of the nematodes underwent disintegration. Moreover, a gene co-expression pattern network was constructed by WGCNA and Cytoscape. The gene co-expression pattern network suggested that metabolic processes, developmental processes, detoxification, biological regulation, and signaling were influential when the B. xylophilus specimens were exposed to A. robusta. Additionally, bZIP transcription factors, ankyrin, ATPases, innexin, major facilitator, and cytochrome P450 played critical roles in the network. This study proposes a model in which mobility improved whenever B. xylophilus entered the traps of A. robusta. The model will provide a solid foundation with which to understand the molecular and evolutionary mechanisms underlying interactions between nematodes and nematophagous fungi. Taken together, these findings contribute in several ways to our understanding of B. xylophilus exposed to microorganisms and provide a basis for establishing an environmentally friendly prevention and control strategy.

Prenylated Benzophenones from Vismia Guianensis Reduced Nematode Growth and Chemotaxis

Biopesticides are generally considered a safer and more environmentally friendly alternative to conventional pesticides. Plant metabolites display a range of pest specific activity ranging from antimicrobial to larvicidal and nematocidal. We herein describe the evaluation of a Guyanese collection of Vismia guianensis (Clusiaceae) for anthelmintic activity. The bioassay-guided evaluation of the hexane extract yielded the new prenylated benzophenone 8,9-epoxyvismiaphenone F (1). The final structures were elucidated based on spectral analysis and comparison to the known metabolite. To evaluate the anthelmintic efficacy of these compounds, Caenorhabditis elegans were exposed to the compounds via a ring assay model. Post-exposure, the numbers of live C. elegans in the compound (middle), bacteria ring were recorded for 3 d, as well as the total number of live worms for each plate. Compound 1 reduced C. elegans' overall growth and reproduction, suggesting that these prenylated benzophenones may hold some promise as natural pesticides.

A Multiplex Gene Expression Assay for Direct Measurement of RNA Transcripts in Crude Lysates of the Nematode Caenorhabditis elegans Used as a Bioanalytical Tool

Gene expression profiling in Caenorhabditis elegans has been demonstrated to be a potential bioanalytical tool to detect the toxic potency of environmental contaminants. The RNA transcripts of genes responding to toxic exposure can be used as biomarkers for detecting these toxins. For routine application in environmental quality monitoring, an easy-to-use multiplex assay is required to reliably quantify expression levels of these biomarkers. In the present study, a bead-based assay was developed to fingerprint gene expression in C. elegans by quantitating messenger RNAs (mRNAs) of multiple target genes directly from crude nematode lysates, circumventing RNA extraction and purification steps. The assay uses signal amplification rather than target amplification for direct measurement of toxin-induced RNA transcripts. Using a 50-gene panel, the expression changes of four candidate reference genes and 46 target mRNAs for various contaminants and wastewaters were successfully measured, and the expression profiles indicated the type of toxin present. Moreover, the multiplex assay response was in line with previous results obtained with more time-consuming reverse-transcription quantitative polymerase chain reaction and microarray analyses. In addition, the transcriptomic profiles of nematodes exposed to wastewater samples and extracts prepared from tissues of swimming crabs were evaluated. The profiles indicated the presence of organic pollutants. The present study illustrates the successful development of a multiplex fluorescent bead-based approach using nematode C. elegans crude lysates for gene expression profiling of target RNAs. This method can be used to routinely fingerprint the presence of toxic contaminants in environmental samples and to identify the most biologically active fraction of the contaminant mixture in a toxicity identification and evaluation approach. Environ Toxicol Chem 2023;42:130-142. © 2022 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

CYP35 family in Caenorhabditis elegans biological processes: fatty acid synthesis, xenobiotic metabolism, and stress responses

With more than 80 cytochrome P450 (CYP) encoding genes found in the nematode Caenorhabditis elegans (C. elegans), the cyp35 genes are one of the important genes involved in many biological processes such as fatty acid synthesis and storage, xenobiotic stress response, dauer and eggshell formation, and xenobiotic metabolism. The C. elegans CYP35 subfamily consisted of A, B, C, and D, which have the closest homolog to human CYP2 family. C. elegans homologs could answer part of the hunt for human disease genes. This review aims to provide an overview of CYP35 in C. elegans and their human homologs, to explore the roles of CYP35 in various C. elegans biological processes, and how the genes of cyp35 upregulation or downregulation are influenced by biological processes, upon exposure to xenobiotics or changes in diet and environment. The C. elegans CYP35 gene expression could be upregulated by heavy metals, pesticides, anti-parasitic and anti-chemotherapeutic agents, polycyclic aromatic hydrocarbons (PAHs), nanoparticles, drugs, and organic chemical compounds. Among the cyp35 genes, cyp-35A2 is involved in most of the C. elegans biological processes regulation. Further venture of cyp35 genes, the closest homolog of CYP2 which is the largest family of human CYPs, may have the power to locate cyps gene targets, discovery of novel therapeutic strategies, and possibly a successful medical regime to combat obesity, cancers, and cyps gene-related diseases.

Xenobiotic-Metabolizing Enzymes in Trematodes

Trematode infections occur worldwide causing considerable deterioration of human health and placing a substantial financial burden on the livestock industry. The hundreds of millions of people afflicted with trematode infections rely entirely on only two drugs (praziquantel and triclabendazole) for treatment. An understanding of anthelmintic biotransformation pathways in parasites should clarify factors that can modulate therapeutic potency of anthelmintics currently in use and may lead to the discovery of synergistic compounds for combination treatments. Despite the pronounced epidemiological significance of trematodes, there is still no adequate understanding of the functionality of their metabolic systems, including xenobiotic-metabolizing enzymes. The review is focused on the structure and functional significance of the xenobiotic-metabolizing system in trematodes. Knowledge in this field can solve practical problems related to the search for new targets for antiparasitic therapy based on a focused action on certain elements of the parasite's metabolic system. Knowledge of the functionality of this system is required to understand the adaptation of the biochemical processes of parasites residing in the host and mechanisms of drug resistance development, as well as to select a promising molecular target for the discovery and development of new anthelmintic drugs.

The induction and inhibition of UDP-glycosyltransferases in Haemonchus contortus and their role in the metabolism of albendazole

Albendazole (ABZ) is an anthelmintic frequently used to treat haemonchosis, a common parasitosis of ruminants caused by the gastrointestinal nematode Haemonchus contortus. This parasite is able to protect itself against ABZ via the formation of inactive ABZ-glycosides. The present study was designed to deepen the knowledge about the role of UDP-glycosyltransferases (UGTs) in ABZ glycosylation in H. contortus. The induction effect of phenobarbital, a classical inducer of UGTs, as well as ABZ and ABZ-sulphoxide (ABZSO, the main active metabolite of ABZ) on UGTs expression and UGT activity toward ABZ was studied ex vivo in isolated adult nematodes. The effect of three potential UGT inhibitors (5-nitrouracil, 4,6-dihydroxy-5-nitropyrimidine and sulfinpyrazone) on ABZ glycosylation was tested. Pre-incubation of nematodes with ABZ and ABZSO led to increased expression of several UGTs as well as ABZ-glycosides formation in subsequent treatment. Phenobarbital also induced UGTs expression, but did not affect ABZ biotransformation. In the nematode's subcellular fraction, sulfinpyrazone inhibited UGT activity toward ABZ, although no effect of other inhibitors was observed. The inhibitory potential of sulfinpyrazone on the formation of ABZ-glycosides was also proved ex vivo in living nematodes. The obtained results confirmed the role of UGTs in ABZ biotransformation in H. contortus adults and revealed sulfinpyrazone as a potent inhibitor of ABZ glycosylation in this parasite. The possible use of sulfinpyrazone with ABZ in combination therapy merits further research.

•

UDP-glycosyl transferases catalyse biotransformation of benzimidazole anthelmintics.

•

Gene expression of several UGTs in H. contortus is affected by albendazole.

•

The inhibition of UGTs leads to a decrease in the amounts of glycosylated metabolites.

Interactions of C. elegans β-tubulins with the microtubule inhibitor and anthelmintic drug albendazole

Parasitic nematodes are major human and agricultural pests, and benzimidazoles are amongst the most important broad spectrum anthelmintic drug class used for their control. Benzimidazole resistance is now widespread in many species of parasitic nematodes in livestock globally and an emerging concern for the sustainable control of human soil transmitted helminths. β-tubulin is the major benzimidazole target, although other genes may influence resistance. Among the six C. elegans β-tubulin genes, loss of ben-1 causes resistance without other apparent defects. Here, we explored the genetics of C. elegans β-tubulin genes in relation to the response to the benzimidazole derivative albendazole. The most highly expressed β-tubulin isotypes, encoded by tbb-1 and tbb-2, were known to be redundant with each other for viability, and their products are predicted not to bind benzimidazoles. We found that tbb-2 mutants, and to a lesser extent tbb-1 mutants, were hypersensitive to albendazole. The double mutant tbb-2 ben-1 is uncoordinated and short, resembling the wild type exposed to albendazole, but the tbb-1 ben-1 double mutant did not show the same phenotypes. These results suggest that tbb-2 is a modifier of ABZ sensitivity. To better understand how BEN-1 mutates to cause benzimidazole resistance, we isolated mutants resistant to albendazole and found that 15 of 16 mutations occurred in the ben-1 coding region. Mutations ranged from likely nulls to hypomorphs, and several corresponded to residues that cause resistance in other organisms. Null alleles of ben-1 are albendazole-resistant and BEN-1 shows high sequence identity with tubulins from other organisms, suggesting that many amino acid changes could cause resistance. However, our results suggest that missense mutations conferring resistance are not evenly distributed across all possible conserved sites. Independent of their roles in benzimidazole resistance, tbb-1 and tbb-2 may have specialized functions as null mutants of tbb-1 or tbb-2 were cold or heat sensitive, respectively.

Transcriptomic analyses implicate neuronal plasticity and chloride homeostasis in ivermectin resistance and response to treatment in a parasitic nematode

The antiparasitic drug ivermectin plays an essential role in human and animal health globally. However, ivermectin resistance is widespread in veterinary helminths and there are growing concerns of sub-optimal responses to treatment in related helminths of humans. Despite decades of research, the genetic mechanisms underlying ivermectin resistance are poorly understood in parasitic helminths. This reflects significant uncertainty regarding the mode of action of ivermectin in parasitic helminths, and the genetic complexity of these organisms; parasitic helminths have large, rapidly evolving genomes and differences in evolutionary history and genetic background can confound comparisons between resistant and susceptible populations. We undertook a controlled genetic cross of a multi-drug resistant and a susceptible reference isolate of Haemonchus contortus, an economically important gastrointestinal nematode of sheep, and ivermectin-selected the F2 population for comparison with an untreated F2 control. RNA-seq analyses of male and female adults of all populations identified high transcriptomic differentiation between parental isolates, which was significantly reduced in the F2, allowing differences associated specifically with ivermectin resistance to be identified. In all resistant populations, there was constitutive upregulation of a single gene, HCON_00155390:cky-1, a putative pharyngeal-expressed transcription factor, in a narrow locus on chromosome V previously shown to be under ivermectin selection. In addition, we detected sex-specific differences in gene expression between resistant and susceptible populations, including constitutive upregulation of a P-glycoprotein, HCON_00162780:pgp-11, in resistant males only. After ivermectin selection, we identified differential expression of genes with roles in neuronal function and chloride homeostasis, which is consistent with an adaptive response to ivermectin-induced hyperpolarisation of neuromuscular cells. Overall, we show the utility of a genetic cross to identify differences in gene expression that are specific to ivermectin selection and provide a framework to better understand ivermectin resistance and response to treatment in parasitic helminths.

Ivermectin-induced gene expression changes in adult Parascaris univalens and Caenorhabditis elegans: a comparative approach to study anthelminthic metabolism and resistance in vitro

Background

The nematode Parascaris univalens is one of the most prevalent parasitic pathogens infecting horses but anthelmintic resistance undermines treatment approaches. The molecular mechanisms underlying drug activity and resistance remain poorly understood in this parasite since experimental in vitro models are lacking. The aim of this study was to evaluate the use of Caenorhabditis elegans as a model for P. univalens drug metabolism/resistance studies by a comparative gene expression approach after in vitro exposure to the anthelmintic drug ivermectin (IVM).

Methods

Twelve adult P. univalens worms in groups of three were exposed to ivermectin (IVM, 10^–13 M, 10^–11 M, 10^–9 M) or left unexposed for 24 h at 37 °C, and total RNA, extracted from the anterior end of the worms, was sequenced using Illumina NovaSeq. Differentially expressed genes (DEGs) involved in metabolism, transportation, or gene expression with annotated Caernorhabditis elegans orthologues were identified as candidate genes to be involved in IVM metabolism/resistance. Similarly, groups of 300 adult C. elegans worms were exposed to IVM (10^–9 M, 10^–8 M and 10^–7 M) or left unexposed for 4 h at 20 °C. Quantitative RT-PCR of RNA extracted from the C. elegans worm pools was used to compare against the expression of selected P. univalens candidate genes after drug treatment.

Results

After IVM exposure, 1085 DEGs were found in adult P. univalens worms but the relative gene expression changes were small and large variabilities were found between different worms. Fifteen of the DEGs were chosen for further characterization in C. elegans after comparative bioinformatics analyses. Candidate genes, including the putative drug target lgc-37, responded to IVM in P. univalens, but marginal to no responses were observed in C. elegans despite dose-dependent behavioral effects observed in C. elegans after IVM exposure. Thus, the overlap in IVM-induced gene expression in this small set of genes was minor in adult worms of the two nematode species.

Conclusion

This is the first time to our knowledge that a comparative gene expression approach has evaluated C. elegans as a model to understand IVM metabolism/resistance in P. univalens. Genes in P. univalens adults that responded to IVM treatment were identified. However, identifying conserved genes in P. univalens and C. elegans involved in IVM metabolism/resistance by comparing gene expression of candidate genes proved challenging. The approach appears promising but was limited by the number of genes studied (n = 15). Future studies comparing a larger number of genes between the two species may result in identification of additional candidate genes involved in drug metabolism and/or resistance.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05260-4.

Metabolism and pharmacokinetics of mebendazole in Japanese pufferfish (Takifugu rubripes)

As a typical and broad-spectrum benzimidazole, mebendazole (MBZ) has long been used in human and veterinary medicine to treat parasitic infestations, and is widely employed in the aquaculture of Japanese pufferfish (Takifugu rubripes). However, there have been no studies examining the pharmacokinetic characteristics of MBZ in Japanese pufferfish. Furthermore, the presence of MBZ and its metabolites in animal-derived raw food represents a notable safety concern. Here, we investigated the metabolism of MBZ using a UPLC-Q-TOF system. Additionally, we evaluated the pharmacokinetics of MBZ and two metabolites, 2-amino-5(6)-benzoylbenzimidazole (MBZ-NH₂) and 5-hydroxymebendazole (MBZ-OH), in Japanese pufferfish following intramuscular injection of 20 mg/kg MBZ. We detected three metabolites of MBZ (M1-M3), among which, 2-amino-5(6)-(a-hydroxybenzyl) benzimidazole (M3) was detected in an aquatic animal for the first time. The plasma dispositions of MBZ, MBZ-NH₂, and MBZ-OH were characterized by low plasma clearance, medium distribution volume, and long terminal half-life. Moreover, these compounds were widely distributed in the muscle, from which they were rapidly cleared. The pharmacokinetics and metabolism of mebendazole in Japanese pufferfish are described for the first time in this study. Our findings provide a basis for the rational application of MBZ in Japanese pufferfish farming and contribute to our understanding of the metabolism of MBZ in cultured fish.

Absence of Polymorphisms in Codons 167, 198 and 200 of All Seven β-tubulin Isotypes of Benzimidazole Susceptible and Resistant Parascaris spp. Specimens from Australia

Benzimidazoles resistance is widespread in strongyle parasitic nematodes and associated with polym orphisms in the codons 167, 198 and 200 of isotype 1 β-tubulin (tbb-1). In ascarids, benzimidazole (BZ) resistance has rarely been reported and in none of these cases were any of these polymorphisms detected. Here, available genome and transcriptome data from WormBase ParaSite were used to compare the complete β-tubulin reservoirs of Parascaris univalens, Ascaris suum and Ascaris lumbricoides. Adult Parascaris spp. specimens collected in Australia from horses after BZ treatment (susceptible, n = 13) or surviving BZ treatment and collected after ivermectin treatment (resistant, n = 10) were genotyped regarding codons 167, 198 and 200 using Sanger sequencing. Phylogenetic analyses clearly showed that there are no one-to-one ascarid orthologs of strongyle tbb-1 genes. In the reference genomes, as well as phenotypically susceptible and resistant Parascaris spp. from Australia, six out of seven β-tubulin genes showed a BZ-susceptible genotype (F167, E198, F200). The only exception were the testis-specific β-tubulin D genes from all three ascarid species that encode tyrosine at codon 200. This was observed independently of the BZ-susceptibility phenotype of Parascaris spp. These data suggest that different mechanisms lead to BZ resistance in ascarid and strongyle nematodes.

Helminth lipidomics: Technical aspects and future prospects

Lipidomics is a relatively recent molecular research field, and explores lipids (fats) and their biology using advanced mass spectrometry technologies. Although this field has expanded significantly in biomedical and biotechnological disciplines, it is still in its infancy in molecular parasitology. Our goal here is to review and discuss technical aspects of MS-based lipidomics and its recent applications to parasitic worms, as well as challenges and future directions for worm lipid research. In a multi-omic paradigm, we expect that the exploration of lipidomic data for parasitic worms will yield important insights into lipid-associated biological pathways and processes, including the regulation of essential signalling pathways, parasite invasion, establishment, adaptation and development.

•

Lipids are involved in critical biological functions in parasitic worms.

•

Lipidomics is an emerging research field in molecular helminthology.

•

This article covers technological advances and applications to parasitic worms.

•

It also discusses challenges and future directions for lipidomic research.

Screening the pandemic response box identified benzimidazole carbamates, Olorofim and ravuconazole as promising drug candidates for the treatment of eumycetoma

Eumycetoma is a chronic subcutaneous neglected tropical disease that can be caused by more than 40 different fungal causative agents. The most common causative agents produce black grains and belong to the fungal orders Sordariales and Pleosporales. The current antifungal agents used to treat eumycetoma are itraconazole or terbinafine, however, their cure rates are low. To find novel drugs for eumycetoma, we screened 400 diverse drug-like molecules from the Pandemic Response Box against common eumycetoma causative agents as part of the Open Source Mycetoma initiative (MycetOS). 26 compounds were able to inhibit the growth of Madurella mycetomatis, Madurella pseudomycetomatis and Madurella tropicana, 26 compounds inhibited Falciformispora senegalensis and seven inhibited growth of Medicopsis romeroi in vitro. Four compounds were able to inhibit the growth of all five species of fungi tested. They are the benzimidazole carbamates fenbendazole and carbendazim, the 8-aminoquinolone derivative tafenoquine and MMV1578570. Minimal inhibitory concentrations were then determined for the compounds active against M. mycetomatis. Compounds showing potent activity in vitro were further tested in vivo. Fenbendazole, MMV1782387, ravuconazole and olorofim were able to significantly prolong Galleria mellonella larvae survival and are promising candidates to explore in mycetoma treatment and to also serve as scaffolds for medicinal chemistry optimisation in the search for novel antifungals to treat eumycetoma.

Comparison of electrophysiological and motility assays to study anthelmintic effects in Caenorhabditis elegans

Currently, only a few chemical drug classes are available to control the global burden of nematode infections in humans and animals. Most of these drugs exert their anthelmintic activity by interacting with proteins such as ion channels, and the nematode neuromuscular system remains a promising target for novel intervention strategies. Many commonly-used phenotypic readouts such as motility provide only indirect insight into neuromuscular function and the site(s) of action of chemical compounds. Electrophysiological recordings provide more specific information but are typically technically challenging and lack high throughput for drug discovery. Because drug discovery relies strongly on the evaluation and ranking of drug candidates, including closely related chemical derivatives, precise assays and assay combinations are needed for capturing and distinguishing subtle drug effects.

Past studies show that nematode motility and pharyngeal pumping (feeding) are inhibited by most anthelmintic drugs. Here we compare two microfluidic devices (“chips”) that record electrophysiological signals from the nematode pharynx (electropharyngeograms; EPGs) ─ the ScreenChip™ and the 8-channel EPG platform ─ to evaluate their respective utility for anthelmintic research. We additionally compared EPG data with whole-worm motility measurements obtained with the wMicroTracker instrument. As references, we used three macrocyclic lactones (ivermectin, moxidectin, and milbemycin oxime), and levamisole, which act on different ion channels. Drug potencies (IC₅₀ and IC₉₅ values) from concentration-response curves, and the time-course of drug effects, were compared across platforms and across drugs. Drug effects on pump timing and EPG waveforms were also investigated. These experiments confirmed drug-class specific effects of the tested anthelmintics and illustrated the relative strengths and limitations of the different assays for anthelmintic research.

•

Anthelmintic drugs inhibit pharyngeal pumping and motility in C. elegans.

•

Two electrophysiological assays and one motility assay were compared.

•

Macrocyclic lactones and levamisole have drug-class-specific effects.

•

A combination of assays most fully reveals anthelmintic effects.

•

Strengths and limitations of the three assays were identified.

Encapsulating Anti-Parasite Benzimidazole Drugs into Lipid-Coated Calcium Phosphate Nanoparticles to Efficiently Induce Skin Cancer Cell Apoptosis

Benzimidazole (BMZ) family of anti-worm drugs has been now repurposed as anti-cancer drugs. However, offering a general reformulation method for these drugs is essential due to their hydrophobicity and low aqueous solubility. In this work, we developed a general approach to load typical BMZ drugs as tiny nanocrystals within lipid-coated calcium phosphate (LCP) nanoparticles. BMZ drug-loaded LCP nanoparticles increased their solubility in PBS by 100–200% and significantly enhanced the anti-cancer efficacy in the treatment of B16F0 melanoma cells. These drug-LCP nanoparticles induced much more cancer cell apoptosis, generated much more reactive oxygen species (ROS) and inhibited Bcl-2 expression of cancer cells. Moreover, BMZ drug-loaded LCP nanoparticles caused morphological change and extension disruption of cancer cells, and significantly reduced migration activity, representing high possibility for inhibition of tumor dissemination and metastasis. Very advantageously, BMZ drug-loaded LCP nanoparticles did not show any obvious toxicity, Bcl-2 inhibition and morphological changes in HEK293T healthy cells. In conclusion, BMZ drug-incorporated LCP nanoformulations may be a valuable nanomedicine that is able to inhibit primary tumors and prevent tumor dissemination with minimum side effects on healthy cells and tissues.

Two novel loci underlie natural differences in Caenorhabditis elegans abamectin responses

Parasitic nematodes cause a massive worldwide burden on human health along with a loss of livestock and agriculture productivity. Anthelmintics have been widely successful in treating parasitic nematodes. However, resistance is increasing, and little is known about the molecular and genetic causes of resistance for most of these drugs. The free-living roundworm Caenorhabditis elegans provides a tractable model to identify genes that underlie resistance. Unlike parasitic nematodes, C. elegans is easy to maintain in the laboratory, has a complete and well annotated genome, and has many genetic tools. Using a combination of wild isolates and a panel of recombinant inbred lines constructed from crosses of two genetically and phenotypically divergent strains, we identified three genomic regions on chromosome V that underlie natural differences in response to the macrocyclic lactone (ML) abamectin. One locus was identified previously and encodes an alpha subunit of a glutamate-gated chloride channel (glc-1). Here, we validate and narrow two novel loci using near-isogenic lines. Additionally, we generate a list of prioritized candidate genes identified in C. elegans and in the parasite Haemonchus contortus by comparison of ML resistance loci. These genes could represent previously unidentified resistance genes shared across nematode species and should be evaluated in the future. Our work highlights the advantages of using C. elegans as a model to better understand ML resistance in parasitic nematodes.

UDP-Glycosyltransferases and Albendazole Metabolism in the Juvenile Stages of Haemonchus contortus

The nematode Haemonchus contortus, a gastrointestinal parasite of ruminants, can severely burden livestock production. Although anthelmintics are the mainstay in the treatment of haemonchosis, their efficacy diminishes due to drug-resistance development in H. contortus. An increased anthelmintics inactivation via biotransformation belongs to a significant drug-resistance mechanism in H. contortus. UDP-glycosyltransferases (UGTs) participate in the metabolic inactivation of anthelmintics and other xenobiotic substrates through their conjugation with activated sugar, which drives the elimination of the xenobiotics due to enhanced solubility. The UGTs family, in terms of the biotransformation of commonly used anthelmintics, has been well described in adults as a target stage. In contrast, the free-living juvenile stages of H. contortus have attracted less attention. The expression of UGTs considerably varies throughout the life cycle of the juvenile nematodes, suggesting their different roles. Furthermore, the constitutive expression in a susceptible strain with two resistant strains shows several resistance-related changes in UGTs expression, and the exposure of juvenile stages of H. contortus to albendazole (ABZ) and ABZ-sulfoxide (ABZSO; in sublethal concentrations) leads to the increased expression of several UGTs. The anthelmintic drug ABZ and its primary metabolite ABZSO biotransformation, tested in the juvenile stages, shows significant differences between susceptible and resistant strain. Moreover, higher amounts of glycosidated metabolites of ABZ are formed in the resistant strain. Our results show similarly, as in adults, the UGTs and glycosidations significant for resistance-related differences in ABZ biotransformation and warrant further investigation in their individual functions.

Correlation of NHR-48 Transcriptional Modulator Expression with Selected CYP Genes’ Expression during Thiabendazole Treatment of Anisakis simplex (s.l.)?—An In Vitro Study

Anisakis simplex (s.l.) is a complex of three sibling (biological) species of parasitic nematodes of marine mammals, including A. berlandi, A. pegreffii and A. simplex (s.s.). It is characterized by a complex life cycle in which humans can become accidental hosts by consuming dishes made of raw or undercooked fish containing L3 larvae, which in many regions of the world is related to the national or regional culinary tradition. This has spurred scientific efforts to develop new methods for treating the disease, called anisakiasis, and to neutralize invasive L3. Thiabendazole (TBZ) is a wide-spectrum anthelminthic with a higher efficacy than albendazole, a drug whose long-term use induces resistance in many parasitic species. Cytochromes P450 participate in TBZ metabolism, and the expression of their genes is controlled by nuclear hormone receptors (NHR). This study aimed to examine the effects of TBZ on the above-described pathway in invasive larvae of A. simplex (s.l.). The efficacy of TBZ against A. simplex (s.l.) larvae was observed for the first time. Larvae were cultured in vitro for 72 h in a medium containing TBZ at five concentrations from 0.5 to 1.5 mM. However, the survival curves did not significantly differ from each other. This means that all of the concentrations of TBZ had a similar effect on the A. simplex (s.l.) L3 larvae during in vitro culture. Nevertheless, TBZ modified the expression of nhr-48, cyp13a3 and cyp1a1 genes in the L3 of A. simplex (s.l.).

Quantitative benzimidazole resistance and fitness effects of parasitic nematode beta-tubulin alleles

Infections by parasitic nematodes inflict a huge burden on the health of humans and livestock throughout the world. Anthelmintic drugs are the first line of defense against these infections. Unfortunately, resistance to these drugs is rampant and continues to spread. To improve treatment strategies, we must understand the genetics and molecular mechanisms that underlie resistance. Studies of the fungus Aspergillus nidulans and the free-living nematode Caenorhabditis elegans discovered that a beta-tubulin gene is mutated in benzimidazole (BZ) resistant strains. In parasitic nematode populations, three beta-tubulin alleles, F167Y, E198A, and F200Y, have long been correlated with resistance. Additionally, improvements in sequencing technologies have identified new alleles - E198V, E198L, E198K, E198I, and E198Stop - also correlated with BZ resistance. However, none of these alleles have been proven to cause resistance. To empirically demonstrate this point, we independently introduced the F167Y, E198A, and F200Y alleles as well as two of the newly identified alleles, E198V and E198L, into the BZ susceptible C. elegans N2 genetic background using the CRISPR-Cas9 system. These genome-edited strains were exposed to both albendazole and fenbendazole to quantitatively measure animal responses to BZs. We used a range of concentrations for each BZ compound to define response curves and found that all five of the alleles conferred resistance to BZ compounds equal to a loss of the entire beta-tubulin gene. These results prove that the parasite beta-tubulin alleles cause resistance. The E198V allele is found at low frequencies along with the E198L allele in natural parasite populations, suggesting that it could affect fitness. We performed competitive fitness assays and demonstrated that the E198V allele reduces animal health, supporting the hypothesis that this allele might be less fit in field populations. Overall, we present a powerful platform to quantitatively assess anthelmintic resistance and effects of specific resistance alleles on organismal fitness in the presence or absence of the drug.

Rapid determination of nematode cell and organ susceptibility to toxic treatments

In research focused on the intestine of parasitic nematodes, we recently identified small molecule inhibitors toxic to intestinal cells of larval Ascaris suum (nematode intestinal toxins/toxicants; “NITs”). Some NITs had anthelmintic activity across the phylogenetic diversity of the Nematoda. The whole-worm motility inhibition assay quantified anthelmintic activity, but worm responses to NITs in relation to pathology or affected molecular pathways was not acquired. In this study we extended this research to more comprehensively determine in whole larval A. suum the cells, organ systems, molecular targets, and potential cellular pathways involved in mechanisms of toxicity leading to cell death. The experimental system utilized fluorescent nuclear probes (bisbenzimide, propidium iodide), NITs, an A. suum larval parasite culture system and transcriptional responses (RNA-seq) to NITs. The approach provides for rapid resolution of NIT-induced cell death among organ systems (e.g. intestine, excretory, esophagus, hypodermis and seam cells, and nervous), discriminates among NITs based on cell death profiles, and identifies cells and organ systems with the greatest NIT sensitivity (e.g. intestine and apparent neuronal cells adjacent to the nerve ring). Application was extended to identify cells and organs sensitive to several existing anthelmintics. This approach also resolved intestinal cell death and irreparable damage induced in adult A. suum by two NITs, establishing a new model to elucidate relevant pathologic mechanisms in adult worms. RNA-seq analysis resolved A. suum genes responsive to treatments with three NITs, identifying dihydroorotate dehydrogenase (uridine synthesis) and RAB GTPase(s) (vesicle transport) as potential targets/pathways leading to cell death. A set of genes induced by all three NITs tested suggest common stress or survival responses activated by NITs. Beyond the presented specific lines of research, elements of the overall experimental system presented in this study have broad application toward systematic development of new anthelmintics.

•

A unique rapid cell death assay was developed for parasitic nematodes.

•

Multiple drug-like molecules cause widespread cell death in many organs of A. suum.

•

Multiple cell and organ systems were validated as targets for anthelmintics.

•

Potential drug targets/pathways were implicated in activating cell death processes.

Toward integrative 'omics of the barber's pole worm and related parasitic nematodes

Advances in nucleic acid sequencing, mass spectrometry and computational biology have facilitated the identification, annotation and analysis of genes, transcripts, proteins and metabolites in model nematodes (Caenorhabditis elegans and Pristionchus pacificus) and socioeconomically important parasitic nematodes (Clades I, III, IV and V). Significant progress has been made in genomics and transcriptomics as well as in the proteomics and lipidomics of Haemonchus contortus (the barber's pole worm) - one of the most pathogenic representatives of the order Strongylida. Here, we review salient aspects of genomics, transcriptomics, proteomics, lipidomics, glycomics and functional genomics, and discuss the rise of integrative 'omics of this economically important parasite. Although our knowledge of the molecular biology, genetics and biochemistry of H. contortus and related species has progressed significantly, much remains to be explored, particularly in areas such as drug resistance, unique/unknown genes, host-parasite interactions, parasitism and the pathogenesis of disease, by integrating the use of multiple 'omics methods. This approach should lead to a better understanding of H. contortus and its relatives at a 'systems biology' level, and should assist in developing new interventions against these parasites.

Sub-lethal doses of albendazole induce drug metabolizing enzymes and increase albendazole deactivation in Haemonchus contortus adults

The efficacy of anthelmintic therapy of farm animals rapidly decreases due to drug resistance development in helminths. In resistant isolates, the increased expression and activity of drug-metabolizing enzymes (DMEs), e.g. cytochromes P450 (CYPs), UDP-glycosyltransferases (UGTs) and P-glycoprotein transporters (P-gps), in comparison to sensitive isolates have been described. However, the mechanisms and circumstances of DMEs induction are not well known. Therefore, the present study was designed to find the changes in expression of CYPs, UGTs and P-gps in adult parasitic nematodes Haemonchus contortus exposed to sub-lethal doses of the benzimidazole anthelmintic drug albendazole (ABZ) and its active metabolite ABZ-sulfoxide (ABZSO). In addition, the effect of ABZ at sub-lethal doses on the ability to deactivate ABZ during consequent treatment was studied. The results showed that contact of H. contortus adults with sub-lethal doses of ABZ and ABZSO led to a significant induction of several DMEs, particularly cyp-2, cyp-3, cyp-6, cyp-7, cyp-8, UGT10B1, UGT24C1, UGT26A2, UGT365A1, UGT366C1, UGT368B2, UGT367A1, UGT371A1, UGT372A1 and pgp-3, pgp-9.1, pgp-9.2, pgp-10. This induction led to increased formation of ABZ metabolites (especially glycosides) and their increased export from the helminths' body into the medium. The present study demonstrates for the first time that contact of H. contortus with sub-lethal doses of ABZ (e.g. during underdose treatment) improves the ability of H. contortus adults to deactivate ABZ in consequent therapy.

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.