Misfolded human tRNA isodecoder binds and neutralizes a 3' UTR-embedded Alu element

tRNA derived small RNAs—Small players with big roles

In the past 2 decades, small non-coding RNAs derived from tRNA (tsRNAs or tRNA derived fragments; tRFs) have emerged as new powerful players in the field of small RNA mediated regulation of gene expression, translation, and epigenetic control. tRFs have been identified from evolutionarily divergent organisms from Archaea, the higher plants, to humans. Recent studies have confirmed their roles in cancers and other metabolic disorders in humans and experimental models. They have been implicated in biotic and abiotic stress responses in plants as well. In this review, we summarize the current knowledge on tRFs including types of tRFs, their biogenesis, and mechanisms of action. The review also highlights recent studies involving differential expression profiling of tRFs and elucidation of specific functions of individual tRFs from various species. We also discuss potential considerations while designing experiments involving tRFs identification and characterization and list the available bioinformatics tools for this purpose.

Machine Learning Algorithms Highlight tRNA Information Content and Chargaff’s Second Parity Rule Score as Important Features in Discriminating Probiotics from Non-Probiotics

Simple Summary

Probiotics are a group of beneficial microorganisms that are symbionts of the human gut microbiome. The identification of new probiotics is therefore of paramount importance from both public health and commercial perspectives. In this study, we show for the first time that Artificial Intelligence algorithms can identify novel probiotics and also discriminate them from pathogenic organisms in the human gut. We were also able to determine the information content within tRNA sequences as the key genomic features capable of characterizing probiotics.

Abstract

Probiotic bacteria are microorganisms with beneficial effects on human health and are currently used in numerous food supplements. However, no selection process is able to effectively distinguish probiotics from non-probiotic organisms on the basis of their genomic characteristics. In the current study, four Machine Learning algorithms were employed to accurately identify probiotic bacteria based on their DNA characteristics. Although the prediction accuracies of all algorithms were excellent, the Neural Network returned the highest scores in all the evaluation metrics, managing to discriminate probiotics from non-probiotics with an accuracy greater than 90%. Interestingly, our analysis also highlighted the information content of the tRNA sequences as the most important feature in distinguishing the two groups of organisms probably because tRNAs have regulatory functions and might have allowed probiotics to evolve faster in the human gut environment. Through the methodology presented here, it was also possible to identify seven promising new probiotics that have a higher information content in their tRNA sequences compared to non-probiotics. In conclusion, we prove for the first time that Machine Learning methods can discriminate human probiotic from non-probiotic organisms underlining information within tRNA sequences as the most important genomic feature in distinguishing them.

Genome-wide interrogation of transfer RNA-derived small RNAs in a mouse model of traumatic brain injury

Transfer RNA (tRNA)-derived small RNAs (tsRNAs) are a recently established family of regulatory small non-coding RNAs that modulate diverse biological processes. Growing evidence indicates that tsRNAs are involved in neurological disorders and play a role in the pathogenesis of neurodegenerative disease. However, whether tsRNAs are involved in traumatic brain injury-induced secondary injury remains poorly understood. In this study, a mouse controlled cortical impact model of traumatic brain injury was established, and integrated tsRNA and messenger RNA (mRNA) transcriptome sequencing were used. The results revealed that 103 tsRNAs were differentially expressed in the mouse model of traumatic brain injury at 72 hours, of which 56 tsRNAs were upregulated and 47 tsRNAs were downregulated. Based on microRNA-like seed matching and Pearson correlation analysis, 57 differentially expressed tsRNA-mRNA interaction pairs were identified, including 29 tsRNAs and 26 mRNAs. Moreover, Gene Ontology annotation of target genes revealed that the significantly enriched terms were primarily associated with inflammation and synaptic function. Collectively, our findings suggest that tsRNAs may be associated with traumatic brain injury-induced secondary brain injury, and are thus a potential therapeutic target for traumatic brain injury. The study was approved by the Beijing Neurosurgical Institute Animal Care and Use Committee (approval No. 20190411) on April 11, 2019.

On the Track of the Missing tRNA Genes: A Source of Non-Canonical Functions?

Cellular tRNAs appear today as a diverse population of informative macromolecules with conserved general elements ensuring essential common functions and different and distinctive features securing specific interactions and activities. Their differential expression and the variety of post-transcriptional modifications they are subject to, lead to the existence of complex repertoires of tRNA populations adjusted to defined cellular states. Despite the tRNA-coding genes redundancy in prokaryote and eukaryote genomes, it is surprising to note the absence of genes coding specific translational-active isoacceptors throughout the phylogeny. Through the analysis of different releases of tRNA databases, this review aims to provide a general summary about those “missing tRNA genes.” This absence refers to both tRNAs that are not encoded in the genome, as well as others that show critical sequence variations that would prevent their activity as canonical translation adaptor molecules. Notably, while a group of genes are universally missing, others are absent in particular kingdoms. Functional information available allows to hypothesize that the exclusion of isodecoding molecules would be linked to: 1) reduce ambiguities of signals that define the specificity of the interactions in which the tRNAs are involved; 2) ensure the adaptation of the translational apparatus to the cellular state; 3) divert particular tRNA variants from ribosomal protein synthesis to other cellular functions. This leads to consider the “missing tRNA genes” as a source of putative non-canonical tRNA functions and to broaden the concept of adapter molecules in ribosomal-dependent protein synthesis.

Hijacking tRNAs From Translation: Regulatory Functions of tRNAs in Mammalian Cell Physiology

Transfer tRNAs (tRNAs) are small non-coding RNAs that are highly conserved in all kingdoms of life. Originally discovered as the molecules that deliver amino acids to the growing polypeptide chain during protein synthesis, tRNAs have been believed for a long time to play exclusive role in translation. However, recent studies have identified key roles for tRNAs and tRNA-derived small RNAs in multiple other processes, including regulation of transcription and translation, posttranslational modifications, stress response, and disease. These emerging roles suggest that tRNAs may be central players in the complex machinery of biological regulatory pathways. Here we overview these non-canonical roles of tRNA in normal physiology and disease, focusing largely on eukaryotic and mammalian systems.

Noncanonical Roles of tRNAs: tRNA Fragments and Beyond

As one of the most abundant and conserved RNA species, transfer RNAs (tRNAs) are well known for their role in reading the codons on messenger RNAs and translating them into proteins. In this review, we discuss the noncanonical functions of tRNAs. These include tRNAs as precursors to novel small RNA molecules derived from tRNAs, also called tRNA-derived fragments, that are abundant across species and have diverse functions in different biological processes, including regulating protein translation, Argonaute-dependent gene silencing, and more. Furthermore, the role of tRNAs in biosynthesis and other regulatory pathways, including nutrient sensing, splicing, transcription, retroelement regulation, immune response, and apoptosis, is reviewed. Genome organization and sequence variation of tRNA genes are also discussed in light of their noncanonical functions. Lastly, we discuss the recent applications of tRNAs in genome editing and microbiome sequencing.

Dimerization confers increased stability to nucleases in 5' halves from glycine and glutamic acid tRNAs

We have previously shown that 5′ halves from tRNA^Gly_GCC and tRNA^Glu_CUC are the most enriched small RNAs in the extracellular space of human cell lines, and especially in the non-vesicular fraction. Extracellular RNAs are believed to require protection by either encapsulation in vesicles or ribonucleoprotein complex formation. However, deproteinization of non-vesicular tRNA halves does not affect their retention in size-exclusion chromatography. Thus, we considered alternative explanations for their extracellular stability. In-silico analysis of the sequence of these tRNA-derived fragments showed that tRNA^Gly 5′ halves can form homodimers or heterodimers with tRNA^Glu 5′ halves. This capacity is virtually unique to glycine tRNAs. By analyzing synthetic oligonucleotides by size exclusion chromatography, we provide evidence that dimerization is possible in vitro. tRNA halves with single point substitutions preventing dimerization are degraded faster both in controlled nuclease digestion assays and after transfection in cells, showing that dimerization can stabilize tRNA halves against the action of cellular nucleases. Finally, we give evidence supporting dimerization of endogenous tRNA^Gly_GCC 5′ halves inside cells. Considering recent reports have shown that 5′ tRNA halves from Ala and Cys can form tetramers, our results highlight RNA intermolecular structures as a new layer of complexity in the biology of tRNA-derived fragments.

History of tRNA research in strasbourg

The tRNA molecules, in addition to translating the genetic code into protein and defining the second genetic code via their aminoacylation by aminoacyl-tRNA synthetases, act in many other cellular functions and dysfunctions. This article, illustrated by personal souvenirs, covers the history of ~60 years tRNA research in Strasbourg. Typical examples point up how the work in Strasbourg was a two-way street, influenced by and at the same time influencing investigators outside of France. All along, research in Strasbourg has nurtured the structural and functional diversity of tRNA. It produced massive sequence and crystallographic data on tRNA and its partners, thereby leading to a deeper physicochemical understanding of tRNA architecture, dynamics, and identity. Moreover, it emphasized the role of nucleoside modifications and in the last two decades, highlighted tRNA idiosyncrasies in plants and organelles, together with cellular and health-focused aspects. The tRNA field benefited from a rich local academic heritage and a strong support by both university and CNRS. Its broad interlinks to the worldwide community of tRNA researchers opens to an exciting future. © 2019 IUBMB Life, 2019 © 2019 IUBMB Life, 71(8):1066-1087, 2019.

Differential expression of human tRNA genes drives the abundance of tRNA-derived fragments

The human genome encodes hundreds of transfer RNA (tRNA) genes but their individual contribution to the tRNA pool is not fully understood. Deep sequencing of tRNA transcripts (tRNA-Seq) can estimate tRNA abundance at single gene resolution, but tRNA structures and posttranscriptional modifications impair these analyses. Here we present a bioinformatics strategy to investigate differential tRNA gene expression and use it to compare tRNA-Seq datasets from cultured human cells and human brain. We find that sequencing caveats affect quantitation of only a subset of human tRNA genes. Unexpectedly, we detect several cases where the differences in tRNA expression among samples do not involve variations at the level of isoacceptor tRNA sets (tRNAs charged with the same amino acid but using different anticodons), but rather among tRNA genes within the same isodecoder set (tRNAs having the same anticodon sequence). Because isodecoder tRNAs are functionally equal in terms of genetic translation, their differential expression may be related to noncanonical tRNA functions. We show that several instances of differential tRNA gene expression result in changes in the abundance of tRNA-derived fragments (tRFs) but not of mature tRNAs. Examples of differentially expressed tRFs include PIWI-associated RNAs, tRFs present in tissue samples but not in cells cultured in vitro, and somatic tissue-specific tRFs. Our data support that differential expression of tRNA genes regulate noncanonical tRNA functions performed by tRFs.

A statistical analysis on transcriptome sequences: The enrichment of Alu-element is associated with subcellular location

The Alu-element plays important roles in mediating alternative splicing, RNA editing and translation regulation. However, the distribution and function of the Alu-element are never analysed at the transcriptome level. This study presents a statistical analysis of the Alu-element on human transcriptome. We found that mRNAs and lncRNAs share the same sequence form for the Alu-element. The Alu-element covers 5.8% of the coding transcripts and 17.1% of the coding genes for mRNAs, and covers 9.3% of the transcripts and 13.6% of the genes for lncRNAs. The Alu-element is preferentially located at the 3' end. Statistical analysis demonstrates that the enrichment of Alu-element is associated with subcellular location. For instance, Alu-inclusive transcripts are overexpressed in nucleus, mitochondrion and Golgi apparatus membrane while under-expressed in cell membrane and extracellular space. We found that genes contain both Alu-element and S- domains of 7SL RNA are all associated with cellular activities carried out in mitochondrion.

Data Mining of Small RNA-Seq Suggests an Association Between Prostate Cancer and Altered Abundance of 5' Transfer RNA Halves in Seminal Fluid and Prostatic Tissues

Extracellular RNAs are gaining clinical interest as biofluid-based noninvasive markers for diseases, especially cancer. In particular, derivatives of transfer RNA (tRNA) are emerging as a new class of small-noncoding RNAs with high biomarker potential. We and others previously reported alterations in serum levels of specific tRNA halves in disease states including cancer. Here, we explored seminal fluid for tRNA halves as potential markers of prostate cancer. We found that 5′ tRNA halves are abundant in seminal fluid and are elevated in prostate cancer relative to noncancer patients. Importantly, most of these tRNA halves are also detectable in prostatic tissues, and a subset were increased in malignant relative to adjacent normal tissue. These findings emphasize the potential of 5′ tRNA halves as noninvasive markers for prostate cancer screening and diagnosis and provide leads for future work to elucidate a putative role of the 5′ tRNA halves in carcinogenesis.

Modifications and functional genomics of human transfer RNA

Transfer RNA (tRNA) is present at tens of millions of transcripts in a human cell and is the most abundant RNA in moles among all cellular RNAs. tRNA is also the most extensively modified RNA with, on an average, 13 modifications per molecule. The primary function of tRNA as the adaptor of amino acids and the genetic code in protein synthesis is well known. tRNA modifications play multi-faceted roles in decoding and other cellular processes. The abundance, modification, and aminoacylation (charging) levels of tRNAs contribute to mRNA decoding in ways that reflect the cell type and its environment; however, how these factors work together to maximize translation efficiency remains to be understood. tRNAs also interact with many proteins not involved in translation and this may coordinate translation activity and other processes in the cell. This review focuses on the modifications and the functional genomics of human tRNA and discusses future perspectives on the explorations of human tRNA biology.

Transfer-RNA-mediated enhancement of ribosomal proteins S6 kinases signaling for cell proliferation

While transfer-RNAs (tRNAs) are known to transport amino acids to ribosome, new functions are being unveiled from tRNAs and their fragments beyond protein synthesis. Here we show that phosphorylation of 90-kDa RPS6K (ribosomal proteins S6 kinase) was enhanced by tRNA^Leu overexpression under amino acids starvation condition. The phosphorylation of 90-kDa RPS6K was decreased by siRNA specific to tRNA^Leu and was independent to mTOR (mammalian target of rapamycin) signaling. Among the 90-kDa RPS6K family, RSK1 (ribosomal S6 kinase 1) and MSK2 (mitogen-and stress-activated protein kinase 2) were the major kinases phosphorylated by tRNA^Leu overexpression. Through SILAC (stable isotope labeling by/with amino acids in cell culture) and combined mass spectrometry analysis, we identified EBP1 (ErbB3-binding protein 1) as the tRNA^Leu-binding protein. We suspected that the overexpression of free tRNA^Leu would reinforce ErbB2/ErbB3 signaling pathway by disturbing the interaction between ErbB3 and EBP1, resulting in RSK1/MSK2 phosphorylation, improving cell proliferation and resistance to death. Analysis of samples from patients with breast cancer also indicated an association between tRNA^Leu overexpression and the ErbB2-positive population. Our results suggested a possible link between tRNA^Leu overexpression and RSK1/MSK2 activation and ErbB2/ErbB3 signaling.

A comprehensive repertoire of tRNA-derived fragments in prostate cancer

Prostate cancer (PCa) is the most common cancer among men in developed countries. Although its genetic background is thoroughly investigated, rather little is known about the role of small non-coding RNAs (sncRNA) in this disease. tRNA-derived fragments (tRFs) represent a new class of sncRNAs, which are present in a broad range of species and have been reported to play a role in several cellular processes. Here, we analyzed the expression of tRFs in fresh frozen patient samples derived from normal adjacent prostate and different stages of PCa by RNA-sequencing. We identified 598 unique tRFs, many of which are deregulated in cancer samples when compared to normal adjacent tissue. Most of the identified tRFs are derived from the 5’- and 3’-ends of mature cytosolic tRNAs, but we also found tRFs produced from other parts of tRNAs, including pre-tRNA trailers and leaders, as well as tRFs from mitochondrial tRNAs. The 5’-derived tRFs comprise the most abundant class of tRFs in general and represent the major class among upregulated tRFs. The 3’-derived tRFs types are dominant among downregulated tRFs in PCa. We validated the expression of three tRFs using qPCR. The ratio of tRFs derived from tRNA^LysCTT and tRNA^PheGAA emerged as a good indicator of progression-free survival and a candidate prognostic marker. This study provides a systematic catalogue of tRFs and their dysregulation in PCa and can serve as the basis for further research on the biomarker potential and functional roles of tRFs in this disease.

Swellix: a computational tool to explore RNA conformational space

Background

The sequence of nucleotides in an RNA determines the possible base pairs for an RNA fold and thus also determines the overall shape and function of an RNA. The Swellix program presented here combines a helix abstraction with a combinatorial approach to the RNA folding problem in order to compute all possible non-pseudoknotted RNA structures for RNA sequences. The Swellix program builds on the Crumple program and can include experimental constraints on global RNA structures such as the minimum number and lengths of helices from crystallography, cryoelectron microscopy, or in vivo crosslinking and chemical probing methods.

Results

The conceptual advance in Swellix is to count helices and generate all possible combinations of helices rather than counting and combining base pairs. Swellix bundles similar helices and includes improvements in memory use and efficient parallelization. Biological applications of Swellix are demonstrated by computing the reduction in conformational space and entropy due to naturally modified nucleotides in tRNA sequences and by motif searches in Human Endogenous Retroviral (HERV) RNA sequences. The Swellix motif search reveals occurrences of protein and drug binding motifs in the HERV RNA ensemble that do not occur in minimum free energy or centroid predicted structures.

Conclusions

Swellix presents significant improvements over Crumple in terms of efficiency and memory use. The efficient parallelization of Swellix enables the computation of sequences as long as 418 nucleotides with sufficient experimental constraints. Thus, Swellix provides a practical alternative to free energy minimization tools when multiple structures, kinetically determined structures, or complex RNA-RNA and RNA-protein interactions are present in an RNA folding problem.

Electronic supplementary material

The online version of this article (10.1186/s12859-017-1910-7) contains supplementary material, which is available to authorized users.

Determination of tRNA aminoacylation levels by high-throughput sequencing

Transfer RNA (tRNA) decodes mRNA codons when aminoacylated (charged) with an amino acid at its 3' end. Charged tRNAs turn over rapidly in cells, and variations in charged tRNA fractions are known to be a useful parameter in cellular responses to stress. tRNA charging fractions can be measured for individual tRNA species using acid denaturing gels, or comparatively at the genome level using microarrays. These hybridization-based approaches cannot be used for high resolution analysis of mammalian tRNAs due to their large sequence diversity. Here we develop a high-throughput sequencing method that enables accurate determination of charged tRNA fractions at single-base resolution (Charged DM-tRNA-seq). Our method takes advantage of the recently developed DM-tRNA-seq method, but includes additional chemical steps that specifically remove the 3'A residue in uncharged tRNA. Charging fraction is obtained by counting the fraction of A-ending reads versus A+C-ending reads for each tRNA species in the same sequencing reaction. In HEK293T cells, most cytosolic tRNAs are charged at >80% levels, whereas tRNASer and tRNAThr are charged at lower levels. These low charging levels were validated using acid denaturing gels. Our method should be widely applicable for investigations of tRNA charging as a parameter in biological regulation.

GtRNAdb 2.0: an expanded database of transfer RNA genes identified in complete and draft genomes

Transfer RNAs represent the largest, most ubiquitous class of non-protein coding RNA genes found in all living organisms. The tRNAscan-SE search tool has become the de facto standard for annotating tRNA genes in genomes, and the Genomic tRNA Database (GtRNAdb) was created as a portal for interactive exploration of these gene predictions. Since its published description in 2009, the GtRNAdb has steadily grown in content, and remains the most commonly cited web-based source of tRNA gene information. In this update, we describe not only a major increase in the number of tRNA predictions (>367000) and genomes analyzed (>4370), but more importantly, the integration of new analytic and functional data to improve the quality and biological context of tRNA gene predictions. New information drawn from other sources includes tRNA modification data, epigenetic data, single nucleotide polymorphisms, gene expression and evolutionary conservation. A richer set of analytic data is also presented, including better tRNA functional prediction, non-canonical features, predicted structural impacts from sequence variants and minimum free energy structural predictions. Views of tRNA genes in genomic context are provided via direct links to the UCSC genome browsers. The database can be searched by sequence or gene features, and is available at http://gtrnadb.ucsc.edu/.

Regulatory roles of Alu transcript on gene expression

Alu element is the most successful transposon and it maintains a high level of content in primate genome. However, despite the fact that the expression level of independent Alu element is rather low under common condition, an increasing number of the observations for the Alu transcripts in cells and tissues have been reported recently. Alu transcripts play key roles in the network of gene expression regulation. The main functions of Alu transcript focus on gene regulation both at transcriptional and post-transcriptional levels. This review summarizes major functions of Alu transcripts on gene expression and highlights molecular mechanisms dependent on conserved sequence or secondary structure in order to unravel a relative ubiquitous way that Alu transcript uses to affect the whole genome.

The RNA-centred view of the synapse: non-coding RNAs and synaptic plasticity

If mRNAs were the only RNAs made by a neuron, there would be a simple mapping of mRNAs to proteins. However, microRNAs and other non-coding RNAs (ncRNAs; endo-siRNAs, piRNAs, BC1, BC200, antisense and long ncRNAs, repeat-related transcripts, etc.) regulate mRNAs via effects on protein translation as well as transcriptional and epigenetic mechanisms. Not only are genes ON or OFF, but their ability to be translated can be turned ON or OFF at the level of synapses, supporting an enormous increase in information capacity. Here, I review evidence that ncRNAs are expressed pervasively within dendrites in mammalian brain; that some are activity-dependent and highly enriched near synapses; and that synaptic ncRNAs participate in plasticity responses including learning and memory. Ultimately, ncRNAs can be viewed as the post-it notes of the neuron. They have no literal meaning of their own, but derive their functions from where (and to what) they are stuck. This may explain, in part, why ncRNAs differ so dramatically from protein-coding genes, both in terms of the usual indicators of functionality and in terms of evolutionary constraints. ncRNAs do not appear to be direct mediators of synaptic transmission in the manner of neurotransmitters or receptors, yet they orchestrate synaptic plasticity—and may drive species-specific changes in cognition.

A parallel implementation of the Wuchty algorithm with additional experimental filters to more thoroughly explore RNA conformational space

We present new modifications to the Wuchty algorithm in order to better define and explore possible conformations for an RNA sequence. The new features, including parallelization, energy-independent lonely pair constraints, context-dependent chemical probing constraints, helix filters, and optional multibranch loops, provide useful tools for exploring the landscape of RNA folding. Chemical probing alone may not necessarily define a single unique structure. The helix filters and optional multibranch loops are global constraints on RNA structure that are an especially useful tool for generating models of encapsidated viral RNA for which cryoelectron microscopy or crystallography data may be available. The computations generate a combinatorially complete set of structures near a free energy minimum and thus provide data on the density and diversity of structures near the bottom of a folding funnel for an RNA sequence. The conformational landscapes for some RNA sequences may resemble a low, wide basin rather than a steep funnel that converges to a single structure.

Probing viral genomic structure: alternative viewpoints and alternative structures for satellite tobacco mosaic virus RNA

Viral RNA structure prediction is a valuable tool for development of drugs against viral disease. This work discusses different approaches to predicting encapsidated viral RNA and highlights satellite tobacco mosaic virus (STMV) RNA as a model system with excellent crystallography data. Fundamentally important issues for debate include thermodynamic versus kinetic control of virus assembly and the possible consequences of quasi-species in the primary structure on RNA secondary structure prediction of a single structure or an ensemble of structures. Multiple computational tools and chemical reagents are now available for improved viral RNA structure prediction. Two different predicted structures for encapsidated STMV RNA result from differences in three main areas: a different approach and philosophy to studying encapsidated viral RNA, an emphasis on different RNA motifs, and technical differences in computational methods and chemical reagents. The experiments with traditional chemical probing and SHAPE reagents are compared in terms of chemistry, results, and interpretation for STMV RNA as well as other RNA protein assemblies, such as the 5'UTR of HIV and the ribosome. This discussion of the challenges of viral RNA structure prediction will lead to new experiments and improved future predictions for viral RNA.

Regulation of translation dynamic and neoplastic conversion by tRNA and their pieces

Research on transfer RNA (tRNA) has gone a long way since the existence of this essential adapter of the genetic code was first hypothesized five decades ago. With the new and fascinating discovering of connections between tRNAs and cellular pathways beyond genetic translation, the field of tRNA research has reached a new era. Here, we review some aspects of the emerging variety of tasks performed by full length tRNAs as well as their fragments generated by specific nuclease cleavage. Topics of special focus include the effect of differential expression of tRNAs in healthy tissues as well as their frequent deregulation observed in cancer cells. We also discuss the central role played by tRNAMet in cell metabolism, proliferation, and response to oxidative stress. Finally we review evidences suggesting that tRNAs are critical sources of short RNAs regulating an ever growing variety of cellular processes including translation initiation, control of genomic retroviral sequences, or RNA interference.

Diversity of human tRNA genes from the 1000-genomes project

The sequence diversity of individual human genomes has been extensively analyzed for variations and phenotypic implications for mRNA, miRNA, and long non-coding RNA genes. TRNA (tRNA) also exhibits large sequence diversity in the human genome, but tRNA gene sequence variation and potential functional implications in individual human genomes have not been investigated. Here we capitalize on the sequencing data from the 1000-genomes project to examine the diversity of tRNA genes in the human population. Previous analysis of the reference human genome indicated an unexpected large number of diverse tRNA genes beyond the necessity of translation, suggesting that some tRNA transcripts may perform non-canonical functions. We found 24 new tRNA sequences in>1% and 76 new tRNA sequences in>0.2% of all individuals, indicating that tRNA genes are also subject to evolutionary changes in the human population. Unexpectedly, two abundant new tRNA genes contain base-pair mismatches in the anticodon stem. We experimentally determined that these two new tRNAs have altered structures in vitro; however, one new tRNA is not aminoacylated but extremely stable in HeLa cells, suggesting that this new tRNA can be used for non-canonical function. Our results show that at the scale of human population, tRNA genes are more diverse than conventionally understood, and some new tRNAs may perform non-canonical, extra-translational functions that may be linked to human health and disease.

Transcriptome-wide expansion of non-coding regulatory switches: evidence from co-occurrence of Alu exonization, antisense and editing

Non-coding RNAs from transposable elements of human genome are gaining prominence in modulating transcriptome dynamics. Alu elements, as exonized, edited and antisense components within same transcripts could create novel regulatory switches in response to different transcriptional cues. We provide the first evidence for co-occurrences of these events at transcriptome-wide scale through integrative analysis of data sets across diverse experimental platforms and tissues. This involved the following: (i) positional anchoring of Alu exonization events in the UTRs and CDS of 4663 transcript isoforms from RefSeq mRNAs and (ii) mapping on to them A→I editing events inferred from ∼7 million ESTs from dbEST and antisense transcripts identified from virtual serial analysis of gene expression tags represented in Cancer Genome Anatomy Project next-generation sequencing data sets across 20 tissues. We observed significant enrichment of these events in the 3'UTR as well as positional preference within the embedded Alus. More than 300 genes had co-occurrence of all these events at the exon level and were significantly enriched in apoptosis and lysosomal processes. Further, we demonstrate functional evidence of such dynamic interactions between Alu-mediated events in a time series data from Integrated Personal Omics Profiling during recovery from a viral infection. Such 'single transcript-multiple fate' opportunity facilitated by Alu elements may modulate transcriptional response, especially during stress.

Beyond microRNA--novel RNAs derived from small non-coding RNA and their implication in cancer

Over the recent years, Next Generation Sequencing (NGS) technologies targeting the microRNA transcriptome revealed the existence of many different RNA fragments derived from small RNA species other than microRNA. Although initially discarded as RNA turnover artifacts, accumulating evidence suggests that RNA fragments derived from small nucleolar RNA (snoRNA) and transfer RNA (tRNA) are not just random degradation products but rather stable entities, which may have functional activity in the normal and malignant cell. This review summarizes new findings describing the detection and alterations in expression of snoRNA-derived (sdRNA) and tRNA-derived (tRF) RNAs. We focus on the possible interactions of sdRNAs and tRFs with the canonical microRNA pathways in the cell and present current hypotheses on the function of these RNAs.

PlantRNA, a database for tRNAs of photosynthetic eukaryotes

PlantRNA database (http://plantrna.ibmp.cnrs.fr/) compiles transfer RNA (tRNA) gene sequences retrieved from fully annotated plant nuclear, plastidial and mitochondrial genomes. The set of annotated tRNA gene sequences has been manually curated for maximum quality and confidence. The novelty of this database resides in the inclusion of biological information relevant to the function of all the tRNAs entered in the library. This includes 5'- and 3'-flanking sequences, A and B box sequences, region of transcription initiation and poly(T) transcription termination stretches, tRNA intron sequences, aminoacyl-tRNA synthetases and enzymes responsible for tRNA maturation and modification. Finally, data on mitochondrial import of nuclear-encoded tRNAs as well as the bibliome for the respective tRNAs and tRNA-binding proteins are also included. The current annotation concerns complete genomes from 11 organisms: five flowering plants (Arabidopsis thaliana, Oryza sativa, Populus trichocarpa, Medicago truncatula and Brachypodium distachyon), a moss (Physcomitrella patens), two green algae (Chlamydomonas reinhardtii and Ostreococcus tauri), one glaucophyte (Cyanophora paradoxa), one brown alga (Ectocarpus siliculosus) and a pennate diatom (Phaeodactylum tricornutum). The database will be regularly updated and implemented with new plant genome annotations so as to provide extensive information on tRNA biology to the research community.

Transposable elements reveal a stem cell-specific class of long noncoding RNAs

Background

Numerous studies over the past decade have elucidated a large set of long intergenic noncoding RNAs (lincRNAs) in the human genome. Research since has shown that lincRNAs constitute an important layer of genome regulation across a wide spectrum of species. However, the factors governing their evolution and origins remain relatively unexplored. One possible factor driving lincRNA evolution and biological function is transposable element (TE) insertions. Here, we comprehensively characterize the TE content of lincRNAs relative to genomic averages and protein coding transcripts.

Results

Our analysis of the TE composition of 9,241 human lincRNAs revealed that, in sharp contrast to protein coding genes, 83% of lincRNAs contain a TE, and TEs comprise 42% of lincRNA sequence. lincRNA TE composition varies significantly from genomic averages - L1 and Alu elements are depleted and broad classes of endogenous retroviruses are enriched. TEs occur in biased positions and orientations within lincRNAs, particularly at their transcription start sites, suggesting a role in lincRNA transcriptional regulation. Accordingly, we observed a dramatic example of HERVH transcriptional regulatory signals correlating strongly with stem cell-specific expression of lincRNAs. Conversely, lincRNAs devoid of TEs are expressed at greater levels than lincRNAs with TEs in all tissues and cell lines, particularly in the testis.

Conclusions

TEs pervade lincRNAs, dividing them into classes, and may have shaped lincRNA evolution and function by conferring tissue-specific expression from extant transcriptional regulatory signals.

Genome-wide identification of novel small RNAs in Pseudomonas aeruginosa

Bacterial small regulatory RNAs (sRNAs) function in post-transcriptional control of gene expression and control a variety of processes including metabolic reactions, stress responses and pathogenesis in response to environmental signals. A variety of approaches have been used previously to identify 44 sRNAs in the opportunistic human pathogen Pseudomonas aeruginosa. In this work, RNA sequencing (RNA-seq) is used to identify novel transcripts in P.aeruginosa involving a combination of three different sequencing libraries. Almost all known sRNAs and over 500 novel intergenic sRNAs are identified with this approach. Although the use of three libraries increased the number of novel transcripts identified, there were significant differences in the subset of transcripts detected in each library, underscoring the importance of library preparation strategy and relative sRNA abundance for successful sRNA detection. Nearly 90% of the novel sRNAs have no orthologous bacterial sequences outside of P.aeruginosa, supporting a limited degree of sequence conservation and rapid evolution of sRNAs at the species level. We anticipate that the data will be useful for the study of regulatory sRNAs in bacteria and that the approach described here may be applied to identify sRNAs in any bacterium under different growth and stress conditions.