Dual Transcriptome Profiling of Leishmania-Infected Human Macrophages Reveals Distinct Reprogramming Signatures

TDR Targets 6: driving drug discovery for human pathogens through intensive chemogenomic data integration

The volume of biological, chemical and functional data deposited in the public domain is growing rapidly, thanks to next generation sequencing and highly-automated screening technologies. These datasets represent invaluable resources for drug discovery, particularly for less studied neglected disease pathogens. To leverage these datasets, smart and intensive data integration is required to guide computational inferences across diverse organisms. The TDR Targets chemogenomics resource integrates genomic data from human pathogens and model organisms along with information on bioactive compounds and their annotated activities. This report highlights the latest updates on the available data and functionality in TDR Targets 6. Based on chemogenomic network models providing links between inhibitors and targets, the database now incorporates network-driven target prioritizations, and novel visualizations of network subgraphs displaying chemical- and target-similarity neighborhoods along with associated target-compound bioactivity links. Available data can be browsed and queried through a new user interface, that allow users to perform prioritizations of protein targets and chemical inhibitors. As such, TDR Targets now facilitates the investigation of drug repurposing against pathogen targets, which can potentially help in identifying candidate targets for bioactive compounds with previously unknown targets. TDR Targets is available at https://tdrtargets.org.

MicroRNAs: Biological Regulators in Pathogen-Host Interactions

An inflammatory response is essential for combating invading pathogens. Several effector components, as well as immune cell populations, are involved in mounting an immune response, thereby destroying pathogenic organisms such as bacteria, fungi, viruses, and parasites. In the past decade, microRNAs (miRNAs), a group of noncoding small RNAs, have emerged as functionally significant regulatory molecules with the significant capability of fine-tuning biological processes. The important role of miRNAs in inflammation and immune responses is highlighted by studies in which the regulation of miRNAs in the host was shown to be related to infectious diseases and associated with the eradication or susceptibility of the infection. Here, we review the biological aspects of microRNAs, focusing on their roles as regulators of gene expression during pathogen-host interactions and their implications in the immune response against Leishmania, Trypanosoma, Toxoplasma, and Plasmodium infectious diseases.

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.

Immuno-metabolic profile of human macrophages after Leishmania and Trypanosoma cruzi infection

Macrophages can reprogram their metabolism in response to the surrounding stimuli, which affects their capacity to kill intracellular pathogens. We have investigated the metabolic and immune status of human macrophages after infection with the intracellular trypanosomatid parasites Leishmania donovani, L. amazonensis and T. cruzi and their capacity to respond to a classical polarizing stimulus (LPS and IFN-γ). We found that macrophages infected with Leishmania preferentially upregulate oxidative phosphorylation, which could be contributed by both host cell and parasite, while T. cruzi infection did not significantly increase glycolysis or oxidative phosphorylation. Leishmania and T. cruzi infect macrophages without triggering a strong inflammatory cytokine response, but infection does not prevent a potent response to LPS and IFN-γ. Infection appears to prime macrophages, since the cytokine response to activation with LPS and IFN-γ is more intense in infected macrophages compared to uninfected ones. Metabolic polarization in macrophages can influence infection and immune evasion of these parasites since preventing macrophage cytokine responses would help parasites to establish a persistent infection. However, macrophages remain responsive to classical inflammatory stimuli and could still trigger inflammatory cytokine secretion by macrophages.

SILAC-based quantitative proteomics reveals pleiotropic, phenotypic modulation in primary murine macrophages infected with the protozoan pathogen Leishmania donovani

Leishmaniases are major vector-borne tropical diseases responsible for great human morbidity and mortality, caused by protozoan, trypanosomatid parasites of the genus Leishmania. In the mammalian host, parasites survive and multiply within mononuclear phagocytes, especially macrophages. However, the underlying mechanisms by which Leishmania spp. affect their host are not fully understood. Herein, proteomic alterations of primary, bone marrow-derived BALB/c macrophages are documented after 72 h of infection with Leishmania donovani insect-stage promastigotes, applying a SILAC-based, quantitative proteomics approach. The protocol was optimised by combining strong anion exchange and gel electrophoresis fractionation that displayed similar depth of analysis (combined total of 6189 mouse proteins). Our analyses revealed 86 differentially modulated proteins (35 showing increased and 51 decreased abundance) in response to Leishmania donovani infection. The proteomics results were validated by analysing the abundance of selected proteins. Intracellular Leishmania donovani infection led to changes in various host cell biological processes, including primary metabolism and catabolic process, with a significant enrichment in lysosomal organisation. Overall, our analysis establishes the first proteome of bona fide primary macrophages infected ex vivo with Leishmania donovani, revealing new mechanisms acting at the host/pathogen interface. SIGNIFICANCE: Little is known on proteome changes that occur in primary macrophages after Leishmania donovani infection. This study describes a SILAC-based quantitative proteomics approach to characterise changes of bone marrow-derived macrophages infected with L. donovani promastigotes for 72 h. With the application of SILAC and the use of SAX and GEL fractionation methods, we have tested new routes for proteome quantification of primary macrophages. The protocols developed here can be applicable to other diseases and pathologies. Moreover, this study sheds important new light on the "proteomic reprogramming" of infected macrophages in response to L. donovani promastigotes that affects primary metabolism, cellular catabolic processes, and lysosomal/vacuole organisation. Thus, our study reveals key molecules and processes that act at the host/pathogen interface that may inform on new immuno- or chemotherapeutic interventions to combat leishmaniasis.

Metabolomic Profile of BALB/c Macrophages Infected with Leishmania amazonensis: Deciphering L-Arginine Metabolism

Background: Leishmaniases are neglected tropical diseases that are caused by Leishmania, being endemic worldwide. L-arginine is an essential amino acid that is required for polyamines production on mammal cells. During Leishmania infection of macrophages, L-arginine is used by host and parasite arginase to produce polyamines, leading to parasite survival; or, by nitric oxide synthase 2 to produce nitric oxide leading to parasite killing. Here, we determined the metabolomic profile of BALB/c macrophages that were infected with L. amazonensis wild type or with L. amazonensis arginase knockout, correlating the regulation of L-arginine metabolism from both host and parasite. Methods: The metabolites of infected macrophages were analyzed by capillary electrophoresis coupled with mass spectrometry (CE-MS). The metabolic fingerprints analysis provided the dual profile from the host and parasite. Results: We observed increased levels of proline, glutamic acid, glutamine, L-arginine, ornithine, and putrescine in infected-L. amazonensis wild type macrophages, which indicated that this infection induces the polyamine production. Despite this, we observed reduced levels of ornithine, proline, and trypanothione in infected-L. amazonensis arginase knockout macrophages, indicating that this infection reduces the polyamine production. Conclusions: The metabolome fingerprint indicated that Leishmania infection alters the L-arginine/polyamines/trypanothione metabolism inside the host cell and the parasite arginase impacts on L-arginine metabolism and polyamine production, defining the infection fate.

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.

Pathogen Evasion of Chemokine Response Through Suppression of CXCL10

Clearance of intracellular pathogens, such as Leishmania (L.) major, depends on an immune response with well-regulated cytokine signaling. Here we describe a pathogen-mediated mechanism of evading CXCL10, a chemokine with diverse antimicrobial functions, including T cell recruitment. Infection with L. major in a human monocyte cell line induced robust CXCL10 transcription without increasing extracellular CXCL10 protein concentrations. We found that this transcriptionally independent suppression of CXCL10 is mediated by the virulence factor and protease, glycoprotein-63 (gp63). Specifically, GP63 cleaves CXCL10 after amino acid A81 at the base of a C-terminal alpha-helix. Cytokine cleavage by GP63 demonstrated specificity, as GP63 cleaved CXCL10 and its homologs, which all bind the CXCR3 receptor, but not distantly related chemokines, such as CXCL8 and CCL22. Further characterization demonstrated that CXCL10 cleavage activity by GP63 was produced by both extracellular promastigotes and intracellular amastigotes. Crucially, CXCL10 cleavage impaired T cell chemotaxis in vitro, indicating that cleaved CXCL10 cannot signal through CXCR3. Ultimately, we propose CXCL10 suppression is a convergent mechanism of immune evasion, as Salmonella enterica and Chlamydia trachomatis also suppress CXCL10. This commonality suggests that counteracting CXCL10 suppression may provide a generalizable therapeutic strategy against intracellular pathogens.

Importance

Leishmaniasis, an infectious disease that annually affects over one million people, is caused by intracellular parasites that have evolved to evade the host's attempts to eliminate the parasite. Cutaneous leishmaniasis results in disfiguring skin lesions if the host immune system does not appropriately respond to infection. A family of molecules called chemokines coordinate recruitment of the immune cells required to eliminate infection. Here, we demonstrate a novel mechanism that Leishmania (L.) spp. employ to suppress host chemokines: a Leishmania-encoded protease cleaves chemokines known to recruit T cells that fight off infection. We observe that other common human intracellular pathogens, including Chlamydia trachomatis and Salmonella enterica, reduce levels of the same chemokines, suggesting a strong selective pressure to avoid this component of the immune response. Our study provides new insights into how intracellular pathogens interact with the host immune response to enhance pathogen survival.

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.

Virulence factor RNA transcript expression in the Leishmania Viannia subgenus: influence of species, isolate source, and Leishmania RNA virus-1

BACKGROUND

Leishmania RNA virus-1 (LRV1) is a double-stranded RNA virus identified in 20-25% of Viannia-species endemic to Latin America, and is believed to accelerate cutaneous to mucosal leishmaniasis over time. Our objective was to quantify known virulence factor (VF) RNA transcript expression according to LRV1 status, causative species, and isolate source.

METHODS

Eight cultured isolates of Leishmania were used, four of which were LRV1-positive (Leishmania Viannia braziliensis [n = 1], L. (V.) guyanensis [n = 1], L. (V.) panamensis [n = 2]), and four were LRV1-negative (L. (V.) panamensis [n = 3], L. (V.) braziliensis [n = 1]). Promastigotes were inoculated into macrophage cultures, and harvested at 24 and 48 h. RNA transcript expression of hsp23, hsp70, hsp90, hsp100, mpi, cpb, and gp63 were quantified by qPCR.

RESULTS

RNA transcript expression of hsp100 (p = 0.012), cpb (p = 0.016), and mpi (p = 0.022) showed significant increases from baseline pure culture expression to 24- and 48-h post-macrophage infection, whereas hsp70 (p = 0.004) was significantly decreased. A trend toward increased transcript expression of hsp100 at baseline in isolates of L. (V.) panamensis was noted. Pooled VF RNA transcript expression by L. (V.) panamensis isolates was lower than that of L. (V.) braziliensis and L. (V.) guyananesis at 24 h (p = 0.03). VF RNA transcript expression did not differ by LRV1 status, or source of cultured isolate at baseline, 24, or 48 h; however, a trend toward increased VF RNA transcript expression of 2.71- and 1.93-fold change of mpi (p = 0.11) and hsp90 (p = 0.11), respectively, in LRV1 negative isolates was noted. Similarly, a trend toward lower levels of overall VF RNA transcript expression in clinical isolates (1.15-fold change) compared to ATCC® strains at 24 h was noted (p = 0.07).

CONCLUSIONS

Our findings suggest that known VF RNA transcript expression may be affected by the process of macrophage infection. We were unable to demonstrate definitively that LRV-1 presence affected VF RNA transcript expression in the species and isolates studied. L. (V.) guyanensis and L. (V.) braziliensis demonstrated higher pooled VF RNA transcript expression than L. (V.) panamensis; however, further analyses of protein expression to corroborate this finding are warranted.

Melatonin and Leishmania amazonensis Infection Altered miR-294, miR-30e, and miR-302d Impacting on Tnf, Mcp-1, and Nos2 Expression

Leishmaniases are neglected diseases that cause a large spectrum of clinical manifestations, from cutaneous to visceral lesions. The initial steps of the inflammatory response involve the phagocytosis of Leishmania and the parasite replication inside the macrophage phagolysosome. Melatonin, the darkness-signaling hormone, is involved in modulation of macrophage activation during infectious diseases, controlling the inflammatory response against parasites. In this work, we showed that exogenous melatonin treatment of BALB/c macrophages reduced Leishmania amazonensis infection and modulated host microRNA (miRNA) expression profile, as well as cytokine production such as IL-6, MCP-1/CCL2, and, RANTES/CCL9. The role of one of the regulated miRNA (miR-294-3p) in L. amazonensis BALB/c infection was confirmed with miRNA inhibition assays, which led to increased expression levels of Tnf and Mcp-1/Ccl2 and diminished infectivity. Additionally, melatonin treatment or miR-30e-5p and miR-302d-3p inhibition increased nitric oxide synthase 2 (Nos2) mRNA expression levels and nitric oxide (NO) production, altering the macrophage activation state and reducing infection. Altogether, these data demonstrated the impact of melatonin treatment on the miRNA profile of BALB/c macrophage infected with L. amazonensis defining the infection outcome.

Host and parasite responses in human diffuse cutaneous leishmaniasis caused by L. amazonensis

Diffuse cutaneous leishmaniasis (DCL) is a rare form of leishmaniasis where parasites grow uncontrolled in diffuse lesions across the skin. Meta-transcriptomic analysis of biopsies from DCL patients infected with Leishmania amazonensis demonstrated an infiltration of atypical B cells producing a surprising preponderance of the IgG4 isotype. DCL lesions contained minimal CD8+ T cell transcripts and no evidence of persistent TH2 responses. Whereas localized disease exhibited activated (so-called M1) macrophage presence, transcripts in DCL suggested a regulatory macrophage (R-Mϕ) phenotype with higher levels of ABCB5, DCSTAMP, SPP1, SLAMF9, PPARG, MMPs, and TM4SF19. The high levels of parasite transcripts in DCL and the remarkable uniformity among patients afforded a unique opportunity to study parasite gene expression in this disease. Patterns of parasite gene expression in DCL more closely resembled in vitro parasite growth in resting macrophages, in the absence of T cells. In contrast, parasite gene expression in LCL revealed 336 parasite genes that were differently upregulated, relative to DCL and in vitro macrophage growth, and these transcripts may represent transcripts that are produced by the parasite in response to host immune pressure.

Model-based Evaluation of Gene Expression Changes in Response to Leishmania Infection

Since the development of high-density microarray technology in the late 1990s, global host gene expression changes in response to various stimuli have been extensively studied. More than a dozen peer-reviewed publications have investigated the effect of Leishmania infection in various models since 2001. This review covers the transcriptional changes in macrophage models induced by various Leishmania species and summarizes the resulting impact these studies have on our understanding of the host response to leishmaniasis in vitro. Characterization of the similarities and differences between various model systems will not only further our understanding of Leishmania-induced changes to macrophage gene expression but also identify potential therapeutic targets in the future.

RNA-Seq Revealed Expression of Many Novel Genes Associated With Leishmania donovani Persistence and Clearance in the Host Macrophage

Host- as well as parasite-specific factors are equally crucial in allowing either the Leishmania parasites to dominate, or host macrophages to resist infection. To identify such factors, we infected murine peritoneal macrophages with either the virulent (vAG83) or the non-virulent (nvAG83) parasites of L. donovani. Then, through dual RNA-seq, we simultaneously elucidated the transcriptomic changes occurring both in the host and the parasites. Through Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of the differentially expressed (DE) genes, we showed that the vAG83-infected macrophages exhibit biased anti-inflammatory responses compared to the macrophages infected with the nvAG83. Moreover, the vAG83-infected macrophages displayed suppression of many important cellular processes, including protein synthesis. Further, through protein-protein interaction study, we showed significant downregulation in the expression of many hubs and hub-bottleneck genes in macrophages infected with vAG83 as compared to nvAG83. Cell signaling study showed that these two parasites activated the MAPK and PI3K-AKT signaling pathways differentially in the host cells. Through gene ontology analyses of the parasite-specific genes, we discovered that the genes for virulent factors and parasite survival were significantly upregulated in the intracellular amastigotes of vAG83. In contrast, genes involved in the immune stimulations, and those involved in negative regulation of the cell cycle and transcriptional regulation, were upregulated in the nvAG83. Collectively, these results depicted a differential regulation in the host and the parasite-specific molecules during in vitro persistence and clearance of the parasites.

Toll-Like Receptor and miRNA-let-7e Expression Alter the Inflammatory Response in Leishmania amazonensis-Infected Macrophages

Parasite recognition by Toll-like receptors (TLRs) contributes to macrophage activation and subsequent control of Leishmania infection through the coordinated production of pro-inflammatory and microbicidal effector molecules. The modulation of microRNA (miRNA) expression by Leishmania infection potentially mediates the post-transcriptional regulation of the expression of genes involved in leishmanicidal activity. Here, the contribution of TLR signaling to the miRNA profile and gene expression was evaluated in Leishmania amazonensis-infected murine macrophages. The infectivity of L. amazonensis was higher in murine bone marrow-derived macrophages from mice knockout for myeloid differentiation factor 88 (MyD88^−/−), TLR2 (TLR2^−/−), or TLR4 (TLR4^−/−) than wild type C57BL/6 (WT). L. amazonensis infection of WT macrophages modulated the expression of 32% of the miRNAs analyzed, while 50% were upregulated. The absence of MyD88, TLR2, and TLR4 altered the percentage of miRNAs modulated during L. amazonensis infection, including the downregulation of let-7e expression. Moreover, the absence of signals mediated by MyD88, TLR2, or TLR4 reduced nitric oxide synthase 2 (Nos2) mRNA expression during infection. Indeed, the inhibition of let-7e increased levels of the Nos2 mRNA and NOS2 (or iNOS) protein and nitric oxide (NO) production in L. amazonensis-infected macrophages (4–24 h). The absence of TLR2 and inhibition of let-7e increased the expression of the arginase 1 (Arg1) mRNA but did not alter the protein level during infection. However, higher levels of the L-arginine transporters Cat2B and Cat1 were detected in the absence of Myd88 signaling during infection but were not altered following let-7e inhibition. The inhibition of let-7e impacted the global expression of genes in the TLR pathway by upregulating the expression of recognition and adaptors molecules, such as Tlr6, Tlr9, Ly96, Tirap, Traf 6, Ticam1, Tollip, Casp8, Map3k1, Mapk8, Nfkbib, Nfkbil1, Ppara, Mapk8ip3, Hspd1, and Ube2n, as well as immunomodulators, such as Ptgs2/Cox2, Csf 2, Csf 3, Ifnb1, Il6ra, and Ilr1, impacting NOS2 expression, NO production and parasite infectiveness. In conclusion, L. amazonensis infection alters the TLR signaling pathways by modulating the expression of miRNAs in macrophages to subvert the host immune responses.

In Search of Biomarkers for Pathogenesis and Control of Leishmaniasis by Global Analyses of Leishmania-Infected Macrophages

Leishmaniasis is a vector-borne, neglected tropical disease with a worldwide distribution that can present in a variety of clinical forms, depending on the parasite species and host genetic background. The pathogenesis of this disease remains far from being elucidated because the involvement of a complex immune response orchestrated by host cells significantly affects the clinical outcome. Among these cells, macrophages are the main host cells, produce cytokines and chemokines, thereby triggering events that contribute to the mediation of the host immune response and, subsequently, to the establishment of infection or, alternatively, disease control. There has been relatively limited commercial interest in developing new pharmaceutical compounds to treat leishmaniasis. Moreover, advances in the understanding of the underlying biology of Leishmania spp. have not translated into the development of effective new chemotherapeutic compounds. As a result, biomarkers as surrogate disease endpoints present several potential advantages to be used in the identification of targets capable of facilitating therapeutic interventions considered to ameliorate disease outcome. More recently, large-scale genomic and proteomic analyses have allowed the identification and characterization of the pathways involved in the infection process in both parasites and the host, and these analyses have been shown to be more effective than studying individual molecules to elucidate disease pathogenesis. RNA-seq and proteomics are large-scale approaches that characterize genes or proteins in a given cell line, tissue, or organism to provide a global and more integrated view of the myriad biological processes that occur within a cell than focusing on an individual gene or protein. Bioinformatics provides us with the means to computationally analyze and integrate the large volumes of data generated by high-throughput sequencing approaches. The integration of genomic expression and proteomic data offers a rich multi-dimensional analysis, despite the inherent technical and statistical challenges. We propose that these types of global analyses facilitate the identification, among a large number of genes and proteins, those that hold potential as biomarkers. The present review focuses on large-scale studies that have identified and evaluated relevant biomarkers in macrophages in response to Leishmania infection.

Transcriptional Profiling of Immune-Related Genes in Leishmania infantum-Infected Mice: Identification of Potential Biomarkers of Infection and Progression of Disease

Leishmania spp. is a protozoan parasite that affects millions of people around the world. At present, there is no effective vaccine to prevent leishmaniases in humans. A major limitation in vaccine development is the lack of precise understanding of the particular immunological mechanisms that allow parasite survival in the host. The parasite-host cell interaction induces dramatic changes in transcriptome patterns in both organisms, therefore, a detailed analysis of gene expression in infected tissues will contribute to the evaluation of drug and vaccine candidates, the identification of potential biomarkers, and the understanding of the immunological pathways that lead to protection or progression of disease. In this large-scale analysis, differential expression of 112 immune-related genes has been analyzed using high-throughput qPCR in spleens of infected and naïve Balb/c mice at four different time points. This analysis revealed that early response against Leishmania infection is characterized by the upregulation of Th1 markers and M1-macrophage activation molecules such as Ifng, Stat1, Cxcl9, Cxcl10, Ccr5, Cxcr3, Xcl1, and Ccl3. This activation doesn't protect spleen from infection, since parasitic burden rises along time. This marked difference in gene expression between infected and control mice disappears during intermediate stages of infection, probably related to the strong anti-inflammatory and immunosuppresory signals that are activated early upon infection (Ctla4) or remain activated throughout the experiment (Il18bp). The overexpression of these Th1/M1 markers is restored later in the chronic phase (8 wpi), suggesting the generation of a classical "protective response" against leishmaniasis. Nonetheless, the parasitic burden rockets at this timepoint. This apparent contradiction can be explained by the generation of a regulatory immune response characterized by overexpression of Ifng, Tnfa, Il10, and downregulation Il4 that counteracts the Th1/M1 response. This large pool of data was also used to identify potential biomarkers of infection and parasitic burden in spleen, on the bases of two different regression models. Given the results, gene expression signature analysis appears as a useful tool to identify mechanisms involved in disease outcome and to establish a rational approach for the identification of potential biomarkers useful for monitoring disease progression, new therapies or vaccine development.

UDP-glycosyltransferase genes in trypanosomatid genomes have diversified independently to meet the distinct developmental needs of parasite adaptations

BACKGROUND

Trypanosomatid parasites such as Trypanosoma spp. and Leishmania spp. are a major source of infectious disease in humans and domestic animals worldwide. Fundamental to the host-parasite interactions of these potent pathogens are their cell surfaces, which are highly decorated with glycosylated proteins and other macromolecules. Trypanosomatid genomes contain large multi-copy gene families encoding UDP-dependent glycosyltransferases (UGTs), the primary role of which is cell-surface decoration. Here we report a phylogenetic analysis of UGTs from diverse trypanosomatid genomes, the aim of which was to understand the origin and evolution of their diversity.

RESULTS

By combining phylogenetics with analyses of recombination, and selection, we compared UGT repertoire, genomic context and sequence evolution across 19 trypanosomatids. We identified a UGT lineage present in stercorarian trypanosomes and a free-living kinetoplastid Bodo saltans that likely represents the ancestral state of this gene family. The phylogeny of parasite-specific genes shows that UGTs repertoire in Leishmaniinae and salivarian trypanosomes has expanded independently and with distinct evolutionary dynamics. In the former, the ancestral UGT repertoire was organised in a tandem array from which sporadic transpositions to telomeric regions occurred, allowing expansion most likely through telomeric exchange. In the latter, the ancestral UGT repertoire was comprised of seven subtelomeric lineages, two of which have greatly expanded potentially by gene transposition between these dynamic regions of the genome.

CONCLUSIONS

The phylogeny of UGTs confirms that they represent a substantial parasite-specific innovation, which has diversified independently in the distinct trypanosomatid lineages. Nonetheless, developmental regulation has been a strong driver of UGTs diversification in both African trypanosomes and Leishmania.

Mycobacterium tuberculosis Exploits a Molecular Off Switch of the Immune System for Intracellular Survival

Mycobacterium tuberculosis (M. tuberculosis) survives and multiplies inside human macrophages by subversion of immune mechanisms. Although these immune evasion strategies are well characterised functionally, the underlying molecular mechanisms are poorly understood. Here we show that during infection of human whole blood with M. tuberculosis, host gene transcriptional suppression, rather than activation, is the predominant response. Spatial, temporal and functional characterisation of repressed genes revealed their involvement in pathogen sensing and phagocytosis, degradation within the phagolysosome and antigen processing and presentation. To identify mechanisms underlying suppression of multiple immune genes we undertook epigenetic analyses. We identified significantly differentially expressed microRNAs with known targets in suppressed genes. In addition, after searching regions upstream of the start of transcription of suppressed genes for common sequence motifs, we discovered novel enriched composite sequence patterns, which corresponded to Alu repeat elements, transposable elements known to have wide ranging influences on gene expression. Our findings suggest that to survive within infected cells, mycobacteria exploit a complex immune “molecular off switch” controlled by both microRNAs and Alu regulatory elements.

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.

Trypanosoma cruzi specific mRNA amplification by in vitro transcription improves parasite transcriptomics in host-parasite RNA mixtures

Background

Trypanosomatids are a group of protozoan parasites that includes the etiologic agents of important human illnesses as Chagas disease, sleeping sickness and leishmaniasis. These parasites have a significant distinction from other eukaryotes concerning mRNA structure, since all mature mRNAs have an identical species-specific sequence of 39 nucleotides at the 5′ extremity, named spliced leader (SL). Considering this peculiar aspect of trypanosomatid mRNA, the aim of the present work was to develop a Trypanosoma cruzi specific in vitro transcription (IVT) linear mRNA amplification method in order to improve parasite transcriptomics analyses.

Methods

We designed an oligonucleotide complementary to the last 21 bases of T. cruzi SL sequence, bearing an upstream T7 promoter (T7SL primer), which was used to direct the synthesis of second-strand cDNA. Original mRNA was then amplified by IVT using T7 RNA polymerase. T7SL-amplified RNA from two distinct T. cruzi stages (epimastigotes and trypomastigotes) were deep sequenced in SOLiD platform. Usual poly(A) + RNA and and T7-oligo(dT) amplified RNA (Eberwine method) were also sequenced. RNA-Seq reads were aligned to our new and improved T. cruzi Dm28c genome assembly (PacBio technology) and resulting transcriptome pattern from these three RNA preparation methods were compared, mainly concerning the conservation of mRNA transcritional levels and DEGs detection between epimastigotes and trypomastigotes.

Results

T7SL IVT method detected more potential differentially expressed genes in comparison to either poly(A) + RNA or T7dT IVT, and was also able to produce reliable quantifications of the parasite transcriptome down to 3 ng of total RNA. Furthermore, amplification of parasite mRNA in HeLa/epimastigote RNA mixtures showed that T7SL IVT generates transcriptome quantification with similar detection of differentially expressed genes when parasite RNA mass was only 0.1% of the total mixture (R = 0.78 when compared to poly(A) + RNA).

Conclusions

The T7SL IVT amplification method presented here allows the detection of more potential parasite differentially expressed genes (in comparison to poly(A) + RNA) in host-parasite mixtures or samples with low amount of RNA. This method is especially useful for trypanosomatid transcriptomics because it produces less bias than PCR-based mRNA amplification. Additionally, by simply changing the complementary region of the T7SL primer, the present method can be applied to any trypanosomatid species.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4163-y) contains supplementary material, which is available to authorized users.

Dual RNA-seq reveals no plastic transcriptional response of the coccidian parasite Eimeria falciformis to host immune defenses

Background

Parasites can either respond to differences in immune defenses that exist between individual hosts plastically or, alternatively, follow a genetically canalized (“hard wired”) program of infection. Assuming that large-scale functional plasticity would be discernible in the parasite transcriptome we have performed a dual RNA-seq study of the lifecycle of Eimeria falciformis using infected mice with different immune status as models for coccidian infections.

Results

We compared parasite and host transcriptomes (dual transcriptome) between naïve and challenge infected mice, as well as between immune competent and immune deficient ones. Mice with different immune competence show transcriptional differences as well as differences in parasite reproduction (oocyst shedding). Broad gene categories represented by differently abundant host genes indicate enrichments for immune reaction and tissue repair functions. More specifically, TGF-beta, EGF, TNF and IL-1 and IL-6 are examples of functional annotations represented differently depending on host immune status. Much in contrast, parasite transcriptomes were neither different between Coccidia isolated from immune competent and immune deficient mice, nor between those harvested from naïve and challenge infected mice. Instead, parasite transcriptomes have distinct profiles early and late in infection, characterized largely by biosynthesis or motility associated functional gene groups, respectively. Extracellular sporozoite and oocyst stages showed distinct transcriptional profiles and sporozoite transcriptomes were found enriched for species specific genes and likely pathogenicity factors.

Conclusion

We propose that the niche and host-specific parasite E. falciformis uses a genetically canalized program of infection. This program is likely fixed in an evolutionary process rather than employing phenotypic plasticity to interact with its host. This in turn might limit the potential of the parasite to adapt to new host species or niches, forcing it to coevolve with its host.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4095-6) contains supplementary material, which is available to authorized users.

Leishmania major large RAB GTPase is highly immunogenic in individuals immune to cutaneous and visceral leishmaniasis

Background

We previously identified a Leishmania (L.) major large RAB GTPase (LmlRAB), a new atypical RAB GTPase protein. It is highly conserved in Leishmania species while displaying low level of homology with mammalian homologues. Leishmania small RAB GTPases proteins have been involved in regulation of exocytic and endocytic pathways whereas the role of large RAB GTPases proteins has not been characterized yet. We report here the immunogenicity of both recombinant rLmlRAB and rLmlRABC, in individuals with immunity against L. major or L. infantum.

Methods

PBMC were isolated from individuals cured of L. major (CCLm) or from healthy individuals. The latter were subdivided into high or low IFN-γ responders. Healthy high IFN-γ responders, considered as asymptomatics, were living in an endemic area for L. major (HHRLm) or L. infantum (HHRLi). Healthy low IFN-γ responders (HLR) were considered as naïve controls. Cells from all volunteers were stimulated with rLmlRAB or rLmlRABC. Cytokines were analysed by CBA and ELISA and phenotypes of IFN-γ-producing cells were analysed by flow cytometry.

Results

Both rLmlRAB and rLmlRABC induced high significant levels of IFN-γ in CCLm, HHRLm and HHRLi groups. Phenotype analysis of rLmlRAB and rLmlRABC-stimulated T cells in CCLm individuals showed a significant increase in the percentage of specific IFN-γ-producing CD4+ and CD8+ T cells. rLmlRAB induced significant granzyme B levels in CCLm and HHRLm. Low but significant granzyme B levels were detected in naïve group. IL-10 was detected in immune and naïve individuals.

Conclusion

We showed that rLmlRAB protein and its divergent carboxy-terminal part induced a predominant Th1 response in individuals immune to L. major or L. infantum. Our results suggest that rLmlRAB and rLmlRABC proteins are potential cross-species vaccine candidates against cutaneous and visceral leishmaniasis.

Transcriptomic analysis reveals metabolic switches and surface remodeling as key processes for stage transition in Trypanosoma cruzi

American trypanosomiasis is a chronic and endemic disease which affects millions of people. Trypanosoma cruzi, its causative agent, has a life cycle that involves complex morphological and functional transitions, as well as a variety of environmental conditions. This requires a tight regulation of gene expression, which is achieved mainly by post-transcriptional regulation. In this work we conducted an RNAseq analysis of the three major life cycle stages of T. cruzi: amastigotes, epimastigotes and trypomastigotes. This analysis allowed us to delineate specific transcriptomic profiling for each stage, and also to identify those biological processes of major relevance in each state. Stage specific expression profiling evidenced the plasticity of T. cruzi to adapt quickly to different conditions, with particular focus on membrane remodeling and metabolic shifts along the life cycle. Epimastigotes, which replicate in the gut of insect vectors, showed higher expression of genes related to energy metabolism, mainly Krebs cycle, respiratory chain and oxidative phosphorylation related genes, and anabolism related genes associated to nucleotide and steroid biosynthesis; also, a general down-regulation of surface glycoprotein coding genes was seen at this stage. Trypomastigotes, living extracellularly in the bloodstream of mammals, express a plethora of surface proteins and signaling genes involved in invasion and evasion of immune response. Amastigotes mostly express membrane transporters and genes involved in regulation of cell cycle, and also express a specific subset of surface glycoprotein coding genes. In addition, these results allowed us to improve the annotation of the Dm28c genome, identifying new ORFs and set the stage for construction of networks of co-expression, which can give clues about coded proteins of unknown functions.

Primary macrophages and J774 cells respond differently to infection with Mycobacterium tuberculosis

Macrophages play an essential role in the early immune response to Mycobacterium tuberculosis and are the cell type preferentially infected in vivo. Primary macrophages and macrophage-like cell lines are commonly used as infection models, although the physiological relevance of cell lines, particularly for host-pathogen interaction studies, is debatable. Here we use high-throughput RNA-sequencing to analyse transcriptome dynamics of two macrophage models in response to M. tuberculosis infection. Specifically, we study the early response of bone marrow-derived mouse macrophages and cell line J774 to infection with live and γ-irradiated (killed) M. tuberculosis. We show that infection with live bacilli specifically alters the expression of host genes such as Rsad2, Ifit1/2/3 and Rig-I, whose potential roles in resistance to M. tuberculosis infection have not yet been investigated. In addition, the response of primary macrophages is faster and more intense than that of J774 cells in terms of number of differentially expressed genes and magnitude of induction/repression. Our results point to potentially novel processes leading to immune containment early during M. tuberculosis infection, and support the idea that important differences exist between primary macrophages and cell lines, which should be taken into account when choosing a macrophage model to study host-pathogen interactions.

Developmental differentiation in Leishmania lifecycle progression: post-transcriptional control conducts the orchestra

The successful progression of Leishmania spp. through their lifecycle entails a series of differentiation processes; the proliferative procyclic promastigote forms become quiescent, human-infective metacyclic promastigotes during metacyclogenesis in the sandfly vector, which then differentiate into amastigotes during amastigogenesis in the mammalian host. The progression to these infective forms requires two components: environmental cues and a coordinated cellular response. Recent studies have shown that the Leishmania cellular transformation into mammalian-infective stages is triggered by broad changes in the absolute and relative RNA and protein levels. In this review, we will discuss the implications of Leishmania transcriptomic and proteomic fluctuations, which adapt the parasitic cell for survival.

Meta-transcriptome Profiling of the Human-Leishmania braziliensis Cutaneous Lesion

Host and parasite gene expression in skin biopsies from Leishmania braziliensis-infected patients were simultaneously analyzed using high throughput RNA-sequencing. Biopsies were taken from 8 patients with early cutaneous leishmaniasis and 17 patients with late cutaneous leishmaniasis. Although parasite DNA was found in all patient lesions at the time of biopsy, the patients could be stratified into two groups: one lacking detectable parasite transcripts (PT^Neg) in lesions, and another in which parasite transcripts were readily detected (PT^Pos). These groups exhibited substantial differences in host responses to infection. PT^Pos biopsies contained an unexpected increase in B lymphocyte-specific and immunoglobulin transcripts in the lesions, and an upregulation of immune inhibitory molecules. Biopsies without detectable parasite transcripts showed decreased evidence for B cell activation, but increased expression of antimicrobial genes and genes encoding skin barrier functions. The composition and abundance of L. braziliensis transcripts in PT^Pos lesions were surprisingly conserved among all six patients, with minimal meaningful differences between lesions from patients with early and late cutaneous leishmaniasis. The most abundant parasite transcripts expressed in lesions were distinct from transcripts expressed in vitro in human macrophage cultures infected with L. amazonensis or L. major. Therefore in vitro gene expression in macrophage monolayers may not be a strong predictor of gene expression in lesions. Some of the most highly expressed in vivo transcripts encoded amastin-like proteins, hypothetical genes, putative parasite virulence factors, as well as histones and tubulin. In summary, RNA sequencing allowed us to simultaneously analyze human and L. braziliensis transcriptomes in lesions of infected patients, and identify unexpected differences in host immune responses which correlated with active transcription of parasite genes.

Leishmania spp are intracellular protozoan parasites that replicate primarily within host tissue macrophages. In this paper we simultaneously query host and parasite gene expression in human cutaneous L. braziliensis lesions. We observe an unexpectedly prominent role for B cells and immunoglobulins in lesions in which actively transcribing parasites reside. We also observe that parasite gene expression is surprisingly conserved among L. braziliensis lesions, and the genes that are expressed in lesions are not those that have been previously associated with parasite growth in vitro. This analysis of parasite and host gene expression in lesions may lead to the identification of new parasite virulence factors and may identify host responses that promote parasite persistence in lesions.