Leishmania: Responding to environmental signals and challenges without regulated transcription

Bromodomain Factor 5 as a Target for Antileishmanial Drug Discovery

Leishmaniases are a collection of neglected tropical diseases caused by kinetoplastid parasites in the genus Leishmania. Current chemotherapies are severely limited, and the need for new antileishmanials is of pressing international importance. Bromodomains are epigenetic reader domains that have shown promising therapeutic potential for cancer therapy and may also present an attractive target to treat parasitic diseases. Here, we investigate Leishmania donovani bromodomain factor 5 (LdBDF5) as a target for antileishmanial drug discovery. LdBDF5 contains a pair of bromodomains (BD5.1 and BD5.2) in an N-terminal tandem repeat. We purified recombinant bromodomains of L. donovani BDF5 and determined the structure of BD5.2 by X-ray crystallography. Using a histone peptide microarray and fluorescence polarization assay, we identified binding interactions of LdBDF5 bromodomains with acetylated peptides derived from histones H2B and H4. In orthogonal biophysical assays including thermal shift assays, fluorescence polarization, and NMR, we showed that BDF5 bromodomains bind to human bromodomain inhibitors SGC-CBP30, bromosporine, and I-BRD9; moreover, SGC-CBP30 exhibited activity against Leishmania promastigotes in cell viability assays. These findings exemplify the potential BDF5 holds as a possible drug target in Leishmania and provide a foundation for the future development of optimized antileishmanial compounds targeting this epigenetic reader protein.

The paradigm of intracellular parasite survival and drug resistance in leishmanial parasite through genome plasticity and epigenetics: Perception and future perspective

Leishmania is an intracellular, zoonotic, kinetoplastid eukaryote with more than 1.2 million cases all over the world. The leishmanial chromosomes are divided into polymorphic chromosomal ends, conserved central domains, and antigen-encoding genes found in telomere-proximal regions. The genome flexibility of chromosomal ends of the leishmanial parasite is known to cause drug resistance and intracellular survival through the evasion of host defense mechanisms. Therefore, in this review, we discuss the plasticity of Leishmania genome organization which is the primary cause of drug resistance and parasite survival. Moreover, we have not only elucidated the causes of such genome plasticity which includes aneuploidy, epigenetic factors, copy number variation (CNV), and post-translation modification (PTM) but also highlighted their impact on drug resistance and parasite survival.

Leishmania parasites exchange drug-resistance genes through extracellular vesicles

Leishmania are eukaryotic parasites that have retained the ability to produce extracellular vesicles (EVs) through evolution. To date, it has been unclear if different DNA entities could be associated with Leishmania EVs and whether these could constitute a mechanism of horizontal gene transfer (HGT). Herein, we investigate the DNA content of EVs derived from drug-resistant parasites, as well as the EVs' potential to act as shuttles for DNA transfer. Next-generation sequencing and PCR assays confirm the enrichment of amplicons carrying drug-resistance genes associated with EVs. Transfer assays of drug-resistant EVs highlight a significant impact on the phenotype of recipient parasites induced by the expression of the transferred DNA. Recipient parasites display an enhanced growth and better control of oxidative stress. We provide evidence that eukaryotic EVs function as efficient mediators in HGT, thereby facilitating the transmission of drug-resistance genes and increasing the fitness of cells when encountering stressful environments.

Expression Profile of Genes Related to the Th17 Pathway in Macrophages Infected by Leishmania major and Leishmania amazonensis: The Use of Gene Regulatory Networks in Modeling This Pathway

Leishmania amazonensis and Leishmania major are the causative agents of cutaneous and mucocutaneous diseases. The infections‘ outcome depends on host–parasite interactions and Th1/Th2 response, and in cutaneous form, regulation of Th17 cytokines has been reported to maintain inflammation in lesions. Despite that, the Th17 regulatory scenario remains unclear. With the aim to gain a better understanding of the transcription factors (TFs) and genes involved in Th17 induction, in this study, the role of inducing factors of the Th17 pathway in Leishmania–macrophage infection was addressed through computational modeling of gene regulatory networks (GRNs). The Th17 GRN modeling integrated experimentally validated data available in the literature and gene expression data from a time-series RNA-seq experiment (4, 24, 48, and 72 h post-infection). The generated model comprises a total of 10 TFs, 22 coding genes, and 16 cytokines related to the Th17 immune modulation. Addressing the Th17 induction in infected and uninfected macrophages, an increase of 2- to 3-fold in 4–24 h was observed in the former. However, there was a decrease in basal levels at 48–72 h for both groups. In order to evaluate the possible outcomes triggered by GRN component modulation in the Th17 pathway. The generated GRN models promoted an integrative and dynamic view of Leishmania–macrophage interaction over time that extends beyond the analysis of single-gene expression.

G Protein-Coupled Receptors as Potential Intercellular Communication Mediators in Trypanosomatidae

Detection and transduction of environmental signals, constitute a prerequisite for successful parasite invasion; i.e., Leishmania transmission, survival, pathogenesis and disease manifestation and dissemination, with diverse molecules functioning as inter-cellular signaling ligands. Receptors [i.e., G protein-coupled receptors (GPCRs)] and their associated transduction mechanisms, well conserved through evolution, specialize in this function. However, canonical GPCR-related signal transduction systems have not been described in Leishmania, although orthologs, with reduced domains and function, have been identified in Trypanosomatidae. These inter-cellular communication means seem to be essential for multicellular and unicellular organism’s survival. GPCRs are flexible in their molecular architecture and may interact with the so-called receptor activity-modifying proteins (RAMPs), which modulate their function, changing GPCRs pharmacology, acting as chaperones and regulating signaling and/or trafficking in a receptor-dependent manner. In the skin, vasoactive- and neuro- peptides released in response to the noxious stimuli represented by the insect bite may trigger parasite physiological responses, for example, chemotaxis. For instance, in Leishmania (V.) braziliensis, sensory [Substance P, SP, chemoattractant] and autonomic [Vasoactive Intestinal Peptide, VIP, and Neuropeptide Y, NPY, chemorepellent] neuropeptides at physiological levels stimulate in vitro effects on parasite taxis. VIP and NPY chemotactic effects are impaired by their corresponding receptor antagonists, suggesting that the stimulated responses might be mediated by putative GPCRs (with essential conserved receptor domains); the effect of SP is blocked by [(D-Pro 2, D-Trp7,9]-Substance P (10^-6 M)] suggesting that it might be mediated by neurokinin-1 transmembrane receptors. Additionally, vasoactive molecules like Calcitonin Gene-Related Peptide [CGRP] and Adrenomedullin [AM], exert a chemorepellent effect and increase the expression of a 24 kDa band recognized in western blot analysis by (human-)-RAMP-2 antibodies. In-silico search oriented towards GPCRs-like receptors and signaling cascades detected a RAMP-2-aligned sequence corresponding to Leishmania folylpolyglutamate synthase and a RAMP-3 aligned protein, a hypothetical Leishmania protein with yet unknown function, suggesting that in Leishmania, CGRP and AM activities may be modulated by RAMP- (-2) and (-3) homologs. The possible presence of proteins and molecules potentially involved in GPCRs cascades, i.e., RAMPs, signpost conservation of ancient signaling systems associated with responses, fundamental for cell survival, (i.e., taxis and migration) and may constitute an open field for description of pharmacophores against Leishmania parasites.

Microevolution of Trypanosoma cruzi reveals hybridization and clonal mechanisms driving rapid genome diversification

Protozoa and fungi are known to have extraordinarily diverse mechanisms of genetic exchange. However, the presence and epidemiological relevance of genetic exchange in Trypanosoma cruzi, the agent of Chagas disease, has been controversial and debated for many years. Field studies have identified both predominantly clonal and sexually recombining natural populations. Two of six natural T. cruzi lineages (TcV and TcVI) show hybrid mosaicism, using analysis of single-gene locus markers. The formation of hybrid strains in vitro has been achieved and this provides a framework to study the mechanisms and adaptive significance of genetic exchange. Using whole genome sequencing of a set of experimental hybrids strains, we have confirmed that hybrid formation initially results in tetraploid parasites. The hybrid progeny showed novel mutations that were not attributable to either (diploid) parent showing an increase in amino acid changes. In long-term culture, up to 800 generations, there was a variable but gradual erosion of progeny genomes towards triploidy, yet retention of elevated copy number was observed at several core housekeeping loci. Our findings indicate hybrid formation by fusion of diploid T. cruzi, followed by sporadic genome erosion, but with substantial potential for adaptive evolution, as has been described as a genetic feature of other organisms, such as some fungi.

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.

Genome wide comparison of Leishmania donovani strains from Indian visceral leishmaniasis and para-kala-azar dermal leishmaniasis patients

Visceral leishmaniasis (VL) or Kala-azar, primarily caused by Leishmania donovani, is a major health concern in many countries including India. Growing unresponsiveness among the parasites toward the available drugs is alarming, and so, it is necessary to decipher the underlying mechanism of such development for designing new therapeutics. Moreover, even after successful treatment, some VL patients develop apparently harmless skin lesions known as post-kala-azar dermal leishmaniasis (PKDL) which may serve as a reservoir of the parasite in the transmission cycle. Furthermore, recent reports of para-kala-azar dermal leishmaniasis (para-KDL) cases having PKDL manifestation with concomitant VL, emphasize the necessity of more attention to address complex nature of the parasite for eradicating the disease effectively. In the present study, whole genome sequencing is performed with sodium stibogluconate (SSG) sensitive and resistant L. donovani strains along with SSG sensitive para-KDL strains, derived from the clinical isolates of Indian patients to identify the genomic variations among them. Notably, the analyses of chromosome somy values and genome wide mutation profile in the coding regions reveal distinct clustering of the para-KDL strains with 24 genes being mutated uniquely in this group. Such distinguishing genomic changes among the para-KDL strains could be significant for the parasites to become dermatotropic. Overall, the study reveals a possible correlation of the development of SSG resistance and the transition towards the manifestation of PKDL with chromosome aneuploidy and non-synonymous genetic variations in the coding sequences of the L. donovani strains from Indian patients.

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.

Life Cycle Stage-Specific Accessibility of Leishmania donovani Chromatin at Transcription Start Regions

Leishmania donovani is a parasitic protist that causes the lethal Kala-azar fever in India and East Africa. Gene expression in Leishmania is regulated by gene copy number variation and inducible translation while RNA synthesis initiates at a small number of sites per chromosome and proceeds through polycistronic transcription units, precluding a gene-specific regulation (C. Clayton and M. Shapira, Mol Biochem Parasitol 156:93–101, 2007, https://doi.org/10.1016/j.molbiopara.2007.07.007). Here, we analyze the dynamics of chromatin structure in both life cycle stages of the parasite and find evidence for an additional, epigenetic gene regulation pathway in this early branching eukaryote. The assay for transposase-accessible chromatin using sequencing (ATAC-seq) analysis (J. D. Buenrostro, P. G. Giresi, L. C. Zaba, H. Y. Chang, and W. J. Greenleaf, Nat Methods 10:1213–1218, 2013, https://doi.org/10.1038/nmeth.2688) predominantly shows euchromatin at transcription start regions in fast-growing promastigotes, but mostly heterochromatin in the slowly proliferating amastigotes, the mammalian stage, reflecting a previously shown increase of histone synthesis in the latter stage.

IMPORTANCELeishmania parasites are important pathogens with a global impact and cause poverty-related illness and death. They are devoid of classic cis- and trans-acting transcription regulators but use regulated translation and gene copy number variations to adapt to hosts and environments. In this work, we show that transcription start regions present as open euchromatin in fast-growing insect stages but as less-accessible heterochromatin in the slowly proliferating amastigote stage, indicating an epigenetic control of gene accessibility in this early branching eukaryotic pathogen. This finding should stimulate renewed interest in the control of RNA synthesis in Leishmania and related parasites.

Downregulation of FeSOD-A expression in Leishmania infantum alters trivalent antimony and miltefosine susceptibility

Background

Superoxide dismutase (SOD), a central component of the antioxidant defence system of most organisms, removes excess superoxide anions by converting them to oxygen and hydrogen peroxide. As iron (Fe) SOD is absent in the human host, this enzyme is a promising molecular target for drug development against trypanosomatids.

Results

We obtained Leishmania infantum mutant clones with lower FeSOD-A expression and investigated their phenotypes. Our attempts to delete this enzyme-coding gene using three different methodologies (conventional allelic replacement or two different CRISPR/methods) failed, as FeSOD-A gene copies were probably retained by aneuploidy or gene amplification. Promastigote forms of WT and mutant parasites were used in quantitative reverse-transcription polymerase chain reaction (RT-qPCR) and western blot analyses, and these parasite forms were also used to assess drug susceptibility. RT-qPCR and western blot analyses revealed that FeSOD-A transcript and protein levels were lower in FeSOD-A^−/−/+L. infantum mutant clones than in the wild-type (WT) parasite. The decrease in FeSOD-A expression in L. infantum did not interfere with the parasite growth or susceptibility to amphotericin B. Surprisingly, FeSOD-A^−/−/+L. infantum mutant clones were 1.5- to 2.0-fold more resistant to trivalent antimony and 2.4- to 2.7-fold more resistant to miltefosine. To investigate whether the decrease in FeSOD-A expression was compensated by other enzymes, the transcript levels of five FeSODs and six enzymes from the antioxidant defence system were assessed by RT-qPCR. The transcript level of the enzyme ascorbate peroxidase increased in both the FeSOD-A^−/−/+ mutants tested. The FeSOD-A^−/−/+ mutant parasites were 1.4- to 1.75-fold less tolerant to oxidative stress generated by menadione. Infection analysis using THP-1 macrophages showed that 72 h post-infection, the number of infected macrophages and their intracellular multiplication rate were lower in the FeSOD-A^−/−/+ mutant clones than in the WT parasite.

Conclusions

The unsuccessful attempts to delete FeSOD-A suggest that this gene is essential in L. infantum. This enzyme plays an important role in the defence against oxidative stress and infectivity in THP-1 macrophages. FeSOD-A-deficient L. infantum parasites deregulate their metabolic pathways related to antimony and miltefosine resistance.

Graphic Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-021-04838-8.