Leishmania major chromosome 3 contains two long convergent polycistronic gene clusters separated by a tRNA gene

Genome-wide quantification of polycistronic transcription in Leishmania major

Leishmania major is a human-pathogenic, obligate parasite and the etiological agent of the most prevalent, cutaneous form of leishmaniasis, which is an important neglected, tropical disease with ~1.2 million new infections per year. Leishmania, and the whole order Trypanosomatida, are early eukaryotes with highly diverged gene expression and regulation pathways, setting them apart from their mammalian hosts and from most other eukaryotes. Using precision run-on sequence analysis, we performed a genome-wide mapping and density analysis of RNA polymerases in isolated nuclei of the protozoan parasite Leishmania major. We map transcription initiation sites at divergent strand switch regions and head-tail regions within the chromosomes and correlate them with known sites of chromatin modifications. We confirm continuous, polycistronic RNA synthesis in all RNA polymerase II-dependent gene arrays but find small varying RNA polymerase activities in polycistronic transcription units (PTUs), excluding gene-specific transcription regulation, but not PTU-specific variations. Lastly, we find evidence for transcriptional pausing of all three RNA polymerase classes, hinting at a possible mechanism of transcriptional regulation.IMPORTANCELeishmania spp. are pathogens of humans and animals and cause one of the most important neglected tropical diseases. Regulation of gene expression in Leishmania but also in the related Trypanosoma is radically different from all eukaryotic model organisms, dispensing with regulated, gene-specific transcription, and relying instead on highly regulated translation. Our work sheds light on the initiation, elongation, and termination of transcription, maps unidirectional, polycistronic transcription units, provides evidence for transcriptional pausing at or near starting points of RNA synthesis, and quantifies the varying transcription rates of the polycistronic transcription units. Our results will further the understanding of these important pathogens and should provide a valuable resource for researchers in the field of eukaryotic microbiology.

Preliminary Information of Iranian Lizard Leishmania Promastigote Transcriptome Sequencing by Next Generation Sequencing (NGS) Method

Background

A lizard Leishmania has been isolated from a lizard (Agama agilis) in Iran. Its genome sequence has not been determined, so far.

Methods

The study was done at Shahid Beheshti University of Medical Sciences, Tehran, Iran in 2017-2023. Leishmania promastigotes were cultured in RPMI1640 culture medium and collected at logarithmic phase by centrigugation. Parasite RNA was extracted by the Qiagene standard kit and its quantity and quality was determined and sequenced by NGS method with Illumina PE machine at BGI Company (China).

Results

The number of 8316 mRNA, 83 tRNA, 63 rRNA, 83 ncRNA, 5 snRNA, 1039 snoRNA, 36 region, and 3 repeat regions, 8343 CDS, 9597 Exon and 9292 Genes were identified in promastigote of Iranian lizard Leishmania.

Conclusion

Genomic elements of Iranian lizards Leishmania (with unique characteristics) were determined and identified by NGS system.

Making the Most of Its Short Reads: A Bioinformatics Workflow for Analysing the Short-Read-Only Data of Leishmania orientalis (Formerly Named Leishmania siamensis) Isolate PCM2 in Thailand

Simple Summary

Leishmaniasis is a parasitic disease caused by flagellated protozoa of the genus Leishmania. Multiple genome sequencing platforms have been employed to complete Leishmania genomes at the expense of high cost. This study proposes an integrative bioinformatic workflow for assembling only the short-read data of Leishmania orientalis isolate PCM2 from Thailand and produce an acceptable-quality genome for further genomic analysis. This workflow gives extensive information required for identifying strain-specific markers and virulence-associated genes useful for drug and vaccine development before a more exhaustive and expensive investigation.

Abstract

Background: Leishmania orientalis (formerly named Leishmania siamensis) has been neglected for years in Thailand. The genomic study of L. orientalis has gained much attention recently after the release of the first high-quality reference genome of the isolate LSCM4. The integrative approach of multiple sequencing platforms for whole-genome sequencing has proven effective at the expense of considerably expensive costs. This study presents a preliminary bioinformatic workflow including the use of multi-step de novo assembly coupled with the reference-based assembly method to produce high-quality genomic drafts from the short-read Illumina sequence data of L. orientalis isolate PCM2. Results: The integrating multi-step de novo assembly by MEGAHIT and SPAdes with the reference-based method using the L. enriettii genome and salvaging the unmapped reads resulted in the 30.27 Mb genomic draft of L. orientalis isolate PCM2 with 3367 contigs and 8887 predicted genes. The results from the integrated approach showed the best integrity, coverage, and contig alignment when compared to the genome of L. orientalis isolate LSCM4 collected from the northern province of Thailand. Similar patterns of gene ratios and frequency were observed from the GO biological process annotation. Fifty GO terms were assigned to the assembled genomes, and 23 of these (accounting for 61.6% of the annotated genes) showed higher gene counts and ratios when results from our workflow were compared to those of the LSCM4 isolate. Conclusions: These results indicated that our proposed bioinformatic workflow produced an acceptable-quality genome of L. orientalis strain PCM2 for functional genomic analysis, maximising the usage of the short-read data. This workflow would give extensive information required for identifying strain-specific markers and virulence-associated genes useful for drug and vaccine development before a more exhaustive and expensive investigation.

Heat Shock Proteins as the Druggable Targets in Leishmaniasis: Promises and Perils

Leishmania, the causative agent of leishmaniasis, is an intracellular pathogen that thrives in the insect gut and mammalian macrophages to complete its life cycle. Apart from temperature difference (26 to 37°C), it encounters several harsh conditions, including oxidative stress, inflammatory reactions, and low pH. Heat shock proteins (HSPs) play essential roles in cell survival by strategically reprogramming cellular processes and signaling pathways. HSPs assist cells in multiple functions, including differentiation, adaptation, virulence, and persistence in the host cell. Due to cyclical epidemiological patterns, limited chemotherapeutic options, drug resistance, and the absence of a vaccine, control of leishmaniasis remains a far-fetched dream. The essential roles of HSPs in parasitic differentiation and virulence and increased expression in drug-resistant strains highlight their importance in combating the disease. In this review, we highlighted the diverse physiological importance of HSPs present in Leishmania, emphasizing their significance in disease pathogenesis. Subsequently, we assessed the potential of HSPs as a chemotherapeutic target and underlined the challenges associated with it. Furthermore, we have summarized a few ongoing drug discovery initiatives that need to be explored further to develop clinically successful chemotherapeutic agents in the future.

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

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

Identification and functional characterization of a bacterial homologue of Zeta toxin in Leishmania donovani

Zeta-toxin is a cognate toxin of epsilon antitoxin of prokaryotic Type II toxin-antitoxin system (TA) and play an important role in cell death. An orthologue of bacterial-zeta-toxin (BzT) was identified in Leishmania donovani with similar structural and functional features. Leishmania zeta-toxin (named Ld_ζ1) harboring similar UNAG and ATP-binding pockets showed UNAG kinase and ATP-binding activity. An active Ld_ζ1 was found to express in infective extracellular promastigotes stage of L. donovani and episomal overexpression of an active Ld_ζ1domain-triggered cell death. This study demonstrates the presence of prokaryotic-like-zeta-toxin in eukaryotic parasite Leishmania and its association with cell death. Conceivably, phosphorylated UNAG or analogues, the biochemical mimics of zeta-toxin function mediating cell death can act as a novel anti-leishmanial chemotherapeutics.

Codon choice directs constitutive mRNA levels in trypanosomes

Selective transcription of individual protein coding genes does not occur in trypanosomes and the cellular copy number of each mRNA must be determined post-transcriptionally. Here, we provide evidence that codon choice directs the levels of constitutively expressed mRNAs. First, a novel codon usage metric, the gene expression codon adaptation index (geCAI), was developed that maximised the relationship between codon choice and the measured abundance for a transcriptome. Second, geCAI predictions of mRNA levels were tested using differently coded GFP transgenes and were successful over a 25-fold range, similar to the variation in endogenous mRNAs. Third, translation was necessary for the accelerated mRNA turnover resulting from codon choice. Thus, in trypanosomes, the information determining the levels of most mRNAs resides in the open reading frame and translation is required to access this information.

Transcriptional Regulation of Telomeric Expression Sites and Antigenic Variation in Trypanosomes

INTRODUCTION

Trypanosoma brucei uses antigenic variation to evade the host antibody clearance by periodically changing its Variant Surface Glycoprotein (VSGs) coat. T. brucei encode over 2,500 VSG genes and pseudogenes, however they transcribe only one VSG gene at time from one of the 20 telomeric Expression Sites (ESs). VSGs are transcribed in a monoallelic fashion by RNA polymerase I from an extranucleolar site named ES body. VSG antigenic switching occurs by transcriptional switching between telomeric ESs or by recombination of the VSG gene expressed. VSG expression is developmentally regulated and its transcription is controlled by epigenetic mechanisms and influenced by a telomere position effect.

CONCLUSION

Here, we discuss 1) the molecular basis underlying transcription of telomeric ESs and VSG antigenic switching; 2) the current knowledge of VSG monoallelic expression; 3) the role of inositol phosphate pathway in the regulation of VSG expression and switching; and 4) the developmental regulation of Pol I transcription of procyclin and VSG genes.

Cell cycle stage-specific transcriptional activation of cyclins mediated by HAT2-dependent H4K10 acetylation of promoters in Leishmania donovani

Chromatin modifications affect several processes. In investigating the Leishmania donovani histone acetyltransferase HAT2, using in vitro biochemical assays and HAT2-heterozygous genomic knockout we found the constitutively nuclear HAT2 acetylated histone H4K10 in vitro and in vivo. HAT2 was essential. HAT2-depleted cells displayed growth and cell cycle defects, and poor survival in host cells. Real time PCR and DNA microarray analyses, as well as rescue experiments, revealed that downregulation of cyclins CYC4 and CYC9 were responsible for S phase and G2/M defects of HAT2-depleted cells respectively. Leishmania genes are arranged in unidirectional clusters, and clustered genes are coordinately transcribed as long polycistronic units, typically from divergent strand switch regions (dSSRs) which initiate transcription bidirectionally on opposite strands. In investigating the mechanism by which CYC4 and CYC9 expression levels are reduced in HAT2-depleted cells without other genes in their polycistronic transcription units being coordinately downregulated, we found using reporter assays that CYC4 and CYC9 have their own specific promoters. Chromatin immunoprecipitation assays with H4acetylK10 antibodies and real time PCR analyses of RNA suggested these gene-specific promoters were activated in cell cycle-dependent manner. Nuclear run-on analyses confirmed that CYC4 and CYC9 were transcriptionally activated from their own promoters at specific cell cycle stages. Thus, there are two tiers of gene regulation. Transcription of polycistronic units primarily initiates at dSSRs, and this most likely occurs constitutively. A subset of genes have their own promoters, at least some of which are activated in a cell-cycle dependent manner. This second tier of regulation is more sensitive to H4K10 acetylation levels, resulting in downregulation of expression in HAT2-depleted cells. This report presents the first data pointing to cell cycle-specific activation of promoters in trypanosomatids, thus uncovering new facets of gene regulation in this parasite family.

Nuclear DNA replication and repair in parasites of the genus Leishmania: Exploiting differences to develop innovative therapeutic approaches

Leishmaniasis is a common tropical disease that affects mainly poor people in underdeveloped and developing countries. This largely neglected infection is caused by Leishmania spp, a parasite from the Trypanosomatidae family. This parasitic disease has different clinical manifestations, ranging from localized cutaneous to more harmful visceral forms. The main limitations of the current treatments are their high cost, toxicity, lack of specificity, and long duration. Efforts to improve treatments are necessary to deal with this infectious disease. Many approved drugs to combat diseases as diverse as cancer, bacterial, or viral infections take advantage of specific features of the causing agent or of the disease. Recent evidence indicates that the specific characteristics of the Trypanosomatidae replication and repair machineries could be used as possible targets for the development of new treatments. Here, we review in detail the molecular mechanisms of DNA replication and repair regulation in trypanosomatids of the genus Leishmania and the drugs that could be useful against this disease.

A mitochondrial HSP70 (HSPA9B) is linked to miltefosine resistance and stress response in Leishmania donovani

Background

Protozoan parasites of the genus Leishmania are responsible for leishmaniasis, a neglected tropical disease affecting millions worldwide. Visceral leishmaniasis (VL), caused by Leishmania donovani, is the most severe form of leishmaniasis with high rates of mortality if left untreated. Current treatments include pentavalent antimonials and amphotericin B. However, high toxicity and emergence of resistance hinder the success of these options. Miltefosine (HePC) is the first oral treatment available for leishmaniasis. While treatment with HePC has proven effective, higher tolerance to the drug has been observed, and experimental resistance is easily developed in an in vitro environment. Several studies, including ours, have revealed that HePC resistance has a multi-factorial origin and this work aims to shed light on this complex mechanism.

Methods

2D-DIGE quantitative proteomics comparing the soluble proteomes of sensitive and HePC resistant L. donovani lines identified a protein of interest tentatively involved in drug resistance. To test this link, we employed a gain-of-function approach followed by mutagenesis analysis. Functional studies were complemented with flow cytometry to measure HePC incorporation and cell death.

Results

We identified a mitochondrial HSP70 (HSPA9B) downregulated in HePC-resistant L. donovani promastigotes. The overexpression of HSPA9B in WT lines confers an increased sensitivity to HePC, regardless of whether the expression is ectopic or integrative. Moreover, the increased sensitivity to HePC is specific to the HSPA9B overexpression since dominant negative mutant lines were able to restore HePC susceptibility to WT values. Interestingly, the augmented susceptibility to HePC did not correlate with an increased HePC uptake. Leishmania donovani promastigotes overexpressing HSPA9B were subjected to different environmental stimuli. Our data suggest that HSPA9B is capable of protecting cells from stressful conditions such as low pH and high temperature. This phenotype was further corroborated in axenic amastigotes overexpressing HSPA9B.

Conclusions

The results from this study provide evidence to support the involvement of a mitochondrial HSP70 (HSPA9B) in experimental HePC resistance, a mechanism that is not yet fully understood, and reveal potential fundamental roles of HSPA9B in the biology of Leishmania. Overall, our findings are relevant for current and future antileishmanial chemotherapy strategies.

Electronic supplementary material

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

Chromosomal Translocations in the Parasite Leishmania by a MRE11/RAD50-Independent Microhomology-Mediated End Joining Mechanism

The parasite Leishmania often relies on gene rearrangements to survive stressful environments. However, safeguarding a minimum level of genome integrity is important for cell survival. We hypothesized that maintenance of genomic integrity in Leishmania would imply a leading role of the MRE11 and RAD50 proteins considering their role in DNA repair, chromosomal organization and protection of chromosomes ends in other organisms. Attempts to generate RAD50 null mutants in a wild-type background failed and we provide evidence that this gene is essential. Remarkably, inactivation of RAD50 was possible in a MRE11 null mutant that we had previously generated, providing good evidence that RAD50 may be dispensable in the absence of MRE11. Inactivation of the MRE11 and RAD50 genes led to a decreased frequency of homologous recombination and analysis of the null mutants by whole genome sequencing revealed several chromosomal translocations. Sequencing of the junction between translocated chromosomes highlighted microhomology sequences at the level of breakpoint regions. Sequencing data also showed a decreased coverage at subtelomeric locations in many chromosomes in the MRE11^-/-RAD50^-/- parasites. This study demonstrates an MRE11-independent microhomology-mediated end-joining mechanism and a prominent role for MRE11 and RAD50 in the maintenance of genomic integrity. Moreover, we suggest the possible involvement of RAD50 in subtelomeric regions stability.

The parasite Leishmania relies on gene rearrangements to survive stressful conditions. However, maintaining a minimum level of genomic integrity is crucial for cell survival. Studies in other organisms have provided evidence that the DNA repair proteins MRE11 and RAD50 are involved in chromosomes organization, protection of chromosomes ends and therefore in the maintenance of genomic integrity. In this manuscript, we present the conditional inactivation of the Leishmania infantum RAD50 gene that was only possible in MRE11 deficient cells and suggest the genetic background is crucial for RAD50 inactivation. We demonstrate the occurrence of chromosomal translocations in the MRE11 and RAD50 deficient cells and described a MRE11-independent microhomology-mediated end-joining mechanism at the level of translocation breakpoints. We also suggest a possible involvement of RAD50 in subtelomeric regions stability. Our results highlight that both MRE11 and RAD50 are important for the maintenance of genomic integrity in Leishmania.

Influence of the Microenvironment in the Transcriptome of Leishmania infantum Promastigotes: Sand Fly versus Culture

Zoonotic visceral leishmaniasis is a vector-borne disease caused by Leishmania infantum in the Mediterranean Basin, where domestic dogs and wild canids are the main reservoirs. The promastigote stage replicates and develops within the gut of blood-sucking phlebotomine sand flies. Mature promastigotes are injected in the dermis of the mammalian host and differentiate into the amastigote stage within parasitophorous vacuoles of phagocytic cells. The major vector of L. infantum in Spain is Phlebotomus perniciosus. Promastigotes are routinely axenized and cultured to mimic in vitro the conditions inside the insect gut, which allows for most molecular, cellular, immunological and therapeutical studies otherwise inviable. Culture passages are known to decrease infectivity, which is restored by passage through laboratory animals. The most appropriate source of promastigotes is the gut of the vector host but isolation of the parasite is technically challenging. In fact, this option is not viable unless small samples are sufficient for downstream applications like promastigote cultures and nucleic acid amplification. In this study, in vitro infectivity and differential gene expression have been studied in cultured promastigotes at the stationary phase and in promastigotes isolated from the stomodeal valve of the sand fly P. perniciosus. About 20 ng RNA per sample could be isolated. Each sample contained L. infantum promastigotes from 20 sand flies. RNA was successfully amplified and processed for shotgun genome microarray hybridization analysis. Most differentially regulated genes are involved in regulation of gene expression, intracellular signaling, amino acid metabolism and biosynthesis of surface molecules. Interestingly, meta-analysis by hierarchical clustering supports that up-regulation of 22.4% of the differentially regulated genes is specifically enhanced by the microenvironment (i.e. sand fly gut or culture). The correlation between cultured and naturally developed promastigotes is strong but not very high (Pearson coefficient R² = 0.727). Therefore, the influence of promastigote culturing should be evaluated case-by-case in experimentation.

The protozoan parasite Leishmania infantum causes visceral leishmaniasis in humans and is responsible for a recent outbreak reported in central Spain. Domestic dogs and wild canids are the main reservoirs. The life cycle of the parasite involves two stages and two hosts. The motile promastigote stage differentiates within the gut of the sand fly vector host and develops into non-motile amastigotes within phagocytes of the mammalian host. Promastigotes are routinely cultured in liquid media because it is assumed that they mimic the conditions within the gut of the insect. Therefore, the culture model is used in most studies about the biology of the parasite, pathogenesis and development of vaccines and new compounds for treatment. Isolating promastigotes from the natural microenvironment (i.e. the vector host) is desirable but technically challenging. We were able to perform a high-throughput analysis of gene expression thanks to mRNA amplification. The over-expressed genes detected may influence life cycle progression depending on the promastigote microenvironment (i.e. culture or vector host). Upcoming studies based on these results may reveal new therapeutic targets or vaccine candidates. Our results suggest that evaluating the influence of cultures in experimentation is convenient.

Serum Removal from Culture Induces Growth Arrest, Ploidy Alteration, Decrease in Infectivity and Differential Expression of Crucial Genes in Leishmania infantum Promastigotes

Leishmania infantum is one of the species responsible for visceral leishmaniasis. This species is distributed basically in the Mediterranean basin. A recent outbreak in humans has been reported in Spain. Axenic cultures are performed for most procedures with Leishmania spp. promastigotes. This model is stable and reproducible and mimics the conditions of the gut of the sand fly host, which is the natural environment of promastigote development. Culture media are undefined because they contain mammalian serum, which is a rich source of complex lipids and proteins. Serum deprivation slows down the growth kinetics and therefore, yield in biomass. In fact, we have confirmed that the growth rate decreases, as well as infectivity. Ploidy is also affected. Regarding the transcriptome, a high-throughput approach has revealed a low differential expression rate but important differentially regulated genes. The most remarkable profiles are: up-regulation of the GINS Psf3, the fatty acyl-CoA synthase (FAS1), the glyoxylase I (GLO1), the hydrophilic surface protein B (HASPB), the methylmalonyl-CoA epimerase (MMCE) and an amastin gene; and down-regulation of the gPEPCK and the arginase. Implications for metabolic adaptations, differentiation and infectivity are discussed herein.

The Dynamics of Lateral Gene Transfer in Genus Leishmania - A Route for Adaptation and Species Diversification

Background

The genome of Leishmania major harbours a comparably high proportion of genes of prokaryote origin, acquired by lateral gene transfer (LGT). Some of these are present in closely related trypanosomatids, while some are detected in Leishmania only. We have evaluated the impact and destiny of LGT in genus Leishmania.

Methodology/Principal Findings

To study the dynamics and fate of LGTs we have performed phylogenetic, as well as nucleotide and amino acid composition analyses within orthologous groups of LGTs detected in Leishmania. A set of universal trypanosomatid LGTs was added as a reference group. Both groups of LGTs have, to some extent, ameliorated to resemble the recipient genomes. However, while virtually all of the universal trypanosomatid LGTs are distributed and conserved in the entire genus Leishmania, the LGTs uniquely present in genus Leishmania are more prone to gene loss and display faster rates of evolution. Furthermore, a PCR based approach has been employed to ascertain the presence of a set of twenty LGTs uniquely present in genus Leishmania, and three universal trypanosomatid LGTs, in ten additional strains of Leishmania. Evolutionary rates and predicted expression levels of these LGTs have also been estimated. Ten of the twenty LGTs are distributed and conserved in all species investigated, while the remainder have been subjected to modifications, or undergone pseudogenization, degradation or loss in one or more species.

Conclusions/Significance

LGTs unique to the genus Leishmania have been acquired after the divergence of Leishmania from the other trypanosomatids, and are evolving faster than their recipient genomes. This implies that LGT in genus Leishmania is a continuous and dynamic process contributing to species differentiation and speciation. This study also highlights the importance of carefully evaluating these dynamic genes, e.g. as LGTs have been suggested as potential drug targets.

Leishmania parasites cause leishmaniasis, a neglected tropical disease, estimated to threaten 350 million people in 88 countries worldwide according to the WHO. The genome of Leishmania major harbours a number of genes, which have been proposed as acquired by lateral gene transfer (LGT) from a broad variety of prokaryote donors. Such genes may prove beneficial for the parasites, e.g. by promoting survival of the parasite in new environments. We have studied orthologs to LGTs previously detected uniquely in L. major as well as LGTs shared also by other trypanosomatids. The universal trypanosomatid LGTs are more conserved within genus Leishmania, as compared to LGTs that are exclusively present in genus Leishmania. One possible explanation to this observation is that these have resided in their host genomes for a longer time period. This explanation strengthens the hypothesis that the LGTs unique to genus Leishmaina were acquired after the divergence from the trypanosomes, rather than before the divergence, and then subsequently lost from the trypanosome lineage. An in-depth analysis of a subset of the LGTs, which are present only in genus Leishmania showed that LGT within genus Leishmania is a dynamic process. LGTs, providing beneficial capabilities to the parasite, are demonstrated to become conserved throughout generic diversification, hence contributing to species differentiation, while LGTs of limited use are degraded and lost.

Conserved Curvature of RNA Polymerase I Core Promoter Beyond rRNA Genes: The Case of the Tritryps

In trypanosomatids, the RNA polymerase I (RNAPI)-dependent promoters controlling the ribosomal RNA (rRNA) genes have been well identified. Although the RNAPI transcription machinery recognizes the DNA conformation instead of the DNA sequence of promoters, no conformational study has been reported for these promoters. Here we present the in silico analysis of the intrinsic DNA curvature of the rRNA gene core promoters in Trypanosoma brucei, Trypanosoma cruzi, and Leishmania major. We found that, in spite of the absence of sequence conservation, these promoters hold conformational properties similar to other eukaryotic rRNA promoters. Our results also indicated that the intrinsic DNA curvature pattern is conserved within the Leishmania genus and also among strains of T. cruzi and T. brucei. Furthermore, we analyzed the impact of point mutations on the intrinsic curvature and their impact on the promoter activity. Furthermore, we found that the core promoters of protein-coding genes transcribed by RNAPI in T. brucei show the same conserved conformational characteristics. Overall, our results indicate that DNA intrinsic curvature of the rRNA gene core promoters is conserved in these ancient eukaryotes and such conserved curvature might be a requirement of RNAPI machinery for transcription of not only rRNA genes but also protein-coding genes.

Ribosome profiling reveals translation control as a key mechanism generating differential gene expression in Trypanosoma cruzi

Background

Due to the absence of transcription initiation regulation of protein coding genes transcribed by RNA polymerase II, posttranscriptional regulation is responsible for the majority of gene expression changes in trypanosomatids. Therefore, cataloging the abundance of mRNAs (transcriptome) and the level of their translation (translatome) is a key step to understand control of gene expression in these organisms.

Results

Here we assess the extent of regulation of the transcriptome and the translatome in the Chagas disease causing agent, Trypanosoma cruzi, in both the non-infective (epimastigote) and infective (metacyclic trypomastigote) insect’s life stages using RNA-seq and ribosome profiling. The observed steady state transcript levels support constitutive transcription and maturation implying the existence of distinctive posttranscriptional regulatory mechanisms controlling gene expression levels at those parasite stages. Meanwhile, the downregulation of a large proportion of the translatome indicates a key role of translation control in differentiation into the infective form. The previously described proteomic data correlate better with the translatomes than with the transcriptomes and translational efficiency analysis shows a wide dynamic range, reinforcing the importance of translatability as a regulatory step. Translation efficiencies for protein families like ribosomal components are diminished while translation of the transialidase virulence factors is upregulated in the quiescent infective metacyclic trypomastigote stage.

Conclusions

A large subset of genes is modulated at the translation level in two different stages of Trypanosoma cruzi life cycle. Translation upregulation of virulence factors and downregulation of ribosomal proteins indicates different degrees of control operating to prepare the parasite for an infective life form. Taking together our results show that translational regulation, in addition to regulation of steady state level of mRNA, is a major factor playing a role during the parasite differentiation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1563-8) contains supplementary material, which is available to authorized users.

RNA interference screen reveals a high proportion of mitochondrial proteins essential for correct cell cycle progress in Trypanosoma brucei

Background

Trypanosomatid parasites possess a single mitochondrion which is classically involved in the energetic metabolism of the cell, but also, in a much more original way, through its single and complex DNA (termed kinetoplast), in the correct progress of cell division. In order to identify proteins potentially involved in the cell cycle, we performed RNAi knockdowns of 101 genes encoding mitochondrial proteins using procyclic cells of Trypanosoma brucei.

Results

A major cell growth reduction was observed in 10 cases and a moderate reduction in 29 other cases. These data are overall in agreement with those previously obtained by a case-by-case approach performed on chromosome 1 genes, and quantitatively with those obtained by “high-throughput phenotyping using parallel sequencing of RNA interference targets” (RIT-seq). Nevertheless, a detailed analysis revealed many qualitative discrepancies with the RIT-seq-based approach. Moreover, for 37 out of 39 mutants for which a moderate or severe growth defect was observed here, we noted abnormalities in the cell cycle progress, leading to increased proportions of abnormal cell cycle stages, such as cells containing more than 2 kinetoplasts (K) and/or more than 2 nuclei (N), and modified proportions of the normal phenotypes (1N1K, 1N2K and 2N2K).

Conclusions

These data, together with the observation of other abnormal phenotypes, show that all the corresponding mitochondrial proteins are involved, directly or indirectly, in the correct progress or, less likely, in the regulation, of the cell cycle in T. brucei. They also show how post-genomics analyses performed on a case-by-case basis may yield discrepancies with global approaches.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1505-5) contains supplementary material, which is available to authorized users.

Genomic analysis of sequence-dependent DNA curvature in Leishmania

Leishmania major is a flagellated protozoan parasite of medical importance. Like other members of the Trypanosomatidae family, it possesses unique mechanisms of gene expression such as constitutive polycistronic transcription of directional gene clusters, gene amplification, mRNA trans-splicing, and extensive editing of mitochondrial transcripts. The molecular signals underlying most of these processes remain under investigation. In order to investigate the role of DNA secondary structure signals in gene expression, we carried out a genome-wide in silico analysis of the intrinsic DNA curvature. The L. major genome revealed a lower frequency of high intrinsic curvature regions as well as inter- and intra- chromosomal distribution heterogeneity, when compared to prokaryotic and eukaryotic organisms. Using a novel method aimed at detecting region-integrated intrinsic curvature (RIIC), high DNA curvature was found to be associated with regions implicated in transcription initiation. Those include divergent strand-switch regions between directional gene clusters and regions linked to markers of active transcription initiation such as acetylated H3 histone, TRF4 and SNAP50. These findings suggest a role for DNA curvature in transcription initiation in Leishmania supporting the relevance of DNA secondary structures signals.

Footprints of a trypanosomatid RNA world: pre-small subunit rRNA processing by spliced leader addition trans-splicing

The addition of a capped mini-exon [spliced leader (SL)] through trans-splicing is essential for the maturation of RNA polymerase (pol) II-transcribed polycistronic pre-mRNAs in all members of the Trypanosomatidae family. This process is an inter-molecular splicing reaction that follows the same basic rules of cis-splicing reactions. In this study, we demonstrated that mini-exons were added to precursor ribosomal RNA (pre-rRNA) are transcribed by RNA pol I, including the 5' external transcribed spacer (ETS) region. Additionally, we detected the SL-5'ETS molecule using three distinct methods and located the acceptor site between two known 5'ETS rRNA processing sites (A' and A1) in four different trypanosomatids. Moreover, we detected a polyadenylated 5'ETS upstream of the trans-splicing acceptor site, which also occurs in pre-mRNA trans-splicing. After treatment with an indirect trans-splicing inhibitor (sinefungin), we observed SL-5'ETS decay. However, treatment with 5-fluorouracil (a precursor of RNA synthesis that inhibits the degradation of pre-rRNA) led to the accumulation of SL-5'ETS, suggesting that the molecule may play a role in rRNA degradation. The detection of trans-splicing in these molecules may indicate broad RNA-joining properties, regardless of the polymerase used for transcription.

Glucosylated hydroxymethyluracil, DNA base J, prevents transcriptional readthrough in Leishmania

Some Ts in nuclear DNA of trypanosomes and Leishmania are hydroxylated and glucosylated to yield base J (β-D-glucosyl-hydroxymethyluracil). In Leishmania, about 99% of J is located in telomeric repeats. We show here that most of the remaining J is located at chromosome-internal RNA polymerase II termination sites. This internal J and telomeric J can be reduced by a knockout of J-binding protein 2 (JBP2), an enzyme involved in the first step of J biosynthesis. J levels are further reduced by growing Leishmania JBP2 knockout cells in BrdU-containing medium, resulting in cell death. The loss of internal J in JBP2 knockout cells is accompanied by massive readthrough at RNA polymerase II termination sites. The readthrough varies between transcription units but may extend over 100 kb. We conclude that J is required for proper transcription termination and infer that the absence of internal J kills Leishmania by massive readthrough of transcriptional stops.

Trypanosomatid comparative genomics: Contributions to the study of parasite biology and different parasitic diseases

In 2005, draft sequences of the genomes of Trypanosoma brucei, Trypanosoma cruzi and Leishmania major, also known as the Tri-Tryp genomes, were published. These protozoan parasites are the causative agents of three distinct insect-borne diseases, namely sleeping sickness, Chagas disease and leishmaniasis, all with a worldwide distribution. Despite the large estimated evolutionary distance among them, a conserved core of ~6,200 trypanosomatid genes was found among the Tri-Tryp genomes. Extensive analysis of these genomic sequences has greatly increased our understanding of the biology of these parasites and their host-parasite interactions. In this article, we review the recent advances in the comparative genomics of these three species. This analysis also includes data on additional sequences derived from other trypanosmatid species, as well as recent data on gene expression and functional genomics. In addition to facilitating the identification of key parasite molecules that may provide a better understanding of these complex diseases, genome studies offer a rich source of new information that can be used to define potential new drug targets and vaccine candidates for controlling these parasitic infections.

Genomic organization of leishmania species

Leishmania is a protozoan parasite belonging to the family Trypanosomatidae, which is found among 88 different countries. The parasite lives as an amastigote in vertebrate macrophages and as a promastigote in the digestive tract of sand fly. It can be cultured in the laboratory using appropriate culture media. Although the sexual cycle of Leishmania has not been observed during the promastigote and amastigote stages, it has been reported by some researchers. Leishmania has eukaryotic cell organization. Cell culture is convenient and cost effective, and because posttranslational modifications are common processes in the cultured cells, the cells are used as hosts for preparing eukaryotic recombinant proteins for research. Several transcripts of rDNA in the Leishmania genome are suitable regions for conducting gene transfer. Old World Leishmania spp. has 36 chromosomes, while New World Leishmania spp. has 34 or 35 chromosomes. The genomic organization and parasitic characteristics have been investigated. Leishmania spp. has a unique genomic organization among eukaryotes; the genes do not have introns, and the chromosomes are smaller with larger numbers of genes confined to a smaller space within the nucleus. Leishmania spp. genes are organized on one or both DNA strands and are transcribed as polycistronic (prokaryotic-like) transcripts from undefined promoters. Regulation of gene expression in the members of Trypanosomatidae differs from that in other eukaryotes. The trans-splicing phenomenon is a necessary step for mRNA processing in lower eukaryotes and is observed in Leishmania spp. Another particular feature of RNA editing in Leishmania spp. is that mitochondrial genes encoding respiratory enzymes are edited and transcribed. This review will discuss the chromosomal and mitochondrial (kinetoplast) genomes of Leishmania spp. as well as the phenomenon of RNA editing in the kinetoplast genome.

Multiple levels of gene regulation mediate differentiation of the intracellular pathogen Leishmania

For many years, mRNA abundance has been used as the surrogate measure of gene expression in biological systems. However, recent genome-scale analyses in both bacteria and eukaryotes have revealed that mRNA levels correlate with steady-state protein abundance for only 50-70% of genes, indicating that translation and post-translation processes also play important roles in determining gene expression. What is not yet clear is whether dynamic processes such as cell cycle progression, differentiation, or response to environmental changes change the relationship between mRNA and protein abundance. Here, we describe a systems approach to interrogate promastigote-to-amastigote differentiation in the obligatory intracellular parasitic protozoan Leishmania donovani. Our results indicate that regulation of mRNA levels plays a major role early in the differentiation process, while translation and post-translational regulation are more important in the latter part. In addition, it appears that the differentiation signal causes a transient global increase in the rate of protein synthesis, which is subsequently down-regulated by phosphorylation of α-subunit of translation initiation factor 2. Thus, Leishmania dynamically changes the relationship between mRNA and protein abundance as it adapts to new environmental circumstances. It is likely that similar mechanisms play a more important role than previously recognized in regulation of gene expression in other organisms.

Regulation of gene expression in protozoa parasites

Infections with protozoa parasites are associated with high burdens of morbidity and mortality across the developing world. Despite extensive efforts to control the transmission of these parasites, the spread of populations resistant to drugs and the lack of effective vaccines against them contribute to their persistence as major public health problems. Parasites should perform a strict control on the expression of genes involved in their pathogenicity, differentiation, immune evasion, or drug resistance, and the comprehension of the mechanisms implicated in that control could help to develop novel therapeutic strategies. However, until now these mechanisms are poorly understood in protozoa. Recent investigations into gene expression in protozoa parasites suggest that they possess many of the canonical machineries employed by higher eukaryotes for the control of gene expression at transcriptional, posttranscriptional, and epigenetic levels, but they also contain exclusive mechanisms. Here, we review the current understanding about the regulation of gene expression in Plasmodium sp., Trypanosomatids, Entamoeba histolytica and Trichomonas vaginalis.

Gene expression in trypanosomatid parasites

The parasites Leishmania spp., Trypanosoma brucei, and Trypanosoma cruzi are the trypanosomatid protozoa that cause the deadly human diseases leishmaniasis, African sleeping sickness, and Chagas disease, respectively. These organisms possess unique mechanisms for gene expression such as constitutive polycistronic transcription of protein-coding genes and trans-splicing. Little is known about either the DNA sequences or the proteins that are involved in the initiation and termination of transcription in trypanosomatids. In silico analyses of the genome databases of these parasites led to the identification of a small number of proteins involved in gene expression. However, functional studies have revealed that trypanosomatids have more general transcription factors than originally estimated. Many posttranslational histone modifications, histone variants, and chromatin modifying enzymes have been identified in trypanosomatids, and recent genome-wide studies showed that epigenetic regulation might play a very important role in gene expression in this group of parasites. Here, we review and comment on the most recent findings related to transcription initiation and termination in trypanosomatid protozoa.

Molecular analysis of surface glycoprotein multigene family TrGP expressed on the plasma membrane of Trypanosoma rangeli epimastigotes forms

Trypanosoma rangeli, a non-pathogenic hemoflagelate that in Central and South America infects humans, shares with Trypanosoma cruzi reservoirs and triatomine vectors, as well as geographical distribution. Recently, we have described in T. rangeli a truncated gene copy belonging to the group II of the trans-sialidase superfamily (TrGP). This superfamily, collectively known in T. cruzi as gp85/TS, includes members that are involved in host cell invasion and infectivity. To confirm the presence of this superfamily in the genome of T. rangeli and obtain a better knowledge of its characteristics, we designed a PCR and RT-PCR cloning strategy to allow sequence analysis of both genomic and transcribed copies. We identified two full-length copies of TrGP, some pseudogenes, and N- and C-terminal sequences of several genes. We also analyzed the expression and cellular localization of these proteins in epimastigote forms of a Venezuelan T. rangeli isolate using polyclonal antibodies made against a recombinant peptide from the N-terminal region of a TrGP member. We confirmed that TrGP is a multigenic family that shares many features with T. cruzi gp85/TS, including the telomeric location of some of its members, and by immunofluorescence analysis that its location is at the surface of the parasite.

Organization and evolution of two SIDER retroposon subfamilies and their impact on the Leishmania genome

Background

We have recently identified two large families of extinct transposable elements termed Short Interspersed DEgenerated Retroposons (SIDERs) in the parasitic protozoan Leishmania major. The characterization of SIDER elements was limited to the SIDER2 subfamily, although members of both subfamilies have been shown to play a role in the regulation of gene expression at the post-transcriptional level. Apparent functional domestication of SIDERs prompted further investigation of their characterization, dissemination and evolution throughout the Leishmania genus, with particular attention to the disregarded SIDER1 subfamily.

Results

Using optimized statistical profiles of both SIDER1 and SIDER2 subgroups, we report the first automated and highly sensitive annotation of SIDERs in the genomes of L. infantum, L. braziliensis and L. major. SIDER annotations were combined to in-silico mRNA extremity predictions to generate a detailed distribution map of the repeat family, hence uncovering an enrichment of antisense-oriented SIDER repeats between the polyadenylation and trans-splicing sites of intergenic regions, in contrast to the exclusive sense orientation of SIDER elements within 3'UTRs. Our data indicate that SIDER elements are quite uniformly dispersed throughout all three genomes and that their distribution is generally syntenic. However, only 47.4% of orthologous genes harbor a SIDER element in all three species. There is evidence for species-specific enrichment of SIDERs and for their preferential association, especially for SIDER2s, with different metabolic functions. Investigation of the sequence attributes and evolutionary relationship of SIDERs to other trypanosomatid retroposons reveals that SIDER1 is a truncated version of extinct autonomous ingi-like retroposons (DIREs), which were functional in the ancestral Leishmania genome.

Conclusion

A detailed characterization of the sequence traits for both SIDER subfamilies unveils major differences. The SIDER1 subfamily is more heterogeneous and shows an evolutionary link with vestigial DIRE retroposons as previously observed for the ingi/RIME and L1Tc/NARTc couples identified in the T. brucei and T. cruzi genomes, whereas no identified DIREs are related to SIDER2 sequences. Although SIDER1s and SIDER2s display equivalent genomic distribution globally, the varying degrees of sequence conservation, preferential genomic disposition, and differential association to orthologous genes allude to an intricate web of SIDER assimilation in these parasitic organisms.

Dicistronic tRNA-5S rRNA genes in Yarrowia lipolytica: an alternative TFIIIA-independent way for expression of 5S rRNA genes

In eukaryotes, genes transcribed by RNA polymerase III (Pol III) carry their own internal promoters and as such, are transcribed as individual units. Indeed, a very few cases of dicistronic Pol III genes are yet known. In contrast to other hemiascomycetes, 5S rRNA genes of Yarrowia lipolytica are not embedded into the tandemly repeated rDNA units, but appear scattered throughout the genome. We report here an unprecedented genomic organization: 48 over the 108 copies of the 5S rRNA genes are located 3' of tRNA genes. We show that these peculiar tRNA-5S rRNA dicistronic genes are expressed in vitro and in vivo as Pol III transcriptional fusions without the need of the 5S rRNA gene-specific factor TFIIIA, the deletion of which displays a viable phenotype. We also report the existence of a novel putative non-coding Pol III RNA of unknown function about 70 nucleotide-long (RUF70), the 13 genes of which are devoid of internal Pol III promoters and located 3' of the 13 copies of the tDNA-Trp (CCA). All genes embedded in the various dicistronic genes, fused 5S rRNA genes, RUF70 genes and their leader tRNA genes appear to be efficiently transcribed and their products correctly processed in vivo.

From stop to start: tandem gene arrangement, copy number and trans-splicing sites in the dinoflagellate Amphidinium carterae

Dinoflagellate genomes present unique challenges including large size, modified DNA bases, lack of nucleosomes, and condensed chromosomes. EST sequencing has shown that many genes are found as many slightly different variants implying that many copies are present in the genome. As a preliminary survey of the genome our goal was to obtain genomic sequences for 47 genes from the dinoflagellate Amphidinium carterae. A PCR approach was used to avoid problems with large insert libraries. One primer set was oriented inward to amplify the genomic complement of the cDNA and a second primer set would amplify outward between tandem repeats of the same gene. Each gene was also tested for a spliced leader using cDNA as template. Almost all (14/15) of the highly expressed genes (i.e. those with high representation in the cDNA pool) were shown to be in tandem arrays with short intergenic spacers, and most were trans-spliced. Only two moderately expressed genes were found in tandem arrays. A polyadenylation signal was found in genomic copies containing the sequence AAAAG/C at the exact polyadenylation site and was conserved between species. Four genes were found to have a high intron density (>5 introns) while most either lacked introns, or had only one to three. Actin was selected for deeper sequencing of both genomic and cDNA copies. Two clusters of actin copies were found, separated from each other by many non-coding features such as intron size and sequence. One intron-rich gene was selected for genomic walking using inverse PCR, and was not shown to be in a tandem repeat. The first glimpse of dinoflagellate genome indicates two general categories of genes in dinoflagellates, a highly expressed tandem repeat class and an intron rich less expressed class. This combination of features appears to be unique among eukaryotes.

Genome diversity in microbial eukaryotes

The genomic peculiarities among microbial eukaryotes challenge the conventional wisdom of genome evolution. Currently, many studies and textbooks explore principles of genome evolution from a limited number of eukaryotic lineages, focusing often on only a few representative species of plants, animals and fungi. Increasing emphasis on studies of genomes in microbial eukaryotes has and will continue to uncover features that are either not present in the representative species (e.g. hypervariable karyotypes or highly fragmented mitochondrial genomes) or are exaggerated in microbial groups (e.g. chromosomal processing between germline and somatic nuclei). Data for microbial eukaryotes have emerged from recent genome sequencing projects, enabling comparisons of the genomes from diverse lineages across the eukaryotic phylogenetic tree. Some of these features, including amplified rDNAs, subtelomeric rDNAs and reduced genomes, appear to have evolved multiple times within eukaryotes, whereas other features, such as absolute strand polarity, are found only within single lineages.

Leishmania: conserved evolution--diverse diseases

The landmark completion of the Leishmania major genome sequence and the recent publication of the L. infantum and L. braziliensis genomes revealed the surprising result that, although separated by 15-50 million years of evolution, the Leishmania genomes are highly conserved and have less than 1% species-specific genes. Yet, these three species of Leishmania cause distinctive and diverse diseases in humans. Here, we discuss these findings together with recent microarray and proteomics studies and highlight their importance in understanding Leishmania disease phenotypes.

Global gene expression in Leishmania

The completion of the genomic sequences of many protozoan pathogens of humans, including species of Leishmania, Trypanosoma and Plasmodium, provide new approaches to study the pattern of gene expression during differentiation and development. Leishmania are a major public health risk in many countries and cause a wide spectrum of clinical disease referred to as leishmaniasis. The Leishmania life cycle consists of two morphologically distinct stages: intracellular amastigotes that reside in the phagolysosome of mammalian macrophages, and extracellular promastigotes that reside within the gut of the sandfly vector. DNA microarray analysis is a powerful method to study global gene expression in terms of quantitation of mRNA levels. This review discusses the application of DNA microarray technology to study the pattern of global gene expression of Leishmania promastigote and amastigote life stages. Results from several studies show that, overall, there is a surprisingly low level of differentially expressed genes, ranging from 0.2% to 5% of total genes, between the amastigote and promastigote life stages. Thus, the Leishmania genome can be considered to be constitutively expressed with a limited number of genes showing stage-specific expression. Comparative genomic analyses of gene expression levels between Leishmania major and Leishmania mexicana show that the majority of differentially expressed genes between amastigotes and promastigotes are species specific with relatively few differentially expressed genes in common between these two Leishmania species. Quantitative proteomic analysis of Leishmania relative protein expression shows there is a weak correlation to gene expression. Therefore, Leishmania protein expression levels are likely regulated at the level of translation or by post transcriptional mechanisms, and differential protein modifications may be more important in development than the regulation of gene expression.

Genomic and proteomic expression analysis of Leishmania promastigote and amastigote life stages: The Leishmania genome is constitutively expressed

Leishmania are protozoan parasites that cause a wide spectrum of clinical diseases in humans and are a major public health risk in several countries. Leishmania life cycle consists of an extracellular flagellated promastigote stage within the midgut of a sandfly vector, and a morphological distinct intracellular amastigote stage within macrophages of a mammalian host. This study reports the use of DNA oligonucleotide genome microarrays representing 8160 genes to analyze the mRNA expression profiles of L. major promastigotes and lesion derived amastigotes. Over 94% of the genes were expressed in both life stages. Advanced statistical analysis identified a surprisingly low degree of differential mRNA expression: 1.4% of the total genes in amastigotes and 1.5% in promastigotes. These microarray results demonstrate that the L. major genome is essentially constitutively expressed in both life stages and suggest that Leishmania is constitutively adapted for survival and replication in either the sandfly vector or macrophage host utilizing an appropriate set of genes for each vastly different environment. Quantitative proteomics, using the isotope coded affinity tag (ICAT) technology and mass spectrometry, was used to identify L. infantum promastigote and axenic amastigote differentially expressed proteins. Of the 91 distinct proteins identified, 8% were differentially expressed in the amastigote stage, 20% were differentially expressed in the promastigote stage, and the remaining 72% were considered constitutively expressed. The differential expression was validated by the identification of previously reported stage specific proteins and identified several amastigote and promastigote novel stage specific proteins.

Compared genomics of the strand switch region of Leishmania chromosome 1 reveal a novel genus-specific gene and conserved structural features and sequence motifs

BACKGROUND

Trypanosomatids exhibit a unique gene organization into large directional gene clusters (DGCs) in opposite directions. The transcription "strand switch region" (SSR) separating the two large DGCs that constitute chromosome 1 of Leishmania major has been the subject of several studies and speculations. Thus, it has been suspected of being the single replication origin of the chromosome, the transcription initiation site for both DGCs or even a centromere. Here, we have used an inter-species compared genomics approach on this locus in order to try to identify conserved features or motifs indicative of a putative function.

RESULTS

We isolated, and compared the structure and nucleotide sequence of, this SSR in 15 widely divergent species of Leishmania and Sauroleishmania. As regards its intrachromosomal position, size and AT content, the general structure of this SSR appears extremely stable among species, which is another demonstration of the remarkable structural stability of these genomes at the evolutionary level. Sequence alignments showed several interesting features. Overall, only 30% of nucleotide positions were conserved in the SSR among the 15 species, versus 74% and 62% in the 5' parts of the adjacent XPP and PAXP genes, respectively. However, nucleotide divergences were not distributed homogeneously along this sequence. Thus, a central fragment of approximately 440 bp exhibited 54% of identity among the 15 species. This fragment actually represents a new Leishmania-specific CDS of unknown function which had been overlooked since the annotation of this chromosome. The encoded protein comprises two trans-membrane domains and is classified in the "structural protein" GO category. We cloned this novel gene and expressed it as a recombinant green fluorescent protein-fused version, which showed its localisation to the endoplasmic reticulum. The whole of these data shorten the actual SSR to an 887-bp segment as compared with the original 1.6 kb. In the rest of the SSR, the percentage of identity was much lower, around 22%. Interestingly, the 72-bp fragment where the putatively single transcription initiation site of chromosome 1 was identified is located in a low-conservation portion of the SSR and is itself highly polymorphic amongst species. Nevertheless, it is highly C-rich and presents a unique poly(C) tract in the same position in all species.

CONCLUSION

This inter-specific comparative study, the first of its kind, (a) allowed to reveal a novel genus-specific gene and (b) identified a conserved poly(C) tract in the otherwise highly polymorphic region containing the putative transcription initiation site. This allows hypothesising an intervention of poly(C)-binding proteins known elsewhere to be involved in transcriptional control.

Analysis of the Leishmania donovani transcriptome reveals an ordered progression of transient and permanent changes in gene expression during differentiation

Leishmania donovani is an intracellular protozoan parasite that causes kala-azar in humans. During infection the extracellular insect forms (promastigotes) undergo rapid differentiation to intracellular amastigotes that proliferates in phagolysosomes of mammalian macrophages. We used microarray-based expression profiling to investigate the time-course of changes in RNA abundance during promastigote-to-amastigote differentiation in a host-free system that mimics this process. These studies revealed that several hundred genes underwent an ordered progression of transient or permanent up- and down-regulation during differentiation. Genes that were permanently up-regulated in amastigotes were enriched for transporters and surface proteins, but under-represented in genes involved in protein and other metabolism. Most of these changes occurred late in the differentiation process, when morphological differentiation was essentially complete. Down-regulated genes were over-represented in those involved in cell motility, growth and/or maintenance, and these changes generally occurred earlier in the process. Genes that were transiently up- or down-regulated during differentiation included those encoding heat shock proteins, ubiquitin hydrolases, RNA binding proteins, protein kinases, a protein phosphatase, and a histone deacetylase. These results suggest that changes in mRNA abundance may be important in signal transduction, as well as protein and mRNA turnover, during differentiation. In addition to these mRNA changes, other transcripts including one or more rRNAs and snoRNAs, and non-coding RNAs from several telomeres, also showed substantial changes in abundance during the differentiation process. This paper provides the first genome-scale quantitative analysis of gene expression during the transition from promastigotes to amastigotes and demonstrates the utility of the host-free differentiation system.

AUG-proximal nucleotides regulate protein synthesis in Leishmania tropica

Gene expression in the Leishmania is controlled post-transcriptionally, and is likely to be impacted by both 5' and 3' untranslated regions (UTRs). We have investigated the effects of trinucleotides in the AUG-proximal region (APR) (i.e. positions -3 to -1 upstream of an AUG) on two reporter genes in the context of an endogenous intergenic region of Leishmania tropica. The effects of APRs on protein expression were determined in stable transfectants in vivo. Three APRs, namely, C(-3)C(-2)C(-1), ACC and GCC, yielded robust translation, whereas GTA produced low amounts of proteins. A purine at -3 of an APR was not crucial for efficient translation. Steady-state level of reporter mRNA did not correlate directly with the amount of protein detected. Polysome analysis revealed that APRs modulate translation, at least in part, by influencing mRNA association with ribosomes. An analysis of genomic UTRs in L. major showed that (i) the consensus APR is N(-3)N(-2)C(-1) (where N = any nucleotide), and (ii) the most frequently used APRs include ACA, ACC, ATC, GCC, GCG, GTC and CAC, some of which were translation enhancers in our experimental studies.

Evolutionary Relationships and Protein Domain Architecture in an Expanded Calpain Superfamily in Kinetoplastid Parasites

Employing whole-genome analysis we have characterized a large family of genes coding for calpain-related proteins in three kinetoplastid parasites. We have defined a total of 18 calpain-like sequences in Trypanosoma brucei, 27 in Leishmania major, and 24 in Trypanosoma cruzi. Sequence characterization revealed a well-conserved protease domain in most proteins, although residues critical for catalytic activity were frequently altered. Many of the proteins contain a novel N-terminal sequence motif unique to kinetoplastids. Furthermore, 24 of the sequences contain N-terminal fatty acid acylation motifs indicating association of these proteins with intracellular membranes. This extended family of proteins also includes a group of sequences that completely lack a protease domain but is specifically related to other kinetoplastid calpain-related proteins by a highly conserved N-terminal domain and by genomic organization. All sequences lack the C-terminal calmodulin-related calcium-binding domain typical of most mammalian calpains. Our analysis emphasizes the highly modular structure of calpains and calpain-like proteins, suggesting that they are involved in diverse cellular functions. The discovery of this surprisingly large family of calpain-like proteins in lower eukaryotes that combines novel and conserved sequence modules contributes to our understanding of the evolution of this abundant protein family.

Parasite genomes

The availability of genome sequences and the associated transcriptome and proteome mapping projects has revolutionised research in the field of parasitology. As more parasite species are sequenced, comparative and phylogenetic comparisons are improving the quality of gene prediction and annotation. Genome sequences of parasites are also providing important data sets for understanding parasite biology and identifying new vaccine candidates and drug targets. We review some of the preliminary conclusions from examination of parasite genome sequences and discuss some of the bioinformatics approaches taken in this analysis.

A computational investigation of kinetoplastid trans-splicing

A novel computational approach is presented and applied to predicting trans-splicing sites in 2 chromosomes of Leishmania major.

Trans-splicing is an unusual process in which two separate RNA strands are spliced together to yield a mature mRNA. We present a novel computational approach which has an overall accuracy of 82% and can predict 92% of known trans-splicing sites. We have applied our method to chromosomes 1 and 3 of Leishmania major, with high-confidence predictions for 85% and 88% of annotated genes respectively. We suggest some extensions of our method to other systems.

The translational efficiencies of the two Leishmania infantum HSP70 mRNAs, differing in their 3'-untranslated regions, are affected by shifts in the temperature of growth through different mechanisms

Exposure of Leishmania promastigotes to the temperature of their mammalian hosts induces a typical heat-shock response. In Leishmania infantum, HSP70 is encoded by two types of genes that differ in their 3'-untranslated regions (3'-UTRs). Previously, we have shown that specific transcripts for each gene are present in promastigotes growing at normal temperature (26 degrees C), but only transcripts with 3'-UTR-type I (3'-UTRI) accumulate in a temperature-dependent manner. Here, we have investigated the translational efficiencies of both types of HSP70 transcripts at the different temperatures that the parasite encounters in the insect (26 degrees C, normal temperature) or in the mammalian host (heat-shock temperatures). Interestingly, 3'-UTRI-bearing transcripts (HSP70-I) were found associated with ribosomes in promastigotes at normal and heat-shock temperatures, whereas the HSP70-II transcripts appear to be preferentially translated at heat-shock temperatures but not at 26 degrees C. We have analyzed the function of these UTRs in the translational control by use of plasmid constructs in which the CAT reporter gene was flanked by UTRs of the HSP70 genes. Unexpectedly, it was found that CAT transcripts with 3'-UTRII bind to ribosomes at 26 degrees C, and, indeed, the CAT protein is synthesized. A valid conclusion of these experiments was that both types of 3'-UTRs are essential for translation of HSP70 mRNAs at heat shock temperatures, although the 3'-UTRII is more efficient during severe heat shock (39 degrees C). In addition, these results suggest that sequence region other than the 3'-UTR of HSP70-II gene is involved in the translational silent state of HSP70-II transcripts at 26 degrees C. Finally, a null mutant has been created by targeted disruption of both HSP70-II alleles. Remarkably, the deltaHSP70 mutant synthesizes HSP70 at a lower rate than the wild-type parasites. Overall, our data suggest that the biological function of the HSP70-II gene is to top up HSP70 levels under conditions of stress.

Pre-mRNA trans-splicing: from kinetoplastids to mammals, an easy language for life diversity

Since the discovery that genes are split into intron and exons, the studies of the mechanisms involved in splicing pointed to presence of consensus signals in an attempt to generalize the process for all living cells. However, as discussed in the present review, splicing is a theme full of variations. The trans-splicing of pre-mRNAs, the joining of exons from distinct transcripts, is one of these variations with broad distribution in the phylogenetic tree. The biological meaning of this phenomenon is discussed encompassing reactions resembling a possible noise to mechanisms of gene expression regulation. All of them however, can contribute to the generation of life diversity.

Comparative genomics of trypanosomatid parasitic protozoa

A comparison of gene content and genome architecture of Trypanosoma brucei, Trypanosoma cruzi, and Leishmania major, three related pathogens with different life cycles and disease pathology, revealed a conserved core proteome of about 6200 genes in large syntenic polycistronic gene clusters. Many species-specific genes, especially large surface antigen families, occur at nonsyntenic chromosome-internal and subtelomeric regions. Retroelements, structural RNAs, and gene family expansion are often associated with syntenic discontinuities that-along with gene divergence, acquisition and loss, and rearrangement within the syntenic regions-have shaped the genomes of each parasite. Contrary to recent reports, our analyses reveal no evidence that these species are descended from an ancestor that contained a photosynthetic endosymbiont.

The Trypanosoma cruzi Proteome

To complement the sequencing of the three kinetoplastid genomes reported in this issue, we have undertaken a whole-organism, proteomic analysis of the four life-cycle stages of Trypanosoma cruzi. Peptides mapping to 2784 proteins in 1168 protein groups from the annotated T. cruzi genome were identified across the four life-cycle stages. Protein products were identified from >1000 genes annotated as "hypothetical" in the sequenced genome, including members of a newly defined gene family annotated as mucin-associated surface proteins. The four parasite stages appear to use distinct energy sources, including histidine for stages present in the insect vectors and fatty acids by intracellular amastigotes.

The genome of the kinetoplastid parasite, Leishmania major

Leishmania species cause a spectrum of human diseases in tropical and subtropical regions of the world. We have sequenced the 36 chromosomes of the 32.8-megabase haploid genome of Leishmania major (Friedlin strain) and predict 911 RNA genes, 39 pseudogenes, and 8272 protein-coding genes, of which 36% can be ascribed a putative function. These include genes involved in host-pathogen interactions, such as proteolytic enzymes, and extensive machinery for synthesis of complex surface glycoconjugates. The organization of protein-coding genes into long, strand-specific, polycistronic clusters and lack of general transcription factors in the L. major, Trypanosoma brucei, and Trypanosoma cruzi (Tritryp) genomes suggest that the mechanisms regulating RNA polymerase II-directed transcription are distinct from those operating in other eukaryotes, although the trypanosomatids appear capable of chromatin remodeling. Abundant RNA-binding proteins are encoded in the Tritryp genomes, consistent with active posttranscriptional regulation of gene expression.

Protozoan genomes: gene identification and annotation

The draft sequence of several complete protozoan genomes is now available and genome projects are ongoing for a number of other species. Different strategies are being implemented to identify and annotate protein coding and RNA genes in these genomes, as well as study their genomic architecture. Since the genomes vary greatly in size, GC-content, nucleotide composition, and degree of repetitiveness, genome structure is often a factor in choosing the methodology utilised for annotation. In addition, the approach taken is dictated, to a greater or lesser extent, by the particular reasons for carrying out genome-wide analyses and the level of funding available for projects. Nevertheless, these projects have provided a plethora of material that will aid in understanding the biology and evolution of these parasites, as well as identifying new targets that can be used to design urgently required drug treatments for the diseases they cause.

Characterization and developmental gene regulation of a large gene family encoding amastin surface proteins in Leishmania spp

The ability of Leishmania amastigotes to survive within the drastic environmental changes encountered in the phagolysosomes of mammalian macrophages is heavily dependent on the developmental regulation of a variety of genes. The identification of genes that are expressed preferentially in the mammalian stage of the parasite should increase our understanding of the molecular mechanisms regulating stage-specific gene expression and of the determinants that control its intracellular survival and contribute to its pathogenesis. We report here detailed sequence characterization and structural organization of the amastin gene family in Leishmania major and Leishmania infantum and the study of their developmental gene regulation throughout the parasite's life cycle. Amastin surface proteins represent the largest developmentally regulated gene family reported so far in Leishmania comprising up to 45 members. All the members of the amastin gene family in both Leishmania and Trypanosoma species share a similar structural organization and contain a highly conserved 11 amino acid extracellular domain, which is unique to amastin proteins. The majority of the amastin gene homologs are specifically expressed in the amastigote stage of the parasite. Three distinct RNA elements were identified in the 3'-untranslated regions (3'UTR) of the amastin transcripts. The majority of these transcripts contain a conserved 450 nt cis-acting 3'UTR element shown previously to regulate stage-specific gene expression at the level of translation, which suggests that several amastin homologs may be regulated by a similar mechanism of translational control inside the macrophage. These findings further highlight the unique features of gene expression control in Leishmania.

Strand asymmetry patterns in trypanosomatid parasites

The genome organization of kinetoplastid parasites is unusual, with chromosomes containing several long regions of polycistronically transcribed genes. The regions where the direction of transcription switches have been hypothesized to contain origins of replication and possibly also centromers and promoters. We report that overall strand asymmetry patterns can be observed in Trypanosoma cruzi and Trypanosoma brucei with optima on strand-switch regions. The base skews of T. cruzi and T. brucei divergent strand-switches show patterns analogous to those for bacterial origins of replication, but they differ from those of Leishmania major. Bias in codon usage and the trypanosomatid unidirectional gene clusters predict most of this skew, but fail to properly explain the same trend in intergenic regions, as does the current knowledge of regulatory sequences.

A survey of Leishmania braziliensis genome by shotgun sequencing

We have carried out a survey of the genome of Leishmania (Viannia) braziliensis by shotgun sequencing. Approximately 15% of the haploid genome of the parasite (5.15 Mb of genomic sequence) was obtained. A large number of known and putative genes, predicted to be involved in several cellular processes, were identified. Some genomic features were investigated, such as the general G + C content, which was found to be lower than L. major (57% versus 63%). BlastN searches revealed that 60.2% of the clusterized GSS sequences displayed similarity to L. major genomic sequences, while a BlastX search showed that 45.3% of the thus obtained predicted protein sequences showed similarity to annotated proteins of L. major. Further comparison of the degree of conservation between L. major and L. braziliensis revealed that coding regions are much more conserved than non-coding ones. The shotgun sequence analysis of Leishmania braziliensis appears to be an efficient and suitable strategy contributing to the search for vaccines and novel drug targets. The sequence data described in this paper have been submitted to the dbGSS database under the following accession numbers (BX530413 to BX530454; BX530456 to BX530718; BX538354 to BX539305; BX539350 to BX540325; BX541002 to BX544869; BX544893 to BX545685; BX897701 to BX897710; BX905184 to BX907797; BX907798 to BX908381; BX908403 to BX908718). All data including sequences are also available at (www.ebi.ac.uk/embl/).