Recombinant SINEs are formed at high frequency during induced retrotransposition in vivo

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.

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.

Possible role played by the SINE2 element in gene regulation, as demonstrated by differential processing and polyadenylation in avirulent strains of E. histolytica

Entamoeba histolytica represents a useful model in parasitic organisms due to its complex genomic organization and survival mechanisms. To counteract pathogenic organisms, it is necessary to characterize their molecular biology to design new strategies to combat them. In this report, we investigated a less-known genetic element, short interspersed nuclear element 2 (SINE2), that is present in this ameba and is highly transcribed and polyadenylated. In this study, we show that in two different nonvirulent strains of E. histolytica, SINE2 is differentially processed into two transcript fragments, that is, a full-length 560-nt fragment and a shorter 393-nt fragment bearing an approximately 18-nt polyadenylation tail. Sequence analysis of the SINE2 transcript showed that a Musashi-like protein may bind to it. Also, two putative Musashi-like sequences were identified on the transcript. Semiquantitative expression analysis of the two Musashi-like proteins identified in the E. histolytica genome (XP_648918 and XP_649094) showed that XP_64094 is overexpressed in the nonvirulent strains tested. The information available in the literature and the results presented in this report indicate that SINE2 may affect other genes, as observed with the epigenetic silencing of the G3 strain, by an antisense mechanism or via RNA-protein interactions that may ultimately be involved in the phenotype of nonvirulent strains of E. histolytica.

The Cassandra retrotransposon landscape in sugar beet (Beta vulgaris) and related Amaranthaceae: recombination and re-shuffling lead to a high structural variability

BACKGROUND AND AIMS

Plant genomes contain many retrotransposons and their derivatives, which are subject to rapid sequence turnover. As non-autonomous retrotransposons do not encode any proteins, they experience reduced selective constraints leading to their diversification into multiple families, usually limited to a few closely related species. In contrast, the non-coding Cassandra terminal repeat retrotransposons in miniature (TRIMs) are widespread in many plants. Their hallmark is a conserved 5S rDNA-derived promoter in their long terminal repeats (LTRs). As sugar beet (Beta vulgaris) has a well-described LTR retrotransposon landscape, we aim to characterize TRIMs in beet and related genomes.

METHODS

We identified Cassandra retrotransposons in the sugar beet reference genome and characterized their structural relationships. Genomic organization, chromosomal localization, and distribution of Cassandra-TRIMs across the Amaranthaceae were verified by Southern and fluorescent in situ hybridization.

KEY RESULTS

All 638 Cassandra sequences in the sugar beet genome contain conserved LTRs and thus constitute a single family. Nevertheless, variable internal regions required a subdivision into two Cassandra subfamilies within B. vulgaris. The related Chenopodium quinoa harbours a third subfamily. These subfamilies vary in their distribution within Amaranthaceae genomes, their insertion times and the degree of silencing by small RNAs. Cassandra retrotransposons gave rise to many structural variants, such as solo LTRs or tandemly arranged Cassandra retrotransposons. These Cassandra derivatives point to an interplay of template switch and recombination processes - mechanisms that likely caused Cassandra's subfamily formation and diversification.

CONCLUSIONS

We traced the evolution of Cassandra in the Amaranthaceae and detected a considerable variability within the short internal regions, whereas the LTRs are strongly conserved in sequence and length. Presumably these hallmarks make Cassandra a prime target for unequal recombination, resulting in the observed structural diversity, an example of the impact of LTR-mediated evolutionary mechanisms on the host genome.

Divergence of 3' ends as a driver of short interspersed nuclear element (SINE) evolution in the Salicaceae

Short interspersed nuclear elements (SINEs) are small, non-autonomous and heterogeneous retrotransposons that are widespread in plants. To explore the amplification dynamics and evolutionary history of SINE populations in representative deciduous tree species, we analyzed the genomes of the six following Salicaceae species: Populus deltoides, Populus euphratica, Populus tremula, Populus tremuloides, Populus trichocarpa, and Salix purpurea. We identified 11 Salicaceae SINE families (SaliS-I to SaliS-XI), comprising 27 077 full-length copies. Most of these families harbor segmental similarities, providing evidence for SINE emergence by reshuffling or heterodimerization. We observed two SINE groups, differing in phylogenetic distribution pattern, similarity and 3' end structure. These groups probably emerged during the 'salicoid duplication' (~65 million years ago) in the Salix-Populus progenitor and during the separation of the genus Salix (45-65 million years ago), respectively. In contrast to conserved 5' start motifs across species and SINE families, the 3' ends are highly variable in sequence and length. This extraordinary 3'-end variability results from mutations in the poly(A) tail, which were fixed by subsequent amplificational bursts. We show that the dissemination of newly evolved 3' ends is accomplished by a displacement of older motifs, leading to various 3'-end subpopulations within the SaliS families.

Becoming a Selfish Clan: Recombination Associated to Reverse-Transcription in LTR Retrotransposons

Transposable elements (TEs) are parasitic DNA bits capable of mobilization and mutagenesis, typically suppressed by host’s epigenetic silencing. Since the selfish DNA concept, it is appreciated that genomes are also molded by arms-races against natural TE inhabitants. However, our understanding of evolutionary processes shaping TEs adaptive populations is scarce. Here, we review the events of recombination associated to reverse-transcription in LTR retrotransposons, a process shuffling their genetic variants during replicative mobilization. Current evidence may suggest that recombinogenic retrotransposons could beneficially exploit host suppression, where clan behavior facilitates their speciation and diversification. Novel refinements to retrotransposons life-cycle and evolution models thus emerge.

The conserved 3' Angio-domain defines a superfamily of short interspersed nuclear elements (SINEs) in higher plants

Repetitive sequences are ubiquitous components of eukaryotic genomes affecting genome size and evolution as well as gene regulation. Among them, short interspersed nuclear elements (SINEs) are non-coding retrotransposons usually shorter than 1000 bp. They contain only few short conserved structural motifs, in particular an internal promoter derived from cellular RNAs and a mostly AT-rich 3' tail, whereas the remaining regions are highly variable. SINEs emerge and vanish during evolution, and often diversify into numerous families and subfamilies that are usually specific for only a limited number of species. In contrast, at the 3' end of multiple plant SINEs we detected the highly conserved 'Angio-domain'. This 37 bp segment defines the Angio-SINE superfamily, which encompasses 24 plant SINE families widely distributed across 13 orders within the plant kingdom. We retrieved 28 433 full-length Angio-SINE copies from genome assemblies of 46 plant species, frequently located in genes. Compensatory mutations in and adjacent to the Angio-domain imply selective restraints maintaining its RNA structure. Angio-SINE families share segmental sequence similarities, indicating a modular evolution with strong Angio-domain preservation. We suggest that the conserved domain contributes to the evolutionary success of Angio-SINEs through either structural interactions between SINE RNA and proteins increasing their transpositional efficiency, or by enhancing their accumulation in genes.

Evolutionary modes of emergence of short interspersed nuclear element (SINE) families in grasses

Short interspersed nuclear elements (SINEs) are non-autonomous transposable elements which are propagated by retrotransposition and constitute an inherent part of the genome of most eukaryotic species. Knowledge of heterogeneous and highly abundant SINEs is crucial for de novo (or improvement of) annotation of whole genome sequences. We scanned Poaceae genome sequences of six important cereals (Oryza sativa, Triticum aestivum, Hordeum vulgare, Panicum virgatum, Sorghum bicolor, Zea mays) and Brachypodium distachyon to examine the diversity and evolution of SINE populations. We comparatively analyzed the structural features, distribution, evolutionary relation and abundance of 32 SINE families and subfamilies within grasses, comprising 11 052 individual copies. The investigation of activity profiles within the Poaceae provides insights into their species-specific diversification and amplification. We found that Poaceae SINEs (PoaS) fall into two length categories: simple SINEs of up to 180 bp and dimeric SINEs larger than 240 bp. Detailed analysis at the nucleotide level revealed that multimerization of related and unrelated SINE copies is an important evolutionary mechanism of SINE formation. We conclude that PoaS families diversify by massive reshuffling between SINE families, likely caused by insertion of truncated copies, and provide a model for this evolutionary scenario. Twenty-eight of 32 PoaS families and subfamilies show significant conservation, in particular either in the 5' or 3' regions, across Poaceae species and share large sequence stretches with one or more other PoaS families.

Cytosine DNA methylation at promoter of non LTR retrotransposon and heat shock protein gene (HSP70) of Entamoeba histolytica and lack of correlation with transcription status

Non LTR retrotransposons (EhLINEs and EhSINEs) occupy 11% of the Entamoeba histolytica genome. Since promoter DNA methylation at cytosines has been correlated with transcriptional silencing of transposable elements in model organisms we checked whether this was the case in EhLINE1. We located promoter activity in a 841bp fragment at 5'-end of this element by luciferase reporter assay. From RNAseq and RT-PCR analyses we selected a transcriptionally active and silent copy to study cytosine DNA methylation of the promoter region by bisulfite sequencing. None of the cytosines were methylated in either copy. Further, we looked at methylation status of a few selected cytosines in all 5'-intact EhLINE1 copies by single nucleotide incorporation opposite cytosine in bisulfite-treated DNA, where dGTP would be incorporated if the cytosine was methylated. Again we did not find evidence of cytosine methylation, indicating that expression status of this element was not correlated with promoter DNA methylation. To test for any role of cytosine methylation in transcriptional regulation of the E. histolytica Hsp70 gene in which the promoter is fully methylated under normal growth conditions, we checked methylation status and found that the promoter remained fully methylated during heat-shock as well, although transcription was greatly enhanced by heat-shock, showing that cytosine methylation is not a repressive mark for EhHsp70. Our data present direct evidence that promoter methylation, a common mode of transposon silencing, is unlikely to be involved in transcriptional regulation of EhLINE1, and reinforce the conclusion that promoter DNA methylation may not be a major contributor to transcriptional regulation in E. histolytica.

Rare horizontal transmission does not hide long-term inheritance of SINE highly conserved domains in the metazoan evolution

Transposable elements (TEs) are self-replicating, mobile DNA sequences which constitute a significant fraction of eukaryotic genomes. They are generally considered selfish DNA, as their replication and random insertion may have deleterious effects on genome functionalities, although some beneficial effects and evolutionary potential have been recognized. Short interspersed elements (SINEs) are non-autonomous TEs with a modular structure: a small RNA-related head, a body, and a long interspersed element-related tail. Despite their high turnover rate and de novo emergence, the body may retain highly conserved domains (HCDs) shared among divergent SINE families: in metazoans, at least nine HCD-SINEs have been recognized. Data mining on public molecular databases allowed the retrieval of 16 new HCD-SINE families from cnidarian, molluscs, arthropods, and vertebrates. Tracking the ancestry of HCDs on the metazoan phylogeny revealed that some of them date back to the Radiata–Bilateria split. Moreover, phylogenetic and age versus divergence analyses of the most ancient HCDs suggested that long-term vertical inheritance is the rule, with few horizontal transfer events. We suggest that the evolutionary conservation of HCDs may be linked to their potential to serve as recombination hotspots. This indirectly affects host genomes by maintaining active and diverse SINE lineages, whose insertions may impact (either positively or negatively) on the evolution of the genome.

Evolutionary genomics and population structure of Entamoeba histolytica

Amoebiasis caused by the gastrointestinal parasite Entamoeba histolytica has diverse disease outcomes. Study of genome and evolution of this fascinating parasite will help us to understand the basis of its virulence and explain why, when and how it causes diseases. In this review, we have summarized current knowledge regarding evolutionary genomics of E. histolytica and discussed their association with parasite phenotypes and its differential pathogenic behavior. How genetic diversity reveals parasite population structure has also been discussed. Queries concerning their evolution and population structure which were required to be addressed have also been highlighted. This significantly large amount of genomic data will improve our knowledge about this pathogenic species of Entamoeba.

A genomewide overexpression screen identifies genes involved in the phosphatidylinositol 3-kinase pathway in the human protozoan parasite Entamoeba histolytica

Entamoeba histolytica is a protozoan parasite that causes amoebic dysentery and liver abscess. E. histolytica relies on motility, phagocytosis, host cell adhesion, and proteolysis of extracellular matrix for virulence. In eukaryotic cells, these processes are mediated in part by phosphatidylinositol 3-kinase (PI3K) signaling. Thus, PI3K may be critical for virulence. We utilized a functional genomics approach to identify genes whose products may operate in the PI3K pathway in E. histolytica. We treated a population of trophozoites that were overexpressing genes from a cDNA library with a near-lethal dose of the PI3K inhibitor wortmannin. This screen was based on the rationale that survivors would be overexpressing gene products that directly or indirectly function in the PI3K pathway. We sequenced the overexpressed genes in survivors and identified a cDNA encoding a Rap GTPase, a protein previously shown to participate in the PI3K pathway. This supports the validity of our approach. Genes encoding a coactosin-like protein, EhCoactosin, and a serine-rich E. histolytica protein (SREHP) were also identified. Cells overexpressing EhCoactosin or SREHP were also less sensitive to a second PI3K inhibitor, LY294002. This corroborates the link between these proteins and PI3K. Finally, a mutant cell line with an increased level of phosphatidylinositol (3,4,5)-triphosphate, the product of PI3K activity, exhibited increased expression of SREHP and EhCoactosin. This further supports the functional connection between these proteins and PI3K in E. histolytica. To our knowledge, this is the first forward-genetics screen adapted to reveal genes participating in a signal transduction pathway in this pathogen.

Stress granule formation in Entamoeba histolytica: cross-talk between EhMLBP, EhRLE3 reverse transcriptase and polyubiquitinated proteins

The Entamoeba histolytica-methylated LINE-binding protein (EhMLBP) binds to methylated repetitive DNA and is a positive regulator of a reverse transcriptase of a long interspersed nucleotide element (LINE). This protein protects trophozoites against heat shock by reducing protein aggregation. The presence of EhMLBP and polyubiquitinated proteins in heat shock-induced protein aggregates raised the question whether these proteins interact. This assumption was confirmed by co-immunoprecipitation experiments: ubiquitinated proteins were detected in the perinuclear region of non-stressed E. histolytica trophozoites, whereas ubiquitinated proteins were detected in the perinuclear region and colocalized with EhMLBP in cytoplasmic granules in heat-shocked trophozoites. We also observed that overexpression of the reverse transcriptase of EhRLE3 induced the upregulation of EhMLBP expression and the formation of these EhMLBP-containing granules. Since (i) these EhMLBP-containing granules in the cytoplasm of heat-shocked E. histolytica trophozoites also contain polyubiquitinated proteins and poly(A)(+) mRNA and (ii) their formation is promoted by sodium arsenate, puromycin, and pateamine A and is inhibited by cycloheximide, we propose that these cytoplasmic EhMLBP-containing granules are stress granules. Our data also suggest that the formation of these granules is dependent upon EhMLBP and LINE.

Homologous Recombination Occurs in Entamoeba and Is Enhanced during Growth Stress and Stage Conversion

Homologous recombination (HR) has not been demonstrated in the parasitic protists Entamoeba histolytica or Entamoeba invadens, as no convenient method is available to measure it. However, HR must exist to ensure genome integrity, and possible genetic exchange, especially during stage conversion from trophozoite to cyst. Here we show the up regulation of mitotic and meiotic HR genes in Entamoeba during serum starvation, and encystation. To directly demonstrate HR we use a simple PCR-based method involving inverted repeats, which gives a reliable read out, as the recombination junctions can be determined by sequencing the amplicons. Using this read out, we demonstrate enhanced HR under growth stress in E. histolytica, and during encystation in E. invadens. We also demonstrate recombination between chromosomal inverted repeats. This is the first experimental demonstration of HR in Entamoeba and will help future investigations into this process, and to explore the possibility of meiosis in Entamoeba.

Conserved domains and SINE diversity during animal evolution

Eukaryotic genomes harbour a number of mobile genetic elements (MGEs); moving from one genomic location to another, they are known to impact on the host genome. Short interspersed elements (SINEs) are well-represented, non-autonomous retroelements and they are likely the most diversified MGEs. In some instances, sequence domains conserved across unrelated SINEs have been identified; remarkably, one of these, called Nin, has been conserved since the Radiata-Bilateria splitting. Here we report on two new domains: Inv, derived from Nin, identified in insects and in deuterostomes, and Pln, restricted to polyneopteran insects. The identification of Inv and Pln sequences allowed us to retrieve new SINEs, two in insects and one in a hemichordate. The diverse structural combination of the different domains in different SINE families, during metazoan evolution, offers a clearer view of SINE diversity and their frequent de novo emergence through module exchange, possibly underlying the high evolutionary success of SINEs.

Genomic distribution of SINEs in Entamoeba histolytica strains: implication for genotyping

BACKGROUND

The major clinical manifestations of Entamoeba histolytica infection include amebic colitis and liver abscess. However the majority of infections remain asymptomatic. Earlier reports have shown that some E. histolytica isolates are more virulent than others, suggesting that virulence may be linked to genotype. Here we have looked at the genomic distribution of the retrotransposable short interspersed nuclear elements EhSINE1 and EhSINE2. Due to their mobile nature, some EhSINE copies may occupy different genomic locations among isolates of E. histolytica possibly affecting adjacent gene expression; this variability in location can be exploited to differentiate strains.

RESULTS

We have looked for EhSINE1- and EhSINE2-occupied loci in the genome sequence of Entamoeba histolytica HM-1:IMSS and searched for homologous loci in other strains to determine the insertion status of these elements. A total of 393 EhSINE1 and 119 EhSINE2 loci were analyzed in the available sequenced strains (Rahman, DS4-868, HM1:CA, KU48, KU50, KU27 and MS96-3382. Seventeen loci (13 EhSINE1 and 4 EhSINE2) were identified where a EhSINE1/EhSINE2 sequence was missing from the corresponding locus of other strains. Most of these loci were unoccupied in more than one strain. Some of the loci were analyzed experimentally for SINE occupancy using DNA from strain Rahman. These data helped to correctly assemble the nucleotide sequence at three loci in Rahman. SINE occupancy was also checked at these three loci in 7 other axenically cultivated E. histolytica strains and 16 clinical isolates. Each locus gave a single, specific amplicon with the primer sets used, making this a suitable method for strain typing. Based on presence/absence of SINE and amplification with locus-specific primers, the 23 strains could be divided into eleven genotypes. The results obtained by our method correlated with the data from other typing methods. We also report a bioinformatic analysis of EhSINE2 copies.

CONCLUSIONS

Our results reveal several loci with extensive polymorphism of SINE occupancy among different strains of E. histolytica and prove the principle that the genomic distribution of SINEs is a valid method for typing of E. histolytica strains.

Genome-wide analysis reveals origin of transfer RNA genes from tRNA halves

Transfer RNAs (tRNAs) play an important role linking mitochondrial RNA and amino acids during protein biogenesis. Four types of tRNA genes have been identified in living organisms. However, the evolutionary origin of tRNAs remains largely unknown. In this article, we conduct a deep sequence analysis of diverse genomes that cover all three domains of life to unveil the evolutionary history of tRNA genes from tRNA halves. tRNA half homologs were detected in diverse organisms, and some of them were expressed in mouse tissues. Continuous tRNA genes have a conserved pattern similar to indels, which is, more closely flanking regions have higher single nucleotide substitution rates, whereas tRNA half homologs do not have this pattern. In addition, tRNAs tend to break into tRNA halves when tissues are incubated in vitro, the tendency of tRNA to break into tRNA halves may be a "side-effect" of tRNA genes evolving from tRNA halves. These results suggest that modern tRNAs originated from tRNA halves through a repeat element-mediated mechanism. These findings provide insight into the evolutionary origin of tRNA genes.