microRNAs: a role in drug resistance in parasitic nematodes?

Drug resistance in parasitic nematodes is an increasing problem worldwide, with resistance reported to all three commonly used classes of anthelmintics. Most studies to date have sought to correlate the resistant phenotype with genotypic changes in putative target molecules. Although this approach has identified mutations in several relevant genes, resistance might result from a complex interaction of different factors. Here we propose an alternative mechanism underlying the development of drug resistance based on functional differences in microRNA activity in resistant parasites. microRNAs play an important role in resistance to chemotherapeutic agents in many tumour cells and here we discuss whether they might also be involved in anthelmintic resistance in parasitic nematodes. Copyright 2010 Elsevier Ltd. All rights reserved.

The problem of drug resistance in parasitic nematodes

Parasitic nematodes of humans and animals remain a major problem in many areas of the world. Control of these parasites relies almost exclusively on the use of chemotherapeutic agents, of which three classes are in common use. These are the tetrahydropyrimidines/imidazothiazoles (e.g. pyrantel, levamisole), the benzimidazoles (e.g. fenbendazole) and the macrocyclic lactones (e.g. ivermectin, moxidectin). Whereas most drugs retain efficacy against human parasitic nematodes, there are indications that resistance to ivermectin might be developing in the filarial nematode Onchocerca volvulus [1]. By contrast, the problem of drug resistance in veterinary nematodes is widespread, particularly in parasites of sheep, goats and horses [2] and, more recently, cattle [3]. For sheep parasites such as Teladorsagia circumcincta and Haemonchus contortus, resistance has been reported to all three classes of drug and, indeed, triple-resistant parasites exist in parts of the world, including Scotland [4]. As well as causing considerable economic loss (an estimated £84 million per annum to the sheep industry in the UK [5]), these pathogens cause significant welfare problems in infected animals. In the short term, the availability of novel anthelmintics, such as the amino-acetonitrile derivative (AAD) [6], monepantel, marketed as Zolvix, offer temporary relief. Yet, by analogy with how quickly resistance has developed to the current compounds, it is probably simply a matter of time before the efficacy of novel compounds is compromised.

Understanding the mode of action of anthelmintics and the mechanisms of resistance is important in terms of attempting to monitor and modulate the resistant phenotype. Much of our knowledge of mechanisms of drug action and the likely genes involved in resistance comes from studies on Caenorhabditis elegans [7], where the relative ease of generating resistant mutants and identifying mutations has facilitated such studies. This has been most clearly demonstrated for the β-tubulin gene ben-1, in which mutations conferred resistance to benzimidazole (BZ) compounds in C. elegans [8]. In H. contortus, mutations in the β-tubulin isotype 1 gene (Hc-iso-1) are also associated with BZ resistance. Transformation of BZ-resistant C. elegans with the wild-type Hc-iso-1 allele restores sensitivity to BZ compounds, whereas transformation with the mutant Hc-iso-1 gene had no such effect [9]. These results helped establish that β-tubulin isotype-1 mutations are a major determinant of BZ resistance in parasitic nematodes. Confirming a role for specific gene mutations in resistance to other anthelmintics has, however, proven more difficult. For example, studies in C. elegans have indicated that resistance to ivermectin might be polygenic and that mutations at any one contributing locus may be insufficient for a fully resistant phenotype [10]. It has been suggested that a similar scenario might exist in parasitic nematodes [11]. The situation is confounded by the significant genetic polymorphism of parasite populations, particularly for trichostrongylid species [12], making it difficult to conclude that resistance is due to a particular polymorphism, or whether other linked mutations might be responsible.

Recent studies on possible mechanisms of drug resistance in parasitic nematodes have identified mutations within the coding regions or the promoters of candidate genes first identified in C. elegans ([11,13-15]). Nevertheless, by analogy with other eukaryotes, drug resistance in nematodes could result from a range of different mechanisms in addition to functional mutations in specific genes. These include altered levels of expression of target or non-target genes, such as drug transporters or detoxification enzymes. Changes in the expression of such genes could arise via a number of different processes such as epigenetic modifications in chromatin or by a variety of transcriptional or post-transcriptional mechanisms, including trans-splicing and capping [16]. In tumour cells, one specific type of post-transcriptional regulation, the alteration of microRNA (miRNA) expression or activity, is increasingly implicated in drug resistance. In this article we discuss whether miRNAs are involved in drug resistance in parasitic nematodes.

miRNAs

miRNAs are one of a family of endogenous small non-coding RNAs, including small interfering RNAs (siRNAs) and piwi-interacting RNAs (piRNAs), that have been identified in a diverse range of organisms [17]. They were first identified in C. elegans as key components of the heterochronic pathway, the mechanism by which developmental timing is regulated in the worm [18,19]. miRNAs are the focus of intense interest given their important roles in the control of gene expression in many biological and pathological processes, such as cell and organ development, differentiation and homeostasis, tumour suppression and stem cell regulation [20]. The biogenesis of miRNAs (Box 1) and the mechanisms by which they exert their regulatory function are increasingly understood [21]. miRNAs regulate the expression of target genes at the post-transcriptional level by base pairing with defined sites often, but not exclusively, located in the 3′-UTR of target genes. In animal cells, specific miRNAs have been shown to repress translation, whereas others appear to act by inducing degradation of the target mRNA. Since their discovery in C. elegans, miRNAs have been identified in a diverse range of plants and animals and appear to be abundant in most genomes; the current version (release 15) of miRBase (http://www.mirbase.org/) identifies 940 miRNAs in humans and 175 in C. elegans. In addition, up to one third of all human mRNAs might be regulated by miRNAs [22]. Now that genome sequence data are available for several species of parasitic nematode [23], it is possible to identify miRNAs in parasitic species and to investigate their roles in the biology of parasites. For example, a number of important miRNAs including let-7 and lin-4 have been identified in the B. malayi genome [24].

miRNAs and drug resistance in tumour cells

miRNAs have been the focus of much interest in tumour cell biology because they appear to play a role in the initiation and progress of some cancers and indeed are considered by some to function as oncogenes or tumour suppressors [25]. In addition, many tumour cells express an altered profile of miRNAs, which might have diagnostic and/or prognostic potential [26]. Furthermore, expression profiling of drug-sensitive and drug-resistant tumour cell lines has suggested a role for miRNAs in the development of drug resistance [27]. For example, in the breast cancer cell line MCF-7, doxyrubicin resistance appears to be associated with significant changes in the levels of specific miRNAs [28]. These cells express very high levels of P-glycoprotein (encoded by the mdr-1 gene), an important mechanism of drug efflux. The 3′-UTR of the mdr-1 gene contains a binding site for miR-451, which negatively regulates mdr-1 expression (Figure 1(a)). It was shown that transfection of doxyrubicin-resistant cells with this miRNA resulted in decreased levels of P-glycoprotein and a 2.5-fold increase in sensitivity to doxyrubicin.

Figure 1

Potential mechanisms for miRNA involvement in drug resistance. (a) In human MCF-7 breast cancer cells, miR-451 negatively regulates translation of the P-glycoprotein gene mdr-1, via a complementary miR-451 site in the 3′-UTR of mdr-1. In MCF-7 cells resistant to doxorubicin (MCF-7/DOX), decreased levels of miR-451 result in an increased level of MDR-1 protein with resultant DOX resistance [28]. (b) In C. elegans, miR-1 negatively regulates the unc-29 nAChR subunit via three complementary sites in the unc-29 3′UTR. In mir-1 null mutants, UNC-29 protein level is increased, which alters nAChR composition, leading to a decreased sensitivity to levamisole [37].

A number of recent studies have further demonstrated that drug metabolising enzymes and transporters can be post-transcriptionally regulated by miRNAs and, importantly, that changes in expression of many miRNAs are induced by drug treatment [29]. For example, miR-328 is downregulated ten-fold in xenograft tumours following gemcitabine treatment [30]. miR-328 is involved in the negative regulation of the ABCG2 efflux transporter [31] and therefore downregulation of miR-328 results in overexpression of the transporter and multidrug resistance [30].

As well as altered expression of miRNAs leading to drug resistance, deletions or mutations within the binding targets of miRNAs can affect miRNA function. A single nucleotide polymorphism (SNP) near the binding site of miR-24 in the 3′-UTR of the human dihyrofolate reductase (DHFR) gene prevents repression by miR-24, resulting in overexpression of DHFR and resistance to methotrexate [32]. Indeed mutations within 3′-UTR binding sites occur at a higher level than SNPs in miRNA sequences themselves [33], probably because miRNAs have multiple targets, and therefore miRNA mutation is likely to lead to more serious consequences. In a different system, several drug resistant cell lines that overexpress ABCG2 were shown to contain truncated 3′-UTRs. The deleted region contains a binding site for miR-519c, and the absence of this binding site in resistant cells leads to overexpression of ABCG2 and drug resistance [34]. The accumulating data on miRNA and mRNA expression and drug resistance in cancer cells are now being integrated to predict responses to drugs and improve treatment efficiency, an area referred to as ‘miRNA pharmacogenomics’ [35].

It is possible that a link between miRNAs and drug resistance is more likely to occur in cancer cells, as the cancerous state can itself be due to miRNA mis-expression. However, the study of miRNAs in cancer has clearly established that miRNA-mediated alterations in levels of drug targets, drug transporters, metabolic enzymes or cell apoptosis proteins can lead to drug resistance. Altered gene expression associated with drug resistance in other systems indicates that examination of miRNA activity and 3′-UTR interactions in parasitic nematodes is warranted to improve our understanding of drug resistance mechanisms.

Could miRNAs have a role in drug resistance in nematodes?

The study of miRNAs and their functions in parasitic nematodes is in its infancy, although much is known about specific miRNAs in C. elegans [36]. Recent work on this nematode indicates that miRNAs might play a role in drug resistance, and indicates that this area is ripe for investigation in parasitic nematodes. In C. elegans, miR-1 negatively regulates the expression of two nicotinic acetylcholine receptor (nAChR) subunits, unc-29 and unc-63 [37]. Intriguingly, in mir-1 mutants the expression of both UNC-29 and UNC-63 subunits is increased, and this corresponds with decreased muscle sensitivity to acetylcholine and levamisole (Figure 1(b)). It is perhaps surprising that an increased level of particular AChR subunits leads to an altered response to levamisole. Nevertheless, it was speculated that altering the subunit composition of the AChR could affect receptor biogenesis or function and that miRNAs can alter the composition and activity of other heteromultimeric receptors [37]. These findings are consistent with changes in kinetics and affinity of mammalian nAChRs that follow changes in subunit composition [38] and are supported by studies on the Ascaris suum nAChR which highlighted the importance of subunit composition for drug sensitivity [39].

Interestingly, in a very different approach using immunoprecipitation to isolate miRNA-target complexes in C. elegans during development, mRNAs encoding ion channels and receptors were particularly enriched in the data set [36]. Since these types of proteins are well represented among known drug targets in nematodes, miRNAs might well be involved in altering their expression and activity. Recently, two studies in parasitic nematode species reported alterations in expression levels of receptors and ion channels in drug resistant isolates. Following pyrantel treatment, changes in transcript levels of the nAChR genes unc-29, unc-38 and unc-63 were found in isolates of the dog hookworm Ancylostoma caninum [40]. In this case a highly resistant isolate showed diminished expression of the three receptor subunit genes compared to an isolate with a low level of pyrantel resistance. Whether the lower transcript level of any of these three genes corresponds with an increase in regulatory miRNA expression or activity has not been examined, but would be of great interest. Similarly, isolates of H. contortus resistant to ivermectin were found to have lowered expression of a novel ligand-gated ion channel gene (Hc-GGR3) [41]. Interestingly, a SNP was identified in the Hc-GGR3 3′-UTR of resistant worms. It was proposed that this polymorphism might be associated with resistance, although functional analysis is required to test this. Given the high degree of conservation of some miRNA sequences across species, including the miR-1 sequence, C. elegans presents a suitable model system in which to functionally test potential miRNA-3′-UTR interactions, which cannot readily be examined directly in parasitic species.

As well as changes in expression levels of potential drug targets, alterations in the expression of drug transporters have been reported in parasitic nematodes following anthelmintic treatment. This might be analogous to the changes in transporter levels in tumour cells exposed to drugs, as discussed previously. In H. contortus and O. volvulus, treatment with ivermectin results in a reduction in allele frequency and overexpression of P-glycoprotein and other ABC transporter genes [42]. Similarly, two independent isolates of H. contortus that were selected for BZ resistance show an elevated frequency of the same P-glycoprotein allele [14]. Although different drug classes can interact with different ABC transporters, these findings indicate that the mechanism underlying this allele selection might be the same and have implications for the development of resistance to new drug classes. A recent study in C. elegans also found that exposure to increasing concentrations of ivermectin resulted in increased expression of several P-glycoprotein and ABC transporter genes [43]. It would obviously be interesting to compare the 3′-UTR region of P-glycoprotein genes from susceptible and resistant nematodes for regulation by miRNAs.

Future perspectives

Herein the possibility that changes in gene expression levels associated with drug resistance in parasitic nematodes are mediated by miRNAs has been discussed. This might be a mechanism of altering drug efficacy; given the link between miRNAs and drug responses in cancer cells, this area is ripe for detailed investigation. With the number of genome sequencing projects underway, miRNAs present in parasitic nematodes can be identified using both in silico approaches [44] (for example by combining BLAST [45] and RNAfold [46] analysis to search for homologues of known miRNAs), and small RNA sequencing [24]. Potential target mRNAs can be identified using available bioinformatics programmes such as MiRanda [47] and PicTar [48], combined with immunoprecipitation of mRNAs in RISC complexes and sequencing [49]. A number of potential target mRNAs relevant to drug resistance in parasitic nematodes have been discussed above; bioinformatic and biochemical approaches can be used to determine whether these are likely to be real targets of miRNA regulation. A potential interaction between a miRNA and the 3′-UTR of a target gene from a parasitic nematode could then be verified using C. elegans as a heterologous expression system.

Currently, the potential correlation between antihelmintic resistance and changes in particular gene sequence or expression level has been limited to the study of a small number of candidate genes. Examining and comparing miRNA levels, as well as mRNA expression in drug resistant and susceptible isolates, through the use of microarrays will produce a global picture of changes that might correlate with drug resistance. In addition, miRNA inhibitors and mimics (Box 2) are available and can be used to examine the upregulation or downregulation of miRNAs and their potential mRNA targets. Could these miRNA inhibitors or mimics be effective in reversing drug resistance in nematodes?