Characterization of the xenobiotic response of Caenorhabditis elegans to the anthelmintic drug albendazole and the identification of novel drug glucoside metabolites

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.

Identification of cytochrome P450 gene family and functional analysis of HgCYP33E1 from Heterodera glycines

The cytochrome P450 (CYP) genes of nematode play a crucial role in the metabolic detoxification of xenobiotics including pesticides. Heterodera glycines, also known as the soybean cyst nematode, is a sedentary endoparasite that infests plant roots, causing high annual economic losses in soybean production regions globally. In this study, we identified 36 CYP genes at a genome-wide level of the H. glycines isolate TN10 using all CYPs from Caenorhabditis elegans as queries. Subsequently, a full-length cDNA of HgCYP33E1 which was significantly up-regulated by the conventional nematicide abamectin was initially cloned from H. glycines. It presented significantly higher expressions in the second-stage juvenile (J2) compared to other parasitic stages of H. glycines. qRT-PCR analysis suggested that the expression of HgCYP33E1 was also xenobiotically induced by soybean root exudate and the metabolites of biocontrol agents. Using RNA interference (RNAi), we investigated the function of HgCYP33E1 in H. glycines parasitism and nematicide selectivity. Compared to the control and dsGFP-treated group, silencing of HgCYP33E1 did not affect the J2 behaviors and the early invasion ability, while it decreased the number of J4s in soybean roots after 18-d inoculation with the dsHgCYP33E1-treated nematodes. In addition, knockdown of HgCYP33E1 in H. glycines resulted in an increase in J2 mortality after 24-h incubation with abamectin compared to the GFP dsRNA-soaked and the control group. These findings revealed the potential role of HgCYP33E1 in the xenobiotic detoxification pathway of H. glycines. Moreover, our data also provided valuable gene information for studying the functions of the CYP family in H. glycines host adaption.

Computational analysis of variation in C. elegans ugts

Caenorhabditis elegans are free-living nematodes with a relatively short life cycle and a wealth of genomic information across multiple databases. Uridine diphosphate-glycosyltransferases (UGTs) are a family of enzymes involved in Phase II modification of xenobiotics in C. elegans , which is the addition of a sizeable water-soluble molecule to a xenobiotic to allow for its excretion out of a cell. Little is known about the variation in UGTs across wild isolates and how that might affect their innate immune response. We analyzed the diversity in ugt genes across C. elegans isolates from different geographical locations from the Caenorhabditis elegans Natural Diversity Resource (CaeNDR) database. This was accomplished using whole genome data and data identifying genome regions as hyper-divergent for each isotype. We implemented three steps to identify ugt genes and make inferences based on their variation. First, we created a catalog of UGTs in the N2 reference strain and used them to create a phylogenetic tree that depicts the relationships between the UGT protein sequences. We then quantified ugt variation using the strains from the CaeNDR database and used their data to remove hyper-divergent ugt genes. The third step was to catalog the occurrence of minor allele frequency (MAF) > 0.05 for all the ugts to compare how that aligned with genes classified as hyper-divergent by CaeNDR. Of the 67 ugt genes analyzed, 18 were hyper-divergent. This research will help improve our understanding of ugt variation in C. elegans .

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.

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.

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.

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.

Parasitic helminths (worms) infect people and animals throughout the world and are largely controlled with mass administration of anthelmintic drugs. There are a very limited number of anthelmintics available and parasitic helminths can rapidly develop resistance to these drugs. Ivermectin is a widely used anthelmintic in both humans and animals, but resistance is now widespread in the veterinary field. We crossed ivermectin resistant and ivermectin susceptible parasitic helminths and treated them with ivermectin or left them as untreated controls. This provided resistant and susceptible populations with a similar genetic background with which to study differences in gene expression associated with ivermectin resistance. We identified upregulation of a gene with no previous association with drug resistance (HCON_00155390:cky-1) in male and female worms in all resistant populations. This gene is thought to be expressed in the helminth pharynx (mouthpart) and, in mammals, plays a role in controlling nerve function and protecting nerves from damage. This is consistent with the known effects of ivermectin in inhibiting helminth feeding through pharyngeal paralysis and may implicate a novel mechanism that allows resistant worms to survive treatment.

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.

Effects of Echinococcus multilocularis metacestodes infection and drug treatment on the activities of biotransformation enzymes in mouse liver

The changes of biotransformation enzymes will substantially affect the host's ability to metabolize drugs and other xenobiotic compounds. In order to further elucidate this process and promote the development in treatment of echinococcosis, we investigated the effects of Echinococcus multilocularis infection and drug treatment on biotransformation enzymes in mouse liver. In microsomal and cytosolic fractions, from the six activities assayed, significant decrease of glutathione S-transferases (GST) activity and significant increase of 7-pentoxyresorufin (PROD) and NADPH-cytochrome P450 reductase (CPR) activity were observed in the mice. infected with E. multilocularis metacestodes. In addition, after six weeks treatment of albendazole in E. multilocularis infected mice, significant decreased GST activity and significant increase of 7- ethoxyresorufin (EROD), PROD, and particularly 3-fold higher 7-methoxyresorufin (MROD) activity were observed. The 3-bromopyruvate treated mice only exhibited significantly lower GST activity. Our results demonstrate that E. multilocularis metacestodes infection can affect the activities of main hepatic biotransformation enzymes and such alterations of activity may further affect the hepatic biotransformation of xenobiotics. Moreover, albendazole and 3-bromopyruvate, the promising potential drug against Echinococcus, affected different hepatic biotransformation enzymes and may affect their metabolism. The findings will help to develop rational treatments with less side effects and promote the development of more efficient treatments against E. multilocularis.

Macrostomum lignano as a model to study the genetics and genomics of parasitic flatworms

Hundreds of millions of people worldwide are infected by various species of parasitic flatworms. Without treatment, acute and chronical infections frequently lead to the development of severe pathologies and even death. Emerging data on a decreasing efficiency of some important anthelmintic compounds and the emergence of resistance to them force the search for alternative drugs. Parasitic flatworms have complex life cycles, are laborious and expensive in culturing, and have a range of anatomic and physiological adaptations that complicate the application of standard molecular-biological methods. On the other hand, free-living flatworm species, evolutionarily close to parasitic flatworms, do not have the abovementioned difficulties, which makes them potential alternative models to search for and study homologous genes. In this review, we describe the use of the basal free-living flatworm Macrostomum lignano as such a model. M. lignano has a number of convenient biological and experimental properties, such as fast reproduction, easy and non-expensive laboratory culturing, optical body transparency, obligatory sexual reproduction, annotated genome and transcriptome assemblies, and the availability of modern molecular methods, including transgenesis, gene knockdown by RNA interference, and in situ hybridization. All this makes M. lignano amenable to the most modern approaches of forward and reverse genetics, such as transposon insertional mutagenesis and methods of targeted genome editing by the CRISPR/Cas9 system. Due to the availability of an increasing number of genome and transcriptome assemblies of different parasitic flatworm species, new knowledge generated by studying M. lignano can be easily translated to parasitic flatworms with the help of modern bioinformatic methods of comparative genomics and transcriptomics. In support of this, we provide the results of our bioinformatics search and analysis of genes homologous between M. lignano and parasitic flatworms, which predicts a list of promising gene targets for subsequent research.

Cytochrome P450 monooxygenase of Acanthamoeba castellanii participates in resistance to polyhexamethylene biguanide treatment

Acanthamoeba spp. are free-living parasites that can cause severe infections such as granulomatous amoebic encephalitis (GAE) and amoebic keratitis (AK). Polyhexamethylene biguanide (PHMB) is a topical application for AK treatment. However, PHMB is not entirely effective against all Acanthamoeba strains or isolates. The mechanisms by which Acanthamoeba protects itself against extreme drug conditions without encystation are still unknown. According to a previous study, cytochrome P450 monooxygenase (CYP450MO) plays an important role in the oxidative biotransformation of numerous drugs related to metabolism. In this study, a CYP450MO fragment was inserted into the pGAPDH-EGFP vector and transfected into Acanthamoeba castellanii. We found that CYP450MO-overexpressing Acanthamoeba had higher survival rates than those of the control cells after PHMB treatment. Moreover, we also found that encystation-related genes such as cellulose synthase I (CSI), encystation-mediating serine proteinase (EMSP), and autophagy-related protein 8 (ATG8) expression levels were not significantly different between Acanthamoeba transfected by pGAPDH-EGFP or pGAPDH-EGFP-CYP450MO. We suggest that Acanthamoeba transfected by pGAPDH-EGFP-CYP450MO may not induce encystation-related genes to resist PHMB treatment. In conclusion, these findings indicate that CYP450MO may be an additional target when PHMB is used for treatment of amoebic keratitis.

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

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

•

First case of fenbendazole resistance in Parascaris univalens in Europe.

•

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

•

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

•

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

Xenobiotic metabolism and transport in Caenorhabditis elegans

Caenorhabditis elegans has emerged as a major model in biomedical and environmental toxicology. Numerous papers on toxicology and pharmacology in C. elegans have been published, and this species has now been adopted by investigators in academic toxicology, pharmacology, and drug discovery labs. C. elegans has also attracted the interest of governmental regulatory agencies charged with evaluating the safety of chemicals. However, a major, fundamental aspect of toxicological science remains underdeveloped in C. elegans: xenobiotic metabolism and transport processes that are critical to understanding toxicokinetics and toxicodynamics, and extrapolation to other species. The aim of this review was to initially briefly describe the history and trajectory of the use of C. elegans in toxicological and pharmacological studies. Subsequently, physical barriers to chemical uptake and the role of the worm microbiome in xenobiotic transformation were described. Then a review of what is and is not known regarding the classic Phase I, Phase II, and Phase III processes was performed. In addition, the following were discussed (1) regulation of xenobiotic metabolism; (2) review of published toxicokinetics for specific chemicals; and (3) genetic diversity of these processes in C. elegans. Finally, worm xenobiotic transport and metabolism was placed in an evolutionary context; key areas for future research highlighted; and implications for extrapolating C. elegans toxicity results to other species discussed.

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.).

Transcriptome analyses of Ditylenchus destructor in responses to cold and desiccation stress

The objective of this study was to identify molecular responses in Ditylenchus destructor to cold and desiccation by means of transcriptomes analyses. A total of 102,517 unigenes were obtained, with an average length of 1,076 bp, in which 58,453 (57%) had a functional annotation. A total of 1154 simple sequence repeats (SSRs) distributed over 1078 unigenes were detected. Gene expression profiles in response to cold and desiccation stress and the expression of specific stress-related genes were compared. Gene ontology analysis and pathway-based analysis were used to further investigate the functions of the differentially expressed genes. The reliability of the sequencing data was verified through quantitative real-time PCR analysis of 19 stress-related genes. RNA interference used to further assess the functions of the cold-related unigenes 15628 and 15596 showed that the knockdown of each of these genes led to decreased cold tolerance of D. destructor. Hence, this study revealed molecular processes and pathways active in cold- or dessication-treated nematodes. The transcriptome profiles presented in this study provide insight into the transcriptome complexity and will contribute to further understand stress tolerance in D. destructor.

Functional analysis of UGT201D3 associated with abamectin resistance in Tetranychus cinnabarinus (Boisduval)

Uridine diphosphate (UDP)-glycosyltransferases (UGTs) are widely distributed within living organisms and share roles in biotransformation of various lipophilic endo- and xenobiotics with activated UDP sugars. In this study, it was found that the activity of UGTs in abamectin-resistant (AbR) strain was significantly higher (2.35-fold) than that in susceptible strain (SS) of Tetranychus cinnabarinus. Further analysis showed that 5-nitrouracil, the inhibitor of UGTs, could enhance the lethal effect of abamectin on mites. From the previous microarray results, we found an UGT gene (UGT201D3) overexpressed in AbR strain. Quantitative PCR analysis showed that UGT201D3 was highly expressed and more inducible with abamectin exposure in the AbR strain. After silencing the transcription of UGT201D3, the activity of UGTs was decreased and the susceptibility to abamectin was increased in AbR strain whereas it was not in SS. Furthermore, UGT201D3 gene was then successfully expressed in Escherichia coli. The recombinant UGT201D3 exhibited α-naphthol activity (2.81 ± 0.43 nmol/mg protein/min), and the enzyme activity could be inhibited by abamectin (inhibitory concentration at 50%: 57.50 ± 3.54 μmol/L). High-performance liquid chromatography analysis demonstrated that the recombinant UGT201D3 could effectively deplete abamectin (15.77% ± 3.72%) incubating with 150 μg protein for 6 h. These results provided direct evidence that UGT201D3 was involved in abamectin resistance in T. cinnabarinus.

A Whole Genome Re-Sequencing Based GWA Analysis Reveals Candidate Genes Associated with Ivermectin Resistance in Haemonchus contortus

The most important and broad-spectrum drug used to control the parasitic worms to date is ivermectin (IVM). Resistance against IVM has emerged in parasites, and preserving its efficacy is now becoming a serious issue. The parasitic nematode Haemonchus contortus (Rudolphi, 1803) is economically an important parasite of small ruminants across the globe, which has a successful track record in IVM resistance. There are growing evidences regarding the multigenic nature of IVM resistance, and although some genes have been proposed as candidates of IVM resistance using lower magnification of genome, the genetic basis of IVM resistance still remains poorly resolved. Using the full magnification of genome, we herein applied a population genomics approach to characterize genome-wide signatures of selection among pooled worms from two susceptible and six ivermectin-resistant isolates of H. contortus, and revealed candidate genes under selection in relation to IVM resistance. These candidates also included a previously known IVM-resistance-associated candidate gene HCON_00148840, glc-3. Finally, an RNA-interference-based functional validation assay revealed the HCON_00143950 as IVM-tolerance-associated gene in H. contortus. The possible role of this gene in IVM resistance could be detoxification of xenobiotic in phase I of xenobiotic metabolism. The results of this study further enhance our understanding on the IVM resistance and continue to provide further evidence in favor of multigenic nature of IVM 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.

Intestinal UDP-glucuronosyltransferase as a potential target for the treatment and prevention of lymphatic filariasis

Lymphatic filariasis (LF), a morbid disease caused by the tissue-invasive nematodes Wuchereria bancrofti, Brugia malayi, and Brugia timori, affects millions of people worldwide. Global eradication efforts have significantly reduced worldwide prevalence, but complete elimination has been hampered by limitations of current anti-filarial drugs and the lack of a vaccine. The goal of this study was to evaluate B. malayi intestinal UDP-glucuronosyltransferase (Bm-UGT) as a potential therapeutic target. To evaluate whether Bm-UGT is essential for adult filarial worms, we inhibited its expression using siRNA. This resulted in a 75% knockdown of Bm-ugt mRNA for 6 days and almost complete suppression of detectable Bm-UGT by immunoblot. Reduction in Bm-UGT expression resulted in decreased worm motility for 6 days, 70% reduction in microfilaria release from adult worms, and significant reduction in adult worm metabolism as detected by MTT assays. Because prior allergic-sensitization to a filarial antigen would be a contraindication for its use as a vaccine candidate, we tested plasma from infected and endemic normal populations for Bm-UGT-specific IgE using a luciferase immunoprecipitation assay. All samples (n = 35) tested negative. We then tested two commercially available medicines known to be broad inhibitors of UGTs, sulfinpyrazone and probenecid, for in vitro activity against B. malayi. There were marked macrofilaricidal effects at concentrations achievable in humans and very little effect on microfilariae. In addition, we observed that probenecid and sulfinpyrazone exhibit a synergistic macrofilaricidal effect when used in combination with albendazole. The results of this study demonstrate that Bm-UGT is an essential protein for adult worm survival. Lack of prior IgE sensitization in infected and endemic populations suggest it may be a feasible vaccine candidate. The finding that sulfinpyrazone and probenecid have in vitro effects against adult B. malayi worms suggests that these medications have promise as potential macrofilaricides in humans.

The transcription factor NHR-8: A new target to increase ivermectin efficacy in nematodes

Resistance to the anthelmintic macrocyclic lactone ivermectin (IVM) has a great impact on the control of parasitic nematodes. The mechanisms by which nematodes adapt to IVM remain to be deciphered. We have identified NHR-8, a nuclear hormone receptor involved in the xenobiotic response in Caenorhabditis elegans, as a new regulator of tolerance to IVM. Loss-of-function nhr-8(ok186) C. elegans mutants subjected to larval development assays and electropharyngeogram measurements, displayed hypersensitivity to IVM, and silencing of nhr-8 in IVM-resistant worms increased IVM efficacy. In addition, compared to wild-type worms, nhr-8 mutants under IVM selection pressure failed to acquire tolerance to the drug. In addition, IVM-hypersensitive nhr-8(ok186) worms displayed low transcript levels of several genes from the xenobiotic detoxification network and a concomitant low Pgp-mediated drug efflux activity. Interestingly, some pgp and cyp genes known to impact IVM tolerance in many nematode species, were down regulated in nhr-8 mutants and inversely upregulated in IVM-resistant worms. Moreover, pgp-6 overexpression in nhr-8(ok186) C. elegans increased tolerance to IVM. Importantly, NHR-8 function was rescued in nhr-8(ok186) C. elegans with the homolog of the parasitic nematode Haemonchus contortus, and silencing of Hco-nhr-8 by RNAi on L2 H. contortus larvae increased IVM susceptibility in both susceptible and resistant H. contortus isolates. Thus, our data show that NHR-8 controls the tolerance and development of resistance to IVM in C. elegans and the molecular basis for this relates to the NHR-8-mediated upregulation of IVM detoxification genes. Since our results show that Hco-nhr-8 functions similarly to Cel-nhr-8, this study helps to better understand mechanisms underlying failure in drug efficacy and open perspectives in finding new compounds with NHR-8 antagonist activity to potentiate IVM efficacy.

IVM is the most important broad-spectrum deworming drug used today but resistance to this drug has appeared in parasites of both animals and humans. This seriously jeopardizes the success of current parasite control. Preserving IVM efficacy is a public health issue, whose outcome depends on the understanding of the molecular basis of selection for resistance to these drugs. We unambiguously show that the nuclear hormone receptor NHR-8, is crucial for protection of the nematode model Caenorhabditis elegans against IVM toxicity. Worms deficient in NHR-8 are hypersensitive to IVM and fail to become resistant to IVM under drug pressure. NHR-8 functions in the parasitic nematode of ruminants Haemonchus contortus and similar mechanisms could occur in other target pathogens. By controlling the xenobiotic detoxification network, NHR-8 may contribute to the biotransformation and elimination of IVM and help to desensitize the worm to the drug. This provides novel molecular targets involved in IVM drug tolerance in parasitic nematodes. Such findings could be exploited for targeted therapeutic intervention to treat parasitic nematode infections and delay the process of resistance development to anthelmintic drugs.

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

UDP-glycosyltransferases (UGT), catalysing conjugation of UDP-activated sugar donors to small lipophilic chemicals, are widespread in living organisms from bacteria to fungi, plant, or animals. The progress of genome sequencing has enabled an assessment of the UGT multigene family in Haemonchus contortus (family Trichostrongylidae, Nematoda), a hematophagous gastrointestinal parasite of small ruminants. Here we report 32 putative UGT genes divided into 15 UGT families. Phylogenetic analysis in comparison with UGTs from Caenorhabditis elegans, a free-living model nematode, revealed several single member homologues, a lack of the dramatic gene expansion seen in C. elegans, but also several families (UGT365, UGT366, UGT368) expanded in H. contortus only. The assessment of constitutive UGT mRNA expression in H. contortus adults identified significant differences between females and males. In addition, we compared the expression of selected UGTs in the drug-sensitive ISE strain to two benzimidazole-resistant strains, IRE and WR, with different genetic backgrounds. Constitutive expression of UGT368B2 was significantly higher in both resistant strains than in the sensitive strain. As resistant strains were able to deactivate benzimidazole anthelmintics via glycosylation more effectively then the sensitive strain, UGT368B2 enhanced constitutive expression might contribute to drug resistance in H. contortus.

The transcription factor SKN-1 and detoxification gene ugt-22 alter albendazole efficacy in Caenorhabditis elegans

Parasitic nematodes infect over 1/4 th of the human population and are a major burden on livestock and crop production. Benzimidazole class anthelmintics are widely used to treat infections, but resistance is a widespread problem. Mutation of genes encoding the benzimidazole target β-tubulin is a well-established mechanism of resistance, but recent evidence suggests that metabolism of the drugs may also occur. Our objective was to investigate contributions of the detoxification-response transcription factor SKN-1 to anthelmintic drug resistance using C. elegans. We find that skn-1 mutations alter EC₅₀ of the common benzimidazole albendazole in motility assays by 1.5–1.7 fold. We also identify ugt-22 as a detoxification gene associated with SKN-1 that influences albendazole efficacy. Mutation and overexpression of ugt-22 alter albendazole EC₅₀ by 2.3–2.5-fold. The influence of a nematode UGT on albendazole efficacy is consistent with recent studies demonstrating glucose conjugation of benzimidazoles.

Metabolism of albendazole, ricobendazole and flubendazole in Haemonchus contortus adults: Sex differences, resistance-related differences and the identification of new metabolites

Haemonchus contortus (family Trichostrongylidae, Nematoda), a hematophagous gastrointestinal parasite found in small ruminants, has a great ability to develop resistance to anthelmintic drugs. We studied the biotransformation of the three benzimidazole anthelmintics: albendazole (ABZ), ricobendazole (albendazole S-oxide; RCB) and flubendazole (FLU) in females and males of H. contortus in both a susceptible ISE strain and resistant IRE strain. The ex vivo cultivation of living nematodes in culture medium with or without the anthelmintics was used. Ultrasensitive UHPLC/MS/MS analysis revealed 9, 7 and 12 metabolites of ABZ, RCB and FLU, respectively, with most of these metabolites now described in the present study for the first time in H. contortus. The structure of certain metabolites shows the presence of biotransformation reactions not previously reported in nematodes. There were significant qualitative and semi-quantitative differences in the metabolites formed by male and female worms. In most cases, females metabolized drugs more extensively than males. Adults of the IRE strain were able to form many more metabolites of all the drugs than adults of the ISE strain. Some metabolites were even found only in adults of the IRE strain. These findings suggest that increased drug metabolism may play a role in resistance to benzimidazole drugs in H. contortus.

Pleiotropic alterations in gene expression in Latin American Fasciola hepatica isolates with different susceptibility to drugs

Background

Fasciola hepatica is the main agent of fasciolosis, a zoonotic disease affecting livestock worldwide, and an emerging food-borne disease in humans. Even when effective treatments are available, drugs are costly and can result in tolerance, liver damage and normally they do not prevent reinfection. Drug-resistant strains in livestock have been reported in various countries and, more worryingly, drug resistance in human cases has emerged in South America. The present study aims to characterize the transcriptome of two South American resistant isolates, the Cajamarca isolate from Peru, resistant to both triclabendazole and albendazole (TCBZR/ABZR) and the Rubino isolate from Uruguay, resistant to ABZ (TCBZS/ABZR), and compare them to a sensitive strain (Cenapa, Mexico, TCBZS/ABZS) to reveal putative molecular mechanisms leading to drug resistance.

Results

We observed a major reduction in transcription in the Cajamarca TCBZR/ABZR isolate in comparison to the other isolates. While most of the differentially expressed genes are still unannotated, several trends could be detected. Specific reduction in the expression levels of cytoskeleton proteins was consistent with a role of tubulins as putative targets of triclabendazole (TCBZ). A marked reduction of adenylate cyclase might be underlying pleiotropic effects on diverse metabolic pathways of the parasite. Upregulation of GST mu isoforms suggests this detoxifying mechanism as one of the strategies associated with resistance.

Conclusions

Our results stress the value of transcriptomic approaches as a means of providing novel insights to advance the understanding of drug mode of action and drug resistance. The results provide evidence for pleiotropic variations in drug-resistant isolates consistent with early observations of TCBZ and ABZ effects and recent proteomic findings.

Electronic supplementary material

The online version of this article (10.1186/s13071-017-2553-2) contains supplementary material, which is available to authorized users.

A molecularly imprinted polymer coupled with high-performance liquid chromatography-UV for the determination of albendazole in plasma and urine samples: CCD-RSM design

The present study attempted to develop a fast and sensitive ultrasound-assisted-dispersive-micro-solid phase extraction method for the separation and preconcentration of albendazole from plasma and urine samples.

Comparison of constitutive and thiabendazole-induced expression of five cytochrome P450 genes in fourth-stage larvae of Haemonchus contortus isolates with different drug susceptibility identifies one gene with high constitutive expression in a multi-resistant isolate

Benzimidazoles (BZs) remain amongst the most widely used anthelmintic drug classes against gastro-intestinal nematode infections, although their efficacy is increasingly compromised by resistance. The primary underlying mechanisms for BZ resistance are single-nucleotide polymorphisms (SNPs) in the isotype 1 β-tubulin gene causing the substitutions F167Y, E198A or F200Y. However, resistance is believed to be multi-genic and previous studies have shown that isolates carrying 90-100% F200Y can vary considerably in their resistance level in the egg hatch assay (EHA). Cytochrome P450 monooxygenases (CYPs) are associated with drug resistance in mammals and arthropods and have been considered as mediators of anthelmintic resistance. In Caenorhabditis elegans, several members of the CYP34/35 and CYP31 families are BZ and/or xenobiotic inducible and thiabendazole (TBZ) is metabolised by CYP35D1. Here, expression of all 5 CYPs closely related to the C. elegans CYP34/35 and CYP31 families was investigated in fourth-stage larvae of two susceptible and three BZ-resistant Haemonchus contortus isolates following in vitro exposure to TBZ for 3 and 6 h using real-time RT-PCR. The resistance status of all isolates was determined using EHAs and quantification of resistance-associated β-tubulin SNPs using pyrosequencing. While none of the CYPs was TBZ inducible, constitutive expression of CYP34/35 family member HCOI100383400 was significantly 2.4-3.7-fold higher in the multi-drug resistant WR isolate with the strongest BZ resistance phenotype compared to susceptible and intermediate-level BZ-resistant isolates. Although this increase is only moderate, HCOI100383400 might still be involved in high-level BZ resistance by further decreasing susceptibility in isolates already carrying 100% of a β-tubulin SNP causing BZ resistance. Lower transcript levels were observed for all CYPs in the intermediately resistant IRE isolate in comparison to the susceptible HcH isolate, which, except for CYP HCOI01579500, were statistically non-significant. This suggests that none of the investigated CYPs may contribute to protection against TBZ in this particular isolate.

L-Asparaginase of Leishmania donovani: Metabolic target and its role in Amphotericin B resistance

Emergence of Amphotericin B (AmB) resistant Leishmania donovani has posed major therapeutic challenge against the parasite. Consequently, combination therapy aimed at multiple molecular targets, based on proteome wise network analysis has been recommended. In this regard we had earlier identified and proposed L-asparaginase of Leishmania donovani (LdAI) as a crucial metabolic target. Here we report that both LdAI overexpressing axenic amastigote and promastigote forms of L. donovani survives better when challenged with AmB as compared to wild type strain. Conversely, qRT-PCR analysis showed an upregulation of LdAI in both forms upon AmB treatment. Our data demonstrates the importance of LdAI in imparting immediate protective response to the parasite upon AmB treatment. In the absence of structural and functional information, we modeled LdAI and validated its solution structure through small angle X-ray scattering (SAXS) analysis. We identified its specific inhibitors through ligand and structure-based approach and characterized their effects on enzymatic properties (K_m, V_max, K_cat) of LdAI. We show that in presence of two of the inhibitors L1 and L2, the survival of L. donovani is compromised whereas overexpression of LdAI in these cells restores viability. Taken together, our results conclusively prove that LdAI is a crucial metabolic enzyme conferring early counter measure against AmB treatment by Leishmania.

Deleterious Consequences of UDP-Galactopyranose Mutase Inhibition for Nematodes

Parasitic nematodes pose a serious threat to agriculture, livestock, and human health. Increasing resistance to antiparasitic agents underscores the need to replenish our anthelmintic arsenal. The nonpathogenic Caenorhabditis elegans, which serves as an effective model of parasitic helminths, has been used to search for new anthelmintic leads. We previously reported small-molecule inhibitors of the essential C. elegans protein UDP-galactopyranose mutase (UGM or Glf). This enzyme is required for the generation of galactofuranose (Galf)-containing glycans and is needed in nematodes for proper cuticle formation. Though our first-generation inhibitors were effective in vitro, they elicited no phenotypic effects. These findings are consistent with the known difficulty of targeting nematodes. C. elegans is recalcitrant to pharmacological modulation; typically, less than 0.02% of small molecules elicit a phenotypic effect, even at 40 μM. We postulated that the lack of activity of the UGM inhibitors was due to their carboxylic acid group, which can be exploited by nematodes for detoxification. We therefore tested whether replacement of the carboxylate with an N-acylsulfonamide surrogate would result in active compounds. UGM inhibitors with the carboxylate mimetic can phenocopy the deleterious consequences of UGM depletion in C. elegans. These findings support the use of UGM inhibitors as anthelmintic agents. They also outline a strategy to render small-molecule carboxylates more effective against nematodes.

Reduced efficacy of albendazole against Ascaris lumbricoides in Rwandan schoolchildren

Control of human soil-transmitted helminths (STHs) relies on preventive chemotherapy of schoolchildren applying the benzimidazoles (BZ) albendazole or mebendazole. Anthelmintic resistance (AR) is a common problem in nematodes of veterinary importance but for human STHs, information on drug efficacy is limited and routine monitoring is rarely implemented. Herein, the efficacy of single dose albendazole (400 mg) was evaluated in 12 schools in the Huye district of Rwanda where Ascaris is the predominant STH. Ascaris eggs were detected by wet mount microscopy and the Mini-FLOTAC method to assess cure rate (CR) and faecal egg count reduction (FECR). Blood and faecal samples were analysed for co-infections with Plasmodium sp. and Giardia duodenalis, respectively. Ascaris positive samples collected before and after treatment were analysed for putatively BZ-resistance associated β-tubulin gene single nucleotide polymorphisms. The overall CR was 69.9% by Mini-FLOTAC and 88.6% by wet mount microscopy. The FECR was 75.4% and the 95% calculated confidence intervals were 50.4–87.8% using sample variance, 55.4–88.8% by bootstrapping, and 75.0–75.7% applying a Markov Chain Monte Carlo Bayesian approach. FECR varied widely between 0 and 96.8% for individual schools. No putative BZ-resistance associated polymorphisms were found in the four Ascaris β-tubulin isotype genes examined. Since FECRs <95% indicate reduced efficacy, these findings raise the suspicion of BZ resistance. In the absence of respective molecular evidence, heritable AR in the local Ascaris populations cannot be formally proven. However, since FECRs <95% indicate reduced efficacy, BZ resistance may be suspected which would be alarming and calls for further analyses and routine monitoring in preventive chemotherapy programs.

Detection of Burkholderia pseudomallei toxin-mediated inhibition of protein synthesis using a Caenorhabditis elegans ugt-29 biosensor

Toxins are believed to play a crucial role in Burkholderia pseudomallei pathogenicity, however to date, only a few have been identified. The discovery of additional toxic molecules is limited by the lack of a sensitive indicator of B. pseudomallei toxicity. Previously, from a whole genome transcriptome analysis of B. pseudomallei-infected Caenorhabditis elegans, we noted significant overexpression of a number of worm genes encoding detoxification enzymes, indicating the host's attempt to clear bacterial toxic molecules. One of these genes, ugt-29, a family member of UDP-glucuronosyltransferases, was the most robustly induced phase II detoxification gene. In this study, we show that strong induction of ugt-29 is restricted to infections by the most virulent species among the pathogens tested. We also noted that ugt-29 is activated upon disruption of host protein synthesis. Hence, we propose that UGT-29 could be a promising biosensor to detect B. pseudomallei toxins that compromise host protein synthesis. The identification of bactobolin, a polyketide-peptide hybrid molecule, as a toxic molecule of B. pseudomallei further verifies the utilization of this surveillance system to search for bacterial toxins. Hence, a ugt-29 based reporter should be useful in screening for other molecules that inhibit host protein synthesis.

Haemonchus contortus: Genome Structure, Organization and Comparative Genomics

One of the first genome sequencing projects for a parasitic nematode was that for Haemonchus contortus. The open access data from the Wellcome Trust Sanger Institute provided a valuable early resource for the research community, particularly for the identification of specific genes and genetic markers. Later, a second sequencing project was initiated by the University of Melbourne, and the two draft genome sequences for H. contortus were published back-to-back in 2013. There is a pressing need for long-range genomic information for genetic mapping, population genetics and functional genomic studies, so we are continuing to improve the Wellcome Trust Sanger Institute assembly to provide a finished reference genome for H. contortus. This review describes this process, compares the H. contortus genome assemblies with draft genomes from other members of the strongylid group and discusses future directions for parasite genomics using the H. contortus model.

The Role of Xenobiotic-Metabolizing Enzymes in Anthelmintic Deactivation and Resistance in Helminths

Xenobiotic-metabolizing enzymes (XMEs) modulate the biological activity and behavior of many drugs, including anthelmintics. The effects of anthelmintics can often be abolished by XMEs when the drugs are metabolized to an inefficient compound. XMEs therefore play a significant role in anthelmintic efficacy. Moreover, differences in XMEs between helminths are reflected by differences in anthelmintic metabolism between target species. Taking advantage of the newly sequenced genomes of many helminth species, progress in this field has been remarkable. The present review collects up to date information regarding the most important XMEs (phase I and phase II biotransformation enzymes; efflux transporters) in helminths. The participation of these XMEs in anthelmintic metabolism and their possible roles in drug resistance are evaluated.

Reliable reference gene selection for quantitative real time PCR in Haemonchus contortus

The aim of this work was to identify reliable reference genes for expression studies in adult Haemonchus contortus. Eleven candidate genes were identified and the stability of their expression was assessed in adult males and females of two genetically divergent H. contortus isolates: drug-susceptible (ISE) and multi-drug-resistant (WR). Five genes with the most stable expression patterns were further assessed for suitability as reference genes in anthelmintic-treated H. contortus adults versus non-treated controls. We identified important differences in the expression of a number of candidate genes in anthelmintic-treated samples, confirming the need for careful validation of control genes for such experiments. We propose the use of multiple reference genes for expression studies in this species and found gpd, ama and far most suitable for adult H. contortus.

NHR-176 regulates cyp-35d1 to control hydroxylation-dependent metabolism of thiabendazole in Caenorhabditis elegans

Knowledge of how drugs are metabolized and excreted is an essential component of understanding their fate within and among target and non-target organisms. Thiabendazole (TBZ) was the first benzimidazole (BZ) to be commercially available and remains one of the most important anthelmintic drugs for medical and veterinary use. We have characterized how Caenorhabditis elegans metabolizes and excretes TBZ. We have shown that TBZ directly binds to the nuclear hormone receptor (NHR)-176 and that this receptor is required for the induction by TBZ of the cytochrome P450 (CYP) encoded by cyp-35d1. Further, RNAi inhibition of cyp-35d1 in animals exposed to TBZ causes a reduction in the quantity of a hydroxylated TBZ metabolite and its glucose conjugate that is detected in C. elegans tissue by HPLC. This final metabolite is unique to nematodes and we also identify two P-glycoproteins (PGPs) necessary for its excretion. Finally, we have shown that inhibiting the metabolism we describe increases the susceptibility of C. elegans to TBZ in wild-type and in resistant genetic backgrounds.

The cytochrome P450 family in the parasitic nematode Haemonchus contortus

Haemonchus contortus, a highly pathogenic and economically important parasitic nematode of sheep, is particularly adept at developing resistance to the anthelmintic drugs used in its treatment and control. The basis of anthelmintic resistance is poorly understood for many commonly used drugs with most research being focused on mechanisms involving drug targets or drug efflux. Altered or increased drug metabolism is a possible mechanism that has yet to receive much attention despite the clear role of xenobiotic metabolism in pesticide resistance in insects. The cytochrome P450s (CYPs) are a large family of drug-metabolising enzymes present in almost all living organisms, but for many years thought to be absent from parasitic nematodes. In this paper, we describe the CYP sequences encoded in the H. contortus genome and compare their expression in different parasite life-stages, sexes and tissues. We developed a novel real-time PCR approach based on partially assembled CYP sequences "tags" and confirmed findings in the subsequent draft genome with RNA-seq. Constitutive expression was highest in larval stages for the majority of CYPs, although higher expression was detected in the adult male or female for a small subset of genes. Many CYPs were expressed in the worm intestine. A number of H. contortus genes share high identity with Caenorhabditis elegans CYPs and the similarity in their expression profiles supports their classification as putative orthologues. Notably, H. contortus appears to lack the dramatic CYP subfamily expansions seen in C. elegans and other species, which are typical of CYPs with exogenous roles. However, a small group of H. contortus genes cluster with the C. elegans CYP34 and CYP35 subfamilies and may represent candidate xenobiotic metabolising genes in the parasite.

Development of an in vitro screen for compound bioaccumulation in Haemonchus contortus

The objective of the current study was to establish an in vitro screen and a highly sensitive analytical assay to delineate key physicochemical properties that favor compound bioaccumulation in the L3 life stage of a Haemonchus contortus isolate. Time-dependent studies revealed that absorption and elimination kinetics during the first 6 hr of exposure were sufficient to achieve maximum bioaccumulation for the majority of compounds tested. In subsequent studies, the larvae were incubated for 6 hr in a medium containing 146 compounds (5 μM initial concentration), including both human and veterinary medicines, characterized by a broad range of physicochemical properties. Bioaccumulation of the compounds by the nematodes was determined, and multiple physicochemical descriptors were selected for correlation. Data analysis using Bayes classification model and partial least-square regression revealed that clogD7.4, rotatable bond, E-state, and hydrogen bond donor each correlated with compound bioaccumulation in H. contortus L3. The finding that lipophilicity was critical for transcuticle compound permeation was consistent with previous studies in other parasitic species and in adult H. contortus . The finding of additional physicochemical properties that contribute to compound conformational flexibility, polarity, and electrotopological state shed light on the mechanisms governing transcuticle permeation. The relatively poor correlation between transcuticle and transmembrane permeation indicated the distinct mechanisms of compound permeation, likely due to the different constituents, and their contributions to overall transport function, of the lipid membranes and the porous collagen barrier of the nematode cuticle. Our study, for the first time, establishes a high-throughput screen for compound bioaccumulation in a parasitic nematode and further elucidates physicochemical factors governing transcuticular permeation of compounds. Application of this methodology will help explain the basis for discrepancies observed in receptor binding and whole organism potency assays and facilitate incorporation of drug delivery principles in the design of candidate anthelmintics.

Phenanthrene bioaccumulation in the nematode Caenorhabditis elegans

The contribution of food to the bioaccumulation of xenobiotics and hence toxicity is still an ambiguous issue. It is becoming more and more evident that universal statements cannot be made, but that the relative contribution of food-associated xenobiotics in bioaccumulation depends on species, substance, and environmental conditions. Yet, small-sized benthic or soil animals such as nematodes have largely been disregarded so far. Bioaccumulation of the polycyclic aromatic hydrocarbon phenanthrene in the absence and presence of bacterial food was measured in the nematode Caenorhabditis elegans. Elimination of phenanthrene in the nematodes was biphasic, suggesting that there was a slowly exchanging pool within the nematodes or that biotransformation of phenanthrene took place. Even with food present, dissolved phenanthrene was still the major contributor to bioaccumulated compound in nematode tissues, whereas the diet only contributed about 9%. Toxicokinetic parameters in the treatment without food were different from the ones of the treatment with bacteria, possibly because nematodes depleted their lipid reserves during starvation.

Phenobarbital induction and chemical synergism demonstrate the role of UDP-glucuronosyltransferases in detoxification of naphthalophos by Haemonchus contortus larvae

We used an enzyme induction approach to study the role of detoxification enzymes in the interaction of the anthelmintic compound naphthalophos with Haemonchus contortus larvae. Larvae were treated with the barbiturate phenobarbital, which is known to induce the activity of a number of detoxification enzymes in mammals and insects, including cytochromes P450 (CYPs), UDP-glucuronosyltransferases (UDPGTs), and glutathione (GSH) S-transferases (GSTs). Cotreatment of larvae with phenobarbital and naphthalophos resulted in a significant increase in the naphthalophos 50% inhibitory concentration (IC50) compared to treatment of larvae with the anthelmintic alone (up to a 28-fold increase). The phenobarbital-induced drug tolerance was reversed by cotreatment with the UDPGT inhibitors 5-nitrouracil, 4,6-dihydroxy-5-nitropyrimidine, probenecid, and sulfinpyrazone. Isobologram analysis of the interaction of 5-nitrouracil with naphthalophos in phenobarbital-treated larvae clearly showed the presence of strong synergism. The UDPGT inhibitors 5-nitrouracil, 4,6-dihydroxy-5-nitropyrimidine, and probenecid also showed synergistic effects with non-phenobarbital-treated worms (synergism ratio up to 3.2-fold). This study indicates that H. contortus larvae possess one or more UDPGT enzymes able to detoxify naphthalophos. In highlighting the protective role of this enzyme group, this study reveals the potential for UDPGT enzymes to act as a resistance mechanism that may develop under drug selection pressure in field isolates of this species. In addition, the data indicate the potential for a chemotherapeutic approach utilizing inhibitors of UDPGT enzymes as synergists to increase the activity of naphthalophos against parasitic worms and to combat detoxification-mediated drug resistance if it arises in the field.

Analysis of putative inhibitors of anthelmintic resistance mechanisms in cattle gastrointestinal nematodes

Effects of the cytochrome P450 inhibitor piperonyl butoxide and the P-glycoprotein inhibitor verapamil on the efficacy of ivermectin and thiabendazole were studied in vitro in susceptible and resistant isolates of the cattle parasitic nematodes Cooperia oncophora and Ostertagia ostertagi. The effects of combined use of drug and piperonyl butoxide/verapamil, respectively, were investigated in the Egg Hatch Assay, the Larval Development Assay and the Larval Migration Inhibition Assay. The effects of piperonyl butoxide and verapamil as inhibitors of thiabendazole and ivermectin responses were particularly marked for larval development, where both inhibitors were able to completely eliminate all differences between susceptible and resistant isolates. Even the lowest concentrations of anthelmintics used in combination with inhibitors caused complete inhibition of development. Differences and/or similarities among responses in different isolates were only obtained in the two other assays: in the Egg Hatch Assay piperonyl butoxide caused a shift in concentration-response curves obtained with thiabendazole to the left for all isolates tested, changing relative differences between isolates. In contrast, an effect of verapamil in the Egg Hatch Assay was only apparent for benzimidazole-resistant isolates. In the Larval Migration Inhibition Assay only ivermectin was tested and piperonyl butoxide shifted the concentration-response curves for all isolates to the left, again eliminating differences in EC50 values between susceptible and resistant isolates. This was not the case using verapamil as an inhibitor, where curves for both susceptible and benzimidazole-resistant isolates shifted to the left in Ostertagia isolates. In Cooperia the picture was more complex with ivermectin-resistant isolates showing a larger shift than the susceptible isolate. Single nucleotide polymorphisms in the β-tubulin isotype 1 gene were investigated. Significantly increased frequencies of resistance-associated alleles were observed for the codons 167 and 200 in one benzimidazole-resistant isolate but not in an isolate selected for benzimidazole resistance at an early stage of selection.

Adaptive and specialised transcriptional responses to xenobiotic stress in Caenorhabditis elegans are regulated by nuclear hormone receptors

Characterisation of the pathways by which xenobiotics are metabolised and excreted in both target and non-target organisms is crucial for the rational design of effective and specific novel bioactive molecules. Consequently, we have investigated the induced responses of the model nematode Caenorhabditis elegans to a variety of xenobiotics which represent a range of putative modes of action. The majority of genes that were specifically induced in preliminary microarray analyses encoded enzymes from Phase I and II metabolism, including cytochrome P450s, short chain dehydrogenases, UDP-glucuronosyl transferases and glutathione transferases. Changes in gene expression were confirmed by quantitative PCR and GFP induction in reporter strains driven by promoters for transcription of twelve induced enzymes was investigated. The particular complement of metabolic genes induced was found to be highly contingent on the xenobiotic applied. The known regulators of responses to applied chemicals ahr-1, hif-1, mdt-15 and nhr-8 were not required for any of these inducible responses and skn-1 regulated GFP expression from only two of the promoters. Reporter strains were used in conjunction with systematic RNAi screens to identify transcription factors which drive expression of these genes under xenobiotic exposure. These transcription factors appeared to regulate specific xenobiotic responses and have no reported phenotypes under standard conditions. Focussing on nhr-176 we demonstrate the role of this transcription factor in mediating the resistance to thiabendazole.

Biotransformation of albendazole and activities of selected detoxification enzymes in Haemonchus contortus strains susceptible and resistant to anthelmintics

The increased activity of drug-metabolizing enzymes can protect helminths against the toxic effect of anthelmintics. The aim of this study was to compare the metabolism of the anthelmintic drug albendazole (ABZ) and the activities of selected biotransformation and antioxidant enzymes in three different strains of Haemonchus contortus: the ISE strain (susceptible to common anthelmintics), the BR strain (resistant to benzimidazole anthelmintics) and the WR strain (multi-resistant). H. contortus adults were collected from the abomasum of experimentally infected lambs. In vitro (subcellular fractions of H. contortus homogenate) as well as ex vivo (living nematodes cultivated in flasks with medium) experiments were performed. HPLC with spectrofluorimetric and mass-spectrometric detection was used in the analysis of ABZ metabolites. The in vitro activities of oxidation/antioxidation and conjugation enzymes toward model substrates were also assayed. The in vitro data showed significant differences between the susceptible (ISE) and resistant (BR, WR) strains regarding the activities of peroxidases, catalase and UDP-glucosyltransferases. S-oxidation of ABZ was significantly lower in BR than in the ISE strain. Ex vivo, four ABZ metabolites were identified: ABZ sulphoxide and three ABZ glucosides. In the resistant strains BR and WR, the ex vivo formation of all ABZ glucosides was significantly higher than in the susceptible ISE strain. The altered activities of certain detoxifying enzymes might partly protect the parasites against the toxic effect of the drugs as well as contribute to drug-resistance in these parasites.

Chemical detoxification of small molecules by Caenorhabditis elegans

Caenorhabditis elegans lives in compost and decaying fruit, eats bacteria and is exposed to pathogenic microbes. We show that C. elegans is able to modify diverse microbial small-molecule toxins via both O- and N-glucosylation as well as unusual 3'-O-phosphorylation of the resulting glucosides. The resulting glucosylated derivatives have significantly reduced toxicity to C. elegans, suggesting that these chemical modifications represent a general mechanism for worms to detoxify their environments.