Major changes in chromosomal somy, gene expression and gene dosage driven by SbIII in Leishmania braziliensis and Leishmania panamensis

Insights into the trypanothione system in antimony-resistant and sensitive Leishmania tropica clinical isolates

Pentavalent antimonials are the mainstay treatment against different clinical forms of leishmaniasis. The emergence of resistant isolates in endemic areas has led to treatment failure. Unraveling the underlying resistance mechanism would assist in improving the treatment strategies against resistant isolates. This study aimed to investigate the RNA expression level of glutathione synthetase (GS), Spermidine synthetase (SpS), trypanothione synthetase (TryS) genes involved in trypanothione synthesis, and thiol-dependent reductase (TDR) implicated in drug reduction, in antimony-sensitive and -resistant Leishmania tropica isolates. We investigated 11 antimony-resistant and 11 antimony-sensitive L. tropica clinical isolates from ACL patients. Drug sensitivity of amastigotes was determined in mouse macrophage cell line J774A.1. The RNA expression level in the promastigote forms was analyzed by quantitative real-time PCR. The results revealed a significant increase in the average expression of GS, SpS, and TrpS genes by 2.19, 1.56, and 2.33-fold in resistant isolates compared to sensitive ones. The average expression of TDR was 1.24-fold higher in resistant isolates, which was insignificant. The highest correlation coefficient between inhibitory concentration (IC50) values and gene expression belonged to the TryS, GS, SpS, and TDR genes. Moreover, the intracellular thiol content was increased 2.17-fold in resistant isolates compared to sensitive ones and positively correlated with IC50 values. Our findings suggest that overexpression of trypanothione biosynthesis genes and increased thiol content might play a key role in the antimony resistance of L. tropica clinical isolates. In addition, the diversity of gene expression in the trypanothione system and thiol content among L. tropica clinical isolates highlighted the phenotypic heterogeneity of antimony resistance among the parasite population.

Antimony resistance and gene expression in Leishmania: spotlight on molecular and proteomic aspects

Leishmaniasis is a vector-borne parasitic disease caused by Leishmania parasites with a spectrum of clinical manifestations, ranging from skin lesions to severe visceral complications. Treatment of this infection has been extremely challenging with the concurrent emergence of drug resistance. The differential gene expression and the discrepancies in protein functions contribute to the appearance of 2 distinct phenotypes: resistant and sensitive, but the current diagnostic tools fail to differentiate between them. The identification of gene expression patterns and molecular mechanisms coupled with antimony (Sb) resistance can be leveraged to prompt diagnosis and select the most effective treatment methods. The present study attempts to use comparative expression of Sb resistance-associated genes in resistant and sensitive Leishmania, to disclose their relative abundance in clinical or in vitro selected isolates to gain an understanding of the molecular mechanisms of Sb response/resistance. Data suggest that the analysis of resistance gene expression would verify the Sb resistance or susceptibility only to a certain extent; however, none of the individual expression patterns of the studied genes was diagnostic as a biomarker of Sb response of Leishmania. The findings highlighted will be useful in bridging the knowledge gap and discovering innovative diagnostic tools and novel therapeutic targets.

The adaptive roles of aneuploidy and polyclonality in Leishmania in response to environmental stress

Aneuploidy is generally considered harmful, but in some microorganisms, it can act as an adaptive mechanism against environmental stress. Here, we use Leishmania-a protozoan parasite with remarkable genome plasticity-to study the early steps of aneuploidy evolution under high drug pressure (using antimony or miltefosine as stressors). By combining single-cell genomics, lineage tracing with cellular barcodes, and longitudinal genome characterization, we reveal that aneuploidy changes under antimony pressure result from polyclonal selection of pre-existing karyotypes, complemented by further and rapid de novo alterations in chromosome copy number along evolution. In the case of miltefosine, early parasite adaptation is associated with independent point mutations in a miltefosine transporter gene, while aneuploidy changes only emerge later, upon exposure to increased drug levels. Therefore, polyclonality and genome plasticity are hallmarks of parasite adaptation, but the scenario of aneuploidy dynamics depends on the nature and strength of the environmental stress as well as on the existence of other pre-adaptive mechanisms.

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

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.

Life in plastic, it's fantastic! How Leishmania exploit genome instability to shape gene expression

Leishmania are kinetoplastid pathogens that cause leishmaniasis, a debilitating and potentially life-threatening infection if untreated. Unusually, Leishmania regulate their gene expression largely post-transcriptionally due to the arrangement of their coding genes into polycistronic transcription units that may contain 100s of functionally unrelated genes. Yet, Leishmania are capable of rapid and responsive changes in gene expression to challenging environments, often instead correlating with dynamic changes in their genome composition, ranging from chromosome and gene copy number variations to the generation of extrachromosomal DNA and the accumulation of point mutations. Typically, such events indicate genome instability in other eukaryotes, coinciding with genetic abnormalities, but for Leishmania, exploiting these products of genome instability can provide selectable substrates to catalyse necessary gene expression changes by modifying gene copy number. Unorthodox DNA replication, DNA repair, replication stress factors and DNA repeats are recognised in Leishmania as contributors to this intrinsic instability, but how Leishmania regulate genome plasticity to enhance fitness whilst limiting toxic under- or over-expression of co-amplified and co-transcribed genes is unclear. Herein, we focus on fresh, and detailed insights that improve our understanding of genome plasticity in Leishmania. Furthermore, we discuss emerging models and factors that potentially circumvent regulatory issues arising from polycistronic transcription. Lastly, we highlight key gaps in our understanding of Leishmania genome plasticity and discuss future studies to define, in higher resolution, these complex regulatory interactions.

Genome-wide analysis reveals allelic variation and chromosome copy number variation in paromomycin-resistant Leishmania donovani

In the absence of adequate diagnosis and treatment, leishmaniasis remains a major public health concern on a global scale. Drug resistance remains a key obstacle in controlling and eliminating visceral leishmaniasis. The therapeutic gap due to lack of target-specific medicine and vaccine can be minimized by obtaining parasite's genomic information. This study compared whole-genome sequence of paromomycin-resistant parasite (K133PMM) developed through in vitro adaptation and selection with sensitive Leishmania clinical isolate (K133WT). We found a large number of upstream and intergenic gene variations in K133PMM. There were 259 single nucleotide polymorphisms (SNPs), 187 insertion-deletion (InDels), and 546 copy number variations (CNVs) identified. Most of the genomic variations were found in the gene's upstream and non-coding regions. Ploidy estimation revealed chromosome 5 in tetrasomy and 6, 9, and 12 in trisomy, uniquely in K133PMM. These contain the genes for protein degradation, parasite motility, autophagy, cell cycle maintenance, and drug efflux membrane transporters. Furthermore, we also observed reduction in ploidy of chromosomes 15, 20, and 23, in the resistant parasite containing mostly the genes for hypothetical proteins and membrane transporters. We chronicled correlated genomic conversion and aneuploidy in parasites and hypothesize that this led to rapid evolutionary changes in response to drug induced pressure, which causes them to become resistant.

Transcriptome Analysis of Intracellular Amastigotes of Clinical Leishmania infantum Lines from Therapeutic Failure Patients after Infection of Human Macrophages

Leishmaniasis is considered to be one of the most neglected tropical diseases affecting humans and animals around the world. Due to the absence of an effective vaccine, current treatment is based on chemotherapy. However, the continuous appearance of drug resistance and therapeutic failure (TF) lead to an early obsolescence of treatments. Identification of the factors that contribute to TF and drug resistance in leishmaniasis will constitute a useful tool for establishing future strategies to control this disease. In this manuscript, we evaluated the transcriptomic changes in the intracellular amastigotes of the Leishmania infantum parasites isolated from patients with leishmaniasis and TF at 96 h post-infection of THP-1 cells. The adaptation of the parasites to their new environment leads to expression alterations in the genes involved mainly in the transport through cell membranes, energy and redox metabolism, and detoxification. Specifically, the gene that codes for the prostaglandin f2α synthase seems to be relevant in the pathogenicity and TF since it appears substantially upregulated in all the L. infantum lines. Overall, our results show that at the late infection timepoint, the transcriptome of the parasites undergoes significant changes that probably improve the survival of the Leishmania lines in the host cells, contributing to the TF phenotype as well as drug therapy evasion.

Genome plasticity driven by aneuploidy and loss of heterozygosity in Trypanosoma cruzi

Trypanosoma cruzi the causative agent of Chagas disease shows a marked genetic diversity and divided into at least six Discrete Typing Units (DTUs). High intra genetic variability has been observed in the TcI DTU, the most widely distributed DTU, where patterns of genomic diversity can provide information on ecological and evolutionary processes driving parasite population structure and genome organization. Chromosomal aneuploidies and rearrangements across multigene families represent an evidence of T. cruzi genome plasticity. We explored genomic diversity among 18 Colombian T. cruzi I clones and 15 T. cruzi I South American strains. Our results confirm high genomic variability, heterozygosity and presence of a clade compatible with the TcI_dom genotype, described for strains from humans in Colombia and Venezuela. TcI showed high structural plasticity across the geographical region studied. Differential events of whole and segmental aneuploidy (SA) along chromosomes even between clones from the same strain were found and corroborated by the depth and allelic frequency. We detected loss of heterozygosity (LOH) events in different chromosomes, however, the size and location of segments under LOH varied between clones. Genes adjacent to breakpoints were evaluated, and retrotransposon hot spot genes flanked the beginning of segmental aneuploidies. Our results suggest that T. cruzi genomes, like those of Leishmania, may have a highly unstable structure and there is now an urgent need to design experiments to explore any potential adaptive role for the plasticity observed.

Characterization of Differentially Abundant Proteins Among Leishmania (Viannia) braziliensis Strains Isolated From Atypical or Typical Lesions

Leishmania (Viannia) braziliensis is the main etiological agent of cutaneous and mucocutaneous leishmaniasis in Latin America. Non-ulcerated atypical tegumentary leishmaniasis cases caused by L. braziliensis have been reported in several regions of the American continent, including the Xacriabá indigenous reserve in São João das Missões/Minas Gerais, Brazil. Parasites isolated from these atypical clinical lesions are resistant to antimony-based therapeutics. In the present study, proteins displaying differential abundance in two strains of L. braziliensis isolated from patients with atypical lesions compared with four strains isolated from patients with typical lesions were identified using a quantitative proteomics approach based on tandem mass tag labeling (TMT) and mass spectrometry. A total of 532 (P<0.05) differentially abundant proteins were identified (298 upregulated and 234 downregulated) in strains from atypical lesions compared to strains from typical lesions. Prominent positively regulated proteins in atypical strains included those that may confer greater survival inside macrophages, proteins related to antimony resistance, and proteins associated with higher peroxidase activity. Additionally, we identified proteins showing potential as new drug and vaccine targets. Our findings contribute to the characterization of these intriguing L. braziliensis strains and provide a novel perspective on Atypical Cutaneous Leishmaniasis (ACL) cases that have been associated with therapeutic failures.

Specific Human ATR and ATM Inhibitors Modulate Single Strand DNA Formation in Leishmania major Exposed to Oxidative Agent

Leishmania parasites are the causative agents of a group of neglected tropical diseases known as leishmaniasis. The molecular mechanisms employed by these parasites to adapt to the adverse conditions found in their hosts are not yet completely understood. DNA repair pathways can be used by Leishmania to enable survival in the interior of macrophages, where the parasite is constantly exposed to oxygen reactive species. In higher eukaryotes, DNA repair pathways are coordinated by the central protein kinases ataxia telangiectasia mutated (ATM) and ataxia telangiectasia and Rad3 related (ATR). The enzyme Exonuclease-1 (EXO1) plays important roles in DNA replication, repair, and recombination, and it can be regulated by ATM- and ATR-mediated signaling pathways. In this study, the DNA damage response pathways in promastigote forms of L. major were investigated using bioinformatics tools, exposure of lineages to oxidizing agents and radiation damage, treatment of cells with ATM and ATR inhibitors, and flow cytometry analysis. We demonstrated high structural and important residue conservation for the catalytic activity of the putative LmjEXO1. The overexpression of putative LmjEXO1 made L. major cells more susceptible to genotoxic damage, most likely due to the nuclease activity of this enzyme and the occurrence of hyper-resection of DNA strands. These cells could be rescued by the addition of caffeine or a selective ATM inhibitor. In contrast, ATR-specific inhibition made the control cells more susceptible to oxidative damage in an LmjEXO1 overexpression-like manner. We demonstrated that ATR-specific inhibition results in the formation of extended single-stranded DNA, most likely due to EXO1 nucleasic activity. Antagonistically, ATM inhibition prevented single-strand DNA formation, which could explain the survival phenotype of lineages overexpressing LmjEXO1. These results suggest that an ATM homolog in Leishmania could act to promote end resection by putative LmjEXO1, and an ATR homologue could prevent hyper-resection, ensuring adequate repair of the parasite DNA.

Nitric Oxide Resistance in Leishmania (Viannia) braziliensis Involves Regulation of Glucose Consumption, Glutathione Metabolism and Abundance of Pentose Phosphate Pathway Enzymes

In American Tegumentary Leishmaniasis production of cytokines, reactive oxygen species and nitric oxide (NO) by host macrophages normally lead to parasite death. However, some Leishmania braziliensis strains exhibit natural NO resistance. NO-resistant strains cause more lesions and are frequently more resistant to antimonial treatment than NO-susceptible ones, suggesting that NO-resistant parasites are endowed with specific mechanisms of survival and persistence. To tests this, we analyzed the effect of pro- and antioxidant molecules on the infectivity in vitro of L. braziliensis strains exhibiting polar phenotypes of resistance or susceptibility to NO. In addition, we conducted a comprehensive quantitative mass spectrometry-based proteomics analysis of those parasites. NO-resistant parasites were more infective to peritoneal macrophages, even in the presence of high levels of reactive species. Principal component analysis of protein concentration values clearly differentiated NO-resistant from NO-susceptible parasites, suggesting that there are natural intrinsic differences at molecular level among those strains. Upon NO exposure, NO-resistant parasites rapidly modulated their proteome, increasing their total protein content and glutathione (GSH) metabolism. Furthermore, NO-resistant parasites showed increased glucose analogue uptake, and increased abundance of phosphotransferase and G6PDH after nitrosative challenge, which can contribute to NADPH pool maintenance and fuel the reducing conditions for the recovery of GSH upon NO exposure. Thus, increased glucose consumption and GSH-mediated redox capability may explain the natural resistance of L. braziliensis against NO.

Impact of Genetic Diversity and Genome Plasticity of Leishmania spp. in Treatment and the Search for Novel Chemotherapeutic Targets

Leishmaniasis is one of the major public health concerns in Latin America, Africa, Asia, and Europe. The absence of vaccines for human use and the lack of effective vector control programs make chemotherapy the main strategy to control all forms of the disease. However, the high toxicity of available drugs, limited choice of therapeutic agents, and occurrence of drug-resistant parasite strains are the main challenges related to chemotherapy. Currently, only a small number of drugs are available for leishmaniasis treatment, including pentavalent antimonials (Sb^V), amphotericin B and its formulations, miltefosine, paromomycin sulphate, and pentamidine isethionate. In addition to drug toxicity, therapeutic failure of leishmaniasis is a serious concern. The occurrence of drug-resistant parasites is one of the causes of therapeutic failure and is closely related to the diversity of parasites in this genus. Owing to the enormous plasticity of the genome, resistance can occur by altering different metabolic pathways, demonstrating that resistance mechanisms are multifactorial and extremely complex. Genetic variability and genome plasticity cause not only the available drugs to have limitations, but also make the search for new drugs challenging. Here, we examined the biological characteristics of parasites that hinder drug discovery.

Development of an Amplicon-Based Next-Generation Sequencing Protocol to Identify Leishmania Species and Other Trypanosomatids in Leishmaniasis Endemic Areas

Trypanosomatid infections are an important public health threat affecting many low-income countries across the tropics, particularly in the Americas. Trypanosomatids can infect many vertebrate, invertebrate, and plant species and play an important role as human pathogens. Among these clinically relevant pathogens are species from the genera Leishmania and Trypanosoma. Mixed trypanosomatid infections remain a largely unexplored phenomenon. Herein, we describe the application of an amplicon-based next-generation sequencing (NGS) assay to detect and identify trypanosomatid species in mammalian reservoirs, human patients, and sand fly vectors throughout regions of Leishmania endemicity. Sixty-five samples from different departments of Colombia, including two samples from Venezuela, were analyzed: 49 samples from cutaneous leishmaniasis (CL) patients, 8 from sand flies, 2 from domestic reservoirs (Canis familiaris), and 6 from wild reservoirs (Phyllostomus hastatus). DNA from each sample served to identify the presence of trypanosomatids through conventional PCR using heat shock protein 70 (HSP70) gene as the target. PCR products underwent sequencing by Sanger sequencing and NGS, and trypanosomatid species were identified by using BLASTn against a reference database built from trypanosomatid-derived HSP70 sequences. The alpha and beta diversity indexes of amplicon sequence variants were calculated for each group. The results revealed the presence of mixed infections with more than two Leishmania species in 34% of CL samples analyzed. Trypanosoma cruzi was identified in samples from wild reservoirs, as well as in sand fly vectors. Coinfection events with three different Leishmania species were identified in domestic reservoirs. These findings depose the traditional paradigm of leishmaniasis as being a single-species-driven infection and redraw the choreography of host-pathogen interaction in the context of multiparasitism. Further research is needed to decipher how coinfections may influence disease progression. This knowledge is key to developing an integrated approach for diagnosis and treatment.

IMPORTANCE Traditionally, there has been a frequent, yet incorrect assumption that phlebotomine vectors, animal reservoirs, and human hosts are susceptible to Leishmania infection by a single parasite species. However, current evidence supports that these new vector-parasite-reservoir associations lend vectors and reservoirs greater permissiveness to certain Leishmania species, thus promoting the appearance of coinfection events, particularly in disease-endemic regions. The present study describes the application of an amplicon-based next-generation sequencing (NGS) assay to detect and identify trypanosomatid species in mammalian reservoirs, human patients, and sand fly vectors from regions of endemicity for leishmaniasis. This changes our understanding of the clinical course of leishmaniasis in areas of endemicity.

Insights from Leishmania (Viannia) guyanensis in vitro behavior and intercellular communication

BACKGROUND

Pentavalent antimonial-based chemotherapy is the first-line approach for leishmaniasis treatment and disease control. Nevertheless antimony-resistant parasites have been reported in some endemic regions. Treatment refractoriness is complex and is associated with patient- and parasite-related variables. Although amastigotes are the parasite stage in the vertebrate host and, thus, exposed to the drug, the stress caused by trivalent antimony in promastigotes has been shown to promote significant modification in expression of several genes involved in various biological processes, which will ultimately affect parasite behavior. Leishmania (Viannia) guyanensis is one of the main etiological agents in the Amazon Basin region, with a high relapse rate (approximately 25%).

METHODS

Herein, we conducted several in vitro analyses with L. (V.) guyanensis strains derived from cured and refractory patients after treatment with standardized antimonial therapeutic schemes, in addition to a drug-resistant in vitro-selected strain. Drug sensitivity assessed through Sb(III) half-maximal inhibitory concentration (IC₅₀) assays, growth patterns (with and without drug pressure) and metacyclic-like percentages were determined for all strains and compared to treatment outcomes. Finally, co-cultivation without intercellular contact was followed by parasitic density and Sb(III) IC₅₀ measurements.

RESULTS

Poor treatment response was correlated with increased Sb(III) IC₅₀ values. The decrease in drug sensitivity was associated with a reduced cell replication rate, increased in vitro growth ability, and higher metacyclic-like proportion. Additionally, in vitro co-cultivation assays demonstrated that intercellular communication enabled lower drug sensitivity and enhanced in vitro growth ability, regardless of direct cell contact.

CONCLUSIONS

Data concerning drug sensitivity in the Viannia subgenus are emerging, and L. (V.) guyanensis plays a pivotal epidemiological role in Latin America. Therefore, investigating the parasitic features potentially related to relapses is urgent. Altogether, the data presented here indicate that all tested strains of L. (V.) guyanensis displayed an association between treatment outcome and in vitro parameters, especially the drug sensitivity. Remarkably, sharing enhanced growth ability and decreased drug sensitivity, without intercellular communication, were demonstrated.

Advances in Understanding Leishmania Pathobiology: What Does RNA-Seq Tell Us?

Leishmaniasis is a vector-borne disease caused by a protozoa parasite from over 20 Leishmania species. The clinical manifestations and the outcome of the disease vary greatly. Global RNA sequencing (RNA-Seq) analyses emerged as a powerful technique to profile the changes in the transcriptome that occur in the Leishmania parasites and their infected host cells as the parasites progresses through their life cycle. Following the bite of a sandfly vector, Leishmania are transmitted to a mammalian host where neutrophils and macrophages are key cells mediating the interactions with the parasites and result in either the elimination the infection or contributing to its proliferation. This review focuses on RNA-Seq based transcriptomics analyses and summarizes the main findings derived from this technology. In doing so, we will highlight caveats in our understanding of the parasite’s pathobiology and suggest novel directions for research, including integrating more recent data highlighting the role of the bacterial members of the sandfly gut microbiota and the mammalian host skin microbiota in their potential role in influencing the quantitative and qualitative aspects of leishmaniasis pathology.

Revisiting the heterogeneous global genomic population structure of Leishmania infantum

Leishmania infantum is the main causative agent responsible for visceral leishmaniasis (VL), a disease with global distribution. The genomic structure and genetic variation of this species have been widely studied in different parts of the world. However, in some countries, this information is still yet unknown, as is the genomic behaviour of the main antigens used in VL diagnosis (rK39 and rK28), which have demonstrated variable sensitivity and specificity in a manner dependent on the geographic region analysed. The objective of this study was to explore the genomic architecture and diversity of four Colombian L. infantum isolates obtained in this study and to compare these results with the genetic analysis of 183 L. infantum isolates from across the world (obtained from public databases), as well as to analyse the whole rK39 and rK28 antigen sequences in our dataset. The results showed that, at the global level, L. infantum has high genetic homogeneity and extensive aneuploidy. Furthermore, we demonstrated that there are distinct populations of L. infantum circulating in various countries throughout the globe and that populations of distant countries have close genomic relationships. Additionally, this study demonstrated the high genetic variability of the rK28 antigen worldwide. In conclusion, our study allowed us to (i) expand our knowledge of the genomic structure of global L. infantum; (ii) describe the intra-specific genomic variability of this species; and (iii) understand the genomic characteristics of the main antigens used in the diagnosis of VL. Additionally, this is the first study to report whole-genome sequences of Colombian L. infantum isolates.

Three different mutations in the DNA topoisomerase 1B in Leishmania infantum contribute to resistance to antitumor drug topotecan

Background

The evolution of drug resistance is one of the biggest challenges in leishmaniasis and has prompted the need for new antileishmanial drugs. Repurposing of approved drugs is a faster and very attractive strategy that is gaining supporters worldwide. Different anticancer topoisomerase 1B (TOP1B) inhibitors have shown strong antileishmanial activity and promising selective indices, supporting the potential repurposing of these drugs. However, cancer cells and Leishmania share the ability to become rapidly resistant. The aim of this study was to complete a whole-genome exploration of the effects caused by exposure to topotecan in order to highlight the potential mechanisms deployed by Leishmania to favor its survival in the presence of a TOP1B inhibitor.

Methods

We used a combination of stepwise drug resistance selection, whole-genome sequencing, functional validation, and theoretical approaches to explore the propensity of and potential mechanisms deployed by three independent clones of L. infantum to resist the action of TOP1B inhibitor topotecan.

Results

We demonstrated that L. infantum is capable of becoming resistant to high concentrations of topotecan without impaired growth ability. No gene deletions or amplifications were identified from the next-generation sequencing data in any of the three resistant lines, ruling out the overexpression of efflux pumps as the preferred mechanism of topotecan resistance. We identified three different mutations in the large subunit of the leishmanial TOP1B (Top1B^F187Y, Top1B^G191A, and Top1B^W232R). Overexpression of these mutated alleles in the wild-type background led to high levels of resistance to topotecan. Computational molecular dynamics simulations, in both covalent and non-covalent complexes, showed that these mutations have an effect on the arrangement of the catalytic pentad and on the interaction of these residues with surrounding amino acids and DNA. This altered architecture of the binding pocket results in decreased persistence of topotecan in the ternary complex.

Conclusions

This work helps elucidate the previously unclear potential mechanisms of topotecan resistance in Leishmania by mutations in the large subunit of TOP1B and provides a valuable clue for the design of improved inhibitors to combat resistance in both leishmaniasis and cancer. Our data highlights the importance of including drug resistance evaluation in drug discovery cascades.

Graphical abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-021-04947-4.

Comparative analysis of the transcriptional responses of five Leishmania species to trivalent antimony

BACKGROUND

Leishmaniasis is a neglected tropical disease caused by several species of Leishmania. The resistance phenotype of these parasites depends on the characteristics of each species, which contributes to increased therapeutic failures. Understanding the mechanism used by the parasite to survive under treatment pressure in order to identify potential common and specific therapeutic targets is essential for the control of leishmaniasis. The aim of this study was to investigate the expression profiles and potential shared and specific resistance markers of the main Leishmania species of medical importance [subgenus L. (Leishmania): L. donovani, L. infantum and L. amazonensis; subgenus L. (Viannia): L. panamensis and L. braziliensis)] resistant and sensitive to trivalent stibogluconate (SbIII).

METHODS

We conducted comparative analysis of the transcriptomic profiles (only coding sequences) of lines with experimentally induced resistance to SbIII from biological replicates of five Leishmania species available in the databases of four articles based on ortholog attribution. Simultaneously, we carried out functional analysis of ontology and reconstruction of metabolic pathways of the resulting differentially expressed genes (DEGs).

RESULTS

Resistant lines for each species had differential responses in metabolic processes, compound binding, and membrane components concerning their sensitive counterpart. One hundred and thirty-nine metabolic pathways were found, with the three main pathways comprising cysteine and methionine metabolism, glycolysis, and the ribosome. Differentially expressed orthologous genes assigned to species-specific responses predominated, with 899 self-genes. No differentially expressed genes were found in common among the five species. Two common upregulated orthologous genes were found among four species (L. donovani, L. braziliensis, L. amazonensis, and L. panamensis) related to an RNA-binding protein and the NAD(P)H cytochrome-B5-oxidoreductase complex, associated with transcriptional control and de novo synthesis of linoleic acid, critical mechanisms in resistance to antimonials.

CONCLUSION

Herein, we identified potential species-specific genes related to resistance to SbIII. Therefore, we suggest that future studies consider a treatment scheme that is species-specific. Despite the limitations of our study, this is the first approach toward unraveling the pan-genus genetic mechanisms of resistance in leishmaniasis.

DNA Double-Strand Breaks: A Double-Edged Sword for Trypanosomatids

For nearly all eukaryotic cells, stochastic DNA double-strand breaks (DSBs) are one of the most deleterious types of DNA lesions. DSB processing and repair can cause sequence deletions, loss of heterozygosity, and chromosome rearrangements resulting in cell death or carcinogenesis. However, trypanosomatids (single-celled eukaryotes parasites) do not seem to follow this premise strictly. Several studies have shown that trypanosomatids depend on DSBs to perform several events of paramount importance during their life cycle. For Trypanosoma brucei, DSBs formation is associated with host immune evasion via antigenic variation. In Trypanosoma cruzi, DSBs play a crucial role in the genetic exchange, a mechanism that is still little explored but appear to be of fundamental importance for generating variability. In Leishmania spp., DSBs are necessary to generate genomic changes by gene copy number variation (CNVs), events that are essential for these organisms to overcome inhospitable conditions. As DSB repair in trypanosomatids is primarily conducted via homologous recombination (HR), most of the events associated with DSBs are HR-dependent. This review will discuss the latest findings on how trypanosomatids balance the benefits and inexorable challenges caused by DSBs.

Distribution, treatment outcome and genetic diversity of Leishmania species in military personnel from Colombia with cutaneous leishmaniasis

BACKGROUND

Leishmaniasis is one of the most important infectious diseases affecting the Colombian National Army due to the high number of reported cases and exposure throughout military operations in endemic areas. The main aim of this study was to estimate the geographical distribution along with the genetic diversity and treatment outcome of Leishmania species in Colombian military personnel.

METHODS

Skin lesion samples by smear and aspirate were collected in 136 patients having parasitological cutaneous leishmaniasis (CL) diagnosis. DNA was extracted, the nuclear marker heat shock protein 70 (HSP70) was amplified by PCR and sequenced. Leishmania species were identified by BLASTn. The geo-spatial distribution of the identified parasites was determined according to the possible site of infection. Gene tree was constructed by maximum likelihood (ML), diversity indices (π, h) were estimated and haplotype network was constructed under the Templeton-Crandall-Sing algorithm in order to determine the geographic relationships of the genetic variants of Leishmania species circulating in Colombian military population.

RESULTS

The species were identified in 77.94% of the samples, with a predominance of L. braziliensis (65.09%), followed by L. panamensis (31.13%), L. naiffi by the first time reported in Colombia in two patients (1.89%) as well as L. lindenbergi in a single patient (0.945%) with possible infection in the municipality of Miraflores, Guaviare and L. infantum in a single patient (0.945%) notified with CL in the municipality of Tumaco, Nariño. The phylogenetic analysis was consistent according to bootstrap, showing four strongly differentiated clades.

CONCLUSIONS

The geo-spatial distribution suggested that L. braziliensis has a greater abundance, while L. panamensis has a greater dispersion. The phylogenetic relationships of Leishmania species in Colombian military personnel was estimated with the confirmation of two new species circulating without prior report in the country and a species with no background for CL in the Colombian army. A substantial genetic diversity of Leishmania braziliensis was defined. This study contributes through the understanding of the molecular epidemiology to the CL transmission in Colombia.

Comparative transcriptomic analysis of antimony resistant and susceptible Leishmania infantum lines

Background

One of the major challenges to leishmaniasis treatment is the emergence of parasites resistant to antimony. To study differentially expressed genes associated with drug resistance, we performed a comparative transcriptomic analysis between wild-type and potassium antimonyl tartrate (Sb^III)-resistant Leishmania infantum lines using high-throughput RNA sequencing.

Methods

All the cDNA libraries were constructed from promastigote forms of each line, sequenced and analyzed using STAR for mapping the reads against the reference genome (L. infantum JPCM5) and DESeq2 for differential expression statistical analyses. All the genes were functionally annotated using sequence similarity search.

Results

The analytical pipeline considering an adjusted p-value < 0.05 and fold change > 2.0 identified 933 transcripts differentially expressed (DE) between wild-type and Sb^III-resistant L. infantum lines. Out of 933 DE transcripts, 504 presented functional annotation and 429 were assigned as hypothetical proteins. A total of 837 transcripts were upregulated and 96 were downregulated in the Sb^III-resistant L. infantum line. Using this DE dataset, the proteins were further grouped in functional classes according to the gene ontology database. The functional enrichment analysis for biological processes showed that the upregulated transcripts in the Sb^III-resistant line are associated with protein phosphorylation, microtubule-based movement, ubiquitination, host–parasite interaction, cellular process and other categories. The downregulated transcripts in the Sb^III-resistant line are assigned in the GO categories: ribonucleoprotein complex, ribosome biogenesis, rRNA processing, nucleosome assembly and translation.

Conclusions

The transcriptomic profile of L. infantum showed a robust set of genes from different metabolic pathways associated with the antimony resistance phenotype in this parasite. Our results address the complex and multifactorial antimony resistance mechanisms in Leishmania, identifying several candidate genes that may be further evaluated as molecular targets for chemotherapy of leishmaniasis.

Genomic Diversification, Structural Plasticity, and Hybridization in Leishmania (Viannia) braziliensis

Leishmania (Viannia) braziliensis is an important Leishmania species circulating in several Central and South American countries. Among Leishmania species circulating in Brazil, Argentina and Colombia, L. braziliensis has the highest genomic variability. However, genomic variability at the whole genome level has been only studied in Brazilian and Peruvian isolates; to date, no Colombian isolates have been studied. Considering that in Colombia, L. braziliensis is a species with great clinical and therapeutic relevance, as well as the role of genetic variability in the epidemiology of leishmaniasis, we analyzed and evaluated intraspecific genomic variability of L. braziliensis from Colombian and Bolivian isolates and compared them with Brazilian isolates. Twenty-one genomes were analyzed, six from Colombian patients, one from a Bolivian patient, and 14 Brazilian isolates downloaded from public databases. The results obtained of Phylogenomic analysis showed the existence of four well-supported clades, which evidenced intraspecific variability. The whole-genome analysis revealed structural variations in the somy, mainly in the Brazilian genomes (clade 1 and clade 3), low copy number variations, and a moderate number of single-nucleotide polymorphisms (SNPs) in all genomes analyzed. Interestingly, the genomes belonging to clades 2 and 3 from Colombia and Brazil, respectively, were characterized by low heterozygosity (~90% of SNP loci were homozygous) and regions suggestive of loss of heterozygosity (LOH). Additionally, we observed the drastic whole genome loss of heterozygosity and possible hybridization events in one genome belonging to clade 4. Unique/shared SNPs between and within the four clades were identified, revealing the importance of some of them in biological processes of L. braziliensis. Our analyses demonstrate high genomic variability of L. braziliensis in different regions of South America, mainly in Colombia and suggest that this species exhibits striking genomic diversity and a capacity of genomic hybridization; additionally, this is the first study to report whole-genome sequences of Colombian L. braziliensis isolates.

Casein kinase 1.2 over expression restores stress resistance to Leishmania donovani HSP23 null mutants

Leishmania donovani is a trypanosomatidic parasite and causes the lethal kala-azar fever, a neglected tropical disease. The Trypanosomatida are devoid of transcriptional gene regulation and rely on gene copy number variations and translational control for their adaption to changing conditions. To survive at mammalian tissue temperatures, L. donovani relies on the small heat shock protein HSP23, the loss of which renders the parasites stress sensitive and impairs their proliferation. Here, we analysed a spontaneous escape mutant with wild type-like in vitro growth. Further selection of this escape strains resulted in a complete reversion of the phenotype. Whole genome sequencing revealed a correlation between stress tolerance and the massive amplification of a six-gene cluster on chromosome 35, with further analysis showing over expression of the casein kinase 1.2 gene as responsible. In vitro phosphorylation experiments established both HSP23 and the related P23 co-chaperone as substrates and modulators of casein kinase 1.2, providing evidence for another crucial link between chaperones and signal transduction protein kinases in this early branching eukaryote.

In-depth quantitative proteomics uncovers specie-specific metabolic programs in Leishmania (Viannia) species

Leishmania species are responsible for a broad spectrum of diseases, denominated Leishmaniasis, affecting over 12 million people worldwide. During the last decade, there have been impressive efforts for sequencing the genome of most of the pathogenic Leishmania spp. as well as hundreds of strains, but large-scale proteomics analyses did not follow these achievements and the Leishmania proteome remained mostly uncharacterized. Here, we report a comprehensive comparative study of the proteomes of strains representing L. braziliensis, L. panamensis and L. guyanensis species. Proteins extracted by SDS-mediated lysis were processed following the multi-enzyme digestion-filter aided sample preparation (FASP) procedure and analysed by high accuracy mass spectrometry. "Total Protein Approach" and "Proteomic Ruler" were applied for absolute quantification of proteins. Principal component analysis demonstrated very high reproducibility among biological replicates and a very clear differentiation of the three species. Our dataset comprises near 7000 proteins, representing the most complete Leishmania proteome yet known, and provides a comprehensive quantitative picture of the proteomes of the three species in terms of protein concentration and copy numbers. Analysis of the abundance of proteins from the major energy metabolic processes allow us to highlight remarkably differences among the species and suggest that these parasites depend on distinct energy substrates to obtain ATP. Whereas L. braziliensis relies the more on glycolysis, L. panamensis and L. guyanensis seem to depend mainly on mitochondrial respiration. These results were confirmed by biochemical assays showing opposite profiles for glucose uptake and O2 consumption in these species. In addition, we provide quantitative data about different membrane proteins, transporters, and lipids, all of which contribute for significant species-specific differences and provide rich substrate for explore new molecules for diagnosing purposes. Data are available via ProteomeXchange with identifier PXD017696.

Application of next generation sequencing (NGS) for descriptive analysis of 30 genomes of Leishmania infantum isolates in Middle-North Brazil

Visceral leishmaniasis (VL) is a life-threatening disease caused by the protozoa Leishmania donovani and L. infantum. Likely, L. infantum was introduced in the New World by the Iberic colonizers. Due to recent introduction, the genetic diversity is low. Access to genomic information through the sequencing of Leishmania isolates allows the characterization of populations through the identification and analysis of variations. Population structure information may reveal important data on disease dynamics. Aiming to describe the genetic diversity of L. infantum from the Middle-North, Brazil, next generation sequencing of 30 Leishmania isolates obtained in the city of Teresina, from where the disease dispersed, was performed. The variations were categorized accordingly to the genome region and impact and provided the basis for chromosomal ploidy and population structure analysis. The results showed low diversity between the isolates and the Iberic reference genome JPCM5. Most variations were seen in non-coding regions, with modifying impact. The ploidy number analysis showed aneuploid profile. The population structure analysis revealed the presence of two L. infantum populations identified in Teresina. Further population genetics studies with a larger number of isolates should be performed in order to identify the genetic background associated with virulence and parasite ecology.

Can We Harness Immune Responses to Improve Drug Treatment in Leishmaniasis?

Leishmaniasis is a vector-borne parasitic disease that has been neglected in priority for control and eradication of malaria, tuberculosis, and HIV/AIDS. Collectively, over one seventh of the world’s population is at risk of being infected with 0.7–1.2 million new infections reported annually. Clinical manifestations range from self-healing cutaneous lesions to fatal visceral disease. The first anti-leishmanial drugs were introduced in the 1950′s and, despite several shortcomings, remain the mainstay for treatment. Regardless of this and the steady increase in infections over the years, particularly among populations of low economic status, research on leishmaniasis remains under funded. This review looks at the drugs currently in clinical use and how they interact with the host immune response. Employing chemoimmunotherapeutic approaches may be one viable alternative to improve the efficacy of novel/existing drugs and extend their lifespan in clinical use.

Experimental Strategies to Explore Drug Action and Resistance in Kinetoplastid Parasites

Kinetoplastids are the causative agents of leishmaniasis, human African trypanosomiasis, and American trypanosomiasis. They are responsible for high mortality and morbidity in (sub)tropical regions. Adequate treatment options are limited and have several drawbacks, such as toxicity, need for parenteral administration, and occurrence of treatment failure and drug resistance. Therefore, there is an urgency for the development of new drugs. Phenotypic screening already allowed the identification of promising new chemical entities with anti-kinetoplastid activity potential, but knowledge on their mode-of-action (MoA) is lacking due to the generally applied whole-cell based approach. However, identification of the drug target is essential to steer further drug discovery and development. Multiple complementary techniques have indeed been used for MoA elucidation. In this review, the different 'omics' approaches employed to define the MoA or mode-of-resistance of current reference drugs and some new anti-kinetoplastid compounds are discussed.

Regulation of Translation in the Protozoan Parasite Leishmania

Leishmaniasis represents a serious health problem worldwide and drug resistance is a growing concern. Leishmania parasites use unusual mechanisms to control their gene expression. In contrast to many other species, they do not have transcriptional regulation. The lack of transcriptional control is mainly compensated by post-transcriptional mechanisms, including tight translational control and regulation of mRNA stability/translatability by RNA-binding proteins. Modulation of translation plays a major role in parasite survival and adaptation to dramatically different environments during change of host; however, our knowledge of fine molecular mechanisms of translation in Leishmania remains limited. Here, we review the current progress in our understanding of how changes in the translational machinery promote parasite differentiation during transmission from a sand fly to a mammalian host, and discuss how translational reprogramming can contribute to the development of drug resistance.

MRPA-independent mechanisms of antimony resistance in Leishmania infantum

Control of both human and canine leishmaniasis is based on a very short list of chemotherapeutic agents, headed by antimonial derivatives (Sb). The utility of these molecules is severely threatened by high rates of drug resistance. The ABC transporter MRPA is one of the few key Sb resistance proteins described to date, whose role in detoxification has been thoroughly studied in Leishmania parasites. Nonetheless, its rapid amplification during drug selection complicates the discovery of other mechanisms potentially involved in Sb resistance. In this study, stepwise drug-resistance selection and next-generation sequencing were combined in the search for novel Sb-resistance mechanisms deployed by parasites when MRPA is abolished by targeted gene disruption. The gene mrpA is not essential in L. infantum, and its disruption leads to an Sb hypersensitive phenotype in both promastigotes and amastigotes. Five independent mrpA^-/- mutants were selected for antimony resistance. These mutants displayed major changes in their ploidy, as well as extrachromosomal linear amplifications of the subtelomeric region of chromosome 23, which includes the genes coding for ABCC1 and ABCC2. Overexpression of ABCC2, but not of ABCC1, resulted in increased Sb tolerance in the mrpA^-/- mutant. SNP analyses revealed three different heterozygous mutations in the gene coding for a serine acetyltransferase (SAT) involved in de novo cysteine synthesis in Leishmania. Overexpression of sat^Q390K, sat^G321R and sat^G325R variants led to a 2-3.2 -fold increase in Sb resistance in mrpA^-/- parasites. Only sat^G321R and sat^G325R induced increased Sb resistance in wild-type parasites. These results reinforce and expand knowledge on the complex nature of Sb resistance in Leishmania parasites.

Intraspecific Genomic Divergence and Minor Structural Variations in Leishmania (Viannia) panamensis

Leishmania (Viannia) panamensis is one of the most important Leishmania species associated with cutaneous leishmaniasis (CL) in Latin America. Despite its wide geographic distribution and pathogenic potential in humans and animals, the genomic variability of this species is low compared with other Leishmania species circulating in the same geographical area. No studies have reported a detailed analysis of the whole genome of L. panamensis from clinical isolates using DNA high-throughput sequencing to clarify its intraspecific genomic variability or plausible divergence. Therefore, this study aimed to evaluate the intraspecific genomic variability of L. panamensis from Colombia and Panama. A total of 22 genomes were analyzed, 19 from Colombian patients with CL and three genomes from Panama obtained from public databases. The phylogenomic analysis revealed the potential existence of three well-supported clades as evidence of intraspecific divergence. Additionally, the whole-genome analysis showed low structural variations in terms of ploidy, copy number variations, and single-nucleotide polymorphisms (SNPs). SNPs shared among all clades were identified, revealing their importance in different biological processes of L. panamensis. The findings not only expand our knowledge of intraspecific genomic variability of one of the most important Leishmania species in South America but also highlights the possible existence of different clades/lineages/subpopulations across a geographic scale.

Trans-Atlantic Spill Over: Deconstructing the Ecological Adaptation of Leishmania infantum in the Americas

Pathogen fitness landscapes change when transmission cycles establish in non-native environments or spill over into new vectors and hosts. The introduction of Leishmania infantum in the Americas into the Neotropics during European colonization represents a unique case study to investigate the mechanisms of ecological adaptation of this important parasite. Defining the evolutionary trajectories that drive L. infantum fitness in this new environment are of great public health importance as they will allow unique insight into pathways of host/pathogen co-evolution and their consequences for region-specific changes in disease manifestation. This review summarizes current knowledge on L. infantum genetic and phenotypic diversity in the Americas and its possible role in the unique epidemiology of visceral leishmaniasis (VL) in the New World. We highlight the importance of appreciating adaptive molecular mechanisms in L. infantum to understand the parasites' successful establishment on the continent.

Environmental Conditions May Shape the Patterns of Genomic Variations in Leishmania panamensis

Due to the absence of transcriptional regulation of gene expression in Leishmania parasites, it is now well accepted that several forms of genomic variations modulate the levels of critical proteins through changes in gene dosage. We previously observed many of these variations in our reference laboratory strain of L. panamensis (PSC-1 strain), including chromosomes with an increased somy and the presence of a putative linear minichromosome derived from chromosome 34. Here, we compared the previously described genomic variations with those occurring after exposure of this strain to increasing concentrations of trivalent antimony (Sb^III), as well as those present in two geographically unrelated clinical isolates of L. panamensis. We observed changes in the somy of several chromosomes, amplifications of several chromosomal regions, and copy number variations in gene arrays after exposure to Sb^III. Occurrence of amplifications potentially beneficial for the Sb-resistant phenotype appears to be associated with the loss of other forms of amplification, such as the linear minichromosome. In contrast, we found no evidence of changes in somy or amplification of relatively large chromosomal regions in the clinical isolates. In these isolates, the predominant amplifications appear to be those that generate genes arrays; however, in many cases, the amplified arrays have a notably higher number of copies than those from the untreated and Sb-treated laboratory samples.