Nuclear and mitochondrial tRNA-lookalikes in the human genome

Ribonuclease Inhibitor 1 (RNH1) Regulates Sperm tsRNA Generation for Paternal Inheritance through Interacting with Angiogenin in the Caput Epididymis

Environmental stressors can induce paternal epigenetic modifications that are a key determinant of the intergenerational inheritance of acquired phenotypes in mammals. Some of them can affect phenotypic expression through inducing changes in tRNA-derived small RNAs (tsRNAs), which modify paternal epigenetic regulation in sperm. However, it is unclear how these stressors can affect changes in the expression levels of tsRNAs and their related endonucleases in the male reproductive organs. We found that Ribonuclease inhibitor 1 (RNH1), an oxidation responder, interacts with ANG to regulate sperm tsRNA generation in the mouse caput epididymis. On the other hand, inflammation and oxidative stress induced by either lipopolysaccharide (LPS) or palmitate (PA) treatments weakened the RNH1-ANG interaction in the epididymal epithelial cells (EEC). Accordingly, ANG translocation increased from the nucleus to the cytoplasm, which led to ANG upregulation and increases in cytoplasmic tsRNA expression levels. In conclusion, as an antioxidant, RNH1 regulates tsRNA generation through targeting ANG in the mouse caput epididymis. Moreover, the tsRNA is an epigenetic factor in sperm that modulates paternal inheritance in offspring via the fertilization process.

The Typical tRNA Co-Expresses Multiple 5' tRNA Halves Whose Sequences and Abundances Depend on Isodecoder and Isoacceptor and Change with Tissue Type, Cell Type, and Disease

Transfer RNA-derived fragments (tRFs) are noncoding RNAs that arise from either mature transfer RNAs (tRNAs) or their precursors. One important category of tRFs comprises the tRNA halves, which are generated through cleavage at the anticodon. A given tRNA typically gives rise to several co-expressed 5'-tRNA halves (5'-tRHs) that differ in the location of their 3' ends. These 5'-tRHs, even though distinct, have traditionally been treated as indistinguishable from one another due to their near-identical sequences and lengths. We focused on co-expressed 5'-tRHs that arise from the same tRNA and systematically examined their exact sequences and abundances across 10 different human tissues. To this end, we manually curated and analyzed several hundred human RNA-seq datasets from NCBI's Sequence Run Archive (SRA). We grouped datasets from the same tissue into their own collection and examined each group separately. We found that a given tRNA produces different groups of co-expressed 5'-tRHs in different tissues, different cell lines, and different diseases. Importantly, the co-expressed 5'-tRHs differ in their sequences, absolute abundances, and relative abundances, even among tRNAs with near-identical sequences from the same isodecoder or isoacceptor group. The findings suggest that co-expressed 5'-tRHs that are produced from the same tRNA or closely related tRNAs have distinct, context-dependent roles. Moreover, our analyses show that cell lines modeling the same tissue type and disease may not be interchangeable when it comes to experimenting with tRFs.

The diverse structural modes of tRNA binding and recognition

tRNAs are short noncoding RNAs responsible for decoding mRNA codon triplets, delivering correct amino acids to the ribosome, and mediating polypeptide chain formation. Due to their key roles during translation, tRNAs have a highly conserved shape and large sets of tRNAs are present in all living organisms. Regardless of sequence variability, all tRNAs fold into a relatively rigid three-dimensional L-shaped structure. The conserved tertiary organization of canonical tRNA arises through the formation of two orthogonal helices, consisting of the acceptor and anticodon domains. Both elements fold independently to stabilize the overall structure of tRNAs through intramolecular interactions between the D- and T-arm. During tRNA maturation, different modifying enzymes posttranscriptionally attach chemical groups to specific nucleotides, which not only affect translation elongation rates but also restrict local folding processes and confer local flexibility when required. The characteristic structural features of tRNAs are also employed by various maturation factors and modification enzymes to assure the selection, recognition, and positioning of specific sites within the substrate tRNAs. The cellular functional repertoire of tRNAs continues to extend well beyond their role in translation, partly, due to the expanding pool of tRNA-derived fragments. Here, we aim to summarize the most recent developments in the field to understand how three-dimensional structure affects the canonical and noncanonical functions of tRNA.

Generating Mitochondrial-Nuclear Exchange (MNX) Mice to Identify Mitochondrial Determinants of Cancer Metastasis

Understanding the contributions of mitochondrial genetics to disease pathogenesis is facilitated by a new and unique model-the mitochondrial-nuclear exchange mouse. Here we report the rationale for their development, the methods used to create them, and a brief summary of how MNX mice have been used to understand the contributions of mitochondrial DNA in multiple diseases, focusing on cancer metastasis. Polymorphisms in mtDNA which distinguish mouse strains exert intrinsic and extrinsic effects on metastasis efficiency by altering epigenetic marks in the nuclear genome, changing production of reactive oxygen species, altering the microbiota, and influencing immune responses to cancer cells. Although the focus of this report is cancer metastasis, MNX mice have proven to be valuable in studying mitochondrial contributions to other diseases as well.

The Evolutionary Dynamics of the Mitochondrial tRNA in the Cichlid Fish Family

Simple Summary

Cichlids are a unique example of fish diversity and species richness which have been explained by sympatric speciation at different freshwater sources in Africa. The mitochondria contribute to cell vitality by providing energy. It contains a circular genome with an established translation system that is spatially independent of the cytosolic counterpart. The current study aimed to investigate the evolutionary dynamics of the mitochondrial tRNA and its role in cichlids’ diversity. The available cichlid mitogenomes in the public database were filtered, in addition to newly sequenced accessions from a specific cichlid group known as the haplotilapiine lineage that is widely distributed in the Egyptian sector of the Nile River. Based on the comparative analysis of mitogenomic data, we identified 22 tRNA genes, in which a single gene was D-armless, while the cloverleaf secondary structure subdivided into stem-loop formations was predicted and used to define the levels of genetic divergence for the remained tRNAs. Peculiarly, in cichlids, the formation known as “T-arm” showed the lowest polymorphism levels among other structures in contrast to other organisms (e.g., scorpions). Comparing the whole family to the specific haplotilapiine lineage showed that the tryptophan tRNA was the most conserved tRNA, with signatures of possible purifying selection.

Abstract

The mitochondrial transfer RNA genes (tRNAs) attract more attention due to their highly dynamic and rapidly evolving nature. The current study aimed to detect and evaluate the dynamics, characteristic patterns, and variations of mitochondrial tRNAs. The study was conducted in two main parts: first, the published mitogenomic sequences of cichlids mt tRNAs have been filtered. Second, the filtered mitochondrial tRNA and additional new mitogenomes representing the most prevalent Egyptian tilapiine were compared and analyzed. Our results revealed that all 22 tRNAs of cichlids folded into a classical cloverleaf secondary structure with four domains, except for trnS^GCU, missing the D domain in all cichlids. When consensus tRNAs were compared, most of the mutations were observed in the trnP at nucleotide levels (substitutions and indels), in contrast to trnL^UAA. From a structural perspective, the anticodon loop and T-loop formations were the most conserved structures among all parts of the tRNA in contrast to the A-stem and D-loop formations. The trnW was the lowest polymorphic unneutral tRNA among all cichlids (both the family and the haplotilapiine lineage), in contrast with the neutral trnD that was extremely polymorphic among and within the haplotilapiine lineage species compared to other cichlids species. From a phylogenetic perspective, the trnC was extremely hypervariable and neutral tRNA in both haplotilapiine lineage and cichlids but was unable to report correct phylogenetic signal for the cichlids. In contrast to trnI and trnY, less variable neutral tRNAs that were able to cluster the haplotilapiine lineage and cichlids species as previously reported. By observing the DNA polymorphism in the coding DNA sequences (CDS), the highest affected amino acid by non-synonymous mutations was isoleucine and was equally mutated to valine and vice versa; no correlation between mutations in CDS and tRNAs was statistically found. The current study provides an insight into the mitochondrial tRNA evolution and its effect on the cichlid diversity and speciation model at the maternal level.

The tRNA-Derived Fragment tRF-24-V29K9UV3IU Functions as a miRNA-like RNA to Prevent Gastric Cancer Progression by Inhibiting GPR78 Expression

Emerging studies have proved that tRNA-derived fragments (tRFs) play vital roles in tumor metastasis; however, the function of tRFs in gastric cancer (GC) remains largely unclear. We investigated the role of tRF-24-V29K9UV3IU in growth and metastasis of GC using a xenograft mouse model. Differential gene expression downstream of tRF-24-V29K9UV3IU was identified by transcriptome sequencing, and interaction was then verified by a dual luciferase reporter and RNA immunoprecipitation. MKN-45 cells were also used to explore the biological functions of tRF-24-V29K9UV3IU in vitro. Here, knockdown of tRF-24-V29K9UV3IU promoted tumor growth and metastasis of GC in vivo. The expression of tRF-24-V29K9UV3IU and E-cadherin (epithelial cell marker) was down-regulated in tumors of mice following tRF-24-V29K9UV3IU knockdown, whereas the mesenchymal cell markers N-cadherin and vimentin displayed an opposite trend. Transcriptome sequencing identified 87 differentially expressed genes (DEGs) down-regulated in the tRF-24-V29K9UV3IU-overexpressed groups compared with the control group. Among them, G-protein–coupled receptor 78 (GPR78), the most significantly down-regulated DEG, was also predicted to be a target of tRF-24-V29K9UV3IU. Moreover, tRF-24-V29K9UV3IU could function as a miRNA-like fragment and bind to AGO2 and directly silence GPR78 expression by complementing with the 3′-untranslated region of the GPR78 mRNA. Functionally, overexpression of tRF-24-V29K9UV3IU significantly suppressed proliferation, migration, and invasion and promoted apoptosis of MKN-45 cells, whereas GPR78 attenuated these effects. Therefore, our data suggest that tRF-24-V29K9UV3IU functions as a miRNA-like fragment to suppress GPR78 expression and thus inhibit GC progression. These observations suggest that the tRF-24-V29K9UV3IU/GPR78 axis serves as a potential therapeutic target in GC.

Bray P, Pliatsika V, Telonis AG, Rigoutsos I. IsoMiRmap: fast, deterministic and exhaustive mining of isomiRs from short RNA-seq datasets

MOTIVATION

MicroRNA (miRNA) precursor arms give rise to multiple isoforms simultaneously called 'isomiRs.' IsomiRs from the same arm typically differ by a few nucleotides at either their 5' or 3' termini or both. In humans, the identities and abundances of isomiRs depend on a person's sex and genetic ancestry as well as on tissue type, tissue state and disease type/subtype. Moreover, nearly half of the time the most abundant isomiR differs from the miRNA sequence found in public databases. Accurate mining of isomiRs from deep sequencing data is thus important.

RESULTS

We developed isoMiRmap, a fast, standalone, user-friendly mining tool that identifies and quantifies all isomiRs by directly processing short RNA-seq datasets. IsoMiRmap is a portable 'plug-and-play' tool, requires minimal setup, has modest computing and storage requirements, and can process an RNA-seq dataset with 50 million reads in just a few minutes on an average laptop. IsoMiRmap deterministically and exhaustively reports all isomiRs in a given deep sequencing dataset and quantifies them accurately (no double-counting). IsoMiRmap comprehensively reports all miRNA precursor locations from which an isomiR may be transcribed, tags as 'ambiguous' isomiRs whose sequences exist both inside and outside of the space of known miRNA sequences and reports the public identifiers of common single-nucleotide polymorphisms and documented somatic mutations that may be present in an isomiR. IsoMiRmap also identifies isomiRs with 3' non-templated post-transcriptional additions. Compared to similar tools, isoMiRmap is the fastest, reports more bona fide isomiRs, and provides the most comprehensive information related to an isomiR's transcriptional origin.

AVAILABILITY AND IMPLEMENTATION

The codes for isoMiRmap are freely available at https://cm.jefferson.edu/isoMiRmap/ and https://github.com/TJU-CMC-Org/isoMiRmap/. IsomiR profiles for the datasets of the 1000 Genomes Project, spanning five population groups, and The Cancer Genome Atlas (TCGA), spanning 33 cancer studies, are also available at https://cm.jefferson.edu/isoMiRmap/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Combining tRNA sequencing methods to characterize plant tRNA expression and post-transcriptional modification

Differences in tRNA expression have been implicated in a remarkable number of biological processes. There is growing evidence that tRNA genes can play dramatically different roles depending on both expression and post-transcriptional modification, yet sequencing tRNAs to measure abundance and detect modifications remains challenging. Their secondary structure and extensive post-transcriptional modifications interfere with RNA-seq library preparation methods and have limited the utility of high-throughput sequencing technologies. Here, we combine two modifications to standard RNA-seq methods by treating with the demethylating enzyme AlkB and ligating with tRNA-specific adapters in order to sequence tRNAs from four species of flowering plants, a group that has been shown to have some of the most extensive rates of post-transcriptional tRNA modifications. This protocol has the advantage of detecting full-length tRNAs and sequence variants that can be used to infer many post-transcriptional modifications. We used the resulting data to produce a modification index of almost all unique reference tRNAs in Arabidopsis thaliana, which exhibited many anciently conserved similarities with humans but also positions that appear to be 'hot spots' for modifications in angiosperm tRNAs. We also found evidence based on northern blot analysis and droplet digital PCR that, even after demethylation treatment, tRNA-seq can produce highly biased estimates of absolute expression levels most likely due to biased reverse transcription. Nevertheless, the generation of full-length tRNA sequences with modification data is still promising for assessing differences in relative tRNA expression across treatments, tissues or subcellular fractions and help elucidate the functional roles of tRNA modifications.

Roles of mitochondria in the hallmarks of metastasis

Although mitochondrial contributions to cancer have been recognised for approximately a century, given that mitochondrial DNA (mtDNA) is dwarfed by the size of the nuclear genome (nDNA), nuclear genetics has represented a focal point in cancer biology, often at the expense of mtDNA and mitochondria. However, genomic sequencing and advances in in vivo models underscore the importance of mtDNA and mitochondria in cancer and metastasis. In this review, we explore the roles of mitochondria in the four defined 'hallmarks of metastasis': motility and invasion, microenvironment modulation, plasticity and colonisation. Biochemical processes within the mitochondria of both cancer cells and the stromal cells with which they interact are critical for each metastatic hallmark. We unravel complex dynamics in mitochondrial contributions to cancer, which are context-dependent and capable of either promoting metastasis or being leveraged to prevent it at various points of the metastatic cascade. Ultimately, mitochondrial contributions to cancer and metastasis are rooted in the capacity of these organelles to tune metabolic and genetic responses to dynamic microenvironmental cues.

Intronic tRNAs of mitochondrial origin regulate constitutive and alternative splicing

BACKGROUND

The presence of nuclear mitochondrial DNA (numtDNA) has been reported within several nuclear genomes. Next to mitochondrial protein-coding genes, numtDNA sequences also encode for mitochondrial tRNA genes. However, the biological roles of numtDNA remain elusive.

RESULTS

Employing in silico analysis, we identify 281 mitochondrial tRNA homologs in the human genome, which we term nimtRNAs (nuclear intronic mitochondrial-derived tRNAs), being contained within introns of 76 nuclear host genes. Despite base changes in nimtRNAs when compared to their mtRNA homologs, a canonical tRNA cloverleaf structure is maintained. To address potential functions of intronic nimtRNAs, we insert them into introns of constitutive and alternative splicing reporters and demonstrate that nimtRNAs promote pre-mRNA splicing, dependent on the number and positioning of nimtRNA genes and splice site recognition efficiency. A mutational analysis reveals that the nimtRNA cloverleaf structure is required for the observed splicing increase. Utilizing a CRISPR/Cas9 approach, we show that a partial deletion of a single endogenous nimtRNALys within intron 28 of the PPFIBP1 gene decreases inclusion of the downstream-located exon 29 of the PPFIBP1 mRNA. By employing a pull-down approach followed by mass spectrometry, a 3'-splice site-associated protein network is identified, including KHDRBS1, which we show directly interacts with nimtRNATyr by an electrophoretic mobility shift assay.

CONCLUSIONS

We propose that nimtRNAs, along with associated protein factors, can act as a novel class of intronic splicing regulatory elements in the human genome by participating in the regulation of splicing.

On the expanding roles of tRNA fragments in modulating cell behavior

The fragments that derive from transfer RNAs (tRNAs) are an emerging category of regulatory RNAs. Known as tRFs, these fragments were reported for the first time only a decade ago, making them a relatively recent addition to the ever-expanding pantheon of non-coding RNAs. tRFs are short, 16-35 nucleotides (nts) in length, and produced through cleavage of mature and precursor tRNAs at various positions. Both cleavage positions and relative tRF abundance depend strongly on context, including the tissue type, tissue state, and disease, as well as the sex, population of origin, and race/ethnicity of an individual. These dependencies increase the urgency to understand the regulatory roles of tRFs. Such efforts are gaining momentum, and comprise experimental and computational approaches. System-level studies across many tissues and thousands of samples have produced strong evidence that tRFs have important and multi-faceted roles. Here, we review the relevant literature on tRF biology in higher organisms, single cell eukaryotes, and prokaryotes.

tRNA Fragments Populations Analysis in Mutants Affecting tRNAs Processing and tRNA Methylation

tRNA fragments (tRFs) are a class of small non-coding RNAs (sncRNAs) derived from tRNAs. tRFs are highly abundant in many cell types including stem cells and cancer cells, and are found in all domains of life. Beyond translation control, tRFs have several functions ranging from transposon silencing to cell proliferation control. However, the analysis of tRFs presents specific challenges and their biogenesis is not well understood. They are very heterogeneous and highly modified by numerous post-transcriptional modifications. Here we describe a bioinformatic pipeline (tRFs-Galaxy) to study tRFs populations and shed light onto tRNA fragments biogenesis in Drosophila melanogaster. Indeed, we used small RNAs Illumina sequencing datasets extracted from wild type and mutant ovaries affecting two different highly conserved steps of tRNA biogenesis: 5'pre-tRNA processing (RNase-P subunit Rpp30) and tRNA 2'-O-methylation (dTrm7_34 and dTrm7_32). Using our pipeline, we show how defects in tRNA biogenesis affect nuclear and mitochondrial tRFs populations and other small non-coding RNAs biogenesis, such as small nucleolar RNAs (snoRNAs). This tRF analysis workflow will advance the current understanding of tRFs biogenesis, which is crucial to better comprehend tRFs roles and their implication in human pathology.

tRNA-Derived Small RNAs and Their Potential Roles in Cardiac Hypertrophy

Transfer RNAs (tRNAs) are abundantly expressed, small non-coding RNAs that have long been recognized as essential components of the protein translation machinery. The tRNA-derived small RNAs (tsRNAs), including tRNA halves (tiRNAs), and tRNA fragments (tRFs), were unexpectedly discovered and have been implicated in a variety of important biological functions such as cell proliferation, cell differentiation, and apoptosis. Mechanistically, tsRNAs regulate mRNA destabilization and translation, as well as retro-element reverse transcriptional and post-transcriptional processes. Emerging evidence has shown that tsRNAs are expressed in the heart, and their expression can be induced by pathological stress, such as hypertrophy. Interestingly, cardiac pathophysiological conditions, such as oxidative stress, aging, and metabolic disorders can be viewed as inducers of tsRNA biogenesis, which further highlights the potential involvement of tsRNAs in these conditions. There is increasing enthusiasm for investigating the molecular and biological functions of tsRNAs in the heart and their role in cardiovascular disease. It is anticipated that this new class of small non-coding RNAs will offer new perspectives in understanding disease mechanisms and may provide new therapeutic targets to treat cardiovascular disease.

Transfer RNA-Derived Small RNAs: Another Layer of Gene Regulation and Novel Targets for Disease Therapeutics

Decades after identification as essential for protein synthesis, transfer RNAs (tRNAs) have been implicated in various cellular processes beyond translation. tRNA-derived small RNAs (tsRNAs), referred to as tRNA-derived fragments (tRFs) or tRNA-derived, stress-induced RNAs (tiRNAs), are produced by cleavage at different sites from mature or pre-tRNAs. They are classified into six major types representing potentially thousands of unique sequences and have been implicated to play a wide variety of regulatory roles in maintaining normal homeostasis, cancer cell viability, tumorigenesis, ribosome biogenesis, chromatin remodeling, translational regulation, intergenerational inheritance, retrotransposon regulation, and viral replication. However, the detailed mechanisms governing these processes remain unknown. Aberrant expression of tsRNAs is found in various human disease conditions, suggesting that a further understanding of the regulatory role of tsRNAs will assist in identifying novel biomarkers, potential therapeutic targets, and gene-regulatory tools. Here, we highlight the classification, biogenesis, and biological role of tsRNAs in regulatory mechanisms of normal and disease states.

OncotRF: an online resource for exploration of tRNA-derived fragments in human cancers

Transfer RNA-derived fragments (tRFs) are a new class of small non-coding RNAs whose biological roles in cancers are not well understood. Emerging evidence suggests that tRFs are involved in gene regulation at multiple levels. In this study, we constructed an integrative database, OncotRF (http://bioinformatics.zju.edu.cn/OncotRF), for in silico exploration of tRF functions, and identification of diagnostic and prognostic biomarkers in cancers. The database contains an analysis pipeline for tRF identification and characterization, analysis results of 11,211 small RNA sequencing samples and 8,776 RNA sequencing samples, and clinicopathologic annotation data from The Cancer Genome Atlas (TCGA). The results include: tRF identification and quantification across 33 cancers, abnormally expressed tRFs and genes, tRF-gene correlations, tRF-gene networks, survival analyses, and tRF-related functional enrichment analyses. Users are also able to identify differentially expressed tRFs, predict their functions, and assess the relevance of the tRF expression levels to the clinical outcome according to user-defined groups. Additionally, an online Kaplan-Meier plotter is available in OncotRF for plotting survival curves according to user-defined groups. OncotRF will be a valuable online database and functional annotation tool for researchers studying the roles, functions, and mechanisms of tRFs in human cancers.

"Mitochondrial Toolbox" - A Review of Online Resources to Explore Mitochondrial Genomics

Mitochondria play a significant role in many biological systems. There is emerging evidence that differences in the mitochondrial genome may contribute to multiple common diseases, leading to an increasing number of studies exploring mitochondrial genomics. There is often a large amount of complex data generated (for example via next generation sequencing), which requires optimised bioinformatics tools to efficiently and effectively generate robust outcomes from these large datasets. Twenty-four online resources dedicated to mitochondrial genomics were reviewed. This 'mitochondrial toolbox' summary resource will enable researchers to rapidly identify the resource(s) most suitable for their needs. These resources fulfil a variety of functions, with some being highly specialised. No single tool will provide all users with the resources they require; therefore, the most suitable tool will vary between users depending on the nature of the work they aim to carry out. Genetics resources are well established for phylogeny and DNA sequence changes, but further epigenetic and gene expression resources need to be developed for mitochondrial genomics.

Ribosomal RNA fragmentation into short RNAs (rRFs) is modulated in a sex- and population of origin-specific manner

Background

The advent of next generation sequencing (NGS) has allowed the discovery of short and long non-coding RNAs (ncRNAs) in an unbiased manner using reverse genetics approaches, enabling the discovery of multiple categories of ncRNAs and characterization of the way their expression is regulated. We previously showed that the identities and abundances of microRNA isoforms (isomiRs) and transfer RNA-derived fragments (tRFs) are tightly regulated, and that they depend on a person’s sex and population origin, as well as on tissue type, tissue state, and disease type. Here, we characterize the regulation and distribution of fragments derived from ribosomal RNAs (rRNAs). rRNAs form a group that includes four (5S, 5.8S, 18S, 28S) rRNAs encoded by the human nuclear genome and two (12S, 16S) by the mitochondrial genome. rRNAs constitute the most abundant RNA type in eukaryotic cells.

Results

We analyzed rRNA-derived fragments (rRFs) across 434 transcriptomic datasets obtained from lymphoblastoid cell lines (LCLs) derived from healthy participants of the 1000 Genomes Project. The 434 datasets represent five human populations and both sexes. We examined each of the six rRNAs and their respective rRFs, and did so separately for each population and sex. Our analysis shows that all six rRNAs produce rRFs with unique identities, normalized abundances, and lengths. The rRFs arise from the 5′-end (5′-rRFs), the interior (i-rRFs), and the 3′-end (3′-rRFs) or straddle the 5′ or 3′ terminus of the parental rRNA (x-rRFs). Notably, a large number of rRFs are produced in a population-specific or sex-specific manner. Preliminary evidence suggests that rRF production is also tissue-dependent. Of note, we find that rRF production is not affected by the identity of the processing laboratory or the library preparation kit.

Conclusions

Our findings suggest that rRFs are produced in a regimented manner by currently unknown processes that are influenced by both ubiquitous as well as population-specific and sex-specific factors. The properties of rRFs mirror the previously reported properties of isomiRs and tRFs and have implications for the study of homeostasis and disease.

Exploration of CCA-added RNAs revealed the expression of mitochondrial non-coding RNAs regulated by CCA-adding enzyme

Post-transcriptional non-template additions of nucleotides to 3'-ends of RNAs play important roles in the stability and function of RNA molecules. Although tRNA nucleotidyltransferase (CCA-adding enzyme) is known to add CCA trinucleotides to 3'-ends of tRNAs, whether other RNA species can be endogenous substrates of CCA-adding enzyme has not been widely explored yet. Herein, we used YAMAT-seq to identify non-tRNA substrates of CCA-adding enzyme. YAMAT-seq captures RNA species that form secondary structures with 4-nt protruding 3'-ends of the sequence 5'-NCCA-3', which is the hallmark structure of RNAs that are generated by CCA-adding enzyme. By executing YAMAT-seq for human breast cancer cells and mining the sequence data, we identified novel candidate substrates of CCA-adding enzyme. These included fourteen 'CCA-RNAs' that only contain CCA as non-genomic sequences, and eleven 'NCCA-RNAs' that contain CCA and other nucleotides as non-genomic sequences. All newly-identified (N)CCA-RNAs were derived from the mitochondrial genome and were localized in mitochondria. Knockdown of CCA-adding enzyme severely reduced the expression levels of (N)CCA-RNAs, suggesting that the CCA-adding enzyme-catalyzed CCA additions stabilize the expression of (N)CCA-RNAs. Furthermore, expression levels of (N)CCA-RNAs were severely reduced by various cellular treatments, including UV irradiation, amino acid starvation, inhibition of mitochondrial respiratory complexes, and inhibition of the cell cycle. These results revealed a novel CCA-mediated regulatory pathway for the expression of mitochondrial non-coding RNAs.

Short RNA regulators: the past, the present, the future, and implications for precision medicine and health disparities

We herein provide a brief review of the trajectory that the field of short RNA research followed in the last 25 years. We place emphasis on the unexpected discoveries and the ramifications of these discoveries for the field, as well as offer some thoughts about what the next 25 years may bring. Arguably, the uncovered dependence of different types of short RNAs on individual attributes such as a person's sex, population origin, race, and on tissue type, tissue state, and disease was most unexpected. This dependence has important ramifications in that it will provide a boost to our understanding of the molecular mechanisms of health disparities as well as pave the way for novel approaches to designing improved and personalized diagnostics and therapeutics.

A threonyl-tRNA synthetase-like protein has tRNA aminoacylation and editing activities

TARS and TARS2 encode cytoplasmic and mitochondrial threonyl-tRNA synthetases (ThrRSs) in mammals, respectively. Interestingly, in higher eukaryotes, a third gene, TARSL2, encodes a ThrRS-like protein (ThrRS-L), which is highly homologous to cytoplasmic ThrRS but with a different N-terminal extension (N-extension). Whether ThrRS-L has canonical functions is unknown. In this work, we studied the organ expression pattern, cellular localization, canonical aminoacylation and editing activities of mouse ThrRS-L (mThrRS-L). Tarsl2 is ubiquitously but unevenly expressed in mouse tissues. Different from mouse cytoplasmic ThrRS (mThrRS), mThrRS-L is located in both the cytoplasm and nucleus; the nuclear distribution is mediated via a nuclear localization sequence at its C-terminus. Native mThrRS-L enriched from HEK293T cells was active in aminoacylation and editing. To investigate the in vitro catalytic properties of mThrRS-L accurately, we replaced the N-extension of mThrRS-L with that of mThrRS. The chimeric protein (mThrRS-L-N^T) has amino acid activation, aminoacylation and editing activities. We compared the activities and cross-species tRNA recognition between mThrRS-L-N^T and mThrRS. Despite having a similar aminoacylation activity, mThrRS-L-N^T and mThrRS exhibit differences in tRNA recognition and editing capacity. Our results provided the first analysis of the aminoacylation and editing activities of ThrRS-L, and improved our understanding of Tarsl2.

Accurate Profiling and Quantification of tRNA Fragments from RNA-Seq Data: A Vade Mecum for MINTmap

There is an increasing interest within the scientific community in identifying tRNA-derived fragments (tRFs) and elucidating the roles they play in the cell. Such endeavors can be greatly facilitated by mining the numerous datasets from many cellular contexts that exist publicly. However, the standard mapping tools cannot be used for the purpose. Several factors complicate this endeavor including: the presence of multiple identical or nearly identical isodecoders at various genomic locations; the presence of identical sequence segments that are shared by isodecoders of the same or even different anticodons; the existence of numerous partial tRNA sequences across the genome; the existence of hundreds of "lookalike" sequences that resemble true tRNAs; and others. This is generating a need for specialized tools that can mine deep sequencing data to identify and quantify tRFs. We discuss the various complicating factors and their ramifications, and how to use and run MINTmap, a tool that addresses these considerations.

tRex: A Web Portal for Exploration of tRNA-Derived Fragments in Arabidopsis thaliana

tRNA-derived fragments (tRFs) constitute a new class of short regulatory RNAs that are a product of nascent or mature tRNA processing. tRF sequences have been identified in all domains of life; however, most published research pertains to human, yeast and some bacterial organisms. Despite growing interest in plant tRFs and accumulating evidence of their function in plant development and stress responses, no public, web-based repository dedicated to these molecules is currently available. Here, we introduce tRex (http://combio.pl/trex)-the first comprehensive data-driven online resource specifically dedicated to tRFs in the model plant Arabidopsis thaliana. The portal is based on verified Arabidopsis tRNA annotation and includes in-house-generated and publicly available small RNA sequencing experiments from various tissues, ecotypes, genotypes and stress conditions. The provided web-based tools are designed in a user-friendly manner and allow for seamless exploration of the data that are presented in the form of dynamic tables and cumulative coverage profiles. The tRex database is connected to external genomic and citation resources, which makes it a one-stop solution for Arabidopsis tRF-related research.

YAMAT-seq: an efficient method for high-throughput sequencing of mature transfer RNAs

Besides translation, transfer RNAs (tRNAs) play many non-canonical roles in various biological pathways and exhibit highly variable expression profiles. To unravel the emerging complexities of tRNA biology and molecular mechanisms underlying them, an efficient tRNA sequencing method is required. However, the rigid structure of tRNA has been presenting a challenge to the development of such methods. We report the development of Y-shaped Adapter-ligated MAture TRNA sequencing (YAMAT-seq), an efficient and convenient method for high-throughput sequencing of mature tRNAs. YAMAT-seq circumvents the issue of inefficient adapter ligation, a characteristic of conventional RNA sequencing methods for mature tRNAs, by employing the efficient and specific ligation of Y-shaped adapter to mature tRNAs using T4 RNA Ligase 2. Subsequent cDNA amplification and next-generation sequencing successfully yield numerous mature tRNA sequences. YAMAT-seq has high specificity for mature tRNAs and high sensitivity to detect most isoacceptors from minute amount of total RNA. Moreover, YAMAT-seq shows quantitative capability to estimate expression levels of mature tRNAs, and has high reproducibility and broad applicability for various cell lines. YAMAT-seq thus provides high-throughput technique for identifying tRNA profiles and their regulations in various transcriptomes, which could play important regulatory roles in translation and other biological processes.

MINTmap: fast and exhaustive profiling of nuclear and mitochondrial tRNA fragments from short RNA-seq data

Transfer RNA fragments (tRFs) are an established class of constitutive regulatory molecules that arise from precursor and mature tRNAs. RNA deep sequencing (RNA-seq) has greatly facilitated the study of tRFs. However, the repeat nature of the tRNA templates and the idiosyncrasies of tRNA sequences necessitate the development and use of methodologies that differ markedly from those used to analyze RNA-seq data when studying microRNAs (miRNAs) or messenger RNAs (mRNAs). Here we present MINTmap (for MItochondrial and Nuclear TRF mapping), a method and a software package that was developed specifically for the quick, deterministic and exhaustive identification of tRFs in short RNA-seq datasets. In addition to identifying them, MINTmap is able to unambiguously calculate and report both raw and normalized abundances for the discovered tRFs. Furthermore, to ensure specificity, MINTmap identifies the subset of discovered tRFs that could be originating outside of tRNA space and flags them as candidate false positives. Our comparative analysis shows that MINTmap exhibits superior sensitivity and specificity to other available methods while also being exceptionally fast. The MINTmap codes are available through https://github.com/TJU-CMC-Org/MINTmap/ under an open source GNU GPL v3.0 license.

Noncoding RNA: Current Deep Sequencing Data Analysis Approaches and Challenges

One of the most significant biological discoveries of the last decade is represented by the reality that the vast majority of the transcribed genomic output comprises diverse classes of noncoding RNAs (ncRNAs) that may play key roles and/or be affected by many biochemical cellular processes (i.e., RNA editing), with implications in human health and disease. With 90% of the human genome being transcribed and novel classes of ncRNA emerging (tRNA-derived small RNAs and circular RNAs among others), the great majority of the human transcriptome suggests that many important ncRNA functions/processes are yet to be discovered. An approach to filling such vast void of knowledge has been recently provided by the increasing application of next-generation sequencing (NGS), offering the unprecedented opportunity to obtain a more accurate profiling with higher resolution, increased throughput, sequencing depth, and low experimental complexity, concurrently posing an increasing challenge in terms of efficiency, accuracy, and usability of data analysis software. This review provides an overview of ncRNAs, NGS technology, and the most recent/popular computational approaches and the challenges they attempt to solve, which are essential to a more sensitive and comprehensive ncRNA annotation capable of furthering our understanding of this still vastly uncharted genomic territory.

Dissecting tRNA-derived fragment complexities using personalized transcriptomes reveals novel fragment classes and unexpected dependencies

We analyzed transcriptomic data from 452 healthy men and women representing five different human populations and two races, and, 311 breast cancer samples from The Cancer Genome Atlas. Our studies revealed numerous constitutive, distinct fragments with overlapping sequences and quantized lengths that persist across dozens of individuals and arise from the genomic loci of all nuclear and mitochondrial human transfer RNAs (tRNAs). Surprisingly, we discovered that the tRNA fragments' length, starting and ending points, and relative abundance depend on gender, population, race and also on amino acid identity, anticodon, genomic locus, tissue, disease, and disease subtype. Moreover, the length distribution of mitochondrially-encoded tRNAs differs from that of nuclearly-encoded tRNAs, and the specifics of these distributions depend on tissue. Notably, tRNA fragments from the same anticodon do not have correlated abundances. We also report on a novel category of tRNA fragments that significantly contribute to the differences we observe across tissues, genders, populations, and races: these fragments, referred to as i-tRFs, are abundant in human tissues, wholly internal to the respective mature tRNA, and can straddle the anticodon. HITS-CLIP data analysis revealed that tRNA fragments are loaded on Argonaute in a cell-dependent manner, suggesting cell-dependent functional roles through the RNA interference pathway. We validated experimentally two i-tRF molecules: the first was found in 21 of 22 tested breast tumor and adjacent normal samples and was differentially abundant between health and disease whereas the second was found in all eight tested breast cancer cell lines.

MINTbase: a framework for the interactive exploration of mitochondrial and nuclear tRNA fragments

Motivation: It has been known that mature transfer RNAs (tRNAs) that are encoded in the nuclear genome give rise to short molecules, collectively known as tRNA fragments or tRFs. Recently, we reported that, in healthy individuals and in patients, tRFs are constitutive, arise from mitochondrial as well as from nuclear tRNAs, and have composition and abundances that depend on a person’s sex, population origin and race as well as on tissue, disease and disease subtype. Our findings as well as similar work by other groups highlight the importance of tRFs and presage an increase in the community’s interest in elucidating the roles of tRFs in health and disease.

Results: We created MINTbase, a web-based framework that serves the dual-purpose of being a content repository for tRFs and a tool for the interactive exploration of these newly discovered molecules. A key feature of MINTbase is that it deterministically and exhaustively enumerates all possible genomic locations where a sequence fragment can be found and indicates which fragments are exclusive to tRNA space, and thus can be considered as tRFs: this is a very important consideration given that the genomes of higher organisms are riddled with partial tRNA sequences and with tRNA-lookalikes whose aberrant transcripts can be mistaken for tRFs. MINTbase is extremely flexible and integrates and presents tRF information from multiple yet interconnected vantage points (‘vistas’). Vistas permit the user to interactively personalize the information that is returned and the manner in which it is displayed. MINTbase can report comparative information on how a tRF is distributed across all anticodon/amino acid combinations, provides alignments between a tRNA and multiple tRFs with which the user can interact, provides details on published studies that reported a tRF as expressed, etc. Importantly, we designed MINTbase to contain all possible tRFs that could ever be produced by mature tRNAs: this allows us to report on their genomic distributions, anticodon/amino acid properties, alignments, etc. while giving users the ability to at-will investigate candidate tRF molecules before embarking on focused experimental explorations. Lastly, we also introduce a new labeling scheme that is tRF-sequence-based and allows users to associate a tRF with a universally unique label (‘tRF-license plate’) that is independent of a genome assembly and does not require any brokering mechanism.

Availability and Implementation: MINTbase is freely accessible at http://cm.jefferson.edu/MINTbase/. Dataset submissions to MINTbase can be initiated at http://cm.jefferson.edu/MINTsubmit/.

Contact:isidore.rigoutsos@jefferson.edu

Supplementary information: Supplementary data are available at Bioinformatics online.

Consequential considerations when mapping tRNA fragments

We examine several of the choices that went into the design of tDRmapper, a recently reported tool for identifying transfer RNA (tRNA) fragments in deep sequencing data, evaluate them in the context of currently available knowledge, and discuss their potential impact on the output that the tool generates.

Four-leaf clover qRT-PCR: A convenient method for selective quantification of mature tRNA

Transfer RNAs (tRNAs) play a central role in translation and also recently appear to have a variety of other functions in biological processes beyond translation. Here we report the development of Four-Leaf clover qRT-PCR (FL-PCR), a convenient PCR-based method, which can specifically quantify individual mature tRNA species. In FL-PCR, T4 RNA ligase 2 specifically ligates a stem-loop adapter to mature tRNAs but not to precursor tRNAs or tRNA fragments. Subsequent TaqMan qRT-PCR amplifies only unmodified regions of the tRNA-adapter ligation products; therefore, FL-PCR quantification is not influenced by tRNA post-transcriptional modifications. FL-PCR has broad applicability for the quantification of various tRNAs in different cell types, and thus provides a much-needed simple method for analyzing tRNA abundance and heterogeneity.

Mitochondrial tRNA-lookalikes in nuclear chromosomes: could they be functional?

The presence in human nuclear chromosomes of multiple sequences that are highly similar to human mitochondrial tRNAs (tRNA-lookalikes) raises intriguing questions about the possible functionality of these genomic loci. In this perspective, we explore the significance of the mitochondrial tRNA-lookalikes based on a series of properties that argue for their non-accidental nature. We particularly focus on the possibility of transcription as well as on potential functional roles for these sequences that can range from their acting as DNA regulatory elements to forming functional mature tRNAs or tRNA-derived fragments. Extension of our analysis to other simians (chimp, gorilla, rhesus, and squirrel monkey), 2 rodents (mouse and rat), a marsupial (opossum) and 3 invertebrates (fruit-fly, worm, and sponge) revealed that mitochondrial tRNA-lookalikes are prevalent in primates and the opossum but absent from the other analyzed organisms.

tRNA-Derived Short Non-coding RNA as Interacting Partners of Argonaute Proteins

The advent of next-generation sequencing technologies has not only accelerated findings on various novel non-coding RNA (ncRNA) species but also led to the revision of the biological significance and versatility of fundamental RNA species with canonical function, such as transfer RNAs (tRNAs). Although tRNAs are best known as adapter components of translational machinery, recent studies suggest that tRNAs are not always end products but can further serve as a source for short ncRNAs. In many organisms, various tRNA-derived ncRNA species are produced from mature tRNAs or their precursor transcripts as functional molecules involved in various biological processes beyond translation. In this review, we focus on the tRNA-derived ncRNAs associated with Argonaute proteins and summarize recent studies on their conceivable biogenesis factors and on their emerging roles in gene expression regulation as regulatory RNAs.

Sex hormone-dependent tRNA halves enhance cell proliferation in breast and prostate cancers

Sex hormones and their receptors play critical roles in the development and progression of the breast and prostate cancers. Here we report that a novel type of transfer RNA (tRNA)-derived small RNA, termed Sex HOrmone-dependent TRNA-derived RNAs (SHOT-RNAs), are specifically and abundantly expressed in estrogen receptor (ER)-positive breast cancer and androgen receptor (AR)-positive prostate cancer cell lines. SHOT-RNAs are not abundantly present in ER(-) breast cancer, AR(-) prostate cancer, or other examined cancer cell lines from other tissues. ER-dependent accumulation of SHOT-RNAs is not limited to a cell culture system, but it also occurs in luminal-type breast cancer patient tissues. SHOT-RNAs are produced from aminoacylated mature tRNAs by angiogenin-mediated anticodon cleavage, which is promoted by sex hormones and their receptors. Resultant 5'- and 3'-SHOT-RNAs, corresponding to 5'- and 3'-tRNA halves, bear a cyclic phosphate (cP) and an amino acid at the 3'-end, respectively. By devising a "cP-RNA-seq" method that is able to exclusively amplify and sequence cP-containing RNAs, we identified the complete repertoire of 5'-SHOT-RNAs. Furthermore, 5'-SHOT-RNA, but not 3'-SHOT-RNA, has significant functional involvement in cell proliferation. These results have unveiled a novel tRNA-engaged pathway in tumorigenesis of hormone-dependent cancers and implicate SHOT-RNAs as potential candidates for biomarkers and therapeutic targets.