Evolution of non-LTR retrotransposons in the trypanosomatid genomes: Leishmania major has lost the active elements

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.

Variability of the Pr77 sequence of L1Tc retrotransposon among six T. cruzi strains belonging to different discrete typing units (DTUs)

All trypanosomatid genomes are colonized by non-LTR retrotransposons which exhibit a highly conserved 77-nt sequence at their 5' ends, known as the Pr77-hallmark (Pr77). The wide distribution of Pr77 is expected to be related to the gene regulation processes in these organisms as it has promoter and HDV-like ribozyme activities at the DNA and RNA levels, respectively. The identification of Pr77 hallmark-bearing retrotransposons and the study of the associations of mobile elements with relevant genes have been analyzed in the genomes of six strains of Trypanosoma cruzi belonging to different discrete typing units (DTUs) and with different geographical origins and host/vectors. The genomes have been sequenced, assembled and annotated. BUSCO analyses indicated a good quality for the assemblies that were used in comparative analyses. The results show differences among the six genomes in the copy number of genes related to virulence processes, the abundance of retrotransposons bearing the Pr77 sequence and the presence of the Pr77 hallmarks not associated with retroelements. The analyses also show frequent associations of Pr77-bearing retrotransposons and single Pr77 hallmarks with genes coding for trans-sialidases, RHS, MASP or hypothetical proteins, showing variable proportion depending on the type of retroelement, gene class and parasite strain. These differences in the genomic distribution of active retroelements and other Pr77-containing elements have shaped the genome architecture of these six strains and might be contributing to the phenotypic variability existing among them.

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.

The Tritryps Comparative Repeatome: Insights on Repetitive Element Evolution in Trypanosomatid Pathogens

The major human pathogens Trypanosoma cruzi, Trypanosoma brucei, and Leishmania major are collectively known as the Tritryps. The initial comparative analysis of their genomes has uncovered that Tritryps share a great number of genes, but repetitive DNA seems to be extremely variable between them. However, the in-depth characterization of repetitive DNA in these pathogens has been in part neglected, mainly due to the well-known technical challenges of studying repetitive sequences from de novo assemblies using short reads. Here, we compared the repetitive DNA repertories between the Tritryps genomes using genome-wide, low-coverage Illumina sequencing coupled to RepeatExplorer analysis. Our work demonstrates that this extensively implemented approach for studying higher eukaryote repeatomes is also useful for protozoan parasites like trypanosomatids, as we recovered previously observed differences in the presence and amount of repetitive DNA families. Additionally, our estimations of repetitive DNA abundance were comparable to those obtained from enhanced-quality assemblies using longer reads. Importantly, our methodology allowed us to describe a previously undescribed transposable element in Leishmania major (TATE element), highlighting its potential to accurately recover distinctive features from poorly characterized repeatomes. Together, our results support the application of this low-cost, low-coverage sequencing approach for the extensive characterization of repetitive DNA evolutionary dynamics in trypanosomatid and other protozoan genomes.

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.

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.

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.

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.

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.

Comparative in-silico genome analysis of Leishmania (Leishmania) donovani: A step towards its species specificity

Comparative genome analysis of recently sequenced Leishmania (L.) donovani was unexplored so far. The present study deals with the complete scanning of L. (L.) donovani genome revealing its interspecies variations. 60 distinctly present genes in L. (L.) donovani were identified when the whole genome was compared with Leishmania (L.) infantum. Similarly 72, 159, and 265 species specific genes were identified in L. (L.) donovani when compared to Leishmania (L.) major, Leishmania (L.) mexicana and Leishmania (Viannia) braziliensis respectively. The cross comparison of L. (L.) donovani in parallel with the other sequenced species of leishmanial led to the identification of 55 genes which are highly specific and expressed exclusively in L. (L.) donovani. We found mainly the discrepancies of surface proteins such as amastins, proteases, and peptidases. Also 415 repeat containing proteins in L. (L.) donovani and their differential distribution in other leishmanial species were identified which might have a potential role during pathogenesis. The genes identified can be evaluated as drug targets for anti-leishmanial treatment, exploring the scope for extensive future investigations.

•

Comparative genome analysis identifies 55 species specific L. (L.) donovani genes.

•

Discrepancies of surface proteins such as amastins, proteases, and peptidases are identified in L. (L.) donovani.

•

Apical Membrane Antigen (AMA1) might be a novel factor which helps L. (L.) donovani invasion.

•

Novel A2 and amastin genes in L. (L.) donovani genome are identified.

•

Our study identifies differential gene distribution in L. (L.) donovani with respect to other leishmanial species.

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.

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.

TSIDER1, a short and non-autonomous Salivarian trypanosome-specific retroposon related to the ingi6 subclade

Retroposons of the ingi clade are the most abundant transposable elements identified in the trypanosomatid genomes. Some are long autonomous elements (ingi, L1Tc) while others, such as RIME and NARTc, are short non-coding elements that parasitize the retrotransposition machinery of the active autonomous ones for their own mobilization. Here, we identified a new family of short non-autonomous retroposons of the ingi clade, called TSIDER1, which are present in the genome of Salivarian (African) trypanosomes, Trypanosoma brucei, T. congolense and T. vivax, but absent in the T. cruzi and Leishmania spp. genomes and, as such, TSIDER1 is the only retroposon subfamily conserved at the nucleotide level between African trypanosome species. We identified three TvSIDER1 families within the genome of T. vivax and the high level of sequence conservation within the TvSIDER1a and TvSIDER1b groups suggests that they are still active. We propose that TvSIDER1a/b elements are using the Tvingi retrotransposition machinery, as they are preceded by the same conserved pattern characteristic of the ingi6 subclade, which corresponds to the retroposon-encoded endonuclease binding site. In contrast, TcoSIDER1, TbSIDER1 and TvSIDER1c are too divergent to be considered as active retroposons. The relatively low number of SIDER elements identified in the T. congolense (70 copies), T. vivax (32 copies) and T. brucei (22 copies) genomes confirms that trypanosomes have not expanded short transposable elements, which is in contrast to Leishmania spp. (∼2000 copies), where SIDER play a role in the regulation of gene expression.

Centromere-associated topoisomerase activity in bloodstream form Trypanosoma brucei

Topoisomerase-II accumulates at centromeres during prometaphase, where it resolves the DNA catenations that represent the last link between sister chromatids. Previously, using approaches including etoposide-mediated topoisomerase-II cleavage, we mapped centromeric domains in trypanosomes, early branching eukaryotes in which chromosome segregation is poorly understood. Here, we show that in bloodstream form Trypanosoma brucei, RNAi-mediated depletion of topoisomerase-IIα, but not topoisomerase-IIβ, results in the abolition of centromere-localized activity and is lethal. Both phenotypes can be rescued by expression of the corresponding enzyme from T. cruzi. Therefore, processes which govern centromere-specific topoisomerase-II accumulation/activation have been functionally conserved within trypanosomes, despite the long evolutionary separation of these species and differences in centromeric DNA organization. The variable carboxyl terminal region of topoisomerase-II has a major role in regulating biological function. We therefore generated T. brucei lines expressing T. cruzi topoisomerase-II truncated at the carboxyl terminus and examined activity at centromeres after the RNAi-mediated depletion of the endogenous enzyme. A region necessary for nuclear localization was delineated to six residues. In other organisms, sumoylation of topoisomerase-II has been shown to be necessary for regulated chromosome segregation. Evidence that we present here suggests that sumoylation of the T. brucei enzyme is not required for centromere-specific cleavage activity.

Rapid decay of unstable Leishmania mRNAs bearing a conserved retroposon signature 3'-UTR motif is initiated by a site-specific endonucleolytic cleavage without prior deadenylation

We have previously shown that the Leishmania genome possess two widespread families of extinct retroposons termed Short Interspersed DEgenerated Retroposons (SIDER1/2) that play a role in post-transcriptional regulation. Moreover, we have demonstrated that SIDER2 retroposons promote mRNA degradation. Here we provide new insights into the mechanism by which unstable Leishmania mRNAs harboring a SIDER2 retroposon in their 3′-untranslated region are degraded. We show that, unlike most eukaryotic transcripts, SIDER2-bearing mRNAs do not undergo poly(A) tail shortening prior to rapid turnover, but instead, they are targeted for degradation by a site-specific endonucleolytic cleavage. The main cleavage site was mapped in two randomly selected SIDER2-containing mRNAs in vivo between an AU dinucleotide at the 5′-end of the second 79-nt signature (signature II), which represents the most conserved sequence amongst SIDER2 retroposons. Deletion of signature II abolished endonucleolytic cleavage and deadenylation-independent decay and increased mRNA stability. Interestingly, we show that overexpression of SIDER2 anti-sense RNA can increase sense transcript abundance and stability, and that complementarity to the cleavage region is required for protecting SIDER2-containing transcripts from degradation. These results establish a new paradigm for how unstable mRNAs are degraded in Leishmania and could serve as the basis for a better understanding of mRNA decay pathways in general.

Trypanosomatid genomes contain several subfamilies of ingi-related retroposons

Retroposons are ubiquitous transposable elements found in the genomes of most eukaryotes, including trypanosomatids. The African and American trypanosomes (Trypanosoma brucei and Trypanosoma cruzi) contain long autonomous retroposons of the ingi clade (Tbingi and L1Tc, respectively) and short nonautonomous truncated versions (TbRIME and NARTc, respectively), as well as degenerate ingi-related retroposons devoid of coding capacity (DIREs). In contrast, Leishmania major contains only remnants of extinct retroposons (LmDIREs) and of short nonautonomous heterogeneous elements (LmSIDERs). We extend this comparative and evolutionary analysis of retroposons to the genomes of two other African trypanosomes (Trypanosoma congolense and Trypanosoma vivax) and another Leishmania sp. (Leishmania braziliensis). Three new potentially functional retroposons of the ingi clade have been identified: Tvingi in T. vivax and Tcoingi and L1Tco in T. congolense. T. congolense is the first trypanosomatid containing two classes of potentially active retroposons of the ingi clade. We analyzed sequences located upstream of these new long autonomous ingi-related elements, which code for the recognition site of the retroposon-encoded endonuclease. The closely related Tcoingi and Tvingi elements show the same conserved pattern, indicating that the Tcoingi- and Tvingi-encoded endonucleases share site specificity. Similarly, the conserved pattern previously identified upstream of L1Tc has also been detected at the same relative position upstream of L1Tco elements. A phylogenetic analysis of all ingi-related retroposons identified so far, including DIREs, clearly shows that several distinct subfamilies have emerged and coexisted, though in the course of trypanosomatid evolution, only a few have been maintained as active elements in modern trypanosomatid (sub)species.

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.

The SIDER2 elements, interspersed repeated sequences that populate the Leishmania genomes, constitute subfamilies showing chromosomal proximity relationship

Background

Protozoan parasites of the genus Leishmania are causative agents of a diverse spectrum of human diseases collectively known as leishmaniasis. These eukaryotic pathogens that diverged early from the main eukaryotic lineage possess a number of unusual genomic, molecular and biochemical features. The completion of the genome projects for three Leishmania species has generated invaluable information enabling a direct analysis of genome structure and organization.

Results

By using DNA macroarrays, made with Leishmania infantum genomic clones and hybridized with total DNA from the parasite, we identified a clone containing a repeated sequence. An analysis of the recently completed genome sequence of L. infantum, using this repeated sequence as bait, led to the identification of a new class of repeated elements that are interspersed along the different L. infantum chromosomes. These elements turned out to be homologues of SIDER2 sequences, which were recently identified in the Leishmania major genome; thus, we adopted this nomenclature for the Leishmania elements described herein. Since SIDER2 elements are very heterogeneous in sequence, their precise identification is rather laborious. We have characterized 54 LiSIDER2 elements in chromosome 32 and 27 ones in chromosome 20. The mean size for these elements is 550 bp and their sequence is G+C rich (mean value of 66.5%). On the basis of sequence similarity, these elements can be grouped in subfamilies that show a remarkable relationship of proximity, i.e. SIDER2s of a given subfamily locate close in a chromosomal region without intercalating elements. For comparative purposes, we have identified the SIDER2 elements existing in L. major and Leishmania braziliensis chromosomes 32. While SIDER2 elements are highly conserved both in number and location between L. infantum and L. major, no such conservation exists when comparing with SIDER2s in L. braziliensis chromosome 32.

Conclusion

SIDER2 elements constitute a relevant piece in the Leishmania genome organization. Sequence characteristics, genomic distribution and evolutionarily conservation of SIDER2s are suggestive of relevant functions for these elements in Leishmania. Apart from a proved involvement in post-trancriptional mechanisms of gene regulation, SIDER2 elements could be involved in DNA amplification processes and, perhaps, in chromosome segregation as centromeric sequences.

Genomic organization and transcription analysis of the 195-bp satellite DNA in Trypanosoma cruzi

The 195-bp satellite DNA is the most abundant Trypanosoma cruzi repetitive sequence. Here we show by RNA blotting and RT-PCR that 195 SAT is intensely transcribed. We observed a positive correlation between the level of satellite RNA and the abundance of the satellite copies in the genome of T. cruzi strains and that the satellite expression is not developmentally regulated. By analyzing CL Brener individual reads, we estimated that 195 SAT corresponds to approximately 5% of the CL Brener genome. 195 SAT elements were found in only 37 annotated contigs, indicating that a large number of satellite copies were not incorporated into the assembled data. The assembled satellite units are distributed in non-syntenic regions with Trypanosoma brucei and Leishmania major genomes, enriched with surface proteins, retroelements, RHS and hypothetical proteins. Satellite repeats were not observed in annotated subtelomeric regions. We report that 12 satellite sequences are truncated by the retroelement VIPER.

Role of transposable elements in trypanosomatids

Transposable elements constitute 2-5% of the genome content in trypanosomatid parasites. Some of them are involved in critical cellular functions, such as the regulation of gene expression in Leishmania spp. In this review, we highlight the remarkable role extinct transposable elements can play as the source of potential new functions.

Telomeric co-localization of the modified base J and contingency genes in the protozoan parasite Trypanosoma cruzi

Base J or beta-d-glucosylhydroxymethyluracil is a modification of thymine residues within the genome of kinetoplastid parasites. In organisms known to contain the modified base, J is located mainly within the telomeric repeats. However, in Trypanosoma brucei, a small fraction of J is also located within the silent subtelomeric variant surface glycoprotein (VSG) gene expression sites, but not in the active expression site, suggesting a role for J in regulating telomeric genes involved in pathogenesis. With the identification of surface glycoprotein genes adjacent to telomeres in the South American Trypanosome, Trypanosoma cruzi, we became interested in the telomeric distribution of base J. Analysis of J and telomeric repeat sequences by J immunoblots and Southern blots following DNA digestion, reveals approximately 25% of J outside the telomeric repeat sequences. Moreover, the analysis of DNA sequences immunoprecipitated with J antiserum, localized J within subtelomeric regions rich in life-stage-specific surface glycoprotein genes involved in pathogenesis. Interestingly, the pattern of J within these regions is developmentally regulated. These studies provide a framework to characterize the role of base J in the regulation of telomeric gene expression/diversity in T. cruzi.

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

Trypanosomatids are unicellular protists that include the human pathogens Leishmania spp. (leishmaniasis), Trypanosoma brucei (sleeping sickness), and Trypanosoma cruzi (Chagas disease). Analysis of their recently completed genomes confirmed the presence of non-long-terminal repeat retrotransposons, also called retroposons. Using the 79-bp signature sequence common to all trypanosomatid retroposons as bait, we identified in the Leishmania major genome two new large families of small elements--LmSIDER1 (785 copies) and LmSIDER2 (1,073 copies)--that fulfill all the characteristics of extinct trypanosomatid retroposons. LmSIDERs are approximately 70 times more abundant in L. major compared to T. brucei and are found almost exclusively within the 3'-untranslated regions (3'UTRs) of L. major mRNAs. We provide experimental evidence that LmSIDER2 act as mRNA instability elements and that LmSIDER2-containing mRNAs are generally expressed at lower levels compared to the non-LmSIDER2 mRNAs. The considerable expansion of LmSIDERs within 3'UTRs in an organism lacking transcriptional control and their role in regulating mRNA stability indicate that Leishmania have probably recycled these short retroposons to globally modulate the expression of a number of genes. To our knowledge, this is the first example in eukaryotes of the domestication and expansion of a family of mobile elements that have evolved to fulfill a critical cellular function.

Comparative genomics: from genotype to disease phenotype in the leishmaniases

Recent progress in sequencing the genomes of several Leishmania species, causative agents of cutaneous, mucocutaneous and visceral leishmaniasis, is revealing unusual features of potential relevance to parasite virulence and pathogenesis in the host. While the genomes of Leishmania major, Leishmania braziliensis and Leishmania infantum are highly similar in content and organisation, species-specific genes and mechanisms distinguish one from another. In particular, the presence of retrotransposons and the components of a putative RNA interference machinery in L. braziliensis suggest the potential for both greater diversity and more tractable experimentation in this Leishmania Viannia species.

Comparative genomic analysis of three Leishmania species that cause diverse human disease

Leishmania parasites cause a broad spectrum of clinical disease. Here we report the sequencing of the genomes of two species of Leishmania: Leishmania infantum and Leishmania braziliensis. The comparison of these sequences with the published genome of Leishmania major reveals marked conservation of synteny and identifies only approximately 200 genes with a differential distribution between the three species. L. braziliensis, contrary to Leishmania species examined so far, possesses components of a putative RNA-mediated interference pathway, telomere-associated transposable elements and spliced leader-associated SLACS retrotransposons. We show that pseudogene formation and gene loss are the principal forces shaping the different genomes. Genes that are differentially distributed between the species encode proteins implicated in host-pathogen interactions and parasite survival in the macrophage.

The L1Tc non-LTR retrotransposon of Trypanosoma cruzi contains an internal RNA-pol II-dependent promoter that strongly activates gene transcription and generates unspliced transcripts

L1Tc is the best represented autonomous LINE of the Trypanosoma cruzi genome, throughout which several functional copies may exist. In this study, we show that the first 77 bp of L1Tc (Pr77) (also present in the T. cruzi non-autonomous retrotransposon NARTc, in the Trypanosoma brucei RIME/ingi elements, and in the T. cruzi, T. brucei and Leishmania major degenerate L1Tc/ingi-related elements [DIREs]) behave as a promoter element that activates gene transcription. The transcription rate promoted by Pr77 is 10–14-fold higher than that mediated by sequences located upstream from the T. cruzi tandemly repeated genes KMP11 and the GAPDH. The Pr77 promoter-derived mRNAs initiate at nucleotide +1 of L1Tc, are unspliced and translated. L1Tc transcripts show a moderate half life and are RNA pol II dependent. The presence of an internal promoter at the 5′ end of L1Tc favors the production of full-length L1Tc RNAs and reinforces the hypothesis that this mobile element may be naturally autonomous in its transposition.