Overexpression of ubiquitin and amino acid permease genes in association with antimony resistance in Leishmania tropica field isolates.

The mainstay therapy against leishmaniasis is still pentavalent antimonial drugs; however, the rate of antimony resistance is increasing in endemic regions such as Iran. Understanding the molecular basis of resistance to antimonials could be helpful to improve treatment strategies. This study aimed to recognize genes involved in antimony resistance of Leishmania tropica field isolates. Sensitive and resistant L. tropica parasites were isolated from anthroponotic cutaneous leishmaniasis patients and drug susceptibility of parasites to meglumine antimoniate (Glucantime®) was confirmed using in vitro assay. Then, complementary DNA-amplified fragment length polymorphism (cDNA-AFLP) and real-time reverse transcriptase-PCR (RT-PCR) approaches were utilized on mRNAs from resistant and sensitive L. tropica isolates. We identified 2 known genes, ubiquitin implicated in protein degradation and amino acid permease (AAP3) involved in arginine uptake. Also, we identified 1 gene encoding hypothetical protein. Real-time RT-PCR revealed a significant upregulation of ubiquitin (2.54-fold), and AAP3 (2.86-fold) (P<0.05) in a resistant isolate compared to a sensitive one. Our results suggest that overexpression of ubiquitin and AAP3 could potentially implicated in natural antimony resistance.

INTRODUCTION

Leishmaniasis is one of the tropical diseases caused by infection with the dimorphic protozoan parasites Leishmania. The clinical manifestations vary from self-healing cutaneous lesion to fatal visceral disease [1]. Cutaneous leishmaniasis (CL) is the most common form of the disease which has been reported from 77 countries with incidence of 1-1.5 million annual new cases [2]. CL is an important public health problem in different parts of Iran and 13-25% of reported CL in some parts of Iran are anthroponotic cutaneous leishmaniasis (ACL) caused by Leishmania tropica [3,4]. Chemotherapy still remains the sole choice for treatment and antimonial compounds such as meglumine antimoniate (Glucantime®) have been in use for more than 70 years [5]. Unfortunately, in recent years antimonial therapy has been imperiled due to an increasing rate of treatment failure attributed to antimony resistant parasites in some endemic foci, especially in India where 37-70% of visceral leishmaniasis patients failed to respond after antimonial therapy [6]. Failure of antimonial treatment was also reported in 10-12% of ACL patients in Iran which was shown to be due to drug resistance [7]. Most of the knowledge regarding antimony resistance process has derived from laboratory generated parasites suggesting presence of different mechanisms such as decreased capacity of converting SbV to SbIII, reduced drug uptake, and elevated drug efflux [8]. In laboratory generated parasites, elevated levels of trypanothione [9] and downregulation of aquaglyceroporin (AQP1) which uptake trivalent metalloid [10] have been demonstrated. Additionally, overexpression of ATP-binding cassette (ABC) protein MRPA which sequestered Sb-thiol complexes has been related to the drug resistance [11].

Recent studies on resistant field isolates have suggested that in addition to the conventional mechanisms other factors such as apoptosis might be implicated in natural antimony resistance [12,13]. Considering that, more insight into the mechanisms responsible for natural antimony resistance could improve current treatment strategies.

Several techniques such as RT-PCR, microarrays, and proteomics are currently used to explore molecular markers of antimony resistance [14,15]. cDNA-amplified fragment length polymorphism (cDNA-AFLP) is one of the most sensitive technologies for detection of differentially expressed gene which does not require any prior information of gene sequences [16, 17]. In this study, cDNA-AFLP was used to identify the potential biomarkers of antimonial drug resistance in L. tropica field isolates.

MATERIALS AND METHODS

Clinical isolates and cell culture

L. tropica field strains used in the present study included meglumine antimoniate sensitive (MHOM/IR/10/Mash-175) and resistant (MHOM/IR/10/Mash-827) isolates which were obtained from patients with ACL in Mashhad, north-east of Iran [7]. Unresponsiveness was defined as patients who regardless of receiving at least 3 full courses of Glucantime® still showed active lesions [7]. The Institutional Ethical Committee of the School of Public Health, Tehran University of Medical Sciences reviewed and approved the study protocol.

The parasites within 6 subcultures after the isolation from the patients were used. Cryopreserved parasites were recovered from liquid nitrogen and cultured in RPMI 1640 medium (Gibco/BRL, Rockville, Maryland, USA) supplemented with 10% fetal bovine serum (FBS; Gibco/BRL), 100 U/ml penicillin, 100 µg/ml streptomycin (Gibco/BRL), in the absence of antimonial drug pressure, and incubated at 26℃.

In vitro drug sensitivity assay

The drug susceptibility of the field isolates in monocyte-macrophage mouse cell-line (J774A.1) was analyzed in triplicate as describe previously with slight modifications [7]. Briefly, J774A.1 (5×104 cells/well) were grown in RPMI with 10% FBS in Labtek 8-well chamber slides (Nunc, New York, USA) and incubated at 37℃ for 24 hr to allow cell attachment, then the cells were infected with stationary phase promastigotes at a parasite to macrophage ratio of 10:1. After 4 hr incubation at 37℃, cells were washed for 3 times and incubated for 72 hr in the presence of serial dilution of Glucantime® (Sanofi-Aventis, Paris, France). Pentavalent antimony concentrations used for the sensitive isolate were 2, 4, 6, 8, 10, and 12 µg/ml, and the for resistant isolate were 35, 40, 45, 50, 55, 60, 65, and 70 µg/ml (the doses were obtained based on previous screening test). The fresh drug was added, and slides were incubated for an additional 72 hr, and then stained using Giemsa. The number of amastigotes in 100 randomly chosen macrophages was counted. The inhibitory concentration 50% (IC50) was calculated using the values for the number of amastigotes/macrophage. IC50 was defined as the effective dose of Glucantime® that reduces the survival of L. tropica by 50%.

RNA extraction

Total RNA extraction was carried out using Tripure isolation reagent (Roch, Mannheim, Germany) according to the manufacture's instructions with minor modifications. Briefly, 1×108 log-phase promastigotes of sensitive and resistant isolates were harvested and lysed using Tripure reagent. The lysate was mixed with chloroform and vortexed vigorously, and then centrifuged at 12,000 g for 15 min at 4℃. Afterwards, the upper colorless phase was precipitated by adding isopropanol and washed with 75% ethanol. Finally, the precipitated RNA was dried at room temperature and dissolved in RNase free water. The quantity and quality of RNA were analyzed using nanodrop (ND-1000, Thermo Scientific Fisher, Waltham, Massachusetts, USA) and gel electrophoresis, respectively. To avoid any genomic contamination, RNA was treated with DNase (Qiagen, Hilden, Germany) according to the manufacture's protocols.

cDNA-AFLP

The cDNA-AFLP assay was conducted as describe previously with some modifications [18]. Single strand cDNA (sscDNA) was synthesized using 10 µg RNA, 20 pmol/µl oligo-dT (Fermentas, Burlington, Canada), 20 pmol/µl random hexamer (Fermentas), 10 mM of dNTP mix (Fermentas) incubated at 65℃ for 10 min, followed by addition of 20 U RNase inhibitor (Fermentas), 4 µL of 5× reverse transcriptase (RT) buffer containing Tris-HCl (pH 8.3) (Fermentas), and 200 U RevertAid premium RT (Fermentas). It was then incubated at 25℃ for 10 min followed by incubation at 50℃ for 60 min. The integrity of cDNA was checked using glyceraldehyde 3-phosphate dehydrogenase (GAPDH) primers as the housekeeping gene (Table 1). The PCR condition was an initial denaturing step of 94℃ for 5 min and 30 repetitions of denaturation at 94℃ for 30 sec, annealing at 60℃ for 30 sec, and extension at 72℃ for 45 sec with a final extension of 72℃ for 7 min.

Table 1

Sequences of adaptors and primer used in cDNA-AFLP and real-time RT-PCR

Double strand cDNA (dscDNA) synthesis was performed using 40 U DNA Polymerase I (Fermentas), 1.5 U Ribonuclease H (Roche, Mannheim, Germany), and 10 mM of dNTP mix at 16℃ for 2 hr, followed by purification with High Pure PCR Clean up Micro Kit (Roche). Double restriction digestion was carried out with 5 U MboI and 5 U EcoRI (Fermentas) on 5 µg purified dscDNA at 37℃ for 2 hr.

Digested dscDNA fragments were ligated to AFLP adaptors (Table 1), 8 µg AD Mbo1, AD EcoR1, 4 µg ad Mbo1, and ad EcoR1 in a final volume of 60 µl. Briefly, it was conducted in the following conditions: 1 min at 55℃, decreasing to 10℃ over 45 min (i.e., 1℃ per 1 min) and then 6 U T4 DNA Ligase (Roche, Mannheim, Germany) was added and incubated at 4℃ overnight. After purification with High Pure PCR Clean up Micro Kit (Roche), pre-amplification with pre-amp primers (Table 1) was performed according to the following conditions: 30 cycles of 94℃, 30 sec, 64℃, 30 sec, and 72℃, 45 sec, and final extension at 72℃ for 7 min.

The sensitive amplification was done using various dilutions of pre-amplified products to get the optimal conditions. PCR amplification was carried out as described above using sensitive primers containing adaptor sequences plus 1 nucleotide at 3' end (Table 1).

The products resulting from sensitive amplifications were resolved on 8% non-denaturating PAGE and stained using silver nitrate. Finally the gels were wrapped and scanned, and then checked for the presence of differentially expressed transcription-derived fragments (TDFs).

Isolation, cloning, and sequencing of TDFs

Selected bands were extracted from the polyacrylamide gels. The eluted DNA was re-amplified with selective primers (Table 1) from resultant profile. After checking PCR products on 1.5% agarose gel, they were cloned into a pGEM-T Vector System I (Promega, Fitchburg, Wisconsin, USA). To identify the isolated TDFs, recombinant plasmid containing an unknown DNA was subjected to sequencing using universal primers (T7 promotor and SP6) (Bioneer, Seoul, South Korea). Homology searches were performed in non-redundant nucleic and protein databases BLAST (http://www.ncbi.nlm.nih.gov/BLAST).

Real-time RT-PCR

In order to confirm the differences in expression patterns of TDF-derived genes between resistant and sensitive L. tropica isolates, real-time RT-PCR was conducted. Target gene primers were designed by Primer 3 software version 0.4.0 (http://frodo.wi.mit.edu) according to the identified genes (Table 1). The cDNA was synthesized from 1 µg of total RNA using QuantiTect® Reverse Transcription Kit (Qiagen, Hilden, Germany) following the manufacturer's instructions. RT-PCR was performed in triplicate with 20 µl volumes using SYBR® Premix Ex Taq™ II (Takara, Tokyo, Japan) in a CFX96™ Real-Time System (C1000™ Thermal Cycler) (BioRad, Hercules, California, USA). The PCR condition was as follows: activation at 95℃ for 30 sec, amplification at 95℃ for 5 sec, 60℃ for 30 sec for 40 cycles. The relative value of the expression level of each gene was calculated by comparing the cycle threshold (CTs) of the target genes with that of the housekeeping gene (GAPDH) using the 2-ΔΔct method [19].

Statistical analysis

In vitro drug susceptibility assay was presented as IC50 which was determined using the linear regression analysis. Each experiment was conducted at least thrice, and the results are expressed as the mean±standard error of the mean (SEM). The significance of differences was determined using the Student's t-test. The level of statistical significance was P<0.05.

RESULTS

Drug susceptibility of Leishmania tropica field isolates

In vitro drug sensitivity of L. tropica to Glucantime® was determined in J774-1 cells. The results showed that IC50 values of Sb (V) for resistant isolate was 50.0±3.5 µg/ml, whereas for sensitive isolate was 6.0±1.5 µg/ml. The IC50 of resistant isolate was about 8-fold higher than the sensitive strain.

The cDNA-AFLP analysis of clinical isolates

To identify differentially-expressed transcripts putatively associated with natural antimony resistance, a cDNA-AFLP approach was utilized on resistant and sensitive L. tropica field isolates. Total RNA was extracted and the sscDNA was synthesized from each sample. The integrity of cDNA was confirmed by specific primers of the housekeeping gene (GAPDH) (Fig. 1).

Fig. 1

Detection of GAPDH as the housekeeping gene (154 bp) by reverse transcriptase-PCR (RT-PCR) on agarose gel (1.5%). M, 50 bp (base pair) molecular weight marker; S, sensitive; R, resistant; NC, negative control.

The cDNA AFLP results revealed differential expressions of TDFs on 8% non-denaturating PAGE by silver staining using different primer combinations (Fig. 2). Differentially-expressed bands were extracted, re-amplified, and cloned into a TA-cloning system. After then, the recombinant plasmids were sequenced. The sequencing results were searched in non-redundant nucleic and protein databases to find homologous sequences and gene identity.

Fig. 2

Expression patterns of TDFs extracted from cDNA-AFLP PAGE. Sensitive amplification of cDNA-AFLP on a PAGE from 3 different primer combinations of S3 MboI/S2 EcoRI, S3 MboI/S1 EcoRI, and S4 MboI/S4 MboI. Arrows point to differentially expressed TDFs which were isolated with code numbers of 1, 2, and 3 (associated genes which derived from these TDFs mentioned in Table 2). M, 50 bp molecular weight marker; S, sensitve; R, resistant.

Blast analysis revealed that 1 clone highly matched with ubiquitine implicated in protein degradation and another clone highly homologous to amino acid permease (AAP3) involved in arginine uptake. Moreover, a clone was found to be similar to a gene which coded a hypothetical protein. The details of sequencing results of identified genes are presented in Table 2.

Table 2

Differentially expressed transcription-derived fragments identified via cDNA-AFLP

Validation of cDNA-AFLP analysis by real-time RT-PCR

In order to validate the cDNA-AFLP results, the mRNA expressions of ubiquitin and AAP3 in resistant and sensitive L. tropica field isolates were analyzed using the real-time RT-PCR. Fig. 3 shows a significant upregulation of ubiquitin (2.54-fold) and AAP3 (2.86-fold) in the resistant isolate compared to the sensitive one (P<0.05).

Fig. 3

The relative expression levels of ubiquitin and amino acid permease (AAP3) genes in sensitive and resistant L. tropica isolates using real-time RT-PCR. The expression of GAPDH was used to normalize the data. Values given are the mean±SD of 3 different experiments (P<0.05). S, sensitive; R, resistant.

DISCUSSION

Emergence of antimony drug resistance has eroded treatment and control of leishmaniasis in endemic regions such as Iran [6,7]. The molecular mechanisms operating in natural antimony resistance are not well known. Different approaches are currently used to elucidate potential factors implicated in this phenomena. In the current study, we utilized cDNA-AFLP to identify genes potentially involved in drug resistance of L. tropica field isolates. We found 1 unknown gene and also 2 known genes including ubiquitin and AAP3 genes. Furthermore, differential expression of known genes was confirmed by real time RT-PCR. Our findings showed a significant upregulation of ubiquitine and AAP3 mRNA expression in the resistant isolate compared to the sensitive one.

Ubiquitin belongs to the class of heat-shock proteins which play important roles in various cellular processes such as degradation of defective proteins, endocytosis, transcription, DNA repair, and apoptosis [20]. One of the well understood function of ubiquitin is protein degradation via ubiquitin-proteasome pathway (UPP) by which it prevents cells from toxic accumulation of abnormal proteins [21]. Through the UPP, ubiquitin attached to lysine residues of the target proteins leading to degradation of the ubiquitin-tagged protein by the 26S proteasome [22]. Several studies have revealed that the wild-type ubiquitin is cytoprotective and vital for cell surveillance [21,22]. Tsirigotis et al. [23] studied the sensitivity of mammalian cells expressing mutant ubiqutin to protein-damaging agents such as cadmium suggested that genetic alteration within ubiqutin coding region inhibits protein degradation and sensitizes cells to toxic effects of aberrant proteins. Another study showed that in mutant ubiquitin-expressing cells the viability was decreased when the cells were exposed to heat shock and hydrogen peroxide [24].

It was demonstrated that metalloid compounds such as arsenic and cadmium induce heat shock and/or oxidative stress responses leading to production of aberrant oxidized proteins [25,26]. Additionally, these compounds inhibit ubiquitination and degradation of aberrant proteins which might lead to oxidative stress-induced apoptosis [27]. It was documented that the elevated expression of ubiqutin in organisms such as Saccharomyces cerevisiae contributes to resistance to oxidative stress through degradation of oxidatively damaged proteins [28]. Moreover, another study demonstrated that overexpression of ubiquitin in yeast increased ubiquitin-dependent degradation of damaged proteins and could prevent heath shock toxicity and cell death [29]. Therefore, our results suggest that the overexpression of ubiquitin in antimony resistant L. tropica isolate could contribute to degradation of oxidized proteins and defend cells from oxidative stress related to antimonial compounds.

Also, our results identified the upregulation of AAP3 in resistant isolate in comparison to the sensitive one. To date several AAPs have been identified in Leishmania [30]. AAP3 is a highly specific arginine transporter that belongs to a large amino acid/auxin permease family (AAAP) which localizes to the surface membrane of Leishmania [31]. Arginine is known as a precursor of polyamine biosynthesis [32] which transferred into Leishmania from external sources through AAP3 [31]. This amino acid is hydrolyzed by arginase producing ornitine which participate in the synthesis of polyamine and trypanothione [32]. The trypanothione is a major reduced thiol of Leishmania parasites and plays a crucial role in detoxification of antimonial drugs [9,33]. Furthermore, the association of increased trypanothione levels with antimony resistance has been observed in Leishmania isolates [9,34]. Haimeur et al. [9] reported that the upregulation of trypanothione in resistant isolates is correlated with overexpression of enzymes involved in the synthesis of polyamines (ornithine decarboxylase; ODC). Based on the above data, it is concluded that an elevated expression of AAP3 in antimony resistant L. tropica isolate might contribute to enhance the level of trypanothione and consequently detoxification of antimonial drugs.

In conclusion, for the first time, we identified ubiquitin and AAP3 genes as potential antimony resistance candidate target genes in Iranian strains of L. tropica through cDNA-AFLP approach. The overexpression of aforementioned genes might be implicated in natural antimony resistance and could be helpful to get more insight into the molecular etiology of antimony resistance in Leishmania spp. However, further studies are required to establish the exact mechanisms by which these genes act in the resistance process.