The role of several ABC transporter genes in ivermectin resistance in Caenorhabditis elegans

Biocontrol Potential of Entomopathogenic Fungi Against Plant-Parasitic Nematodes: A Caenorhabditis elegans-Based Screening and Mechanistic Study

Plant-parasitic nematodes and insect pests critically threaten agricultural productivity, but chemical pesticides face limitations due to resistance and environmental concerns, necessitating eco-friendly biopesticides targeting both pests and nematodes. Here, we developed a high-throughput screening platform using Caenorhabditis elegans to identify entomopathogenic fungi exhibiting nematocidal activity against Meloidogyne incognita. Among 32 tested strains, nine Metarhizium spp. and one Beauveria strain demonstrated dual efficacy against C. elegans and M. incognita. Metarhizium anisopliae CQMa421 showed the highest virulence, suppressing nematode reproduction by 42.7% and inducing >80% mortality. Pot experiments revealed a 50% reduction in the root galling index and 50.3% fewer root galls in Solanum lycopersicum. The CQMa421 filtrate caused irreversible locomotor deficits and reduced egg hatching rates by 28%. Concurrently, intestinal damage, elevated oxidative stress and autophagy were observed in C. elegans. This was accompanied by a transcriptome-wide modulation of genes involved in detoxification and immune defense pathways. These findings demonstrate the efficacy of our C. elegans-based screening method for identifying fungi with nematocidal potential. CQMa421's virulence against M. incognita suggests its promise for pest management, while molecular insights highlight pathways that may contribute to the future design of future nematicides. This study advances fungal biocontrol agents and offers a sustainable strategy for agriculture.

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.

Expression of transporter genes in anthelmintic resistant isolates of Haemonchus contortus

ATP-binding cassette (ABC) transporters, including P-glycoproteins (PGP), have been implicated in drug resistance in different organisms including Haemonchus contortus. This study confirmed the resistance status of H. contortus isolates selected for ivermectin (IVM) and oxfendazole (OXF) resistances using the fecal egg count reduction test and evaluated the gene expression of seven ABC transporters using RT-qPCR for two biological scenarios: the effect of selection for anthelmintic resistance and the effect of drug exposure on gene expression. Gene expression results showed that selection for IVM resistance led to the significant upregulation of Hco-pgp-9a (1.5-fold), Hco-pgp-11 (3-fold) and Hco-haf-9 (1.5-fold) (p < 0.05). Similarly, selection for OXF resistance led to the significant upregulation of Hco-pgp-9a (3-fold), Hco-pgp-11 (4-fold) and Hco-haf-9 (2-fold) when comparing with the unselected ISE isolate (p < 0.05). Exposure of selected isolates to anthelmintics lead to no significant upregulation of the studied transporter genes. We also observed instances where there was strong intragroup variation regarding samples originating from parasites obtained from different individual hosts pointing that the interactions of the animal host with the tested anthelmintics may also play a role in the expression of the studied nematode genes.

Importance of ABC Transporters in the Survival of Parasitic Nematodes and the Prospect for the Development of Novel Control Strategies

ABC transporters, a family of ATP-dependent transmembrane proteins, are responsible for the active transport of a wide range of molecules across cell membranes, including drugs, toxins, and nutrients. Nematodes possess a great diversity of ABC transporters; however, only P-glycoproteins have been well-characterized compared to other classes. The ABC transport proteins have been implicated in developing resistance to various classes of anthelmintic drugs in parasitic nematodes; their role in plant and human parasitic nematodes still needs further investigation. Therefore, ABC transport proteins offer a potential opportunity to develop nematode control strategies. Multidrug resistance inhibitors are becoming more attractive for controlling nematodes due to their potential to increase drug efficacy in two ways: (i) by limiting drug efflux from nematodes, thereby increasing the amount of drug that reaches its target site, and (ii) by reducing drug excretion by host animals, thereby enhancing drug bioavailability. This article reviews the role of ABC transporters in the survival of parasitic nematodes, including the genes involved, their regulation and physiological roles, as well as recent developments in their characterization. It also discusses the association of ABC transporters with anthelmintic resistance and the possibility of targeting them with next-generation inhibitors or nutraceuticals (e.g., polyphenols) to control parasitic infections.

Transcriptional response of Meloidogyne incognita to non-fumigant nematicides

There is limited research about the impacts of new nematicides, including fluazaindolizine, fluopyram, and fluensulfone, on the plant-parasitic nematode Meloidogyne incognita, despite it being a pervasive agricultural pest. In this study, M. incognita second-stage juveniles were exposed for 24-h to fluensulfone, fluazaindolizine, fluopyram, and oxamyl and total RNA was extracted and sequenced using next-generation sequencing to determine gene expression. The effects of nematicide exposure on cellular detoxification pathways, common differentially expressed (DE) genes, and fatty acid and retinol-binding genes were examined. Fluopyram and oxamyl had the smallest impacts on the M. incognita transcriptome with 48 and 151 genes that were DE, respectively. These compounds also elicited a weak response in the cellular detoxification pathway and fatty acid and retinol-binding (FAR) genes. Fluensulfone and fluazaindolizine produced robust transcriptional responses with 1208 and 2611 DE genes, respectively. These compounds had strong impacts on cellular detoxification, causing differential regulation of transcription factors and genes in the detox pathway. These compounds strongly down-regulated FAR genes between 52-85%. Having a greater understanding of how these compounds function at a molecular level will help to promote proper stewardship, aid with nematicide discovery, and help to stay a step ahead of nematicide resistance.

Stage-specific miRNAs regulate gene expression associated with growth, development and parasite-host interaction during the intra-mammalian migration of the zoonotic helminth parasite Fasciola hepatica

Background

MiRNAs are small non-coding RNAs that post-transcriptionally regulate gene expression in organisms ranging from viruses to mammals. There is great relevance in understanding how miRNAs regulate genes involved in the growth, development, and maturation of the many parasitic worms (helminths) that together afflict more than 2 billion people.

Results

Here, we describe the miRNAs expressed by each of the predominant intra-mammalian development stages of Fasciola hepatica, a foodborne flatworm that infects a wide range of mammals worldwide, most importantly humans and their livestock. A total of 124 miRNAs were profiled, 72 of which had been previously reported and three of which were conserved miRNA sequences described here for the first time. The remaining 49 miRNAs were novel sequences of which, 31 were conserved with F. gigantica and the remaining 18 were specific to F. hepatica. The newly excysted juveniles express 22 unique miRNAs while the immature liver and mature bile duct stages each express 16 unique miRNAs. We discovered several sequence variant miRNAs (IsomiRs) as well as miRNA clusters that exhibit strict temporal expression paralleling parasite development. Target analysis revealed the close association between miRNA expression and stage-specific changes in the transcriptome; for example, we identified specific miRNAs that target parasite proteases known to be essential for intestinal wall penetration (cathepsin L3). Moreover, we demonstrate that miRNAs fine-tune the expression of genes involved in the metabolic pathways that allow the parasites to move from an aerobic external environment to the anerobic environment of the host.

Conclusions

These results provide novel insight into the regulation of helminth parasite development and identifies new genes and miRNAs for therapeutic development to limit the virulence and pathogenesis caused by F. hepatica.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08644-z.

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.

Transporter gene expression and Wolbachia quantification in adults of Dirofilaria immitis treated in vitro with ivermectin or moxidectin alone or in combination with doxycycline for 12 hours

Due to their marked larvicidal activity, macrocyclic lactones (MLs) are used for the prevention of heartworm disease ( Dirofilaria immitis) in dogs. They have also been shown to eliminate adult parasites after long-term administration, with a so-called "slow-kill" effect. In addition, recent studies have established that a combination of doxycycline, which eliminates the endosymbiont Wolbachia, and MLs has superior adulticide effects when compared to MLs alone. It has been hypothesized that the apparent synergism between doxycycline/MLs may be due to interaction with drug efflux transport proteins. The aim of the present study was to evaluate gene expression of several transport proteins in D. immitis adults treated in vitro either with doxycycline alone, ivermectin alone, moxidectin alone, or a combination of ivermectin or moxidectin with doxycycline for 12h. Quantitative PCR analysis showed a sex-dependent response to treatments. In female worms, Dim-pgp-10, Dim-haf-1 and Dim-haf-5 were upregulated compared to controls with doxycycline alone and when combined with ivermectin. Moxidectin did not induce any changes in gene expression. In males, moxidectin administered alone induced a slight increase in Dim-pgp-10, Dim-pgp-11and Di-avr-14, while ivermectin in combination with doxycycline produced significant upregulation of the ML receptor Di-avr-14. These results suggest possible synergism between the two drug classes and different susceptibility of males vs. females to adulticide effects.

Crowdsourced Identification of Potential Target Genes for CTV Induced Gene Silencing for Controlling the Citrus Greening Vector Diaphorina citri

Citrus Greening or Huanglongbing (HLB) is a disease of citrus, causing high reduction in citrus production and is transmitted by the Asian citrus psyllid Diaphorina citri Kuwayama vectoring a phloem-limited bacterium Candidatus Liberibacter sp. We report research results using crowdsourcing challenge strategy identifying potential gene targets in D. citri to control the insect using RNA interference (RNAi). From 63 submitted sequences, 43 were selected and tested by feeding them to D. citri using artificial diet assays. After feeding on artificial diet, the three most effective dsRNAs causing 30% mortality above control silenced genes expressing iron-sulfur cluster subunit of the mitochondrial electron transport chain complex (Rieske), heme iron-binding terminal oxidase enzyme (Cytochrome P450) and tetrahydrobiopterin (BH₄) pathway enzyme (Pterin 4α-Carbinolamine Dehydratase). These sequences were cloned into a citrus phloem-limited virus (Citrus tristeza virus, CTV T36) expressing dsRNA against these target genes in citrus. The use of a viral mediated “para-transgenic” citrus plant system caused higher mortality to adult D. citri than what was observed using artificial diet, reaching 100% when detached citrus leaves with the engineered CTV expressing dsRNA were fed to adult D. citri. Using this approach, a virus-induced gene silencing (VIGS) can be used to test future transgenic cultivars before genetically engineering citrus. RNA Seq analysis after feeding D. citri CTV-RIE on infected leaves identified transcriptionally modified genes located upstream and downstream of the targeted RIE gene. These genes were annotated showing that many are associated with the primary function of the Rieske gene that was targeted by VIGS.

Genomic Epidemiology in Filarial Nematodes: Transforming the Basis for Elimination Program Decisions

Onchocerciasis and lymphatic filariasis are targeted for elimination, primarily using mass drug administration at the country and community levels. Elimination of transmission is the onchocerciasis target and global elimination as a public health problem is the end point for lymphatic filariasis. Where program duration, treatment coverage, and compliance are sufficiently high, elimination is achievable for both parasites within defined geographic areas. However, transmission has re-emerged after apparent elimination in some areas, and in others has continued despite years of mass drug treatment. A critical question is whether this re-emergence and/or persistence of transmission is due to persistence of local parasites-i.e., the result of insufficient duration or drug coverage, poor parasite response to the drugs, or inadequate methods of assessment and/or criteria for determining when to stop treatment-or due to re-introduction of parasites via human or vector movement from another endemic area. We review recent genetics-based research exploring these questions in Onchocerca volvulus, the filarial nematode that causes onchocerciasis, and Wuchereria bancrofti, the major pathogen for lymphatic filariasis. We focus in particular on the combination of genomic epidemiology and genome-wide associations to delineate transmission zones and distinguish between local and introduced parasites as the source of resurgence or continuing transmission, and to identify genetic markers associated with parasite response to chemotherapy. Our ultimate goal is to assist elimination efforts by developing easy-to-use tools that incorporate genetic information about transmission and drug response for more effective mass drug distribution, surveillance strategies, and decisions on when to stop interventions to improve sustainability of elimination.

Differential ABC transporter gene expression in adult Dirofilaria immitis males and females following in vitro treatment with ivermectin, doxycycline or a combination of both

Background

Combination doxycycline/macrocyclic lactone (ML) protocols have been shown to provide a more rapid adulticidal and microfilaricidal effect than either MLs or doxycycline alone, although female worms were reported to have a higher tolerance to treatments compared to male worms. The present study aimed to evaluate how ABC transporters may be involved in the synergic effect of the combination treatment. Adult worms of D. immitis were treated in vitro for 24 hours with doxycycline (DOXY), ivermectin (IVM) and a combination of both, and changes in the modulation of ABC transporter genes were measured. Levels of doxycycline inside different treatment media, post-treatment, were determined through HPLC analysis.

Results

Quantitative RT-PCR analysis showed the presence of changes in the modulation of ABC transporter genes evaluated in this study. In particular, in female worms, the combination treatment induced a substantial increase in gene expressions, especially of Dim-pgp-10 and Dim-haf-4; whereas in male worms, the greatest increase in gene expression was observed for Dim-pgp-10 and Dim-pgp-11 when treated with DMSO + IVM and DMSO + DOXY/IVM. HPLC analysis of the DOXY concentrations in the media after in vitro treatments of male worms showed a slight difference between the DMSO + DOXY samples and the combination (DMSO + DOXY + IVM), while no difference was observed among females.

Conclusions

Further studies are required to explain whether the modulation of cellular efflux plays a role, even partially, in the adulticide effect of doxycycline/macrocyclic lactone combinations in heartworm-infected dogs. To the authors’ knowledge, this is the first study to evaluate P-gp expression in adult D. immitis.

Perspectives on the utility of moxidectin for the control of parasitic nematodes in the face of developing anthelmintic resistance

Macrocyclic lactone (ML) anthelmintics are the most important class of anthelmintics because of our high dependence on them for the control of nematode parasites and some ectoparasites in livestock, companion animals and in humans. However, resistance to MLs is of increasing concern. Resistance is commonplace throughout the world in nematode parasites of small ruminants and is of increasing concern in horses, cattle, dogs and other animals. It is suspected in Onchocerca volvulus in humans. In most animals, resistance first arose to the avermectins, such as ivermectin (IVM), and subsequently to moxidectin (MOX). Usually when parasite populations are ML-resistant, MOX is more effective than avermectins. MOX may have higher intrinsic potency against some parasites, especially filarial nematodes, than the avermectins. However, it clearly has a significantly different pharmacokinetic profile. It is highly distributed to lipid tissues, less likely to be removed by ABC efflux transporters, is poorly metabolized and has a long half-life. This results in effective concentrations persisting for longer in target hosts. It also has a high safety index. Limited data suggest that anthelmintic resistance may be overcome, at least temporarily, if a high concentration can be maintained at the site of the parasites for a prolonged period of time. Because of the properties of MOX, there are reasonable prospects that strains of parasites that are resistant to avermectins at currently recommended doses will be controlled by MOX if it can be administered at sufficiently high doses and in formulations that enhance its persistence in the host. This review examines the properties of MOX that support this contention and compares them with the properties of other MLs. The case for using MOX to better control ML-resistant parasites is summarised and some outstanding research questions are presented.

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

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

Anthelmintic drug actions in resistant and susceptible C. elegans revealed by electrophysiological recordings in a multichannel microfluidic device

Many anthelmintic drugs used to treat parasitic nematode infections target proteins that regulate electrical activity of neurons and muscles: ion channels (ICs) and neurotransmitter receptors (NTRs). Perturbation of IC/NTR function disrupts worm behavior and can lead to paralysis, starvation, immune attack and expulsion. Limitations of current anthelmintics include a limited spectrum of activity across species and the threat of drug resistance, highlighting the need for new drugs for human and veterinary medicine. Although ICs/NTRs are valuable anthelmintic targets, electrophysiological recordings are not commonly included in drug development pipelines. We designed a medium-throughput platform for recording electropharyngeograms (EPGs)-the electrical signals emitted by muscles and neurons of the pharynx during pharyngeal pumping (feeding)-in Caenorhabditis elegans and parasitic nematodes. The current study in C. elegans expands previous work in several ways. Detecting anthelmintic bioactivity in drugs, compounds or natural products requires robust, sustained pharyngeal pumping under baseline conditions. We generated concentration-response curves for stimulating pumping by perfusing 8-channel microfluidic devices (chips) with the neuromodulator serotonin, or with E. coli bacteria (C. elegans' food in the laboratory). Worm orientation in the chip (head-first vs. tail-first) affected the response to E. coli but not to serotonin. Using a panel of anthelmintics-ivermectin, levamisole and piperazine-targeting different ICs/NTRs, we determined the effects of concentration and treatment duration on EPG activity, and successfully distinguished control (N2) and drug-resistant worms (avr-14; avr-15; glc-1, unc-38 and unc-49). EPG recordings detected anthelmintic activity of drugs that target ICs/NTRs located in the pharynx as well as at extra-pharyngeal sites. A bus-8 mutant with enhanced permeability was more sensitive than controls to drug treatment. These results provide a useful framework for investigators who would like to more easily incorporate electrophysiology as a routine component of their anthelmintic research workflow.

Dominance of P-glycoprotein 12 in phenotypic resistance conversion against ivermectin in Caenorhabditis elegans

While diseases caused by nematodes remains a considerable drawback for the livestock, agriculture and public health, anthelmintics drug resistance has been observed over the past years and is a major concern for parasite control. Ivermectin, initially considered as a highly potent drug, currently presents a reduced anti-helminthic efficacy, which is influenced by expression of several ATP-binding cassette transporters (ABC), among them the P-glycoproteins (Pgps). Here we present some evidences of Pgps dominance during Ivermectin resistance/susceptibility using Pgps double silencing in C. elegans and the phylogenetic relationship of Pgps among nematodes, which strengthen the use of this model for study of drug resistance in nematodes. Firstly, we evaluated the quantitative gene expression of 12 out the 15 known Pgps from resistant and WT strains of C. elegans, we demonstrated the upregulation of Pgps 12 and 13 and downregulation of all remaining Pgps in ivermectin resistant strain. By using an RNAi loss-of-function approach we observed that Pgp 12 gene silencing reverts the resistance phenotype to ivermectin, while Pgp 4 gene silencing does not alter the resistance phenotype but induces a resistance in wild type strain. Interestingly, the dual silencing of Pgp 12 and Pgp 4 expression demonstrates the dominance of phenotype promoted by Pgp 12 silencing. Finally, in silico analysis reveals a close relationship between Pgps from C. elegans and several nematodes parasites. Taken together, our results indicate that Pgp 12 is crucial for the resistance to ivermectin and thus a good candidate for further studies aiming to develop specific inhibitors to this transporter, allowing the continuous use of ivermectin to control the burden on animal and human health inflicted by nematode parasites globally.

Polymorphism in ABC transporter genes of Dirofilaria immitis

Dirofilaria immitis, a filarial nematode, causes dirofilariasis in dogs, cats and occasionally in humans. Prevention of the disease has been mainly by monthly use of the macrocyclic lactone (ML) endectocides during the mosquito transmission season. Recently, ML resistance has been confirmed in D. immitis and therefore, there is a need to find new classes of anthelmintics. One of the mechanisms associated with ML resistance in nematodes has been the possible role of ATP binding cassette (ABC) transporters in reducing drug concentrations at receptor sites. ABC transporters, mainly from sub-families B, C and G, may contribute to multidrug resistance (MDR) by active efflux of drugs out of the cell. Gene products of ABC transporters may thus serve as the targets for agents that may modulate susceptibility to drugs, by inhibiting drug transport. ABC transporters are believed to be involved in a variety of physiological functions critical to the parasite, such as sterol transport, and therefore may also serve as the target for drugs that can act as anthelmintics on their own. Knowledge of polymorphism in these ABC transporter genes in nematode parasites could provide useful information for the process of drug design. We have identified 15 ABC transporter genes from sub-families A, B, C and G, in D. immitis, by comparative genomic approaches and analyzed them for polymorphism. Whole genome sequencing data from four ML susceptible (SUS) and four loss of efficacy (LOE) pooled populations were used for single nucleotide polymorphism (SNP) genotyping. Out of 231 SNPs identified in those 15 ABC transporter genes, 89 and 75 of them were specific to the SUS or LOE populations, respectively. A few of the SNPs identified may affect gene expression, protein function, substrate specificity or resistance development and may be useful for transporter inhibitor/anthelmintic drug design, or in order to anticipate resistance development.

•

In the D. immitis genome, all ABC-A, -B, -C and -G transporter genes were identified.

•

Within 15 ABC transporter genes identified in D. immitis, 231 SNP loci were found.

•

Four exonic SNPs caused changes in predicted secondary structure of ABC proteins.

•

D. immitis populations have low genetic variability among ABC transporter genes.

In silico analysis of the binding of anthelmintics to Caenorhabditis elegansP-glycoprotein 1

Macrocyclic lactones (ML) are important anthelmintics used in animals and humans against parasite nematodes, but their therapeutic success is compromised by the spread of ML resistance. Some ABC transporters, such as p-glycoproteins (Pgps), are selected and overexpressed in ML-resistant nematodes, supporting a role for some drug efflux proteins in ML resistance. However, the role of such proteins in ML transport remains to be clarified at the molecular level. Recently, Caenorhabditis elegans Pgp-1 (Cel-Pgp-1) has been crystallized, and its drug-modulated ATPase function characterized in vitro revealed Cel-Pgp-1 as a multidrug transporter. Using this crystal structure, we have developed an in silico drug docking model in order to study the binding of ML and other anthelmintic drugs to Cel-Pgp-1. All tested ML bound with high affinity in a unique site, within the inner chamber of the protein, supporting that ML may be transported by Cel-Pgp-1. Interestingly, interacting residues delineate a ML specific fingerprint involving H-bonds, including T1028. In particular, benzofurane and spiroketal moieties bound to specific sub-sites. When compared with the aglycone ML, such as moxidectin and ivermectin aglycone, avermectin anthelmintics have significant higher affinity for Cel-Pgp-1, likely due to the sugar substituent(s) that bind to a specific area involving H-bonds at Y771. Triclabendazole, closantel and emodepside bound with good affinities to different sub-sites in the inner chamber, partially overlapping with the ML binding site, suggesting that they could compete for Cel-Pgp-1-mediated ML transport. In conclusion, this work provides novel information on the role of nematode Pgps in transporting anthelmintics, and a valuable tool to predict drug-drug interactions and to rationally design new competitive inhibitors of clinically-relevant nematode Pgps, to improve anthelmintic therapeutics.

Acquired Tolerance to Ivermectin and Moxidectin after Drug Selection Pressure in the Nematode Caenorhabditis elegans

Ivermectin and moxidectin are the most widely administered anthelmintic macrocyclic lactones (MLs) to treat human and animal nematode infections. Their widespread and frequent use has led to a high level of resistance to these drugs. Although they have the same mode of action, differences in terms of selection for drug resistance have been reported. Our objective was to study and compare changes occurring upon ivermectin or moxidectin selection in the model nematode Caenorhabditis elegans C. elegans worms were submitted to stepwise exposure to increasing doses of moxidectin. The sensitivity of moxidectin-selected worms to MLs was determined in a larval development assay and compared with those of wild-type and ivermectin-selected strains. Selection with either ivermectin or moxidectin led to acquired tolerance to ivermectin, moxidectin, and eprinomectin. Importantly, moxidectin was the most potent ML in both ivermectin- and moxidectin-selected strains. Interestingly, this order of potency was also observed in a resistant Haemonchus contortus isolate. In addition, ivermectin- and moxidectin-selected strains displayed constitutive overexpression of several genes involved in xenobiotic metabolism and transport. Moreover, verapamil potentiated sensitivity to ivermectin and moxidectin, demonstrating that ABC transporters play a role in ML sensitivity in ML-selected C. elegans strains. Finally, both ivermectin- and moxidectin-selected strains displayed a dye-filling-defective phenotype. Overall, this work demonstrated that selection with ivermectin or moxidectin led to cross-resistance to several MLs in nematodes and that the induction of detoxification systems and defects in the integrity of amphidial neurons are two mechanisms that appear to affect the responsiveness of worms to both ivermectin and moxidectin.

ABC-B transporter genes in Dirofilaria immitis

Dirofilaria immitis is a filarial nematode causing infection and heartworm disease in dogs and other canids, cats, and occasionally in humans. Prevention with macrocyclic lactones (ML) is recommended during the mosquito transmission season. Recently, ML resistance has been reported. ABC-B transporter genes are thought to be involved in the mechanism of ML resistance in other nematodes. This study aimed to identify all the ABC-B transporter genes in D. immitis using as a reference the nDi.2.2 D. immitis whole genome, which is not completely annotated. Using bioinformatic tools and PCR amplification on pooled D. immitis genomic DNA and on pooled cDNA, nine ABC transporter genes including one pseudogene were characterized. Bioinformatic and phylogenetic analyses allowed identification of three P-glycoproteins (Pgps) (Dim-pgp-3 Dim-pgp-10, Dim-pgp-11), of two ABC-B half transporter genes (one ortholog of Cel-haf-4 and Cel-haf-9; and one ortholog of Cel-haf-1 and Cel-haf-3), of one ABC half transporter gene (ortholog of Cel-haf-5) that contained an ABC-C motif, and of one additional half transporter that would require functional study for characterization. The number of ABC-B transporter genes identified was lower than in Caenorhabditis elegans and Haemonchuscontortus. Further studies are needed to understand their possible role in ML resistance in D. immitis. These ABC transporters constitute a base for ML resistance investigation in D. immitis and advance our understanding of the molecular biology of this parasite.

•

Identification of ABC-B full and half transporter genes in Dirofilaria immitis.

•

Phylogenetic analysis of the D. immitis ABC-B transporter genes.

•

Lower number of ABC-B transporter genes in D. immitis compared with Clade V nematodes.

Effects of in vitro exposure to ivermectin and levamisole on the expression patterns of ABC transporters in Haemonchus contortus larvae

This study investigated the interaction of ATP binding cassette (ABC) transport proteins with ivermectin (IVM) and levamisole (LEV) in larvae of susceptible and resistant isolates of Haemonchus contortus in vitro by measuring transcription patterns following exposure to these anthelmintics. Furthermore, we studied the consequences of drug exposure by measuring the sensitivity of L₃ to subsequent exposure to higher drug concentrations using larval migration assays. The most highly transcribed transporter genes in both susceptible and resistant L₃ were pgp-9.3, abcf-1, mrp-5, abcf-2, pgp-3, and pgp-10. The resistant isolate showed significantly higher transcription of pgp-1, pgp-9.1 and pgp-9.2 compared to the susceptible isolate. Five P-gp genes and the haf-6 gene showed significantly higher transcription (up to 12.6-fold) after 3 h exposure to IVM in the resistant isolate. Similarly, five P-gp genes, haf-6 and abcf-1 were transcribed at significantly higher levels (up to 10.3-fold) following 3 h exposure to LEV in this isolate. On the other hand, there were no significant changes in transcriptional patterns of all transporter genes in the susceptible isolate following 3 and 6 h exposure to IVM or LEV. In contrast to these isolate-specific transcription changes, both isolates showed an increase in R-123 efflux following exposure to the drugs, suggesting that the drugs stimulated activity of existing transporter proteins in both isolates. Exposure of resistant larvae to IVM or LEV resulted, in some instances, in an increase in the proportion of the population able to migrate at the highest IVM concentrations in subsequent migration assays. The significant increase in transcription of some ABC transporter genes following 3 h exposure to both IVM and LEV in the resistant isolate only, suggests that an ability to rapidly upregulate protective pathways in response to drugs may be a component of the resistance displayed by this isolate.

•

We studied interaction of anthelmintics with ABC transporters in Haemonchus contortus.

•

pgp-1, 2, -9.1, -10, and -11 and haf-6 up-regulated after 3 h exposure to ivermectin.

•

Up-regulation occurred in a drug-resistant isolate but not in a -susceptible isolate.

•

A proportion of the drug exposed larvae showed tolerance to subsequent drug treatment.

•

Rapid up-regulation of transporters may be component of resistance in parasitic nematodes.

Investigating the Role of the Host Multidrug Resistance Associated Protein Transporter Family in Burkholderia cepacia Complex Pathogenicity Using a Caenorhabditis elegans Infection Model

This study investigated the relationship between host efflux system of the non-vertebrate nematode Caenorhabditis elegans and Burkholderia cepacia complex (Bcc) strain virulence. This is the first comprehensive effort to profile host-transporters within the context of Bcc infection. With this aim, two different toxicity tests were performed: a slow killing assay that monitors mortality of the host by intestinal colonization and a fast killing assay that assesses production of toxins. A Virulence Ranking scheme was defined, that expressed the toxicity of the Bcc panel members, based on the percentage of surviving worms. According to this ranking the 18 Bcc strains were divided in 4 distinct groups. Only the Cystic Fibrosis isolated strains possessed profound nematode killing ability to accumulate in worms’ intestines. For the transporter analysis a complete set of isogenic nematode single Multidrug Resistance associated Protein (MRP) efflux mutants and a number of efflux inhibitors were interrogated in the host toxicity assays. The Bcc pathogenicity profile of the 7 isogenic C. elegans MRP knock-out strains functionality was classified in two distinct groups. Disabling host transporters enhanced nematode mortality more than 50% in 5 out of 7 mutants when compared to wild type. In particular mrp-2 was the most susceptible phenotype with increased mortality for 13 out 18 Bcc strains, whereas mrp-3 and mrp-4 knock-outs had lower mortality rates, suggesting a different role in toxin-substrate recognition. The use of MRP efflux inhibitors in the assays resulted in substantially increased (>40% on average) mortality of wild-type worms.

Macrocyclic lactones differ in interaction with recombinant P-glycoprotein 9 of the parasitic nematode Cylicocylus elongatus and ketoconazole in a yeast growth assay

Macrocyclic lactones (MLs) are widely used parasiticides against nematodes and arthropods, but resistance is frequently observed in parasitic nematodes of horses and livestock. Reports claiming resistance or decreased susceptibility in human nematodes are increasing. Since no target site directed ML resistance mechanisms have been identified, non-specific mechanisms were frequently implicated in ML resistance, including P-glycoproteins (Pgps, designated ABCB1 in vertebrates). Nematode genomes encode many different Pgps (e.g. 10 in the sheep parasite Haemonchus contortus). ML transport was shown for mammalian Pgps, Pgps on nematode egg shells, and very recently for Pgp-2 of H. contortus. Here, Pgp-9 from the equine parasite Cylicocyclus elongatus (Cyathostominae) was expressed in a Saccharomyces cerevisiae strain lacking seven endogenous efflux transporters. Pgp was detected on these yeasts by flow cytometry and chemiluminescence using the monoclonal antibody UIC2, which is specific for the active Pgp conformation. In a growth assay, Pgp-9 increased resistance to the fungicides ketoconazole, actinomycin D, valinomycin and daunorubicin, but not to the anthelmintic fungicide thiabendazole. Since no fungicidal activity has been described for MLs, their interaction with Pgp-9 was investigated in an assay involving two drugs: Yeasts were incubated with the highest ketoconazole concentration not affecting growth plus increasing concentrations of MLs to determine competition between or modulation of transport of both drugs. Already equimolar concentrations of ivermectin and eprinomectin inhibited growth, and at fourfold higher ML concentrations growth was virtually abolished. Selamectin and doramectin did not increase susceptibility to ketoconazole at all, although doramectin has been shown previously to strongly interact with human and canine Pgp. An intermediate interaction was observed for moxidectin. This was substantiated by increased binding of UIC2 antibodies in the presence of ivermectin, moxidectin, daunorubicin and ketoconazole but not selamectin. These results demonstrate direct effects of MLs on a recombinant nematode Pgp in an ML-specific manner.

Macrocyclic lactones (MLs) are widely used drugs against parasitic nematodes, but drug resistance is rapidly increasing in prevalence and spatial distribution in parasites of ruminants and horses, and is suspected in human nematodes after mass drug applications. Changes in expression levels or the amino acid sequences of P-glycoprotein (Pgp) transporters have frequently been implicated in ML resistance, but direct evidence for transport of MLs by nematode Pgps is still missing. Here, cloning of pgp-9 of the equine parasite Cylicocyclus elongatus and its functional recombinant expression in a Saccharomyces cerevisiae yeast strain deficient in seven endogenous ABC transporters is described. Expression decreased susceptibility to several fungicidal mammalian Pgp substrates including e.g. actinomycin D and ketoconazole, but had no influence on susceptibility to the benzimidazole thiabendazole, which is active against both, yeasts and nematodes. Addition of some MLs strongly increased ketoconazole susceptibility in yeasts expressing C. elongatus Pgp-9, while other MLs had no effect. These interactions are a strong hint that some MLs act as substrates or at least as inhibitors of Pgp-9 mediated drug transport.

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

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

Increased susceptibility of a triclabendazole (TCBZ)-resistant isolate of Fasciola hepatica to TCBZ following co-incubation in vitro with the P-glycoprotein inhibitor, R(+)-verapamil

A study was carried out to investigate whether the action of triclabendazole sulphoxide (TCBZ.SO) against the liver fluke, Fasciola hepatica is altered by inhibition of P-glycoprotein (Pgp)-linked drug efflux pumps. The Oberon TCBZ-resistant and Cullompton TCBZ-susceptible fluke isolates were used for this in vitro study and the Pgp inhibitor selected was R(+)-verapamil [R(+)-VPL]. For experiments with the Oberon isolate, flukes were incubated for 24 h with either R(+)-VPL (1×10-4 m) on its own, TCBZ.SO (15 μg mL-1) alone, a combination of R(+)-VPL (1×10-4 m) plus TCBZ.SO (15 μg mL-1), TCBZ.SO (50 μg mL-1) on its own, or a combination of TCBZ.SO (50 μg mL-1) plus R(+)-VPL (1×10-4 m). They were also incubated in TCBZ.SO (50 μg mL-1) alone or in combination with R(+)-VPL (1×10-4 m) until they became inactive; and in TCBZ.SO (50 μg mL-1) alone for a time to match that of the combination inactivity time. Flukes from the Cullompton isolate were treated with either TCBZ.SO (50 μg mL-1) alone or in combination with R(+)-VPL (1×10-4 m) until they became inactive, or with TCBZ.SO (50 μg mL-1) alone time-matched to the combination inactivity time. Morphological changes resulting from drug treatment and following Pgp inhibition were assessed by means of scanning electron microscopy. Incubation in R(+)-VPL alone had a minimal effect on either isolate. TCBZ.SO treatment had a relatively greater impact on the TCBZ-susceptible Cullompton isolate. When R(+)-VPL was combined with TCBZ.SO in the incubation medium, however, the surface disruption to both isolates was more severe than that seen after TCBZ.SO treatment alone; also, the time taken to reach inactivity was shorter. More significantly, though, the potentiation of drug activity was greater in the Oberon isolate; also, it was more distinct at the higher concentration of TCBZ.SO. So, the Oberon isolate appears to be particularly sensitive to efflux pump inhibition. The results of this study suggest that enhanced drug efflux in the Oberon isolate may be involved in the mechanism of resistance to TCBZ.

NeuroChip: a microfluidic electrophysiological device for genetic and chemical biology screening of Caenorhabditis elegans adult and larvae

Genetic and chemical biology screens of C. elegans have been of enormous benefit in providing fundamental insight into neural function and neuroactive drugs. Recently the exploitation of microfluidic devices has added greater power to this experimental approach providing more discrete and higher throughput phenotypic analysis of neural systems. Here we make a significant addition to this repertoire through the design of a semi-automated microfluidic device, NeuroChip, which has been optimised for selecting worms based on the electrophysiological features of the pharyngeal neural network. We demonstrate this device has the capability to sort mutant from wild-type worms based on high definition extracellular electrophysiological recordings. NeuroChip resolves discrete differences in excitatory, inhibitory and neuromodulatory components of the neural network from individual animals. Worms may be fed into the device consecutively from a reservoir and recovered unharmed. It combines microfluidics with integrated electrode recording for sequential trapping, restraining, recording, releasing and recovering of C. elegans. Thus mutant worms may be selected, recovered and propagated enabling mutagenesis screens based on an electrophysiological phenotype. Drugs may be rapidly applied during the recording thus permitting compound screening. For toxicology, this analysis can provide a precise description of sub-lethal effects on neural function. The chamber has been modified to accommodate L2 larval stages showing applicability for small size nematodes including parasitic species which otherwise are not tractable to this experimental approach. We also combine NeuroChip with optogenetics for targeted interrogation of the function of the neural circuit. NeuroChip thus adds a new tool for exploitation of C. elegans and has applications in neurogenetics, drug discovery and neurotoxicology.

Genetic variants and increased expression of Parascaris equorum P-glycoprotein-11 in populations with decreased ivermectin susceptibility

Macrocyclic lactones (MLs) represent the major drug class for control of parasitic infections in humans and animals. However, recently reports of treatment failures became more frequent. In addition to human and ruminant parasitic nematodes this also is the case for the horse-nematode Parascaris equorum. Nevertheless, to date the molecular basis of ML resistance is still not understood. Unspecific resistance mechanisms involving transporters such as P-glycoproteins (Pgps) are expected to contribute to ML resistance in nematodes. Here, complete sequences of two P. equorum Pgps were cloned and identified as orthologs of Caenorhabditis elegans Ppg-11 and an unnamed Caenorhabditis briggsae Pgp designated as Pgp-16 using phylogenetic analysis. Quantitative real-time PCR was used to compare expression between tissues. Significantly higher PeqPgp-11 expression was found in the gut for both genders, whereas for PeqPgp-16 the body wall was identified as predominant expression site. Furthermore, Pgps were analyzed regarding their participation in resistance development. Using SeqDoC analyses, Pgp-sequences of P. equorum populations with different ML susceptibility were compared. This approach revealed three single nucleotide polymorphisms (SNPs) causing missense mutations in the PeqPgp-11 sequence which correlated with decreased ML susceptibility. However, no resistance associated differences in mRNA expression levels were detected between embryonated eggs of these populations. In contrast, comparison of two pre-adult groups with different ivermectin (IVM) susceptibility revealed the presence of the three SNPs and in addition statistically significant PeqPgp-11 overexpression in the group of worms with reduced susceptibility. These results indicate that Pgp-11 might be involved in IVM resistance in P. equorum as it shows increased expression in an IVM exposed life-cycle stage of an IVM resistant population as well as occurrence of putatively resistance associated SNPs in populations with reduced IVM susceptibility. These SNPs are promising diagnostic candidates for detection of ML resistance with potential also for other parasitic nematode species.

Gene expression analysis of ABC transporters in a resistant Cooperia oncophora isolate following in vivo and in vitro exposure to macrocyclic lactones

Members of the ATP-binding cassette (ABC) transporter family (P-glycoproteins, Half-transporters and Multidrug Resistant Proteins) potentially play a role in the development of anthelmintic resistance. The aim of this study was to investigate the possible involvement of ABC transporters in anthelmintic resistance in the bovine parasite, Cooperia oncophora. Partial sequences of 15 members of the ABC transporter protein family were identified, by mining a transcriptome dataset combined with a degenerate PCR approach. Reverse transcriptase PCR showed that most of the ABC transporters identified were constitutively transcribed throughout the life cycle of C. oncophora. Constitutive differences in gene transcript levels between a susceptible and resistant isolate were only observed for Con-haf-9 and Con-mrp-1 in eggs of the resistant isolate, while no differences were observed in L3 or the adult life stage. Analysis of resistant adult worms, collected from calves 14 days after treatment with either ivermectin or moxidectin, showed a significant 3- to 5-fold increase in the transcript levels of Con-pgp-11 compared to non-exposed worms. Interestingly, a 4-fold transcriptional up-regulation of Con-pgp-11 was also observed in L3 of the resistant isolate, after in vitro exposure to different concentrations of ivermectin, whereas this effect was not observed in exposed L3 of the susceptible isolate. The results suggest that the worms of this particular resistant isolate have acquired the ability to up-regulate Con-pgp-11 upon exposure to macrocyclic lactones. Further work is needed to understand the genetic basis underpinning this process and the functional role of PGP-11.

Caenorhabditis elegans: modest increase of susceptibility to ivermectin in individual P-glycoprotein loss-of-function strains

P-glycoproteins (Pgps) are members of the ABC transporter superfamily and are involved in detoxification mechanisms of single- and multicellular organisms. Their importance for survival of organisms in the presence of harmful drug concentrations has been widely studied in cancer cells but Pgp-dependent drug resistance of parasites has also been demonstrated. Ivermectin (IVM), a widely used anthelmintic in human and veterinary medicine, is a known substrate at least of mammalian Pgps and resistance against IVM is proposed to be associated with Pgps. The consequences of loss of Pgp function for the development of the model nematode Caenorhabditis elegans were analysed in the presence of IVM. Either strains missing only a single Pgp were used or Pgp activity generally was inhibited using verapamil (VPL). Loss-of-function of individual Pgp resulted in a statistically significant increase in IVM susceptibility in terms of impaired development with decreases in EC₅₀ values between 1.5- and 4.3-fold. Absence of seven Pgps resulted in a higher impact on IVM susceptibility of C. elegans since it resulted in EC₅₀ values decreased by 2.4- to 4.3-fold. This increase in IVM susceptibility was even more pronounced than that observed when Pgp function was blocked in general by VPL (approximately 2.5-fold). This study demonstrates clearly that Pgps are of importance for IVM detoxification in the model organism C. elegans and that some Pgps obviously have a higher impact than others.