Translational reprogramming as a driver of antimony-drug resistance in Leishmania

Gutierrez Guarnizo, Sneider Alexander; Tikhonova, Elena B.; Karamyshev, Andrey L.; Muskus, Carlos E.; Karamysheva, Zemfira N.

doi:10.1038/s41467-023-38221-1

Download PDF

Article
Open access
Published: 05 May 2023

Translational reprogramming as a driver of antimony-drug resistance in Leishmania

Nature Communications volume 14, Article number: 2605 (2023) Cite this article

2377 Accesses
1 Citations
77 Altmetric
Metrics details

Subjects

Abstract

Leishmania is a unicellular protozoan that has a limited transcriptional control and mostly uses post-transcriptional regulation of gene expression, although the molecular mechanisms of the process are still poorly understood. Treatments of leishmaniasis, pathologies associated with Leishmania infections, are limited due to drug resistance. Here, we report dramatic differences in mRNA translation in antimony drug-resistant and sensitive strains at the full translatome level. The major differences (2431 differentially translated transcripts) were demonstrated in the absence of the drug pressure supporting that complex preemptive adaptations are needed to efficiently compensate for the loss of biological fitness once they are exposed to the antimony. In contrast, drug-resistant parasites exposed to antimony activated a highly selective translation of only 156 transcripts. This selective mRNA translation is associated with surface protein rearrangement, optimized energy metabolism, amastins upregulation, and improved antioxidant response. We propose a novel model that establishes translational control as a major driver of antimony-resistant phenotypes in Leishmania.

tRNA modification reprogramming contributes to artemisinin resistance in Plasmodium falciparum

Article Open access 17 April 2024

Transcriptional profiling of macrophages reveals distinct parasite stage-driven signatures during early infection by Leishmania donovani

Article Open access 16 April 2022

Analysis by RNA-seq of transcriptomic changes elicited by heat shock in Leishmania major

Article Open access 06 May 2019

Introduction

Leishmaniasis is a neglected tropical disease caused by protozoan parasites of the genus Leishmania. This infection remains a grave public health dilemma, causing ~0.7–1 million new cases per year, with 350 million people at risk of infection¹. Depending on the host immune system status and Leishmania species/strains involved, the disease can take four different clinical forms: cutaneous leishmaniasis (CL), which causes local or diffuse skin ulcers; mucocutaneous leishmaniasis (ML) that affects the mouth and nasal mucosa; visceral leishmaniasis (VL), a systemic infection that affects internal organs containing macrophages; and post-kala-azar dermal leishmaniasis (PKDL), a dermal sequela of VL. VL is lethal without treatment, causing between 20,000 and 30,000 deaths per year².

In the absence of an approved vaccine, the treatment for leishmaniasis relies on chemotherapy. Pentavalent antimonials, commercially administered as sodium stibogluconate (Pentostam) and meglumine antimoniate (Glucantime), have been used as the drug of choice to treat leishmaniasis during the last seven decades³. Nonetheless, the efficacy of antimonials is progressively diminishing, resulting in treatment failure rates of 65% in Bihar, India^4,5,6. The phenomenon of resistance to pentavalent antimonials has been observed in different regions, including Asia, North Africa, and Latin America^4,5,6,7,8. Even though some treatment failure could partially be explained by the host’s immune system and pharmacodynamic or pharmacokinetic factors, the major loss of antimonial efficacy is consistently associated with the parasite’s ability to develop drug resistance⁹. With limited antileishmanial alternatives, there is an urgent need to understand antimonial drug-resistance mechanisms.

There are indications that resistance to antimony is a complex process that involves changes in the genome, metabolome, and lipid remodeling^10,11,12,13. Antimony-resistant parasites have been shown to combat the drug in several ways: reducing drug uptake by downregulation of the aquaglyceroporin 1 transporter; optimizing the antioxidant response by thiol metabolism activation; tuning the energy metabolism to fuel the high energy demanded by the antioxidant response; using the thiols as drug sequesters to form thiol-metal complexes; and increasing the efflux of thiol-metal complexes via exocytosis by overexpression of ABC transporters^14,15,16. Nonetheless, the molecular mechanisms orchestrating the coordinated responses required to combat antimony remain unknown. Understanding these processes is crucial to improve the design of new therapeutic strategies by developing combination therapies to rescue antimonials efficacies, and further exploring the basic biology of Leishmania parasites.

Previous studies have shown that antimony-resistant parasites have differential expression of ribosomal proteins and proteins involved in the regulation of mRNA translation^{17,18,19,20,21}. Nevertheless, there is no clear association between translational control and the parasite’s resistance mechanisms. Interestingly, in the absence of transcriptional control, the parasites mediate changes in gene expression at the post-transcriptional level, whereby the regulation of mRNA translation assumes a central role in influencing the parasite’s phenotypes^22,23.

Recent studies indicate that Leishmania parasites use translational regulation to modulate their phenotype in response to environmental changes²⁴. It was also demonstrated that the protein kinase CDPK1 plays a role in translation and coordinates the Leishmania resistance to paromomycin, an antileishmanial alternative²⁵. However, these studies are mostly rare, and translational regulation of Leishmania’s drug resistance, especially to the main antileishmanial drugs, is still very poorly studied. In this study, we have used a translatomic approach that couples polysome profiling and deep RNA sequencing to assess the changes in the translational efficiency of antimony-resistant parasites growing with or without drug challenge.

In this work, we show that antimony-resistant parasites use translational regulation to coordinate a complex response that mediates the antimony resistance. We have discovered that antimony-resistant parasites exhibit a dramatic remodeling of the translatome, even in the absence of the drug, providing a preemptive adaptation at the level of translation. We have also provided evidence that in retaliation to the drug, a highly selective translation is activated. This response involves alteration of the expression of genes directly required to combat the drug. Translatome analysis allows us to observe changes that would be undetectable by conventional transcriptomic analysis. Our study also reveals that putative translational modulators represent the largest group of detected gene variants, supporting the importance of translational control in the development of drug resistance. We propose a new model of the role of translational control in the coordination of antimony-resistant phenotypes.

Results

Development of antimony-resistant strains

The resistance to antimony is considered a limiting factor for leishmaniasis treatment; however, the parasites’ mechanisms for antimony resistance are poorly understood. Since Leishmania lacks conventional eukaryotic transcriptional control due to its polycistronic nature, gene expression is primarily regulated by post-transcriptional mechanisms. Although the drug resistance in Leishmania has been studied for the past 20 years, the role of translational regulation has remained essentially neglected. In this study, we explored the potential role of translational regulation in the modulation of antimony-resistant phenotypes.

To identify changes at the translational level that are associated with antimony-resistant phenotypes, the drug-resistant strains were derived from the wild-type (WT) Leishmania tropica strain, one of the Leishmania species that can induce both cutaneous and viscerotropic leishmaniasis^26,27,28. Briefly, the antimony-sensitive wild-type strain was treated with an increasing concentration of trivalent antimony (Sb^III), the active form of antimonial drugs, for six months. Moderately resistant (MR) and highly resistant (HR) strains were selected by their half-maximal effective concentration (EC50), estimated by the colorimetric MTT (3-(4,5-dimethylthiazol-2-yl)−2,5-diphenyltetrazolium bromide) cell viability assay (Fig. 1a). WT, MR, and HR strains showed significantly different EC50 values: 10.46 ± 0,65; 260.37 ± 34.93; and 631.73 ± 73.7 µg/mL of Sb^III, respectively (Fig. 1b). The curves of growth for the wild type and the derived resistant strains did not show significant differences (Supplementary Fig. 1). Consistently, the MTT viability assay showed a similar signal intensity of formazan for sensitive and resistant parasites (see left wells in the panels of Fig. 1a). In addition, to test whether the resistance phenotype was stable, we submitted the highly resistant strain to 25 consecutive passages in absence of the drug (Supplementary Fig. 2). Notably, after these passages in culture the HR strain still showed a significantly higher EC50 value when compared to the wild type (WT: 9.34 ± 2.35; HR: 278.6 ± 80 µg/mL of Sb^III). We noticed that even though the EC50 value for the resistant strain was reduced, it still showed about 30-fold higher EC50 than the wild type. Thus, our data reveal that the HR strain has a stable resistant phenotype. An infectivity assay was completed to test whether the selected HR strain shows changes in the infection to a human phagocytic cell line in vitro. The results indicate that the derived resistant strain has a significant reduction in infectivity (Supplementary Fig. 3). As the reduction in the infectivity correlates with the level of resistance (HR versus MR and WT), it is possible that it represents parasite’s trade-off adaptations to acquire drug resistance at the expense of infectivity. Consistently, various studies have previously reported a decreased infectivity in resistant parasites^29,30,31.

Since WT and HR strains showed clear phenotypic differences in terms of antimony susceptibility, they were chosen to conduct polysome profiling experiments followed by deep RNA-seq.

**Fig. 1: Effective stepwise selection of Sb^III-resistant parasites.**

Polysome profiling experimental design and bioinformatic workflow for translatome analysis

Upon completion of drug-resistance validation, Sb^III-sensitive and resistant L. tropica strains were subjected to polysome profiling followed by deep RNA-seq for comparative translatome analyses. mRNA translation is usually conducted by several ribosomes termed polysomes; more ribosomes associated with particular mRNA means mRNA is more involved in translation and likely translated more efficiently (Supplementary Fig. 4A). Polysome profiling is a technique that includes the separation of mRNA transcripts associated with a different number of ribosomes. Polysome profiling followed by deep RNA-seq is a good approach to evaluate mRNA translation status on the full translatome level. It is an adequate technique to study gene expression regulation, especially in Leishmania because this organism has very limited transcriptional control and its gene expression is primarily regulated during translation. The major goal of this work is to delineate the selective translation of mRNAs associated with drug resistance through analysis of mRNA engagement with ribosome and polysomes of the drug-resistant Leishmania strain vs WT and in response to the drug. We wanted to examine two types of changes: changes at the translational level in drug-resistant strain occurring in the absence of the drug (general basal changes), and those only detectable when the resistant parasites are growing under Sb^III exposure (changes to combat the drug). Both of these scenarios were considered in the experimental design (Fig. 2). The highly resistant HR strain was grown with and without the drug, while the WT strain was only grown without Sb^III because of high sensitivity to the drug; the cells were lysed, and lysates were used in polysome profiling experiments. Ribosome subunits (40 S, 60 S), monosomes (80 S), and polysomes were well separated by ultracentrifugation in sucrose gradient as demonstrated by their characteristic A₂₆₀ absorbance peaks (Supplementary Fig. 4B, C). A₂₆₀ spectra and polysome-to-monosome (P/M) ratio for WT and HR strain grown with and without the drug were very similar, indicating that the resistance to Sb^III is not associated with a significant change in overall translational capacity (Supplementary Fig. 5 and Supplementary Data 1). A P/M mean score of 0.54 indicates that during the log phase of growth, L. tropica contains twice as many monosomes as polysomes, independently of the Sb^III-resistance phenotypes or drug challenge.

**Fig. 2: Schematic representation of the experimental design.**

Total mRNAs and fractions corresponding to monosomes (one ribosome or 80 S), light (2–5 ribosomes), and heavy (6 and more ribosomes) polysomes were selected for further analysis (Supplementary Fig. 4C). mRNAs in heavy polysomes were considered to be more engaged in translation and consequently translated more efficiently than in light polysomes or monosomes. Each experimental condition was evaluated in three independent biological replicates, generating a total of 36 samples analyzed by deep RNA-seq (Fig. 2a). The produced reads were mapped against the Leishmania major genome of reference (TriTrypDB-51 Friedlin strain) since it is the most completely annotated Leishmania’s genome, and it is phylogenetically closely related to L. tropica. In consequence, the genes are reported using L. major IDs³². The mapped reads were used for differential translational analysis (DTA) and variant calling analysis (VCA). DTA was used to detect differentially translated transcripts (DTTs) based on the number of reads mapped per transcript and experimental condition. VCA was included to estimate the proportion of DTTs with associated gene variants. Both analyses were compared in terms of the genes and biological processes involved (Fig. 2a). The workflows for DTA and VCA were based on DESeq2 algorithm which uses negative binomial generalized linear models and the FreeBayes algorithm which uses the Bayesian statistic approach, respectively^33,34. Both bioinformatics workflows are described in Supplementary Fig. 6.

To estimate the basal changes in the translatome after the development of drug resistance, the HR strain was compared with the WT strain (control) in the absence of the Sb^III (Fig. 2b). To evaluate the changes in the translatome associated with active drug resistance, the HR growing under Sb^III challenge was compared with the HR without Sb^III treatment (control) (Fig. 2b). More specifically, monosomes were compared against monosomes, and the same type of comparison was done for light and heavy polysomes. As result, a total of six dual comparisons were included to detect DTTs (Fig. 2b). Raw and processed data are available in GEO database repository (GEO ID: GSE173848, BioProject: PRJNA727294).

Next, in order to distinguish changes in translatomes from those derived from differences in mRNA abundance, we compared the total transcriptome (total mRNA used as input for polysome profiling) and heavy polysome mRNA (more efficiently translated transcripts) as shown in Fig. 2c. DESeq2 algorithm was used for both “basal conditions” and “drug challenge”. After DESeq2 analysis and the respective matching, the genes were classified into 4 groups. Group 1: genes that were detected as differentially expressed only in heavy polysomes. Group 2: genes that were detected as differentially expressed only in the total transcriptome. Group 3: genes that were detected as differentially expressed in both total transcriptome and heavy polysomes. Group 4: genes that were not considered as differentially expressed in any of the previously mentioned groups.

Sb^III-resistant strain has a dramatically different translatome in the absence of drug challenge

To detect mRNA-specific translational changes, we first examined changes observed in the absence of drug challenge (basal level, HR compared to WT). The analysis showed that the transcript composition per polysome fraction was drastically different in the Sb^III-resistant strain even in the absence of drug (Fig. 3a and Supplementary Data 2). 49.88% were downregulated and 50.11% upregulated in total polysomes, specifically in light polysomes (906 down; 951 up), and heavy polysomes (1096 down; 1064 up). Monosomes contained only four upregulated DTTs and no downregulated ones (Fig. 3a and Supplementary Data 2). Several transcripts can be detected as a DTT in more than one polysome fraction, therefore, a Venn diagram analysis was performed to identify the distribution of unique transcripts (Fig. 3b and Supplementary Data 3). In total, 2431 differentially translated unique transcripts were detected. 1586 DTTs (65.24%) were identified in both light and heavy polysomes.

**Fig. 3: Sb^III-resistant parasites exhibit dramatic translational reprogramming.**

Next, we performed a comparative analysis of the total mRNA to distinguish changes in translatome from those derived from differences in mRNA abundance. The transcriptomic analysis under basal conditions demonstrates that the vast majority of changes in resistant parasites are observed in translatome (2055 transcripts—translation only, heavy polysomes) rather than transcriptome (62 transcripts—transcription only) and provides strong evidence that translational reprogramming is major driver of antimony-resistant phenotypes in Leishmania (Fig. 3e and Supplementary Data 4). Thus, changes in mRNAs abundance are very limited in the resistant strain, however, engagement of mRNAs in translation is drastically different from the sensitive parasites.

The antimony challenge triggers a highly selective translation in resistant parasites

Next, we sought to determine if resistant parasites modulate the population of transcripts under translation in response to drug challenges. Compared with the dramatic changes detected at the basal level (Fig. 3), a relatively small number of DTTs was detected under the drug challenge, with the majority of DTTs exhibiting an increase (Fig. 3c and Supplementary Data 5). They were mostly associated with polysome fractions: light polysomes (up: 57, down: 9), and heavy polysomes (up: 134, down: 30), while only two DTTs were detected in monosomes (up: 2, down: 0), (Fig. 3c and Supplementary Data 5). Furthermore, the Venn diagram analysis demonstrated that the DTTs correspond to 189 unique transcripts. 122 (64.5%) of these transcripts were exclusively associated with heavy polysomes, where the highest level of translational efficiency is expected (Fig. 3d and Supplementary Data 6). These results suggest that in response to the drug challenge, Sb^III-resistant parasites exhibit a shift towards a highly selective translation to prioritize the efficient translation of a specific population of mRNAs.

Remarkably, even though monosomes are the most abundant type of ribosomes based on Leishmania’s polysome profiles (Supplementary Figs. 1C and 2A), only a small number of DTTs were associated with them (Fig. 3a, c), and the majority of the changes in mRNAs were associated with polysomes, indicating engagement of specific mRNA population in highly efficient translation in response to the drug.

The mRNA abundance analysis identified 104 transcripts differentially expressed in transcriptome only, while analysis of heavy polysome fraction revealed 131 differentially translated transcripts (Fig. 3f and Supplementary Data 7), supporting that translational control plays a predominant role under the drug challenge. However, it appears that changes in mRNA abundance are more pronounced under the drug challenge in comparison with basal conditions (Fig. 3e, f).

Translatome data validation by RT-qPCR and proteomic analyses

Before proceeding with the functional analysis of the translational remodeling, the translatome data were validated by both quantitative reverse transcription-polymerase chain reaction (RT-qPCR) and proteome analyses. The identification of an internal control for RT-qPCR analysis was based on the comparison of normalized read counts. Particularly, a one-way ANOVA analysis showed that mRNA of the protein corresponding to LmjF.35.1945 in L. major had a similar expression level in all analyzed samples, therefore it was selected as an internal control for comparative analyses in RT-qPCR (Fig. 4a and Supplementary Data 8). In addition, OmpA was added as a synthetic mRNA Spike-In during polysome profiling fractionation as a reference for internal control identification (Fig. 4b). To simplify the validation by RT-qPCR, we focused on the changes in heavy polysomes occurring in response to the drug (Fig. 4). A total of 14 genes with different expression profiles were considered for validation and the primer design is summarized in Supplementary Data 9. A Spearman test showed a statistically significant positive correlation (r²: 0.9165) between the fold changes estimated from RNA-Seq and RT-qPCR data, thus validating the translatome analysis (Fig. 4c).

**Fig. 4: Validation of translatome data by RT-qPCR and proteome analyses.**

Next, we asked if the DTTs detected after polysome analysis are also differentially expressed at the total RNA level (input for polysome profiling). For some gene targets, such as TRYR (LmjF.05.0350), RING-P (LmjF.31.1920), and amastin (LmjF.34.1720), the expression level was compared between total RNA input and RNA exclusively isolated from heavy polysome fractions (Fig. 4d). The analysis showed that even though some transcripts do not change at the total RNA level, they were significantly enriched in heavy polysomes in response to drug challenge, indicating a substantially increased translation despite similar total mRNA levels. Widespread uncoupling between transcriptome and translatome is well documented, supporting that many genes are regulated at the translational level³⁵. Translational control is especially important in Leishmania and other trypanosomatids with its polycistronic transcription^36,37, and our data support this.

Since translatome analysis detects mRNA involved in translation, we also explored whether the changes in translational efficiency are in agreement with protein levels using tandem mass spectrometry. While the DTA showed 189 DTTs in resistant parasites growing under drug challenge (Fig. 3d), the proteomic analysis identified only 56 differentially expressed proteins (DEPs). DEPs were identified by statistical filtering using the same parameters as DTTs (absolute fold change ≥1.5 and P value corrected by Benjamini–Hochberg FDR ≤ 0.05). Of these DEPs, 20 were uncharacterized or hypothetical proteins. Therefore, we focused on the remaining 36 DEPs to match proteome and translatome analyses.

The majority of the proteins showed a similar expression pattern across translatome and proteome analyses. The expression pattern of 20 DEPs matched DTAs of Sb^III-resistant parasites grown under drug challenge (Fig. 4e). In addition, 6 DEPs (homocysteine S-methyltransferase, ATP-binding cassette protein subfamily A, putative calpain-like cysteine peptidase family C2, fatty acyl-CoA reductase, Sm protein F, DNA polymerase kappa) showed similar expression in the translatome of resistant parasites growing in the absence of drug challenge. The remaining 10 DEPs did not match the DTA. The comparison between proteome and translatome is summarized in Supplementary Data 10.

The translatome analysis identifies markers associated with antimony resistance

Although the majority of the DTTs identified in response to the drug have still an unknown association with Sb^III-resistance, several of our findings are consistent with previously documented antimony resistance markers. ABC transporters are consistently linked with antimony efflux while the S-adenosylmethionine synthetase is a methyl donor generator that improves the antioxidant response^16,38,39. These genes are frequently upregulated in Sb^III-resistant parasites and consistently were detected as more efficiently translated in response to the drug challenge in our experiments as well (Fig. 5a, c and Supplementary Data 8). The Kinetoplastid membrane protein-11, which has been reported to be downregulated in antimony-resistant parasites⁴⁰, exhibited a significant reduction in the translational efficiency in our study (Fig. 5b and Supplementary Data 8). In addition, several studies have reported the upregulation of several amastins in Sb^III-resistant parasites^41,42,43,44. Interestingly, here we detected a group of 32 amastins displaying a dramatic increase in translation in Sb^III-resistant parasites, both with and without drug challenge (Fig. 5d and Supplementary Data 8). The shift in translatomic profiles was also validated across RT-qPCR analysis for ABC-H transporter and δ-amastin (Fig. 5e, f). Unexpectedly, the aquaglyceroporin 1 (AQP1) that modulates Sb^III uptake and is typically reported as downregulated in antimony-resistant strains¹⁴, showed an unusual tendency of being highly translated and displayed a higher association with light polysomes in Sb^III-resistant parasites challenged with the drug (Supplementary Data 5). The potentially new and previously reported antimony-resistant markers are summarized in Supplementary Data 11. Interestingly, the identified antimony modulators are involved in specific biological processes, such as calcium metabolism, energy metabolism, flagellar component, lipid metabolism, surface protein remodeling, translational control, vesicle transport, and drug efflux.

**Fig. 5: Transcripts associated with the mechanisms of Sb^III drug resistance show a shift in translational efficiency.**

Sb^III-resistant parasites growing without drug challenge display a preemptive adaptation to antimony drug through dramatic changes in translatome

After translatome data validation, we focused on the biological interpretation of the translational remodeling and identification of key regulators potentially contributing to drug resistance. First, we concentrated on the role of the dramatic changes in the translatome observed in Sb^III-resistant parasites growing without drug challenge. Gene ontology (GO) and metabolic pathway enrichment analysis were accompanied by protein–protein interaction networks based on both homology and co-expression to identify key modulators of Sb^III-resistance. While networks based on homology allow for the clustering of structurally and functionally related proteins, the networks based on co-expression analysis can be used to identify genes with a regulatory role in drug resistance with important therapeutic implications⁴⁵.

The enrichment analysis by functional categories, biological processes, metabolic pathways, protein homology, and protein co-expression is shown in Fig. 6. As more than two thousand DTTs associated with polysomes were detected in resistant parasites in comparison with sensitive ones in the absence of Sb^III (basal changes), the significantly enriched GO categories were diverse and complex. However, two types of processes were notably enriched—energy metabolism (green asterisk) and putative translational modulators (red asterisk, Fig. 6).

**Fig. 6: Global translatome analysis of resistant parasites growing without drug challenge reveals complex preemptive adaptations to facilitate drug resistance.**

Several molecular functions and metabolic pathways associated with energy metabolism (Fig. 6a and Supplementary Data 12) were predominantly upregulated and therefore more efficiently translated in resistant parasites (purine metabolism, carbon metabolism, pyrimidine metabolism, biosynthesis of secondary metabolism, and cellular amino acid metabolic process). In addition, various ion-binding proteins and transmembrane transporters were also associated with an increase in drug resistance (Fig. 6b, c). Interestingly, among numerous implicated mechanisms, our analysis showed the enrichment of proteins that can potentially participate in mRNA translation, including “mRNA binding proteins” and “kinase activity” (significantly enriched GO categories associated with gene molecular function); “mRNA processing”, and “aminoacyl t-RNA biosynthesis” (Fig. 6b, c and Supplementary Data 12).

Complementarily, the interactome based on protein homology showed four enriched clusters (Fig. 6d and Supplementary Data 13). Three of these clusters involve potential translational modulators. A first cluster grouped several GTPases. In this cluster, 27 DTTs were identified and most of them were decreased in resistant parasites. Interestingly, besides cytoskeleton rearrangement, GTPases can modulate different stages of protein biosynthesis by guanosine triphosphate (GTP) hydrolysis⁴⁶. A second cluster included a noteworthy group of 62 kinases that were mostly increased in resistant parasites, several of them correspond to serine/threonine-protein kinases that can participate in multiple processes including parasite proliferation, response to stimulus, and translational regulation⁴⁷. A third cluster grouped 15 DEAD/DEAH box helicases that were predominantly increased in resistant parasites. These proteins are typically involved in ribosome biogenesis^48,49,50. Furthermore, other helicase domain-containing proteins, such as eukaryotic translation initiation factor 4E-1 (LmjF.27.1620, downregulated), 4E-2 (LmjF.19.1500, downregulated), and 4E-3 (LmjF.28.2500, upregulated) were also differentially translated. Finally, a fourth cluster contained mostly downregulated DTTs in resistant parasites growing without drug challenge. It included 15 oxidoreductases, mainly constituted by quinonoid dihydropteridine reductases (QDPR), enzymes that participate in the pathway that recycles the tetrahydrobiopterin (BH4), an antioxidant compound that mediates the H₂O₂ susceptibility⁵¹. Interestingly, the interactome based on protein co-expression showed that a group of at least 13 proteins involved in translation was significantly enriched representing multiple interactions in the network (Fig. 6e and Supplementary Data 14). Particularly, the changes in ribosomal proteins and translation initiation factors could have a key role in the observed translational remodeling, possibly as a preemptive adaptation to drug challenges.

Thus, our analysis revealed that even in the absence of the drug the translatome of resistant parasites is dramatically different from the sensitive strain. As a result, many important processes in the cell, such as energy and amino acid metabolism, transmembrane transport, mRNA binding, mRNA processing, and putative translational modulators are substantially modified.

Drug-resistant parasites activate a highly selective translation in response to antimony drug

As described above the resistant parasites displayed a dramatic remodeling of translatome in comparison with sensitive parasites even in the absence of drug thus supporting that preemptive adaptation is necessary to combat the drug efficiently. Interestingly, their response to the drug challenge was very selective and included only 189 unique DTTs (Fig. 3d). To understand the biological significance of translational changes detected in Sb^III-resistant parasites in response to the drug further analysis was performed. As described before, GO and protein–protein network interaction analyses were carried out to identify key players in Sb^III-resistant parasites in response to drug challenges.

The significantly enriched GO categories associated with cellular components showed notable changes in the translational efficiency of mRNAs encoding proteins associated with microtubule organization (cytoplasmic vesicle, flagellum, and vacuoles) (Fig. 7a and Supplementary Data 12). The significantly enriched GO categories associated with molecular functions identified a highly efficient synthesis of proteins involved in antioxidant response (methionine adenosyltransferase activity, glutathione peroxidase activity, peroxidase activity) and ion-binding proteins (metal ion binding, and cation binding) (Fig. 7b and Supplementary Data 12). In addition, components of the antioxidant response were also enriched in biological process GO categories (S-adenosylmethionine metabolism, response to reactive oxygen species, response to oxidative stress, response to stress, etc.) (Fig. 7c and Supplementary Data 15). Finally, the STRING interaction network showed significant enrichment of proteins associated with glutathione metabolism, equivalent to trypanothione metabolism of Leishmania parasites, the main antioxidant response in trypanosomatids. These 11 nodes are all highly interconnected (Fig. 7d and Supplementary Data 15).

**Fig. 7: Resistant parasites modulate the synthesis of proteins involved in interconnected processes to combat Sb^III challenge.**

Since GO analysis is limited by gene annotation, we manually grouped the DTTs detected under antimony challenge in general functional categories based on gene annotation. The DTTs with known or putative functions were grouped into 9 general functional categories which were clearly interconnected. Balanced “calcium metabolism” (9 genes) and “antioxidant response” (11 genes) can participate in mediating the oxidative stress typically induced by the antimony and promote the drug inactivation via thiol metabolism⁵². Optimized “energy metabolism” (9 genes) could supply the high-energy cost required by the antioxidant response activation and other drug depuration mechanisms¹⁰. Since the flagellar movement is associated with a high-energy cost, the decrease in “flagellar components” (8 genes) can contribute to optimizing energy metabolism by conserving ATP²⁴. Complementary experiments showed that decreased flagellar components correlate with a smaller flagellum in resistant parasites growing under drug challenge (Supplementary Fig. 7). “Vesicle transport and drug efflux” (5 genes) can mediate drug removal¹⁶. The largest group included “Surface protein” (37 genes) that can modulate the parasite’s membrane composition and drug transport⁵³. “Lipid metabolism” can contribute to membrane remodeling and energy metabolism optimization¹⁰. “Translational control” (12 genes), can potentially participate in the activation of highly selective translation in response to the drug (Fig. 7e). These results support the idea that the mechanisms to combat the Sb^III require complex, coordinated processes. The DTTs manually grouped are summarized in Supplementary Data 16.

Differential translational and variant calling analyses detected different genes involved in common biological processes

Finally, we asked if the observed translatome remodeling correlated with changes at the genomic level by the detection of gene variants exclusively observed in resistant parasites. After filtering, a total of 1626 variants were detected. Of them, 623 (38.31% of the total) were homozygous. In terms of the type of variant, most of them correspond to single nucleotide polymorphisms or SNPs (n = 1254 or 77.12% of the total). Interestingly, concerning the genomic region, the variants were mostly distributed in regulatory regions, upstream of the annotated genes (1195 or 73.49% of the total). This observation is consistent with recent studies demonstrating that mutations occur less often in coding sequences of eukaryotic organisms⁵⁴. Moreover, only 186 (11.44% of the total) missense variants were detected. Particularly, 57 missense variants were homozygous, suggesting a higher probability to impact the parasite’s phenotype and have a protagonist role in Sb^III resistance (Supplementary Fig. 8 and Supplementary Data 17).

Then, we estimated the percentage of DTTs with associated gene variants to analyze the correspondence between DTA and VCA. The results showed that most of the DTTs do not contain associated variants. At the basal level, 309 DTTs (12.71% of the total of 2431) showed associated variants (Fig. 8a and Supplementary Data 18). Under the drug challenge, 19 DTTs (10.05% of the total of 189) showed associated variants, and most of them correspond to amastins (Fig. 8b and Supplementary Data 18). These results indicate that genomic changes can only partially explain the dramatic remodeling detected at the translational level. Interestingly, the manual grouping of variants by general functional categories showed that the drug-resistant strain exhibits variants in genes that are associated with functional categories that were also enriched in the translatome analyses (Figs. 7 and 8). These categories include putative translational control, surface proteins, energy metabolism, and kinase/phosphatase proteins (Fig. 8c and Supplementary Data 18). This evidence indicates that during the development of resistance, continuous drug exposure induces selective pressures on some parasite’s biological processes that gradually leads to a coordinated rearrangement at both genomic and translatomic levels.

**Fig. 8: Only a minor group of DTTs showed associated variants but the functional categories involved are consistent between DTA and VCA.**

Discussion

Leishmaniasis is a complex disease, with cutaneous, mucocutaneous, and visceral forms, that is caused by Leishmania parasites and is difficult to treat due to the widespread emergence of drug-resistant strains. Whereas cutaneous leishmaniasis is the most common form of the disease affecting humans, visceral leishmaniasis is the most life-threatening, ranking second to only malaria as the most fatal of parasitic infections. Drug resistance in Leishmania parasites remains a restrictive factor in leishmaniasis treatment; however, the basic molecular mechanisms orchestrating the development of resistance remain poorly understood. In the absence of a vaccine, leishmaniasis disease treatment relies heavily on chemotherapy. Current chemotherapy drugs used to treat leishmaniasis are limited. These include pentavalent antimonials (Sb^V), as the primary treatment, and amphotericin B, miltefosine, pentamidine, and paromomycin as alternative treatments⁵⁵. However, treatment of leishmaniasis is obstructed by drug toxicity, high cost, and treatment failures caused by drug resistance. Furthermore, cross-resistance phenomena have been documented, suggesting a high capacity of Leishmania to adapt to different stressors and drugs^56,57. Resistance to antimony, the primary antileishmanial chemotherapy drug, is related to genomic rearrangements that can modulate the dosage and function of key proteins^11,43,58. However, genomic rearrangements and mutations do not explain all possible drivers of resistance mechanisms. Leishmania, and other trypanosomatids, have polycistronic transcription and the use of transcriptional control is limited⁵⁹. As a result, Leishmania, and other trypanosomatids, use primarily translational control to regulate gene expression^60,61,62. Therefore, besides the changes at the genomic level, translational control could represent the major route to combat the drugs. Nevertheless, the role of translational regulation of drug resistance in Leishmania has not been investigated yet. The major advantage of translational regulation is that it can be easily reversed and readjusted under different environments⁶³. This study describes a new perspective on drug resistance in Leishmania parasites, highlighting translational reprogramming as a crucial driver of antimony-resistant phenotypes.

In this study, a L. tropica strain that is highly resistant to antimony drug was generated on the base of a sensitive strain upon a gradual increase in drug concentration (Fig. 1). The changes in the translatome were studied by coupling polysome profiling and deep RNA-Seq (Supplementary Fig. 4). Two types of comparison were used to identify changes in the translatome exhibited by the derived resistant strain. First, the resistant strain was compared with the sensitive strain to identify basal translatome changes in the absence of drug challenge. Then, the translatome of resistant parasites treated with the antimony drug was compared with the translatome of untreated resistant parasites to identify potential drug-resistance modulators (Fig. 2).

Although the overall translational capacity remains unchanged (Supplementary Fig. 5) and the population of transcripts associated with monosomes remained mostly unaffected (Fig. 3), the population of efficiently translated transcripts (associated with light and heavy polysomes) was significantly different in antimony-resistant parasites under both, basal changes and in response to the drug. Notably, the basal changes of the translatome in the absence of the drug included 2431 differentially translated transcripts in comparison with sensitive parasites (Fig. 3a, b). This suggests that the development of antimony resistance involves dramatic reprogramming of mRNA translation. Our comparative analysis of transcriptome versus translatome reveals very limited changes in mRNA abundance and is in a strong favor of our model supporting that translational reprogramming plays a major role in drug resistance in Leishmania parasites (Figs. 3e, f and 9). The basal changes of the translatome observed in resistant parasites could represent the parasite’s strategies to efficiently compensate for the biological fitness once they are exposed to the antimony (Fig. 9a). In our previous work using the same resistant and sensitive strains, we observe that lipid and metabolite composition is different in the resistant strain even in the absence of the drug correlating well with dramatic changes in translatome observed under basal conditions^10,13. In future studies, it would be interesting to examine several resistant cell lines derived from the same parental sensitive strain and determine if parasites always use the same route of translational reprogramming to develop drug resistance or different solutions are possible to achieve the same outcome by means of translational regulation. The association between drug-resistant phenotypes and reprogrammed mRNA translation has been previously observed in refractory cancer cells. Translational reprogramming is commonly recognized as a source of adaptive plasticity that allows cancer cells to become resistant to the new therapies^45,64,65,66. The identification of key translational modulators associated with drug resistance in Leishmania parasites is essential for the development of new therapeutic strategies.

**Fig. 9: Model of the translational regulation as a central driver of *Leishmania*’s antimony-resistant phenotypes.**

Resistant parasites that are grown while exposed to the antimony induced changes in the translation of a specific population of only 189 transcripts. These drug-resistance modulators were grouped by interconnected biological processes, such as surface protein remodeling, improved energy metabolism, oxidant response, drug inactivation, and drug efflux (Fig. 9b). Together, our results suggest that Leishmania parasites use multifaceted mechanisms to combat the drug’s effect in a coordinated effort. In this study, we have shown that antimony-resistant parasites optimize energy metabolism through changes in translational efficiency (Figs. 6c, 7e, and 8c). In agreement, our previous findings indicated that an optimized carbon-energy metabolism achieved through lipidome and metabolome changes is essential to fuel the high-energy demands required to cope with antioxidant-induced stress and drug efflux via ABC transporters^10,13 (Fig. 8). In addition, resistant parasites downregulate flagellar components, a strategy that could promote an optimized energy expenditure (Fig. 7a). We observed that this event correlates with the resistant parasite’s phenotype of possessing a shorter flagellum, which is associated with lower flagellum activity and lower energy expenditure (Supplementary Fig. 8). Based on the increase of proteins argininosuccinate synthase (LmjF.23.0260), ascorbate peroxidase (LmjF.34.0070), gamma-glutamylcysteine synthetase (LmjF.18.1660), methylenetetrahydrofolate reductase, putative, peptide methionine sulfoxide reductase-like, S-adenosylmethionine synthetase (LmjF.30.3520, LmjF.30.3500), SelR domain-containing protein (LmjF.28.2660), tryparedoxin peroxidase (LmjF.26.0800, LmjF.26.0820), it is evident that the resistant parasites preferentially utilize trypanothione metabolism as an antioxidant response (Fig. 8d). These results have been consistently verified by several studies claiming that trypanothione metabolism is involved in drug inactivation through the formation of thiol-metal complexes, which are necessary for drug efflux^15,67,68,69. We also identified proteins involved in drug efflux via exocytosis, such as ATP-binding cassette protein subfamily H (LmjF.11.0040), QA-SNARE (LmjF.29.0070, LmjF.19.0120), Qb-SNARE (LmjF.19.0010), and vesicle coat proteins (LmjF.29.1810)⁷⁰. Downstream, drug efflux requires remodeling of the parasite’s membrane to allow for the fusion of vesicles containing the thiol-metal complexes with the plasma membrane¹⁰. In this study, we identified multiple proteins that were highly translated in response to the drug are involved in the alteration of the membrane composition, such as aminophospholipid translocase (LmjF.34.3220), lipase (LmjF.31.2460), lipase (LmjF.13.0200), several amastins located in the chromosome 8 and 34, surface antigen-like protein (LmjF.09.0580), and the GRAM domain-containing protein (LmjF.04.0240). Particularly, the translatome analysis showed that the amastin surface proteins are preferentially translated in response to drug challenge. This observation indicates that beyond its role in vertebrate host infection, amastins are important drug-resistance modulators. Though the literature supports the overexpression of amastins in antimony-resistant parasites, their role during drug modulation remains poorly investigated^41,42,43,44.

The gene variant analysis showed that most of the differentially translated transcripts do not have changes at the genomic level. However, genes with detected variants are associated with similar functional categories as those enriched by translatome analysis. These results suggest that the antimony-resistant phenotype is predominantly regulated by translational control.

Throughout this study, we utilized different strategies of analysis, such as gene ontology, protein clustering by homology, and co-expression network to identify putative translational modulators. Overall, we have shown that among several differences, resistant parasites exhibit consistent changes in the expression of putative translational regulators, including the enrichment of protein kinases, helicases, mRNA binding proteins, translation initiation factors, and ribosomal proteins (Figs. 6, 7 and Supplementary Data 15). Interestingly, recent studies have shown that Leishmania parasites can use translational master regulators, such as the calcium-dependent kinase CDKP1 to modify the landscape of translating mRNAs modulating the resistance to paromomycin²⁵. In addition, ribosomal proteins, such as eIF2alpha, are crucial for the translational reprogramming that occurs during parasite stage differentiation⁷¹. New studies in human cancer cells highlight the role of altering translational factors, such as eIF4E or eIF3, to induce dramatic translatomic changes that modulate drug-resistant phenotypes^72,73. The activation of transporters and the antioxidant response demands a lot of energy in Leishmania parasites. The optimization of energy metabolism of resistant parasites observed during basal changes supports the idea of a preemptive adaptation.

Overall, our study demonstrates that during the development of drug resistance, resistant parasites undergo a dramatic translatome remodeling. After stepwise selection, resistant parasites exhibit a plethora of changes directed by translational regulators. Translational regulators themselves undergo changes occurring at both genomic and translational levels and influencing a complex remodeling. The fact that resistant parasites show dramatic changes in the translatome even in the absence of the drug challenge suggests that this remodeling could work as a parasite’s strategy to be prepared for drug exposure. Therefore, we hypothesize that these changes are necessary to compensate for the loss of biological fitness of the parasites exposed to the drug for both immediate survival and as a preemptive adaptation to the drug challenge (Fig. 9a). Our recent studies uncovered dramatic lipidome and metabolome remodeling even in the absence of antimony drug in the same drug-resistant Leishmania tropica parasites supporting the idea that those changes could be essential preemptive adaptations^10,13. Our current findings suggest that translatome reprogramming could direct changes in lipidome and metabolome in resistant parasites. These coordinated changes preemptively prepare parasites to counteract the antimony drug very efficiently. Resistant parasites growing under drug challenge show targeted changes in the translational efficiency of only a few specific modulators of drug resistance (Fig. 9b). As a result, 33 transcripts downregulate their translation, and a specific population of 156 transcripts shifts towards a highly efficient translation by forming heavy polysomes. The produced proteins act in interconnected biological processes that modulate the effective response to the drug. An optimized energy metabolism supplies the energy required to fuel the antioxidant response via thiol metabolism; in turn, thiols further inactivate the drug by the thiol-metal complex conformation; ABC transporters and the traffic by vesicles induce the drug efflux. Drug efflux via exocytosis can be favored by the remodeling of membrane and surface protein. When transcriptional control is limited, translational remodeling offers an energy-saving alternative to efficiently respond to the drugs. Thus, highly selective translation activated by the drug leads to a coordinated response to combat the drug that involves interconnected biological processes.

Deciphering molecular mechanisms of selective mRNA translation to combat stress induced by antimony holds a great promise in the development of innovative strategies for the treatment of leishmaniasis based on targeting master regulators of drug-induced selective translation and is a matter of future investigation.

Methods

Reagents

Schneider insect medium (Sigma-Aldrich, #S0146), fetal bovine serum or FBS (Sigma-Aldrich, #F0926), Dulbecco′s Phosphate Buffered Saline (Sigma-Aldrich, #D8537), potassium antimony (III) tartrate trihydrate (Sigma-Aldrich, #230057), penicillin/streptomycin (Sigma-Aldrich, #P0781), MTT or 3-[4,5-dimethylthiazol-2-yl]-dephenyltetrazolium bromide (Sigma-Aldrich, #M2128), cycloheximide (Sigma-Aldrich, C7698), TRIzol™ LS Reagent (Invitrogen, #10296028), Power SYBR Green PCR Master Mix (Applied Biosystems, #4367659), Direct-zolTM RNA MiniPrep (Zymo Research, #R2052); poly(A) tailing Kit (Thermo Fisher Scientific, AM1350), HEPES (Fisher Scientific, BP310), Dithiotreitol or DTT (Fisher Scientific, BP172), Nonidet P 40 or NP-40 (Thermo Fisher Scientific, #J610055), potassium Chloride (KCl) (Sigma-Aldrich, P9541), magnesium chloride hexahydrate (Acros Organics, #413415000), recombinant Rnasin Ribonuclease Inhibitor (Promega, #N251B), cOmplete Mini EDTA-free protease inhibitors (Roche Diagnostics, #11873580001).

Leishmania parasites culture

Leishmania tropica promastigotes were grown at a density of 5 × 10⁵ cells/mL, in Schneider Insect medium, supplemented with 10% Fetal Bovine Serum (FBS), 10 units/mL Penicillin and 0.1 mg/mL streptomycin.

Antimony resistance stepwise selection

The trivalent antimonial (Sb^III), considered to be the active form of the drug, was administered as potassium antimony (III) tartrate trihydrate⁷⁴. The stepwise Sb^III resistance selection was started at 10 µg/mL of Sb^III, equivalent to the EC50 estimated for the wild-type strain. The derived parasites were sequentially treated with 10, 20, 30, 40, 50, and 100 µg/mL of Sb^III. The MTT cell viability assay was used to estimate the response to the drug at three different points in the kinetics of resistance development: initial point (Sb^III-sensitive strain or WT), middle point i.e., parasites derived from treatment with 50 µg/mL of Sb^III (Sb^III-moderately resistant strain or MR), and final point i.e., parasites derived from treatment with 100 µg/mL of Sb^III (Sb^III-highly resistant strain or HR)⁷⁵. In addition, response to the drug was measured in a similar manner in a highly resistant strain after culturing HR parasites for 25 passages in the absence of the drug. EC50 values were estimated using Probit regression or three-parameter dose-response model. CellTrace^TM CFSE (Thermo Fisher Scientific) was used in cell proliferation assay following manufacture instructions. To estimate differences in parasite growth curves after drug-resistance selection, we used a Neubauer chamber to count the parasite density for 8 days post-seeded in 5 ml of Schneider Insect medium using T25 flasks.

Infectivity assay

The infectivity assay with Leishmania parasites was based on Palacios et al., 2017 protocol⁷⁶. The U937 cell line exhibiting monocyte morphology (ATCC: CRL-1593.2™) was stimulated with phorbol 12-myristate 13-acetate (Sigma-Aldrich®, # P8139) and seeded at 3 × 10⁵ cells per well in 24 well plates (Corning®, # 3514). The cells were incubated at 37 °C and 5% CO2 for 72 h to promote differentiation to phagocytic cells. For in vitro infection, Leishmania promastigotes were cultured in Schneider’s Insect Medium and harvested after 5 days of culture (stationary phase). The infection was initiated by mixing parasites and phagocytic cells in a ratio 1:15 and resuspended in RPMI 1640 medium (Gibco, # 11875093) without FBS. The infection was completed at 34 °C and 5% CO2 for 3 h. Next, the cells were washed with PBS to remove free-moving parasites. Infected cells were resuspended in RPMI supplemented with 10% FBS and incubated at 37 °C for 24 h. Finally, the cells were washed with PBS, fixed with 100% methanol, and stained with Giemsa following the commercial instructions (Sigma-Aldrich, # 64571-95). The number of infected cells was determined by light microscopy using ×1000 total magnification. A total of 100 cells were counted per condition. The experiments were completed in four biological replicates with two technical replicates.

Polysome profiling

Three experimental conditions were analyzed, (1) the WT strain growing without drug treatment, (2) the HR strain growing without drug treatment, and (3) the HR strain treated with 100 µg/mL of Sb^III for 12 h after 36 h of culture to induce the mechanisms necessary for drug resistance. In each case, the promastigotes were seeded at a density of 5 × 10⁵ promastigotes/mL in a final volume of 60 mL of Schneider´s Insect medium supplemented with FBS.

Parasite’s lysate preparation was based on our previous studies⁷⁷. After 48 h of growth at 26 °C (mid-log phase), the cultures were treated with cycloheximide at a final concentration of 100 µg/mL for 10 min at 26 °C to arrest the ribosomes on translated mRNAs. The number of parasites was estimated by hemocytometer using 3.5% formaldehyde solution (Electron Microscopy Sciences, #15710) for cell fixing. Based on preliminary experiments, an optimal parasite concentration in terms of the polysome profiling quality ranged from 1.5 to 2 × 10⁷ promastigotes/mL in 60 mL of medium.

After cell counting, a total of 7.5 × 10⁸ promastigotes were harvested per experimental condition by centrifugation at 1.800×g and 4 °C for 8 min and washed in 30 mL of DBPS containing 100 µg/mL cycloheximide. The cells were centrifuged at 1.800×g and 4 °C for 5 min and resuspended in 1 mL of lysis buffer (20 mM HEPES-KOH (pH 7.5), 10 mM MgCl₂, 100 mM KCl, 2 mM DTT, 1% NP-40, 1× protease inhibitor cocktail (EDTA-free), 80 units/ml RNasin, and 100 µg/ml of cycloheximide). Details of cell lysate preparation are described previously⁷⁷. The lysate was centrifuged at 11,200×g and 4 °C for 10 min and the supernatant was kept on ice. The supernatant was normalized based on the absorbance at 260 nm prior to sucrose-gradient ultracentrifugation.

To separate polysomes, ~500 µl of the clarified lysates was loaded on the top of 10–50% sucrose gradient containing 20 mM HEPES-KOH (pH 7.5), 100 mM KCl, 10 mM MgCl₂, 1 mM DTT, 80 units/ml RNasin, and subjected to ultracentrifugation using SW41 rotor for 2 h at 260,000×g and 4 °C. After centrifugation, ~500 µl of fractions were collected using a Piston Gradient Fractionator (BioComp Instruments). The polysome profiles were monitored by the absorbance at 260 and 280 nm using the Triax™ Flow Cell 1.56 A software. The absorbance at 260 nm Vs. fractions number plot was generated to identify the ribosome distribution. Trizol LS was added to fractions immediately and samples were kept at −80 °C until RNA extraction⁷⁷. The fractions covering the peaks associated with monosomes (one ribosome), light polysomes (2–4 ribosomes), and heavy polysomes (5–9 ribosomes) were analyzed and pooled per each independent experiment.

RNA isolation

Monosome, light, and heavy polysome pools were used for RNA isolation by Direct-zolTM RNA MiniPrep following the manufacturer’s instructions. For RT-qPCR analysis, synthetic and polyadenylated outer membrane protein A (OmpA) mRNA was spiked into all samples before RNA purification and used for normalization as previously described⁷⁷. RNA was recovered in 50 µL nuclease-free water, quantified by Qubit 3.0 (Life Technologies), and compared to reference samples. A total of 2.5 ng of polyadenylated OmpA were added per 1 µg of total RNA before sequencing.

High-throughput sequencing

A total of 36 samples were prepared for sequencing. Sample quality was assessed by Agilent RNA 6000 Nano Reagent on Agilent 2100 Bioanalyzer (Agilent Technologies) and quantified by Qubit RNA HS assay (Thermo Fisher Scientific, #Q32857). Enrichment of poly-A RNA was performed by NEB Next® Poly(A) mRNA Magnetic Isolation Module (New England BioLabs Inc., #7490). Subsequently, library preparation was performed with 2 μg poly-A RNA by NEB Next Ultra II RNA library prep non-directional (New England BioLabs Inc.) following manufacturer instructions. Paired RNA-seq was carried out in an Illumina Hiseq 2000 platform using HCS2.2.68 software, with read length of 150 nucleotides, and an estimated of 40 M PE reads per sample (20 M in each direction).

Bioinformatics analysis of translatome

All bioinformatic analyses were completed using the Galaxy platform (https://usegalaxy.org/)⁷⁸. The paired reads were initially evaluated by the FastQC v0.73 (Babraham Bioinformatics) to verify the quality and the presence of adapters⁷⁹. Then, the adapters were removed and the reads were filtered (Phred quality score >30, read length ≥50 nucleotides) using Trim Galore wrapper script (Version 0.6.3)⁸⁰. After the processing of raw reads, the read quality was confirmed and summarized using the MultiQC algorithm v1.11 (Babraham Bioinformatics)⁸¹. Then, the paired reads were mapped to the L. major Friedlin strain genome (available in tritrypdb.org/ and released on 2019-11-04) using BWA-MEM (Version 0.7.17.1) mapper⁸², the derived.bam files were submitted to coordinated sort using the function SortSAM of Picard tools kit (Version 2.18.2.1)⁸³. MarkDuplicates of the Picard tool kit (2.18.2.2) were used for PCR duplicate removal. The counting matrices (raw number of reads mapped per transcript) were generated using the feature Counts algorithm⁸⁴. The counting matrices were normalized using Deseq2 v1.34.0 normalization for differential expression analysis³³. To detect statistically differentially expressed genes per comparison, a fold change cutoff of ≥1.5 (change in gene expression greater than 50%) and a Benjamini–Hochberg adjusted P value <0.05 were fixed as threshold. DTTs with statistically significant changes were filtered by two parameters: the significance or P value score estimated by a Wald test and corrected by the Benjamini–Hochberg method, and the fold change score estimated by the normalized read counts ratio.

Validation of DTTs by RT-qPCR

cDNA samples were prepared using High-Capacity cDNA Reverse Transcription Kit (Applied Biosystems, #4368813). Reverse transcription-quantitative polymerase chain reactions (RT-qPCR) were done on Quant Studio 12 K Flex Real-Time PCR System using Power SYBR Green PCR Master Mix (Applied Biosystems, #4367659) according to the manufacturer’s protocol. The differential gene expression was estimated using the ΔΔCT method⁸⁵. For this purpose, mRNA for the hypothetical protein LmjF.35.1945/LTRL590_350024600.1 (forward primer: GCGAAGCTGAGGCGGGAGAACGAG; reverse primer: AGCTTCTCCGCATCCGCAGCGAG) was used as the internal control. In addition, 13 additional mRNAs were analyzed by RT-qPCR (Supplementary Data 9).

Validation of DTTs by proteomic analysis

Promastigotes were first mixed with extraction buffer (50 mM ammonium bicarbonate (Fisher Scientific) buffer (ABC), 5% sodium deoxycholate (SDC), (Fisher Scientific). The samples were homogenized by adding 3-mm zirconium beads using a beads beater. The homogenization was done at 4 °C at 4000 rpm for 30 s, followed by a 30-s pause. This step was repeated four times. After centrifuging for 15 min, the supernatant containing extracted proteins was collected and diluted ten times by adding 50 mM ABC buffer. Protein concentration was determined by BCA protein assay kit (ThermoScientific Fisher, #23227) following the manufacturer’s instructions.

A 20-µg aliquot of extracted proteins from each sample was then subjected to reduction, alkylation, and tryptic digestion. ABC (50 mM)-SDC (0.5%) solution was first added to samples to maintain a volume of 50 µL. Proteins were thermally denatured at 80 °C for 10 min. The reduction of proteins was performed by adding a 1.25-µL aliquot of 200 mM dithiothreitol solution and incubating at 60 °C for 45 min. The reduced proteins were then alkylated by adding a 5-µL aliquot of 200 mM iodoacetamide solution and incubation at 37 °C in the dark for 45 min. To quench the excessive IAA, a 1.25-µL aliquot of DTT solution was added again and samples were incubated at 37 °C for 30 min. Following the reduction and alkylation of proteins, trypsin (Promega) was added at a ratio of 1:25 (enzyme: proteins, w/w) into samples and incubated at 37 °C for 18 h. After incubation, formic acid was added at a final concentration of 0.5% (v/v) for both quenching the enzymatic reaction and removing the SDC detergent. Samples were then mixed thoroughly and centrifuged at 21,100×g for 10 min. The supernatant was collected, speed-vac dried, and resuspended in an aqueous solution containing 2% acetonitrile (ACN) and 0.1% formic acid (FA) before LC-MS/MS analysis using ThermoScientific XCalibur 4.2 at Proteomics/Metabolomics core facility, Texas Tech University.

The raw data files obtained from LC-MS/MS analysis were processed with the Proteome Discoverer software version 2.2 (Thermo Fisher Scientific). Database search was performed against UniProtKB/Swiss-Prot L. major database. The search included cysteine carbamidomethylation as a fixed modification and variable modifications, including methionine oxidation and acetylation of protein N-terminal. Trypsin was specified as the proteolysis enzyme and a maximum of two missed cleavages were allowed. For identification, the mass tolerance of 10 ppm and ±0.02 Da was chosen for MS and MS/MS analysis, respectively. The false discovery rate (FDR) was set to be 0.01 at both peptide and proteins and only scans with ≤0.01 FDR as determined by Percolator were used for protein identification. Additionally, only proteins with at least two identified peptides were considered for further analysis. Detailed workflow is described in Supplementary Fig. 6.

Analysis of the differentially expressed genes detected at the translatome versus transcriptome level

After DESeq2 analysis, two groups of samples were used to compare the changes in the translatome Vs. transcriptome. The DTTs detected by analyzing heavy polysome fractions were used to indicate changes in highly translated transcripts. The DEG detected in the total mRNA (total RNA before polysome profiling fractionation) was considered to identify changes in the total transcriptome. In both cases, the statistically significant changes were identified by using a cutoff of absolute fold change ≥1.5 and a Benjamini–Hochberg adjusted P value ≤0.05. The resulting differentially translated/expressed genes were matched by gene ID and these genes were classified into 4 groups. Group 1: genes that were detected as differentially expressed only in heavy polysomes (DTT). Group 2: genes that were detected as differentially expressed only in the total transcriptome (transcriptionally regulated or mRNA stability). Group 3: genes that were detected as differentially expressed in both total transcriptome and heavy polysomes (mRNA abundance). Group 4: the remaining genes that did not show significant changes in either analysis. The analysis was performed using three biological replicates and considered both resistant parasites growing without and with drug challenges.

Functional enrichment analysis

Gene ontology analysis was performed using the TriTrypDB web browser⁸⁶. Networks based on protein homology were based on the interactome published in ref. ⁸⁷. Networks based on protein co-expression were built in (https://string-db.org/). All networks were visualized on Cytoscape 3.7.2.

Bioinformatics approach for call variants analysis

The variants were identified utilizing FreeBayes algorithm (Version 1.3.1), which has been successfully used in Leishmania that commonly exhibits chromosome and gene copy number variations^34,88,89,90. VCFfilter (Version 1.0) was employed to filter FreeBayes VCF for strand bias (Strand balance probability for the reference allele or SRP >30 and Strand balance probability for the reference allele or SAP > 30), placement bias (End Placement Probability or EPP > 30), variant quality (QUAL > 30), and depth of coverage (Total read depth at the locus or DP > 100)⁹¹. The.vcf outputs from VCFfilter were annotated with SnpEff eff function using a custom database built with the.GFF dataset using the SnpEff v4.3 build function⁹². The identified variants were visually and manually inspected for accuracy using the Integrative Genomics Viewer (IGV). The variants exclusively detected in HR strain were considered for downstream analysis. Detailed workflow is described in Supplementary Fig. 6. The variants were summarized in a circus plot using the shinyCircos v1 function⁹³.

Statistical analyses

Statistical analyses summarized in bar and box plots were performed using GraphPad Prizm software v9.1.2. Statistical analyses for RNA-Seq and proteomic data were performed using R v4.1 and proteome Discoverer software version 2.2, respectively.

Reporting summary

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

Data availability

The deep RNA-Seq datasets supporting the conclusions of this article are available in the Gene Expression Omnibus database from the National Center for Biotechnology under accession number GSE173848. The proteomics data supporting the findings of this study are available through Mass Spectrometry Interactive Virtual Environment under accession number MSV000089617 (MassIVE, https://massive.ucsd.edu/ProteoSAFe/static/massive.jsp?redirect=auth). The project data are grouped and available in the BioStudies database under accession number S-BSST850 (https://www.ebi.ac.uk/biostudies/). Source data are provided with this paper.

References

WHO. http://www.who.int/mediacentre/factsheets/fs375/en/ (2020).
Torres-Guerrero, E., Quintanilla-Cedillo, M. R., Ruiz-Esmenjaud, J. & Arenas, R. Leishmaniasis: a review. F1000Research 6, 750–750 (2017).
Article PubMed PubMed Central Google Scholar
Sneader, W. Chemical Medicines. In: Drug Discovery, (ed. W. Sneader). John Wiley & Sons Ltd. (2005).
Rojas, R. et al. Resistance to antimony and treatment failure in human Leishmania (Viannia) infection. J. Infect. Dis. 193, 1375–1383 (2006).
Article CAS PubMed Google Scholar
Unger, A. et al. Association of treatment of American cutaneous leishmaniasis prior to ulcer development with high rate of failure in northeastern Brazil. Am. J. Trop. Med. Hyg. 80, 574–579 (2009).
Article PubMed Google Scholar
Sundar, S. et al. Failure of pentavalent antimony in visceral leishmaniasis in India: report from the center of the Indian epidemic. Clin. Infect. Dis. 31, 1104–1107 (2000).
Article CAS PubMed Google Scholar
Eddaikra, N. et al. Antimony susceptibility of Leishmania isolates collected over a 30-year period in Algeria. PLoS Negl. Trop. Dis. 12, e0006310 (2018).
Article PubMed PubMed Central Google Scholar
Llanos-Cuentas, A. et al. Clinical and parasite species risk factors for pentavalent antimonial treatment failure in cutaneous leishmaniasis in Peru. Clin. Infect. Dis. 46, 223–231 (2008).
Article PubMed Google Scholar
Chakravarty, J. & Sundar, S. Drug resistance in leishmaniasis. J. Glob. Infect. Dis. 2, 167–176 (2010).
Article PubMed PubMed Central Google Scholar
Gutierrez Guarnizo, S. A. et al. Drug-induced lipid remodeling in Leishmania parasites. Microorganisms 9, 790 (2021).
Article PubMed PubMed Central Google Scholar
Laffitte, M. N., Leprohon, P., Papadopoulou, B. & Ouellette, M. Plasticity of the Leishmania genome leading to gene copy number variations and drug resistance. F1000Res 5, 2350 (2016).
Article PubMed PubMed Central Google Scholar
Patino, L. H., Muskus, C. & Ramírez, J. D. Transcriptional responses of Leishmania (Leishmania) amazonensis in the presence of trivalent sodium stibogluconate. Parasites Vectors 12, 348 (2019).
Article PubMed PubMed Central Google Scholar
Gutierrez Guarnizo, S. A., Karamysheva, Z. N., Galeano, E. & Muskus, C. E. Metabolite biomarkers of leishmania antimony resistance. Cells 10, https://doi.org/10.3390/cells10051063 (2021).
Maharjan, M., Singh, S., Chatterjee, M. & Madhubala, R. Role of aquaglyceroporin (AQP1) gene and drug uptake in antimony-resistant clinical isolates of Leishmania donovani. Am. J. Trop. Med. Hyg. 79, 69–75 (2008).
Article CAS PubMed Google Scholar
Singh, N., Chatterjee, M. & Sundar, S. The overexpression of genes of thiol metabolism contribute to drug resistance in clinical isolates of visceral leishmaniasis (kala azar) in India. Parasites Vectors 7, 596 (2014).
Article PubMed PubMed Central Google Scholar
Mookerjee Basu, J. et al. Inhibition of ABC transporters abolishes antimony resistance in Leishmania Infection. Antimicrob. Agents Chemother. 52, 1080–1093 (2008).
Article PubMed Google Scholar
Andrade, J. M. et al. Comparative transcriptomic analysis of antimony resistant and susceptible Leishmania infantum lines. Parasites Vectors 13, 600 (2020).
Article CAS PubMed PubMed Central Google Scholar
Das, S. et al. Over-expression of 60s ribosomal L23a is associated with cellular proliferation in SAG resistant clinical isolates of Leishmania donovani. PLOS Negl. Trop. Dis. 7, e2527 (2013).
Article PubMed PubMed Central Google Scholar
Matrangolo, F. S. V. et al. Comparative proteomic analysis of antimony-resistant and -susceptible Leishmania braziliensis and Leishmania infantum chagasi lines. Mol. Biochem. Parasitol. 190, 63–75 (2013).
Article CAS PubMed Google Scholar
Brotherton, M. C. et al. Proteomic and genomic analyses of antimony resistant Leishmania infantum mutant. PLoS ONE 8, e81899 (2013).
Article ADS PubMed PubMed Central Google Scholar
Medina, J., Cruz-Saavedra, L., Patiño, L. H., Muñoz, M. & Ramírez, J. D. Comparative analysis of the transcriptional responses of five Leishmania species to trivalent antimony. Parasites Vectors 14, 419 (2021).
Article CAS PubMed PubMed Central Google Scholar
Cohen-Freue, G., Holzer, T. R., Forney, J. D. & McMaster, W. R. Global gene expression in Leishmania. Int J. Parasitol. 37, 1077–1086 (2007).
Article CAS PubMed Google Scholar
Karamysheva, Z. N., Gutierrez Guarnizo, S. A. & Karamyshev, A. L. Regulation of translation in the protozoan parasite Leishmania. Int. J. Mol. Sci. 21, https://doi.org/10.3390/ijms21082981 (2020).
Piel, L. et al. Experimental evolution links post-transcriptional regulation to Leishmania fitness gain. PLOS Pathogens. 18, e1010375 https://doi.org/10.1371/journal.ppat.1010375 (2022).
Bhattacharya, A. et al. Coupling chemical mutagenesis to next generation sequencing for the identification of drug resistance mutations in Leishmania. Nat. Commun. 10, 5627 (2019).
Article ADS CAS PubMed PubMed Central Google Scholar
Magill, A. J., Grögl, M., Gasser, R. A. Jr., Sun, W. & Oster, C. N. Visceral infection caused by Leishmania tropica in veterans of operation desert storm. N. Engl. J. Med. 328, 1383–1387 (1993).
Article CAS PubMed Google Scholar
Alborzi, A., Rasouli, M. & Shamsizadeh, A. Leishmania tropica-isolated patient with visceral leishmaniasis in southern Iran. Am. J. Trop. Med. Hyg. 74, 306–307 (2006).
Article PubMed Google Scholar
Weiss, F., Vogenthaler, N., Franco-Paredes, C. & Parker, S. R. Leishmania tropica-induced cutaneous and presumptive concomitant viscerotropic leishmaniasis with prolonged incubation. Arch. Dermatol. 145, 1023–1026 (2009).
Article PubMed Google Scholar
Gazola, K. C. et al. Cell surface carbohydrates and in vivo infectivity of glucantime-sensitive and resistant Leishmania (Viannia) guyjanensis cell lines. Parasitol. Res 87, 935–940 (2001).
Article CAS PubMed Google Scholar
Silva, N., Camacho, N., Figarella, K. & Ponte-Sucre, A. Cell differentiation and infectivity of Leishmania mexicana are inhibited in a strain resistant to an ABC-transporter blocker. Parasitology 128, 629–634 (2004).
Article CAS PubMed Google Scholar
Liarte, D. B. & Murta, S. M. Selection and phenotype characterization of potassium antimony tartrate-resistant populations of four New World Leishmania species. Parasitol. Res. 107, 205–212 (2010).
Article PubMed Google Scholar
Marcili, A. et al. Phylogenetic relationships of Leishmania species based on trypanosomatid barcode (SSU rDNA) and gGAPDH genes: taxonomic revision of Leishmania (L.) infantum chagasi in South America. Infect. Genet. Evol. 25, 44–51 (2014).
Article CAS PubMed Google Scholar
Love, M. I., Huber, W. & Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 550–550 (2014).
Article PubMed PubMed Central Google Scholar
Garrison, E. & Marth, G. Haplotype-based variant detection from short-read sequencing. Preprint at https://arxiv.org/abs/1207.3907 (2012).
Tebaldi, T. et al. Widespread uncoupling between transcriptome and translatome variations after a stimulus in mammalian cells. BMC Genomics 13, 220 (2012).
Article CAS PubMed PubMed Central Google Scholar
Clayton, C. E. Life without transcriptional control? From fly to man and back again. EMBO J. 21, 1881–1888 (2002).
Clayton, C. Regulation of gene expression in trypanosomatids: living with polycistronic transcription. Open Biol. 9, 190072 (2019).
Walker, J. et al. Discovery of factors linked to antimony resistance in Leishmania panamensis through differential proteome analysis. Mol. Biochem. Parasitol. 183, 166–176 (2012).
Article CAS PubMed Google Scholar
Perea, A., Manzano, J. I., Castanys, S. & Gamarro, F. The LABCG2 transporter from the protozoan parasite Leishmania is involved in antimony resistance. Antimicrob. Agents Chemother. 60, 3489–3496 (2016).
Article CAS PubMed PubMed Central Google Scholar
El Fadili, K., Drummelsmith, J., Roy, G., Jardim, A. & Ouellette, M. Down regulation of KMP-11 in Leishmania infantum axenic antimony resistant amastigotes as revealed by a proteomic screen. Exp. Parasitol. 123, 51–57 (2009).
Article CAS PubMed Google Scholar
Rastrojo, A. et al. Genomic and transcriptomic alterations in Leishmania donovani lines experimentally resistant to antileishmanial drugs. Int J. Parasitol. Drugs Drug Resist. 8, 246–264 (2018).
Article PubMed PubMed Central Google Scholar
Verma, A. et al. Transcriptome profiling identifies genes/pathways associated with experimental resistance to paromomycin in Leishmania donovani. Int. J. Parasitol. Drugs Drug Resist. 7, 370–377 (2017).
Article PubMed PubMed Central Google Scholar
Patino, L. H. et al. Major changes in chromosomal somy, gene expression and gene dosage driven by SbIII in Leishmania braziliensis and Leishmania panamensis. Sci. Rep. 9, 9485 (2019).
Article ADS PubMed PubMed Central Google Scholar
Downing, T. et al. Whole genome sequencing of multiple Leishmania donovani clinical isolates provides insights into population structure and mechanisms of drug resistance. Genome Res. 21, 2143–2156 (2011).
Article CAS PubMed PubMed Central Google Scholar
van Dam, S., Võsa, U., van der Graaf, A., Franke, L. & de Magalhães, J. P. Gene co-expression analysis for functional classification and gene-disease predictions. Brief. Bioinforma. 19, 575–592 (2018).
Google Scholar
Wallin, G., Kamerlin, S. C. L. & Åqvist, J. Energetics of activation of GTP hydrolysis on the ribosome. Nat. Commun. 4, 1733 (2013).
Article ADS PubMed Google Scholar
Kolupaeva, V. Serine-threonine protein phosphatases: lost in translation. Biochim Biophys. Acta Mol. Cell Res. 1866, 83–89 (2019).
Article CAS PubMed Google Scholar
Martin, R., Straub, A., Doebele, C. & Bohnsack, M. DExD/H-box RNA helicases in ribosome biogenesis. RNA Biol. 9, https://doi.org/10.4161/rna.21879 (2012).
Bourgeois, C. F., Mortreux, F. & Auboeuf, D. The multiple functions of RNA helicases as drivers and regulators of gene expression. Nat. Rev. Mol. Cell Biol. 17, 426–438 (2016).
Article CAS PubMed Google Scholar
Günzl, A. The pre-mRNA splicing machinery of trypanosomes: complex or simplified? Eukaryot. Cell 9, 1159–1170 (2010).
Article PubMed PubMed Central Google Scholar
Nare, B., Garraway, L. A., Vickers, T. J. & Beverley, S. M. PTR1-dependent synthesis of tetrahydrobiopterin contributes to oxidant susceptibility in the trypanosomatid protozoan parasite Leishmania major. Curr. Genet. 55, 287–299 (2009).
Article CAS PubMed PubMed Central Google Scholar
Görlach, A., Bertram, K., Hudecova, S. & Krizanova, O. Calcium and ROS: a mutual interplay. Redox Biol. 6, 260–271 (2015).
Article PubMed PubMed Central Google Scholar
Bapi, B., Bhakti, L., Pratyay, S., Koushik, D. & Monidipa, G. Characterization of a novel amastin-like surface protein (ALSP) of Leishmania donovani, a probable lipase. Preprint at bioRxiv https://doi.org/10.1101/2020.07.23.218107 (2021).
Monroe, J. G. et al. Mutation bias reflects natural selection in Arabidopsis thaliana. Nature 602, 101–105 (2022).
Article ADS CAS PubMed PubMed Central Google Scholar
Bhargava, P. & Singh, R. Developments in diagnosis and antileishmanial drugs. Interdiscip. Perspect. Infect. Dis. 2012, 626838 (2012).
Article PubMed PubMed Central Google Scholar
Purkait, B. et al. Mechanism of amphotericin B resistance in clinical isolates of Leishmania donovani. Antimicrob. Agents Chemother. 56, 1031–1041 (2012).
Article CAS PubMed PubMed Central Google Scholar
Coelho, A. C., Beverley, S. M. & Cotrim, P. C. Functional genetic identification of PRP1, an ABC transporter superfamily member conferring pentamidine resistance in Leishmania major. Mol. Biochem. Parasitol. 130, 83–90 (2003).
Article CAS PubMed Google Scholar
Dumetz, F. et al. Molecular preadaptation to antimony resistance in Leishmania donovani on the Indian subcontinent. mSphere 3, https://doi.org/10.1128/mSphere.00548-17 (2018).
Kramer, S. Developmental regulation of gene expression in the absence of transcriptional control: the case of kinetoplastids. Mol. Biochem. Parasitol. 181, 61–72 (2012).
Article CAS PubMed Google Scholar
Vasquez, J. J., Hon, C. C., Vanselow, J. T., Schlosser, A. & Siegel, T. N. Comparative ribosome profiling reveals extensive translational complexity in different Trypanosoma brucei life cycle stages. Nucleic Acids Res. 42, 3623–3637 (2014).
Article CAS PubMed PubMed Central Google Scholar
Smircich, P. et al. Ribosome profiling reveals translation control as a key mechanism generating differential gene expression in Trypanosoma cruzi. BMC Genomics 16, 443 (2015).
Article PubMed PubMed Central Google Scholar
Bifeld, E. et al. Ribosome profiling reveals HSP90 inhibitor effects on stage-specific protein synthesis in Leishmania donovani. mSystems 3, e00214–e00218 (2018).
Article CAS PubMed PubMed Central Google Scholar
Bajak, K. & Clayton, C. Polysome profiling and metabolic labeling methods to measure translation in Trypanosoma brucei. Methods Mol. Biol. 2116, 99–108 (2020).
Article CAS PubMed Google Scholar
Lee, L. J. et al. Cancer plasticity: the role of mRNA translation. Trends Cancer 7, 134–145 (2021).
Article CAS PubMed Google Scholar
Lelong, E. I. J. et al. Prostate cancer resistance leads to a global deregulation of translation factors and unconventional translation. NAR Cancer. 4, https://doi.org/10.1093/narcan/zcac034 (2022).
McMahon, M. & Ruggero, D. A wobbly road to drug resistance in melanoma: tRNA-modifying enzymes in translation reprogramming. EMBO J. 37, e99978 (2018).
Mukhopadhyay, R. et al. Trypanothione overproduction and resistance to antimonials and arsenicals in Leishmania. Proc. Natl Acad. Sci. USA 93, 10383–10387 (1996).
Article ADS CAS PubMed PubMed Central Google Scholar
Mohebali, M. et al. Gene expression analysis of antimony resistance in Leishmania tropica using quantitative real-time PCR focused on genes involved in trypanothione metabolism and drug transport. Arch. Dermatol Res. 311, 9–17 (2019).
Article CAS PubMed Google Scholar
Légaré, D. et al. The Leishmania ATP-binding cassette protein PGPA is an intracellular metal-thiol transporter ATPase. J. Biol. Chem. 276, 26301–26307 (2001).
Article PubMed Google Scholar
Leprohon, P., Légaré, D., Girard, I., Papadopoulou, B. & Ouellette, M. Modulation of Leishmania ABC protein gene expression through life stages and among drug-resistant parasites. Eukaryot. cell 5, 1713–1725 (2006).
Article CAS PubMed PubMed Central Google Scholar
Cloutier, S. et al. Translational control through eIF2alpha phosphorylation during the Leishmania differentiation process. PLoS ONE 7, e35085 (2012).
Article ADS CAS PubMed PubMed Central Google Scholar
Bertorello, J. et al. Translation reprogramming by eIF3 linked to glioblastoma resistance. NAR Cancer 2, https://doi.org/10.1093/narcan/zcaa020 (2020).
Gong, C. et al. Phosphorylation independent eIF4E translational reprogramming of selective mRNAs determines tamoxifen resistance in breast cancer. Oncogene 39, 3206–3217 (2020).
Article CAS PubMed PubMed Central Google Scholar
Shaked-Mishan, P., Ulrich, N., Ephros, M. & Zilberstein, D. Novel Intracellular SbV reducing activity correlates with antimony susceptibility in Leishmania donovani. J. Biol. Chem. 276, 3971–3976 (2001).
Article CAS PubMed Google Scholar
Riss, T. L., Moravec, R. A. & Niles, A. L., et al. Cell Viability Assays. In: Assay Guidance Manual [Internet]. National Library of Medicine (eds Markossian, S. et al.) (2013).
Palacios, G. et al. Studies in vitro on infectivity and sensitivity to antileishmanial drugs in New World Leishmania species transfected with the green fluorescent protein [pIR3(-)-eGFP]. Parasitology 144, 1718–1725 (2017).
Article CAS PubMed Google Scholar
Karamysheva, Z. N. et al. Polysome profiling in Leishmania, human cells and mouse testis. J. Vis. Exp. https://doi.org/10.3791/57600 (2018).
Afgan, E. et al. The Galaxy platform for accessible, reproducible and collaborative biomedical analyses: 2018 update. Nucleic Acids Res. 46, W537–W544 (2018).
Article CAS PubMed PubMed Central Google Scholar
Andrews, S. Babraham bioinformatics—FastQC a quality control tool for high throughput sequence data. https://www.bioinformatics.babraham.ac.uk/projects/fastqc/ (2015).
Krueger, F. Brahan bioinformatics trim galore. https://www.bioinformatics.babraham.ac.uk/projects/trim_galore/ (2015).
Ewels, P., Magnusson, M., Lundin, S. & Käller, M. MultiQC: summarize analysis results for multiple tools and samples in a single report. Bioinformatics 32, 3047–3048 (2016).
Article CAS PubMed PubMed Central Google Scholar
Li, H. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. Preprint at https://arxiv.org/abs/1303.3997 (2013).
Picard toolkit. Broad Institute. Available at: https://broadinstitute.github.io/picard/. Accessed: 2nd May 2023. Broad Institute, GitHub repository (2019).
Liao, Y., Smyth, G. K. & Shi, W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 30, 923–930 (2014).
Article CAS PubMed Google Scholar
Schmittgen, T. D. & Livak, K. J. Analyzing real-time PCR data by the comparative C(T) method. Nat. Protoc. 3, 1101–1108 (2008).
Article CAS PubMed Google Scholar
Aslett, M. et al. TriTrypDB: a functional genomic resource for the Trypanosomatidae. Nucleic Acids Res. 38, D457–D462 (2010).
Article CAS PubMed Google Scholar
Flórez, A. F. et al. Protein network prediction and topological analysis in Leishmania major as a tool for drug target selection. BMC Bioinforma. 11, 484 (2010).
Article Google Scholar
Wang, J. et al. Investigation of rare and low-frequency variants using high-throughput sequencing with pooled DNA samples. Sci. Rep. 6, 33256 (2016).
Article ADS CAS PubMed PubMed Central Google Scholar
Camacho, E. et al. Complete assembly of the Leishmania donovani (HU3 strain) genome and transcriptome annotation. Sci. Rep. 9, 6127 (2019).
Article ADS PubMed PubMed Central Google Scholar
Rogers, M. B. et al. Chromosome and gene copy number variation allow major structural change between species and strains of Leishmania. Genome Res. 21, 2129–2142 (2011).
Article CAS PubMed PubMed Central Google Scholar
Olson, N. D. et al. Best practices for evaluating single nucleotide variant calling methods for microbial genomics. Front. Genet. 6, 235–235 (2015).
Article PubMed PubMed Central Google Scholar
Cingolani, P. et al. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly 6, 80–92 (2012).
Article CAS PubMed PubMed Central Google Scholar
Yu, Y., Ouyang, Y. & Yao, W. shinyCircos: an R/Shiny application for interactive creation of Circos plot. Bioinformatics 34, 1229–1231 (2018).
Article CAS PubMed Google Scholar

Download references

Acknowledgements

We thank all members of the lab for helpful discussion and support, especially to Sarah Hernandez for the critical reading of the manuscript. We are thankful to Admera Health for deep RNA-Seq and core facilities at Texas Tech University for proteomics services. We acknowledge to the Galaxy Project (https://usegalaxy.org) for the web resource, computational server, and ToolShed repository used in this study for bioinformatic analyses. This study was supported by funds from “Minciencias and Gobernación del Tolima-Colombia as part of the program 755–2016” to S.A.GG and C.E.M; the Start-up fund from Texas Tech University Health Sciences Center to A.L.K.; and in part by the National Institute Of General Medical Sciences of the National Institutes of Health under Award Number R15GM146171 to Z.N.K. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.

Author information

These authors contributed equally: Andrey L. Karamyshev, Carlos E. Muskus, Zemfira N. Karamysheva.

Authors and Affiliations

Department of Cell Biology and Biochemistry, Texas Tech University Health Sciences Center, Lubbock, TX, 79430, USA
Sneider Alexander Gutierrez Guarnizo, Elena B. Tikhonova & Andrey L. Karamyshev
Programa de Estudio y Control de Enfermedades Tropicales, Universidad de Antioquia. Medellín, Medellín, 050010, Colombia
Sneider Alexander Gutierrez Guarnizo & Carlos E. Muskus
Department of Biological Sciences, Texas Tech University, Lubbock, TX, 79409, USA
Zemfira N. Karamysheva

Authors

Sneider Alexander Gutierrez Guarnizo
View author publications
You can also search for this author in PubMed Google Scholar
Elena B. Tikhonova
View author publications
You can also search for this author in PubMed Google Scholar
Andrey L. Karamyshev
View author publications
You can also search for this author in PubMed Google Scholar
Carlos E. Muskus
View author publications
You can also search for this author in PubMed Google Scholar
Zemfira N. Karamysheva
View author publications
You can also search for this author in PubMed Google Scholar

Contributions

Z.N.K., C.M., and S.A.G.G. designed experiments and analyzed the results. S.A.G.G., Z.N.K., E.B.T., and A.L.K. performed experiments. S.A.G.G. performed bioinformatics analysis. S.A.G.G. and Z.N.K. wrote the manuscript. All the authors reviewed the final manuscript version.

Corresponding authors

Correspondence to Andrey L. Karamyshev, Carlos E. Muskus or Zemfira N. Karamysheva.

Ethics declarations

Competing interests

The authors declare no competing interests.

Peer review

Peer review information

Nature Communications thanks the anonymous reviewer(s) for their contribution to the peer review of this work. A peer review file is available.

Additional information

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Supplementary Information file

Peer Review File

Description of additional supplementary files

Supplementary Data 1

Supplementary Data 2

Supplementary Data 3

Supplementary Data 4

Supplementary Data 5

Supplementary Data 6

Supplementary Data 7

Supplementary Data 8

Supplementary Data 9

Supplementary Data 10

Supplementary Data 11

Supplementary Data 12

Supplementary Data 13

Supplementary Data 14

Supplementary Data 15

Supplementary Data 16

Supplementary Data 17

Supplementary Data 18

Reporting Summary

Source data

Source Data

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Cite this article

Gutierrez Guarnizo, S.A., Tikhonova, E.B., Karamyshev, A.L. et al. Translational reprogramming as a driver of antimony-drug resistance in Leishmania. Nat Commun 14, 2605 (2023). https://doi.org/10.1038/s41467-023-38221-1

Download citation

Received: 25 May 2022
Accepted: 20 April 2023
Published: 05 May 2023
DOI: https://doi.org/10.1038/s41467-023-38221-1

Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.