The Leishmania ribosome: more than passive mRNA translating machinery.

While the central dogma of molecular biology describes how genetic information flows, gene expression is also affected by epigenetic and epitranscriptomic processes. A recent report by Rajan et al. demonstrates how pseudouridylation of a Leishmania ribosomal rRNA affects the expression of particular proteins: an example of epitranslatomic control.

A putative Leishmania DNA polymerase theta protects the parasite against oxidative damage.

Leishmania infantum is a protozoan parasite that is phagocytized by human macrophages. The host macrophages kill the parasite by generating oxidative compounds that induce DNA damage. We have identified, purified and biochemically characterized a DNA polymerase θ from L. infantum (LiPolθ), demonstrating that it is a DNA-dependent DNA polymerase involved in translesion synthesis of 8oxoG, abasic sites and thymine glycol lesions. Stably transfected L. infantum parasites expressing LiPolθ were significantly more resistant to oxidative and interstrand cross-linking agents, e.g. hydrogen peroxide, cisplatin and mitomycin C. Moreover, LiPolθ-overexpressing parasites showed an increased infectivity toward its natural macrophage host. Therefore, we propose that LiPolθ is a translesion synthesis polymerase involved in parasite DNA damage tolerance, to confer resistance against macrophage aggression.

Complete and de novo assembly of the Leishmania braziliensis (M2904) genome.

Leishmania braziliensis is the etiological agent of American mucosal leishmaniasis, one of the most severe clinical forms of leishmaniasis. Here, we report the assembly of the L. braziliensis (M2904) genome into 35 continuous chromosomes. Also, the annotation of 8395 genes is provided. The public availability of this information will contribute to a better knowledge of this pathogen and help in the search for vaccines and novel drug targets aimed to control the disease caused by this Leishmania species.

Datasets of Iso-Seq transcripts for decoding transcriptome complexity in four Leishmania species.

The Iso-Seq technology, based on PacBio sequencing, enables the generation of high-quality, full-length transcripts, providing insights into transcriptome complexity. In this study, total RNA from promastigotes of four Leishmania species (Leishmania braziliensis, Leishmania donovani, Leishmania infantum and Leishmania major) was sequenced using Single Molecule, Real-Time (SMRT) Sequencing (PacBio) methodology. The Iso-seq transcripts were categorized as either complete or truncated according to the presence or absence of the Spliced-Leader (SL) sequence at their 5'-end, respectively. Moreover, only transcripts having a poly-A+ at their 3'-end were considered. Supplied datasets represent valuable information that may help to uncover novel transcripts and alternative splicing events in a parasite that regulates its gene expression at the post-transcriptional level. A better knowledge of gene expression regulation in Leishmania will open avenues for the development of new drugs to treat leishmaniasis, a devastating disease that has worldwide distribution. Additionally, the bioinformatics pipeline followed here may guide the analysis of Iso-Seq data derived from related trypanosomatids like Trypanosoma cruzi (Chagas disease agent) and Trypanosoma brucei (sleeping disease). © 2023 The Authors. Published by Elsevier Inc. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).

A Proteogenomic Approach to Unravel New Proteins Encoded in the Leishmania donovani (HU3) Genome.

The high-throughput proteomics data generated by increasingly more sensible mass spectrometers greatly contribute to our better understanding of molecular and cellular mechanisms operating in live beings. Nevertheless, proteomics analyses are based on accurate genomic and protein annotations, and some information may be lost if these resources are incomplete. Here, we show that most proteomics data may be recovered by interconnecting genomics and proteomics approaches (i.e., following a proteogenomic strategy), resulting, in turn, in an improvement of gene/protein models. In this study, we generated proteomics data from Leishmania donovani (HU3 strain) promastigotes that allowed us to detect 1908 proteins in this developmental stage on the basis of the currently annotated proteins available in public databases. However, when the proteomics data were searched against all possible open reading frames existing in the L. donovani genome, twenty new protein-coding genes could be annotated. Additionally, 43 previously annotated proteins were extended at their N-terminal ends to accommodate peptides detected in the proteomics data. Also, different post-translational modifications (phosphorylation, acetylation, methylation, among others) were found to occur in a large number of Leishmania proteins. Finally, a detailed comparative analysis of the L. donovani and Leishmania major experimental proteomes served to illustrate how inaccurate conclusions can be raised if proteomes are compared solely on the basis of the listed proteins identified in each proteome. Finally, we have created data entries (based on freely available repositories) to provide and maintain updated gene/protein models. Raw data are available via ProteomeXchange with the identifier PXD051920.

Exploring the functionality and conservation of Alba proteins in Trypanosoma cruzi: A focus on biological diversity and RNA binding ability.

Trypanosoma cruzi is the causative agent of Chagas disease, as well as a trypanosomatid parasite with a complex biological cycle that requires precise mechanisms for regulating gene expression. In Trypanosomatidae, gene regulation occurs mainly at the mRNA level through the recognition of cis elements by RNA-binding proteins (RBPs). Alba family members are ubiquitous DNA/RNA-binding proteins with representatives in trypanosomatid parasites functionally related to gene expression regulation. Although T. cruzi possesses two groups of Alba proteins (Alba1/2 and Alba30/40), their functional role remains poorly understood. Thus, herein, a characterization of T. cruzi Alba (TcAlba) proteins was undertaken. Physicochemical, structural, and phylogenetic analysis of TcAlba showed features compatible with RBPs, such as hydrophilicity, RBP domains/motifs, and evolutionary conservation of the Alba-domain, mainly regarding other trypanosomatid Alba. However, in silico RNA interaction analysis of T. cruzi Alba proteins showed that TcAlba30/40 proteins, but not TcAlba1/2, would directly interact with the assayed RNA molecules, suggesting that these two groups of TcAlba proteins have different targets. Given the marked differences existing between both T. cruzi Alba groups (TcAlba1/2 and TcAlba30/40), regarding sequence divergence, RNA binding potential, and life-cycle expression patterns, we suggest that they would be involved in different biological processes.

Alpha tubulin genes from Leishmania braziliensis: genomic organization, gene structure and insights on their expression.

BACKGROUND: Alpha tubulin is a fundamental component of the cytoskeleton which is responsible for cell shape and is involved in cell division, ciliary and flagellar motility and intracellular transport. Alpha tubulin gene expression varies according to the morphological changes suffered by Leishmania in its life cycle. However, the objective of studying the mechanisms responsible for the differential expression has resulted to be a difficult task due to the complex genome organization of tubulin genes and to the non-conventional mechanisms of gene regulation operating in Leishmania. RESULTS: We started this work by analyzing the genomic organization of α-tubulin genes in the Leishmania braziliensis genome database. The genomic organization of L. braziliensis α-tubulin genes differs from that existing in the L. major and L. infantum genomes. Two loci containing α-tubulin genes were found in the chromosomes 13 and 29, even though the existence of sequence gaps does not allow knowing the exact number of genes at each locus. Southern blot assays showed that α-tubulin locus at chromosome 13 contains at least 8 gene copies, which are tandemly organized with a 2.08-kb repetition unit; the locus at chromosome 29 seems to contain a sole α-tubulin gene. In addition, it was found that L. braziliensis α-tubulin locus at chromosome 13 contains two types of α-tubulin genes differing in their 3' UTR, each one presumably containing different regulatory motifs. It was also determined that the mRNA expression levels of these genes are controlled by post-transcriptional mechanisms tightly linked to the growth temperature. Moreover, the decrease in the α-tubulin mRNA abundance observed when promastigotes were cultured at 35°C was accompanied by parasite morphology alterations, similar to that occurring during the promastigote to amastigote differentiation. CONCLUSIONS: Information found in the genome databases indicates that α-tubulin genes have been reorganized in a drastic manner along Leishmania speciation. In the L. braziliensis genome database, two loci containing α-tubulin sequences were found, but only the locus at chromosome 13 contains the prototypic α-tubulin genes, which are repeated in a head-to-tail manner. Also, we determined that the levels of α-tubulin mRNAs are down-regulated drastically in response to heat shock by a post-transcriptional mechanism which is dependent upon active protein synthesis.

The transcriptome of Leishmania major in the axenic promastigote stage: transcript annotation and relative expression levels by RNA-seq.

BACKGROUND: Although the genome sequence of the protozoan parasite Leishmania major was determined several years ago, the knowledge of its transcriptome was incomplete, both regarding the real number of genes and their primary structure. RESULTS: Here, we describe the first comprehensive transcriptome analysis of a parasite from the genus Leishmania. Using high-throughput RNA sequencing (RNA-seq), a total of 10285 transcripts were identified, of which 1884 were considered novel, as they did not match previously annotated genes. In addition, our data indicate that current annotations should be modified for many of the genes. The detailed analysis of the transcript processing sites revealed extensive heterogeneity in the spliced leader (SL) and polyadenylation addition sites. As a result, around 50% of the genes presented multiple transcripts differing in the length of the UTRs, sometimes in the order of hundreds of nucleotides. This transcript heterogeneity could provide an additional source for regulation as the different sizes of UTRs could modify RNA stability and/or influence the efficiency of RNA translation. In addition, for the first time for the Leishmania major promastigote stage, we are providing relative expression transcript levels. CONCLUSIONS: This study provides a concise view of the global transcriptome of the L. major promastigote stage, providing the basis for future comparative analysis with other development stages or other Leishmania species.

Refinement of Leishmania donovani Genome Annotations in the Light of Ribosome-Protected mRNAs Fragments (Ribo-Seq Data).

Advances in next-generation sequencing methodologies have facilitated the assembly of an ever-increasing number of genomes. Gene annotations are typically conducted via specialized software, but the most accurate results require additional manual curation that incorporates insights derived from functional and bioinformatic analyses (e.g., transcriptomics, proteomics, and phylogenetics). In this study, we improved the annotation of the Leishmania donovani (strain HU3) genome using publicly available data from the deep sequencing of ribosome-protected mRNA fragments (Ribo-Seq). As a result of this analysis, we uncovered 70 previously non-annotated protein-coding genes and improved the annotation of around 600 genes. Additionally, we present evidence for small upstream open reading frames (uORFs) in a significant number of transcripts, indicating their potential role in the translational regulation of gene expression. The bioinformatics pipelines developed for these analyses can be used to improve the genome annotations of other organisms for which Ribo-Seq data are available. The improvements provided by these studies will bring us closer to the ultimate goal of a complete and accurately annotated L. donovani genome and will enhance future transcriptomics, proteomics, and genetics studies.

Leishmania infantum (JPCM5) Transcriptome, Gene Models and Resources for an Active Curation of Gene Annotations.

Leishmania infantum is one of the causative agents of visceral leishmaniases, the most severe form of leishmaniasis. An improved assembly for the L. infantum genome was published five years ago, yet delineation of its transcriptome remained to be accomplished. In this work, the transcriptome annotation was attained by a combination of both short and long RNA-seq reads. The good agreement between the results derived from both methodologies confirmed that transcript assembly based on Illumina RNA-seq and further delimitation according to the positions of spliced leader (SAS) and poly-A (PAS) addition sites is an adequate strategy to annotate the transcriptomes of Leishmania, a procedure previously used for transcriptome annotation in other Leishmania species and related trypanosomatids. These analyses also confirmed that the Leishmania transcripts boundaries are relatively slippery, showing extensive heterogeneity at the 5'- and 3'-ends. However, the use of RNA-seq reads derived from the PacBio technology (referred to as Iso-Seq) allowed the authors to uncover some complex transcription patterns occurring at particular loci that would be unnoticed by the use of short RNA-seq reads alone. Thus, Iso-Seq analysis provided evidence that transcript processing at particular loci would be more dynamic than expected. Another noticeable finding was the observation of a case of allelic heterozygosity based on the existence of chimeric Iso-Seq reads that might be generated by an event of intrachromosomal recombination. In addition, we are providing the L. infantum gene models, including both UTRs and CDS regions, that would be helpful for undertaking whole-genome expression studies. Moreover, we have built the foundations of a communal database for the active curation of both gene/transcript models and functional annotations for genes and proteins.

The Astonishing Large Family of HSP40/DnaJ Proteins Existing in Leishmania.

Abrupt environmental changes are faced by Leishmania parasites during transmission from a poikilothermic insect vector to a warm-blooded host. Adaptation to harsh environmental conditions, such as nutrient deprivation, hypoxia, oxidative stress and heat shock needs to be accomplished by rapid reconfiguration of gene expression and remodeling of protein interaction networks. Chaperones play a central role in the maintenance of cellular homeostasis, and they are responsible for crucial tasks such as correct folding of nascent proteins, protein translocation across different subcellular compartments, avoiding protein aggregates and elimination of damaged proteins. Nearly one percent of the gene content in the Leishmania genome corresponds to members of the HSP40 family, a group of proteins that assist HSP70s in a variety of cellular functions. Despite their expected relevance in the parasite biology and infectivity, little is known about their functions or partnership with the different Leishmania HSP70s. Here, we summarize the structural features of the 72 HSP40 proteins encoded in the Leishmania infantum genome and their classification into four categories. A review of proteomic data, together with orthology analyses, allow us to postulate cellular locations and possible functional roles for some of them. A detailed study of the members of this family would provide valuable information and opportunities for drug discovery and improvement of current treatments against leishmaniasis.

Assembly of a Large Collection of Maxicircle Sequences and Their Usefulness for Leishmania Taxonomy and Strain Typing.

Parasites of medical importance, such as Leishmania and Trypanosoma, are characterized by the presence of thousands of circular DNA molecules forming a structure known as kinetoplast, within the mitochondria. The maxicircles, which are equivalent to the mitochondrial genome in other eukaryotes, have been proposed as a promising phylogenetic marker. Using whole-DNA sequencing data, it is also possible to assemble maxicircle sequences as shown here and in previous works. In this study, based on data available in public databases and using a bioinformatics workflow previously reported by our group, we assembled the complete coding region of the maxicircles for 26 prototypical strains of trypanosomatid species. Phylogenetic analysis based on this dataset resulted in a robust tree showing an accurate taxonomy of kinetoplastids, which was also able to discern between closely related Leishmania species that are usually difficult to discriminate by classical methodologies. In addition, we provide a dataset of the maxicircle sequences of 60 Leishmania infantum field isolates from America, Western Europe, North Africa, and Eastern Europe. In agreement with previous studies, our data indicate that L. infantum parasites from Brazil are highly homogeneous and closely related to European strains, which were transferred there during the discovery of America. However, this study showed the existence of different L. infantum populations/clades within the Mediterranean region. A maxicircle signature for each clade has been established. Interestingly, two L. infantum clades were found coexisting in the same region of Spain, one similar to the American strains, represented by the Spanish JPCM5 reference strain, and the other, named "non-JPC like", may be related to an important leishmaniasis outbreak that occurred in Madrid a few years ago. In conclusion, the maxicircle sequence emerges as a robust molecular marker for phylogenetic analysis and species typing within the kinetoplastids, which also has the potential to discriminate intraspecific variability.

Live attenuated vaccines, a favorable strategy to provide long-term immunity against protozoan diseases.

The control of diseases caused by protozoan parasites is one of the United Nations' Sustainable Development Goals. In recent years much research effort has gone into developing a new generation of live attenuated vaccines (LAVs) against malaria, Chagas disease and leishmaniasis. However, there is a bottleneck related to their biosafety, production, and distribution that slows downs further development. The success of irradiated or genetically attenuated sporozoites against malaria, added to the first LAV against leishmaniasis to be evaluated in clinical trials, is indicative that the drawbacks of LAVs are gradually being overcome. However, whether persistence of LAVs is a prerequisite for sustained long-term immunity remains to be clarified, and the procedures necessary for clinical evaluation of vaccine candidates need to be standardized.

Evidence of RNA editing in Leishmania braziliensis promastigotes.

RNA editing in trypanosomatids is an elaborate form of post-transcriptional processing that inserts and deletes uridines in many mitochondrial pre-mRNAs, providing the genetic information needed to create functional transcripts. The process has been extensively analyzed in Trypanosoma brucei, Crithidia fasciculata, and Leishmania tarentolae. However, few data exist on this mechanism in pathogenic Leishmania species. Here, we show evidence that this process also operates in Leishmania braziliensis, being the first time that RNA editing has been described in a species of the Viannia subgenus. A partially edited transcript corresponding to the NADH dehydrogenase subunit 8 (ND8) gene was identified in L. braziliensis promastigotes. Sequence analysis allowed the identification of the maxicircle-encoded cryptogene, which shows a high degree of sequence conservation with the corresponding cryptogenes in other Leishmania species. Although an edition pattern could be postulated for the ND8 transcripts in L. braziliensis, attempts to isolate completely edited transcripts by RT-PCR were not fruitful; instead, many transcripts with partial and unexpected editing patterns were isolated. This data, together with our inability to detect full-size transcripts by Northern blotting in promastigotes of L. braziliensis, led us to the suggestion that the strain used in this study (M2904) lacks of critical RNA guides for a complete edition of ND8 transcripts.

Identification of the interactomes associated with SCD6 and RBP42 proteins in Leishmania braziliensis.

Leishmania are protozoan parasites responsible for leishmaniasis. These parasites present a precise gene regulation that allows them to survive different environmental conditions during their digenetic life cycle. This adaptation depends on the regulation of the expression of a wide variety of genes, which occurs, mainly at the post-transcriptional level. This differential gene expression is achieved by mechanisms based mainly in RNA binding proteins that regulate the translation and/or stability of mRNA targets by interaction with cis elements principally located in the untranslated regions (UTR). In recent studies, our group identified and characterized two proteins, SCD6 and RBP42, as RNA binding proteins in Leishmania braziliensis. To find clues about the cellular processes in which these proteins are involved, this work was aimed to determine the SCD6- and RBP42-interacting proteins (interactome) in L. braziliensis promastigotes. For this purpose, after an in vivo UV cross-linking, cellular extracts were used to immunoprecipitated, by specific antibodies, protein complexes in which SCD6 or RBP42 were present. Protein mass spectrometry analysis of the immunoprecipitated proteins identified 96 proteins presumably associated with SCD6 and 173 proteins associated with RBP42. Notably, a significant proportion of the identified proteins were shared in both interactomes, indicating a possible functional relationship between SCD6 and RBP42. Remarkably, many of the proteins identified in the SCD6 and RBP42 interactomes are related to RNA metabolism and translation processes, and many of them have been described as components of ribonucleoprotein (RNP) granules in Leishmania and related trypanosomatids. Thus, these results support a role of SCD6 and RBP42 in the assembly and/or function of mRNA-protein complexes, participating in the fate (decay/accumulation/translation) of L. braziliensis transcripts. SIGNIFICANCE: Parasites of the Leishmania genus present a particular regulation of gene expression, operating mainly at the post-transcriptional level, surely aimed to modulate quickly both mRNA and protein levels to survive the sudden environmental changes that occur during a parasite's life cycle as it moves from one host to another. This regulation of gene expression processes would be governed by the interaction of mRNA with RNA binding proteins. Nevertheless, the entirety of protein networks involved in these regulatory processes is far from being understood. In this regard, our work is contributing to stablish protein networks in which the L. braziliensis SCD6 and RBP42 proteins are involved; these proteins, in previous works, have been described as RNA binding proteins and found to participate in gene regulation in different cells and organisms. Additionally, our data point out a possible functional relationship between SCD6 and RBP42 proteins as constituents of mRNA granules, like processing bodies or stress granules, which are essential structures in the regulation of gene expression. This knowledge could provide a new approach for the development of therapeutic targets to control Leishmania infections.

Gene Annotation and Transcriptome Delineation on a De Novo Genome Assembly for the Reference Leishmania major Friedlin Strain.

Leishmania major is the main causative agent of cutaneous leishmaniasis in humans. The Friedlin strain of this species (LmjF) was chosen when a multi-laboratory consortium undertook the objective of deciphering the first genome sequence for a parasite of the genus Leishmania. The objective was successfully attained in 2005, and this represented a milestone for Leishmania molecular biology studies around the world. Although the LmjF genome sequence was done following a shotgun strategy and using classical Sanger sequencing, the results were excellent, and this genome assembly served as the reference for subsequent genome assemblies in other Leishmania species. Here, we present a new assembly for the genome of this strain (named LMJFC for clarity), generated by the combination of two high throughput sequencing platforms, Illumina short-read sequencing and PacBio Single Molecular Real-Time (SMRT) sequencing, which provides long-read sequences. Apart from resolving uncertain nucleotide positions, several genomic regions were reorganized and a more precise composition of tandemly repeated gene loci was attained. Additionally, the genome annotation was improved by adding 542 genes and more accurate coding-sequences defined for around two hundred genes, based on the transcriptome delimitation also carried out in this work. As a result, we are providing gene models (including untranslated regions and introns) for 11,238 genes. Genomic information ultimately determines the biology of every organism; therefore, our understanding of molecular mechanisms will depend on the availability of precise genome sequences and accurate gene annotations. In this regard, this work is providing an improved genome sequence and updated transcriptome annotations for the reference L. major Friedlin strain.

Protective Efficacy in a Hamster Model of a Multivalent Vaccine for Human Visceral Leishmaniasis (MuLeVaClin) Consisting of the KMP11, LEISH-F3+, and LJL143 Antigens in Virosomes, Plus GLA-SE Adjuvant.

Visceral leishmaniasis (VL) is the most severe clinical form of leishmaniasis, fatal if untreated. Vaccination is the most cost-effective approach to disease control; however, to date, no vaccines against human VL have been made available. This work examines the efficacy of a novel vaccine consisting of the Leishmania membrane protein KMP11, LEISH-F3+ (a recombinant fusion protein, composed of epitopes of the parasite proteins nucleoside hydrolase, sterol-24-c-methyltransferase, and cysteine protease B), and the sand fly salivary protein LJL143, in two dose ratios. The inclusion of the TLR4 agonist GLA-SE as an adjuvant, and the use of virosomes (VS) as a delivery system, are also examined. In a hamster model of VL, the vaccine elicited antigen-specific immune responses prior to infection with Leishmania infantum. Of note, the responses were greater when higher doses of KMP11 and LEISH-F3+ proteins were administered along with the GLA-SE adjuvant and/or when delivered within VS. Remarkably, hamsters immunized with the complete combination (i.e., all antigens in VS + GLA-SE) showed significantly lower parasite burdens in the spleen compared to those in control animals. This protection was underpinned by a more intense, specific humoral response against the KMP11, LEISH-F3+, and LJL143 antigens in vaccinated animals, but a significantly less intense antibody response to the pool of soluble Leishmania antigens (SLA). Overall, these results indicate that this innovative vaccine formulation confers protection against L. infantum infection, supporting the advancement of the vaccine formulation into process development and manufacturing and the conduction of toxicity studies towards future phase I human clinical trials.

The Experimental Proteome of Leishmania infantum Promastigote and Its Usefulness for Improving Gene Annotations.

Leishmania infantum causes visceral leishmaniasis (kala-azar), the most severe form of leishmaniasis, which is lethal if untreated. A few years ago, the re-sequencing and de novo assembling of the L. infantum (JPCM5 strain) genome was accomplished, and now we aimed to describe and characterize the experimental proteome of this species. In this work, we performed a proteomic analysis from axenic cultured promastigotes and carried out a detailed comparison with other Leishmania experimental proteomes published to date. We identified 2352 proteins based on a search of mass spectrometry data against a database built from the six-frame translated genome sequence of L. infantum. We detected many proteins belonging to organelles such as glycosomes, mitochondria, or flagellum, as well as many metabolic enzymes and many putative RNA binding proteins and molecular chaperones. Moreover, we listed some proteins presenting post-translational modifications, such as phosphorylations, acetylations, and methylations. On the other hand, the identification of peptides mapping to genomic regions previously annotated as non-coding allowed for the correction of annotations, leading to the N-terminal extension of protein sequences and the uncovering of eight novel protein-coding genes. The alliance of proteomics, genomics, and transcriptomics has resulted in a powerful combination for improving the annotation of the L. infantum reference genome.

A postgenomic view of the heat shock proteins in kinetoplastids.

The kinetoplastids Leishmania major, Trypanosoma brucei and Trypanosoma cruzi are causative agents of a diverse spectrum of human diseases: leishmaniasis, sleeping sickness and Chagas' disease, respectively. These protozoa possess digenetic life cycles that involve development in mammalian and insect hosts. It is generally accepted that temperature is a triggering factor of the developmental programme allowing the adaptation of the parasite to the mammalian conditions. The heat shock response is a general homeostatic mechanism that protects cells from the deleterious effects of environmental stresses, such as heat. This response is universal and includes the synthesis of the heat-shock proteins (HSPs). In this review, we summarize the salient features of the different HSP families and describe their main cellular functions. In parallel, we analyse the composition of these families in kinetoplastids according to literature data and our understanding of genome sequence data. The genome sequences of these parasites have been recently completed. The HSP families described here are: HSP110, HSP104, group I chaperonins, HSP90, HSP70, HSP40 and small HSPs. All these families are widely represented in these parasites. In particular, kinetoplastids possess an unprecedented number of members of the HSP70, HSP60 and HSP40 families, suggesting key roles for these HSPs in their biology.

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

Besides their medical relevance, Leishmania is an adequate model for studying post-transcriptional mechanisms of gene expression. In this microorganism, mRNA degradation/stabilization mechanisms together with translational control and post-translational modifications of proteins are the major drivers of gene expression. Leishmania parasites develop as promastigotes in sandflies and as amastigotes in mammalians, and during host transmission, the parasite experiences a sudden temperature increase. Here, changes in the transcriptome of Leishmania major promastigotes after a moderate heat shock were analysed by RNA-seq. Several of the up-regulated transcripts code for heat shock proteins, other for proteins previously reported to be amastigote-specific and many for hypothetical proteins. Many of the transcripts experiencing a decrease in their steady-state levels code for transporters, proteins involved in RNA metabolism or translational factors. In addition, putative long noncoding RNAs were identified among the differentially expressed transcripts. Finally, temperature-dependent changes in the selection of the spliced leader addition sites were inferred from the RNA-seq data, and particular cases were further validated by RT-PCR and Northern blotting. This study provides new insights into the post-transcriptional mechanisms by which Leishmania modulate gene expression.