PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of ijmbInternational Journal of Microbiology
 
Int J Microbiol. 2010; 2010: 819060.
Published online 2010 July 12. doi:  10.1155/2010/819060
PMCID: PMC2913627

Frequency of Drug Resistance Gene Amplification in Clinical Leishmania Strains

Abstract

Experimental studies about Leishmania resistance to metal and antifolates have pointed out that gene amplification is one of the main mechanisms of drug detoxification. Amplified genes code for adenosine triphosphate-dependent transporters (multidrug resistance and P-glycoproteins P), enzymes involved in trypanothione pathway, particularly gamma glutamyl cysteine synthase, and others involved in folates metabolism, such as dihydrofolate reductase and pterine reductase. The aim of this study was to detect and quantify the amplification of these genes in clinical strains of visceral leishmaniasis agents: Leishmania infantum, L. donovani, and L. archibaldi. Relative quantification experiments by means of real-time polymerase chain reaction showed that multidrug resistance gene amplification is the more frequent event. For P-glycoproteins P and dihydrofolate reductase genes, level of amplification was comparable to the level observed after in vitro selection of resistant clones. Gene amplification is therefore a common phenomenon in wild strains concurring to Leishmania genomic plasticity. This finding, which corroborates results of experimental studies, supports a better understanding of metal resistance selection and spreading in endemic areas.

1. Introduction

Visceral leishmaniasis (VL) constitutes a public health problem in East Africa, Asia, Mediterranean basin, and Central and South America. It is estimated that about 500 000 new cases occur each year. In East Africa and India, VL cause large-scale and tenacious epidemics with high case-fatality rates [1]. In numerous countries that cannot afford the cost of liposomal amphotericin B, first line treatment still relies on sodium stibogluconate or N-methyl glucamine [2]. Cases refractory to antimony treatment have been described for a long time in humans [3] as well as in dogs [4]. During the last decades, incidence of resistance to antimony drugs in Leishmania spp. increased dramatically in some VL foci [5]. The mechanisms contributing to drug resistance in vivo are poorly understood. In vitro studies of strains selected for metal resistance showed that several independent pathways could concur to the resistance phenotype [6]. Among them, overexpression of energy dependent transporters seems to play a major role in resistance to antimonials (reviewed in [7]). Genes coding for “ATP binding cassette” (ABC) transporters [8] and multidrug resistance (MDR) genes [9] have been shown to be amplified as extrachromosomal elements in strains selected in vitro for resistance to heavy-metals. Some works showed that genes that concur to glutathione and trypanothione synthesis can also be amplified in resistant parasites [10]. Gene amplification has been also described in Leishmania strains selected for resistance to methotrexate [11]. Other molecules are involved in trivalent antimony or arsenite compounds transport like aquaglyceroporins whose down regulation provides resistance to trivalent antimony and mutation could affect specifically metal transport [12].

However, all these findings derived from studies of strains whose resistance was induced in vitro and data about natural conditions are lacking. The aim of the present study was therefore to evaluate both the frequency of gene amplification and the amplification level of genes involved in drug resistance among Leishmania clinical strains.

2. Materials and Methods

2.1. Leishmania  Strains

The study included 90 Leishmania strains isolated from French and Sudanese VL foci.

Sixty L. infantum strains were collected in France from 16 HIV positive patients presenting with autochthonous VL. When relapses occurred, we performed successive parasite isolations: for each patient, the number of isolates varied from two to seven and the duration of the survey was 10 to 1800 days (Table 1).

Table 1
Main characteristics of the isolates and results of relative quantification of resistance genes for French VL patients. NA: nonamplified.

Twenty seven strains (fifteen L. donovani, six L. archibaldi and six L. infantum), were isolated from 27 Sudanese patients. Any of these patients were infected by HIV. Three other strains (two L. donovani and one L. infantum) were isolated from Sudanese post-Kala-Azar dermatitidis leishmaniasis (PKDL) patients (Table 2).

Table 2
Main characteristics of the isolates and results of relative quantification of resistance genes for Sudanese patients. All patients were visceral leishmaniasis patients except S28, S29, and S30 who presented with post kala-azar dermal leishmaniasis. NA: ...

All strains were identified, typed, and cryopreserved at the National Reference Centre for Leishmaniasis (Montpellier, France).

A Leishmania infantum strain, isolated from a dog in an area without any therapeutic pressure (MCAN/82/GR/MON497), was used as reference (calibrator) for the relative quantification experiments. In vitro sensitivity assays have confirmed the sensitivity of this strain to antimonials [13]. It was grown using the same conditions as the studied strains.

2.2. DNA Extraction

Parasites were grown in RPMI medium (containing 10% foetal calf serum, gentamycin, all from Sigma); cells were harvested by centrifugation, and washed twice with PBS before DNA extraction. DNA was purified from 5 106 promastigotes by means of a QiaAmp DNA mini kit (Qiagen ref 51034) according to the manufacturer specifications. Elution was performed with 100 μL of DNAse free water and the concentration was measured by spectrophotometry. DNA solutions were kept frozen (−40°C) until use.

2.3. Relative Quantifications of the Genes Under Study

We focused on the following genes: two ATP-dependant drug transporters, MDR 1 and P-glycoproteins P A (PGP A); one gene involved in glutathione and trypanothione synthesis, the G-glutamylcysteine synthetase (GSH1); two genes involved in the folate metabolism pathways, dihydrofolate reductase-thymidlate synthase (DHFR-TS) and pterine reductase (PTR1).

We measured their amplification levels by means of real-time PCR, relatively to invariant genes. We used two single copy genes as reference genes: the DNA polymerase [14] and the ascorbate-dependant peroxidase [15]. Primers and probe targetting the DNA polymerase gene were chosen according to Bretagne et al. [16]. For the ascorbate-dependant peroxidase and the genes under study, specific sequences were extracted from Genbank (http://www.ncbi.nlm.nih.gov/) and primers and TaqMan probes were designed by means of Primer3 software (http://frodo.wi.mit.edu/primer3). Table 3 summarizes the sequences of all primers and probes used in this study.

Table 3
Primers and probes designed from the reference sequences of the studied genes.

All assays were duplex PCR, performed in duplicate experiments, on a Stratagene MX 4000 thermal cycler with a multiplex QPCR ready-to-use mix (Stratagene, ref 600549-51). Each assay was performed by using 0.1 ng of sample DNA as template. Thermal profile was 94°C for 10 minutes for polymerase activation, 94°C 30 s, 54°C 30 s, 45 cycles. All the samples results were compared to those of the reference strain by means of the delta (delta Ct) method; this method calculate, for each sample, the level of the genes of interest relatively to the invariants genes and then standardize the results by means of the calibrator [17].

2.4. Statistical Analysis

Statistical analysis was computed using the SAS Stat 9.1.3 (SAS, Cary, NC, USA) statistical software. Whether the occurrence of gene amplification was associated with Leishmania species was tested for each studied gene using the Pearson Chi-squared test.

3. Results

3.1. Thresholds

For each gene studied, we need to estimate a threshold for amplification. In this aim, we determined the gene amplification levels of the sensitive reference strain from ten biological replicates. We considered this level as a sensitive reference value. The threshold defining significant gene amplification was fixed by calculating the mean value of the 10 replicates augmented of 3 standard deviations to avoid over estimation of gene amplification. Table 4 presents means, standard deviations, and thresholds for the genes of interest; all amplification thresholds were therefore lower than 2, indicating that our assays were able to discriminate at least a two fold gene amplification relatively to the chosen reference.

Table 4
Reference values (means, standard deviations, thresholds) from the sensitive reference strain and amplification coefficients of the clinical primo-isolates displaying amplification ratios above the threshold values calculated from the reference strain. ...

3.2. Evidence of Drug Resistance Genes Amplification in Field Strains

Data concerning the intensity and the frequency of gene amplification in the strains under study are presented in Tables Tables1,1, ,2,2, and and4.4. Statistics took only in account results of the initial strain isolated from the French patients.

MDR 1 gene was at least twice amplified in 30 of the 46 primitive strains, representing 65% of the cases. GSH was amplified in only three Sudanese strains (mean amplification level = 2.41). In the first strain, GSH amplification was associated with PGP A, MDR 1, and PTR1 amplification; in the second, GSH amplification was associated with PGP A and MDR 1 amplification; in the third, GSH amplification was associated with MDR 1 amplification.

Five out of the 16 French patients primitive strains showed DHFR gene amplification; in four of these strains MDR 1 gene was also amplified.

3.3. Stability of Gene Amplification Level among Successive Strains from the Same Patient

We analysed gene copy number variations in the sequential strains from the 16 French patients. Variations of the PGP A, MDR 1, GSH, DHFR and PTR genes amplification with respect to the initial level were less than or equal to 28%, 34%, 25%, 33%, and 10%, respectively, except for the increase of PGP A amplification level observed in one patient (F8). This case corresponded to a patient initially treated with N methyl glucamine and he relapsed 22 months after initial disease. For this patient, the two last strains, collected during relapse, presented with PGP A amplification levels of 7.84 and 7.49 instead of an initial level of 0.45 before the first treatment course (Figure 1).

Figure 1
Variations of PGP A gene amplification among the successive Leishmania isolates from patient F8.

3.4. Gene Amplification Frequency in Leishmania Species

We tested the eventual association between the species L. archibaldi, L. donovani (absent in the Mediterranean focus), L. infantum (present in the two studied foci), and gene amplification. Only MDR 1 amplification, which was observed in 83% of the L. infantum strains, was statistically (P  =  .04) associated with the species (Table 5).

Table 5
Distribution of the isolates displaying amplified genes according to the Leishmania species. *Amplification observed only with L. infantum.

4. Discussion

Several studies aimed to document gene amplification among Leishmania genome of strains selected under in vitro drug pressure and its involvement in resistance mechanisms [18]. The first approach used hybridisation experiments after standardization of the samples by means of a housekeeping gene [8]. More recently, micro array technology was developed; a large number of genes could be analysed into a single experiment, improving the screening capacity. This technology confirmed the involvement of gene amplification in Leishmania resistance to antimony and was complemented by other techniques for finest quantification of the variable genes [10, 19]. A better accuracy is obtained by means of comparative real-time PCR, quantifying the number of copy of a target gene to one or more housekeeping genes.

We analysed amplification of five genes involved in drug resistance on the basis of experimental studies. Our approach was reproducible enough to appreciate a two fold increase in gene amplification but presents some limiting aspects. A clear correlation between clinical response to therapy and gene amplification was impossible to establish because most of the French patients were treated with liposomal amphotericin B and we lack precise information about response to meglumine antimoniate administration in Sudanese cases. Other mechanisms than gene amplification are involved in resistance to treatment such as over expression, the two phenomenons being sometimes linked [20]. In this study, promastigotes were tested after 5 to 8 in vitro passages; we hypothesized that the studied characters were stable, as described by others authors in vivo [21] and in vitro [22]. Nevertheless, we believe our results highlight some significant aspects of the phenomenon of drug resistance of Leishmania that need further discussion.

Our main findings show that gene amplification, initially described in experimental studies about drug resistance of laboratory mutants, can also be encountered in Leishmania field strains.

The most frequently amplified gene was the MDR 1 gene, which was found amplified in 65% of the strains with an average amplification level of 6.98. The MDR 1 gene has been characterized in several Leishmania species. This gene was initially found to be amplified in mutants selected for vinblastine (V circles) [23] resistance and transfection experiments showed that the product of this gene can induce drug resistance by active extrusion of the drug. Experimental studies showed that efflux could be inhibited by verapamil [24]. This gene is highly homologous to the mammalian MDR genes involved in resistance to anticancer drugs. The amplification level of MDR 1 we observed in Leishmania was comparable to those reported for Plasmodium falciparum in experimental studies [25]. Resistance of Plasmodium falciparum to quinine has been associated with PfMDR1 over expression, was directly correlated to gene amplification, and confers cross resistance to other drugs [26]. The selection of MDR amplicons could be explained by the intensive, or perhaps exclusive, use of meglumine antimoniate for VL treatment in Sudan and in veterinary practice in France. So the frequency of MDR amplification may reflects the therapeutic pressure exerted on Leishmania. Katakura et al. [27] obtained stable transfectants for L. amazonensis MDR 1 gene at a copy number over 500 under selective pressure and, for PGP A, Guimond et al. [10] showed that a 2.2-fold increase detected by means of microarray technology was measured as a 12.7-fold increase by Southern blot. Our results showed that high copy number of putatively resistance genes can also be observed in natural conditions; one strain (S30) showed a very high copy number both for MDR 1 and PGP A, respectively, 63- and 22-fold. Performing duplex PCR kept the results clear of a technical error related to differences in sample amounts introduced in the reaction tubes.

The PGP A was found amplified in five strains with an average number of copy of 8.42, similar to the ratios reported by Guimond after selection of metal resistant strains [10]. PGP A belongs to the ATP binding cassette transporters family, a subfamily of multidrug resistance proteins (MRP). PGP A recognizes metal-thiol conjugates and confers resistance to trivalent metals (AS III or Sb III compounds) by accumulating the conjugates in intracellular vesicles close to the flagellar pocket [28]. The level of resistance conferred by transfection of this gene is dependent on the Leishmania species, suggesting that others factors concur to resistance. In another study, Haimeur found that the overproduction of PGP A was associated with GSH 1 gene amplification [8]. Our findings showed that such a concomitant amplification of these genes also occurs in field strains in 3 out of 5 cases of PGP A amplification.

Legare showed that an increase of trypanothione synthesis concur to the detoxification of the reduced antimony compound [29]. GSH1 codes for the heavy subunit of gamma-glutamylcysteine synthase [30], one of the enzymes involved in trypanothione pathway and Grondin et al. showed that GSH1 is overexpressed consecutively to gene amplification in metal resistant strains [31]. In the strains under study, GSH amplification seems to be a minor event and further work has to be made for testing the hypothesis of an overexpression without gene amplification.

Amplification of the DHFR gene was only found in five French L. infantum strains. DHFR-TS is a bifunctional enzyme, involved in the main pathway of folates synthesis; a duplication of its genomic region was first described in methotrexate resistant mutants [32]. One can hypothesise that drug pressure was most important in the Mediterranean area, due to the extended use of trimethoprime in veterinary practice [33]. Guimond et al. found a 7.1 increase for DHFR gene amplification level in strains selected for resistance to methotrexate [10]. We show here that some natural strains presented with the same amplification magnitude (average = 10.6). DHFR amplification and mutations were also found in malaria parasites under natural or experimental drug pressure [34]. These mutations allowed to describe DHFR alleles involved in pyrimethamine resistance and to perform genetic studies about the spreading of resistance from Asia to Africa [35].

Amplification of the PTR1 gene appears uncommon in our collection of natural strains. This finding is consistent with the fact that antifolate chemotherapy is unsuccessful in trypanosomatids and that amplification of this gene would be a physiological consequence of inhibition of DHFR [36].

All these genes are known to be amplified as extrachromosomal elements. The R region for DHFR, the H region for PTR1 and PGP A were described after selection for resistance to methotrexate in L. tarentolae [37] and L. major [38]. In the present study, we do not observe such linkage, even in the F8 patient, showing the selection of amplified PGP A parasites. MDR genes are coded by V circles or D circles, depending on the drug used for their selection. LD1 region was the first extrachromosomal amplicon described in Leishmania and was present in about 15% of the wild-type strains [39]. These elements participate to the plasticity of the Leishmania genome.

Except for MDR 1, statistics showed that gene amplification was not associated with Leishmania species in a significant manner. This might be explained by relative low number of strains displaying genes amplification in the study samples, leading to an insufficient power to detect such association.

Our findings illustrate the relevance of quantitative molecular tools for studying gene amplification potentially involved in Leishmania drug resistance. The high frequency and spreading of amplified genes among field strains is crucial for understanding drug resistance spreading and therapeutic failure in Leishmania species. However, these results, obtained from analysis of the promastigote stage, need to be confirmed on tissue amastigotes and clinical significance of such data warrants to be evaluated in further studies that should compare gene amplification to clinical outcome of antimony therapy. In addition, the fact that other mechanisms which could participate to antimony resistance in natural conditions, such as variable gene expression or altered proteins functionality, remains to be documented.

Acknowledgment

This paper was supported by internal funding from Assistance Publique, Hopitaux de Marseille.

References

1. World Health Organization. Programmes and projects. Leishmaniasis: the disease and its epidemiology. October 2009, http://www.who.int/leishmaniasis/disease_epidemiology/en/index.html.
2. Alvar J, Croft S, Olliaro P. Chemotherapy in the treatment and control of leishmaniasis. Advances in Parasitology. 2006;61:223–274. [PubMed]
3. Bryceson A. A policy for leishmaniasis with respect to the prevention and control of drug resistance. Tropical Medicine and International Health. 2001;6(11):928–934. [PubMed]
4. Gramiccia M, Gradoni L, Orsini S. Decreased sensitivity to meglumine antimoniate (Glucantime) of Leishmania infantum isolated from dogs after several courses of drug treatment. Annals of Tropical Medicine and Parasitology. 1992;86(6):613–620. [PubMed]
5. Sundar S. Drug resistance in Indian visceral leishmaniasis. Tropical Medicine and International Health. 2001;6(11):849–854. [PubMed]
6. Ouellette M, Papadopoulou B. Mechanisms of drug resistance in Leishmania . Parasitology Today. 1993;9(5):150–153. [PubMed]
7. Pérez-Victoria JM, Parodi-Talice A, Torres A, Gamarro F, Castanys S. ABC transporters in the protozoan parasite Leishmania . International Microbiology. 2001;4(3):159–166. [PubMed]
8. Haimeur A, Brochu C, Genest P-A, Papadopoulou B, Ouellette M. Amplification of the ABC transporter gene PGPA and increased trypanothione levels in potassium antimonyl tartrate (SbIII) resistant Leishmania tarentolae . Molecular and Biochemical Parasitology. 2000;108(1):131–135. [PubMed]
9. Henderson DM, Sifri CD, Rodgers M, Wirth DF, Hendrickson N, Ullman B. Multidrug resistance in Leishmania donovani is conferred by amplification of a gene homologous to the mammalian mdr1 gene. Molecular and Cellular Biology. 1992;12(6):2855–2865. [PMC free article] [PubMed]
10. Guimond C, Trudel N, Brochu C, et al. Modulation of gene expression in Leishmania drug resistant mutants as determined by targeted DNA microarrays. Nucleic Acids Research. 2003;31(20):5886–5896. [PMC free article] [PubMed]
11. Grondin K, Kündig C, Roy G, Ouellette M. Linear amplicons as precursors of amplified circles in methotrexate-resistant Leishmania tarentolae . Nucleic Acids Research. 1998;26(14):3372–3378. [PMC free article] [PubMed]
12. Uzcategui NL, Zhou Y, Figarella K, Ye J, Mukhopadhyay R, Bhattacharjee H. Alteration in glycerol and metalloid permeability by a single mutation in the extracellular C-loop of Leishmania major aquaglyceroporin LmAQP1. Molecular Microbiology. 2008;70(6):1477–1486. [PMC free article] [PubMed]
13. Ogunkolade BW, Vouldoukis I, Frommel D, Davoust B, Rhodes-Feuillette A, Monjour L. Immunization of dogs with a Leishmania infantum-derived vaccine. Veterinary Parasitology. 1988;28(1-2):33–41. [PubMed]
14. Croan DG, Morrison DA, Ellis JT. Evolution of the genus Leishmania revealed by comparison of DNA and RNA polymerase gene sequences. Molecular and Biochemical Parasitology. 1997;89(2):149–159. [PubMed]
15. Dolai S, Yadav RK, Pal S, Adak S. Leishmania major ascorbate peroxidase overexpression protects cells against reactive oxygen species-mediated cardiolipin oxidation. Free Radical Biology and Medicine. 2008;45(11):1520–1529. [PubMed]
16. Bretagne S, Durand R, Olivi M, et al. Real-time PCR as a new tool for quantifying Leishmania infantum in liver in infected mice. Clinical and Diagnostic Laboratory Immunology. 2001;8(4):828–831. [PMC free article] [PubMed]
17. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCt method. Methods. 2001;25(4):402–408. [PubMed]
18. Beverley SM. Gene amplification in Leishmania . Annual Review of Microbiology. 1991;45:417–444. [PubMed]
19. Ubeda J-M, Légaré D, Raymond F, et al. Modulation of gene expression in drug resistant Leishmania is associated with gene amplification, gene deletion and chromosome aneuploidy. Genome Biology. 2008;9(7, article R115) [PMC free article] [PubMed]
20. Ashutosh, Sundar S, Goyal N. Molecular mechanisms of antimony resistance in Leishmania . Journal of Medical Microbiology. 2007;56(2):143–153. [PubMed]
21. Dube A, Singh N, Sundar S, Singh N. Refractoriness to the treatment of sodium stibogluconate in Indian kala-azar field isolates persist in in vitro and in vivo experimental models. Parasitology Research. 2005;96(4):216–223. [PubMed]
22. Sereno D, Guilvard E, Maquaire S, et al. Experimental studies on the evolution of antimony-resistant phenotype during the in vitro life cycle of Leishmania infantum: implications for the spread of chemoresistance in endemic areas. Acta Tropica. 2001;80(3):195–205. [PubMed]
23. Gueiros-Filho FJ, Viola JP, Gomes FC, et al. Leishmania amazonensis: multidrug resistance in vinblastine-resistant promastigotes is associated with rhodamine 123 efflux, DNA amplification, and RNA overexpression of a Leishmania mdr1 gene. Experimental Parasitology. 1995;81(4):480–490. [PubMed]
24. Prasad V, Kaur J, Dey CS. Arsenite-resistant Leishmania donovani promastigotes express an enhanced membrane P-type adenosine triphosphatase activity that is sensitive to verapamil treatment. Parasitology Research. 2000;86(8):661–664. [PubMed]
25. Ferreira ID, Do Rosário VE, Cravo PVL. Real-time quantitative PCR with SYBR Green I detection for estimating copy numbers of nine drug resistance candidate genes in Plasmodium falciparum . Malaria Journal. 2006;5, article 1 [PMC free article] [PubMed]
26. Cowman AF, Galatis D, Thompson JK. Selection for mefloquine resistance in Plasmodium falciparum is linked to amplification of the pfmdr1 gene and cross-resistance to halofantrine and quinine. Proceedings of the National Academy of Sciences of the United States of America. 1994;91(3):1143–1147. [PubMed]
27. Katakura K, Iwanami M, Ohtomo H, Fujise H, Hashiguchi Y. Structural and functional analysis of the LaMDR1 multidrug resistance gene in Leishmania amazonensis . Biochemical and Biophysical Research Communications. 1999;255(2):289–294. [PubMed]
28. Légaré D, Richard D, Mukhopadhyay R, et al. The Leishmania ATP-binding cassette protein PGPA is an intracellular metal-thiol transporter ATPase. Journal of Biological Chemistry. 2001;276(28):26301–26307. [PubMed]
29. Légaré D, Papadopoulou B, Roy G, et al. Efflux systems and increased trypanothione levels in arsenite-resistant Leishmania . Experimental Parasitology. 1997;87(3):275–282. [PubMed]
30. Carter KC, Hutchison S, Henriquez FL, et al. Resistance of Leishmania donovani to sodium stibogluconate is related to the expression of host and parasite γ-glutamylcysteine synthetase. Antimicrobial Agents and Chemotherapy. 2006;50(1):88–95. [PMC free article] [PubMed]
31. Grondin K, Haimeur A, Mukhopadhyay R, Rosen BP, Ouellette M. Co-amplification of the γ-glutamylcysteine synthetase gene gsh1 and of the ABC transporter gene pgpA in arsenite-resistant Leishmania tarentolae . EMBO Journal. 1997;16(11):3057–3065. [PubMed]
32. Kündig C, Leblanc E, Papadopoulou B, Ouellette M. Role of the locus and of the resistance gene on gene amplification frequency in methotrexate resistant Leishmania tarentolae . Nucleic Acids Research. 1999;27(18):3653–3659. [PMC free article] [PubMed]
33. Pedersen K, Pedersen K, Jensen H, Finster K, Jensen VF, Heuer OE. Occurrence of antimicrobial resistance in bacteria from diagnostic samples from dogs. Journal of Antimicrobial Chemotherapy. 2007;60(4):775–781. [PubMed]
34. Thaithong S, Ranford-Cartwright LC, Siripoon N, et al. Plasmodium falciparum: gene mutations and amplification of dihydrofolate reductase genes in parasites grown in vitro in presence of pyrimethamine. Experimental Parasitology. 2001;98(2):59–70. [PubMed]
35. Roper C, Pearce R, Nair S, Sharp B, Nosten F, Anderson T. Intercontinental spread of pyrimethamine-resistant malaria. Science. 2004;305(5687):p. 1124. [PubMed]
36. Nare B, Hardy LW, Beverley SM. The roles of pteridine reductase 1 and dihydrofolate reductase-thymidylate synthase in pteridine metabolism in the protozoan parasite Leishmania major . Journal of Biological Chemistry. 1997;272(21):13883–13891. [PubMed]
37. Quellette M, Fase-Fowler F, Borst P. The amplified H circle of methotrexate-resistant Leishmania tarentolae contains a novel P-glycoprotein gene. EMBO Journal. 1990;9(4):1027–1033. [PubMed]
38. Callahan HL, Beverley SM. Heavy metal resistance: a new role for P-glycoproteins in Leishmania . Journal of Biological Chemistry. 1991;266(28):18427–18430. [PubMed]
39. Sunkin SM, McDonagh P, Cunningham ML, Beverley SM, Stuart K, Myler PJ. Conservation of the LD1 region in Leishmania includes DNA implicated in LD1 amplification. Molecular and Biochemical Parasitology. 2001;113(2):315–321. [PubMed]

Articles from International Journal of Microbiology are provided here courtesy of Hindawi