UDP-glycosyltransferase family in Haemonchus contortus: Phylogenetic analysis, constitutive expression, sex-differences and resistance-related differences

Study on the effect of GCY-12 gene on albendazole sensitivity of Haemonchus contortus by RNA interference

Introduction

Haemonchus contortus (H. contortus) is a common gastrointestinal nematode in small ruminants, posing a significant threat to the livestock industry. The control of H. contortus often depends on drugs such as albendazole. However, the prolonged and improper use of these drugs by livestock producers has led to widespread resistance among ruminant populations, posing a major challenge to parasite management. It is reported that the resistance of H. contortus to albendazole is related to the single nucleotide polymorphism (SNP) of β-tubulin homologous type I gene, but whether other genes are involved has not been reported.

Methods

Based on the comparative analysis of the transcriptome sequencing data of albendazole-sensitive and albendazole-resistant strains of H. contortus, HCON_00043720 (Receptor-type guanylate cyclase GCY-12, GCY-12) gene was selected as the research object from the 23 differential genes in the cyclic guanosine monophosphate (cGMP) signaling pathway where the growth and development of H. contortus dauer stage is located. The GCY-12 gene of H. contortus eggs was silenced by RNA interference (RNAi) test. The reaction temperature was optimized, the expression of silenced eggs was detected by quantitative real-time polymerase chain reaction (qRT-PCR) technology, and the drug resistance of silenced eggs was detected by egg hatch assay (EHA).

Results

The results showed that the gene expression level decreased significantly after GCY-12 interference, and the sensitivity of H. contortus to albendazole increased.

Discussion

This study highlights the potential role of GCY-12 in modulating albendazole resistance in H. contortus, offering new insights for developing effective therapeutic approaches.

Environment-like exposure of free-living stages of Haemonchus contortus to albendazole: The intra- and intergeneration stability of altered expression of UDP-glycosyl transferases and ABC-transporters

Anthelmintic drugs are widespread environmental contaminants, but their impact is still poorly understood. Although contact of parasitic nematode Haemonchus contortus with traces of anthelmintic drug albendazole (ABZ) altered the expression and activity of several UDP-glycosyl transferases (UGTs) and P-glycoproteins (Pgps, belonging to ABC-transporters), key enzymes in endogenous and xenobiotic metabolism, it is not known whether these changes will last during the life cycle and pass to the next generations. In the present study simulating the environmental-like exposure, free-living stages of H. contortus were exposed or unexposed to a sub-lethal dose of ABZ and its transformation products (ABZs) during L3 development. The L3 served for lambs' infection and obtaining of H. contortus adults and eggs, which were again exposed or unexposed to ABZs during L3 development. The expression pattern of UGTs and Pgps was analysed and compared in the first generation of L3, in the adults, and in the second generation of L3. The results showed that ABZs exposition during larvae development altered the expression of several ugt and pgp genes in L3 and adults. The intrageneration stability of ABZs-evoked changes was observed in the case of three genes, four genes maintained the intergeneration stability. Interestingly, ABZs-induced changes in the expression of some genes became apparent only in the second generation of L3. Taking together, contact of free-living stages of H. contortus with traces of ABZs in the environment evokes changes in the expression of certain UGTs and Pgps, with some of these changes being intra- and inter-generation stable.

Molecular interaction assays in silico of crotapotin from Crotalus durissus terrificus against the molecular target trypanothione reductase from Leishmania braziliensis

Background

Leishmaniasis is a neglected disease that mainly affects impoverished populations and receives limited attention from governments and research institutions. Current treatments are based on antimonial therapies, which present high toxicity and cause significant side effects, such as cardiotoxicity and hepatotoxicity. This study proposes using crotapotin, isolated from Crotalus durissus terrificus venom, as a potential inhibitor of the enzyme trypanothione reductase from Leishmania braziliensis (LbTR).

Methods

In silico assays were conducted to evaluate the interaction of crotapotin with LbTR using molecular docking and molecular dynamics techniques. Recombinant LbTR was expressed in E. coli, and its enzymatic activity was confirmed. The inhibitory action of crotapotin on LbTR was then tested in enzymatic assays.

Results

The stability of these interactions was confirmed over 200 ns molecular dynamics simulations, with a clustering analysis using the GROMACS method revealing a total of 12 distinct clusters. The five most representative clusters showed low RMSD values, indicating high structural stability of the LbTR-crotapotin complex. In particular, cluster 1, with 3,398 frames and an average RMSD of 0.189 nm from the centroid, suggests a dominant stable conformation of the complex. Additional clusters maintained average RMSD values between 0.173 nm and 0.193 nm, further reinforcing the robustness of the complex under physiological conditions. Recombinant LbTR expression was successful, yielding 4.8 mg/L with high purity, as verified by SDS-PAGE. In the enzymatic assays, crotapotin partially inhibited LbTR activity, with an IC50 of 223.4 μM.

Conclusion

The in silico findings suggest a stable and structured interaction between crotapotin and LbTR, with low structural fluctuation, although the inhibition observed in in vitro assays was moderate. These results indicate the potential of crotapotin as a promising basis for developing specific LbTR inhibitors, contributing to the bioprospecting of new antiparasitic agents.

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.

A systematic review of the molecular mechanisms related to anthelmintic resistance in Haemonchus contortus: A contemporary narrative

Haemonchus contortus is a gastrointestinal parasite that affects ruminants (cattle, sheep, etc.), having a significant welfare impact worldwide. The rise of anthelmintic resistance poses a growing challenge to adequate control, compromising the success of treatments. This study presents a systematic review of the molecular mechanisms involved in the resistance of H. contortus to anthelmintic drugs. Following an extensive literature search (9075 total articles/excluding duplications), 61 articles were examined. From these, benzimidazoles (BZD) and macrocyclic lactone (ML) were the most reported drug classes (17 and 29, respectively). The mutations in the β-tubulin gene were the primary mechanism of BZD resistance. Important comparisons from early reports of resistance mechanisms to ML (published before 2020) mainly based on parasite-population gene expression (e.g., ligand-gated chloride channels, LGCC, and P-glycoproteins, PGP) with more recent genomic and transcriptomic data (e.g., transcription factor, cky-1 gene) are highlighted. Additionally, resistance mechanisms to levamisole (LEV) and monepantel are discussed, showing evidence of polymorphisms in genes related to the nicotinic acetylcholine receptors (nAChR). Considering the available data, it is possible to divide the reports into two technological periods, evidencing that the availability of a chromosome-level genome for H. contortus in association with experiments based on controlled genetic crosses and transcriptome-wide data allowed for the visualization of genes and polymorphisms that were previously indistinguishable from unrelated genetic variation (i.e., genetic noise). Therefore, the study of anthelmintic resistance mechanisms is facing new possibilities, reflecting the large data banks and the speed at which this information is being processed. We suggest that new publications on drug resistance should adopt the approaches and refer to this new era of scientific discoveries. Consistent data interpretation, including artificial intelligence (AI) support, will help us to suggest novel biological mechanisms involved in drug resistance and predict its evolution, allowing a more comprehensible approach toward sustainable parasite control strategies.

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.

Chromosome-contiguous genome for the Haecon-5 strain of Haemonchus contortus reveals marked genetic variability and enables the discovery of essential gene candidates

Millions of livestock animals worldwide are infected with the haematophagous barber's pole worm, Haemonchus contortus, the aetiological agent of haemonchosis. Despite the major significance of this parasite worldwide and its widespread resistance to current treatments, the lack of a high-quality genome for the well-defined strain of this parasite from Australia, called Haecon-5, has constrained research in a number of areas including host-parasite interactions, drug discovery and population genetics. To enable research in these areas, we report here a chromosome-contiguous genome (∼280 Mb) for Haecon-5 with high-quality models for 19,234 protein-coding genes. Comparative genomic analyses show significant genomic similarity (synteny) with a UK strain of H. contortus, called MHco3(ISE).N1 (abbreviated as "ISE"), but we also discover marked differences in genomic structure/gene arrangements, distribution of nucleotide variability (single nucleotide polymorphisms (SNPs) and indels) and orthology between Haecon-5 and ISE. We used the genome and extensive transcriptomic resources for Haecon-5 to predict a subset of essential single-copy genes employing a "cross-species" machine learning (ML) approach using a range of features from nucleotide/protein sequences, protein orthology, subcellular localisation, single-cell RNA-seq and/or histone methylation data available for the model organisms Caenorhabditis elegans and Drosophila melanogaster. From a set of 1,464 conserved single copy genes, transcribed in key life-cycle stages of H. contortus, we identified 232 genes whose homologs have critical functions in C. elegans and/or D. melanogaster, and prioritised 10 of them for further characterisation; nine of the 10 genes likely play roles in neurophysiological processes, germline, hypodermis and/or respiration, and one is an unknown (orphan) gene for which no detailed functional information exists. Future studies of these genes/gene products are warranted to elucidate their roles in parasite biology, host-parasite interplay and/or disease. Clearly, the present Haecon-5 reference genome and associated resources now underpin a broad range of fundamental investigations of H. contortus and could assist in accelerating the discovery of novel intervention targets and drug candidates to combat haemonchosis.

Spirocerca lupi draft genome, vaccine and anthelmintic targets

Spirocerca lupi is a parasitic nematode affecting predominantly domestic dogs. It causes spirocercosis, a disease that is often fatal. The assembled draft genome of S. lupi consists of 13,627 predicted protein-coding genes and is approximately 150 Mb in length. Several known anthelmintic gene targets such as for β-Tubulin, glutamate, and GABA receptors as well as known vaccine gene targets such as cysteine protease inhibitor and cytokines were identified in S. lupi by comparing orthologs of C. elegans anthelmintic gene targets as well as orthologs to known vaccine candidates. New anthelmintic targets were predicted through an inclusion-exclusion strategy and new vaccine targets were predicted through an immunoinformatics approach. New anthelminthic targets include DNA-directed RNA polymerases, chitin synthase, polymerases, and other enzymes. New vaccine targets include cuticle collagens. These gene targets provide a starting platform for new drug identification and vaccine design.

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.

Analysis of lncRNA-related studies of ivermectin-sensitive and -resistant strains of Haemonchus contortus

In this study, 858 novel long non-coding RNAs (lncRNAs) were predicted as sensitive and resistant strains of Haemonchus contortus to ivermectin. These lncRNAs underwent bioinformatic analysis. In total, 205 lncRNAs significantly differed using log2 (difference multiplicity) > 1 or log2 (difference multiplicity) <  - 1 and FDR < 0.05 as the threshold for significant difference analysis. We selected five lncRNAs based on significant differences in expression, cis-regulation, and their association with the Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathways. These expressions of lncRNAs, namely MSTRG.12610.1, MSTRG.8169.1, MSTRG.6355.1, MSTRG.980.1, and MSTRG.9045.1, were significantly downregulated. These findings were consistent with the results of transcriptomic sequencing. We further investigated the relative expression of target gene mRNAs and the regulation of mRNA and miRNA, starting with lncRNA cis-regulation of mRNA, and constructed a lncRNA-mRNA-miRNA network regulation. After a series of statistical analyses, we finally screened out UGT8, Unc-116, Fer-related kinase-1, GGPP synthase 1, and sart3, which may be involved in developing drug resistance under the regulation of their corresponding lncRNAs. The findings of this study provide a novel direction for future studies on drug resistance targets.

Transcriptome reveals the roles and potential mechanisms of lncRNAs in the regulation of albendazole resistance in Haemonchus contortus

BACKGROUND

Haemonchus contortus (H. contortus) is the most common parasitic nematode in ruminants and is prevalent worldwide. H. contortus resistance to albendazole (ABZ) hinders the efficacy of anthelmintic drugs, but little is known about the molecular mechanisms that regulate this of drug resistance. Recent research has demonstrated that long noncoding RNAs (lncRNAs) can exert significant influence as pivotal regulators of the emergence of drug resistance.

RESULTS

In this study, transcriptome sequencing was conducted on both albendazole-sensitive (ABZ-sensitive) and albendazole-resistant (ABZ-resistant) H. contortus strains, with three biological replicates for each group. The analysis of lncRNA in the transcriptomic data revealed that there were 276 differentially expressed lncRNA (DElncRNA) between strains with ABZ-sensitive and ABZ-resistant according to the criteria of |log2Foldchange|≥ 1 and FDR < 0.05. Notably, MSTRG.12969.2 and MSTRG.9827.1 exhibited the most significant upregulation and downregulation, respectively, in the resistant strains. The potential roles of the DElncRNAs included catalytic activity, stimulus response, regulation of drug metabolism, and modulation of the immune response. Moreover, we investigated the interactions between DElncRNAs and other RNAs, specifically MSTRG.12741.1, MSTRG.11848.1, MSTRG.5895.1, and MSTRG.14070.1, involved in regulating drug stimulation through cis/trans/antisense/lncRNA‒miRNA-mRNA interaction networks. This regulation leads to a decrease (or increase) in the expression of relevant genes, consequently enhancing the resistance of H. contortus to albendazole. Furthermore, through comprehensive analysis of competitive endogenous RNAs (ceRNAs) involved in drug resistance-related pathways, such as the mTOR signalling pathway and ABC transporter signalling pathway, the relevance of the MSTRG.2499.1-novel-m0062-3p-HCON_00099610 interaction was identified to mainly involve the regulation of catalytic activity, metabolism, ubiquitination and transcriptional regulation of gene promoters. Additionally, quantitative real-time polymerase chain reaction (qRT-PCR) validation indicated that the transcription profiles of six DElncRNAs and six DEmRNAs were consistent with those obtained by RNA-seq.

CONCLUSIONS

The results of the present study allowed us to better understand the changes in the lncRNA expression profile of ABZ-resistant H. contortus. In total, these results suggest that the lncRNAs MSTRG.963.1, MSTRG.12741.1, MSTRG.11848.1 and MSTRG.2499.1 play important roles in the development of ABZ resistance and can serve as promising biomarkers for further study.

Short-chain dehydrogenases in Haemonchus contortus: changes during life cycle and in relation to drug-resistance

Short-chain dehydrogenases/reductases (SDRs) regulate the activities of many hormones and other signaling molecules and participate in the deactivation of various carbonyl-bearing xenobiotics. Nevertheless, knowledge about these important enzymes in helminths remains limited. The aim of our study was to characterize the SDR superfamily in the parasitic nematode Haemonchus contortus. Genome localization of SDRs was explored, and phylogenetic analysis in comparison with SDRs from free-living nematode Caenorhabditis elegans and the domestic sheep (Ovis aries, a typical host of H. contortus) was constructed. The expression profile of selected SDRs during the life cycle along with differences between the drug-susceptible and drug-resistant strains, were also studied. Genome sequencing enabled the identification of 46 members of the SDR family in H. contortus. A number of genes have no orthologue in the sheep genome. In all developmental stages of H. contortus, SDR1, SDR3, SDR5, SDR6, SDR14, and SDR18 genes were the most expressed, although in individual stages, huge differences in expression levels were observed. A comparison of SDRs expression between the drug-susceptible and drug-resistant strains of H. contortus revealed several SDRs with changed expression in the resistant strain. Specifically, SDR1, SDR12, SDR13, SDR16 are SDR candidates related to drug-resistance, as the expression of these SDRs is consistently increased in most stages of the drug-resistant H. contortus. These findings revealing several SDR enzymes of H. contortus warrant further investigation.

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.

Identification and characterization of UDP-glycosyltransferase genes and the potential role in response to insecticides exposure in Bactrocera dorsalis

BACKGROUND

The oriental fruit fly, Bactrocera dorsalis (Hendel) is a worldwide pest damaging a wide range of hosts. Due to the long-term indiscriminate use of insecticides, B. dorsalis has developed serious resistance to several insecticides. UDP-glycosyltransferases (UGTs) are secondary metabolic enzymes involved in biotransformation and play an important role in the metabolism of plant secondary metabolites and synthetic insecticides in insects. Thus, we suspect that UGTs in B. dorsalis play an important role in insecticide tolerance.

RESULTS

In this study, 31 UGT genes were identified in the genome of B. dorsalis, belonging to 13 subfamilies. Real-time quantitative polymerase chain reaction (RT-qPCR) results revealed that 12 UGT genes were highly expressed in the antennae, midgut, Malpighian tubule and fat body. The mRNA expressions of 17 UGT genes were up-regulated upon exposure to λ-cyhalothrin, imidacloprid, abamectin and chlorpyrifos. Knockdown of the selected five UGT genes (BdUGT301D2, BdUGT35F2, BdUGT36K2, BdUGT49D2, BdUGT50B5) by RNA interference increased the mortality of B. dorsalis from 9.29% to 27.22% upon exposure to four insecticides.

CONCLUSION

The abundance of UGTs in B. dorsalis is similar to other insect species, and 12 out of 31 UGTs were specifically expressed in metabolic tissues, suggesting a key role in detoxification. Down-regulation of five selected UGT genes increased the susceptibility of B. dorsalis to various insecticides, indicating that UGTs may play an important role in tolerance of B. dorsalis to multiple insecticides. © 2022 Society of Chemical Industry.

Genomic signatures of selection associated with benzimidazole drug treatments in Haemonchus contortus field populations

Genome-wide methods offer a powerful approach to detect signatures of drug selection. However, limited availability of suitable reference genomes and the difficulty of obtaining field populations with well-defined, distinct drug treatment histories mean there is little information on the signatures of selection in parasitic nematodes and on how best to detect them. This study addresses these knowledge gaps by using field populations of Haemonchus contortus with well-defined benzimidazole treatment histories, leveraging a recently completed chromosomal-scale reference genome assembly. We generated a panel of 49,393 genomic markers to genotype 20 individual adult worms from each of four H. contortus populations: two from closed sheep flocks with an approximate 20 year history of frequent benzimidazole treatment, and two populations with a history of little or no treatment. Sampling occurred in the same geographical region to limit genetic differentiation and maximise the detection sensitivity. A clear signature of selection was detected on chromosome I, centred on the isotype-1 β-tubulin gene. Two additional, but weaker, signatures of selection were detected; one near the middle of chromosome I spanning 3.75 Mbp and 259 annotated genes, and one on chromosome II spanning a region of 3.3 Mbp and 206 annotated genes, including the isotype-2 β-tubulin locus. We also assessed how sensitivity was impacted by sequencing depth, worm number, and pooled versus individual worm sequence data. This study provides the first known direct genome-wide evidence for any parasitic nematode, that the isotype-1 β-tubulin gene is quantitatively the single most important benzimidazole resistance locus. It also identified two additional genomic regions that likely contain benzimidazole resistance loci of secondary importance. This study provides an experimental framework to maximise the power of genome-wide approaches to detect signatures of selection driven by anthelmintic drug treatments in field populations of parasitic nematodes.

Transgenic Expression of Haemonchus contortus Cytochrome P450 Hco-cyp-13A11 Decreases Susceptibility to Particular but Not All Macrocyclic Lactones in the Model Organism Caenorhabditis elegans

The number of reported macrocyclic lactones (ML) resistance cases across all livestock hosts is steadily increasing. Different studies in the parasitic nematode Haemonchus contortus assume the participation of cytochrome P450s (Cyps) enzymes in ML resistance. Still, functional data about their individual contribution to resistance or substrate specificity is missing. Via microinjection, transgenic Caenorhabditis elegans expressing HCON_00141052 (transgene-Hco-cyp-13A11) from extrachromosomal arrays were generated. After 24 h of exposure to different concentrations of ivermectin (IVM), ivermectin aglycone (IVMa), selamectin (SEL), doramectin (DRM), eprinomectin (EPR), and moxidectin (MOX), motility assays were performed to determine the impact of the H. contortus Cyp to the susceptibility of the worms against each ML. While transgene-Hco-cyp-13A11 significantly decreased susceptibility to IVM (four-fold), IVMa (2-fold), and SEL (3-fold), a slight effect for DRM and no effect for MOX, and EPR was observed. This substrate specificity of Hco-cyp-13A11 could not be explained by molecular modeling and docking studies. Hco-Cyp-13A11 molecular models were obtained for alleles from isolates with different resistance statuses. Although 14 amino acid polymorphisms were detected, none was resistance specific. In conclusion, Hco-cyp-13A11 decreased IVM, IVMa, and SEL susceptibility to a different extent, but its potential impact on ML resistance is not driven by polymorphisms.

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.

UGT440A1 Is Associated With Motility, Reproduction, and Pathogenicity of the Plant-Parasitic Nematode Bursaphelenchus xylophilus

Pine wilt disease (PWD) caused by Bursaphelenchus xylophilus is considered a major threat to pine forests worldwide. Uridine diphosphate (UDP)-glycosyltransferases (UGTs) catalyze the conjugation of small lipophilic compounds with sugars and play crucial roles in the detoxification and homeostatic processes in all living organisms. We investigated the molecular characteristics and biological functions of the gene UGT440A1 that encodes UGTs in B. xylophilus. The in situ hybridization results indicated that UGT440A1 is expressed in all developmental stages of B. xylophilus, particularly in the head, intestine, and hypodermis of the second-stage of juveniles (J2), third-stage of juveniles (J3) and fourth-stage of juveniles (J4) females and in almost the whole body of J4 males and adults. Recombinant UGT440A1 was observed mainly in the inclusion bodies, and the enzyme activity assay revealed that UGT440A1 could catalyze the glycosylation reaction of two types of flavonols (kaempferol and quercetin). RNA interference (RNAi) of UGT440A1 suppressed motility, feeding, and reproduction of B. xylophilus. Furthermore, UGT440A1 knockdown caused a delay in the development of PWD symptoms in the pine seedlings inoculated with the nematodes. These results suggest that UGT440A1 is involved in the pathogenic process of B. xylophilus and the information may facilitate a better understanding of the molecular mechanism of PWD.

Environmental circulation of the anthelmintic drug albendazole affects expression and activity of resistance-related genes in the parasitic nematode Haemonchus contortus

Veterinary anthelmintics excreted from treated animals pass to soil, subsequently to plants and then to their consumers. This circulation might have various consequences, including drug-resistance promotion in helminths. The present study was designed to follow the effect of the environmental circulation of the common anthelmintic drug albendazole (ABZ) in real farm conditions on the parasitic nematode Haemonchus contortus in vivo. Two fields with fodder plants (clover and alfalfa) were fertilized, the first with dung from ABZ-treated sheep (at the recommended dosage), the second with dung from non-treated sheep (controls). After a 10-week growth period, the fresh fodder from both fields was used to feed two groups of sheep, which were infected with H. contortus. Eggs and adult nematodes from the animals of both groups were isolated, and various parameters were compared. No significant changes in the eggs' sensitivity to ABZ and thiabendazole were observed. However, significantly increased expression of several cytochromes P450 and UDP-glycosyl transferases as well as increased oxidation and glycosylation of ABZ and ABZ-sulfoxide (ABZ-SO) was found in the exposed nematodes. These results show that ABZ environmental circulation improves the ability of the helminths to deactivate ABZ.

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.

UGT1A Gene Family Members Serve as Potential Targets and Prognostic Biomarkers for Pancreatic Cancer

Background

Pancreatic cancer (PC) is one of the most common cancers worldwide, with high mortality. The UGT1A gene family plays important roles in pharmacology and toxicology, contributing to interindividual differences in drug disposition. However, mRNA expression and prognostic value of the UGT1A gene family in PC have not been identified.

Methods

Oncomine, GEPIA2, DAVID 6.8, Metascape, Kaplan-Meier plotter, cBioPortal, GeneMANIA, TRRUST v2, TIMER, and R software were used in our study.

Results

The transcriptional levels of UGT1A1/3/6/8/9/10 in PC tissues were significantly higher than those in normal tissues. These results were further validated using five pairs of PC tumor tissues and adjacent nontumor tissues. A significant correlation was found between the expression of UGT1A1/6/10 and the pathological stage of PC. PC patients with lower transcriptional levels of UGT1A1/4/5/6/10 were associated with a better prognosis. The differentially expressed UGT1A gene family functions were primarily related to the glucuronidation pathway, cytokine-cytokine receptor interactions, and the ILK signaling pathway. Our data suggest that HNF1A, AHR, and CDX2 are key transcription factors for the UGT1A gene family. Furthermore, the expression levels of UGT1A1/3/8/9/10 were positively correlated with the activities of tumor-infiltrating immune cells, especially B cells. The expression levels of UGT1A6/9 were negatively correlated with macrophage infiltration levels.

Conclusions

These results suggest that the UGT1A gene family could serve as a potential prognostic biomarker and target for PC. However, future studies are required to validate our findings and promote the clinical utility of the UGT1A gene family in PC.

Molecular and evolutionary basis for survival, its failure, and virulence factors of the zoonotic nematode Anisakis pegreffii

Parasitism is a highly successful life strategy and a driving force in genetic diversity that has evolved many times over. Accidental infections of non-targeted hosts represent an opportunity for lateral host switches and parasite niche expansion. However, if directed toward organisms that are phylogenetically distant from parasite's natural host, such as humans, it may present a dead-end environment where the parasite fails to mature or is even killed by host immunity. One example are nematodes of Anisakidae family, genus Anisakis, that through evolution have lost the ability to propagate in terrestrial hosts, but can survive for a limited time in humans causing anisakiasis. To scrutinize versatility of Anisakis to infect an evolutionary-distant host, we performed transcriptomic profiling of larvae successfully migrating through the rat, a representative model of accidental human infection and compared it to that of larvae infecting an evolutionary-familiar, paratenic host (fish). In a homeothermic accidental host Anisakis upregulated ribosome-related genes, cell division, cuticle constituents, oxidative phosphorylation, in an unsuccessful attempt to molt to the next stage. In contrast, in the paratenic poikilothermic host where metabolic pathways were moderately upregulated or silenced, larvae prepared for dormancy by triggering autophagy and longevity pathways. Identified differences and the modelling of handful of shared transcripts, provide the first insights into evolution of larval nematode virulence, warranting their further investigation as potential drug therapy targets.

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.

Caenorhabditis elegans in anthelmintic research - Old model, new perspectives

For more than four decades, the free-living nematode Caenorhabditis elegans has been extensively used in anthelmintic research. Classic genetic screens and heterologous expression in the C. elegans model enormously contributed to the identification and characterization of molecular targets of all major anthelmintic drug classes. Although these findings provided substantial insights into common anthelmintic mechanisms, a breakthrough in the treatment and control of parasitic nematodes is still not in sight. Instead, we are facing increasing evidence that the enormous diversity within the phylum Nematoda cannot be recapitulated by any single free-living or parasitic species and the development of novel broad-spectrum anthelmintics is not be a simple goal. In the present review, we summarize certain milestones and challenges of the C. elegans model with focus on drug target identification, anthelmintic drug discovery and identification of resistance mechanisms. Furthermore, we present new perspectives and strategies on how current progress in C. elegans research will support future anthelmintic research.

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.

A putative UDP-glycosyltransferase from Heterorhabditis bacteriophora suppresses antimicrobial peptide gene expression and factors related to ecdysone signaling

Insect pathogens have adopted an array of mechanisms to subvert the immune pathways of their respective hosts. Suppression may occur directly at the level of host-pathogen interactions, for instance phagocytic capacity or phenoloxidase activation, or at the upstream signaling pathways that regulate these immune effectors. Insect pathogens of the family Baculoviridae, for example, are known to produce a UDP-glycosyltransferase (UGT) that negatively regulates ecdysone signaling. Normally, ecdysone positively regulates both molting and antimicrobial peptide production, so the inactivation of ecdysone by glycosylation results in a failure of host larvae to molt, and probably a reduced antimicrobial response. Here, we examine a putative ecdysteroid glycosyltransferase, Hba_07292 (Hb-ugt-1), which was previously identified in the hemolymph-activated transcriptome of the entomopathogenic nematode Heterorhabditis bacteriophora. Injection of recombinant Hb-ugt-1 (rHb-ugt-1) into Drosophila melanogaster flies resulted in diminished upregulation of antimicrobial peptides associated with both the Toll and Immune deficiency pathways. Ecdysone was implicated in this suppression by a reduction in Broad Complex expression and reduced pupation rates in r Hb-ugt-1-injected larvae. In addition to the finding that H. bacteriophora excreted-secreted products contain glycosyltransferase activity, these results demonstrate that Hb-ugt-1 is an immunosuppressive factor and that its activity likely involves the inactivation of ecdysone.

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.

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

We have undertaken a detailed analysis of the biotransformation of five of the most therapeutically important benzimidazole anthelmintics - albendazole (ABZ), mebendazole (MBZ), thiabendazole (TBZ), oxfendazole (OxBZ) and fenbendazole (FBZ) - in Caenorhabditis elegans and the ruminant parasite Haemonchus contortus. Drug metabolites were detected by LC-MS/MS analysis in supernatants of C. elegans cultures with a hexose conjugate, most likely glucose, dominating for all five drugs. This work adds to a growing body of evidence that glucose conjugation is a major pathway of xenobiotic metabolism in nematodes and may be a target for enhancement of anthelmintic potency. Consistent with this, we found that biotransformation of albendazole by C. elegans reduced drug potency. Glucose metabolite production by C. elegans was reduced in the presence of the pharmacological inhibitor chrysin suggesting that UDP-glucuronosyl/glucosyl transferase (UGT) enzymes may catalyze benzimidazole glucosidation. Similar glucoside metabolites were detected following ex vivo culture of adult Haemonchus contortus. As a step towards identifying nematode enzymes potentially responsible for benzimidazole biotransformation, we characterised the transcriptomic response to each of the benzimidazole drugs using the C. elegans resistant strain CB3474 ben-1(e1880)III. In the case of albendazole, mebendazole, thiabendazole, and oxfendazole the shared transcriptomic response was dominated by the up-regulation of classical xenobiotic response genes including a shared group of UGT enzymes (ugt-14/25/33/34/37/41/8/9). In the case of fenbendazole, a much greater number of genes were up-regulated, as well as developmental and brood size effects suggesting the presence of secondary drug targets in addition to BEN-1. The transcriptional xenobiotic response of a multiply resistant H. contortus strain UGA/2004 was essentially undetectable in the adult stage but present in the L3 infective stage, albeit more muted than C. elegans. This suggests that xenobiotic responses may be less efficient in stages of parasitic nematodes that reside in the host compared with the free-living stages.

•

C. e. & H. c. display hexose conjugation (likely glucose) and excretion of 5 BZs.

•

C. elegans (C.e.) biotransformation of ABZ reduces drug potency.

•

UGT inhibitor chrysin reduces ABZ biotransformation by C. elegans.

•

Transcriptomic response of C. e. (ben-1) to 5 BZs dominated by xenobiotic response and additional targets for FBZ.

•

Minimal transcriptomic response of H. contortus to ABZ exposure.