Next Article in Journal
Plant Acyl-CoA-Binding Proteins—Their Lipid and Protein Interactors in Abiotic and Biotic Stresses
Previous Article in Journal
Weapons Evolve Faster Than Sperm in Bovids and Cervids
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cells
Volume 10
Issue 5
10.3390/cells10051063
cells-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Metabolite Biomarkers of Leishmania Antimony Resistance

by
Sneider Alexander Gutierrez Guarnizo
1,
Zemfira N. Karamysheva
2,
Elkin Galeano
3 and
Carlos E. Muskus
1,*
1
Programa de Estudio y Control de Enfermedades Tropicales, Facultad de Medicina, Universidad de Antioquia, Medellín 050010, Colombia
2
Department of Biological Sciences, Texas Tech University, Lubbock, TX 79409, USA
3
Grupo de Investigación en Sustancias Bioactivas-GISB, Universidad de Antioquia, Medellín 050010, Colombia
*
Author to whom correspondence should be addressed.
Cells 2021, 10(5), 1063; https://0-doi-org.brum.beds.ac.uk/10.3390/cells10051063
Submission received: 8 April 2021 / Revised: 20 April 2021 / Accepted: 22 April 2021 / Published: 30 April 2021
Browse Figures
Review Reports Versions Notes

Abstract

:
Leishmania parasites cause leishmaniasis, one of the most epidemiologically important neglected tropical diseases. Leishmania exhibits a high ability of developing drug resistance, and drug resistance is one of the main threats to public health, as it is associated with increased incidence, mortality, and healthcare costs. The antimonial drug is the main historically implemented drug for leishmaniasis. Nevertheless, even though antimony resistance has been widely documented, the mechanisms involved are not completely understood. In this study, we aimed to identify potential metabolite biomarkers of antimony resistance that could improve leishmaniasis treatment. Here, using L. tropica promastigotes as the biological model, we showed that the level of response to antimony can be potentially predicted using 1H-NMR-based metabolomic profiling. Antimony-resistant parasites exhibited differences in metabolite composition at the intracellular and extracellular levels, suggesting that a metabolic remodeling is required to combat the drug. Simple and time-saving exometabolomic analysis can be efficiently used for the differentiation of sensitive and resistant parasites. Our findings suggest that changes in metabolite composition are associated with an optimized response to the osmotic/oxidative stress and a rearrangement of carbon-energy metabolism. The activation of energy metabolism can be linked to the high energy requirement during the antioxidant stress response. We also found that metabolites such as proline and lactate change linearly with the level of resistance to antimony, showing a close relationship with the parasite’s efficiency of drug resistance. A list of potential metabolite biomarkers is described and discussed.

1. Introduction

Leishmaniasis is a group of vector-borne tropical diseases considered to be one of the most epidemiologically important as it affects approximately 2 million people per year [1,2]. Leishmaniasis is caused by Leishmania parasites, a genus of kinetoplastids that exhibit two different stages during their life cycle: the promastigote stage adapted to survive in the gut of phlebotomine sandfly vectors, and the amastigote stage, evolutionarily adapted to survive inside the vertebrate host’s macrophages [3].
Once Leishmania is established in the vertebrate host, the clinical outcome is the result of the integration of multiple factors including the Leishmania species and/or strain, the host’s immune response, and the presence of coinfections. Consequently, the disease can be present as a wide spectrum of clinical outcomes grouped in four clinical forms: simple or diffuse cutaneous leishmaniasis (CL) characterized by local or diffuse skin ulcers; mucocutaneous leishmaniasis (ML) affecting mainly nasal mucosa; visceral leishmaniasis (VL), a systemic infection affecting organs containing macrophages such as the bone marrow, liver, spleen, and lymph nodes; and Post-kala-azar dermal leishmaniasis (PKDL) considered a dermal sequela of VL. Notably, the VL can be lethal in the absence of treatment [4].
For the last seven decades, pentavalent antimonial (SbV) has been used as the primary leishmaniasis therapy. It is commonly accepted that SbV acts as a prodrug, requiring reduction to the active trivalent antimony (SbIII). Evidence shows that SbIII induces oxidative and osmotic stress [5,6], inhibits the glycolytic pathway and fatty acid β-oxidation [7], interferes with the purine salvage pathway [8], inhibits DNA topoisomerase I [9], and competes with zinc (ZnII) for its binding to the CCHC and CCCH zinc finger domains [10].
Unfortunately, Leishmania has developed resistance to antimonials, leading to over sixty percent inefficacy in Bihar, India [11], a phenomenon that has also been seen in North Africa and Latin America [12,13,14,15,16,17]. Treatment failures due to drug resistant organisms are very common for L. donovani, L. braziliensis/L. guyanensis, L. tropica, and L. major, species responsible for visceral and cutaneous leishmaniasis. Hence, it is important for the clinician to know the isolate’s drug resistance status in order to make the correct therapeutic choice. Other drugs such as paromomycin, amphotericin b, and miltefosine, have been employed, but Leishmania has also shown clinical resistance against all of them [18,19,20]. Furthermore, there is an evidence of cross-resistance between antimonial and other drugs such as amphotericin b and paromomycin [19,21]. As a result, antimonials are still the main therapeutic option in several endemic countries including Colombia, Venezuela, and Brazil [22]. Together, these facts suggest that drug resistance can be considered one of the main challenges to efficiently combat leishmaniasis, and antimony resistance is a phenomenon with high clinical impact that needs to be further studied. However, no molecular methods are currently validated to track drug-resistant isolates in clinical setting. Therefore, there is an urgent need to identify biomarkers suitable to use in clinical settings and develop standardized molecular methods to measure drug resistance in clinical isolates [23].
“Omics” technologies constitute a powerful approach to analyze the molecular remodeling involved in different biological processes including drug resistance by offering the opportunity to study massively, qualitatively, and quantitatively different types of biomolecules such as genes, transcripts, proteins, lipids, or metabolites at a global cell level [24].
Studies of resistance to antimony using omics technologies have mostly focused on the genomic, transcriptomic, and proteomic levels [25,26]. Thus far, the evidence suggests that SbIII-resistant parasites combat antimony in several main ways: (1) decreasing drug uptake by downregulation of the aquaglyceroporin (AQP1) transporter [27]; (2) inhibiting drug by formation of inactive thiol–SbIII complex and activating the thiol metabolism leading to the production of trypanothione, an important molecule involved in antioxidant response [28]; (3) increasing drug efflux by overexpression of ABC transporters that transport the thiol–metal complex toward vesicles that are secreted via exocytosis [29,30,31]; and finally (4) remodeling of carbon and lipid metabolism, which has been consistently documented in SbIII-resistant parasites, probably to optimize energy metabolism and modify the cell membrane composition [26,32,33,34].
However, few studies address the metabolome of Leishmania’s resistant phenotypes. Metabolomics is one of the latest omics technologies that has been applied successfully in many areas of life sciences, being particularly useful for phenotypic analysis since metabolites are the endpoints of the active pathways and therefore they are optimal indicators of the cell’s physiology [35]. Because of this, metabolites could potentially be used as biomarkers to differentiate antimony-resistant parasites.
Furthermore, studies have mainly focused on intracellular metabolism rather than changes at the extracellular level [36,37,38]. Interestingly, the exometabolome or metabolic footprinting is not only informative of the metabolic changes in response to different environmental conditions but also it requires an easier sample preparation than intracellular metabolomic analysis, offering the opportunity to simplify the phenotypic characterization [39].
Mass spectrometry analysis (MS) has shown that antimony-sensitive and -resistant clinical isolates exhibit dramatic differences in metabolomic profiles, particularly in lipid metabolism [40], sulfur-containing amino acids, and polyamine biosynthetic pathways [36]. 13C-isotope-labeling and MS confirmed that the amino acid metabolism is remodeled by involving activation of the redox system in antimony-resistant parasites [37].
However, due to Leishmania’s complexity, more studies are required to identify robust and reliable biomarkers of SbIII resistance. For instance, it remains unexplored as to whether the abundance of differentially expressed metabolites change proportionally to the level of SbIII resistance. Indeed, these types of compounds could be more closely related with the response required to combat the drug and more informative in predicting the level of resistance to SbIII.
Additionally, metabolomics is based on two analytical methods, MS and nuclear magnetic resonance (NMR), which are considered complementary techniques with different advantages and disadvantages; yet, in the Leishmania antimony resistance field, most of the studies have used MS. Although 1H-NMR is less sensitive, it is more robust and reproducible, the sample preparation can be easier, and the sample is not lost during the spectra acquisition, making this a valuable technique for the search of biomarkers [41,42].
This study aimed to identify the metabolites that can be potentially used as biomarkers of the SbIII resistance level across 1H-NMR analyses from both intracellular and extracellular extracts of Leishmania promastigotes.

2. Materials and Methods

2.1. Reagents

Phosphate buffered saline (Santa Cruz Biotechnology, Dallas, TX, USA, product: 362182), Schneider insect medium (Sigma Aldrich, St. Louis, MI, USA, product: S0146), potassium antimony (III) tartrate trihydrate (Sigma Aldrich, St. Louis, MI, USA, product: 28300-74-5), fetal bovine serum (Eurobio Scientific, Paris, France, product: CVFSVF0001), penicillin and streptomycin (Thermo Fischer Scientific, Waltham, MA, USA, product: 15140122), deuterium oxide (Merck, Darmstadt, Hesse, Germany, product: 7789-20-0), sodium azide (NaN3) (Merck, Darmstadt, Hesse, Germany, product: 26628-22-8), 3-(trimethylsilyl) propionic-2,2,3,3-d4 acid sodium salt or TSP (Sigma Aldrich, St. Louis, MI, USA, product: 24493-21-8), and 2′,7′-dichlorodihydrofluorescein diacetate or H2DCFDA (Thermo Fischer Scientific, Waltham, MA, USA, product: D399).

2.2. Parasites, Culture, and Drug Treatment

Three levels of SbIII-resistance were considered for the metabolomic analysis by 1H-NMR. The Leishmania tropica wild type or antimony sensitive strain (WT), moderately resistant (MR) isogenic derivative of WT, and highly resistant isogenic derivative (HR) were used in this study. The WT, MR, and HR strains have a growing half-maximal effective concentration (EC50) to SbIII of 10.4 ± 0.6, 377 ± 34.93, and 631.7 ± 73.7 µg/mL, respectively, and thus they were chosen to represent a low, moderate, and high level of SbIII resistance, respectively.
Both sensitive and resistant strains were independently grown in culture medium. In each case, the promastigotes were seeded at 1 × 106 promastigotes/mL in a final volume of 10 mL of Schneider insect medium supplemented with 10% of fetal bovine serum (FBS), 10 units/mL penicillin, and 0.1 mg/mL of streptomycin. The MR and HR strains were treated with the drug at 5 and 10 times the estimated EC50 value for WT strain (EC50 for WT strain = 10 µg/mL SbIII), respectively. These drug concentrations were used since they correspond to the drug concentration used to select each strain during the drug resistance stepwise selection. The WT strain was grown in drug-free medium.
Growth was monitored by evaluating cell count using a Neubauer chamber to compare the growth phases per strain in order to synchronize the cell cultures and to avoid technical variability.

2.3. Extraction of Intracellular Metabolites

During the stationary phase, the promastigotes were diluted at 15–20 × 106 parasites/mL using Schneider insect medium without fetal bovine serum. To isolate truly resistant parasite population and avoid leakage of metabolites, we enriched the population of live parasites or parasites without membrane damage by centrifugation of 3 mL of Ficoll-Histopaque®-1077 and 6 mL of diluted parasite’s solution in 15 mL conical tubes at 1800× g at room temperature for 10 min. After Ficoll gradient centrifugation, the parasite population without membrane damage was enriched between middle and upper phase of Ficoll gradient, and 2 mL of these phases were carefully recovered in a 15 mL conical tube (Figure S1).
The parasite’s metabolism was quenched by incubating the tube containing the recovered parasites in −80 °C ethanol until they reached 4 °C (≈10 s). After that, any manipulation was performed at 4 °C to avoid sample degradation [43].
The quenched parasites were harvested and washed twice in 1 mL of cooled PBS-1X at 1800× g, 4 °C for 5 min. The number of parasites between different samples was adjusted to a total of 1 × 107 promastigotes, and the final pellet was obtained by centrifugation at 1800× g, 4 °C for 5 min, and immediately resuspended in 200 µL of methanol (−20 °C) for metabolite extraction. To homogeneously release the intracellular metabolites, we treated the parasites with 5 cycles of 1 min freezing in liquid nitrogen and thawed them at room temperature with 20 s of 3400 rpm vortex (Labnet vortex mixer 120 V).
The tubes containing the lysate were shaken overnight at 250 rpm and 4 °C. Then, the samples were incubated at −80 °C for 20 min and centrifuged at 20,000× g, 4 °C, for 30 min to precipitate proteins. Finally, the samples were dried using a SpeedVac concentrator (Thermo Fisher Scientific, Waltham, MA, USA, SPD111V P1) and used for 1H-NMR analysis [44].

2.4. Extraction of Extracellular Metabolites

During the stationary phase, 5 mL of medium was collected and filtered using a 0.22 µm acrodisc syringe filter, previously washed with PBS (10 times) to remove membrane additives. The collected medium was mixed with 2 parts of methanol and incubated overnight at −80 °C for protein precipitation. Then, the samples were centrifuged at 1600× g for 30 min, and each supernatant was recovered for lyophilization using a VirTis Bench Top 4K SP scientific followed by 1H-NMR analysis [44].

2.5. Proton Nuclear Magnetic Resonance Spectroscopy 1H-NMR Analysis

The dried samples were reconstituted in 550 μL of deuterated phosphate buffer (50 mM NaH2PO4, 0.02% NaN3, and 0.04 mM or 1.5 mM TSP for intracellular or extracellular extracts, respectively, pH = 7.0 ± 0.04, and D2O).
After this the samples were vortexed for 30 s and centrifuged at ≈16,000× g at 4 °C for 15 min to remove any insoluble material. Supernatants were then transferred to NMR tubes with 5 mm of diameter and 7 inches of length (14219-032, VWR).
High resolution one-dimensional 1H-NMR spectra were obtained on a 600 MHz Bruker Avance III spectrometer (Bruker BioSpin Ltd., Billerica, MA, USA), coupled to a 5 mm Prodigy® TCI cryoprobe. Applying a standard Bruker 1D spectroscopy pre-saturation pulse sequence (noesypr1d) with optimal water suppression and 256 and 128 scans for intracellular and extracellular metabolome analysis, respectively [45]. D2O was used for the internal lock signal and TSP-d4 as the internal standard with a chemical shift of δ 0.0 at 300 K. All raw 1H-NMR spectra are available under request. The samples were randomized before analysis to avoid progressive bias.

2.6. Data Processing

The 1H-NMR data processing was adapted from Mickiewicz et al. (2019) [46]. The spectra were processed in 2 different ways. First, the spectra were annotated to reveal the type of metabolites that can be identified at the intracellular and extracellular level in L. tropica’s methanolic extracts using Chenomx software 8.5 (www.chenomx.com, accessed on 23 April 2021, Edmonton, AB, Canada) in combination with the public database, human metabolite database (HMDB, www.hmdb.ca; accessed on 13 March 2021).
Then, the spectra were processed in parallel using the software MestReNova 14.2. This included the referencing to the TSP peak at 0.0 ppm, apodization along t1 with an exponential of 0.5 Hz, phasing by Metabonomics algorithm, and baseline correction applying the Whittaker Smoother algorithm. The phasing and baseline quality were verified manually. Finally, the spectra were binned using the method average sum with a width of each integral region of 0.04 ppm, from 0 to 10 ppm. Bins covering signals that are not of interest, such as water and residual methanol signals, were manually removed before statistical analysis.

2.7. Statistical Analysis

The datasets generated after spectra binning were normalized by probabilistic quotient normalization (PQN) [47] and imported to SIMCA Umetrics V14.1 software (Umetrics, Umeå, Sweden) for multivariate statistical analysis, adapting the strategy used by Mickiewicz et al. (2018) [48]. First, the datasets were scaled by the Pareto method [49]. Then, principal component analysis (PCA) and orthogonal partial least squares discriminant analysis (OPLS-DA) were completed as unsupervised and supervised multivariate methods, respectively.
The PCA was carried out to summarize the source of variation in each dataset and explore the sample grouping and the presence of outliers. The OPLS-DA method was used to extract maximum information on discriminant signals (compounds) from the spectra.
The OPLS-DA models were validated by sevenfold cross-validation (CV), calculating the R2Y (the percentage of variation explained by the model), Q2 (the predictive ability of the model), and CV-ANOVA (cross-validated analysis of variance) p-value with a cutoff of 0.05 [50,51]. The significant metabolites were selected from the OPLS-DA regression coefficients and variable influence on projection (VIP) scores higher than 1 using the Jackknife technique as bias estimator (p-value ≤ 0.05).
S-line diagrams were performed to summarize the contribution of each VIP to the OPLS-DA models. S-line plot visualizes the signals (chemical shift) distributed by p(ctr) [1] loading and colored according to the absolute value of the correlation loading, p(corr) [1], or modeled correlation. The top end of the color scale visualizes the NMR shifts that influence the separation of the groups [52].
To identify metabolites linearly affected by the antimony resistance level, which can be more precise and reliable for the prediction of resistant phenotypes, we calculated Pearson correlation scores and filtered them by false discovered rate (FDR ≤ 0.05). The metabolites with a higher absolute Pearson coefficient were prioritized and analyzed using an open source receiver operating curve characteristics (ROC) analysis tool in MetaboAnalyst 5.0 [53].

2.8. Measurement of Reactive Oxygen Species (ROS) Levels

To measure the ROS production in antimony-sensitive and -resistant parasites, we used hydrogen peroxide (H2O2) as an oxidative stress inducer, adapting the methodology used by Karampetsou et al. (2019) [54]. Parasites were harvested by centrifugation at 1800× g for 5 min and washed twice in 1 mL of PBS 1X. Then, the cell density was normalized at 1 × 106 parasites/mL. In each case, a group of parasites was incubated in PBS without any additional treatment (control group). On the other hand, an additional group was treated with 0.5 mM hydrogen peroxide (H2O2) for 4 h in a final volume of 1 mL. Once the period of exposure was completed, the pelleted parasites were resuspended in 100 µL of H2DCFDA at 0.4 µg/mL. H2DCFDA is an indicator for reactive oxygen species (ROS) production in live cells. The samples were incubated for 1 h at room temperature in darkness. Then, the parasites were washed in 1 mL of PBS and resuspended in 300 µL of fresh PBS. Staining with propidium iodide (0.125 µg/mL) was used as a cell viability control. Finally, the fluorescence of PI and H2DCFDA were acquired using 610/20 and 530/30 as detection filters, respectively, in a LSRFortessa™ Becton Dickinson cell flow cytometer.

3. Results

3.1. Experimental Design

Several reports demonstrated that the most significant changes in the metabolite’s composition occur during stationary phase [40,55]. Therefore, we collected samples from this stage of parasite growth. To detect metabolomic changes in SbIII-resistant parasites, we analyzed both intracellular and extracellular extracts of Leishmania’s promastigote in at least three independent biological replicates. The intracellular metabolome reflects the cell physiology, while the extracellular metabolome represents the metabolites uptake and excretion in a particular condition.
In this study, we analyzed three strains (WT, MR, and HR) representing increasing SbIII resistance levels because of a progressive stepwise resistance selection in vitro. Since the strains show different phenotypes (SbIII resistance levels), they were considered independent for statistical purposes.
The two resistant strains (MR and HR) were treated with SbIII to stimulate their acquired mechanism of resistance, while the WT was grown without drug treatment and used as a control since the sensitive parasite cannot grow until stationary phase in the presence of the drug. Under the described experimental conditions, all strains showed a similar growth curve, eliminating growth stage bias and allowing for the comparison of 1H-NMR profiling between groups (Figure S2).
Typically, a 1H-NMR spectrum shows the proton signals distribution based on their chemical shift (a proton with a particular chemical environment) and their intensity (abundance), offering the opportunity to characterize the type of compound (fingerprint or signals pattern) and estimate its concentration (signal intensity). On this basis, the workflow represented in Figure 1 included sample preparation followed by 1H-HMR analysis. Spectra annotation was performed assigning the corresponding compound to each detected signal in order to define the Leishmania’s metabolome coverage detected by 1H-NMR. After spectra binning, the multivariable analysis was used to identify the group of signals (bins) significantly differentiating SbIII-resistant phenotypes or VIPs. The signals detected as VIPs were matched with the previously annotated spectra, and the respective compounds were suggested as potential biomarkers of SbIII resistance phenotypes, prioritizing those that exhibited a significant linear correlation between metabolite concentration (signal intensity) and SbIII resistance level. Finally, we interactively analyzed the intracellular and extracellular metabolome to discuss the metabolic pathways potentially affected (Figure 1).

3.2. Metabolomic Coverage in Intracellular and Extracellular Methanolic Extracts Using 1H-NMR

After 1H-NMR analysis, both intracellular and extracellular spectra were manually annotated to define the metabolome coverage. Usually, the metabolome coverage is limited by several factors including the metabolite concentration and the solvent. In our study here, we employed a methanol solvent for extraction since it is recognized as having a better metabolite coverage [56].
Excluding the reference compound and the solvent’s signals, we annotated 26 and 29 metabolites in the spectra from intracellular and extracellular extracts, respectively. The same compounds were detected in both SbIII-sensitive and -resistant parasites (Figure 2 and Figure 3, Tables S1 and S2).
Considering that some metabolites are commonly detected at the intracellular and extracellular level, a total of 40 different compounds were detected by 1H-NMR, 15 of which were commonly detected in both intracellular and extracellular extracts. A total of 11 and 14 compounds were exclusively identified at the intracellular or extracellular levels, respectively (Table S3).
The type of metabolites detected were mainly distributed in three metabolite classes: amino acids (total matching compounds: 18), tricarboxylic acids (total matching compounds: 3), and purines (total matching compounds: 2) (Table S4). Additionally, these compounds were enriched significantly in six pathways of the KEGG database, including aminoacyl-tRNA biosynthesis (total hits: 16); valine, leucine, and isoleucine biosynthesis (total hits: 4); arginine biosynthesis (total hits: 4); alanine, aspartate, and glutamate metabolism (total hits: 5); glyoxylate and dicarboxylate metabolism (total hits: 5); pantothenate and CoA biosynthesis (total hits: 2); nicotinate; and nicotinamide metabolism (total hits: 3) (Table S5).

3.3. 1H-NMR Spectra Efficiently Differentiated the Antimony-Resistant Phenotypes Both at the Intracellular and Extracellular Levels

The PCA models were based on two principal components that in total summarized 74.2% and 43.2% of the variance in the dataset for intracellular and extracellular analysis, respectively (Figure 4A and Figure 5A). PCA score plots showed that resistant strains were clustered away from the sensitive strain, suggesting that the metabolomic profiles were different in both intracellular and extracellular extracts. Additionally, no outliers were detected since all samples were distributed within the ellipse, representing Hotelling’s T-squared 95% confidence interval (Figure 4A and Figure 5A).
The OPLS-DA analysis was carried out to reveal metabolic differences between SbIII-sensitive and SbIII-resistant parasites. Models were obtained using one predictive and one orthogonal component. The OPLS-DA score scatter plots showed a clear separation between SbIII-resistant (MR and HR strain) and SbIII-sensitive phenotypes (WT) at both the intracellular and extracellular levels (Figure 4B and Figure 5B). Furthermore, the models were correctly validated, obtaining good correlation (R2Y) and predictive (Q2) scores and significant CV-ANOVA p-values (R2Y > 0.6, Q2 > 0.6, R2Y ≠ Q2 < 0.3, p-value < 0.05) (Table S6).
Major contributing metabolites for group separation were identified after the validation of the estimators (coefficient value and VIP score) using the Jackknife method as shown in an S-line plot. A total of seven and five metabolites at the intracellular and extracellular levels, respectively, were selected as potential biomarkers of antimony resistance. At the intracellular level, SbIII-resistant parasites showed an upregulation of alanine, proline, arginine, and lysine, and a downregulation of sn-glycero-3-phosphocholine, betaine, and acetate (Figure 4C). Additionally, at the extracellular level, SbIII-resistant parasites showed lower levels of proline, valine, lactate, and threonine (Figure 5C). Some of the signals with high correlation values (color bars) are not shown because they were not significant by the Jackknife method or because the associated metabolites could not be identified by signal overlapping (Figure 4C and Figure 5C). The VIP score and the coefficient size are summarized in the Figure S3.

3.4. Proline and Lactate Changed Linearly with the Antimony Resistance Level

Once the potentially important metabolites differentiating SbIII-sensitive and SbIII-resistant parasites were identified from the OPLS-DA models, we then explored whether the expression of those metabolites exhibited a significant linear correlation between the SbIII resistance level and the metabolite abundance. Interestingly, proline metabolite abundance showed a positive correlation with the SbIII resistance level in intracellular extracts, while both proline and lactate metabolite abundance showed a negative correlation with SbIII resistance level in extracellular extracts (Figure 6 and Figure 7).
To evaluate the importance of discriminating metabolites, we applied the receiver operating characteristic (ROC) curve analysis to the metabolites with a higher correlation value as a binary classifier (sensitive vs. resistant strains). ROC analyses showed optimal sensitivity and specificity scores with areas under the curve (AUCs) of 1 for all evaluated metabolites in a univariate approach (Figure 6 and Figure 7). However, these AUC scores should be carefully interpreted since here we analyzed only one Leishmania species using an in vitro approach, and consequently the metabolites were described as potential biomarkers.
Even though glycerol-phosphocholine and betaine did not show a significant Pearson correlation, they still showed a high capability to differentiate SbIII-sensitive and SbIII-resistant parasites on the basis of the AUC score, and thus they were also included (Figure 6).

4. Discussion

Antimony resistance remains an important determinant in leishmaniasis treatment, and biomarkers of resistant phenotypes can help to select an optimal medication. Here, we applied metabolomics to L. tropica strains with different resistance levels to SbIII as biological model. Clinical isolates of L. tropica are known to exhibit progressive reduction of antimony susceptibility, demanding strategies to identify and prevent the dissemination of drug-resistant strains [17]. Through the multivariate statistical analysis, we were able to identify 10 different metabolites (proline, alanine, arginine, lysine, betaine, glycerol-phosphocholine, lactate, acetate, threonine, and valine) differentiating the metabolic profiles of SbIII-resistant and -sensitive parasites at the intracellular and/or extracellular levels (Figure 4C and Figure 5C).
Proline was the only significant metabolite differentiating antimony-sensitive and -resistant parasites in both intracellular and extracellular extracts. Specifically, a decreased proline abundance at the extracellular level and an increased proline abundance at the intracellular level were representative features of the SbIII-resistant phenotypes. A similar observation was reported also for another Leishmania species: L. infantum promastigotes selected for SbIII resistance, and L. donovani promastigotes derived from clinical isolates; both of which showed higher proline concentration at the intracellular level [37,40]. Additionally, proline upregulation was also observed during both logarithmic and stationary phases of SbIII-resistant Leishmania growing under drug pressure [32]. However, these studies did not include the exometabolomic analysis.
There was also the observation that the parasite showed better tolerance to SbIII in the medium with high proline concentration but not against other drugs such as amphotericin or miltefosine [57]. Remarkably, in our study, we found that proline concentration was differentially increased in SbIII-resistant parasites at the intracellular level. Additionally, the exometabolomic analysis confirmed that resistant parasites depleted proline from the medium at higher level. The fact that the changes of proline metabolism can be detected at the extracellular level makes the exometabolomic analysis useful for SbIII resistance profiling. Interestingly, in this study, both intracellular and extracellular analysis showed that proline concentration changes linearly with the drug resistance level. This evidence indicates that proline levels are closely associated with the parasite’s ability to efficiently combat the drug and could be considered a robust metabolite biomarker of antimony resistance in Leishmania (Figure 6 and Figure 7).
In different types of eukaryotic cells including Leishmania, proline acts as an osmolyte and regulates reactive oxygen species (ROS). Is commonly accepted that antimonials induce a ROS imbalance [5,6]. Interestingly, here we also verified that the level of resistance to SbIII was associated with better response to combat the oxidative stress induced in vitro by hydrogen peroxide (H2O2). This free radical is commonly produced by the host macrophages during phagocytosis, inducing promastigote death in an apoptosis-like manner [58] (Figure S4).
Notably, proline analogues have been explored as promising therapeutic alternatives in trypanosomes [59]. In Leishmania parasites, ω-amino acid 4-amino-butiric (GABA), a proline analogue, inhibits the proline/alanine transporter and causes toxicity in both promastigotes and amastigotes [60]. This indicates that proline transport inhibitors could potentially be used to reverse SbIII-resistant phenotypes. More studies are needed to evaluate the effect of combined proline analogue and antimonial therapy.
Alanine, another amino acid contributing to the osmotic balance in Leishmania, was also found to be significantly increased in resistant parasites [61]. Both alanine and proline have common transport in Leishmania and it has been hypothesized that the regulation of the proline/alanine transporters such as LdAAP24 are key to osmotic stress responses in environments such as those induced by SbIII [62]. Another study has also detected alanine to be differentially expressed in sensitive parasites growing under SbIII pressure [36], suggesting that the alanine increase may not be exclusive of resistant phenotypes.
At the intracellular level, betaine and glycerol-phosphocholine were significantly decreased in SbIII-resistant parasites at the intracellular level. Betaine is a donor of methyl groups while glycerol-phosphocholine is a precursor for choline biosynthesis [63]. Choline can feed into phosphatidylcholine biosynthesis but also can be oxidized to produce betaine [64]. There is an evidence that phosphatidylcholine biosynthesis is decreased in Leishmania SbIII-resistant parasites growing under drug challenge, supporting the idea that glycerol-phosphocholine degradation taking place in SbIII-resistant parasites is necessary to increase methyl donor availability [40,65].
Together, these observations suggest that the decrease detected in betaine and sn-glycero-3-phosphocholine could be associated with their ability to donate or release methyl groups. Methyl groups are highly required to feed the methionine cycle and then the polyamine biosynthesis, two pathways interconnected by the decarboxylated S-adenosyl-L-methionine. Consistently, polyamine biosynthesis and the downstream trypanothione pathway have been widely recognized to be activated in SbIII-resistant parasites challenged with SbIII [5,66].
Other amino acids were also differentially detected. Arginine was increased in SbIII-resistant parasites at the intracellular level. This profile has been also detected in L. infantum promastigotes [37], while an arginine downregulation has been associated with sensitive parasites growing under antimony pressure [36]. It has been suggested that Leishmania parasites can sense arginine availability regulating the expression of Leishmania arginine transporter (AAP3) under stress conditions [67]. Interestingly, there is evidence that the AAP3 transporter is upregulated at the mRNA level in SbIII-resistant parasites derived from clinical isolates [68], which could explain the tendency toward a higher arginine intracellular concentration in SbIII-resistant phenotypes. Arginine is also interconnected with polyamine metabolism, producing trypanothione downstream via arginase, which catalyzes the enzymatic hydrolysis of L-arginine (L-Arg) to L-ornithine and urea [69].
We observed that lysine was increased in SbIII-resistant parasites at the intracellular level during the stationary phase. Conversely, other studies found that lysine was significantly increased in SbIII-resistant parasites during logarithmic phase [70]. It is possible that these opposing findings could be explained by the fact that the analyses were completed in different Leishmania species during different growth phases. However, more studies are required to clarify the role of lysine in antimony resistance.
Besides amino acids, we also selected other metabolites that could be important to differentiate the phenotypes of SbIII resistance. In Leishmania, acetate is produced from acetyl-CoA, a product of the glycolytic pathway [71]. In trypanosomes, acetate has been shown to be an important precursor for lipid biosynthesis [72]. It is possible that the observed decline in acetate is associated with the energy and lipid production. This is supported by evidence that suggests that energy metabolism and lipid remodeling play an important role in SbIII resistance [40,73].
Additionally, other metabolites involved in energy metabolism were detected. Lactate is a product of glucose catabolism. In Leishmania, glycolysis can take place in aerobic as well as in anaerobic conditions, and the lactate product is secreted [71]. We detected that lactate was increased in sensitive parasites at the extracellular level, suggesting that in response to SbIII, there is a change in energy metabolism. Consistent with this, it is commonly known that while the glycolysis is inhibited in SbIII-sensitive parasites, resistant parasites overexpress this pathway [26,74,75].
Furthermore, our findings suggest that resistant parasites more efficiently utilize threonine and valine since these metabolites were detected in lower concentrations in extracellular extracts. Threonine is another exogenous amino acid that can be used by Leishmania to supply protein synthesis and as a carbon source to feed the tricarboxylic acid cycle (TCA) [76]. Valine has been shown to be critical for Leishmania viability, and L-valine is degraded into pathways that lead to the production of succinyl coenzyme A, an important precursor that can be used in TCA cycle or fatty acid β-oxidation [77].
In summary, the metabolomic profiling by 1H-HMR allows for a clear differentiation between SbIII-resistant and -sensitive parasites (Figure 8). The metabolites identified in this study were associated with two general processes. The first group of metabolites included metabolites responding to osmotic and oxidative stress (proline, alanine, arginine, betaine, glycerol-phosphocholine). Some of these metabolites also contribute to the activation of thiol metabolism, which induces drug inactivation by production of the thiol–SbIII complex and promotes the oxidative stress balance via trypanothione. Trypanothione system, the parasite’s equivalent to the mammalian glutathione, is the main detoxifying system against oxidative damage. Since Leishmania lacks catalases, trypanothione regeneration relies on nicotinamide adenine dinucleotide phosphate (NADPH) [78]. As consequence, thiol metabolism requires a large amount of energy, demanding continuous production of reductant in the form of nicotinamide adenine dinucleotide phosphate (NADPH) [79].
The second group of metabolites (acetate, threonine, valine, and lactate) can contribute to energy production via pathways such as the TCA cycle [80]. TCA cycle, malic enzyme, or pentose phosphate pathway can fuel the NADPH required for thiol metabolism activation (first group of metabolites). The coupling between the activation of thiol metabolism and energy production is essential for resistant parasites to combat the drug efficiently (Figure 8).
Currently, there are no reliable molecular methods to determine drug susceptibility in clinical settings, and treatment failure due to drug resistance is a big problem in Leishmania species [23]. In the absence of validated standardized methods and markers to track drug resistance in clinical isolates, the further investigation and validation of identified drug resistance biomarkers in a clinical setting is crucial. More studies are also required to validate the proposed biomarkers in other Leishmania species due to the high complexity of the model. Likewise, studies using gene-editing techniques are needed to clarify the metabolic adaptations under drug-induced stress conditions.

5. Conclusions

SbIII-resistant and -sensitive parasites showed significant changes in the metabolite composition and these phenotypes can be potentially predicted using 1H-NMR profiling. The differences in metabolite composition suggest that the development of resistance to SbIII involves an optimized response to the osmotic/oxidative stress induced by the drug as well as a rearrangement of the carbon metabolism possibly to compensate the production of energy in form of reductant, typically required to keep the antioxidant response active.
At the intracellular and extracellular levels, metabolites such as proline and lactate are informative not only of the resistant phenotype but also of the level of resistance, suggesting that the metabolic pathways in which these compounds are involved are closely related with the mechanism activated to combat the drug, and consequently these metabolites can be prioritized as a more robust biomarker of SbIII resistance.
In many studies, metabolomic analysis has been performed at the intracellular level; however, our study showed that simpler and time-saving exometabolomic analysis can be also efficiently used for differentiation of sensitive and resistant parasites.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/cells10051063/s1, Figure S1: Enrichment of parasites without membrane damage for intracellular metabolomic analysis. Figure S2: Comparison of the curves of growth between the 3 L. tropica strains. Table S1: Chemical shifts per compound detected in 1H-NMR spectra of Leishmania’s intracellular extracts. Table S2: Chemical shifts per compound detected in a 1H-NMR spectrum of Leishmania’s extracellular extracts. Table S3: Raw data from Venn diagram comparing the different and common compounds detected at the intracellular and extracellular level in Leishmania parasites by 1H-NMR. Table S4: Total compounds detected by 1H-NMR in L. tropica distributed by the metabolites main class. Table S5: Total compounds distributed by KEGG metabolic pathway. Table S6: Comparison of statistical measures calculated for the supervised OPLS-DA models for the differentiation of SbIII-sensitive and -resistant parasites at the intracellular and extracellular levels. Figure S3: Regression coefficients and VIP scores based on the metabolomic profiling dataset for the OPLS-DA models comparing SbIII-resistant parasites versus SbIII-sensitive parasites. Figure S4: SbIII-resistant parasites showed better tolerance to oxidative stress induced by hydrogen peroxide.

Author Contributions

S.A.G.G., Z.N.K., E.G., and C.E.M. contributed to conceptualization; S.A.G.G. wrote the manuscript, analyzed the data, and designed and prepared the figures; S.A.G.G., Z.N.K., E.G., and C.E.M. contributed to the experimental design. E.G. and S.A.G.G. completed the 1H-NMR spectrum acquisition. All authors discussed and edited the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by Minciencias and Gobernación del Tolima-Colombia as part of the program 755-2016.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding authors.

Acknowledgments

We acknowledge the laboratory “Programa de Estudio y Control de Enfermedades Tropicales” (PECET) of the University of Antioquia at Medellín, Colombia, for providing the Leishmania strains used in the study. We would like to thank A.L. Karamyshev for the critical reading of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. WHO. Leishmaniasis Fact Sheet. 2020. Available online: http://www.who.int/mediacentre/factsheets/fs375/en/ (accessed on 13 March 2021).
  2. PAHO/WHO. Leishmaniasis Fact Sheet. 2017. Available online: https://www.paho.org/hq/dmdocuments/2017/2017-cha-leishmaniasis-factsheet-work.pdf (accessed on 13 March 2021).
  3. Chappuis, F.; Sundar, S.; Hailu, A.; Ghalib, H.; Rijal, S.; Peeling, R.W.; Alvar, J.; Boelaert, M. Visceral leishmaniasis: What are the needs for diagnosis, treatment and control? Nat. Rev. Microbiol. 2007, 5, 873–882. [Google Scholar] [CrossRef] [PubMed]
  4. Torres-Guerrero, E.; Quintanilla-Cedillo, M.R.; Ruiz-Esmenjaud, J.; Arenas, R. Leishmaniasis: A review. F1000Research 2017, 6, 750. [Google Scholar] [CrossRef]
  5. Mandal, G.; Wyllie, S.; Singh, N.; Sundar, S.; Fairlamb, A.H.; Chatterjee, M. Increased levels of thiols protect antimony unresponsive Leishmania donovani field isolates against reactive oxygen species generated by trivalent antimony. Parasitology 2007, 134, 1679–1687. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Tirmenstein, M.A.; Plews, P.I.; Walker, C.V.; Woolery, M.D.; Wey, H.E.; Toraason, M.A. Antimony-Induced Oxidative Stress and Toxicity in Cultured Cardiac Myocytes. Toxicol. Appl. Pharmacol. 1995, 130, 41–47. [Google Scholar] [CrossRef]
  7. Berman, J.D.; Gallalee, J.V.; Best, J.M. Sodium stibogluconate (Pentostam) inhibition of glucose catabolism via the glycolytic pathway, and fatty acid beta-oxidation in Leishmania mexicana amastigotes. Biochem. Pharmacol. 1987, 36, 197–201. [Google Scholar] [CrossRef]
  8. Ferreira, C.D.S.; Pimenta, A.M.D.C.; Demicheli, C.; Frézard, F. Characterization of reactions of antimoniate and meglumine antimoniate with a guanine ribonucleoside at different pH. Biometals 2006, 19, 573–581. [Google Scholar] [CrossRef] [PubMed]
  9. Walker, J.; Saravia, N.G. Inhibition of Leishmania donovani promastigote DNA topoisomerase I and human monocyte DNA topoisomerases I and II by antimonial drugs and classical antitopoisomerase agents. J. Parasitol. 2004, 90, 1155–1162. [Google Scholar] [CrossRef]
  10. Frézard, F.; Silva, H.; Pimenta, A.M.; Farrell, N.; Demicheli, C. Greater binding affinity of trivalent antimony to a CCCH zinc finger domain compared to a CCHC domain of kinetoplastid proteins. Metallomics 2012, 4, 433–440. [Google Scholar] [CrossRef] [PubMed]
  11. Sundar, S.; More, D.K.; Singh, M.K.; Singh, V.P.; Sharma, S.; Makharia, A.; Kumar, P.C.K.; Murray, H.W. Failure of Pentavalent Antimony in Visceral Leishmaniasis in India: Report from the Center of the Indian Epidemic. Clin. Infect. Dis. 2000, 31, 1104–1107. [Google Scholar] [CrossRef] [Green Version]
  12. Rojas, R.; Valderrama, L.; Valderrama, M.; Varona, M.X.; Ouellette, M.; Saravia, N.G. Resistance to Antimony and Treatment Failure in Human Leishmania (Viannia) Infection. J. Infect. Dis. 2006, 193, 1375–1383. [Google Scholar] [CrossRef] [Green Version]
  13. Llanos-Cuentas, A.; Tulliano, G.; Araujo-Castillo, R.; Miranda-Verastegui, C.; Santamaria-Castrellon, G.; Ramirez, L.; Lazo, M.; De Doncker, S.; Boelaert, M.; Robays, J.; et al. Clinical and parasite species risk factors for pentavalent antimonial treatment failure in cutaneous leishmaniasis in Peru. Clin. Infect. Dis. 2008, 46, 223–231. [Google Scholar] [CrossRef] [PubMed]
  14. Arevalo, J.; Ramirez, L.; Adaui, V.; Zimic, M.; Tulliano, G.; Miranda-Verástegui, C.; Lazo, M.; Loayza-Muro, R.; De Doncker, S.; Maurer, A.; et al. Influence of Leishmania (Viannia) species on the response to antimonial treatment in patients with American tegumentary leishmaniasis. J. Infect. Dis. 2007, 195, 1846–1851. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Unger, A.; O’Neal, S.; Machado, P.R.; Guimarães, L.H.; Morgan, D.J.; Schriefer, A.; Bacellar, O.; Glesby, M.J.; Carvalho, E.M. Association of treatment of American cutaneous leishmaniasis prior to ulcer development with high rate of failure in northeastern Brazil. Am. J. Trop. Med. Hyg. 2009, 80, 574–579. [Google Scholar] [CrossRef] [Green Version]
  16. Soto, J.; Toledo, J.; Vega, J.; Berman, J. Short report: Efficacy of pentavalent antimony for treatment of colombian cutaneous leishmaniasis. Am. J. Trop. Med. Hyg. 2005, 72, 421–422. [Google Scholar] [CrossRef] [PubMed]
  17. Eddaikra, N.; Ait-Oudhia, K.; Kherrachi, I.; Oury, B.; Moulti-Mati, F.; Benikhlef, R.; Harrat, Z.; Sereno, D. Antimony susceptibility of Leishmania isolates collected over a 30-year period in Algeria. PLoS Negl. Trop. Dis. 2018, 12, e0006310. [Google Scholar] [CrossRef] [Green Version]
  18. Kumar, D.; Kulshrestha, A.; Singh, R.; Salotra, P. In vitro susceptibility of field isolates of Leishmania donovani to Miltefosine and amphotericin B: Correlation with sodium antimony gluconate susceptibility and implications for treatment in areas of endemicity. Antimicrob. Agents Chemother. 2009, 53, 835–838. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Purkait, B.; Kumar, A.; Nandi, N.; Sardar, A.H.; Das, S.; Kumar, S.; Pandey, K.; Ravidas, V.; Kumar, M.; De, T.; et al. Mechanism of Amphotericin B Resistance in Clinical Isolates of Leishmania donovani. Antimicrob. Agents Chemother. 2012, 56, 1031–1041. [Google Scholar] [CrossRef] [Green Version]
  20. Srivastava, S.; Mishra, J.; Gupta, A.K.; Singh, A.; Shankar, P.; Singh, S. Laboratory confirmed miltefosine resistant cases of visceral leishmaniasis from India. Parasites Vectors 2017, 10, 1–11. [Google Scholar] [CrossRef] [Green Version]
  21. 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. 2003, 130, 83–90. [Google Scholar] [CrossRef] [Green Version]
  22. PAHO. Manual de Diagnóstico y Tratamiento de la Leishmaniasis. 2018. Available online: https://www.paho.org/par/index.php?option=com_docman&view=download&alias=575-manual-de-diagnostico-y-tratamiento-de-las-leishmaniasis&category_slug=publicaciones-con-contrapartes&Itemid=253 (accessed on 10 March 2021).
  23. Sereno, D.; Harrat, Z.; Eddaikra, N. Meta-analysis and discussion on challenges to translate Leishmania drug resistance phenotyping into the clinic. Acta Trop. 2019, 191, 204–211. [Google Scholar] [CrossRef]
  24. Karczewski, K.J.; Snyder, M.P. Integrative omics for health and disease. Nat. Rev. Genet. 2018, 19, 299–310. [Google Scholar] [CrossRef] [PubMed]
  25. Laffitte, M.N.; Leprohon, P.; Papadopoulou, B.; Ouellette, M. Plasticity of the Leishmania genome leading to gene copy number variations and drug resistance. F1000Research 2016, 5, 2350. [Google Scholar] [CrossRef] [Green Version]
  26. Patino, L.H.; Muskus, C.; Ramírez, J.D. Transcriptional responses of Leishmania (Leishmania) amazonensis in the presence of trivalent sodium stibogluconate. Parasites Vectors 2019, 12, 348. [Google Scholar] [CrossRef] [Green Version]
  27. 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. 2008, 79, 69–75. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Singh, N.; Almeida, R.; Kothari, H.; Kumar, P.; Mandal, G.; Chatterjee, M.; Venkatachalam, S.; Govind, M.K.; Mandal, S.K.; Sundar, S. Differential gene expression analysis in antimony-unresponsive Indian kala azar (visceral leishmaniasis) clinical isolates by DNA microarray. Parasitology 2007, 134, 777–787. [Google Scholar] [CrossRef]
  29. Mookerjee Basu, J.; Mookerjee, A.; Banerjee, R.; Saha, M.; Singh, S.; Naskar, K.; Tripathy, G.; Sinha, P.K.; Pandey, K.; Sundar, S.; et al. Inhibition of ABC transporters abolishes antimony resistance in Leishmania Infection. Antimicrob. Agents Chemother. 2008, 52, 1080–1093. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Jeddi, F.; Piarroux, R.; Mary, C. Antimony Resistance in Leishmania, Focusing on Experimental Research. J. Trop. Med. 2011, 2011. [Google Scholar] [CrossRef] [Green Version]
  31. Karamysheva, Z.N.; Gutierrez Guarnizo, S.A.; Karamyshev, A.L. Regulation of Translation in the Protozoan Parasite Leishmania. Int. J. Mol. Sci. 2020, 21, 2981. [Google Scholar] [CrossRef]
  32. Berg, M.; Vanaerschot, M.; Jankevics, A.; Cuypers, B.; Maes, I.; Mukherjee, S.; Khanal, B.; Rijal, S.; Roy, S.; Opperdoes, F.; et al. Metabolic adaptations of Leishmania donovani in relation to differentiation, drug resistance, and drug pressure. Mol. Microbiol. 2013, 90, 428–442. [Google Scholar] [CrossRef] [PubMed]
  33. Brotherton, M.C.; Bourassa, S.; Leprohon, P.; Légaré, D.; Poirier, G.G.; Droit, A.; Ouellette, M. Proteomic and genomic analyses of antimony resistant Leishmania infantum mutant. PLoS ONE 2013, 8, e81899. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Gutierrez Guarnizo, S.A.; Tikhonova, E.B.; Zabet-Moghaddam, M.; Zhang, K.; Muskus, C.; Karamyshev, A.L.; Karamysheva, Z.N. Drug-Induced Lipid Remodeling in Leishmania Parasites. Microorganisms 2021, 9, 790. [Google Scholar] [CrossRef]
  35. Pinu, F.R.; Goldansaz, S.A.; Jaine, J. Translational Metabolomics: Current Challenges and Future Opportunities. Metabolites 2019, 9, 108. [Google Scholar] [CrossRef] [Green Version]
  36. Canuto, G.A.; Castilho-Martins, E.A.; Tavares, M.; López-Gonzálvez, A.; Rivas, L.; Barbas, C. CE-ESI-MS metabolic fingerprinting of Leishmania resistance to antimony treatment. Electrophoresis 2012, 33, 1901–1910. [Google Scholar] [CrossRef] [PubMed]
  37. Rojo, D.; Canuto, G.A.B.; Castilho-Martins, E.A.; Tavares, M.F.M.; Barbas, C.; López-Gonzálvez, Á.; Rivas, L. A Multiplatform Metabolomic Approach to the Basis of Antimonial Action and Resistance in Leishmania infantum. PLoS ONE 2015, 10, e0130675. [Google Scholar] [CrossRef]
  38. t’Kindt, R.; Jankevics, A.; Scheltema, R.A.; Zheng, L.; Watson, D.G.; Dujardin, J.C.; Breitling, R.; Coombs, G.H.; Decuypere, S. Towards an unbiased metabolic profiling of protozoan parasites: Optimisation of a Leishmania sampling protocol for HILIC-orbitrap analysis. Anal. Bioanal. Chem. 2010, 398, 2059–2069. [Google Scholar] [CrossRef] [Green Version]
  39. Pinu, F.R.; Villas-Boas, S.G. Extracellular Microbial Metabolomics: The State of the Art. Metabolites 2017, 7, 43. [Google Scholar] [CrossRef] [PubMed]
  40. t’Kindt, R.; Scheltema, R.A.; Jankevics, A.; Brunker, K.; Rijal, S.; Dujardin, J.-C.; Breitling, R.; Watson, D.G.; Coombs, G.H.; Decuypere, S. Metabolomics to Unveil and Understand Phenotypic Diversity between Pathogen Populations. PLoS Negl. Trop. Dis. 2010, 4, e904. [Google Scholar] [CrossRef]
  41. Djukovic, D.; Nagana Gowda, G.A.; Raftery, D. Chapter 18—Mass Spectrometry and NMR Spectroscopy–Based Quantitative Metabolomics. In Proteomic and Metabolomic Approaches to Biomarker Discovery; Issaq, H.J., Veenstra, T.D., Eds.; Academic Press: Boston, MA, USA, 2013; pp. 279–297. [Google Scholar]
  42. Emwas, A.H.; Roy, R.; McKay, R.T.; Tenori, L.; Saccenti, E.; Gowda, G.A.N.; Raftery, D.; Alahmari, F.; Jaremko, L.; Jaremko, M.; et al. NMR Spectroscopy for Metabolomics Research. Metabolites 2019, 9, 123. [Google Scholar] [CrossRef] [Green Version]
  43. Ser, Z.; Liu, X.; Tang, N.N.; Locasale, J.W. Extraction parameters for metabolomics from cultured cells. Anal. Biochem. 2015, 475, 22–28. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Kostidis, S.; Addie, R.D.; Morreau, H.; Mayboroda, O.A.; Giera, M. Quantitative NMR analysis of intra- and extracellular metabolism of mammalian cells: A tutorial. Anal. Chim. Acta 2017, 980, 1–24. [Google Scholar] [CrossRef]
  45. McKay, R.T. How the 1D-NOESY suppresses solvent signal in metabonomics NMR spectroscopy: An examination of the pulse sequence components and evolution. Concepts Magn. Reson. Part A 2011, 38, 197–220. [Google Scholar] [CrossRef]
  46. Mickiewicz, B.; Hyndman, M.; Vogel, H. Metabolite Profiling of Clinical Cancer Biofluid Samples by NMR Spectroscopy: Methods and Protocols; Humana Press: New York, NY, USA, 2019; Volume 1928, pp. 251–274. [Google Scholar]
  47. Dieterle, F.; Ross, A.; Schlotterbeck, G.; Senn, H. Probabilistic Quotient Normalization as Robust Method to Account for Dilution of Complex Biological Mixtures. Application in 1H NMR Metabonomics. Anal. Chem. 2006, 78, 4281–4290. [Google Scholar] [CrossRef] [PubMed]
  48. Mickiewicz, B.; Thompson, G.C.; Blackwood, J.; Jenne, C.N.; Winston, B.W.; Vogel, H.J.; Joffe, A.R. Biomarker Phenotype for Early Diagnosis and Triage of Sepsis to the Pediatric Intensive Care Unit. Sci. Rep. 2018, 8, 16606. [Google Scholar] [CrossRef]
  49. Yang, J.; Zhao, X.; Lu, X.; Lin, X.; Xu, G. A data preprocessing strategy for metabolomics to reduce the mask effect in data analysis. Front. Mol. Biosci. 2015, 2, 4. [Google Scholar] [CrossRef] [Green Version]
  50. Eriksson, L.; Johansson, E.; Kettaneh-Wold, N.; Trygg, J.; Wikstrom, C.; Wold, S. Multi- and Megavariate Data Analysis Part I: Basic Principles and Applications; Umetrics AB: Umeå, Sweden, 2006. [Google Scholar]
  51. Worley, B.; Powers, R. PCA as a practical indicator of OPLS-DA model reliability. Curr. Metab. 2016, 4, 97–103. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Cloarec, O.; Dumas, M.-E.; Craig, A.; Barton, R.H.; Trygg, J.; Hudson, J.; Blancher, C.; Gauguier, D.; Lindon, J.C.; Holmes, E.; et al. Statistical Total Correlation Spectroscopy:  An Exploratory Approach for Latent Biomarker Identification from Metabolic 1H NMR Data Sets. Anal. Chem. 2005, 77, 1282–1289. [Google Scholar] [CrossRef] [PubMed]
  53. Metz, C.E. Basic principles of ROC analysis. Semin. Nucl. Med. 1978, 8, 283–298. [Google Scholar] [CrossRef]
  54. Karampetsou, K.; Koutsoni, O.S.; Dotsika, E. Quantification of Nitric Oxide and Reactive Oxygen Species in Leishmania-infected J774A.1 Macrophages as a Response to the in vitro treatment with a Natural Product Compound. Bio Protoc 2019, 9, e3442. [Google Scholar] [CrossRef] [PubMed]
  55. Arjmand, M.; Madrakian, A.; Khalili, G.; Najafi Dastnaee, A.; Zamani, Z.; Akbari, Z. Metabolomics-Based Study of Logarithmic and Stationary Phases of Promastigotes in Leishmania major by 1H NMR Spectroscopy. Iran. Biomed. J. 2016, 20, 77–83. [Google Scholar] [CrossRef]
  56. Want, E.J.; O’Maille, G.; Smith, C.A.; Brandon, T.R.; Uritboonthai, W.; Qin, C.; Trauger, S.A.; Siuzdak, G. Solvent-dependent metabolite distribution, clustering, and protein extraction for serum profiling with mass spectrometry. Anal. Chem. 2006, 78, 743–752. [Google Scholar] [CrossRef] [PubMed]
  57. Berg, M.; García-Hernández, R.; Cuypers, B.; Vanaerschot, M.; Manzano, J.I.; Poveda, J.A.; Ferragut, J.A.; Castanys, S.; Dujardin, J.-C.; Gamarro, F. Experimental Resistance to Drug Combinations in Leishmania donovani: Metabolic and Phenotypic Adaptations. Antimicrob. Agents Chemother. 2015, 59, 2242–2255. [Google Scholar] [CrossRef] [Green Version]
  58. Das, M.; Mukherjee, S.B.; Shaha, C. Hydrogen peroxide induces apoptosis-like death in Leishmania donovani promastigotes. J. Cell Sci. 2001, 114, 2461–2469. [Google Scholar] [PubMed]
  59. Sayé, M.; Gauna, L.; Valera-Vera, E.; Reigada, C.; Miranda, M.R.; Pereira, C.A. Crystal violet structural analogues identified by in silico drug repositioning present anti-Trypanosoma cruzi activity through inhibition of proline transporter TcAAAP069. PLoS Negl. Trop. Dis. 2020, 14, e0007481. [Google Scholar] [CrossRef] [Green Version]
  60. Prati, F.; Goldman-Pinkovich, A.; Lizzi, F.; Belluti, F.; Koren, R.; Zilberstein, D.; Bolognesi, M.L. Quinone-Amino Acid Conjugates Targeting Leishmania Amino Acid Transporters. PLoS ONE 2014, 9, e107994. [Google Scholar] [CrossRef]
  61. Burrows, C.; Blum, J.J. Effect of hyper-osmotic stress on alanine content of Leishmania major promastigotes. J. Protozool. 1991, 38, 47–52. [Google Scholar] [CrossRef]
  62. Inbar, E.; Schlisselberg, D.; Suter Grotemeyer, M.; Rentsch, D.; Zilberstein, D. A versatile proline/alanine transporter in the unicellular pathogen Leishmania donovani regulates amino acid homoeostasis and osmotic stress responses. Biochem. J. 2012, 449, 555–566. [Google Scholar] [CrossRef]
  63. Obeid, R. The metabolic burden of methyl donor deficiency with focus on the betaine homocysteine methyltransferase pathway. Nutrients 2013, 5, 3481–3495. [Google Scholar] [CrossRef] [Green Version]
  64. Fagone, P.; Jackowski, S. Phosphatidylcholine and the CDP-choline cycle. Biochim. Biophys. Acta 2013, 1831, 523–532. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Shaw, C.D.; Lonchamp, J.; Downing, T.; Imamura, H.; Freeman, T.M.; Cotton, J.A.; Sanders, M.; Blackburn, G.; Dujardin, J.C.; Rijal, S.; et al. In vitro selection of miltefosine resistance in promastigotes of Leishmania donovani from Nepal: Genomic and metabolomic characterization. Mol. Microbiol. 2016, 99, 1134–1148. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. 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 2014, 7, 596. [Google Scholar] [CrossRef] [Green Version]
  67. Goldman-Pinkovich, A.; Kannan, S.; Nitzan-Koren, R.; Puri, M.; Pawar, H.; Bar-Avraham, Y.; McDonald, J.; Sur, A.; Zhang, W.-W.; Matlashewski, G.; et al. Sensing Host Arginine Is Essential for Leishmania Parasites’ Intracellular Development. MBio 2020, 11. [Google Scholar] [CrossRef] [PubMed]
  68. Kazemi-Rad, E.; Mohebali, M.; Khadem-Erfan, M.B.; Hajjaran, H.; Hadighi, R.; Khamesipour, A.; Rezaie, S.; Saffari, M.; Raoofian, R.; Heidari, M. Overexpression of ubiquitin and amino acid permease genes in association with antimony resistance in Leishmania tropica field isolates. Korean J. Parasitol. 2013, 51, 413–419. [Google Scholar] [CrossRef] [PubMed]
  69. Colotti, G.; Ilari, A. Polyamine metabolism in Leishmania: From arginine to trypanothione. Amino Acids 2011, 40, 269–285. [Google Scholar] [CrossRef] [PubMed]
  70. Dumetz, F.; Cuypers, B.; Imamura, H.; Zander, D.; D’Haenens, E.; Maes, I.; Domagalska, M.A.; Clos, J.; Dujardin, J.C.; De Muylder, G. Molecular Preadaptation to Antimony Resistance in Leishmania donovani on the Indian Subcontinent. MSphere 2018, 3, e00548-17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  71. Saunders, E.C.; Souza, D.P.; Naderer, T.; Sernee, M.F.; Ralton, J.E.; Doyle, M.A.; Macrae, J.I.; Chambers, J.L.; Heng, J.; Nahid, A.; et al. Central carbon metabolism of Leishmania parasites. Parasitology 2010, 137, 1303–1313. [Google Scholar] [CrossRef] [PubMed]
  72. Rivière, L.; Moreau, P.; Allmann, S.; Hahn, M.; Biran, M.; Plazolles, N.; Franconi, J.-M.; Boshart, M.; Bringaud, F. Acetate produced in the mitochondrion is the essential precursor for lipid biosynthesis in procyclic trypanosomes. Proc. Natl. Acad. Sci. USA 2009, 106, 12694–12699. [Google Scholar] [CrossRef] [Green Version]
  73. Cuypers, B.; Berg, M.; Imamura, H.; Dumetz, F.; De Muylder, G.; Domagalska, M.A.; Rijal, S.; Bhattarai, N.R.; Maes, I.; Sanders, M.; et al. Integrated genomic and metabolomic profiling of ISC1, an emerging Leishmania donovani population in the Indian subcontinent. Infect. Genet. Evol. 2018, 62, 170–178. [Google Scholar] [CrossRef]
  74. Matrangolo, F.S.V.; Liarte, D.B.; Andrade, L.C.; de Melo, M.F.; Andrade, J.M.; Ferreira, R.F.; Santiago, A.S.; Pirovani, C.P.; Silva-Pereira, R.A.; Murta, S.M.F. Comparative proteomic analysis of antimony-resistant and -susceptible Leishmania braziliensis and Leishmania infantum chagasi lines. Mol. Biochem. Parasitol. 2013, 190, 63–75. [Google Scholar] [CrossRef]
  75. Zheng, Z.; Chen, J.; Ma, G.; Satoskar, A.R.; Li, J. Integrative genomic, proteomic and phenotypic studies of Leishmania donovani strains revealed genetic features associated with virulence and antimony-resistance. Parasites Vectors 2020, 13, 1–14. [Google Scholar] [CrossRef]
  76. Nayak, A. Comparative Study of Amino Acids Utilisation in Leishmania Mexicana. Ph.D. Thesis, University of Glasgow, Glasgow, UK, 2018. Available online: http://theses.gla.ac.uk/9136/ (accessed on 10 March 2021).
  77. Kohlmeier, M. Isoleucine. In Nutrient Metabolism; Kohlmeier, M., Ed.; Academic Press: London, UK, 2003; pp. 377–383. [Google Scholar]
  78. Murray, H.W. Susceptibility of Leishmania to oxygen intermediates and killing by normal macrophages. J. Exp. Med. 1981, 153, 1302–1315. [Google Scholar] [CrossRef] [Green Version]
  79. McConville, M.J.; Naderer, T. Metabolic pathways required for the intracellular survival of Leishmania. Annu. Rev. Microbiol. 2011, 65, 543–561. [Google Scholar] [CrossRef] [PubMed]
  80. Doyle, M.A.; MacRae, J.I.; De Souza, D.P.; Saunders, E.C.; McConville, M.J.; Likić, V.A. LeishCyc: A biochemical pathways database for Leishmania major. BMC Syst. Biol. 2009, 3, 57. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Cells 10 01063 g001 550
Figure 1. Experimental design to differentiate the metabolomic profiles in antimony-resistant and -sensitive L. tropica parasites. The workflow for intracellular and extracellular analysis is shown in parallel. (1) After parasite growth, the cell pellet and supernatant were collected and used for intracellular and extracellular analysis, respectively. (2) The methanol extracts were analyzed by 1H-HMR using a Noesygppr1d pulse. (3) Spectra processing included the annotation, apodization, alignment, phasing, baseline correction, bucketing, data normalization, and data scaling, integrating the following software: Chenomx v8.5, MestReNova v14.2, and SIMCA Umetrics v14.1. (4) Unsupervised multivariate analysis by PCA for preliminary grouping visualization and outlier detection based on the Hotelling t-squared statistic with 95% confidence limit. (5) Supervised multivariate analysis by OPLS-DA for the identification of statistically significant “variable importance in projection” (VIP) or potential biomarker, filtering by the Jackknife method to correct estimation bias. (6) Finally, the VIPs were analyzed to discuss the possible biological meaning.
Figure 1. Experimental design to differentiate the metabolomic profiles in antimony-resistant and -sensitive L. tropica parasites. The workflow for intracellular and extracellular analysis is shown in parallel. (1) After parasite growth, the cell pellet and supernatant were collected and used for intracellular and extracellular analysis, respectively. (2) The methanol extracts were analyzed by 1H-HMR using a Noesygppr1d pulse. (3) Spectra processing included the annotation, apodization, alignment, phasing, baseline correction, bucketing, data normalization, and data scaling, integrating the following software: Chenomx v8.5, MestReNova v14.2, and SIMCA Umetrics v14.1. (4) Unsupervised multivariate analysis by PCA for preliminary grouping visualization and outlier detection based on the Hotelling t-squared statistic with 95% confidence limit. (5) Supervised multivariate analysis by OPLS-DA for the identification of statistically significant “variable importance in projection” (VIP) or potential biomarker, filtering by the Jackknife method to correct estimation bias. (6) Finally, the VIPs were analyzed to discuss the possible biological meaning.
Cells 10 01063 g001
Cells 10 01063 g002 550
Figure 2. 1H-NMR spectra of L. tropica intracellular extracts. (A) Full 1H-NMR spectrum from δ 0.0 to δ 9.3. (B) Enlarged spectrum from δ 0.0 to δ 2.3. (C) Enlarged spectrum from δ 2.3 to δ 4.5. (D) Enlarged spectrum from δ 5.0 to δ 7.2. (E) Enlarged spectrum from δ 7.1 to δ 9.3. Peaks of 26 compounds: 1, valine; 2, isoleucine; 3, propylene glycol; 4, 3-hydroxyisovalerate; 5, lactate; 6, alanine; 7, arginine; 8, acetate; 9, proline; 10, methionine; 11, succinate; 12, β-alanine; 13, malate; 14, lysine; 15, glycine; 16, glycerophosphocholine; 17, betaine; 18, IMP; 19, fumarate; 20, ∂-methylhistidine; 21, desaminotyrosine; 22, xanthine; 23, hypoxanthine; 24, Formate; 25, AMP; 26, NAD. TSP-d4 as an internal standard with chemical shift at δ 0.0. Only one spectrum is represented for simplification. The spectrum corresponds to the HR strain.
Figure 2. 1H-NMR spectra of L. tropica intracellular extracts. (A) Full 1H-NMR spectrum from δ 0.0 to δ 9.3. (B) Enlarged spectrum from δ 0.0 to δ 2.3. (C) Enlarged spectrum from δ 2.3 to δ 4.5. (D) Enlarged spectrum from δ 5.0 to δ 7.2. (E) Enlarged spectrum from δ 7.1 to δ 9.3. Peaks of 26 compounds: 1, valine; 2, isoleucine; 3, propylene glycol; 4, 3-hydroxyisovalerate; 5, lactate; 6, alanine; 7, arginine; 8, acetate; 9, proline; 10, methionine; 11, succinate; 12, β-alanine; 13, malate; 14, lysine; 15, glycine; 16, glycerophosphocholine; 17, betaine; 18, IMP; 19, fumarate; 20, ∂-methylhistidine; 21, desaminotyrosine; 22, xanthine; 23, hypoxanthine; 24, Formate; 25, AMP; 26, NAD. TSP-d4 as an internal standard with chemical shift at δ 0.0. Only one spectrum is represented for simplification. The spectrum corresponds to the HR strain.
Cells 10 01063 g002
Cells 10 01063 g003 550
Figure 3. 1H-NMR spectra of L. tropica extracellular extract. (A) Full 1H-NMR spectrum from δ 0.0 to δ 9.5 (TSP-d4 as an internal standard with chemical shift at δ 0.0). (B) Enlarged spectrum from δ 0.0 to δ 3.1. (C) Enlarged spectrum from δ 3.1 to δ 6.3. (D) Enlarged spectrum from δ 6.3 to δ 9.5. Peaks of 29 compounds: 1, leucine; 2, isoleucine; 2, valine; 4, lactate; 5, threonine; 6, alanine; 7, lysine; 8, proline; 9, glutamine; 10, succinate, 11, ß-alanine; 12, malate; 13, aspartate; 14, arginine; 15, glycine; 16, l-serine; 17, myoinositol; 18, trehalose; 19, uracil; 20, fumarate; 21, N-acetyl tyrosine; 22, tyrosine; 23, phenylalanine; 24, tryptophan, 25, methylhistidine; 26, hypoxanthine; 27, formate; 28, imidazole, 29, nicotinate. Only one spectrum is represented for simplification. The spectrum corresponds to the HR strain.
Figure 3. 1H-NMR spectra of L. tropica extracellular extract. (A) Full 1H-NMR spectrum from δ 0.0 to δ 9.5 (TSP-d4 as an internal standard with chemical shift at δ 0.0). (B) Enlarged spectrum from δ 0.0 to δ 3.1. (C) Enlarged spectrum from δ 3.1 to δ 6.3. (D) Enlarged spectrum from δ 6.3 to δ 9.5. Peaks of 29 compounds: 1, leucine; 2, isoleucine; 2, valine; 4, lactate; 5, threonine; 6, alanine; 7, lysine; 8, proline; 9, glutamine; 10, succinate, 11, ß-alanine; 12, malate; 13, aspartate; 14, arginine; 15, glycine; 16, l-serine; 17, myoinositol; 18, trehalose; 19, uracil; 20, fumarate; 21, N-acetyl tyrosine; 22, tyrosine; 23, phenylalanine; 24, tryptophan, 25, methylhistidine; 26, hypoxanthine; 27, formate; 28, imidazole, 29, nicotinate. Only one spectrum is represented for simplification. The spectrum corresponds to the HR strain.
Cells 10 01063 g003
Cells 10 01063 g004 550
Figure 4. Multivariate analysis showing significant differences between the metabolomic profiles of sensitive and resistant Leishmania parasites at the intracellular level by 1H-NMR. (A) Score plot of PCA. (B) OPLS-DA score plot comparing antimony-resistant vs. antimony-sensitive parasites. (C) S-line plot highlighting the metabolites with major contributions to group separation. Each individual dot in panel (A) or (B) represents an observation or sample. Some dots can be overlapped. Blue dots: WT strain; green dots: MR strain; red dots: HR strain. The dotted ellipse of score plots describe the 95% confidence interval of the Hotelling’s T-squared distribution. The color bar in the S-line plot corresponds to the absolute value of the correlation loading in the discrimination model. A total of 9 observations are shown distributed in 3 independent biological replicates per experimental condition (WT, MR, and HR).
Figure 4. Multivariate analysis showing significant differences between the metabolomic profiles of sensitive and resistant Leishmania parasites at the intracellular level by 1H-NMR. (A) Score plot of PCA. (B) OPLS-DA score plot comparing antimony-resistant vs. antimony-sensitive parasites. (C) S-line plot highlighting the metabolites with major contributions to group separation. Each individual dot in panel (A) or (B) represents an observation or sample. Some dots can be overlapped. Blue dots: WT strain; green dots: MR strain; red dots: HR strain. The dotted ellipse of score plots describe the 95% confidence interval of the Hotelling’s T-squared distribution. The color bar in the S-line plot corresponds to the absolute value of the correlation loading in the discrimination model. A total of 9 observations are shown distributed in 3 independent biological replicates per experimental condition (WT, MR, and HR).
Cells 10 01063 g004
Cells 10 01063 g005 550
Figure 5. Multivariate analysis showing significant differences in metabolomic profiles of sensitive and resistant Leishmania parasites at the extracellular level by 1H-NMR. (A) Score plot of PCA. (B) OPLS-DA score plot comparing antimony-resistant vs. antimony-sensitive parasites. (C) S-line plot highlighting the metabolites with major contributions to group separation. Each individual dot in panel (A) or (B) represents an observation or sample. Some dots can be overlapped. Blue dots: WT strain; green dots: MR strain; red dots: HR strain. The dotted ellipse of score plots describe the 95% confidence interval of the Hotelling’s T-squared distribution. The color bar in the S-line plot corresponds to the absolute value of the correlation loading in the discrimination model. A total of 20 observations or independent biological replicates are shown, 8 of them correspond to WT, 6 to MR, and 6 to HR group.
Figure 5. Multivariate analysis showing significant differences in metabolomic profiles of sensitive and resistant Leishmania parasites at the extracellular level by 1H-NMR. (A) Score plot of PCA. (B) OPLS-DA score plot comparing antimony-resistant vs. antimony-sensitive parasites. (C) S-line plot highlighting the metabolites with major contributions to group separation. Each individual dot in panel (A) or (B) represents an observation or sample. Some dots can be overlapped. Blue dots: WT strain; green dots: MR strain; red dots: HR strain. The dotted ellipse of score plots describe the 95% confidence interval of the Hotelling’s T-squared distribution. The color bar in the S-line plot corresponds to the absolute value of the correlation loading in the discrimination model. A total of 20 observations or independent biological replicates are shown, 8 of them correspond to WT, 6 to MR, and 6 to HR group.
Cells 10 01063 g005
Cells 10 01063 g006 550
Figure 6. Proline, glycerol-phosphocholine, and betaine as potential biomarkers of antimony resistance at the intracellular level. (A) Bar plot of metabolite concentration as the raw intensity values. (B) Box plot of metabolite concentration as the normalized intensity values and showing the Pearson correlation scores. (C) Univariate ROC analysis comparing SbIII-resistant (MR and HR) vs. -sensitive strains (WT). Coefficient of determination by Pearson correlation (R2). False discovery rate of the p-values obtained for the R2 score (FDR). Significant high correlations (FDR ≤ 0.05) are highlighted with a yellow star. Sensitivity or true positive rate (TPR). Specificity or false positive rate (FPR). Area under the curve for univariate ROC analysis (AUC).
Figure 6. Proline, glycerol-phosphocholine, and betaine as potential biomarkers of antimony resistance at the intracellular level. (A) Bar plot of metabolite concentration as the raw intensity values. (B) Box plot of metabolite concentration as the normalized intensity values and showing the Pearson correlation scores. (C) Univariate ROC analysis comparing SbIII-resistant (MR and HR) vs. -sensitive strains (WT). Coefficient of determination by Pearson correlation (R2). False discovery rate of the p-values obtained for the R2 score (FDR). Significant high correlations (FDR ≤ 0.05) are highlighted with a yellow star. Sensitivity or true positive rate (TPR). Specificity or false positive rate (FPR). Area under the curve for univariate ROC analysis (AUC).
Cells 10 01063 g006
Cells 10 01063 g007 550
Figure 7. Lactate, proline, and threonine as potential biomarkers of antimony resistance at the extracellular level. (A) Bar plot of metabolite concentration as the raw intensity values. (B) Box plot of metabolite concentration as the normalized intensity values and showing the Pearson correlation scores. (C) Univariate ROC analysis comparing SbIII-resistant (MR and HR) vs. -sensitive strains (WT). Coefficient of determination by Pearson correlation (R2). False discovery rate of the p-values obtained for the R2 score (FDR). Significant high correlations (FDR ≤ 0.05) are highlighted with a yellow star. Sensitivity or true positive rate (TPR). Specificity or false positive rate (FPR). Area under the curve for univariate ROC analysis (AUC).
Figure 7. Lactate, proline, and threonine as potential biomarkers of antimony resistance at the extracellular level. (A) Bar plot of metabolite concentration as the raw intensity values. (B) Box plot of metabolite concentration as the normalized intensity values and showing the Pearson correlation scores. (C) Univariate ROC analysis comparing SbIII-resistant (MR and HR) vs. -sensitive strains (WT). Coefficient of determination by Pearson correlation (R2). False discovery rate of the p-values obtained for the R2 score (FDR). Significant high correlations (FDR ≤ 0.05) are highlighted with a yellow star. Sensitivity or true positive rate (TPR). Specificity or false positive rate (FPR). Area under the curve for univariate ROC analysis (AUC).
Cells 10 01063 g007
Cells 10 01063 g008 550
Figure 8. SbIII-resistant parasites showed a coordinated metabolomic remodeling to combat the stress caused by the drug. Yellow boxes highlight some of the effects typically caused by SbIII in Leishmania parasites: glycolysis and β-oxidation inhibition, osmotic and oxidative stress. Green boxes highlight the group of metabolites contributing to the osmotic/oxidative balance and thiol metabolism activation. Blue boxes show the metabolites potentially contributing to the energy metabolism required to fuel thiol metabolism in the form of NADPH. This coordinated metabolic response is essential to develop resistance and combat the drug efficiently in Leishmania parasites. Metabolites increased in SbIII-resistant parasites (↑). Metabolites decreased in SbIII-resistant parasites (↓). Metabolites differentially detected at the intracellular (In) and extracellular levels (Ex).
Figure 8. SbIII-resistant parasites showed a coordinated metabolomic remodeling to combat the stress caused by the drug. Yellow boxes highlight some of the effects typically caused by SbIII in Leishmania parasites: glycolysis and β-oxidation inhibition, osmotic and oxidative stress. Green boxes highlight the group of metabolites contributing to the osmotic/oxidative balance and thiol metabolism activation. Blue boxes show the metabolites potentially contributing to the energy metabolism required to fuel thiol metabolism in the form of NADPH. This coordinated metabolic response is essential to develop resistance and combat the drug efficiently in Leishmania parasites. Metabolites increased in SbIII-resistant parasites (↑). Metabolites decreased in SbIII-resistant parasites (↓). Metabolites differentially detected at the intracellular (In) and extracellular levels (Ex).
Cells 10 01063 g008
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gutierrez Guarnizo, S.A.; Karamysheva, Z.N.; Galeano, E.; Muskus, C.E. Metabolite Biomarkers of Leishmania Antimony Resistance. Cells 2021, 10, 1063. https://0-doi-org.brum.beds.ac.uk/10.3390/cells10051063

AMA Style

Gutierrez Guarnizo SA, Karamysheva ZN, Galeano E, Muskus CE. Metabolite Biomarkers of Leishmania Antimony Resistance. Cells. 2021; 10(5):1063. https://0-doi-org.brum.beds.ac.uk/10.3390/cells10051063

Chicago/Turabian Style

Gutierrez Guarnizo, Sneider Alexander, Zemfira N. Karamysheva, Elkin Galeano, and Carlos E. Muskus. 2021. "Metabolite Biomarkers of Leishmania Antimony Resistance" Cells 10, no. 5: 1063. https://0-doi-org.brum.beds.ac.uk/10.3390/cells10051063

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop