http://orcid.org/0000-0002-4713-575X
Figures
Abstract
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.
Citation: Figueiredo LA, Rebouças TF, Ferreira SR, Rodrigues-Luiz GF, Miranda RC, Araujo RN, et al. (2018) Dominance of P-glycoprotein 12 in phenotypic resistance conversion against ivermectin in Caenorhabditis elegans. PLoS ONE 13(2): e0192995. https://doi.org/10.1371/journal.pone.0192995
Editor: Rui Lu, Louisiana State University, UNITED STATES
Received: May 30, 2017; Accepted: February 3, 2018; Published: February 23, 2018
Copyright: © 2018 Figueiredo et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All relevant data are within the paper.
Funding: LAF received a Master Degree Fellowship from FAPEMIG. TFR received fellowships from CAPES and FAPEMIG. RTF is a Research Fellow from CNPq (Bolsa de Produtividade em Pesquisa).
Competing interests: The authors have declared that no competing interests exist.
Introduction
Diseases caused by nematodes remains a considerable drawback for the livestock, agriculture and public health, contributing to increasing production costs, decreased food supply and deficits in host health [1–2]. Programs for parasite control are culturally performed through employment of anthelmintics and the improper use of these drugs for treatment of human and animal populations, like underdosing, contributes to selection of resistant parasites [3–5]. Ivermectin is an important macrocyclic lactone (ML) and the most commonly used drug to treat parasitic nematode infections of the gastrointestinal tract, the lungs, filarial infections and infestations with arthropods [6]. The ubiquitous use of this drug has led to the selection of resistance in several target species including Cooperia oncophora [7], Haemonchus contortus [8], Ostertagia spp. [7], Trichostrongylus spp. [8], Onchocerca volvulus [4, 9–10] and Dirofilaria immitis [11–12].
Previous studies have shown that ML anthelmintics can be transported by several ATP-binding cassette transporters (ABC), whose expression levels might be regulated in ML resistant nematodes [13–15]. Among these ABC transporters, the P-glycoproteins (Pgps) are known to perform cell efflux of several foreign substances [13–16] and their overexpression are related to multidrug resistance in mammalian cells and human tumor cells [17–18]. Although the relationship between Pgps and ivermectin resistance in nematodes has been shown [13, 14, 16, 19–20], here we present some evidences of Pgps dominance 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.
Currently it is known that Caenorhabditis elegans genome encodes 15 Pgps [21]. Due to the simplicity of maintenance, innocuous character and the simple life cycle that last less than 4 days [22], the C. elegans display an excellent and feasible model compared to many other nematodes, providing the opportunity to conduct different assays focusing on viability, behavior, survival, reproduction and development influenced by anthelmintic [23–25]. The aims of this study were focused to further understand the role of Pgps in ivermectin resistance by using an RNAi loss-of-function approach and determine the impact of our findings by use of in silico analyses to evaluate similarity among nematodes Pgps protein sequences.
Materials an methods
Strains and selection of ivermectin-resistant C. elegans
The wild-type Bristol N2 strain of C. elegans was kindly provided by Dr. Carlos Eduardo Winter (University of São Paulo), and cultured on nematode growth medium (NGM) with NA22 Escherichia coli under standard conditions[22]. Resistant worms were stepwise in vitro selection by addition of increased doses (1 ng/ml) of ivermectin (Ourofino, Brazil) to a Petri dish plate. After 1 week at 1 ng/ml ivermectin, all the survivor worms were transferred onto plates containing 2 ng/ml. Thereafter, C. elegans were transferred onto NGM plates containing higher concentrations of drug. After 48 weeks of culture, resistant C. elegans were able to grow at 30 ng/mL ivermectin (IVR30), a concentration that is lethal to the wild type.
C. elegans L3 production
IVR30 and wild-type (WT) C. elegans were cultured over seven days in BOD incubator at 22°C, washed with M9 medium [26] and filtered through three consecutive sieves with pores of 40 μm, 30 μm, and 20 μm, respectively. L3 larvae retained in the 20 μm strainer were collected by backwashing and were washed by centrifugation at 700 g for 4 minutes, followed by two washes with M9 medium.
Pgp gene expression in ivermectin-resistant strain
Worms were harvested from NGM plates and washed using M9 buffer [26]. Total RNA from worm pellets (IVR30 and WT) were extracted using RNA NucleoSpin® RNA II (Macherey-Nagel, Germany) and converted into cDNA by using a Transcriptor High Fidelity cDNA Synthesis Kit (Invitrogen, USA). Quantitative PCR was performed using 50ng of cDNA, 1 μM of gene specific primers and 2X Power Sybr®-Green ABI (Thermo Fisher, USA) in 10μL final volume. Reactions were performed in the Applied Biosystems 7500 Real Time PCR System (Life Technologies, USA) using the following PCR cycling conditions: 95°C 10 min; 95°C 15s and 60°C 1min, for 40 cycles. Analyses were performed at 7500 Software v2.0.1 and results were determined by using 2-ΔΔCt method [27]. Specific Pgps and housekeeping primers were designed using a web-based tool (Primique). Primer sequences as well as sequence confirmation PCR products are presented in Table 1.
RNAi production and assay
Fragments designated for RNAi of Pgps 4 and 12 were obtained by polymerase chain reaction (PCR) from cDNA and were cloned into pGEM®-T Easy vector (Promega, USA). The resulting plasmids were transformed into DH5α E. coli strain. The positive DNA clones were extracted with Wizard® Plus SV Minipreps DNA Purification System (Promega, USA) and digested with NOT I enzyme restriction (Promega, USA). The digested amplicons were cloned into L4440 feeding vector (Addgene, USA). The resulting plasmids were transformed into HT115(DE3) RNase III-deficient E. coli strain. Single colonies of HT115 bacteria containing the cloned L4440 plasmid were picked and grown in LB culture with 100 μg/mL ampicillin and 12.5 μg/mL tetracycline for approximately 2 hours at 37°C and 150 rpm until reaching OD600 (optical density) between 0.5–0.6. Further, 1 nM IPTG was added and bacteria was seed into NGM-agar plates previously supplemented with 100 μg/mL ampicillin, 12.5 μg/mL tetracycline and 0.4mM IPTG. Seeded plates were dried at room temperature and IVR30 and WT C. elegans L3 were transferred onto plates. All plates were kept in BOD incubator at 22°C for 7 days until C. elegans have consumed all bacteria. Additional control using transformed E. coli strain cultured with regular (non-supplemented with IPTG) NGM agar were used in the experiment (control-non silenced).
Motility test
RNAi silenced (IVR30 and WT) and control-non silenced (IVR30 and WT) C. elegans L3 were resuspended in M9, and 1000 larvae in 100 μL of suspension were distributed on 96-well microplate (1000 larvae/well). Ivermectin was added at a final concentration of 30 ng/mL (34mM). Microplates containing ivermectin and larvae were incubated in BOD incubator at 22°C. After 48 hours, 15 μL of culture suspension containing 150 larvae was removed for further analysis and quantification of the number of paralyzed larvae using an optical microscope at 100x magnification. RNAi silenced and control-non silenced worms were cultured only in M9 medium, for negative control group (demonstrating the natural paralyses observed in a culture), and cultured in M9 medium with 1000mM of ivermectin for positive control (where worms were always 100% paralyzed). Larvae were considered paralyzed when presenting with straight body and absence of any motility [28–29].
Propidium iodide staining
After 48 hours of incubation with ivermectin (30 ng/mL), propidium iodide (Sigma-Aldrich, USA) was added to the microplates at a final concentration of 20 μM [29–30]. Again, RNAi silenced and control-non silenced worms were cultured only in M9 medium, for negative control group, and cultured in M9 medium with 1000mM of ivermectin for positive control (The controls fluorescence emission were used to calibrate the software analysis). Microplates were incubated for 15 minutes at room temperature in a horizontal shaker at 120 rpm followed by reading at LAS ImageQuant™ GE 4000 with excitation in white light and emission at 605 nm. Densitometric analyses of the images were performed using the software GE Image Quant TL 8.1.
P-glycoproteins diversity
For P-glycoprotein sequences, a search for protein nematodes available in NCBI database (National Center for Biotechnology Information) was performed using the following keywords: P-glycoprotein, PGP, multidrug resistance protein. Hypothetical and partial protein sequences were excluded. The sequences were aligned using ClustalW 2.0 software with the default parameters. To identify the clusters formed by P-glycoprotein sequences, the pairwise distance was calculated and the distance matrixes were further generated. The distances between the sequences were generated using the package PHYLIP. To provide a visual representation of each distance matrix, we used the multidimensional scaling (MDS) plot with two dimensions. The MDS and graphing were performed using the R software platform.
Statistical analyses
The statistical analysis and graphs were made using GraphPad Prism 7.0 software (GraphPad Inc, USA) using one-way ANOVA test with Turkey’s multiple comparisons test to compare groups and One sample t test to compare the sample value to the reference number (WT basal expression). The significance level was p <0.05.
Results
Expression pattern of Pgps in IVR30 strain
Using resistant and WT strains of C. elegans we evaluated the quantitative gene expression of 12 out the 15 known Pgps. Only the Pgps 12 and 13 were upregulated in the IVR30 in comparison with WT strain. All the other Pgps were downregulated (Fig 1).
Pgps gene expression for WT strains are represented by dashed line. P values were obtained with One sample t test, comparing sample value to WT expression reference value.
To assess the role of Pgps during drug resistance/susceptibility to ivermectin, the expression of all Pgps were knocked down by RNAi. However, consistent results were obtained only for Pgp-4 and Pgp-12.
RNAi-mediated Pgp 12 gene silencing reverts the resistance phenotype to ivermectin
RNAi-mediated impairment of Pgp 12 expression in the IVR30 strain induced a significant increase in worm paralysis (Fig 2A) and viability (Fig 2B), when compared to the control-non silenced worms (p = 0,0207 and p = 0,0009, respectively), after ivermectin challenge. Indeed, the paralysis and viability values of Pgp 12 RNAi-silenced worms approached those of WT strain silenced and not silenced. Gene silencing of Pgp 12 did not affect paralysis and viability of WT worms.
(A) Percentage of L3 paralysis after 48 hours of ivermectin exposure for groups Pgp 12 silenced [RNAi+], Pgp 12 non silenced [RNAi-] and negative control [CT], cultured only in M9 medium and not exposed to ivermectin treatment. (B) Assessment L3 viability cell using propidium iodide marker after 48 hours ivermectin exposure in Pgp 12 silenced [RNAi+] and Pgp 12 non silenced [RNAi-] groups. P value refers to one-way ANOVA test with Turkey’s multiple comparisons test.
RNAi-mediated Pgp 4 gene silencing does not alter the resistance phenotype but induces a resistance in wild type strain
Following the challenge with ivermectin, no significant changes in worm paralysis (Fig 3A) and viability (Fig 3B) were observed in Pgp 4 silenced worms in comparison to the control-non silenced worms of the IVR30 strain. Interestingly, in WT strains RNAi-mediated silencing of Pgp 4 and ivermectin challenge promoted a significant decrease in the paralysis and viability of worms, therefore inducing a phenotype similar to that observed in resistant worms. The paralysis and viability between the both silenced groups (WT and IVR30) did not present any significant differences. Similarly, the results observed for Pgp 4 gene silencing in WT was equivalent to observed for the non-silenced resistant strain.
(A) Percentage of L3 paralysis after 48 hours of ivermectin exposure for groups Pgp 4 silenced [RNAi+], Pgp 4 non silenced [RNAi-] and negative control [CT], cultured only in M9 medium and not exposed to ivermectin treatment. (B) Assessment of L3 viability cell using propidium iodide marker after 48 hours ivermectin exposure in Pgp 4 silenced [RNAi+] and Pgp 4 non silenced [RNAi-] groups. P value refers to one-way ANOVA test with Turkey’s multiple comparisons test.
Dual silencing of Pgp 12 and Pgp 4 expression leads to the dominance of phenotype promoted by Pgp 12 silencing
Dual silencing of Pgp 12 and Pgp 4 was performed to investigate whether these share a role in the ivermectin resistant phenotype. Following ivermectin challenge, we found a significant increase in paralysis and viability of worms in the IVR30 strain when compared to the control-non silenced group (p = 0.0395 and p = 0.0293, respectively). Significant differences in the paralysis (p = 0.0190) but not in the cellular viability was detected when the dual silenced groups (WT and IVR30) were compared. It is worth noticing that although Pgp 4 silencing affected the paralysis and viability of WT worms, dual silencing did not produce significant changes after ivermectin challenge (Fig 4).
(A) Percentage of L3 paralysis percentage after 48 hours ivermectin exposure for groups Pgps 4 and 12 silenced [RNAi+], Pgps 4 and 12 non silenced [RNAi-] and negative control [CT], cultured only in M9 medium and not exposed to ivermectin treatment. (B) Assessement of L3 viability cell using propidium iodide marker after 48 hours ivermectin exposure in Pgps 4 and 12 silenced [RNAi+] and Pgps 4 and 12 non silenced [RNAi-] groups. P value refers to one-way ANOVA test with Turkey’s multiple comparisons test.
In silico analysis reveals a close relationship between Pgps from C. elegans and nematodes parasites
A phylogenetic approach was performed using protein sequences of Pgps from different organisms including helminthes parasites with medical relevance to investigate the level of similarity among them (Fig 5). Protein sequence cluster analysis indicate that C. elegans Pgps are highly similar to the Pgps of nematodes parasites such as Ascaris suum (A), Brugia malayi (B), Cooperia oncophora (C), Haemonchus contortus (H), Loa loa (L), Onchocerca volvulus (O), Parascaris equorum (P), Strongyloides ratti (S), Trichuris trichiura (T) and Wuchereria bancrofti (W).
As follow: Ascaris suum (A), Brugia malayi (B), Cooperia oncophora (C), Haemonchus contortus (H), Loa loa (L), Onchocerca volvulus (O), Parascaris equorum (P), Strongyloides ratti (S), Trichuris trichiura (T), Wuchereria bancrofti (W), C. brenneri (1), C. briggsae (2), C. remanei (4), Pgp 4 of C. elegans (Δ), Pgp 12 of C. elegans (•), all remaining Pgps from C. elegans (2).
Discussion
Initially considered as a highly potent drug, ivermectin has now a reduced anti-helminthic efficacy and therefore research on the possible mechanisms of resistance are important, especially in the light of the negative effects of diseases caused by parasitic nematodes in animal [31–32] and in human [33–34] healthcare. Previous studies revealed the involvement of Pgps in the development of ivermectin resistance in various nematodes [16, 19–20], including the free-living nematode, C. elegans [35–37]. We have confirmed some of these results and extended the study on the role of Pgps in ivermectin resistance developed by C. elegans to show that Pgp 4 might be directly related to susceptible phenotype and that Pgp 12 directly related to resistant phenotype in this worm. In addition, we show that Pgps 1 to 13 of C. elegans have a close phylogenetic relationship to Pgps of other parasitic nematodes.
In our results we observed that expression of most Pgps were downregulated when compared with WT whereas Pgps 12 and 13 were upregulated in an ivermectin-resistant C. elegans strain (IVR30). In contrast, Yan et al. [38] found that most of the 15 Pgps they analyzed were upregulated in their ivermectin-resistant strain (IVR10). We hypothesized that the divergence between the results may be due to difference of resistant strains selected, as in the study cited above they obtained a 10 ng/mL ivermectin-resistant strain and we selected a 30 ng/mL ivermectin-resistant C. elegans strain. In both WT and IVR30 we silenced Pgp 4 and Pgp 12 that were characteristically downregulated and upregulated, respectively. The RNAi method has been widely used in C. elegans to study gene function [39]. Interestingly, the results for paralysis and cellular viability observed for the Pgp 12 gene silencing in IVR30 strain were similar to those observed in both silenced or non-silenced WT strains. Therefore, it indicates that RNAi-mediated Pgp 12 gene silencing in IVR30 reverted the resistant phenotype to the WT, suggesting that this Pgp may be directly related to ivermectin resistance phenotype. Two studies using several C. elegans Pgps knockout strains observed that deletion of Pgp 12 resulted in increased sensitivity to ivermectin [35–36], corroborating our results.
The Pgp 4 gene silencing in IVR30 strain did not lead to a phenotypic changes after ivermectin exposure. However, such gene silencing promoted a significant reduction of paralysis and increased cellular viability in the WT strain, which rendered a phenotype observed in the resistant IVR30 strain. Furthermore, the RNAi-mediated Pgp 4 gene silencing in IVR30 did not alter the resistance phenotype but induced resistance in the WT strain, suggesting that Pgp 4 may be directly related to ivermectin susceptibility. In contrast, previous results with C. elegans Pgp 4 knockout strains showed a slight increase in sensitivity to the drug in the WT. However, other authors [38] using a high concentration of ivermectin (20 ng/mL) found that the RNAi-mediated silencing of Pgp 4 led to downregulation of this gene and a concomitant increase in the motility of the WT strain of C. elegans, a finding that is similar to the effect reported in this study.
To date no study have describe the double-silencing of nematode Pgps. In our results, the dual silencing of Pgp 12 and Pgp 4 expression in the IVR30 strain lead to increase of worm paralysis and reduction of cellular viability after exposure to ivermectin. However, the results obtained for the dual Pgp 12 and Pgp 4 silencing were similar to those obtained when Pgp 12 was silenced alone, suggesting that dual silencing of Pgp 12 and Pgp 4 expression in the IVR30 strain led to the dominance of the phenotype promoted by Pgp 12 silencing. Moreover, dual silencing of these genes did not alter the susceptibility status of the WT strain. It is possible that silencing of both Pgps in the WT strain generated a compensation phenotype effect so that the resistant phenotype, generated by Pgp 4 silencing was suppressed by Pgp 12 silencing phenotype.
To assess whether the roles ascribed for the Pgps in C. elegans resistance to ivermectin could be extrapolated to parasite nematodes, one first step is to evaluate similarity of these Pgps within the phylum Nematoda. The spatial projection allows a global view of the closest sequences (closest in spatial projection points) and the most distant (farthest points in the projection), allowing the analysis of the distribution (diversity) and the formation of groups. The distances between the closest and farthest points can be used to infer phylogenetic trees, diversity and divergence between sequences [40]. The protein sequence cluster analysis performed indicates that all C. elegans Pgps studied herein are highly similar to the Pgps of various nematodes parasites (Fig 5).
Our study is the first to perform dual silencing and reveal a compensation mechanism of one Pgp over another to modulate C. elegans resistance status. We are also the first to analyze phylogenetically the Pgps and show that they are highly similar in many of the nematodes of medical and economical relevance.
Conclusion
In conclusion, here we show that Pgp 12 is a good candidate for further studies aiming to develop specific inhibitors to this transporter so that ivermectin can still be used to control the burden on animal and human health inflicted by nematode parasites globally.
References
- 1. Moncayo AL, Cooper PJ. Geohelminth infections: impact on allergic diseases. The international journal of biochemistry & cell biology. 2006; 38: 1031–1035.
- 2.
Nicol JM, Turner SJ, Coyne DL, den Nijs L H S TM Z. Current Nematode Threats to World Agriculture, In: John J, Gheysen G, Fenoll C, Ankara, Genomics and Molecular Genetics of Plant-Nematode Interactions. 2011; Springer.
- 3. Wolstenholme AJ, Fairweather I, Prichard R, Samson-Himmelstjerna G, Sangster NC., 2004. Drug resistance in veterinary helminthes. Trends Parasitol. 2004: Vol.20 No.10.
- 4. Osei-Atweneboana MY, Awadzi K, Attah SK, Boakye DA, Gyapong JO, Prichard RK. Phenotypic evidence of emerging ivermectin resistance in Onchocerca volvulus. PLoS Negl Trop Dis. 2011; pmid:21468315
- 5.
World Health Organization, 2015.
- 6. Geary TG. Ivermectin 20 years on: maturation of a wonder drug. Trends in parasitology. 2005; v. 21, n. 11, p. 530–2.
- 7. Areskog M, Sollenberg S, Engström A, Samson-Himmelstjerna G, Höglund J. A controlled study on gastrointestinal nematodes from two Swedish cattle farms showing field evidence of ivermectin resistance. Parasites & Vectors. 2014; 7:13.
- 8. Holsback L, da Silva MA, Patelli THC,de Jesus AP, Sanches JR. Resistance of Haemonchus, Cooperia, Trichostrongylus, and Oesophagostomum to ivermectin in dairy cattle in Paraná. Semina: Ciências Agrárias, Londrina. 2015; v. 36, n. 3, suplemento 1, p. 2031–2036.
- 9. Awadzi K, Boakye DA, Edwards G, Opoku NO, Attah SK, Osei-Atweneboana MY, et al. An investigation of persistent microfilaridermias despite multiple treatments with ivermectin, in two onchocerciasis-endemic foci in Ghana. Ann Trop Med Parasitol. 2004a; 98: 231–149.
- 10. Awadzi K, Attah SK, Addy ET, Opoku NO, Quartey BT, Lazdins-Helds JK, et al. Thirty month follow-up of suboptimal responders to multiple treatments with ivermectin, in two onchocerciasis-endemic foci in Ghanda. Ann Trop Med Parasitol. 2004b; 98: 359.
- 11. Bourguinat C, Keller K, Blagburn B, Schenker R, Geary TG, Prichard RK. Correlation between loss of efficacy of macrocyclic lactone heartworm anthelmintics and P-glycoprotein genotype. Veterinary Parasitol. 2011;176:374–381.
- 12. Pulaski CN, Malone JB, Bourguinat C, Prichard R, Geary T, Ward D, et al. Establishment of macrocyclic lactone resistant Dirofilaria immitis isolates in experimentally infected laboratory dogs. Parasites Vectors. 2014; 7:494. pmid:25376278
- 13. Xu M, Molento M, Blackhall WJ, Ribeiro P, Beech RN, Prichard RK. Ivermectin resistance in nematodes may be caused by alteration of P-glycoprotein homolog. Mol. Biochem. Parasitol. 1998; 91, 327–335. pmid:9566525
- 14. Kerboeuf D, Guégnard F, Le Vern Y.Detection of P-glycoprotein-mediated multidrug resistance against anthelmintics in Haemonchus contortus using anti-human mdr1 monoclonal antibodies. Parasitology Research. 2003; Volume 91, Issue 1, pp 79–85.
- 15. Jones PM, George AM. The ABC transporter structure and mechanism: perspectives on recent research. Cell Mol Life Sci. 2004; 61:682–699. pmid:15052411
- 16. Blackhall W,Liu HY, Xu M, Prichard RK, Beech RN. Selection at a P-glycoprotein gene in ivermectin- and moxidectin-selected strains of Haemonchus contortus. Mol Biochem Parasitol. 1998; Sep 15;95(2):193–201 pmid:9803412
- 17. Sharom FJ. The P-glycoprotein multidrug transporter. Biochemical Society Essays Biochem. 2011; 50, 161–178.
- 18. Lespine A, Menez C, Bourguinat C, Prichard RK. P-glycoproteins and other multidrug resistance transporters in the pharmacology of anthelmintics: prospects for reversing transport-dependent anthelmintic resistance. Int. J. Parasitol. Drugs Drug Resist. 2012; 2:58–75. pmid:24533264
- 19. Bartley DJ. McAllister H, Bartley Y, Dupuy J, Ménez C, Alvinerie M. P-glycoprotein interfering agents potentiate ivermectin susceptibility in ivermectin sensitive and resistant isolates of Teladorsagia circumcincta and Haemonchus contortus. Parasitology. 2009; v. 136, n. 9, p. 1081.
- 20. Dicker AJ, Nisbet AJ, Skuce PJ. Gene expression changes in a P-glycoprotein (Tci-pgp-9) putatively associated with ivermectin resistance in Teladorsagia circumcincta. International Journal for Parasitology. 2011; 41, 935–942. pmid:21683705
- 21. Sheps JA, Ralph S, Zhao Z, Baillie DL, Ling V. The ABC transporter gene family of Caenorhabditis elegans has implications for the evolutionary dynamics of multidrug resistance in eukaryotes. Genome Biol. 2004;5:R15. pmid:15003118.
- 22. Brenner S. The Genetics of Caenorhabditis elegans. Genetics. 1974; 77(1), 71–94. pmid:4366476
- 23. Simpkin KG, Coles GC. The use of Caenorhabditis elegans for anthelmintic screening. J. Chem. Technol. Biotechnol.1981; 31: 66–69.
- 24. Jorgensen EM, Mango SE. The art and design of genetic screens: Caenorhabditis elegans. Nature Reviews Genetics. 2002; 3, 356–369. pmid:11988761
- 25. James CE, Davey MW. A rapid colorimetric assay for the quantitation of the viability of free-living larvae of nematodes in vitro. Parasitol. Res. 2007; 101, 975–980. pmid:17492469
- 26.
Stiernagle T. Maintenance of C. elegans. In: WormBook: The Online Review of C. elegans Biology [Internet]. Pasadena (CA): WormBook; 2005. 2006 http://www.ncbi.nlm.nih.gov/books/NBK19649/.
- 27. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001; Dec;25(4):402–8. pmid:11846609
- 28. Park IK, Kim J, Lee SG, Shin SC. Nematicidal activity of plant essential oils and components from Ajowan (Trachyspermum ammi), Allspice (Pimenta dioica) and Litsea (Litsea cubeba) essential oils against pine wood nematode (Bursaphelenchus Xylophilus). Journal of Nematology. 2007; vol. 39, no. 3, pp. 275–279.
- 29. Ferreira SR, Mendes TAO, Bueno LL, Araújo JV, Bartholomeu DC, Fujiwara RT. A New Methodology for Evaluation of Nematode Viability. BioMed Research International.2015; pmid:25866820
- 30. Hunt PR, Olejnik N, Sprando RL. Toxicity ranking of heavy metals with screening method using adult Caenorhabditis elegans and propidium iodide replicates toxicity ranking in rat. Food and Chemical Toxicology. 2012; vol. 50, no. 9, pp. 3280–3290.
- 31. Amarante AF, Bricarello PA, Rocha RA, Gennari SM. Resistance of Santa Ines, Suffolk and Ile de France sheep to naturally acquired gastrointestinal nematode infections. Veterinary parasitology. 2004; 120, 91–106. pmid:15019147
- 32. Graef J, Demeler J, Skuce P, Mitreva M, Samson-Himmelstjerna G, Vercruysse J, Claerebout E, Geldhof P. Gene expression analysis of ABC transporters in a resistant Cooperia oncophora isolate following in vivo and in vitro exposure to macrocyclic lactones. Parasitology. 2013; 140, 499–508. pmid:23279803
- 33. Brooker S. Estimating the global distribution and disease burden of intestinal nematode infections: Adding up the numbers—A review. Int J Parasitol. 2010; 40(10): 1137–1144. pmid:20430032
- 34. Murray CJ, Vos T, Lozano R, Naghavi M, Flaxman AD, Michaud C. Disability-adjusted life years (DALYs) for 291 diseases and injuries in 21 regions, 1990–2010: a systematic analysis for the Global Burden of Disease Study 2010. Lancet. 2013; 380: 2197–2223.
- 35. James CE, Davey MW. Increased expression of ABC transport proteins is associated with ivermectina resistance in the model nematode Caenorhabditis elegans. International Journal for Parasitology. 2009; 39, 213–220. pmid:18708066
- 36. Ardelli BF, Prichard RK. Inhibition of P-glycoprotein enhances sensitivity of Caenorhabditis elegans to ivermectin. Veterinary parasitology. 2013; v. 191, n. 3–4, p. 264–75.
- 37. Janssen IJI, Krücken J, Demeler J, von Samson-Himmelstjerna G. Caenorhabditis elegans: modest increase of susceptibility to ivermectin in individual P-glycoprotein loss-of-function strains. Experimental parasitology. 2013; v. 134, n. 2, p. 171–7.
- 38. Yan R, Urdaneta-Marquez L, Keller K, James CE, Davey MW, Prichard RK. The role of several ABC transporter genes in ivermectin resistance in Caenorhabditis elegans. Veterinary Parasitology. 2012; 190, 519–529. pmid:22840641
- 39. Fire A, Xu S, Montgomery MK, Kostas SA, Driver SE, Mello CC. Potent and specific genetic interference by double-stranded RNA in Caenorhabditis elegans. Nature. 1998; vol 391.
- 40. Tamura K, Sudhir Kumar S. Evolutionary Distance Estimation Under Heterogeneous Substitution Pattern Among Lineages. Mol. Biol. Evol. 2002; 19(10):1727–1736. pmid:12270899