Evolution of tRNA-like sequences and genome variability

Protein-coding tRNA sequences?

Transfer RNAs (tRNAs) are ancient molecules likely predating the translation machinery. These extremely conserved RNA molecules transfer amino acids to the ribosome for the synthesis of proteins encoded by mRNAs, but canonical tRNAs are not protein-coding RNAs. Surprisely, when virtually translated, I observed that peptides derived from tRNA sequences match thousands of protein entries in databases. The analysis of these sequences indicates that the vast majority of these tRNA-derived proteins are annotated as small hypothetical peptides, likely arising from sequencing, prediction and/or annotation errors. But life often surpasses fiction. Importantly, tRNA-encoded amino acid domains were also found embedded in large functional proteins. Phylogenetic analysis of representative tRNA-derived protein domains may provide new insights into the origin, plasticity, and evolution of protein-coding genes.

Seryl-tRNA synthetase specificity for tRNASec in Bacterial Sec biosynthesis

tRNA synthetases are responsible for decoding the molecular information, from codons to amino acids. Seryl-tRNA synthetase (SerRS), besides the five isoacceptors of tRNA^Ser, recognizes tRNA^[Ser]Sec for the incorporation of selenocysteine (Sec, U) into selenoproteins. The selenocysteine synthesis pathway is known and is dependent on several protein-protein and protein-RNA interactions. Those interactions are not fully described, in particular, involving tRNA^[Ser]Sec and SerRS. Here we describe the molecular interactions between the Escherichia coli Seryl-tRNA synthetase (EcSerRS) and tRNA^[Ser]Sec in order to determine their specificity, selectivity and binding order, leading to tRNA aminoacylation. The dissociation constant of EcSerRS and tRNA^[Ser]Sec was determined as (126 ± 20) nM. We also demonstrate that EcSerRS binds initially to tRNA^[Ser]Sec in the presence of ATP for further recognition by E. coli selenocysteine synthetase (EcSelA) for Ser to Sec conversion. The proposed studies clarify the mechanism of tRNA^[Ser]Sec incorporation in Bacteria as well as of other domains of life.

tRNA Genes Affect Chromosome Structure and Function via Local Effects

The genome is packaged and organized in an ordered, nonrandom manner, and specific chromatin segments contact nuclear substructures to mediate this organization. tRNA genes (tDNAs) are binding sites for transcription factors and architectural proteins and are thought to play an important role in the organization of the genome. In this study, we investigate the roles of tDNAs in genomic organization and chromosome function by editing a chromosome so that it lacked any tDNAs. Surprisingly our analyses of this tDNA-less chromosome show that loss of tDNAs does not grossly affect chromatin architecture or chromosome tethering and mobility. However, loss of tDNAs affects local nucleosome positioning and the binding of SMC proteins at these loci. The absence of tDNAs also leads to changes in centromere clustering and a reduction in the frequency of long-range HML-HMR heterochromatin clustering with concomitant effects on gene silencing. We propose that the tDNAs primarily affect local chromatin structure, which results in effects on long-range chromosome architecture.

The Ribosome as a Missing Link in Prebiotic Evolution III: Over-Representation of tRNA- and rRNA-Like Sequences and Plieofunctionality of Ribosome-Related Molecules Argues for the Evolution of Primitive Genomes from Ribosomal RNA Modules

We propose that ribosomal RNA (rRNA) formed the basis of the first cellular genomes, and provide evidence from a review of relevant literature and proteonomic tests. We have proposed previously that the ribosome may represent the vestige of the first self-replicating entity in which rRNAs also functioned as genes that were transcribed into functional messenger RNAs (mRNAs) encoding ribosomal proteins. rRNAs also encoded polymerases to replicate itself and a full complement of the transfer RNAs (tRNAs) required to translate its genes. We explore here a further prediction of our “ribosome-first” theory: the ribosomal genome provided the basis for the first cellular genomes. Modern genomes should therefore contain an unexpectedly large percentage of tRNA- and rRNA-like modules derived from both sense and antisense reading frames, and these should encode non-ribosomal proteins, as well as ribosomal ones with key cell functions. Ribosomal proteins should also have been co-opted by cellular evolution to play extra-ribosomal functions. We review existing literature supporting these predictions. We provide additional, new data demonstrating that rRNA-like sequences occur at significantly higher frequencies than predicted on the basis of mRNA duplications or randomized RNA sequences. These data support our “ribosome-first” theory of cellular evolution.

SMORE: Synteny Modulator of Repetitive Elements

Several families of multicopy genes, such as transfer ribonucleic acids (tRNAs) and ribosomal RNAs (rRNAs), are subject to concerted evolution, an effect that keeps sequences of paralogous genes effectively identical. Under these circumstances, it is impossible to distinguish orthologs from paralogs on the basis of sequence similarity alone. Synteny, the preservation of relative genomic locations, however, also remains informative for the disambiguation of evolutionary relationships in this situation. In this contribution, we describe an automatic pipeline for the evolutionary analysis of such cases that use genome-wide alignments as a starting point to assign orthology relationships determined by synteny. The evolution of tRNAs in primates as well as the history of the Y RNA family in vertebrates and nematodes are used to showcase the method. The pipeline is freely available.

The origin and evolution of six miniature inverted-repeat transposable elements in Bombyx mori and Rhodnius prolixus

Miniature inverted-repeat transposable elements (MITEs) are a specific group of nonautonomous DNA transposons, and they are distributed in a wide range of hosts. However, the origin and evolutionary history of MITEs in eukaryotic genomes remain unclear. In this study, six MITEs were identified in the silkworm (Bombyx mori). Five elements are grouped into four known superfamilies of DNA transposons, and one represents a novel class of MITEs. Unexpectedly, six similar MITEs are also present in the triatomine bug (Rhodnius prolixus) that diverged from the common ancestor with the silkworm about 370 Ma. However, they show different lengths in two species, suggesting that they are different derivatives of progenitor transposons. Three direct progenitor transposons (Sola1, hobo/Ac/Tam [hAT], and Ginger2) are also identified in some other organisms, and several lines of evidence suggested that these autonomous elements might have been independently and horizontally transferred into their hosts. Furthermore, it is speculated that the twisted-wing parasites may be the candidate vectors for these horizontal transfers. The data presented in this study provide some new insights into the origin and evolutionary history of MITEs in the silkworm and triatomine bug.

Evolutionary patterns of non-coding RNAs

A plethora of new functions of non-coding RNAs (ncRNAs) have been discovered in past few years. In fact, RNA is emerging as the central player in cellular regulation, taking on active roles in multiple regulatory layers from transcription, RNA maturation, and RNA modification to translational regulation. Nevertheless, very little is known about the evolution of this "Modern RNA World" and its components. In this contribution, we attempt to provide at least a cursory overview of the diversity of ncRNAs and functional RNA motifs in non-translated regions of regular messenger RNAs (mRNAs) with an emphasis on evolutionary questions. This survey is complemented by an in-depth analysis of examples from different classes of RNAs focusing mostly on their evolution in the vertebrate lineage. We present a survey of Y RNA genes in vertebrates and study the molecular evolution of the U7 snRNA, the snoRNAs E1/U17, E2, and E3, the Y RNA family, the let-7 microRNA (miRNA) family, and the mRNA-like evf-1 gene. We furthermore discuss the statistical distribution of miRNAs in metazoans, which suggests an explosive increase in the miRNA repertoire in vertebrates. The analysis of the transcription of ncRNAs suggests that small RNAs in general are genetically mobile in the sense that their association with a hostgene (e.g. when transcribed from introns of a mRNA) can change on evolutionary time scales. The let-7 family demonstrates, that even the mode of transcription (as intron or as exon) can change among paralogous ncRNA.

tRNA-like elements in Haloferax volcanii

All functional RNAs are generated from precursor molecules by a plethora of processing steps. The generation of mature RNA molecules by processing is an important layer of gene expression regulation catalysed by ribonucleases. Here, we analysed 5S rRNA processing in the halophilic Archaeon Haloferax volcanii. Earlier experiments showed that the 5S rRNA is cleaved at its 5' end by the endonuclease tRNase Z. Interestingly, a tRNA-like structure was identified upstream of the 5S rRNA that might be used as a processing signal. Here, we show that this tRNA-like element is indeed recognised as a processing signal by tRNase Z. Substrates containing mutations in the tRNA-like sequence are no longer processed, whereas a substrate containing a deletion in the 5S rRNA sequence is still cleaved. Therefore, an intact 5S rRNA structure is not required for processing. Further, we used bioinformatics analyses to identify additional sequences in Haloferax containing tRNA-like structures. This search resulted in the identification of all tRNAs, the tRNA-like structure upstream of the 5S RNA and 47 new tRNA-like structural elements. However, the in vitro processing of selected examples showed no cleavage of these newly identified elements. Thus, tRNA-like elements are not a general processing signal in Haloferax.

Human tRNA genes function as chromatin insulators

Insulators help separate active chromatin domains from silenced ones. In yeast, gene promoters act as insulators to block the spread of Sir and HP1 mediated silencing while in metazoans most insulators are multipartite autonomous entities. tDNAs are repetitive sequences dispersed throughout the human genome and we now show that some of these tDNAs can function as insulators in human cells. Using computational methods, we identified putative human tDNA insulators. Using silencer blocking, transgene protection and repressor blocking assays we show that some of these tDNA-containing fragments can function as barrier insulators in human cells. We find that these elements also have the ability to block enhancers from activating RNA pol II transcribed promoters. Characterization of a putative tDNA insulator in human cells reveals that the site possesses chromatin signatures similar to those observed at other better-characterized eukaryotic insulators. Enhanced 4C analysis demonstrates that the tDNA insulator makes long-range chromatin contacts with other tDNAs and ETC sites but not with intervening or flanking RNA pol II transcribed genes.

A novel class of miniature inverted repeat transposable elements (MITEs) that contain hitchhiking (GTCY)(n) microsatellites

The movement of miniature inverted repeat transposable elements (MITEs) modifies genome structure and function. We describe the microsatellite-associated interspersed nuclear element 2 (MINE-2), that integrates at consensus WTTTT target sites, creates dinucleotide TT target site duplications (TSDs), and forms predicted MITE-like secondary structures; a 5' subterminal inverted repeat (SIR; AGGGTTCCGTAG) that is partially complementary to a 5' inverted repeat (IR; ACGAAGCCCT) and 3'-SIRs (TTACGGAACCCT). A (GTCY)(n) microsatellite is hitchhiking downstream of conserved 5'MINE-2 secondary structures, causing flanking sequence similarity amongst mobile microsatellite loci. Transfection of insect cell lines indicates that MITE-like secondary structures are sufficient to mediate genome integration, and provides insight into the transposition mechanism used by MINE-2s.

Genomic organization of eukaryotic tRNAs

BACKGROUND

Surprisingly little is known about the organization and distribution of tRNA genes and tRNA-related sequences on a genome-wide scale. While tRNA gene complements are usually reported in passing as part of genome annotation efforts, and peculiar features such as the tandem arrangements of tRNA gene in Entamoeba histolytica have been described in some detail, systematic comparative studies are rare and mostly restricted to bacteria. We therefore set out to survey the genomic arrangement of tRNA genes and pseudogenes in a wide range of eukaryotes to identify common patterns and taxon-specific peculiarities.

RESULTS

In line with previous reports, we find that tRNA complements evolve rapidly and tRNA gene and pseudogene locations are subject to rapid turnover. At phylum level, the distributions of the number of tRNA genes and pseudogenes numbers are very broad, with standard deviations on the order of the mean. Even among closely related species we observe dramatic changes in local organization. For instance, 65% and 87% of the tRNA genes and pseudogenes are located in genomic clusters in zebrafish and stickleback, resp., while such arrangements are relatively rare in the other three sequenced teleost fish genomes. Among basal metazoa, Trichoplax adherens has hardly any duplicated tRNA gene, while the sea anemone Nematostella vectensis boasts more than 17000 tRNA genes and pseudogenes. Dramatic variations are observed even within the eutherian mammals. Higher primates, for instance, have 616 +/- 120 tRNA genes and pseudogenes of which 17% to 36% are arranged in clusters, while the genome of the bushbaby Otolemur garnetti has 45225 tRNA genes and pseudogenes of which only 5.6% appear in clusters. In contrast, the distribution is surprisingly uniform across plant genomes. Consistent with this variability, syntenic conservation of tRNA genes and pseudogenes is also poor in general, with turn-over rates comparable to those of unconstrained sequence elements. Despite this large variation in abundance in Eukarya we observe a significant correlation between the number of tRNA genes, tRNA pseudogenes, and genome size.

CONCLUSIONS

The genomic organization of tRNA genes and pseudogenes shows complex lineage-specific patterns characterized by an extensive variability that is in striking contrast to the extreme levels of sequence-conservation of the tRNAs themselves. The comprehensive analysis of the genomic organization of tRNA genes and pseudogenes in Eukarya provides a basis for further studies into the interplay of tRNA gene arrangements and genome organization in general.

A sequence similar to tRNA 3 Lys gene is embedded in HIV-1 U3-R and promotes minus-strand transfer

We identified a sequence embedded in the U3/R region of HIV-1 RNA that is highly complementary to human tRNA₃^Lys. The free energy of annealing to tRNA₃^Lys is significantly lower for this sequence and the primer-binding site than for other similar length viral sequences. The only interruption in complementarity is a 29-nucleotide segment inserted where a tRNA intron would be expected. The insert contains the TATA box for viral RNA transcription. The embedded sequence includes a nine-nucleotide segment previously reported to aid minus strand transfer by binding the primer tRNA₃^Lys. Reconstituting transfer in vitro, we show that including segments from the embedded sequence in the acceptor template, beyond the nine nucleotides, further increases transfer efficiency. We propose that a tRNA₃^Lys gene was incorporated during HIV-1 evolution and retained largely intact because of its roles in transcription and strand transfer.

Using triplet periodicity of nucleotide sequences for finding potential reading frame shifts in genes

We introduce a novel approach for the detection of possible mutations leading to a reading frame (RF) shift in a gene. Deletions and insertions of DNA coding regions are considerable events for genes because an RF shift results in modifications of the extensive region of amino acid sequence coded by a gene. The suggested method is based on the phenomenon of triplet periodicity (TP) in coding regions of genes and its relative resistance to substitutions in DNA sequence. We attempted to extend 326 933 regions of continuous TP found in genes from the KEGG databank by considering possible insertions and deletions. We revealed totally 824 genes where such extension was possible and statistically significant. Then we generated amino acid sequences according to active (KEGG's) and hypothetically ancient RFs in order to find confirmation of a shift at a protein level. Consequently, 64 sequences have protein similarities only for ancient RF, 176 only for active RF, 3 for both and 581 have no protein similarity at all. We aimed to have revealed lower bound for the number of genes in which a shift between RF and TP is possible. Further ways to increase the number of revealed RF shifts are discussed.

Search and Classification of Potential Minisatellite Sequences from Bacterial Genomes

We used the method of Information Decomposition developed by us to identify the latent dinucleotide periodicity regions in bacterial genomes. The number of potential minisatellite sequences obtained at high level of statistical significance was 454. Then we classified the periodicity matrices and obtained 45 classes. We used the other new method developed by us--Modified Profile Analysis--to reveal more periodic sequences in the presence of indels using the classes obtained. The number of sequences found by combination of these two methods was 3949. Most of them cannot be revealed by other methods including dynamic programming and Fourier transformation.

Evaluation of tRNA gene PCR for identification of mollicutes

We evaluated the applicability of tRNA gene PCR in combination with fluorescent capillary electrophoresis with an ABI310 genetic analyzer (Applied Biosystems, Calif.) for the identification of different mollicute species. A total of 103 strains and DNA extracts of 30 different species belonging to the genera Acholeplasma, Mycoplasma, and Ureaplasma were studied. Reproducible peak profiles were generated for all samples, except for one M. genitalium isolate, the three M. gallisepticum isolates, and 8 of the 24 Ureaplasma cultures, where no amplification could be obtained. Clustering revealed numerous discrepancies compared to the identifications that had been previously obtained by means of biochemical and serological tests. Final identification was obtained by 16S rRNA gene amplification followed by sequence analysis and/or restriction digestion. This confirmed the identification obtained by tRNA gene PCR in all cases. Seven samples yielded an unexpected tRNA gene PCR profile. Sequence analysis of the 16S rRNA genes showed that six of these samples were mixed and that one had a unique sequence that did not match any of the published sequences, pointing to the existence of a not-yet-described species. In conclusion, we found tRNA gene PCR to be a rapid and discriminatory method to correctly identify a large collection of different species of the class of Mollicutes and to recognize not-yet-described groups.