Transcriptomic profiling of gene expression and RNA processing during Leishmania major differentiation

Structural and mechanistic insights into the function of Leishmania ribosome lacking a single pseudouridine modification

Leishmania is the causative agent of cutaneous and visceral diseases affecting millions of individuals worldwide. Pseudouridine (Ψ), the most abundant modification on rRNA, changes during the parasite life cycle. Alterations in the level of a specific Ψ in helix 69 (H69) affected ribosome function. To decipher the molecular mechanism of this phenotype, we determine the structure of ribosomes lacking the single Ψ and its parental strain at ∼2.4-3 Å resolution using cryo-EM. Our findings demonstrate the significance of a single Ψ on H69 to its structure and the importance for its interactions with helix 44 and specific tRNAs. Our study suggests that rRNA modification affects translation of mRNAs carrying codon bias due to selective accommodation of tRNAs by the ribosome. Based on the high-resolution structures, we propose a mechanism explaining how the ribosome selects specific tRNAs.

Evaluation of expression variations in virulence-related genes of Leishmania major after several culture passages compared with Phlebotomus papatasi isolated promastigotes

Cutaneous leishmaniasis (CL) is a prevalent infectious disease with considerable morbidity annually. Here, we aimed to investigate the likely variations in gene expression of glycoprotein63 (gp63), heat shock protein 70 (HSP70), histone, arginase, cysteine protease B (CPB), Leishmania homologue of receptors for activated C kinase (LACK), small hydrophilic endoplasmic reticulum-associated protein (SHERP) in metacyclic promastigotes of L. major isolated from Phlebotomus papatasi sand flies and promastigotes excessively cultured in culture medium. The parasites were collected from suspected CL cases in Pasteur Institute of Iran, cultured and inoculated into the female BALB/c mice (2×106 promastigotes). Sand flies were trapped in Qom province, fed with the blood of euthanized infected mice and subsequently dissected in order to isolate the midgut including stomodeal valve. The metacyclic promastigotes were isolated from Ph. papatasi (Pro-Ppap) using peanut agglutinin test (PNA), then continuously cultured in RPMI-1640 medium enriched with fetal bovine serum, penicillin (100 U/ml) and streptomycin (100 mg/ml) to reach stationary phase (Pro-Stat). The gene expression was evaluated in both parasitic stages (Pro-Ppap and Pro-Stat) using qRT-PCR. Out results showed a significant increased gene expression at Pro-Ppap stage for gp63 (P = 0.002), SHERP (P = 0.001) and histone (P = 0.026) genes, in comparison with Pro-Stat stage. Noticeably, significant changes were, also, demonstrated in 10th to 15th passages [gp63 (P = 0.041), arginase (P = 0.016), LACK (P = 0.025)] and in 5th to 20th passage (SHERP) (P = 0.029). In conclusion, the findings of the present study seem to be essential in designing Leishmania studies, in particular regarding host-parasite interaction, immunization and infectivity studies.

Omics Approaches in Drug Development against Leishmaniasis: Current Scenario and Future Prospects

Leishmaniasis is a zoonotic disease transmitted in humans by the bite of Leishmania-infected phlebotomine sandflies. Each year approximately 58,500 cases of leishmaniasis are diagnosed across the globe, with a mortality rate of nearly seven percent. There are over 20 parasitic strains of Leishmania which are known to cause distinct types of leishmaniasis and pose an endemic threat to humans worldwide. Therefore, it is crucial to develop potential medications and vaccines to combat leishmaniasis. However, the task of developing therapeutic solutions is challenging due to Leishmania's digenetic lifecycle. The challenge is further intensified by cases of resistance against the available drugs. Owing to these challenges, the conventional drug development regimen is further limited by target discovery and ligand suitability for the targets. On the other hand, as an added advantage, the emergence of omics-based tools, such as high-end proteomics, transcriptomics and genomics, has hastened the pace of target discovery and target-based drug development. It is now becoming apparent that multi-omics convergence and an inter-connected systems approach is less time-consuming and more cost-effective for any drug-development process. This comprehensive review is an attempt to summarize the current knowledge on the muti-omics approach in drug development against leishmaniasis. In particular, it elaborates the potential target identification from secreted proteins in various stages of Leishmania infection and also illustrates the convergence of transcriptomic and genomic data towards the collective goal of drug discovery. This review also provides an understanding of the potential parasite's drug targets and drug resistance characteristics of the parasite, which can be used in designing effective and specific therapeutics.

Structural analysis of PpSP15 and PsSP9 sand fly salivary proteins designed with a self-cleavable linker as a live vaccine candidate against cutaneous leishmaniasis

Background

Leishmania parasites are deposited in the host through sand fly bites along with sand fly saliva. Therefore, salivary proteins are promising vaccine candidates for controlling leishmaniasis. Herein, two immunogenic salivary proteins, PpSP15 from Phlebotomus papatasi and PsSP9 from Phlebotomus sergenti, were selected as vaccine candidates to be delivered by live Leishmania tarentolae as vector. The stepwise in silico protocol advantaged in this study for multi-protein design in L. tarentolae is then described in detail.

Methods

All possible combinations of two salivary proteins, PpSP15 and PsSP9, with or without T2A peptide were designed at the mRNA and protein levels. Then, the best combination for the vaccine candidate was selected based on mRNA and protein stability along with peptide analysis.

Results

At the mRNA level, the most favored secondary structure was PpSP15-T2A-PsSP9. At the protein level, the refined three-dimensional models of all combinations were structurally valid; however, local quality estimation showed that the PpSp15-T2A-PsSP9 fusion had higher stability for each amino acid position, with low root-mean-square deviation (RMSD), compared with the original proteins. In silico evaluation confirmed the PpSP15-T2A-PsSP9 combination as a good Th1-polarizing candidate in terms of high IFN-γ production and low IL-10/TGF-β ratio in response to three consecutive immunizations. Potential protein expression was then confirmed by Western blotting.

Conclusions

The approach presented herein is among the first studies to have privileged protein homology modeling along with mRNA analysis for logical live vaccine design-coding multi-proteins.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05437-x.

Four layer multi-omics reveals molecular responses to aneuploidy in Leishmania

Aneuploidy causes system-wide disruptions in the stochiometric balances of transcripts, proteins, and metabolites, often resulting in detrimental effects for the organism. The protozoan parasite Leishmania has an unusually high tolerance for aneuploidy, but the molecular and functional consequences for the pathogen remain poorly understood. Here, we addressed this question in vitro and present the first integrated analysis of the genome, transcriptome, proteome, and metabolome of highly aneuploid Leishmania donovani strains. Our analyses unambiguously establish that aneuploidy in Leishmania proportionally impacts the average transcript- and protein abundance levels of affected chromosomes, ultimately correlating with the degree of metabolic differences between closely related aneuploid strains. This proportionality was present in both proliferative and non-proliferative in vitro promastigotes. However, as in other Eukaryotes, we observed attenuation of dosage effects for protein complex subunits and in addition, non-cytoplasmic proteins. Differentially expressed transcripts and proteins between aneuploid Leishmania strains also originated from non-aneuploid chromosomes. At protein level, these were enriched for proteins involved in protein metabolism, such as chaperones and chaperonins, peptidases, and heat-shock proteins. In conclusion, our results further support the view that aneuploidy in Leishmania can be adaptive. Additionally, we believe that the high karyotype diversity in vitro and absence of classical transcriptional regulation make Leishmania an attractive model to study processes of protein homeostasis in the context of aneuploidy and beyond.

Promastigote-to-Amastigote Conversion in Leishmania spp.—A Molecular View

A key factor in the successful infection of a mammalian host by Leishmania parasites is their conversion from extracellular motile promastigotes into intracellular amastigotes. We discuss the physical and chemical triggers that induce this conversion and the accompanying changes at the molecular level crucial for the survival of these intracellular parasites. Special emphasis is given to the reliance of these trypanosomatids on the post-transcriptional regulation of gene expression but also to the role played by protein kinases, chaperone proteins and proteolytic enzymes. Lastly, we offer a model to integrate the transduction of different stress signals for the induction of stage conversion.

New Insights on Heme Uptake in Leishmania spp

The protozoan parasite Leishmania, responsible for leishmaniasis, is one of the few aerobic organisms that cannot synthesize the essential molecule heme. Therefore, it has developed specialized pathways to scavenge it from its host. In recent years, some proteins involved in the import of heme, such as LHR1 and LFLVCRB, have been identified, but relevant aspects regarding the process remain unknown. Here, we characterized the kinetics of the uptake of the heme analogue Zn(II) Mesoporphyrin IX (ZnMP) in Leishmania major promastigotes as a model of a parasite causing cutaneous leishmaniasis with special focus on the force that drives the process. We found that ZnMP uptake is an active, inducible, and pH-dependent process that does not require a plasma membrane proton gradient but requires the presence of the monovalent cations Na⁺ and/or K⁺. In addition, we demonstrated that this parasite can efflux this porphyrin against a concentration gradient. We also found that ZnMP uptake differs among different dermotropic or viscerotropic Leishmania species and does not correlate with LHR1 or LFLVCRB expression levels. Finally, we showed that these transporters have only partially overlapping functions. Altogether, these findings contribute to a deeper understanding of an important process in the biology of this parasite.

Virulence factors of Leishmania parasite: Their paramount importance in unraveling novel vaccine candidates and therapeutic targets

Leishmaniasis is a neglected tropical disease and remains a global concern for healthcare. It is caused by an opportunistic protozoan parasite belonging to the genus Leishmania and affects millions worldwide. This disease is mainly prevalent in tropical and subtropical regions and is associated with a high risk of public morbidity and mortality if left untreated. Transmission of this deadly disease is aggravated by the bite of female sand-fly vectors (Phlebotomus and Lutzomyia). With time, significant advancement in leishmaniasis-related research has been carried out to cope with the disease burden. Still, the Leishmania parasite has also co-evolved with its host and adapted successfully within the host's lethal milieu/environment. Thus, understanding and knowledge of various leishmanial virulence factors responsible for the parasitic infection are essential for exploring drug targets and vaccine candidates. The present review elucidates the importance of virulence factors in pathogenesis and summarizes the major leishmanial virulence molecules contributing to the parasitic infection during host-pathogen interaction. Furthermore, we have also elaborated on the potential contribution of leishmanial virulence proteins in developing vaccine candidates and exploring novel therapeutics against this parasitic disease. We aim to represent a clearer picture of parasite pathogenesis within the human host that can further aid in unraveling new strategies to fight against the deadly infection of leishmaniasis.

Finding Correlations Between mRNA and Protein Levels in Leishmania Development: Is There a Discrepancy?

Multiple genes and proteins have been identified as differentially expressed in the stages of the Leishmania life cycle. The differentiation processes are implicated in specific transcriptional and proteomic adjustments driven by gene expression regulation mechanisms. Leishmania parasites lack gene-specific transcriptional control, and gene expression regulation mostly depends on posttranscriptional mechanisms. Due to the lack of transcriptional regulation, criticism regarding the relevance of transcript quantification as a possible and efficient prediction of protein levels is recurrent in studies that use transcriptomic information. The advent of high-throughput technologies has improved the analysis of genomes, transcriptomes and proteomes for different organisms under several conditions. Nevertheless, defining the correlation between transcriptional and proteomic profiles requires arduous and expensive work and remains a challenge in Leishmania. In this review, we analyze transcriptomic and proteomic data for several Leishmania species in two different stages of the parasite life cycle: metacyclogenesis and amastigogenesis (amastigote differentiation). We found a correlation between mRNA and protein levels of 60.9% and 69.8% for metacyclogenesis and amastigogenesis, respectively; showing that majority mRNA and protein levels increase or decrease concomitantly. Among the analyzed genes that did not present correlation indicate that transcriptomic data should be carefully interpreted as protein expression. We also discuss possible explanations and mechanisms involved for this lack of correlation.

Polycistronic Expression of Multi-Subunit Complexes in the Eukaryotic Environment: A Narrative Review

Protein complexes are involved in many vital biological processes. Therefore, researchers need these protein complexes for biochemical and biophysical studies. Several methods exist for expressing multi-subunit proteins in eukaryotic cells, such as 2A sequences, IRES, or intein. Nevertheless, each of these elements has several disadvantages that limit their usage. In this article, we suggest a new system for expressing multi-subunit proteins, which have several advantages over existing methods meanwhile it, lacks most of their disadvantages. Leishmania is a unicellular eukaryote and member of the Trypanosomatidae family. In the expression system of Leishmania, pre-long RNAs that contain several protein sequences transcribe. Then these long RNAs separate into mature mRNAs in the process named trans splicing. For producing multi-subunit protein, Leishmania transformed with a vector containing the sequences of all subunits. Therefore, those subunits translate and form the complex under eukaryotic cell conditions. The sequence of each protein must separate by the spatial sequence needed for trans splicing. Based on a Leishmania expression pattern, not only is it possible to produce the complexes with the correct structures and post-translational modifications, but also it is possible to overcome previous method problems.

Immune Responses in Leishmaniases: An Overview

Leishmaniasis is a parasitic, widespread, and neglected disease that affects more than 90 countries in the world. More than 20 Leishmania species cause different forms of leishmaniasis that range in severity from cutaneous lesions to systemic infection. The diversity of leishmaniasis forms is due to the species of parasite, vector, environmental and social factors, genetic background, nutritional status, as well as immunocompetence of the host. Here, we discuss the role of the immune system, its molecules, and responses in the establishment, development, and outcome of Leishmaniasis, focusing on innate immune cells and Leishmania major interactions.

Stress conditions promote Leishmania hybridization in vitro marked by expression of the ancestral gamete fusogen HAP2 as revealed by single-cell RNA-seq

Leishmania are protozoan parasites transmitted by the bite of sand fly vectors producing a wide spectrum of diseases in their mammalian hosts. These diverse clinical outcomes are directly associated with parasite strain and species diversity. Although Leishmania reproduction is mainly clonal, a cryptic sexual cycle capable of producing hybrid genotypes has been inferred from population genetic studies and directly demonstrated by laboratory crosses. Experimentally, mating competence has been largely confined to promastigotes developing in the sand fly midgut. The ability to hybridize culture promastigotes in vitro has been limited so far to low-efficiency crosses between two Leishmania tropica strains, L747 and MA37, that mate with high efficiency in flies. Here, we show that exposure of promastigote cultures to DNA damage stress produces a remarkably enhanced efficiency of in vitro hybridization of the L. tropica strains and extends to other species, including Leishmania donovani, Leishmania infantum, and Leishmania braziliensis, a capacity to generate intra- and interspecific hybrids. Whole-genome sequencing and total DNA content analyses indicate that the hybrids are in each case full genome, mostly tetraploid hybrids. Single-cell RNA sequencing of the L747 and MA37 parental lines highlights the transcriptome heterogeneity of culture promastigotes and reveals discrete clusters that emerge post-irradiation in which genes potentially involved in genetic exchange are expressed, including the ancestral gamete fusogen HAP2. By generating reporter constructs for HAP2, we could select for promastigotes that could either hybridize or not in vitro. Overall, this work reveals that there are specific populations involved in Leishmania hybridization associated with a discernible transcriptomic signature, and that stress facilitated in vitro hybridization can be a transformative approach to generate large numbers of hybrid genotypes between diverse species and strains.

Possible Involvement of Hsp90 in the Regulation of Telomere Length and Telomerase Activity During the Leishmania amazonensis Developmental Cycle and Population Proliferation

The Leishmania developmental cycle comprises three main life forms in two hosts, indicating that the parasite is continually challenged due to drastic environmental changes. The disruption of this cycle is critical for discovering new therapies to eradicate leishmaniasis, a neglected disease that affects millions worldwide. Telomeres, the physical ends of chromosomes, maintain genome stability and cell proliferation and are potential antiparasitic drug targets. Therefore, understanding how telomere length is regulated during parasite development is vital. Here, we show that telomeres form clusters spread in the nucleoplasm of the three parasite life forms. We also observed that amastigotes telomeres are shorter than metacyclic and procyclic promastigotes and that in parasites with continuous in vitro passages, telomere length increases over time. These observed differences in telomere length among parasite’s life stages were not due to lack/inhibition of telomerase since enzyme activity was detected in all parasite life stages, although the catalysis was temperature-dependent. These data led us to test if, similar to other eukaryotes, parasite telomere length maintenance could be regulated by Hsp83, the ortholog of Hsp90 in trypanosomatids, and Leishmania (LHsp90). Parasites were then treated with the Hsp90 inhibitor 17AAG. The results showed that 17AAG disturbed parasite growth, induced accumulation into G2/M phases, and telomere shortening in a time-dependent manner. It has also inhibited procyclic promastigote’s telomerase activity. Besides, LHsp90 interacts with the telomerase TERT component as shown by immunoprecipitation, strongly suggesting a new role for LHsp90 as a parasite telomerase component involved in controlling telomere length maintenance and parasite life span.

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.

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.

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.

Sclerosing Sialadenitis Is Associated With Salivary Gland Hypofunction and a Unique Gene Expression Profile in Sjögren's Syndrome

Purpose

To develop a novel method to quantify the amount of fibrosis in the salivary gland and to investigate the relationship between fibrosis and specific symptoms associated with Sjögren’s syndrome (SS) using this method.

Materials and Methods

Paraffin-embedded labial salivary gland (LSG) slides from 20 female SS patients and their clinical and LSG pathology data were obtained from the Sjögren’s International Collaborative Clinical Alliance. Relative interstitial fibrosis area (RIFA) in Masson’s trichrome-stained LSG sections was quantified from digitally scanned slides and used for correlation analysis. Gene expression levels were assessed by microarray analysis. Core promoter accessibility for RIFA-correlated genes was determined using DNase I hypersensitive sites sequencing analysis.

Results

RIFA was significantly correlated with unstimulated whole saliva flow rate in SS patients. Sixteen genes were significantly and positively correlated with RIFA. In a separate analysis, a group of differentially expressed genes was identified by comparing severe and moderate fibrosis groups. This combined set of genes was distinct from differentially expressed genes identified in lung epithelium from idiopathic pulmonary fibrosis patients compared with controls. Single-cell RNA sequencing analysis of salivary glands suggested most of the RIFA-correlated genes are expressed by fibroblasts in the gland and are in a permissive chromatin state.

Conclusion

RIFA quantification is a novel method for assessing interstitial fibrosis and the impact of fibrosis on SS symptoms. Loss of gland function may be associated with salivary gland fibrosis, which is likely to be driven by a unique set of genes that are mainly expressed by fibroblasts.

Involvement of Leishmania Phosphatases in Parasite Biology and Pathogeny

In the Leishmania lifecycle, the motile promastigote form is transmitted from the sand fly vector to a mammalian host during a blood meal. Inside vertebrate host macrophages, the parasites can differentiate into the amastigote form and multiply, causing leishmaniasis, one of the most significant neglected tropical diseases. Leishmania parasites face different conditions throughout their development inside sand flies. Once in the mammalian host, the parasites have to overcome the microbicide repertoire of the cells of the immune system to successfully establish the infection. In this context, the expression of protein phosphatases is of particular interest. Several members of the serine/threonine-specific protein phosphatase (STP), protein tyrosine phosphatase (PTP), and histidine acid phosphatase (HAcP) families have been described in different Leishmania species. Although their physiological roles have not been fully elucidated, many studies suggest they have an involvement with parasite biology and pathogeny. Phosphatases play a role in adaptation to nutrient starvation during parasite passage through the sand fly midgut. They are also important to parasite virulence, mainly due to the modulation of host cytokine production and impairment of the microbiocidal potential of macrophages. Furthermore, recent whole-genome expression analyses have shown that different phosphatases are upregulated in metacyclic promastigotes, the infective form of the mammalian host. Leishmania phosphatases are also upregulated in drug-resistant strains, probably due to the increase in drug efflux related to the activation of ABC transporters. Throughout this review, we will describe the physiological roles that have been attributed to Leishmania endogenous phosphatases, including their involvement in the adaptation, survival, and proliferation of the parasites inside their hosts.

The Leishmania donovani LDBPK_220120.1 Gene Encodes for an Atypical Dual Specificity Lipid-Like Phosphatase Expressed in Promastigotes and Amastigotes; Substrate Specificity, Intracellular Localizations, and Putative Role(s)

The intracellular protozoan parasites of the Leishmania genus are responsible for Leishmaniases, vector borne diseases with a wide range of clinical manifestations. Leishmania (L.) donovani causes visceral leishmaniasis (kala azar), the most severe of these diseases. Along their biological cycle, Leishmania parasites undergo distinct developmental transitions including metacyclogenesis and differentiation of metacyclic promastigotes (MPs) to amastigotes. Metacyclogenesis inside the phlebotomine sandfly host's midgut converts the procyclic dividing promastigotes to non-dividing infective MPs eventually injected into the skin of mammalian hosts and phagocytosed by macrophages where the MPs are converted inside modified phagolysosomes to the intracellular amastigotes. These developmental transitions involve dramatic changes in cell size and shape and reformatting of the flagellum requiring thus membrane and cytoskeleton remodeling in which phosphoinositide (PI) signaling and metabolism must play central roles. This study reports on the LDBPK_220120.1 gene, the L. donovani ortholog of LmjF.22.0250 from L. major that encodes a phosphatase from the "Atypical Lipid Phosphatases" (ALPs) enzyme family. We confirmed the expression of the LDBPK_220120.1 gene product in both L. donovani promastigotes and axenic amastigotes and showed that it behaves in vitro as a Dual Specificity P-Tyr and monophosphorylated [PI(3)P and PI(4)P] PI phosphatase and therefore named it LdTyrPIP_22 (Leishmaniad onovani Tyrosine PI Phosphatase, gene locus at chromosome 22). By immunofluorescence confocal microscopy we localized the LdTyrPIP_22 in several intracellular sites in the cell body of L. donovani promastigotes and amastigotes and in the flagellum. A temperature and pH shift from 25°C to 37°C and from pH 7 to 5.5, induced a pronounced recruitment of LdTyrPIP_22 epitopes to the flagellar pocket and a redistribution around the nucleus. These results suggest possible role(s) for this P-Tyr/PI phosphatase in the regulation of processes initiated or upregulated by this temperature/pH shift that contribute to the developmental transition from MPs to amastigotes inside the mammalian host macrophages.

Catalase impairs Leishmania mexicana development and virulence

Catalase is one of the most abundant enzymes on Earth. It decomposes hydrogen peroxide, thus protecting cells from dangerous reactive oxygen species. The catalase-encoding gene is conspicuously absent from the genome of most representatives of the family Trypanosomatidae. Here, we expressed this protein from the Leishmania mexicana Β-TUBULIN locus using a novel bicistronic expression system, which relies on the 2A peptide of Teschovirus A. We demonstrated that catalase-expressing parasites are severely compromised in their ability to develop in insects, to be transmitted and to infect mice, and to cause clinical manifestation in their mammalian host. Taken together, our data support the hypothesis that the presence of catalase is not compatible with the dixenous life cycle of Leishmania, resulting in loss of this gene from the genome during the evolution of these parasites.

Time-resolved transcriptional profiling of Trichinella-infected murine myocytes helps to elucidate host-pathogen interactions in the muscle stage

Background

Parasites of the genus Trichinella are the pathogenic agents of trichinellosis, which is a widespread and severe foodborne parasitic disease. Trichinella spiralis resides primarily in mammalian skeletal muscle cells. After invading the cells of the host organism, T. spiralis must elude or invalidate the host’s innate and adaptive immune responses to survive. It is necessary to characterize the pathogenesis of trichinellosis to help to prevent the occurrence and further progression of this disease. The aims of this study were to elucidate the mechanisms of nurse cell formation, pathogenesis and immune evasion of T. spiralis, to provide valuable information for further research investigating the basic cell biology of Trichinella-infected muscle cells and the interaction between T. spiralis and its host.

Methods

We performed transcriptome profiling by RNA sequencing to identify global changes at 1, 3, 7, 10 and 15 days post-infection (dpi) in gene expression in the diaphragm after the parasite entered and persisted within the murine myocytes; the mice were infected by intravenous injection of newborn larvae. Gene expression analysis was based on the alignment results. Differentially expressed genes (DEGs) were identified based on their expression levels in various samples, and functional annotation and enrichment analysis were performed.

Results

The most extensive and dynamic gene expression responses in host diaphragms were observed during early infection (1 dpi). The number of DEGs and genes annotated in the Kyoto Encyclopedia of Genes and Genomes and Gene Ontology databases decreased significantly in the infected mice compared to the uninfected mice at 3 and 7 dpi, suddenly increased sharply at 10 dpi, and then decreased to a lower level at 15 dpi, similar to that observed at 3 and 7 dpi. The massive initial reaction of the murine muscle cells to Trichinella infection steadied in the later stages of infection, with little additional changes detected for the remaining duration of the studied process. Although there were hundreds of DEGs at each time point, only 11 genes were consistently up- or downregulated at all 5 time points.

Conclusions

The gene expression patterns identified in this study can be employed to characterize the coordinated response of T. spiralis-infected myocytes in a time-resolved manner. This comprehensive dataset presents a distinct and sensitive picture of the interaction between host and parasite during intracellular infection, which can help to elucidate how pathogens evade host defenses and coordinate the biological functions of host cells to survive in the mammalian environment.

Epitranscriptome machinery in Trypanosomatids: New players on the table?

Trypanosoma and Leishmania parasites cause devastating tropical diseases resulting in serious global health consequences. These organisms have complex life cycles with mammalian hosts and insect vectors. The parasites must, therefore, survive in different environments, demanding rapid physiological and metabolic changes. These responses depend upon regulation of gene expression, which primarily occurs posttranscriptionally. Altering the composition or conformation of RNA through nucleotide modifications is one posttranscriptional mechanism of regulating RNA fate and function, and modifications including N6-methyladenosine (m6A), N1-methyladenosine (m1A), N5-methylcytidine (m5C), N4-acetylcytidine (ac4C), and pseudouridine (Ψ), dynamically regulate RNA stability and translation in diverse organisms. Little is known about RNA modifications and their machinery in Trypanosomatids, but we hypothesize that they regulate parasite gene expression and are vital for survival. Here, we identified Trypanosomatid homologs for writers of m1A, m5C, ac4C, and Ψ and analyze their evolutionary relationships. We systematically review the evidence for their functions and assess their potential use as therapeutic targets. This work provides new insights into the roles of these proteins in Trypanosomatid parasite biology and treatment of the diseases they cause and illustrates that Trypanosomatids provide an excellent model system to study RNA modifications, their molecular, cellular, and biological consequences, and their regulation and interplay.

Nuclear mRNA maturation and mRNA export control: from trypanosomes to opisthokonts

The passage of mRNAs through the nuclear pores into the cytoplasm is essential in all eukaryotes. For regulation, mRNA export is tightly connected to the full machinery of nuclear mRNA processing, starting at transcription. Export competence of pre-mRNAs gradually increases by both transient and permanent interactions with multiple RNA processing and export factors. mRNA export is best understood in opisthokonts, with limited knowledge in plants and protozoa. Here, I review and compare nuclear mRNA processing and export between opisthokonts and Trypanosoma brucei. The parasite has many unusual features in nuclear mRNA processing, such as polycistronic transcription and trans-splicing. It lacks several nuclear complexes and nuclear-pore-associated proteins that in opisthokonts play major roles in mRNA export. As a consequence, trypanosome mRNA export control is not tight and export can even start co-transcriptionally. Whether trypanosomes regulate mRNA export at all, or whether leakage of immature mRNA to the cytoplasm is kept to a low level by a fast kinetics of mRNA processing remains to be investigated. mRNA export had to be present in the last common ancestor of eukaryotes. Trypanosomes are evolutionary very distant from opisthokonts and a comparison helps understanding the evolution of mRNA export.

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.

Comparative analysis of gene expression between Babesia bovis blood stages and kinetes allowed by improved genome annotation

Throughout their life cycle, Babesia parasites alternate between a mammalian host, where they cause babesiosis, and the tick vector. Transition between hosts results in distinct environmental signals that influence patterns of gene expression, consistent with the morphological and functional changes operating in the parasites during their life stages. In addition, comparing differential patterns of gene expression among mammalian and tick parasite stages can provide clues for developing improved methods of control. Hereby, we upgraded the genome assembly of Babesia bovis, a bovine hemoparasite, closing a 139 kbp gap, and used RNA-Seq datasets derived from mammalian blood and tick kinete stages to update the genome annotation. Of the originally annotated genes, 1,254 required structural changes, and 326 new genes were identified, leading to a different predicted proteome compared to the original annotation. Next, the RNA-Seq data was used to identify B. bovis genes that were differentially expressed in the vertebrate and arthropod hosts. In blood stages, 28% of the genes were upregulated up to 300 fold, whereas 26% of the genes in kinetes, a tick stage, were upregulated up to >19,000 fold. We thus discovered differentially expressed genes that may play key biological roles and serve as suitable targets for the development of vaccines to control bovine babesiosis.

Macrophages as host, effector and immunoregulatory cells in leishmaniasis: Impact of tissue micro-environment and metabolism

Leishmania are protozoan parasites that predominantly reside in myeloid cells within their mammalian hosts. Monocytes and macrophages play a central role in the pathogenesis of all forms of leishmaniasis, including cutaneous and visceral leishmaniasis. The present review will highlight the diverse roles of macrophages in leishmaniasis as initial replicative niche, antimicrobial effectors, immunoregulators and as safe hideaway for parasites persisting after clinical cure. These multiplex activities are either ascribed to defined subpopulations of macrophages (e.g., Ly6C^highCCR2⁺ inflammatory monocytes/monocyte-derived dendritic cells) or result from different activation statuses of tissue macrophages (e.g., macrophages carrying markers of of classical [M1] or alternative activation [M2]). The latter are shaped by immune- and stromal cell-derived cytokines (e.g., IFN-γ, IL-4, IL-10, TGF-β), micro milieu factors (e.g., hypoxia, tonicity, amino acid availability), host cell-derived enzymes, secretory products and metabolites (e.g., heme oxygenase-1, arginase 1, indoleamine 2,3-dioxygenase, NOS2/NO, NOX2/ROS, lipids) as well as by parasite products (e.g., leishmanolysin/gp63, lipophosphoglycan). Exciting avenues of current research address the transcriptional, epigenetic and translational reprogramming of macrophages in a Leishmania species- and tissue context-dependent manner.

Dual transcriptome analysis reveals differential gene expression modulation influenced by Leishmania arginase and host genetic background

The outcome of Leishmania infection is strongly influenced by the host's genetic background. BALB/c mice are susceptible to Leishmania infection, while C57BL/6 mice show discrete resistance. Central to the fate of the infection is the availability of l-arginine and the related metabolic processes in the host and parasite. Depending on l-arginine availability, nitric oxide synthase 2 (NOS2) of the host cell produces nitric oxide (NO) controlling the parasite growth. On the other hand, Leishmania can also use host l-arginine for the production of polyamines through its own arginase activity, thus favouring parasite replication. Considering RNA-seq data, we analysed the dual modulation of host and parasite gene expression of BALB/c or C57BL/6 mouse bone marrow-derived macrophages (BMDMs) after 4 h of infection with Leishmania amazonensis wild-type (La-WT) or L. amazonensis arginase knockout (La-arg-). We identified 12 641 host transcripts and 8282 parasite transcripts by alignment analysis with the respective Mus musculus and L. mexicana genomes. The comparison of BALB/c_La-arg-versus BALB/c_La-WT revealed 233 modulated transcripts, with most related to the immune response and some related to the amino acid transporters and l-arginine metabolism. In contrast, the comparison of C57BL/6_La-arg-vs. C57BL/6_La-WT revealed only 30 modulated transcripts, including some related to the immune response but none related to amino acid transport or l-arginine metabolism. The transcriptome profiles of the intracellular amastigote revealed 94 modulated transcripts in the comparison of La-arg-_BALB/c vs. La-WT_BALB/c and 45 modulated transcripts in the comparison of La-arg-_C57BL/6 vs. La-WT_C57BL/6. Taken together, our data present new insights into the impact of parasite arginase activity on the orchestration of the host gene expression modulation, including in the immune response and amino acid transport and metabolism, mainly in susceptible BALB/c-infected macrophages. Moreover, we show how parasite arginase activity affects parasite gene expression modulation, including amino acid uptake and amastin expression.

Purine-responsive expression of the Leishmania donovani NT3 purine nucleobase transporter is mediated by a conserved RNA stem-loop

The ability to modulate gene expression in response to changes in the host environment is essential for survival of the kinetoplastid parasite Leishmania Unlike most eukaryotes, gene expression in kinetoplastids is predominately regulated posttranscriptionally. Consequently, RNA-binding proteins and mRNA-encoded sequence elements serve as primary determinants of gene regulation in these organisms; however, few have defined roles in specific stress response pathways. Leishmania species cannot synthesize purines de novo and must scavenge these essential nutrients from the host. Leishmania have evolved a robust stress response to withstand sustained periods of purine scarcity during their life cycle. The purine nucleobase transporter LdNT3 is among the most substantially up-regulated proteins in purine-starved Leishmania donovani parasites. Here we report that the posttranslational stability of the LdNT3 protein is unchanged in response to purine starvation. Instead, LdNT3 up-regulation is primarily mediated by a 33-nucleotide-long sequence in the LdNT3 mRNA 3' UTR that is predicted to adopt a stem-loop structure. Although this sequence is highly conserved within the mRNAs of orthologous transporters in multiple kinetoplastid species, putative stem-loops from L. donovani and Trypanosoma brucei nucleobase transporter mRNAs were not functionally interchangeable for purine-responsive regulation. Through mutational analysis of the element, we demonstrate that species specificity is attributable to just three variant bases within the predicted loop. Finally, we provide evidence that the abundance of the trans-acting factor that binds the LdNT3 stem-loop in vivo is substantially higher than required for regulation of LdNT3 alone, implying a potential role in regulating other purine-responsive genes.

Uncovering the key dimensions of high-throughput biomolecular data using deep learning

Recent advances in high-throughput single-cell RNA-seq have enabled us to measure thousands of gene expression levels at single-cell resolution. However, the transcriptomic profiles are high-dimensional and sparse in nature. To address it, a deep learning framework based on auto-encoder, termed DeepAE, is proposed to elucidate high-dimensional transcriptomic profiling data in an encode-decode manner. Comparative experiments were conducted on nine transcriptomic profiling datasets to compare DeepAE with four benchmark methods. The results demonstrate that the proposed DeepAE outperforms the benchmark methods with robust performance on uncovering the key dimensions of single-cell RNA-seq data. In addition, we also investigate the performance of DeepAE in other contexts and platforms such as mass cytometry and metabolic profiling in a comprehensive manner. Gene ontology enrichment and pathology analysis are conducted to reveal the mechanisms behind the robust performance of DeepAE by uncovering its key dimensions.

Development and Characterization of an Avirulent Leishmania major Strain

Leishmania major causes cutaneous leishmaniasis. An antileishmanial vaccine for humans is unavailable. In this study, we report development of two attenuated L. major strains-5ASKH-HP and LV39-HP-by continuous culture (high passage) of the corresponding virulent strains (low passage). Both avirulent strains showed similar changes in proteome profiles when analyzed by surface-enhanced laser desorption ionization mass spectrometry. Liquid chromatography-mass spectrometry and microarray characterization of 5ASKH strains revealed substantially altered gene and protein expression profiles, respectively. Both virulent and avirulent L. major strains grew comparably in culture, but the avirulent strain survived significantly less in BALB/c-derived peritoneal macrophages. Both attenuated strains failed to infect BALB/c mice and elicited IFN-γ, but not IL-4 and IL-10, responses. 5ASKH-HP parasites failed to induce significant infection even in severely immunocompromised- SCID or inducible NO synthase-, CD40-, or IL-12-deficient mice, indicating attenuation. The avirulent strain induced less IL-10, but higher IL-12, in macrophages. The avirulent strain failed to reduce CD40 relocation to the detergent-resistant membrane domain and to inhibit CD40-induced phosphorylation of the kinases Lyn and protein kinase C-β and MAPKs MKK-3/6 and p38MAPK or to upregulate MEK-1/2 and ERK-1/2 in BALB/c-derived peritoneal macrophages. The virulent and the avirulent strains reciprocally modulated CD40-induced Ras-mediated signaling through PI-3K and Raf-1. Avirulent 5ASKH-primed BALB/c mice were protected against virulent L. major challenge infection. The loss of virulence accompanied by substantially altered proteome profiles and the elicitation of host-protective immune responses indicate plausibly irreversible attenuation of the L. major strain and its potential use as a vaccine strain.

Genome-Wide Proteomics and Phosphoproteomics Analysis of Leishmania spp. During Differentiation

Determining variations in protein abundance and/or posttranslational modification as a function of time or upon induction by a signal in a particular cell type is central to quantitative proteomics. Isobaric labeling methodologies now allow for parallel quantification of proteins at various conditions concurrently or multiplexing in relatively quantitative proteomics workflows. Hence, mapping the protein expression profiles of various developmental stages of Leishmania parasites is possible with high-resolution mass spectrometry. To analyze global changes in protein expression and cellular signaling pathways during Leishmania differentiation and development is possible with a quantitative proteomics approach. The tandem mass tags (TMT) approach provides a chemical labeling method based on the principle of amine reactive tags; the maximum number of conditions that can be multiplexed is 10-plex. We describe herein a detailed method for sample preparation, TMT-labeling, mass spectrometry and data analysis of different developmental stages of Leishmania donovani parasites. This quantitative proteomic approach is useful to study dynamic changes in protein expression levels during L. donovani differentiation, and also allows in-depth analysis of signaling pathways via phosphoproteomics.

Leishmania braziliensis prostaglandin F2α synthase impacts host infection

Background

Prostaglandins (PG) are lipid mediators derived from arachidonic acid metabolism. They are involved in cellular processes such as inflammation and tissue homeostasis. PG production is not restricted to multicellular organisms. Trypanosomatids also synthesize several metabolites of arachidonic acid. Nevertheless, their biological role in these early-branching parasites and their role in host-parasite interaction are not well elucidated. Prostaglandin F_2α synthase (PGF2S) has been observed in the Leishmania braziliensis secreted proteome and in L. donovani extracellular vesicles. Furthermore, we previously reported a positive correlation between L. braziliensis PGF2S (LbrPGF2S) expression and pathogenicity in mice.

Methods

LbrPGF2S gene expression and PGF2α synthesis in promastigotes were detected and quantified by western blotting and EIA assay kit, respectively. To investigate LbrPGF2S localization in amastigotes during bone marrow-derived macrophage infection, parasites expressing mCherry-LbrPGF2S were generated and followed by time-lapse imaging for 48 h post-infection. PGF2S homolog sequences from Leishmania and humans were analyzed in silico using ClustalW on Geneious v6 and EMBOSS Needle.

Results

Leishmania braziliensis promastigotes synthesize prostaglandin F_2α in the presence of arachidonic acid, with peak production in the stationary growth phase under heat stress. LbrPGF2S is a cytoplasmic protein enriched in the secretory site of the parasite cell body, the flagellar pocket. It is an enzyme constitutively expressed throughout promastigote development, but overexpression of LbrPGF2S leads to an increase of infectivity in vitro. The data suggest that LbrPGF2S may be released from intracellular amastigotes into the cytoplasm of bone marrow-derived macrophages over a 48-hour infection period, using time-lapse microscopy and mCherry-PGF2S (mChPGF2S)-expressing parasites.

Conclusions

LbrPGF2S, a parasite-derived protein, is targeted to the host cell cytoplasm. The putative transfer of this enzyme, involved in pro-inflammatory lipid mediator synthesis, to the host cell suggests a potential role in host-parasite interaction and may partially explain the increased pathogenicity associated with overexpression of LbrPGF2S in L. braziliensis. Our data provide valuable insights to help understand the importance of parasite-derived lipid mediators in pathogenesis.

Differential immune response modulation in early Leishmania amazonensis infection of BALB/c and C57BL/6 macrophages based on transcriptome profiles

The fate of Leishmania infection can be strongly influenced by the host genetic background. In this work, we describe gene expression modulation of the immune system based on dual global transcriptome profiles of bone marrow-derived macrophages (BMDMs) from BALB/c and C57BL/6 mice infected with Leishmania amazonensis. A total of 12,641 host transcripts were identified according to the alignment to the Mus musculus genome. Differentially expressed genes (DEGs) profiling revealed a differential modulation of the basal genetic background between the two hosts independent of L. amazonensis infection. In addition, in response to early L. amazonensis infection, 10 genes were modulated in infected BALB/c vs. non-infected BALB/c macrophages; and 127 genes were modulated in infected C57BL/6 vs. non-infected C57BL/6 macrophages. These modulated genes appeared to be related to the main immune response processes, such as recognition, antigen presentation, costimulation and proliferation. The distinct gene expression was correlated with the susceptibility and resistance to infection of each host. Furthermore, upon comparing the DEGs in BMDMs vs. peritoneal macrophages, we observed no differences in the gene expression patterns of Jun, Fcgr1 and Il1b, suggesting a similar activation trends of transcription factor binding, recognition and phagocytosis, as well as the proinflammatory cytokine production in response to early L. amazonensis infection. Analysis of the DEG profile of the parasite revealed only one DEG among the 8,282 transcripts, indicating that parasite gene expression in early infection does not depend on the host genetic background.

Analysis of mRNA processing at whole transcriptome level, transcriptomic profile and genome sequence refinement of Trypanosoma cruzi

The genomic sequence of Trypanosoma cruzi, the protozoan causative of Chagas disease was published more than a decade ago. However, due to their complexity, its complete haploid predicted sequence and therefore its genetic repertoire remains unconfirmed. In this work, we have used RNAseq data to improve the previous genome assembly of Sylvio X10 strain and to define the complete transcriptome at trypomastigote stage (mammalian stage). A total of 22,977 transcripts were identified, of which more than half could be considered novel as they did not match previously annotated genes. Moreover, for the first time in T. cruzi, we are providing their relative abundance levels. We have identified that Sylvio X10 trypomastigotes exhibit a predominance of surface protein genes, specifically those encoding trans-sialidase and mucin-like proteins. On the other hand, detailed analysis of the pre-mRNA processing sites revealed some similarities but also some differences in the spliced leader and different polyadenylation addition sites compared to close related kinetoplastid parasites. Our results also confirm that transcription is bidirectional as occur in other kinetoplastids and the proportion of forward-sense and reverse-sense transcripts is almost equivalent, demonstrating that a strand-specificity does not exist.

Heme synthesis through the life cycle of the heme auxotrophic parasite Leishmania major

Heme is an essential molecule synthetized through a broadly conserved 8-step route that has been lost in trypanosomatid parasites. Interestingly, Leishmania reacquired by horizontal gene transfer from γ-proteobacteria the genes coding for the last 3 enzymes of the pathway. Here we show that intracellular amastigotes of Leishmania major can scavenge heme precursors from the host cell to fulfill their heme requirements, demonstrating the functionality of this partial pathway. To dissect its role throughout the L. major life cycle, the significance of L. major ferrochelatase (LmFeCH), the terminal enzyme of the route, was evaluated. LmFeCH expression in a heterologous system demonstrated its activity. Knockout promastigotes lacking lmfech were not able to use the ferrochelatase substrate protoporphyrin IX as a source of heme. In vivo infection of Phlebotomus perniciosus with knockout promastigotes shows that LmFeCH is not required for their development in the sandfly. In contrast, the replication of intracellular amastigotes was hampered in vitro by the deletion of lmfech. However, LmFeCH^-/- parasites produced disease in a cutaneous leishmaniasis murine model in a similar way as control parasites. Therefore, although L. major can synthesize de novo heme from macrophage precursors, this activity is dispensable being an unsuited target for leishmaniasis treatment.-Orrego, L. M., Cabello-Donayre, M., Vargas, P., Martínez-García, M., Sánchez, C., Pineda-Molina, E., Jiménez, M., Molina, R., Pérez-Victoria, J. M. Heme synthesis through the life cycle of the heme auxotrophic parasite Leishmania major.

Transcriptional responses of Leishmania (Leishmania) amazonensis in the presence of trivalent sodium stibogluconate

Background

In the last decade, resistance to antimonials has become a serious problem due to the emergence of drug-resistant strains. Therefore, understanding the mechanisms used by Leishmania parasites to survive under drug pressure is essential, particularly for species of medical-veterinary importance such as L. amazonensis.

Methods

Here, we used RNA-seq technology to analyse transcriptome profiles and identify global changes in gene expression between antimony-resistant and -sensitive L. amazonensis promastigotes.

Results

A total of 723 differentially expressed genes were identified between resistant and sensitive lines. Comparative transcriptomic analysis revealed that genes encoding proteins involved in metabolism (fatty acids) and stress response, as well as those associated with antimony resistance in other Leishmania species, were upregulated in the antimony-resistant line. Most importantly, we observed upregulation of genes encoding autophagy proteins, suggesting that in the presence of trivalent stibogluconate (Sb^III) L. amazonensis can activate these genes either as a survival strategy or to induce cell death, as has been observed in other parasites.

Conclusions

This work identified global transcriptomic changes in an in vitro-adapted strain in response to Sb^III. Our results provide relevant information to continue understanding the mechanism used by parasites of the subgenus Leishmania (L. amazonensis) to generate an antimony-resistant phenotype.

Electronic supplementary material

The online version of this article (10.1186/s13071-019-3603-8) contains supplementary material, which is available to authorized users.

Genome and transcriptome analyses of Leishmania spp.: opening Pandora's box

In the last 30 years, significant advances in genetic manipulation tools along with complete genome and transcriptome sequencing have advanced our understanding of the biology of Leishmania parasites and their interplay with the sand fly and mammalian hosts. High-throughput sequencing in association with CRISPR/Cas9 have prepared the ground for significant advances. Given the richness of the progress made over the last decade, in this article, we focused on the most recent contributions of genome-wide and transcriptome analyses of Leishmania spp., which permit the comparison of life cycle stages, the evaluation of different strains and species in their natural niches and in the field and the simultaneously comparison of the gene expression profiles of parasites and hosts.

Regulation of gene expression in trypanosomatids: living with polycistronic transcription

In trypanosomes, RNA polymerase II transcription is polycistronic and individual mRNAs are excised by trans-splicing and polyadenylation. The lack of individual gene transcription control is compensated by control of mRNA processing, translation and degradation. Although the basic mechanisms of mRNA decay and translation are evolutionarily conserved, there are also unique aspects, such as the existence of six cap-binding translation initiation factor homologues, a novel decapping enzyme and an mRNA stabilizing complex that is recruited by RNA-binding proteins. High-throughput analyses have identified nearly a hundred regulatory mRNA-binding proteins, making trypanosomes valuable as a model system to investigate post-transcriptional regulation.

Functional genomics in sand fly-derived Leishmania promastigotes

Background

Leishmania development in the sand fly gut leads to highly infective forms called metacyclic promastigotes. This process can be routinely mimicked in culture. Gene expression–profiling studies by transcriptome analysis have been performed with the aim of studying promastigote forms in the sand fly gut, as well as differences between sand fly–and culture-derived promastigotes.

Findings

Transcriptome analysis has revealed the crucial role of the microenvironment in parasite development within the sand fly gut because substantial differences and moderate correlation between the transcriptomes of cultured and sand fly–derived promastigotes have been found. Sand fly–derived metacyclics are more infective than metacyclics in culture. Therefore, some caution should be exercised when using cultured promastigotes, depending on the experimental design. The most remarkable examples are the hydrophilic acidic surface protein/small endoplasmic reticulum protein (HASP/SHERP) cluster, the glycoprotein 63 (gp63), and autophagy genes, which are up-regulated in sand fly–derived promastigotes compared with cultured promastigotes. Because HASP/SHERP genes are up-regulated in nectomonad and metacyclic promastigotes in the sand fly, the encoded proteins are not metacyclic specific. Metacyclic promastigotes are distinguished by morphology and high infectivity. Isolating them from the sand fly gut is not exempt from technical difficulty, because other promastigote forms remain in the gut even 15 days after infection. Leishmania major procyclic promastigotes within the sand fly gut up-regulate genes involved in cell cycle regulation and glucose catabolism, whereas metacyclics increase transcript levels of fatty acid biosynthesis and ATP-coupled proton transport genes. Most parasite's signal transduction pathways remain uncharacterized. Future elucidation may improve understanding of parasite development, particularly signaling molecule-encoding genes in sand fly versus culture and between promastigote forms in the sand fly gut.

Conclusions

Transcriptome analysis has been demonstrated to be technically efficacious to study differential gene expression in sand fly gut promastigote forms. Transcript and protein levels are not well correlated in these organisms (approximately 25% quantitative coincidences), especially under stress situations and at differentiation processes. However, transcript and protein levels behave similarly in approximately 60% of cases from a qualitative point of view (increase, decrease, or no variation). Changes in translational efficiency observed in other trypanosomatids strongly suggest that the differences are due to translational regulation and regulation of the steady-state protein levels. The lack of low-input sample strategies does not allow translatome and proteome analysis of sand fly–derived promastigotes so far.

Comparative transcriptomics in Leishmania braziliensis: disclosing differential gene expression of coding and putative noncoding RNAs across developmental stages

Leishmaniasis is a worldwide public health problem caused by protozoan parasites of the genus Leishmania. Leishmania braziliensis is the most important species responsible for tegumentary leishmaniases in Brazil. An understanding of the molecular mechanisms underlying the success of this parasite is urgently needed. An in-depth study on the modulation of gene expression across the life cycle stages of L. braziliensis covering coding and noncoding RNAs (ncRNAs) was missing and is presented herein. Analyses of differentially expressed (DE) genes revealed that most prominent differences were observed between the transcriptomes of insect and mammalian proliferative forms (6,576 genes). Gene ontology (GO) analysis indicated stage-specific enriched biological processes. A computational pipeline and 5 ncRNA predictors allowed the identification of 11,372 putative ncRNAs. Most of the DE ncRNAs were found between the transcriptomes of insect and mammalian proliferative stages (38%). Of the DE ncRNAs, 295 were DE in all three stages and displayed a wide range of lengths, chromosomal distributions and locations; many of them had a distinct expression profile compared to that of their protein-coding neighbors. Thirty-five putative ncRNAs were submitted to northern blotting analysis, and one or more hybridization-positive signals were observed in 22 of these ncRNAs. This work presents an overview of the L. braziliensis transcriptome and its adjustments throughout development. In addition to determining the general features of the transcriptome at each life stage and the profile of protein-coding transcripts, we identified and characterized a variety of noncoding transcripts. The novel putative ncRNAs uncovered in L. braziliensis might be regulatory elements to be further investigated.

UTRme: A Scoring-Based Tool to Annotate Untranslated Regions in Trypanosomatid Genomes

Most signals involved in post-transcriptional regulatory networks are located in the untranslated regions (UTRs) of the mRNAs. Therefore, to deepen our understanding of gene expression regulation, delimitation of these regions with high accuracy is needed. The trypanosomatid lineage includes a variety of parasitic protozoans causing a significant worldwide burden on human health. Given their peculiar mechanisms of gene expression, these organisms depend on post-transcriptional regulation as the main level of gene expression control. In this context, the definition of the UTR regions becomes of key importance. We have developed UTR-mini-exon (UTRme), a graphical user interface (GUI) stand-alone application to identify and annotate 5′ and 3′ UTR regions in a highly accurate way. UTRme implements a multiple scoring system tailored to address the issue of false positive UTR assignment that frequently arise because of the characteristics of the intergenic regions. Even though it was developed for trypanosomatids, the tool can be used to predict 3′ sites in any eukaryote and 5′ UTRs in any organism where trans-splicing occurs (such as the model organism C. elegans). UTRme offers a way for non-bioinformaticians to precisely determine UTRs from transcriptomic data. The tool is freely available via the conda and github repositories.

Prediction of Leishmania major Key Proteins Via Topological Analysis of Protein-Protein Interaction Network

Background:

Although leishmaniasis is regarded as a public health problem, no effective vaccine or decisive treatment has been introduced for this disease. Therefore, representing novel therapeutic proteins is essential. Protein-protein Interaction network analysis is a suitable tool to discover novel drug targets for leishmania major. To this aim, gene and protein expression data is used for instructing protein network and the key proteins are highlighted.

Materials and Methods:

In this computational and bioinformatics study, the protein/gene expression data related to leishmania major were studied, and 252 candidate proteins were extracted. Then, the protein networks of these proteins were explored and visualized by using String database and Cytoscape software. Finally, clustering and gene ontology were performed by MCODE and PANTHER databases, respectively.

Results:

Based on gene ontology analysis, most of the leishmania major proteins were located in cell compartments and membrane. Catalytic activity and binding were regarded as the relevant molecular functions and metabolic and cellular processes were the significant biological process. In this network analysis, UB-EP52, EF-2, chaperonin, Hsp70.4, Hsp60, tubulin alpha and beta chain, and ENOL and LACK were introduced as hub-bottleneck proteins. Based on clustering analysis, Lmjf.32.3270, ENOL and Lmjf.13.0290 were determined as seed proteins in each cluster.

Conclusion:

The results indicated that hub proteins play a significant role in pathogenesis and life cycle of leishmania major. Further studies of hubs will provide a better understanding of leishmaniasis mechanisms. Finally, these key hub proteins could be a suitable and helpful potential for drug targets and treating leishmaniasis by considering their validation.

An ex vivo abomasal ovine model to study the immediate immune response in the context of Haemonchus contortus larval-stage

We have set up an ex vivo ovine abomasal model, which can mimic the multicellular process to explore the early steps in haemonchine nematode infection using RNA-seq technology. Ovine abomasal explants were collected for histological and transcriptional analysis and supernatants collected to quantitate lactate dehydrogenase (LDH) enzymes. Atotal of 233 were substantially induced genes between L₄-inoculated and uninoculated-control tissues, respectively. However, a total of 14 were considerably down-regulated genes between the 51 aforementioned tissues. Fifteen pathways were annotated by Kyoto Encyclopedia of Genes, and Genomes pathway analysis accounted for the significant percentage in immediate response to larval-stage of H. contortus. Key genes upregulated in response to the addition of L₄-inoculum of H. contortus were IL-6, IL-8, C1q, Atypical chemokine receptor-3, chemokine ligand-2, manganese superoxide dismutase, integrin alpha-7, -8, -9, integrin subunit beta-1, integrin subunit beta 6, intercellular adhesion molecule-1 and actin alpha-1. This study shows for the first time that galectin-1 is up-regulated in an ex vivo abomasal segment model exposed to L₄-inoculum of H. contortus following 6 h of incubation. The abomasal segment model has been shown to be a suitable tool to study the haemonchine larval-stage effects on the ovine abomasal tissues prior to in vivo assessment.

Function of Macrophage and Parasite Phosphatases in Leishmaniasis

The kinetoplastid protozoan parasites belonging to the genus Leishmania are the causative agents of different clinical forms of leishmaniasis, a vector-borne infectious disease with worldwide prevalence. The protective host immune response against Leishmania parasites relies on myeloid cells such as dendritic cells and macrophages in which upon stimulation by cytokines (e.g., interferon-γ) a complex network of signaling pathways is switched on leading to strong antimicrobial activities directed against the intracellular parasite stage. The regulation of these pathways classically depends on post-translational modifications of proteins, with phosphorylation events playing a cardinal role. Leishmania parasites deactivate their phagocytic host cells by inducing specific mammalian phosphatases that are capable to impede signaling. On the other hand, there is now also evidence that Leishmania spp. themselves express phosphatases that might target host cell molecules and thereby facilitate the intracellular survival of the parasite. This review will present an overview on the modulation of host phosphatases by Leishmania parasites as well as on the known families of Leishmania phosphatases and their possible function as virulence factors. A more detailed understanding of the role of phosphatases in Leishmania–host cell interactions might open new avenues for the treatment of non-healing, progressive forms of leishmaniasis.

Characterisation of Casein Kinase 1.1 in Leishmania donovani Using the CRISPR Cas9 Toolkit

The recent adaptation of CRISPR Cas9 genome editing to Leishmania spp. has opened a new era in deciphering Leishmania biology. The method was recently improved using a PCR-based CRISPR Cas9 approach, which eliminated the need for cloning. This new approach, which allows high-throughput gene deletion, was successfully validated in L. mexicana and L. major. In this study, we validated the toolkit in Leishmania donovani targeting the flagellar protein PF16, confirming that the tagged protein localizes to the flagellum and that null mutants lose their motility. We then used the technique to characterise CK1.1, a member of the casein kinase 1 family, which is involved in the regulation of many cellular processes. We showed that CK1.1 is a low-abundance protein present in promastigotes and in amastigotes. We demonstrated that CK1.1 is not essential for promastigote and axenic amastigote survival or for axenic amastigote differentiation, although it may have a role during stationary phase. Altogether, our data validate the use of PCR-based CRISPR Cas9 toolkit in L. donovani, which will be crucial for genetic modification of hamster-derived, disease-relevant parasites.

RNA-seq transcriptional profiling of Leishmania amazonensis reveals an arginase-dependent gene expression regulation

Background

Leishmania is a protozoan parasite that alternates its life cycle between the sand-fly vector and the mammalian host. This alternation involves environmental changes and leads the parasite to dynamic modifications in morphology, metabolism, cellular signaling and regulation of gene expression to allow for a rapid adaptation to new conditions. The L-arginine pathway in L. amazonensis is important during the parasite life cycle and interferes in the establishment and maintenance of the infection in mammalian macrophages. Host arginase is an immune-regulatory enzyme that can reduce the production of nitric oxide by activated macrophages, directing the availability of L-arginine to the polyamine pathway, resulting in parasite replication. In this work, we performed transcriptional profiling to identify differentially expressed genes in L. amazonensis wild-type (La-WT) versus L. amazonensis arginase knockout (La-arg^-) promastigotes and axenic amastigotes.

Methodology/Principal findings

A total of 8253 transcripts were identified in La-WT and La-arg^- promastigotes and axenic amastigotes, about 60% of them codifying hypothetical proteins and 443 novel transcripts, which did not match any previously annotated genes. Our RNA-seq data revealed that 85% of genes were constitutively expressed. The comparison of transcriptome and metabolome data showed lower levels of arginase and higher levels of glutamate-5-kinase in La-WT axenic amastigotes compared to promastigotes. The absence of arginase activity in promastigotes increased the levels of pyrroline 5-carboxylate reductase, but decreased the levels of arginosuccinate synthase, pyrroline 5-carboxylate dehydrogenase, acetylornithine deacetylase and spermidine synthase transcripts levels. These observations can explain previous metabolomic data pointing to the increase of L-arginine, citrulline and L-glutamate and reduction of aspartate, proline, ornithine and putrescine. Altogether, these results indicate that arginase activity is important in Leishmania gene expression modulation during differentiation and adaptation to environmental changes. Here, we confirmed this hypothesis with the identification of differential gene expression of the enzymes involved in biosynthesis of amino acids, arginine and proline metabolism and arginine biosynthesis.

Conclusions/Significance

All data provided information about the transcriptomic profiling and the expression levels of La-WT and La-arg^- promastigotes and axenic amastigotes. These findings revealed the importance of arginase in parasite survival and differentiation, and indicated the existence of a coordinated response in the absence of arginase activity related to arginine and polyamine pathways.

Leishmania are auxotrophic for many essential nutrients, including amino acids. In this way, the parasite needs to uptake the amino acids from the environment. The uptake of amino acids is mediated by amino acid transporters that are unique for Leishmania. As part of polyamine pathway, the arginase converts L-arginine to ornithine and furthermore to putrescine, products which are essential for parasite growth. On the other hand, the absence of arginase activity could alter the metabolism of the parasite to surpass the external signals during the life cycle and the fate of infection. The transcriptional profiling of La-WT and La-arg^- promastigotes and axenic amastigotes revealed 8253 transcripts, 60% encoding hypothetical proteins and 443 novel transcripts. In addition, our data revealed that 85% of the genes were constitutively expressed. Among the 15% (1268 genes) of the differentially expressed genes, we identified genes up- and down-regulated comparing the transcript abundance from different life cycle stages of the parasite and in the presence or absence of arginase. We also combined the transcriptional with metabolic profile that revealed a proportional correlation between enzyme and metabolites in the polyamine pathway. The differentiation of promastigotes to amastigotes alters the expression of enzymes from polyamines biosynthesis, which modulates ornithine, L-glutamate, proline and putrescine levels. In addition, the absence of arginase activity increased the levels of L-arginine, citrulline and L-glutamate and decreased the levels of aspartate, proline, ornithine and putrescine in promastigotes by differential modulation of genes involved in its metabolism. Altogether these data provided additional insights into how Leishmania is able to modulate its biological functions in the presence or absence of arginase activity to survive during environmental changes.

Characterization of the Protein Tyrosine Phosphatase LmPRL-1 Secreted by Leishmania major via the Exosome Pathway

Similar to other intracellular pathogens, Leishmania parasites are known to evade the antimicrobial effector functions of host immune cells. To date, however, only a few virulence factors have been described for Leishmania major, one of the causative agents of cutaneous leishmaniasis. Here, we have characterized the expression and function of an L. major phosphatase, which we termed LmPRL-1. This enzyme shows a strong structural similarity to the human phosphatases of regenerating liver (PRL-1, -2, and -3) that regulate the proliferation, differentiation, and motility of cells. The biochemical characterization of the L. major phosphatase revealed that the enzyme is redox sensitive. When analyzing the subcellular localization of LmPRL-1 in promastigotes, amastigotes, and infected macrophages, we found that the phosphatase was predominantly expressed and secreted by promastigotes via the exosome route. Finally, we observed that ectopic expression of LmPRL-1 in L. major led to an increased number of parasites in macrophages. From these data, we conclude that the L. major phosphatase LmPRL-1 contributes to the intracellular survival of the parasites in macrophages.

Multiplexed Spliced-Leader Sequencing: A high-throughput, selective method for RNA-seq in Trypanosomatids

High throughput sequencing techniques are poorly adapted for in vivo studies of parasites, which require prior in vitro culturing and purification. Trypanosomatids, a group of kinetoplastid protozoans, possess a distinctive feature in their transcriptional mechanism whereby a specific Spliced Leader (SL) sequence is added to the 5'end of each mRNA by trans-splicing. This allows to discriminate Trypansomatid RNA from mammalian RNA and forms the basis of our new multiplexed protocol for high-throughput, selective RNA-sequencing called SL-seq. We provided a proof-of-concept of SL-seq in Leishmania donovani, the main causative agent of visceral leishmaniasis in humans, and successfully applied the method to sequence Leishmania mRNA directly from infected macrophages and from highly diluted mixes with human RNA. mRNA profiles obtained with SL-seq corresponded largely to those obtained from conventional poly-A tail purification methods, indicating both enumerate the same mRNA pool. However, SL-seq offers additional advantages, including lower sequencing depth requirements, fast and simple library prep and high resolution splice site detection. SL-seq is therefore ideal for fast and massive parallel sequencing of parasite transcriptomes directly from host tissues. Since SLs are also present in Nematodes, Cnidaria and primitive chordates, this method could also have high potential for transcriptomics studies in other organisms.

Early and late gene expression profiles of the ovine mucosa in response to Haemonchus contortus infection employing Illumina RNA-seq technology

We conducted herein transcriptome sequencing of the ovine abomasal tissues using the Illumina HiSeq 4000 platform to segregate early and late H. contortus-infected sheep (7 and 50days post-infected groups, respectively) from the control naive ones. A total of 548, 357 and 7 were substantially induced genes in 7days post-infection versus uninfected-control group, 50days post-infection versus 7days post-infection (7dpi), and 50days post-infection (50dpi) versus uninfected-control group, respectively. However, a total of 301, 355 and 11 were significantly repressed genes between 7dpi versus uninfected-control group, 50dpi versus 7dpi, and 50dpi versus uninfected-control group, correspondingly. This indicates that H. contortus infection induced a more potent activation of abomasal gene expression in the early stage of infection as compared to the late stage. Seven pathways were annotated by Kyoto Encyclopedia of Genes, and Genomes pathway analysis accounted for the significant percentage in early H. contortus infection. This study shows for the first time that both galectin-11 and matricellular protein osteopontin are up-regulated in abomasal tissue after chronic H. contortus infection, while galectin-4 is specifically down-regulated in the early infection. Additionally, our results showed that the induction or repression of these molecules is likely to determine the infection progression.

The Transcriptome of Leishmania major Developmental Stages in Their Natural Sand Fly Vector

The life cycle of the Leishmania parasite in the sand fly vector involves differentiation into several distinctive forms, each thought to represent an adaptation to specific microenvironments in the midgut of the fly. Based on transcriptome sequencing (RNA-Seq) results, we describe the first high-resolution analysis of the transcriptome dynamics of four distinct stages of Leishmania major as they develop in a natural vector, Phlebotomus duboscqi. The early transformation from tissue amastigotes to procyclic promastigotes in the blood-fed midgut was accompanied by the greatest number of differentially expressed genes, including the downregulation of amastins, and upregulation of multiple cell surface proteins, sugar and amino acid transporters, and genes related to glucose metabolism and cell cycle progression. The global changes accompanying post-blood meal differentiation of procyclic promastigotes to the nectomonad and metacyclic stages were less extensive, though each displayed a unique signature. The transcriptome of nectomonads, which has not been studied previously, revealed changes consistent with cell cycle arrest and the upregulation of genes associated with starvation and stress, including autophagic pathways of protein recycling. Maturation to the infective, metacyclic stage was accompanied by changes suggesting preadaptation to the intracellular environment of the mammalian host, demonstrated by the amastigote-like profiles of surface proteins and metabolism-related genes. Finally, a direct comparison between sand fly-derived and culture-derived metacyclics revealed a reassuring similarity between the two forms, with the in vivo forms distinguished mainly by a stronger upregulation of transcripts associated with nutrient stress.

The life cycle of Leishmania parasites in the sand fly vector includes their growth and development as morphologically distinct forms of extracellular promastigotes found within the different microenvironments of the gut. Based on RNA-Seq, we provide here the first high-resolution, transcriptomic analysis of Leishmania insect stages during their cyclical development in vivo, from tissue amastigotes ingested with the blood meal to infective, metacyclic promastigotes that initiate infection in the mammalian host. The most extensive genetic reprogramming occurred during the early transformation of amastigotes to rapidly dividing procyclic promastigotes in the blood-fed midgut, with major changes in the abundance of mRNAs for surface proteins and metabolism. The post-blood meal-adapted nectomonad stage was characterized by the downregulation of cell cycle-related genes and the upregulation of stress- and starvation-related genes. Finally, the transcriptome of metacyclic promastigotes shifted to a more amastigote-like profile, suggesting their preadaptation to the intracellular host environment.