Members of a large retroposon family are determinants of post-transcriptional gene expression in Leishmania

R-loops acted on by RNase H1 influence DNA replication timing and genome stability in Leishmania

Genomes in eukaryotes normally undergo DNA replication in a choreographed temporal order, resulting in early and late replicating chromosome compartments. Leishmania, a human protozoan parasite, displays an unconventional DNA replication program in which the timing of DNA replication completion is chromosome size-dependent: larger chromosomes complete replication later then smaller ones. Here we show that both R-loops and RNase H1, a ribonuclease that resolves RNA-DNA hybrids, accumulate in Leishmania major chromosomes in a pattern that reflects their replication timing. Furthermore, we demonstrate that such differential organisation of R-loops, RNase H1 and DNA replication timing across the parasite's chromosomes correlates with size-dependent differences in chromatin accessibility, G quadruplex distribution and sequence content. Using conditional gene excision, we show that loss of RNase H1 leads to transient growth perturbation and permanently abrogates the differences in DNA replication timing across chromosomes, as well as altering levels of aneuploidy and increasing chromosome instability in a size-dependent manner. This work provides a link between R-loop homeostasis and DNA replication timing in a eukaryotic parasite and demonstrates that orchestration of DNA replication dictates levels of genome plasticity in Leishmania.

Evolutionary divergent clusters of transcribed extinct truncated retroposons drive low mRNA expression and developmental regulation in the protozoan Leishmania

BACKGROUND

The Leishmania genome harbors formerly active short interspersed degenerated retroposons (SIDERs) representing the largest family of repetitive elements among trypanosomatids. Their substantial expansion in Leishmania is a strong predictor of important biological functions. In this study, we combined multilevel bioinformatic predictions with high-throughput genomic and transcriptomic analyses to gain novel insights into the diversified roles retroposons of the SIDER2 subfamily play in Leishmania genome evolution and expression.

RESULTS

We show that SIDER2 retroposons form various evolutionary divergent clusters, each harboring homologous SIDER2 sequences usually located nearby in the linear sequence of chromosomes. This intriguing genomic organization underscores the importance of SIDER2 proximity in shaping chromosome dynamics and co-regulation. Accordingly, we show that transcripts belonging to the same SIDER2 cluster can display similar levels of expression. SIDER2 retroposons are mostly transcribed as part of 3'UTRs and account for 13% of the Leishmania transcriptome. Genome-wide expression profiling studies underscore SIDER2 association generally with low mRNA expression. The remarkable link of SIDER2 retroposons with downregulation of gene expression supports their co-option as major regulators of mRNA abundance. SIDER2 sequences also add to the diversification of the Leishmania gene expression repertoire since ~ 35% of SIDER2-containing transcripts can be differentially regulated throughout the parasite development, with a few encoding key virulence factors. In addition, we provide evidence for a functional bias of SIDER2-containing transcripts with protein kinase and transmembrane transporter activities being most represented.

CONCLUSIONS

Altogether, these findings provide important conceptual advances into evolutionary innovations of transcribed extinct retroposons acting as major RNA cis-regulators.

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

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

Life in plastic, it's fantastic! How Leishmania exploit genome instability to shape gene expression

Leishmania are kinetoplastid pathogens that cause leishmaniasis, a debilitating and potentially life-threatening infection if untreated. Unusually, Leishmania regulate their gene expression largely post-transcriptionally due to the arrangement of their coding genes into polycistronic transcription units that may contain 100s of functionally unrelated genes. Yet, Leishmania are capable of rapid and responsive changes in gene expression to challenging environments, often instead correlating with dynamic changes in their genome composition, ranging from chromosome and gene copy number variations to the generation of extrachromosomal DNA and the accumulation of point mutations. Typically, such events indicate genome instability in other eukaryotes, coinciding with genetic abnormalities, but for Leishmania, exploiting these products of genome instability can provide selectable substrates to catalyse necessary gene expression changes by modifying gene copy number. Unorthodox DNA replication, DNA repair, replication stress factors and DNA repeats are recognised in Leishmania as contributors to this intrinsic instability, but how Leishmania regulate genome plasticity to enhance fitness whilst limiting toxic under- or over-expression of co-amplified and co-transcribed genes is unclear. Herein, we focus on fresh, and detailed insights that improve our understanding of genome plasticity in Leishmania. Furthermore, we discuss emerging models and factors that potentially circumvent regulatory issues arising from polycistronic transcription. Lastly, we highlight key gaps in our understanding of Leishmania genome plasticity and discuss future studies to define, in higher resolution, these complex regulatory interactions.

The non-LTR retrotransposons of Entamoeba histolytica: genomic organization and biology

Genome sequence analysis of Entamoeba species revealed various classes of transposable elements. While E. histolytica and E. dispar are rich in non-long terminal repeat (LTR) retrotransposons, E. invadens contains predominantly DNA transposons. Non-LTR retrotransposons of E. histolytica constitute three families of long interspersed nuclear elements (LINEs), and their short, nonautonomous partners, SINEs. They occupy ~ 11% of the genome. The EhLINE1/EhSINE1 family is the most abundant and best studied. EhLINE1 is 4.8 kb, with two ORFs that encode functions needed for retrotransposition. ORF1 codes for the nucleic acid-binding protein, and ORF2 has domains for reverse transcriptase (RT) and endonuclease (EN). Most copies of EhLINEs lack complete ORFs. ORF1p is expressed constitutively, but ORF2p is not detected. Retrotransposition could be demonstrated upon ectopic over expression of ORF2p, showing that retrotransposition machinery is functional. The newly retrotransposed sequences showed a high degree of recombination. In transcriptomic analysis, RNA-Seq reads were mapped to individual EhLINE1 copies. Although full-length copies were transcribed, no full-length 4.8 kb transcripts were seen. Rather, sense transcripts mapped to ORF1, RT and EN domains. Intriguingly, there was strong antisense transcription almost exclusively from the RT domain. These unique features of EhLINE1 could serve to attenuate retrotransposition in E. histolytica.

TBP and SNAP50 transcription factors bind specifically to the Pr77 promoter sequence from trypanosomatid non-LTR retrotransposons

Background

Trypanosomatid genomes are colonized by active and inactive mobile DNA elements, such as LINE, SINE-like, SIDER and DIRE retrotransposons. These elements all share a 77-nucleotide-long sequence at their 5′ ends, known as Pr77, which activates transcription, thereby generating abundant unspliced and translatable transcripts. However, transcription factors that mediates this process have still not been reported.

Methods

TATA-binding protein (TBP) and small nuclear RNA-activating protein 50 kDa (SNAP50) recombinant proteins and specific antibodies raised against them were generated. Protein capture assay, electrophoretic mobility-shift assays (EMSA) and EMSA competition assays carried out using these proteins and nuclear proteins of the parasite together to specific DNA sequences used as probes allowed detecting direct interaction of these transcription factors to Pr77 sequence.

Results

This study identified TBP and SNAP50 as part of the DNA-protein complex formed by the Pr77 promoter sequence and nuclear proteins of Trypanosoma cruzi. TBP establishes direct and specific contact with the Pr77 sequence, where the DPE and DPE downstream regions are docking sites with preferential binding. TBP binds cooperatively (Hill coefficient = 1.67) to Pr77 and to both strands of the Pr77 sequence, while the conformation of this highly structured sequence is not involved in TBP binding. Direct binding of SNAP50 to the Pr77 sequence is weak and may be mediated by protein–protein interactions through other trypanosomatid nuclear proteins.

Conclusions

Identification of the transcription factors that mediate Pr77 transcription may help to elucidate how these retrotransposons are mobilized within the trypanosomatid genomes and their roles in gene regulation processes in this human parasite.

Graphic abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-021-04803-5.

Epitranscriptome machinery in Trypanosomatids: New players on the table?

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

Transcriptomic analysis of Entamoeba histolytica reveals domain-specific sense strand expression of LINE-encoded ORFs with massive antisense expression of RT domain

LINEs are retrotransposable elements found in diverse organisms. Their activity is kept in check by several mechanisms, including transcriptional silencing. Here we have analyzed the transcription status of LINE1 copies in the early-branching parasitic protist Entamoeba histolytica. Full-length EhLINE1 encodes ORF1, and ORF2 with reverse transcriptase (RT) and endonuclease (EN) domains. RNA-Seq analysis of EhLINE1 copies (both truncated and full-length) showed unique features. Firstly, although 20/41 transcribed copies were full-length, we failed to detect any full-length transcripts. Rather, sense-strand transcripts mapped to the functional domains- ORF1, RT and EN. Secondly, there was strong antisense transcription specifically from RT domain. No antisense transcripts were seen from ORF1. Antisense RT transcripts did not encode known functional peptides. They could possibly be involved in attenuating translation of RT domain, as we failed to detect ORF2p, whereas ORF1p was detectable. Lack of full-length transcripts and strong antisense RT expression may serve to limit EhLINE1 retrotransposition.

Leishmania: Responding to environmental signals and challenges without regulated transcription

Here we describe the non-canonical control of gene expression in Leishmania, a single-cell parasite that is responsible for one of the major neglected tropical diseases. We discuss the lack of regulated RNA synthesis, the post-transcriptional gene regulation including RNA stability and regulated translation. We also show that genetic adaptations such as mosaic aneuploidy, gene copy number variations and DNA sequence polymorphisms are important means for overcoming drug challenge and environmental diversity. These mechanisms are discussed in the context of the unique flow of genetic information found in Leishmania and related protists.

The Leishmania donovani SENP Protease Is Required for SUMO Processing but Not for Viability

The protozoan parasite Leishmania donovani is part of an early eukaryotic branch and depends on post-transcriptional mechanisms for gene expression regulation. This includes post-transcriptional protein modifications, such as protein phosphorylation. The presence of genes for protein SUMOylation, i.e., the covalent attachment of small ubiquitin-like modifier (SUMO) polypeptides, in the Leishmania genomes prompted us to investigate the importance of the sentrin-specific protease (SENP) and its putative client, SUMO, for the vitality and infectivity of Leishmania donovani. While SENP null mutants are viable with reduced vitality, viable SUMO null mutant lines could not be obtained. SUMO C-terminal processing is disrupted in SENP null mutants, preventing SUMO from covalent attachment to proteins and nuclear translocation. Infectivity in vitro is not affected by the loss of SENP-dependent SUMO processing. We conclude that SENP is required for SUMO processing, but that functions of unprocessed SUMO are critical for Leishmania viability.

Casein kinase 1.2 over expression restores stress resistance to Leishmania donovani HSP23 null mutants

Leishmania donovani is a trypanosomatidic parasite and causes the lethal kala-azar fever, a neglected tropical disease. The Trypanosomatida are devoid of transcriptional gene regulation and rely on gene copy number variations and translational control for their adaption to changing conditions. To survive at mammalian tissue temperatures, L. donovani relies on the small heat shock protein HSP23, the loss of which renders the parasites stress sensitive and impairs their proliferation. Here, we analysed a spontaneous escape mutant with wild type-like in vitro growth. Further selection of this escape strains resulted in a complete reversion of the phenotype. Whole genome sequencing revealed a correlation between stress tolerance and the massive amplification of a six-gene cluster on chromosome 35, with further analysis showing over expression of the casein kinase 1.2 gene as responsible. In vitro phosphorylation experiments established both HSP23 and the related P23 co-chaperone as substrates and modulators of casein kinase 1.2, providing evidence for another crucial link between chaperones and signal transduction protein kinases in this early branching eukaryote.

The mRNA-bound Proteome of Leishmania mexicana: Novel Genetic Insight into an Ancient Parasite

Leishmania parasite infections, termed the leishmaniases, cause significant global infectious disease burden. The lifecycle of the parasite embodies three main stages that require precise coordination of gene regulation to survive environmental shifts between sandfly and mammalian hosts. Constitutive transcription in kinetoplastid parasites means that gene regulation is overwhelmingly reliant on post-transcriptional mechanisms, yet strikingly few Leishmania trans-regulators are known. Using optimized crosslinking and deep, quantified mass spectrometry, we present a comprehensive analysis of 1400 mRNA binding proteins (mRBPs) and whole cell proteomes from the three main Leishmania lifecycle stages. Supporting the validity, although the crosslinked RBPome is magnitudes more enriched, the protein identities of the crosslinked and non-crosslinked RBPomes were nearly identical. Moreover, multiple candidate RBPs were endogenously tagged and found to associate with discrete mRNA target pools in a stage-specific manner. Results indicate that in L. mexicana parasites, mRNA levels are not a strong predictor of the whole cell expression or RNA binding potential of encoded proteins. Evidence includes a low correlation between transcript and corresponding protein expression and stage-specific variation in protein expression versus RNA binding potential. Unsurprisingly, RNA binding protein enrichment correlates strongly with relative replication efficiency of the specific lifecycle stage. Our study is the first to quantitatively define and compare the mRBPome of multiple stages in kinetoplastid parasites. It provides novel, in-depth insight into the trans-regulatory mRNA:Protein (mRNP) complexes that drive Leishmania parasite lifecycle progression.

A Leishmania infantum genetic marker associated with miltefosine treatment failure for visceral leishmaniasis

Background

Miltefosine has been used successfully to treat visceral leishmaniasis (VL) in India, but it was unsuccessful for VL in a clinical trial in Brazil.

Methods

To identify molecular markers that predict VL treatment failure whole genome sequencing of 26 L. infantum isolates, from cured and relapsed patients allowed a GWAS analysis of SNPs, gene and chromosome copy number variations.

Findings

A strong association was identified (p = 0·0005) between the presence of a genetically stable L. infantumMiltefosine Sensitivity Locus (MSL), and a positive response to miltefosine treatment. The risk of treatment failure increased 9·4-fold (95% CI 2·11–53·54) when an isolate did not have the MSL. The complete absence of the MSL predicted miltefosine failure with 0·92 (95% CI 0·65–0·996) sensitivity and 0·78 (95% CI 0·52–0·92) specificity. A genotyping survey of L. infantum (n = 157) showed that the frequency of MSL varies in a cline from 95% in North East Brazil to <5% in the South East. The MSL was found in the genomes of all L. infantum and L. donovani sequenced isolates from the Old World (n = 671), where miltefosine can have a cure rate higher than 93%.

Interpretation

Knowledge on the presence or absence of the MSL in L. infantum will allow stratification of patients prior to treatment, helping to establish better therapeutic strategies for VL treatment.

Fund

CNPq, FAPES, GCRF MRC and Wellcome Trust.

Biology of Trypanosoma cruzi Retrotransposons: From an Enzymatic to a Structural Point of View

INTRODUCTION

An important portion of the Trypanosoma cruzi genome is composed of mobile genetic elements, which are interspersed with genes on all chromosomes. The L1Tc non-LTR retrotransposon and its truncated version NARTc are the most highly represented and best studied of these elements. L1Tc is actively transcribed in all three forms of the Trypanosoma parasite and encodes the proteins that enable it to autonomously mobilize. This mini review discusses the enzymatic properties of L1Tc that enable its mobilization and possibly the mobilization of other non-autonomous retrotransposons in Trypanosoma. We also briefly review the Hepatitis Delta Virus-like autocatalytic and 2A self-cleaving viral-like sequences contained in L1Tc that regulate post-transcriptional properties such as relative protein abundance and mRNA stability. Special emphasis is placed on the Pr77 dual system, which is based on the RNA pol II-dependent internal promoter of L1Tc and NARTc and the HDV-like ribozyme activity encoded by the first 77 nucleotides of the element's DNA and RNA. The high degree of conservation of the Pr77 sequence, referred to as the "Pr77-hallmark", among different trypanosomatid retroelements suggests that these mobile elements are responsible for the distribution of regulatory sequences within the genome they inhabit.

CONCLUSION

We also discuss how the involvement of L1Tc and NARTc in the gene regulatory processes of these parasites could justify their domestication and long-term coexistence in these ancient organisms.

Gene and Chromosomal Copy Number Variations as an Adaptive Mechanism Towards a Parasitic Lifestyle in Trypanosomatids

Trypanosomatids are a group of kinetoplastid parasites including some of great public health importance, causing debilitating and life-long lasting diseases that affect more than 24 million people worldwide. Among the trypanosomatids, Trypanosoma cruzi, Trypanosoma brucei and species from the Leishmania genus are the most well studied parasites, due to their high prevalence in human infections. These parasites have an extreme genomic and phenotypic variability, with a massive expansion in the copy number of species-specific multigene families enrolled in host-parasite interactions that mediate cellular invasion and immune evasion processes. As most trypanosomatids are heteroxenous, and therefore their lifecycles involve the transition between different hosts, these parasites have developed several strategies to ensure a rapid adaptation to changing environments. Among these strategies, a rapid shift in the repertoire of expressed genes, genetic variability and genome plasticity are key mechanisms. Trypanosomatid genomes are organized into large directional gene clusters that are transcribed polycistronically, where genes derived from the same polycistron may have very distinct mRNA levels. This particular mode of transcription implies that the control of gene expression operates mainly at post-transcriptional level. In this sense, gene duplications/losses were already associated with changes in mRNA levels in these parasites. Gene duplications also allow the generation of sequence variability, as the newly formed copy can diverge without loss of function of the original copy. Recently, aneuploidies have been shown to occur in several Leishmania species and T. cruzi strains. Although aneuploidies are usually associated with debilitating phenotypes in superior eukaryotes, recent data shows that it could also provide increased fitness in stress conditions and generate drug resistance in unicellular eukaryotes. In this review, we will focus on gene and chromosomal copy number variations and their relevance to the evolution of trypanosomatid parasites.

The Pumilio-domain protein PUF6 contributes to SIDER2 retroposon-mediated mRNA decay in Leishmania

Leishmania and other trypanosomatid protozoa lack control at the level of transcription initiation and regulate gene expression exclusively post-transcriptionally. We have reported previously that Leishmania harbors a unique class of short interspersed degenerate retroposons (SIDERs) that are predominantly located within 3'UTRs and play a major role in post-transcriptional control. We have shown that members of the SIDER2 subfamily initiate mRNA decay through endonucleolytic cleavage within the second conserved 79-nt signature sequence of SIDER2 retroposons. Here, we have developed an optimized MS2 coat protein tethering system to capture trans-acting factor(s) regulating SIDER2-mediated mRNA decay. Tethering of the MS2 coat protein to a reporter RNA harboring two MS2 stem-loop aptamers and the cognate SIDER2-containing 3'UTR in combination with immunoprecipitation and mass spectrometry analysis led to the identification of RNA-binding proteins with known functions in mRNA decay. Among the candidate SIDER2-interacting proteins that were individually tethered to a SIDER2 reporter RNA, the Pumilio-domain protein PUF6 was shown to enhance degradation and reduce transcript half-life. Furthermore, we showed that PUF6 binds to SIDER2 sequences that include the regulatory 79-nt signature motif, hence contributing to the mRNA decay process. Consistent with a role of PUF6 in SIDER2-mediated decay, genetic inactivation of PUF6 resulted in increased accumulation and higher stability of endogenous SIDER2-bearing transcripts. Overall, these studies provide new insights into regulated mRNA decay pathways in Leishmania controlled by SIDER2 retroposons and propose a broader role for PUF proteins in mRNA decay within the eukaryotic kingdom.

Identification of the centromeres of Leishmania major: revealing the hidden pieces

Leishmania affects millions of people worldwide. Its genome undergoes constitutive mosaic aneuploidy, a type of genomic plasticity that may serve as an adaptive strategy to survive distinct host environments. We previously found high rates of asymmetric chromosome allotments during mitosis that lead to the generation of such ploidy. However, the underlying molecular events remain elusive. Centromeres and kinetochores most likely play a key role in this process, yet their identification has failed using classical methods. Our analysis of the unconventional kinetochore complex recently discovered in Trypanosoma brucei (KKTs) leads to the identification of a Leishmania KKT gene candidate (LmKKT1). The GFP-tagged LmKKT1 displays "kinetochore-like" dynamics of intranuclear localization throughout the cell cycle. By ChIP-Seq assay, one major peak per chromosome is revealed, covering a region of 4 ±2 kb. We find two largely conserved motifs mapping to 14 of 36 chromosomes while a higher density of retroposons are observed in 27 of 36 centromeres. The identification of centromeres and of a kinetochore component of Leishmania chromosomes opens avenues to explore their role in mosaic aneuploidy.

Disclosing 3' UTR cis-elements and putative partners involved in gene expression regulation in Leishmania spp

To identify putative cis-elements involved in gene expression regulation in Leishmania, we previously conducted an in silico investigation to find conserved intercoding sequences (CICS) in the genomes of L. major, L. infantum, and L. braziliensis. Here, the CICS databank was explored to search for sequences that were present in the untranslated regions (UTRs) of groups of genes showing similar expression profiles during in vitro differentiation. Using a selectable marker as a reporter gene, flanked by either an intact 3' UTR or a UTR lacking the conserved element, the regulatory role of a CICS was confirmed. We observed that the pattern of modulation of the mRNA levels was altered in the absence of the CICS. We also identified putative CICS RNA-binding proteins. This study suggests that the publicly available CICS database is a useful tool for identifying regulatory cis-elements for Leishmania genes and suggests the existence of post-transcriptional regulons in Leishmania.

RNA secondary structure and nucleotide composition of the conserved hallmark sequence of Leishmania SIDER2 retroposons are essential for endonucleolytic cleavage and mRNA degradation

We have reported previously that Short Interspersed Degenerate Retroposons of the SIDER2 subfamily, largely located within 3'UTRs of Leishmania transcripts, promote rapid turnover of mRNAs through endonucleolytic cleavage within the highly conserved second tandem 79-nt hallmark sequence (79-nt SII). Here, we used site-directed mutagenesis and in silico RNA structural studies to delineate the cis-acting requirements within 79-nt SII for cleavage and mRNA degradation. The putative cleavage site(s) and other nucleotides predicted to alter the RNA secondary structure of 79-nt SII were either deleted or mutated and their effect on mRNA turnover was monitored using a gene reporter system. We found that short deletions of 8-nt spanning the two predicted cleavage sites block degradation of SIDER2-containing transcripts, leading to mRNA accumulation. Furthermore, single or double substitutions of the dinucleotides targeted for cleavage as well as mutations altering the predicted RNA secondary structure encompassing both cleavage sites also prevent mRNA degradation, confirming that these dinucleotides are the bona fide cleavage sites. In line with these results, we show that stage-regulated SIDER2 inactivation correlates with the absence of endonucleolytic cleavage. Overall, these data demonstrate that both cleavage sites within the conserved 79-nt SII as well as RNA folding in this region are essential for SIDER2-mediated mRNA decay, and further support that SIDER2-harboring transcripts are targeted for degradation by endonucleolytic cleavage.

SIDER2 retroposon-mediated mRNA decay in Leishmania is coupled to translation

We previously reported that Short Interspersed Degenerate Retroposons of the SIDER2 subfamily predominantly located within 3' untranslated regions (UTRs) of Leishmania transcripts promote rapid turnover that is initiated by endonucleolytic cleavage. Here, we investigated whether SIDER2-mediated mRNA decay is linked to translation. We show that preventing translation initiation by inserting a hairpin structure at the 5'-end of a SIDER2-containing mRNA blocks degradation. Similarly, global inhibition of translation elongation by cycloheximide or termination by puromycin causes stabilisation of SIDER2-containing transcripts. Altogether, these findings support that the mechanism of SIDER2-mediated decay is coupled to translation, possibly through the recruitment of decay factors to elongating ribosomes.

Dataset for distribution of SIDER2 elements in the Leishmania major genome and transcriptome

This paper contains data related to the research article entitled “Genomic cartography and proposal of nomenclature for the repeated, interspersed elements of the Leishmania major SIDER2 family and identification of SIDER2-containing transcripts” [1]. SIDER2 elements are repeated sequences, derived from, nowadays, extinct retrotransposons, that populate the genomes of protist of the genera Leishmania. This dataset (Supplementary file 1), an inventory of 1100 SIDER2 elements, was generated by surveying the L. major complete genome using bioinformatics tools with further manual refinements. In addition to the genomic distribution of these elements (summarized in Fig. 1), this dataset contains information regarding their association with specific transcripts, based on the recently established transcriptome for L. major[2].

Genomic cartography and proposal of nomenclature for the repeated, interspersed elements of the Leishmania major SIDER2 family and identification of SIDER2-containing transcripts

The genomes of most eukaryotic organisms contain a large number of transposable elements that are able to move from one genomic site to another either by transferring of DNA mobile elements (transposons) or transpose via reverse transcription of an RNA intermediate (retroposons). An exception to this rule is found in protists of the subgenus Leishmania, in which active retroposons degenerated after a flourishing era, leaving only retroposon remains; these have been classified into two families: SIDER1 and SIDER2. In this work, we have re-examined the elements belonging to the family SIDER2 present in the genome of Leishmania major with the aim of providing a nomenclature that will facilitate a future reference to particular elements. According to sequence conservation, the 1100 SIDER2 elements have been grouped into subfamilies, and the inferred taxonomic relationships have also been incorporated into the nomenclature. Additionally, we are providing detailed data regarding the genomic distribution of these elements and their association with specific transcripts, based on the recently established transcriptome for L. major. Thus, the presented data can help to study and better understand the roles played by these degenerated retroposons in both regulation of gene expression and genome plasticity.

Genomic and phylogenetic evidence of VIPER retrotransposon domestication in trypanosomatids

Transposable elements are important residents of eukaryotic genomes and eventually the host can domesticate them to serve cellular functions. We reported here a possible domestication event of the vestigial interposed retroelement (VIPER) in trypanosomatids. We found a large gene in a syntenic location in Leishmania braziliensis, L. panamensis, Leptomanas pyrrhocoris, and Crithidia fasciculata whose products share similarity in the C-terminal portion with the third protein of VIPER. No remnants of other VIPER regions surrounding the gene sequence were found. We hypothesise that the domestication event occurred more than 50 mya and the conservation of this gene suggests it might perform some function in the host species.

Plasticity of the Leishmania genome leading to gene copy number variations and drug resistance

Leishmania has a plastic genome, and drug pressure can select for gene copy number variation (CNV). CNVs can apply either to whole chromosomes, leading to aneuploidy, or to specific genomic regions. For the latter, the amplification of chromosomal regions occurs at the level of homologous direct or inverted repeated sequences leading to extrachromosomal circular or linear amplified DNAs. This ability of Leishmania to respond to drug pressure by CNVs has led to the development of genomic screens such as Cos-Seq, which has the potential of expediting the discovery of drug targets for novel promising drug candidates.

The Trypanosomatid Pr77-hallmark contains a downstream core promoter element essential for transcription activity of the Trypanosoma cruzi L1Tc retrotransposon

Background

Trypanosomatid genomes are highly colonized by non-LTR retroelements that make up to 5 % of the nuclear genome. These elements are mainly accumulated in the strand switch regions (SSRs) where polycistronic transcription is initiated and have a 77 nt-long sequence - Pr77 - at their 5′ ends. L1Tc is the best represented retrotransposon in the Trypanosoma cruzi genome and is a potentially functional autonomous element that encodes its own retrotransposition machinery. The Pr77 of the T. cruzi L1Tc element activates gene transcription via RNA polymerase II, generating abundant, unspliced transcripts which are translated.

Results

The present manuscript describes the identification of a downstream core promoter element (DPE) in the L1Tc Pr77 sequence. Just four nucleotides long (CGTG), it covers in Pr77 positions +25 to +28 of the described L1Tc transcription start site. The Pr77-DPE motif is conserved in terms of sequence composition and position in the Pr77 of most trypanosomatid non-LTR retrotransposons, independent of the coding or non-coding capacity of these retroelements. Transcription assays in T. cruzi stable transfectants with vector containing point mutations at 17 locations of the Pr77 nucleotide sequence evidence that the DPE motif is essential for the promoter function of Pr77. Furthermore, the obtained data show that other nucleotides also contributed to the promoter function of Pr77. In addition, the presented results indicate that parasite nuclear proteins specifically bind to different regions of the Pr77 sequence although the strongest binding is to the DPE motif. Moreover, it is shown that the DPE sense single-stranded sequence is being required in DNA-protein recognition of nuclear factors.

Conclusions

The Pr77 sequence present in most of non-LTR retrotransposons of trypanosomatids contains a downstream core promoter element (DPE) which is conserved in terms of nucleotide composition and location. The Pr77-DPE motif is essential for the transcriptional activity of Pr77 although other nucleotides are also involved. DPE has a high affinity binding for nuclear proteins in T. cruzi. The wide retroelement-mediated distribution of Pr77 suggests that it may represent an important tool for regulating gene expression in trypanosomatids.

Electronic supplementary material

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

Resequencing and assembly of seven complex loci to improve the Leishmania major (Friedlin strain) reference genome

Background

Leishmania parasites cause severe human diseases known as leishmaniasis. These eukaryotic microorganisms possess an atypical chromosomal architecture and the regulation of gene expression occurs almost exclusively at post-transcriptional levels. Accordingly, sequencing of the genome of Leishmania major, and subsequently the genome of other related species, was paramount for highlighting these peculiar molecular aspects. Recently, we carried out an analysis of gene expression by massive sequencing of RNA in the L. major promastigote, and data derived from that analysis were suggestive of possible errors in the current genome assembly for this Leishmania species.

Results

During the analysis by RNA-Seq of the transcriptome for L. major Friedlin strain, 163,714 reads could not be aligned with the reference genome. Thus, de novo assembly with these reads was carried out and the resulting contigs were further analyzed. After detailed homology searches using available databases, it was postulated that 15 contigs might correspond to genomic sequences lost during the initial genome assembly of the L. major Friedlin strain. This was experimentally confirmed by PCR amplification, cloning and sequencing of the new genomic regions. As a result, we have identified seven regions of the L. major (Friedlin) genome that were lost during the sequence assembly. This led to the uncovering of six new genes (LmjF.15.1475, LmjF.15.0285, LmjF.24.0765, LmjF.14.0860, LmjF.19.0305, and LmjF.27.2035), and correction of the annotation for two others (LmjF.15.1480 and LmjF.27.2030). Our data suggest that these genomic regions probably collapsed during the genome assembly due to the existence of gene duplications and/or repeated regions surrounding the missed genes.

Conclusion

RNA-seq data helped to reconstruct some genomic regions misassembled during the L. major Friedlin genome assembly, which is otherwise quite robust. On the other hand, this study shows that data derived from massive sequencing approaches, including RNA-Seq, should be carefully inspected to improve current genome definition and gene annotations.

Electronic supplementary material

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

Integrating ribosomal promoter vectors that offer a choice of constitutive expression profiles in Leishmania donovani

We have designed a novel series of integrating ribosomal RNA promoter vectors with five incrementally different constitutive expression profiles, covering a 250-fold range. Differential expression was achieved by placing different combinations of synthetic or leishmanial DNA sequences upstream and downstream of the transgene coding sequence in order to modulate pre-mRNA processing efficiency and mRNA stability, respectively. All of the vectors have extensive multiple cloning sites, and versions are available for producing N- or C- terminal GFP fusions at each of the possible relative expression levels. In addition, the modular configuration of the vectors allows drug resistance cassettes and other components to be readily exchanged. In toto, these vectors should be useful additions to the toolkit available for molecular and genetic studies of Leishmania donovani.

Reverse transcriptase genes are highly abundant and transcriptionally active in marine plankton assemblages

Genes encoding reverse transcriptases (RTs) are found in most eukaryotes, often as a component of retrotransposons, as well as in retroviruses and in prokaryotic retroelements. We investigated the abundance, classification and transcriptional status of RTs based on Tara Oceans marine metagenomes and metatranscriptomes encompassing a wide organism size range. Our analyses revealed that RTs predominate large-size fraction metagenomes (>5 μm), where they reached a maximum of 13.5% of the total gene abundance. Metagenomic RTs were widely distributed across the phylogeny of known RTs, but many belonged to previously uncharacterized clades. Metatranscriptomic RTs showed distinct abundance patterns across samples compared with metagenomic RTs. The relative abundances of viral and bacterial RTs among identified RT sequences were higher in metatranscriptomes than in metagenomes and these sequences were detected in all metatranscriptome size fractions. Overall, these observations suggest an active proliferation of various RT-assisted elements, which could be involved in genome evolution or adaptive processes of plankton assemblage.

Species-specific antimonial sensitivity in Leishmania is driven by post-transcriptional regulation of AQP1

Leishmania is a digenetic protozoan parasite causing leishmaniasis in humans. The different clinical forms of leishmaniasis are caused by more than twenty species of Leishmania that are transmitted by nearly thirty species of phlebotomine sand flies. Pentavalent antimonials (such as Pentostam or Glucantime) are the first line drugs for treating leishmaniasis. Recent studies suggest that pentavalent antimony (Sb(V)) acts as a pro-drug, which is converted to the more active trivalent form (Sb(III)). However, sensitivity to trivalent antimony varies among different Leishmania species. In general, Leishmania species causing cutaneous leishmaniasis (CL) are more sensitive to Sb(III) than the species responsible for visceral leishmaniasis (VL). Leishmania aquaglyceroporin (AQP1) facilitates the adventitious passage of antimonite down a concentration gradient. In this study, we show that Leishmania species causing CL accumulate more antimonite, and therefore exhibit higher sensitivity to antimonials, than the species responsible for VL. This species-specific differential sensitivity to antimonite is directly proportional to the expression levels of AQP1 mRNA. We show that the stability of AQP1 mRNA in different Leishmania species is regulated by their respective 3’-untranslated regions. The differential regulation of AQP1 mRNA explains the distinct antimonial sensitivity of each species.

The degree of response to antimonial drugs varies widely between species and even among strains of the same species of the protozoan parasite Leishmania. However, the molecular mechanism(s) is unknown. In this study, we show that Leishmania aquaglyceroporin AQP1 drives this species-specific antimonial resistance. Aquaglyceroporins are channel proteins that facilitate the passage of small uncharged molecules, such as glycerol and water, across the biological membranes. AQP1 helps the parasite cope with the osmotic challenges it faces during its life cycle. Additionally, AQP1 is an adventitious facilitator of antimonite, the active form of pentavalent antimonial drugs. We show that AQP1 expression level is species-specific, and less AQP1 in visceral species compared to the cutaneous species results in increased resistance to antimonials. We also demonstrate that the 3’-untranslated regions (3’-UTR) of the AQP1 mRNA is a major determining factor of species-specific regulation of AQP1. Along with water homeostasis, aquaglyceroporins are also involved in directed cell migration. The variable levels of AQP1 in different Leishmania species may enable them to find their appropriate niches in vertebrate hosts and cope with the species-specific osmotic challenges during their life cycles.

Identification and analysis of ingi-related retroposons in the trypanosomatid genomes

Transposable elements (TE), defined as discrete pieces of DNA that can move from one site to another site in genomes, represent significant components of eukaryotic genomes, including trypanosomatids. Up to 5% of the trypanosomatid genome content is composed of retroposons of the ingi clade, further divided into subclades and subfamilies ranging from short extinct truncated elements (SIDER) to long active elements (ingi). Important differences in ingi-related retroposon content have been reported between trypanosomatid species. For instance, Leishmania spp. have expanded and recycled a whole SIDER family to fulfill an important biological pathway, i.e., regulation of gene expression, while trypanosome genomes are primarily composed of active elements. Here, we present an overview of the computational methods used to identify, annotate, and analyze ingi-related retroposons for providing a comprehensive picture of all these TE families in newly available trypanosomatid genome sequences.

Pr77 and L1TcRz: A dual system within the 5'-end of L1Tc retrotransposon, internal promoter and HDV-like ribozyme

The sequence corresponding to the first 77 nucleotides of the L1Tc and NARTc non-LTR retrotransposons from Trypanosoma cruzi is an internal promoter (Pr77) that generates abundant, although poorly translatable, un-spliced transcripts. It has been recently described that L1TcRz, an HDV-like ribozyme, resides within the 5'-end of the RNA from the L1Tc and NARTc retrotransposons. Remarkably, the same first 77 nucleotides of L1Tc/NARTc elements comprise both the Pr77 internal promoter and the HDV-like L1TcRz. The L1TcRz cleaves on the 5'-side of the +1 nucleotide of the L1Tc element insuring that the promoter and the ribozyme functions travel with the transposon during retrotransposition. The ribozyme activity would prevent the mobilization of upstream sequences and insure the individuality of the L1Tc/NARTc copies transcribed from associated tandems. The Pr77/L1TcRz sequence is also found in other trypanosomatid's non-LTR retrotransposons and degenerated retroposons. The possible conservation of the ribozyme activity in a widely degenerated retrotransposon, as the Leishmania SIDERs, could indicate that the presence of this element and the catalytic activity could play some favorable genetic regulation. The functional implications of the Pr77/L1TcRz dual system in the regulation of the L1Tc/NARTc retrotransposons and in the gene expression of trypanosomatids are also discussed in this paper.

Genome evolution in trypanosomatid parasites

A decade of genome sequencing has transformed our understanding of how trypanosomatid parasites have evolved and provided fresh impetus to explaining the origins of parasitism in the Kinetoplastida. In this review, I will consider the many ways in which genome sequences have influenced our view of genomic reduction in trypanosomatids; how species-specific genes, and the genomic domains they occupy, have illuminated the innovations in trypanosomatid genomes; and how comparative genomics has exposed the molecular mechanisms responsible for innovation and adaptation to a parasitic lifestyle.

Genome-wide stochastic adaptive DNA amplification at direct and inverted DNA repeats in the parasite Leishmania

The human parasite Leishmania uses adaptive gene rearrangements and amplification involving repeated sequences on a genome-wide scale as one strategy to adapt to a changing environment.

Gene amplification of specific loci has been described in all kingdoms of life. In the protozoan parasite Leishmania, the product of amplification is usually part of extrachromosomal circular or linear amplicons that are formed at the level of direct or inverted repeated sequences. A bioinformatics screen revealed that repeated sequences are widely distributed in the Leishmania genome and the repeats are chromosome-specific, conserved among species, and generally present in low copy number. Using sensitive PCR assays, we provide evidence that the Leishmania genome is continuously being rearranged at the level of these repeated sequences, which serve as a functional platform for constitutive and stochastic amplification (and deletion) of genomic segments in the population. This process is adaptive as the copy number of advantageous extrachromosomal circular or linear elements increases upon selective pressure and is reversible when selection is removed. We also provide mechanistic insights on the formation of circular and linear amplicons through RAD51 recombinase-dependent and -independent mechanisms, respectively. The whole genome of Leishmania is thus stochastically rearranged at the level of repeated sequences, and the selection of parasite subpopulations with changes in the copy number of specific loci is used as a strategy to respond to a changing environment.

Variations in the copy number of DNA segments account for a substantial amount of genome diversity of most organisms. DNA amplification, a contributor to copy number variation, can occur in response to various stresses or after altered growth conditions, leading to extensive and often reversible genetic variation. DNA amplification in the parasite Leishmania occurs outside the normal chromosomes and arises by DNA rearrangements involving homologous repeated sequences. We show here that such repeated sequences are widespread in the Leishmania genome and that most of the Leishmania genome is subject to stochastic gene rearrangements mediated by these low-copy repeat sequences. Thus, although cells in the population have a common core genome, many individual cells will differ from the rest of the population by carrying one or more distinct extrachromosomal amplicon. Upon selection with either drugs or culture conditions, a subpopulation can emerge where the amplicon copy number per cell increases, and this clone of cells can then expand to dominate the population. We propose that Leishmania uses adaptive gene amplification at a genome-wide scale as one strategy to adapt to a changing environment.

A short RNA stem-loop is necessary and sufficient for repression of gene expression during early logarithmic phase in trypanosomes

We have compared the transcriptomes of cultured procyclic Trypanosoma brucei cells in early and late logarithmic phases and found that ∼200 mRNAs were differentially regulated. In late log phase cells, the most upregulated mRNA encoded the nucleobase transporter NT8. The 3' untranslated region (UTR) of NT8 contains a short stem-loop cis-element that is necessary for the regulation of NT8 expression in response to external purine levels. When placed in the 3'-UTR of an unregulated transcript, the cis-element is sufficient to confer regulation in response to purines. To our knowledge, this is the first example of a discrete RNA element that can autonomously regulate gene expression in trypanosomes in response to an external factor and reveals an unprecedented purine-dependent signaling pathway that controls gene expression in eukaryotes.

The wide expansion of hepatitis delta virus-like ribozymes throughout trypanosomatid genomes is linked to the spreading of L1Tc/ingi clade mobile elements

Background

Hepatitis Delta Virus (HDV)-like ribozymes have recently been found in many mobile elements in which they take part in a mechanism that releases intermediate RNAs from cellular co-transcripts. L1Tc in Trypanosoma cruzi is one of the elements in which such a ribozyme is located. It lies in the so-called Pr77-hallmark, a conserved region shared by retrotransposons belonging to the trypanosomatid L1Tc/ingi clade. The wide distribution of the Pr77-hallmark detected in trypanosomatid retrotransposons renders the potential catalytic activity of these elements worthy of study: their distribution might contribute to host genetic regulation at the mRNA level. Indeed, in Leishmania spp, the pervasive presence of these HDV-like ribozyme-containing mobile elements in certain 3′-untranslated regions of protein-coding genes has been linked to mRNA downregulation.

Results

Intensive screening of publicly available trypanosomatid genomes, combined with manual folding analyses, allowed the isolation of putatively Pr77-hallmarks with HDV-like ribozyme activity. This work describes the conservation of an HDV-like ribozyme structure in the Pr77 sequence of retrotransposons in a wide range of trypanosomatids, the catalytic function of which is maintained in the majority.

These results are consistent with the previously suggested common phylogenetic origin of the elements that belong to this clade, although in some cases loss of functionality appears to have occurred and/or perhaps molecular domestication by the host.

Conclusions

These HDV-like ribozymes are widely distributed within retrotransposons across trypanosomatid genomes. This type of ribozyme was once thought to be rare in nature, but in fact it would seem to be abundant in trypanosomatid transcripts. It can even form part of the pool of mRNA 3′-untranslated regions, particularly in Leishmania spp. Its putative regulatory role in host genetic expression is discussed.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-340) contains supplementary material, which is available to authorized users.

Metabolic reprogramming during purine stress in the protozoan pathogen Leishmania donovani

The ability of Leishmania to survive in their insect or mammalian host is dependent upon an ability to sense and adapt to changes in the microenvironment. However, little is known about the molecular mechanisms underlying the parasite response to environmental changes, such as nutrient availability. To elucidate nutrient stress response pathways in Leishmania donovani, we have used purine starvation as the paradigm. The salvage of purines from the host milieu is obligatory for parasite replication; nevertheless, purine-starved parasites can persist in culture without supplementary purine for over three months, indicating that the response to purine starvation is robust and engenders parasite survival under conditions of extreme scarcity. To understand metabolic reprogramming during purine starvation we have employed global approaches. Whole proteome comparisons between purine-starved and purine-replete parasites over a 6–48 h span have revealed a temporal and coordinated response to purine starvation. Purine transporters and enzymes involved in acquisition at the cell surface are upregulated within a few hours of purine removal from the media, while other key purine salvage components are upregulated later in the time-course and more modestly. After 48 h, the proteome of purine-starved parasites is extensively remodeled and adaptations to purine stress appear tailored to deal with both purine deprivation and general stress. To probe the molecular mechanisms affecting proteome remodeling in response to purine starvation, comparative RNA-seq analyses, qRT-PCR, and luciferase reporter assays were performed on purine-starved versus purine-replete parasites. While the regulation of a minority of proteins tracked with changes at the mRNA level, for many regulated proteins it appears that proteome remodeling during purine stress occurs primarily via translational and/or post-translational mechanisms.

Leishmania, the cause of a deadly spectrum of diseases in humans, surmounts a number of environmental challenges, including changes in the availability of salvageable nutrients, to successfully colonize its host. Adaptation to environmental stress is clearly of significance in parasite biology, but the underlying mechanisms are not well understood. To simulate the response to periodic nutrient scarcity in vivo, we have induced purine starvation in vitro. Purines are essential for growth and viability, and serve as the major energy currency of cells. Leishmania cannot synthesize purines and must salvage them from the surroundings. Extracellular purine depletion in culture induces a robust survival response in Leishmania, whereby growth arrests, but parasites persist for months. To profile the events that enable endurance of purine starvation, we used shotgun proteomics. Our data suggest that purine starvation induces extensive proteome remodeling, tailored to enhance purine capture and recycling, reduce energy expenditures, and maintain viability of the metabolically active, non-dividing population. Through global and targeted approaches, we reveal that proteome remodeling is multifaceted, and occurs through an array of responses at the mRNA, translational, and post-translational level. Our data provide one of the most inclusive views of adaptation to microenvironmental stress in Leishmania.

The streamlined genome of Phytomonas spp. relative to human pathogenic kinetoplastids reveals a parasite tailored for plants

Members of the family Trypanosomatidae infect many organisms, including animals, plants and humans. Plant-infecting trypanosomes are grouped under the single genus Phytomonas, failing to reflect the wide biological and pathological diversity of these protists. While some Phytomonas spp. multiply in the latex of plants, or in fruit or seeds without apparent pathogenicity, others colonize the phloem sap and afflict plants of substantial economic value, including the coffee tree, coconut and oil palms. Plant trypanosomes have not been studied extensively at the genome level, a major gap in understanding and controlling pathogenesis. We describe the genome sequences of two plant trypanosomatids, one pathogenic isolate from a Guianan coconut and one non-symptomatic isolate from Euphorbia collected in France. Although these parasites have extremely distinct pathogenic impacts, very few genes are unique to either, with the vast majority of genes shared by both isolates. Significantly, both Phytomonas spp. genomes consist essentially of single copy genes for the bulk of their metabolic enzymes, whereas other trypanosomatids e.g. Leishmania and Trypanosoma possess multiple paralogous genes or families. Indeed, comparison with other trypanosomatid genomes revealed a highly streamlined genome, encoding for a minimized metabolic system while conserving the major pathways, and with retention of a full complement of endomembrane organelles, but with no evidence for functional complexity. Identification of the metabolic genes of Phytomonas provides opportunities for establishing in vitro culturing of these fastidious parasites and new tools for the control of agricultural plant disease.

Some plant trypanosomes, single-celled organisms living in phloem sap, are responsible for important palm diseases, inducing frequent expensive and toxic insecticide treatments against their insect vectors. Other trypanosomes multiply in latex tubes without detriment to their host. Despite the wide range of behaviors and impacts, these trypanosomes have been rather unceremoniously lumped into a single genus: Phytomonas. A battery of molecular probes has been used for their characterization but no clear phylogeny or classification has been established. We have sequenced the genomes of a pathogenic phloem-specific Phytomonas from a diseased South American coconut palm and a latex-specific isolate collected from an apparently healthy wild euphorb in the south of France. Upon comparison with each other and with human pathogenic trypanosomes, both Phytomonas revealed distinctive compact genomes, consisting essentially of single-copy genes, with the vast majority of genes shared by both isolates irrespective of their effect on the host. A strong cohort of enzymes in the sugar metabolism pathways was consistent with the nutritional environments found in plants. The genetic nuances may reveal the basis for the behavioral differences between these two unique plant parasites, and indicate the direction of our future studies in search of effective treatment of the crop disease parasites.

Dosage compensation via transposable element mediated rewiring of a regulatory network

Transposable elements (TEs) may contribute to evolutionary innovations through the rewiring of networks by supplying ready-to-use cis regulatory elements. Genes on the Drosophila X chromosome are coordinately regulated by the male specific lethal (MSL) complex to achieve dosage compensation in males. We show that the acquisition of dozens of MSL binding sites on evolutionarily new X chromosomes was facilitated by the independent co-option of a mutant helitron TE that attracts the MSL complex (TE domestication). The recently formed neo-X recruits helitrons that provide dozens of functional, but suboptimal, MSL binding sites, whereas the older XR chromosome has ceased acquisition and appears to have fine-tuned the binding affinities of more ancient elements for the MSL complex. Thus, TE-mediated rewiring of regulatory networks through domestication and amplification may be followed by fine-tuning of the cis-regulatory element supplied by the TE and erosion of nonfunctional regions.

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.

The structure and repertoire of small interfering RNAs in Leishmania (Viannia) braziliensis reveal diversification in the trypanosomatid RNAi pathway

Among trypanosomatid protozoa the mechanism of RNA interference (RNAi) has been investigated in Trypanosoma brucei and to a lesser extent in Leishmania braziliensis. Although these two parasitic organisms belong to the same family, they are evolutionarily distantly related raising questions about the conservation of the RNAi pathway. Here we carried out an in-depth analysis of small interfering RNAs (siRNAs) associated with L. braziliensis Argonaute1 (LbrAGO1). In contrast to T. brucei, Leishmania siRNAs are sensitive to 3' end oxidation, indicating the absence of blocking groups, and the Leishmania genome does not code for a HEN1 RNA 2'-O-methyltransferase, which modifies small RNA 3' ends. Consistent with this observation, ~20% of siRNA 3' ends carry non-templated uridines. Thus siRNA biogenesis, and most likely their metabolism, is different in these organisms. Similarly to T. brucei, putative mobile elements and repeats constitute the major Leishmania siRNA-producing loci and AGO1 ablation leads to accumulation of long transcripts derived from putative mobile elements. However, contrary to T. brucei, no siRNAs were detected from other genomic regions with the potential to form double-stranded RNA, namely sites of convergent transcription and inverted repeats. Thus, our results indicate that organism-specific diversification has occurred in the RNAi pathway during evolution of the trypanosomatid lineage.

Novel features of a PIWI-like protein homolog in the parasitic protozoan Leishmania

In contrast to nearly all eukaryotes, the Old World Leishmania species L. infantum and L. major lack the bona fide RNAi machinery genes. Interestingly, both Leishmania genomes code for an atypical Argonaute-like protein that possesses a PIWI domain but lacks the PAZ domain found in Argonautes from RNAi proficient organisms. Using sub-cellular fractionation and confocal fluorescence microscopy, we show that unlike other eukaryotes, the PIWI-like protein is mainly localized in the single mitochondrion in Leishmania. To predict PIWI function, we generated a knockout mutant for the PIWI gene in both L. infantum (Lin) and L. major species by double-targeted gene replacement. Depletion of PIWI has no effect on the viability of insect promastigote forms but leads to an important growth defect of the mammalian amastigote lifestage in vitro and significantly delays disease pathology in mice, consistent with a higher expression of the PIWI transcript in amastigotes. Moreover, amastigotes lacking PIWI display a higher sensitivity to apoptosis inducing agents than wild type parasites, suggesting that PIWI may be a sensor for apoptotic stimuli. Furthermore, a whole-genome DNA microarray analysis revealed that loss of LinPIWI in Leishmania amastigotes affects mostly the expression of specific subsets of developmentally regulated genes. Several transcripts encoding surface and membrane-bound proteins were found downregulated in the LinPIWI^(−/−) mutant whereas all histone transcripts were upregulated in the null mutant, supporting the possibility that PIWI plays a direct or indirect role in the stability of these transcripts. Although our data suggest that PIWI is not involved in the biogenesis or the stability of small noncoding RNAs, additional studies are required to gain further insights into the role of this protein on RNA regulation and amastigote development in Leishmania.

Transposable elements: an abundant and natural source of regulatory sequences for host genes

The fact that transposable elements (TEs) can influence host gene expression was first recognized more than 50 years ago. However, since that time, TEs have been widely regarded as harmful genetic parasites-selfish elements that are rarely co-opted by the genome to serve a beneficial role. Here, we survey recent findings that relate to TE impact on host genes and remind the reader that TEs, in contrast to other noncoding parts of the genome, are uniquely suited to gene regulatory functions. We review recent studies that demonstrate the role of TEs in establishing and rewiring gene regulatory networks and discuss the overall ubiquity of exaptation. We suggest that although individuals within a population can be harmed by the deleterious effects of new TE insertions, the presence of TE sequences in a genome is of overall benefit to the population.

Gene expression regulation in trypanosomatids

Trypanosomatids are protozoan micro-organisms that cause serious health problems in humans and domestic animals. In addition to their medical relevance, these pathogens have novel biological structures and processes. From nuclear DNA transcription to mRNA translation, trypanosomes use unusual mechanisms to control gene expression. For example, transcription by RNAPII (RNA polymerase II) is polycistronic, and only a few transcription initiation sites have been identified so far. The sequences present in the polycistronic units code for proteins having unrelated functions, that is, not involved in a similar metabolic pathway. Owing to these biological constraints, these micro-organisms regulate gene expression mostly by post-transcriptional events. Consequently, the function of proteins that recognize RNA elements preferentially at the 3' UTR (untranslated region) of transcripts is central. It was recently shown that mRNP (messenger ribonucleoprotein) complexes are organized within post-transcriptional operons to co-ordinately regulate gene expression of functionally linked transcripts. In the present chapter we will focus on particular characteristics of gene expression in the so-called TriTryp parasites: Trypanosoma cruzi, Trypanosoma brucei and Leishmania major.

A 43-nucleotide U-rich element in 3'-untranslated region of large number of Trypanosoma cruzi transcripts is important for mRNA abundance in intracellular amastigotes

Trypanosoma cruzi, the agent of Chagas disease, does not seem to control gene expression through regulation of transcription initiation and makes use of post-transcriptional mechanisms. We report here a 43-nt U-rich RNA element located in the 3'-untranslated region (3'-UTR) of a large number of T. cruzi mRNAs that is important for mRNA abundance in the intracellular amastigote stage of the parasite. Whole genome scan analysis, differential display RT-PCR, Northern blot, and RT-PCR analyses were used to determine the transcript levels of more than 900 U-rich-containing mRNAs of large gene families as well as single and low copy number genes. Our results indicate that the 43-nt U-rich mRNA element is preferentially present in amastigotes. The cis-element of a protein kinase 3'-UTR but not its mutated version promoted the expression of the green fluorescent protein reporter gene in amastigotes. The regulatory cis-element, but not its mutated version, was also shown to interact with the trypanosome-specific RNA-binding protein (RBP) TcUBP1 but not with other related RBPs. Co-immunoprecipitation experiments of TcUBP1-containing ribonucleoprotein complexes formed in vivo validated the interaction with representative endogenous RNAs having the element. These results suggest that this 43-nt U-rich element together with other yet unidentified sequences might be involved in the modulation of abundance and/or translation of subsets of transcripts in the amastigote stage.

Evolution of genome content: population dynamics of transposable elements in flies and humans

Recent research is starting to shed light on the factors that influence the population and evolutionary dynamics of transposable elements (TEs) and TE life cycles. Genomes differ sharply in the number of TE copies, in the level of TE activity, in the diversity of TE families and types, and in the proportion of old and young TEs. In this chapter, we focus on two well-studied genomes with strikingly different architectures, humans and Drosophila, which represent two extremes in terms of TE diversity and population dynamics. We argue that some of the answers might lie in (1) the larger population size and consequently more effective selection against new TE insertions due to ectopic recombination in flies compared to humans; and (2) in the faster rate of DNA loss in flies compared to humans leading to much faster removal of fixed TE copies from the fly genome.