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

In silico approach for the identification of tRNA-derived small non-coding RNAs in SARS-CoV infection

tsRNAs (tRNA-derived small non-coding RNAs), including tRNA halves (tiRNAs) and tRNA fragments (tRFs), have been implicated in some viral infections, such as respiratory viral infections. However, their involvement in SARS-CoV infection is completely unknown. A comprehensive analysis was performed to determine tsRNA populations in a mouse model of SARS-CoV-infected samples containing the wild-type and attenuated viruses. Data from the Gene Expression Omnibus (GEO) dataset at NCBI (accession ID GSE90624 ) was used for this study. A count matrix was generated for the tRNAs. Differentially expressed tRNAs, followed by tsRNAs derived from each significant tRNAs at different conditions and time points between the two groups WT(SARS-CoV-MA15-WT) vs Mock and ΔE (SARS-CoV-MA15-ΔE) vs Mock were identified. Notably, significantly differentially expressed tRNAs at 2dpi but not at 4dpi. The tsRNAs originating from differentially expressed tRNAs across all the samples belonging to each condition (WT, ΔE, and Mock) were identified. Intriguingly, tRFs (tRNA-derived RNA fragments) exhibited higher levels compared to tiRNAs (tRNA-derived stress-induced RNAs) across all samples associated with WT SARS-CoV strain compared to ΔE and mock-infected samples. This discrepancy suggests a non-random formation of tsRNAs, hinting at a possible involvement of tsRNAs in SARS-CoV viral infection.

A novel tRNA-derived fragment tRF-3023b suppresses inflammation in RAW264.7 cells by targeting Cul4a through NF-κB signaling

The role of transfer RNA (tRNA)-derived fragment (tRF) in various diseases has been established. However, the effect of tRF-3023b on inflammation remains unclear. Inflammation was imitated in RAW264.7 cells by adding Lipopolysaccharide (LPS). Cells were first divided into control, LPS, and LPS + Bulleyaconitine A (BLA) groups. The contents of TNF-α, IL-6, and MCP-1 were quantified using ELISA. The levels of cyclooxygenase-2 (COX2), inducible nitric oxide synthase (iNOS), and the phosphorylation of nuclear factor-kappa B (NF-κB)-P65 (p-P65) were detected by Western blotting. RNA sequencing was utilized to find differentially expressed tRFs (DE-tRFs) among three groups. The levels of various tRFs were checked by quantitative real-time PCR (qRT-PCR). Cell cycle and apoptosis were checked by flow cytometry. Dluciferase reporter assay was applied to predict and confirm the interaction between tRF-3023b and Cullin 4A (Cul4a), subsequently RNA pull-down followed by mass spectrometry analysis were conducted. BLA treatment decreased the contents of TNF-α, IL-6, MCP-1, and the expression levels of COX2, iNOS, p-P65. We found 6 DE-tRFs in LPS + BLA group compared to LPS group, tRF-3023b was high expression in control and BLA groups, and the lowest in LPS group. Cul4a was a direct target of tRF-3023b. tRF-3023b mimic affected the cell cycle distribution, promoted cells apoptosis, and suppressed the TNF-α, IL-6, MCP-1, COX2, iNOS and p-P65. The suppression of Cul4a affected the cell cycle distribution, resulted in an increase of cell apoptosis while a decrease of TNF-α, IL-6, MCP-1, COX2, iNOS and p-P65. Furthermore, Cul4a overexpression reversed the effect of tRF-3023b mimic. Cul4a knockdown reversed the effect of tRF-3023b inhibitor. Our study positions tRF-3023b as a compelling candidate, through its interaction with Cul4a, the underlying mechanism on inflammation maybe related to NF-κB pathway. The study provides a basis for exploring new therapeutic strategies for inflammation.

The role of microRNAs in responses to drought and heat stress in peanut (Arachis hypogaea)

MicroRNAs (miRNAs) are 21-24 nt small RNAs (sRNAs) that negatively regulate protein-coding genes and/or trigger phased small-interfering RNA (phasiRNA) production. Two thousand nine hundred miRNA families, of which ∼40 are deeply conserved, have been identified in ∼80 different plant species genomes. miRNA functions in response to abiotic stresses is less understood than their roles in development. Only seven peanut MIRNA families are documented in miRBase, yet a reference genome assembly is now published and over 480 plant-like MIRNA loci were predicted in the diploid peanut progenitor Arachis duranensis genome. We explored by computational analysis of a leaf sRNA library and publicly available sRNA, degradome, and transcriptome datasets the miRNA and phasiRNA space associated with drought and heat stresses in peanut. We characterized 33 novel candidate and 33 ancient conserved families of MIRNAs and present degradome evidence for their cleavage activities on mRNA targets, including several noncanonical targets and novel phasiRNA-producing noncoding and mRNA loci with validated novel targets such as miR1509 targeting serine/threonine-protein phosphatase7 and miRc20 and ahy-miR3514 targeting penta-tricopeptide repeats (PPRs), in contradistinction to other claims of miR1509/173/7122 superfamily miRNAs indirectly targeting PPRs via TAS-like noncoding RNA loci. We characterized the inverse correlations of significantly differentially expressed drought- and heat-regulated miRNAs, assayed by sRNA blots or transcriptome datasets, with target mRNA expressions in the same datasets. Meta-analysis of an expression atlas and over representation of miRNA target genes in co-expression networks suggest that miRNAs have functions in unique aspects of peanut gynophore development. Genome-wide MIRNA annotation of the published allopolyploid peanut genome can facilitate molecular breeding of value-added traits.

Differentially expressed tRNA-derived fragments in bovine fetuses with assisted reproduction induced congenital overgrowth syndrome

Background: As couples struggle with infertility and livestock producers wish to rapidly improve genetic merit in their herd, assisted reproductive technologies (ART) have become increasingly popular in human medicine as well as the livestock industry. Utilizing ART can cause an increased risk of congenital overgrowth syndromes, such as Large Offspring Syndrome (LOS) in ruminants. A dysregulation of transcripts has been observed in bovine fetuses with LOS, which is suggested to be a cause of the phenotype. Our recent study identified variations in tRNA expression in LOS individuals, leading us to hypothesize that variations in tRNA expression can influence the availability of their processed regulatory products, tRNA-derived fragments (tRFs). Due to their resemblance in size to microRNAs, studies suggest that tRFs target mRNA transcripts and regulate gene expression. Thus, we have sequenced small RNA isolated from skeletal muscle and liver of day 105 bovine fetuses to elucidate the mechanisms contributing to LOS. Moreover, we have utilized our previously generated tRNA sequencing data to analyze the contribution of tRNA availability to tRF abundance. Results: 22,289 and 7,737 unique tRFs were predicted in the liver and muscle tissue respectively. The greatest number of reads originated from 5' tRFs in muscle and 5' halves in liver. In addition, mitochondrial (MT) and nuclear derived tRF expression was tissue-specific with most MT-tRFs and nuclear tRFs derived from LysUUU and iMetCAU in muscle, and AsnGUU and GlyGCC in liver. Despite variation in tRF abundance within treatment groups, we identified differentially expressed (DE) tRFs across Control-AI, ART-Normal, and ART-LOS groups with the most DE tRFs between ART-Normal and ART-LOS groups. Many DE tRFs target transcripts enriched in pathways related to growth and development in the muscle and tumor development in the liver. Finally, we found positive correlation coefficients between tRNA availability and tRF expression in muscle (R = 0.47) and liver (0.6). Conclusion: Our results highlight the dysregulation of tRF expression and its regulatory roles in LOS. These tRFs were found to target both imprinted and non-imprinted genes in muscle as well as genes linked to tumor development in the liver. Furthermore, we found that tRNA transcription is a highly modulated event that plays a part in the biogenesis of tRFs. This study is the first to investigate the relationship between tRNA and tRF expression in combination with ART-induced LOS.

The role of extracellular vesicle miRNAs and tRNAs in synovial fibroblast senescence

Extracellular vesicles are mediators of intercellular communication with critical roles in cellular senescence and ageing. In arthritis, senescence is linked to the activation of a pro-inflammatory phenotype contributing to chronic arthritis pathogenesis. We hypothesised that senescent osteoarthritic synovial fibroblasts induce senescence and a pro-inflammatory phenotype in non-senescent osteoarthritic fibroblasts, mediated through extracellular vesicle cargo. Small RNA-sequencing and mass spectrometry proteomics were performed on extracellular vesicles isolated from the secretome of non-senescent and irradiation-induced senescent synovial fibroblasts. β-galactosidase staining confirmed senescence in SFs. RNA sequencing identified 17 differentially expressed miRNAs, 11 lncRNAs, 14 tRNAs and one snoRNA and, 21 differentially abundant proteins were identified by mass spectrometry. Bioinformatics analysis of miRNAs identified fibrosis, cell proliferation, autophagy, and cell cycle as significant pathways, tRNA analysis was enriched for signaling pathways including FGF, PI3K/AKT and MAPK, whilst protein analysis identified PAX3-FOXO1, MYC and TFGB1 as enriched upstream regulators involved in senescence and cell cycle arrest. Finally, treatment of non-senescent synovial fibroblasts with senescent extracellular vesicles confirmed the bystander effect, inducing senescence in non-senescent cells potentially through down regulation of NF-κβ and cAMP response element signaling pathways thus supporting our hypothesis. Understanding the exact composition of EV-derived small RNAs of senescent cells in this way will inform our understanding of their roles in inflammation, intercellular communication, and as active molecules in the senescence bystander effect.

An Expanded Landscape of Unusually Short RNAs in 11 Samples from Six Eukaryotic Organisms

Small RNA sequencing (sRNA-Seq) approaches unveiled sequences derived from longer non-coding RNAs, such as transfer RNA (tRNA) and ribosomal RNA (rRNA) fragments, known as tRFs and rRFs, respectively. However, rRNAs and RNAs shorter than 16 nt are often depleted from library preparations/sequencing analyses, although they may be functional. Here, we sought to obtain a complete repertoire of small RNAs by sequencing the total RNA from 11 samples of 6 different eukaryotic organisms, from yeasts to human, in an extended 8- to 30-nt window of RNA length. The 8- to 15-nt window essentially contained fragments of longer non-coding RNAs, such as microRNAs, PIWI-associated RNAs (piRNAs), small nucleolar RNAs (snoRNAs), tRNAs and rRNAs. Notably, unusually short RNAs < 16 nt were more abundant than those >16 nt in bilaterian organisms. A new RT-qPCR method confirmed that two unusually short rRFs of 12 and 13 nt were more overly abundant (~3-log difference) than two microRNAs. We propose to not deplete rRNA and to reduce the lower threshold of RNA length to include unusually short RNAs in sRNA-Seq analyses and datasets, as their abundance and diversity support their potential role and importance as biomarkers of disease and/or mediators of cellular function.

Making Invisible RNA Visible: Discriminative Sequencing Methods for RNA Molecules with Specific Terminal Formations

Next generation sequencing of RNA molecules (RNA-seq) has become a common tool to characterize the expression profiles of RNAs and their regulations in normal physiological processes and diseases. Although increasingly accumulating RNA-seq data are widely available through publicly accessible sites, most of the data for short non-coding RNAs (sncRNAs) have been obtained for microRNA (miRNA) analyses by standard RNA-seq, which only capture the sncRNAs with 5'-phosphate (5'-P) and 3'-hydroxyl (3'-OH) ends. The sncRNAs with other terminal formations such as those with a 5'-hydroxyl end (5'-OH), a 3'-phosphate (3'-P) end, or a 2',3'-cyclic phosphate end (2',3'-cP) cannot be efficiently amplified and sequenced by standard RNA-seq. Due to the invisibility in standard RNA-seq data, these non-miRNA-sncRNAs have been a hidden component in the transcriptome. However, as the functional significances of these sncRNAs have become increasingly apparent, specific RNA-seq methods compatible with various terminal formations of sncRNAs have been developed and started shedding light on the previously unrecognized sncRNAs that lack 5'-P/3'-OH ends. In this review, we summarize the expanding world of sncRNAs with various terminal formations and the strategic approaches of specific RNA-seq methods to distinctively characterize their expression profiles.

Deciphering the RNA universe in sperm in its role as a vertical information carrier

The inheritance of neurophysiologic and neuropsychologic complex diseases can only partly be explained by the Mendelian concept of genetic inheritance. Previous research showed that both psychological disorders like post-traumatic stress disorder and metabolic diseases are more prevalent in the progeny of affected parents. This could suggest an epigenetic mode of transmission. Human studies give first insight into the scope of intergenerational influence of stressors but are limited in exploring the underlying mechanisms. Animal models have elucidated the mechanistic underpinnings of epigenetic transmission. In this review, we summarize progress on the mechanisms of paternal intergenerational transmission by means of sperm RNA in mouse models. We discuss relevant details for the modelling of RNA-mediated transmission, point towards currently unanswered questions and propose experimental considerations for tackling these questions.

Transfer RNA-Derived Small RNAs in the Pathogenesis of Parasitic Protozoa

Transfer RNA (tRNA)-derived small RNAs (tsRNAs) are newly identified non-coding small RNAs that have recently attracted attention due to their functional significance in both prokaryotes and eukaryotes. tsRNAs originated from the cleavage of precursor or mature tRNAs by specific nucleases. According to the start and end sites, tsRNAs can be broadly divided into tRNA halves (31-40 nucleotides) and tRNA-derived fragments (tRFs, 14-30 nucleotides). tsRNAs have been reported in multiple organisms to be involved in gene expression regulation, protein synthesis, and signal transduction. As a novel regulator, tsRNAs have also been identified in various protozoan parasites. The conserved biogenesis of tsRNAs in early-branching eukaryotes strongly suggests the universality of this machinery, which requires future research on their shared and potentially disparate biological functions. Here, we reviewed the recent studies of tsRNAs in several representative protozoan parasites including their biogenesis and the roles in parasite biology and intercellular communication. Furthermore, we discussed the remaining questions and potential future works for tsRNAs in this group of organisms.

CLASH Analyst: A Web Server to Identify In Vivo RNA–RNA Interactions from CLASH Data

Non-coding RNAs, such as miRNAs and piRNAs, play critical roles in gene regulation through base-pairing interactions with their target molecules. The recent development of the crosslinking, ligation, and sequencing of hybrids (CLASH) method has allowed scientists to map transcriptome-wide RNA–RNA interactions by identifying chimeric reads consisting of fragments from regulatory RNAs and their targets. However, analyzing CLASH data requires scientists to use advanced bioinformatics, and currently available tools are limited for users with little bioinformatic experience. In addition, many published CLASH studies do not show the full scope of RNA–RNA interactions that were captured, highlighting the importance of reanalyzing published data. Here, we present CLASH Analyst, a web server that can analyze raw CLASH data within a fully customizable and easy-to-use interface. CLASH Analyst accepts raw CLASH data as input and identifies the RNA chimeras containing the regulatory and target RNAs according to the user’s interest. Detailed annotation of the captured RNA–RNA interactions is then presented for the user to visualize within the server or download for further analysis. We demonstrate that CLASH Analyst can identify miRNA- and piRNA-targeting sites reported from published CLASH data and should be applicable to analyze other RNA–RNA interactions. CLASH Analyst is freely available for academic use.

Transfer RNA-Derived Fragments, the Underappreciated Regulatory Small RNAs in Microbial Pathogenesis

In eukaryotic organisms, transfer RNA (tRNA)-derived fragments have diverse biological functions. Considering the conserved sequences of tRNAs, it is not surprising that endogenous tRNA fragments in bacteria also play important regulatory roles. Recent studies have shown that microbes secrete extracellular vesicles (EVs) containing tRNA fragments and that the EVs deliver tRNA fragments to eukaryotic hosts where they regulate gene expression. Here, we review the literature describing microbial tRNA fragment biogenesis and how the fragments secreted in microbial EVs suppress the host immune response, thereby facilitating chronic infection. Also, we discuss knowledge gaps and research challenges for understanding the pathogenic roles of microbial tRNA fragments in regulating the host response to infection.

Integrated genomic analysis reveals regulatory pathways and dynamic landscapes of the tRNA transcriptome

tRNAs and tRNA-derived RNA fragments (tRFs) play various roles in many cellular processes outside of protein synthesis. However, comprehensive investigations of tRNA/tRF regulation are rare. In this study, we used new algorithms to extensively analyze the publicly available data from 1332 ChIP-Seq and 42 small-RNA-Seq experiments in human cell lines and tissues to investigate the transcriptional and posttranscriptional regulatory mechanisms of tRNAs. We found that histone acetylation, cAMP, and pluripotency pathways play important roles in the regulation of the tRNA gene transcription in a cell-specific manner. Analysis of RNA-Seq data identified 950 high-confidence tRFs, and the results suggested that tRNA pools are dramatically distinct across the samples in terms of expression profiles and tRF composition. The mismatch analysis identified new potential modification sites and specific modification patterns in tRNA families. The results also show that RNA library preparation technologies have a considerable impact on tRNA profiling and need to be optimized in the future.

The tRNA-Derived Fragment-3017A Promotes Metastasis by Inhibiting NELL2 in Human Gastric Cancer

tRNA-derived fragments (tRFs) are a new classification of small non-coding RNAs (sncRNAs) derived from the specific cleavage of precursors and mature tRNAs. Accumulating recent evidence has shown that tRFs are frequently abnormal in several cancers. Nevertheless, the role of tRFs in gastric cancer and its mechanism remain unclear. In this study, we found abnormal expression of tRF-3017A (derived from tRNA-Val-TAC) in gastric cancer tissues and cell lines and confirmed its effect on promoting the invasion and migration of gastric cancer cells through functional experiments in vitro. Analysis of clinicopathologic data showed patients with higher tRF-3017A were associated with significantly higher lymph node metastasis. Mechanistic investigation implies that tRF-3017A regulates the tumor suppressor gene NELL2 through forming the RNA-induced silencing complex (RISC) with Argonaute (AGO) proteins. In this study, we found that higher tRF-3017A were associated with significantly higher lymph node metastasis in gastric cancer patients and the tRF-3017A may play a role in promoting the migration and invasion of gastric cancer cells by silencing tumor suppressor NELL2.

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.

Infection-induced 5'-half molecules of tRNAHisGUG activate Toll-like receptor 7

Toll-like receptors (TLRs) play a crucial role in the innate immune response. Although endosomal TLR7 recognizes single-stranded RNAs, their endogenous RNA ligands have not been fully explored. Here, we report 5'-tRNA half molecules as abundant activators of TLR7. Mycobacterial infection and accompanying surface TLR activation up-regulate the expression of 5'-tRNA half molecules in human monocyte-derived macrophages (HMDMs). The abundant accumulation of 5'-tRNA halves also occur in HMDM-secreted extracellular vehicles (EVs); the abundance of EV-5'-tRNAHisGUG half molecules is >200-fold higher than that of the most abundant EV-microRNA (miRNA). Sequence identification of the 5'-tRNA halves using cP-RNA-seq revealed abundant and selective packaging of specific 5'-tRNA half species into EVs. The EV-5'-tRNAHisGUG half was experimentally demonstrated to be delivered into endosomes in recipient cells and to activate endosomal TLR7. Up-regulation of the 5'-tRNA half molecules was also observed in the plasma of patients infected with Mycobacterium tuberculosis. These results unveil a novel tRNA-engaged pathway in the innate immune response and assign the role of "immune activators" to 5'-tRNA half molecules.

Integrative analysis of transcriptomes highlights potential functions of transfer-RNA-derived small RNAs in experimental intracerebral hemorrhage

Transfer-RNA-derived small RNAs (tsRNAs) are a novel class of short non-coding RNAs, that possess regulatory functions. However, their biological roles in hemorrhagic stroke are not understood. In this study, by RNA sequencing, we investigated the tsRNA expression profiles of intracerebral hemorrhagic rat brains in the chronic phase. A total of 331 tsRNAs were identified (308 in sham and 309 in intracerebral hemorrhage). Among them, the validation revealed that 7 tsRNAs (1 up-regulated and 6 down-regulated) were significantly changed. Subsequently, we predicted the target mRNAs of the 7 tsRNAs. Through integrative analysis, the predicted targets were validated by mRNA microarray data. Moreover, we confirmed the functions of tsRNAs targeting mRNAs in vitro. Furthermore, using bioinformatics tools and databases, we developed a tsRNA-mRNA-pathway interaction network to visualize their potential functions. Bioinformatics analyses and confirmatory experiments indicated that the altered genes were mainly enriched in several signaling pathways. These pathways were interrelated with intracerebral hemorrhage, such as response to oxidative stress, endocytosis, and regulation of G protein-coupled receptor signaling pathway. In summary, this study systematically revealed the profiles of tsRNAs after an experimental intracerebral hemorrhage. These results may provide novel therapeutic targets following a hemorrhagic stroke in the chronic phase.

Modulation of Small RNA Signatures in Schwann-Cell-Derived Extracellular Vesicles by the p75 Neurotrophin Receptor and Sortilin

Schwann cells (SCs) are the main glial cells of the peripheral nervous system (PNS) and are known to be involved in various pathophysiological processes, such as diabetic neuropathy and nerve regeneration, through neurotrophin signaling. Such glial trophic support to axons, as well as neuronal survival/death signaling, has previously been linked to the p75 neurotrophin receptor (p75NTR) and its co-receptor Sortilin. Recently, SC-derived extracellular vesicles (EVs) were shown to be important for axon growth and nerve regeneration, but cargo of these glial cell-derived EVs has not yet been well-characterized. In this study, we aimed to characterize signatures of small RNAs in EVs derived from wild-type (WT) SCs and define differentially expressed small RNAs in EVs derived from SCs with genetic deletions of p75NTR (Ngfr-/-) or Sortilin (Sort1-/-). Using RNA sequencing, we identified a total of 366 miRNAs in EVs derived from WT SCs of which the most highly expressed are linked to the regulation of axonogenesis, axon guidance and axon extension, suggesting an involvement of SC EVs in axonal homeostasis. Signaling of SC EVs to non-neuronal cells was also suggested by the presence of several miRNAs important for regulation of the endothelial cell apoptotic process. Ablated p75NTR or sortilin expression in SCs translated into a set of differentially regulated tRNAs and miRNAs, with impact in autophagy and several cellular signaling pathways such as the phosphatidylinositol signaling system. With this work, we identified the global expression profile of small RNAs present in SC-derived EVs and provided evidence for a regulatory function of these vesicles on the homeostasis of other cell types of the PNS. Differentially identified miRNAs can pave the way to a better understanding of p75NTR and sortilin roles regarding PNS homeostasis and disease.

A low-bias and sensitive small RNA library preparation method using randomized splint ligation

Small RNAs are important regulators of gene expression and are involved in human development and disease. Next generation sequencing (NGS) allows for scalable, genome-wide studies of small RNA; however, current methods are challenged by low sensitivity and high bias, limiting their ability to capture an accurate representation of the cellular small RNA population. Several studies have shown that this bias primarily arises during the ligation of single-strand adapters during library preparation, and that this ligation bias is magnified by 2′-O-methyl modifications (2′OMe) on the 3′ terminal nucleotide. In this study, we developed a novel library preparation process using randomized splint ligation with a cleavable adapter, a design which resolves previous challenges associated with this ligation strategy. We show that a randomized splint ligation based workflow can reduce bias and increase the sensitivity of small RNA sequencing for a wide variety of small RNAs, including microRNA (miRNA) and tRNA fragments as well as 2′OMe modified RNA, including Piwi-interacting RNA and plant miRNA. Finally, we demonstrate that this workflow detects more differentially expressed miRNA between tumorous and matched normal tissues. Overall, this library preparation process allows for highly accurate small RNA sequencing and will enable studies of 2′OMe modified RNA with new levels of detail.

Osteoblastogenesis Alters Small RNA Profiles in EVs Derived from Bone Marrow Stem Cells (BMSCs) and Adipose Stem Cells (ASCs)

Multipotent stem cells (MSCs) are used in various therapeutic applications based on their paracrine secretion activity. Here, we set out to identify and characterize the paracrine factors released during osteoblastogenesis, with a special focus on small non-coding RNAs released in extracellular vesicles (EVs). Bone marrow stem cells (BMSCs) and adipose stem cells (ASCs) from healthy human donors were used as representatives of MSCs. We isolated EVs secreted before and after induction of osteoblastic differentiation and found that the EVs contained a specific subset of microRNAs (miRNAs) and tRNA-derived small RNAs (tsRNA) compared to their parental cells. Osteoblastic differentiation had a larger effect on the small RNA profile of BMSC-EVs relative to ASC-EVs. Our data showed that EVs from different MSC origin exhibited distinct expression profiles of small RNA profiles when undergoing osteoblastogenesis, a factor that should be taken into consideration for stem cell therapy.

Lactate Induces Production of the tRNAHis Half to Promote B-lymphoblastic Cell Proliferation

High plasma lactate is emerging as a critical regulator in development and progression of many human malignancies. Small RNAs derived from cleavage of mature tRNAs have been implicated in many cellular stresses, but the detailed mechanisms that respond to lactic acid (LA; acidic lactate) are not well defined. Here, using an Epstein-Barr virus (EBV)-immortalized B lymphoblastic cell line (LCL) as a model, we report that LA induces cleavage of mature tRNA at the anticodon loop, particularly production of three 5'-tRNA halves (5'-HisGUG, 5'-ValAAC, and 5'-GlyGCC), along with increased expression of RNA polymerase III and angiogenin (ANG). Of these, only the 5'-HisGUG half binds to the chromatin regulator argonaute-2 (AGO2) instead of the AGO1 protein for stability. Notably, the levels of ANG and 5'-HisGUG half expression in peripheral blood mononuclear cells from B cell lymphoma patients are tightly correlated with lactate dehydrogenase (LDH; a lactate indicator) in plasma. Silencing production of the 5'-HisGUG half by small interfering RNA or inhibition of ANG significantly reduces colony formation and growth of LA-induced tumor cells in vitro and in vivo using a murine xenograft model. Overall, our findings identify a novel molecular therapeutic target for the diagnosis and treatment of B cell lymphoma.

Transfer RNA-derived small RNAs: potential applications as novel biomarkers for disease diagnosis and prognosis

Transfer RNA-derived small RNA (tsRNA)s are novel non-coding RNAs, expressed in a variety of tissues and organs. Two subtypes of tsRNAs have been reported: tRNA-derived stress-induced RNA (tiRNA)s and tRNA-derived fragment (tRF)s. tsRNAs have been reported to play essential roles and possess different biological functions in a variety of physiological activities. Recently, tsRNAs have been implicated in a large number of diseases, such as cancers (including breast cancer, ovarian cancer, lung cancer, prostate cancer, colorectal cancer, etc.), neurological disorders, viral infections, metabolic diseases and angiogenesis-related diseases. Although the biological functions of tsRNAs are still poorly understood, correlations between dysregulated tsRNA expression and disease development have been recently reported. Additionally, their capabilities as potential biomarkers for disease diagnosis and prognosis have been revealed in clinical studies. In this review, we summarize the current knowledge of tsRNAs, and discuss their potential clinical applications as biomarkers in different diseases. Although the regulation of tsRNAs is similar to miRNAs in regards to the related physiological and pathological processes, the higher stability and expression levels of tsRNAs place them as ideal biomarkers for the diagnosis and prognosis in cancer and other diseases. Therefore, it is worth to verify the possibility and reliability of these reported tsRNAs as potential biomarkers for clinical applications in disease diagnosis and prognosis.

Oxidative stress enhances the expression of 2',3'-cyclic phosphate-containing RNAs

Eukaryotic cells equip robust systems to respond to stress conditions. In stressed mammalian cells, angiogenin endoribonuclease cleaves anticodon-loops of tRNAs to generate tRNA halves termed tRNA-derived stress-induced RNAs (tiRNAs), which promote stress granule formation and regulate translation. The 5'-tiRNAs (5'-tRNA halves) contain a 2',3'-cyclic phosphate (cP) and thus belong to cP-containing RNAs (cP-RNAs). The cP-RNAs form a hidden layer of the transcriptome because standard RNA-seq cannot amplify and sequence them. In this study, we performed genome-wide analyses of short cP-RNA transcriptome in oxidative stress-exposed human cells. Using cP-RNA-seq that can specifically sequence cP-RNAs, we identified tiRNAs and numerous other cP-RNAs that are mainly derived from rRNAs and mRNAs. Although tiRNAs were produced from a wide variety of tRNA species, abundant species of tiRNAs were derived from a focal-specific subset of tRNAs. Regarding rRNA- and mRNA-derived cP-RNAs, determination of the processing sites of substrate RNAs revealed highly specific RNA cleavage events between pyrimidines and adenosine in generation of those cP-RNAs. Those cP-RNAs were derived from specific loci of substrate RNAs rather than from the overall region, implying that cP-RNAs are produced by regulated biogenesis pathways and not by random degradation events. We experimentally confirmed the identified sequences to be expressed as cP-RNAs in the cells, and their expressions were upregulated upon induction of oxidative stress. These analyses of the cP-RNA transcriptome unravel an abundant class of short ncRNAs that accumulate in cells under oxidative stress.

tRNA 2'-O-methylation by a duo of TRM7/FTSJ1 proteins modulates small RNA silencing in Drosophila

2′-O-Methylation (Nm) represents one of the most common RNA modifications. Nm affects RNA structure and function with crucial roles in various RNA-mediated processes ranging from RNA silencing, translation, self versus non-self recognition to viral defense mechanisms. Here, we identify two Nm methyltransferases (Nm-MTases) in Drosophila melanogaster (CG7009 and CG5220) as functional orthologs of yeast TRM7 and human FTSJ1. Genetic knockout studies together with MALDI-TOF mass spectrometry and RiboMethSeq mapping revealed that CG7009 is responsible for methylating the wobble position in tRNA^Phe, tRNA^Trp and tRNA^Leu, while CG5220 methylates position C32 in the same tRNAs and also targets additional tRNAs. CG7009 or CG5220 mutant animals were viable and fertile but exhibited various phenotypes such as lifespan reduction, small RNA pathways dysfunction and increased sensitivity to RNA virus infections. Our results provide the first detailed characterization of two TRM7 family members in Drosophila and uncover a molecular link between enzymes catalyzing Nm at specific tRNAs and small RNA-induced gene silencing pathways.

AGO unchained: Canonical and non-canonical roles of Argonaute proteins in mammals

Argonaute (AGO) proteins play key roles in animal physiology by binding to small RNAs and regulating the expression of their targets. In mammals, they do so through two distinct pathways: the miRNA pathway represses genes through a multiprotein complex that promotes both decay and translational repression; the siRNA pathway represses transcripts through direct Ago2-mediated cleavage. Here, we review our current knowledge of mechanistic details and physiological requirements of both these pathways and briefly discuss their implications to human disease.

Small RNA Sequencing Reveals Transfer RNA-derived Small RNA Expression Profiles in Retinal Neovascularization

Retinal neovascularization (RNV) is characterized in retinopathy of prematurity (ROP), diabetic retinopathy (DR), and retinal vein occlusion (RVO), which leads to severe vision loss and even blindness. To reveal the altered transfer RNA-derived small RNA (tsRNA)s in RNV, and to investigate the underlying mechanisms of the altered tsRNAs involved in RNV, we carried out a small RNA sequencing to profile tsRNA expressions in the retinas of mice with oxygen-induced retinopathy (OIR) and control mice. A total of 45 tsRNAs were significantly changed (fold change ≥ 1.5 and P < 0.05) in the retinas of OIR mice compared with controls. Validation by quantitative real-time reverse transcription polymerase chain reaction (qRT-PCR) in four selected tsRNAs was consistent with the results of small RNA sequencing. Bioinformatics analyses identified 153 altered target genes of the four validated tsRNAs. These altered target genes were largely enriched in developmental process, cell periphery and protein binding, as well as Th1 and Th2 cell differentiation pathway. Our study suggests tsRNAs play key roles in the pathogenesis of RNV, indicating their therapeutic potential to treat patients with RNV. Moreover, small RNA sequencing is a useful tool to identify changes in tsRNA expression, an important indicator of the progress of retinal diseases.

Cell Cycle-Dependent Regulation and Function of ARGONAUTE1 in Plants

Regulated gene expression is key to the orchestrated progression of the cell cycle. Many genes are expressed at specific points in the cell cycle, including important cell cycle regulators, plus factors involved in signal transduction, hormonal regulation, and metabolic control. We demonstrate that post-embryonic depletion of Arabidopsis (Arabidopsis thaliana) ARGONAUTE1 (AGO1), the main effector of plant microRNAs (miRNAs), impairs cell division in the root meristem. We utilized the highly synchronizable tobacco (Nicotiana tabacum) Bright yellow 2 (BY2) cell suspension to analyze mRNA, small RNAs, and mRNA cleavage products of synchronized BY2 cells at S, G2, M, and G1 phases of the cell cycle. This revealed that in plants, only a few miRNAs show differential accumulation during the cell cycle, and miRNA-target pairs were only identified for a small proportion of the more than 13,000 differentially expressed genes during the cell cycle. However, this unique set of miRNA-target pairs could be key to attenuate the expression of several transcription factors and disease resistance genes. We also demonstrate that AGO1 binds to a set of 19-nucleotide, tRNA-derived fragments during the cell cycle progression.

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.

tRNA-Derived Small Non-Coding RNAs as Novel Epigenetic Molecules Regulating Adipogenesis

tRNA-derived fragments (tRFs), a novel type of non-coding RNA derived from tRNAs, play an important part in governing gene expressions at a post-transcriptional level. To date, the regulatory mechanism of tRFs governing fat deposition and adipogenesis is completely unknown. In this study, high fat diet was employed to induce an obese rat model, and tRFs transcriptome sequencing was conducted to identify differentially expressed tRFs that response to obesity. We found out that tRF^GluTTC, which promoted preadipocyte proliferation by increasing expressions of cell cycle regulatory factors, had the highest fold change in the 296 differentially expressed tRFs. Moreover, tRF^GluTTC also suppressed preadipocyte differentiation by reducing triglyceride content and lipid accumulation, and by decreasing expressions of genes that related to fatty acid synthesis. According to results of luciferase activity analysis, tRF^GluTTC directly targeted Kruppel-like factor (KLF) 9, KLF11, and KLF12, thus significantly suppressing mRNA expressions of these target genes. Moreover, tRF^GluTTC suppressed adipogenesis, accompanying by suppressing expressions of adipogenic transcription factors (aP2, PPARγ, and C/EBPα). In conclusion, these results imply that tRF^GluTTC may act as a novel epigenetic molecule regulating adipogenesis and could provide a new strategy for the intervention treatment of obesity.

tRNA deamination by ADAT requires substrate-specific recognition mechanisms and can be inhibited by tRFs

Adenosine deaminase acting on transfer RNA (ADAT) is an essential eukaryotic enzyme that catalyzes the deamination of adenosine to inosine at the first position of tRNA anticodons. Mammalian ADATs modify eight different tRNAs, having increased their substrate range from a bacterial ancestor that likely deaminated exclusively tRNA^Arg Here we investigate the recognition mechanisms of tRNA^Arg and tRNA^Ala by human ADAT to shed light on the process of substrate expansion that took place during the evolution of the enzyme. We show that tRNA recognition by human ADAT does not depend on conserved identity elements, but on the overall structural features of tRNA. We find that ancestral-like interactions are conserved for tRNA^Arg, while eukaryote-specific substrates use alternative mechanisms. These recognition studies show that human ADAT can be inhibited by tRNA fragments in vitro, including naturally occurring fragments involved in important regulatory pathways.

tRNA Fragments Show Intertwining with mRNAs of Specific Repeat Content and Have Links to Disparities

tRNA-derived fragments (tRF) are a class of potent regulatory RNAs. We mined the datasets from The Cancer Genome Atlas (TCGA) representing 32 cancer types with a deterministic and exhaustive pipeline for tRNA fragments. We found that mitochondrial tRNAs contribute disproportionally more tRFs than nuclear tRNAs. Through integrative analyses, we uncovered a multitude of statistically significant and context-dependent associations between the identified tRFs and mRNAs. In many of the 32 cancer types, these associations involve mRNAs from developmental processes, receptor tyrosine kinase signaling, the proteasome, and metabolic pathways that include glycolysis, oxidative phosphorylation, and ATP synthesis. Even though the pathways are common to multiple cancers, the association of specific mRNAs with tRFs depends on and differs from cancer to cancer. The associations between tRFs and mRNAs extend to genomic properties as well; specifically, tRFs are positively correlated with shorter genes that have a higher density in repeats, such as ALUs, MIRs, and ERVLs. Conversely, tRFs are negatively correlated with longer genes that have a lower repeat density, suggesting a possible dichotomy between cell proliferation and differentiation. Analyses of bladder, lung, and kidney cancer data indicate that the tRF-mRNA wiring can also depend on a patient's sex. Sex-dependent associations involve cyclin-dependent kinases in bladder cancer, the MAPK signaling pathway in lung cancer, and purine metabolism in kidney cancer. Taken together, these findings suggest diverse and wide-ranging roles for tRFs and highlight the extensive interconnections of tRFs with key cellular processes and human genomic architecture. SIGNIFICANCE: Across 32 TCGA cancer contexts, nuclear and mitochondrial tRNA fragments exhibit associations with mRNAs that belong to concrete pathways, encode proteins with particular destinations, have a biased repeat content, and are sex dependent.

The AGO proteins: an overview

Small RNAs govern almost every biological process in eukaryotes associating with the Argonaute (AGO) proteins to form the RNA-induced silencing complex (mRISC). AGO proteins constitute the core of RISCs with different members having variety of protein-binding partners and biochemical properties. This review focuses on the AGO subfamily of the AGOs that are ubiquitously expressed and are associated with small RNAs. The structure, function and role of the AGO proteins in the cell is discussed in detail.

Expression analysis of transfer RNA‑derived fragments in the blood of patients with moyamoya disease: A preliminary study

Moyamoya disease (MMD) is a rare chronic cerebrovascular disease mainly found in individuals of East Asian ethnicity, and its pathogenesis is largely unknown. Transfer RNA-derived fragments (tRFs) are novel biological entities involved in many biological processes; however, whether tRFs contribute towards MMD pathogenesis remains unexplored. In the present study, deep sequencing technology was used to identify alterations in tRF expression profiles between patients with MMD and healthy controls. The sequencing findings were validated using reverse transcription-quantitative polymerase chain reaction (RT-qPCR). Subsequently, the putative target genes of tRFs were predicted using miRNA target prediction algorithms. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed to further evaluate potential functions of tRFs. The sequencing results demonstrated that 38 tRFs were differentially expressed between patients and controls, of which 22 were upregulated and 16 were downregulated. RT-qPCR analysis confirmed the validity of the sequencing results. GO and KEGG pathway enrichment analyses indicated that 15 pathways were associated with the selected tRFs. These pathways were mainly involved in angiogenesis and metabolism, both of which are physiopathological fundamentals of MMD. The results provided a novel insight into the mechanisms underlying MMD pathogenesis, and demonstrated that tRFs may serve as potential therapeutic targets for the future treatment of MMD.

Generation of 2',3'-Cyclic Phosphate-Containing RNAs as a Hidden Layer of the Transcriptome

Cellular RNA molecules contain phosphate or hydroxyl ends. A 2′,3′-cyclic phosphate (cP) is one of the 3′-terminal forms of RNAs mainly generated from RNA cleavage by ribonucleases. Although transcriptome profiling using RNA-seq has become a ubiquitous tool in biological and medical research, cP-containing RNAs (cP-RNAs) form a hidden transcriptome layer, which is infrequently recognized and characterized, because standard RNA-seq is unable to capture them. Despite cP-RNAs’ invisibility in RNA-seq data, increasing evidence indicates that they are not accumulated simply as non-functional degradation products; rather, they have physiological roles in various biological processes, designating them as noteworthy functional molecules. This review summarizes our current knowledge of cP-RNA biogenesis pathways and their catalytic enzymatic activities, discusses how the cP-RNA generation affects biological processes, and explores future directions to further investigate cP-RNA biology.

Transfer RNA-derived fragments and tRNA halves: biogenesis, biological functions and their roles in diseases

The number of studies on non-coding RNAs has increased substantially in recent years owing to their importance in gene regulation. However, the biological functions of small RNAs from abundant species of housekeeping non-coding RNAs (rRNA, tRNA, etc.) remain a highly studied topic. tRNA-derived small RNAs (tsRNAs) refer to the specific cleavage of tRNAs by specific nucleases [e.g., Dicer and angiogenin (ANG)] in particular cells or tissues or under certain conditions such as stress and hypoxia. tsRNAs are a type of non-coding small RNA that are widely found in the prokaryotic and eukaryotic transcriptomes and are generated from mature tRNAs or precursor tRNAs at different sites. There are two main types of tsRNAs, tRNA-derived fragments (tRFs) and tRNA halves. tRFs are 14-30 nucleotides (nt) long and mainly consist of three subclasses: tRF-5, tRF-3, and tRF-1. tRNA halves, which are 31-40 nt long, are generated by specific cleavage in the anticodon loops of mature tRNAs. There are two types of tRNA halves, 5'-tRNA halves and 3'-tRNA halves. tsRNAs have multiple biological functions including acting as signaling molecules in stress responses and as regulators of gene expression. Additionally, they have been considered to be involved in RNA processing, cell proliferation, translation suppression, the modulation of DNA damage response, and neurodegeneration. More importantly, they are closely related to the occurrence of many human diseases such as tumors, infectious diseases, metabolic diseases, and neurological diseases. Moreover, tsRNAs have the potential to become new biomarkers for disease diagnosis. Continuous investigations will help us to understand their generation and regulatory mechanisms as well as the possible roles of tRFs and tRNA halves.

Human tRNA-Derived Small RNAs Modulate Host-Oral Microbial Interactions

Coevolution of the human host and its associated microbiota has led to sophisticated interactions to maintain a delicate homeostasis. Emerging evidence suggests that in addition to small molecules, peptides, and proteins, small regulatory noncoding RNAs (sRNAs) might play an important role in cross-domain interactions. In this study, we revealed the presence of diverse host transfer RNA-derived small RNAs (tsRNAs) among human salivary sRNAs. We selected 2 tsRNAs (tsRNA-000794 and tsRNA-020498) for further study based on their high sequence similarity to specific tRNAs from a group of Gram-negative oral bacteria, including Fusobacterium nucleatum, a key oral commensal and opportunistic pathogen. We showed that the presence of F. nucleatum triggers exosome-mediated release of tsRNA-000794 and tsRNA-020498 by human normal oral keratinocyte cells. Furthermore, both tsRNA candidates exerted a growth inhibition effect on F. nucleatum, likely through interference with bacterial protein biosynthesis, but did not affect the growth of Streptococcus mitis, a health-associated oral Gram-positive bacterium whose genome does not carry sequences bearing high similarity to either tsRNA. Our data provide the first line of evidence for the modulatory role of host-derived tsRNAs in the microbial-host interaction.

MINTbase v2.0: a comprehensive database for tRNA-derived fragments that includes nuclear and mitochondrial fragments from all The Cancer Genome Atlas projects

MINTbase is a repository that comprises nuclear and mitochondrial tRNA-derived fragments ('tRFs') found in multiple human tissues. The original version of MINTbase comprised tRFs obtained from 768 transcriptomic datasets. We used our deterministic and exhaustive tRF mining pipeline to process all of The Cancer Genome Atlas datasets (TCGA). We identified 23 413 tRFs with abundance of ≥ 1.0 reads-per-million (RPM). To facilitate further studies of tRFs by the community, we just released version 2.0 of MINTbase that contains information about 26 531 distinct human tRFs from 11 719 human datasets as of October 2017. Key new elements include: the ability to filter tRFs on-the-fly by minimum abundance thresholding; the ability to filter tRFs by tissue keywords; easy access to information about a tRF's maximum abundance and the datasets that contain it; the ability to generate relative abundance plots for tRFs across cancer types and convert them into embeddable figures; MODOMICS information about modifications of the parental tRNA, etc. Version 2.0 of MINTbase contains 15x more datasets and nearly 4x more distinct tRFs than the original version, yet continues to offer fast, interactive access to its contents. Version 2.0 is available freely at http://cm.jefferson.edu/MINTbase/.

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.

Race Disparities in the Contribution of miRNA Isoforms and tRNA-Derived Fragments to Triple-Negative Breast Cancer

Triple-negative breast cancer (TNBC) is a breast cancer subtype characterized by marked differences between White and Black/African-American women. We performed a systems-level analysis on datasets from The Cancer Genome Atlas to elucidate how the expression patterns of mRNAs are shaped by regulatory noncoding RNAs (ncRNA). Specifically, we studied isomiRs, that is, isoforms of miRNAs, and tRNA-derived fragments (tRF). In normal breast tissue, we observed a marked cohesiveness in both the ncRNA and mRNA layers and the associations between them. This cohesiveness was widely disrupted in TNBC. Many mRNAs become either differentially expressed or differentially wired between normal breast and TNBC in tandem with isomiR or tRF dysregulation. The affected pathways included energy metabolism, cell signaling, and immune responses. Within TNBC, the wiring of the affected pathways with isomiRs and tRFs differed in each race. Multiple isomiRs and tRFs arising from specific miRNA loci (e.g., miR-200c, miR-21, the miR-17/92 cluster, the miR-183/96/182 cluster) and from specific tRNA loci (e.g., the nuclear tRNA^Gly and tRNA^Leu, the mitochondrial tRNA^Val and tRNA^Pro) were strongly associated with the observed race disparities in TNBC. We highlight the race-specific aspects of transcriptome wiring by discussing in detail the metastasis-related MAPK and the Wnt/β-catenin signaling pathways, two of the many key pathways that were found differentially wired. In conclusion, by employing a data- and knowledge-driven approach, we comprehensively analyzed the normal and cancer transcriptomes to uncover novel key contributors to the race-based disparities of TNBC.Significance: This big data-driven study comparing normal and cancer transcriptomes uncovers RNA expression differences between Caucasian and African-American patients with triple-negative breast cancer that might help explain disparities in incidence and aggressive character. Cancer Res; 78(5); 1140-54. ©2017 AACR.

The biogenesis pathway of tRNA-derived piRNAs in Bombyx germ cells

Transfer RNAs (tRNAs) function in translational machinery and further serves as a source of short non-coding RNAs (ncRNAs). tRNA-derived ncRNAs show differential expression profiles and play roles in many biological processes beyond translation. Molecular mechanisms that shape and regulate their expression profiles are largely unknown. Here, we report the mechanism of biogenesis for tRNA-derived Piwi-interacting RNAs (td-piRNAs) expressed in Bombyx BmN4 cells. In the cells, two cytoplasmic tRNA species, tRNA^AspGUC and tRNA^HisGUG, served as major sources for td-piRNAs, which were derived from the 5′-part of the respective tRNAs. cP-RNA-seq identified the two tRNAs as major substrates for the 5′-tRNA halves as well, suggesting a previously uncharacterized link between 5′-tRNA halves and td-piRNAs. An increase in levels of the 5′-tRNA halves, induced by BmNSun2 knockdown, enhanced the td-piRNA expression levels without quantitative change in mature tRNAs, indicating that 5′-tRNA halves, not mature tRNAs, are the direct precursors for td-piRNAs. For the generation of tRNA^HisGUG-derived piRNAs, BmThg1l-mediated nucleotide addition to −1 position of tRNA^HisGUG was required, revealing an important function of BmThg1l in piRNA biogenesis. Our study advances the understanding of biogenesis mechanisms and the genesis of specific expression profiles for tRNA-derived ncRNAs.

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.

5'-Terminal nucleotide variations in human cytoplasmic tRNAHisGUG and its 5'-halves

Transfer RNAs (tRNAs) are fundamental adapter components of translational machinery. tRNAs can further serve as a source of tRNA-derived noncoding RNAs that play important roles in various biological processes beyond translation. Among all species of tRNAs, tRNA^HisGUG has been known to uniquely contain an additional guanosine residue at the -1 position (G_-1) of its 5'-end. To analyze this -1 nucleotide in detail, we developed a TaqMan qRT-PCR method that can distinctively quantify human mature cytoplasmic tRNA^HisGUG containing G_-1, U_-1, A_-1, or C_-1 or lacking the -1 nucleotide (starting from G₁). Application of this method to the mature tRNA fraction of BT-474 breast cancer cells revealed the presence of tRNA^HisGUG containing U_-1 as well as the one containing G_-1 Moreover, tRNA lacking the -1 nucleotide was also detected, thus indicating the heterogeneous expression of 5'-tRNA^HisGUG variants. A sequence library of sex hormone-induced 5'-tRNA halves (5'-SHOT-RNAs), identified via cP-RNA-seq of a BT-474 small RNA fraction, also demonstrated the expression of 5'-tRNA^HisGUG halves containing G_-1, U_-1, or G₁ as 5'-terminal nucleotides. Although the detected 5'-nucleotide species were identical, the relative abundances differed widely between mature tRNA and 5'-half from the same BT-474 cells. The majority of mature tRNAs contained the -1 nucleotide, whereas the majority of 5'-halves lacked this nucleotide, which was biochemically confirmed using a primer extension assay. These results reveal the novel identities of tRNA^HisGUG molecules and provide insights into tRNA^HisGUG maturation and the regulation of tRNA half production.

Nucleoside modifications in the regulation of gene expression: focus on tRNA

Both, DNA and RNA nucleoside modifications contribute to the complex multi-level regulation of gene expression. Modified bases in tRNAs modulate protein translation rates in a highly dynamic manner. Synonymous codons, which differ by the third nucleoside in the triplet but code for the same amino acid, may be utilized at different rates according to codon-anticodon affinity. Nucleoside modifications in the tRNA anticodon loop can favor the interaction with selected codons by stabilizing specific base pairs. Similarly, weakening of base pairing can discriminate against binding to near-cognate codons. mRNAs enriched in favored codons are translated in higher rates constituting a fine-tuning mechanism for protein synthesis. This so-called codon bias establishes a basic protein level, but sometimes it is necessary to further adjust the production rate of a particular protein to actual requirements, brought by, e.g., stages in circadian rhythms, cell cycle progression or exposure to stress. Such an adjustment is realized by the dynamic change of tRNA modifications resulting in the preferential translation of mRNAs coding for example for stress proteins to facilitate cell survival. Furthermore, tRNAs contribute in an entirely different way to another, less specific stress response consisting in modification-dependent tRNA cleavage that contributes to the general down-regulation of protein synthesis. In this review, we summarize control functions of nucleoside modifications in gene regulation with a focus on recent findings on protein synthesis control by tRNA base modifications.

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.

Role of Non-Coding RNAs in the Transgenerational Epigenetic Transmission of the Effects of Reprotoxicants

Non-coding RNAs (ncRNAs) are regulatory elements of gene expression and chromatin structure. Both long and small ncRNAs can also act as inductors and targets of epigenetic programs. Epigenetic patterns can be transmitted from one cell to the daughter cell, but, importantly, also through generations. Diversity of ncRNAs is emerging with new and surprising roles. Functional interactions among ncRNAs and between specific ncRNAs and structural elements of the chromatin are drawing a complex landscape. In this scenario, epigenetic changes induced by environmental stressors, including reprotoxicants, can explain some transgenerationally-transmitted phenotypes in non-Mendelian ways. In this review, we analyze mechanisms of action of reprotoxicants upon different types of ncRNAs and epigenetic modifications causing transgenerationally transmitted characters through germ cells but affecting germ cells and reproductive systems. A functional model of epigenetic mechanisms of transgenerational transmission ncRNAs-mediated is also proposed.