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

Transfer RNA-derived fragment tRF-23-Q99P9P9NDD promotes progression of gastric cancer by targeting ACADSB

Gastric cancer (GC) is one of the most common gastrointestinal tumors. As a newly discovered type of non-coding RNAs, transfer RNA (tRNA)‍-derived small RNAs (tsRNAs) play a dual biological role in cancer. Our previous studies have demonstrated the potential of tRF-23-Q99P9P9NDD as a diagnostic and prognostic biomarker for GC. In this work, we confirmed for the first time that tRF-23-Q99P9P9NDD can promote the proliferation, migration, and invasion of GC cells in vitro. The dual luciferase reporter gene assay confirmed that tRF-23-Q99P9P9NDD could bind to the 3' untranslated region (UTR) site of acyl-coenzyme A dehydrogenase short/branched chain (ACADSB). In addition, ACADSB could rescue the effect of tRF-23-Q99P9P9NDD on GC cells. Next, we used Gene Ontology (GO), the Kyoto Encyclopedia of Genes and Genomes (KEGG), and Gene Set Enrichment Analysis (GSEA) to find that downregulated ACADSB in GC may promote lipid accumulation by inhibiting fatty acid catabolism and ferroptosis. Finally, we verified the correlation between ACADSB and 12 ferroptosis genes at the transcriptional level, as well as the changes in reactive oxygen species (ROS) levels by flow cytometry. In summary, this study proposes that tRF-23-Q99P9P9NDD may affect GC lipid metabolism and ferroptosis by targeting ACADSB, thereby promoting GC progression. It provides a theoretical basis for the diagnostic and prognostic monitoring value of GC and opens up new possibilities for treatment.

Altered expression of transfer RNAs and their possible roles in brain white matter injury

Transfer RNAs (tRNAs) can regulate cell behavior and are associated with neurological disorders. Here, we aimed to investigate the expression levels of tRNAs in oligodendrocyte precursor cells (OPCs) and their possible roles in the regulation of brain white matter injury (WMI). Newborn Sprague-Dawley rats (postnatal day 5) were used to establish a model that mimicked neonatal brain WMI. RNA-array analysis was performed to examine the expression of tRNAs in OPCs. psRNAtarget software was used to predict target mRNAs of significantly altered tRNAs. Gene ontology (GO) and KEGG were used to analyze the pathways for target mRNAs. Eighty-nine tRNAs were changed after WMI (fold change absolute ≥1.5, P  < 0.01), with 31 downregulated and 58 upregulated. Among them, three significantly changed tRNAs were identified, with two being significantly increased (chr10.trna1314-ProTGG and chr2.trna2771-ProAGG) and one significantly decreased (chr10.trna11264-GlyTCC). Further, target mRNA prediction and GO/KEGG pathway analysis indicated that the target mRNAs of these tRNAs are mainly involved in G-protein coupled receptor signaling pathways and beta-alanine metabolism, which are both related to myelin formation. In summary, the expression of tRNAs in OPCs was significantly altered after brain WMI, suggesting that tRNAs may play important roles in regulating WMI. This improves the knowledge about WMI pathophysiology and may provide novel treatment targets for WMI.

Cell-type-specific expression of tRNAs in the brain regulates cellular homeostasis

Defects in tRNA biogenesis are associated with multiple neurological disorders, yet our understanding of these diseases has been hampered by an inability to determine tRNA expression in individual cell types within a complex tissue. Here, we developed a mouse model in which RNA polymerase III is conditionally epitope tagged in a Cre-dependent manner, allowing us to accurately profile tRNA expression in any cell type in vivo. We investigated tRNA expression in diverse nervous system cell types, revealing dramatic heterogeneity in the expression of tRNA genes between populations. We found that while maintenance of levels of tRNA isoacceptor families is critical for cellular homeostasis, neurons are differentially vulnerable to insults to distinct tRNA isoacceptor families. Cell-type-specific translatome analysis suggests that the balance between tRNA availability and codon demand may underlie such differential resilience. Our work provides a platform for investigating the complexities of mRNA translation and tRNA biology in the brain.

Novel insights into transfer RNA-derived small RNA (tsRNA) in cardio-metabolic diseases

Cardio-metabolic diseases, including a cluster of metabolic disorders and their secondary affections on cardiovascular physiology, are gradually brought to the forefront by researchers due to their high prevalence and mortality, as well as an unidentified pathogenesis. tRNA-derived small RNAs (tsRNAs), cleaved by several specific enzymes and once considered as some "metabolic junks" in the past, have been proved to possess numerous functions in human bodies. More interestingly, such a potential also seems to influence the progression of cardio-metabolic diseases to some extent. In this review, the biogenesis, classification and mechanisms of tsRNAs will be discussed based on some latest studies, and their relations with several cardio-metabolic diseases will be highlighted in sequence. Lastly, some future prospects, such as their clinical applications as biomarkers and therapeutic targets will also be mentioned, in order to provide researchers with a comprehensive understanding of the research status of tsRNAs as well as its association with cardio-metabolic diseases, thus presenting as a beacon to indicate directions for the next stage of study.

Multiple regulatory roles of the transfer RNA-derived small RNAs in cancers

With the development of sequencing technology, transfer RNA (tRNA)-derived small RNAs (tsRNAs) have received extensive attention as a new type of small noncoding RNAs. Based on the differences in the cleavage sites of nucleases on tRNAs, tsRNAs can be divided into two categories, tRNA halves (tiRNAs) and tRNA-derived fragments (tRFs), each with specific subcellular localizations. Additionally, the biogenesis of tsRNAs is tissue-specific and can be regulated by tRNA modifications. In this review, we first elaborated on the classification and biogenesis of tsRNAs. After summarizing the latest mechanisms of tsRNAs, including transcriptional gene silencing, post-transcriptional gene silencing, nascent RNA silencing, translation regulation, rRNA regulation, and reverse transcription regulation, we explored the representative biological functions of tsRNAs in tumors. Furthermore, this review summarized the clinical value of tsRNAs in cancers, thus providing theoretical support for their potential as novel biomarkers and therapeutic targets.

RNA modifications in physiology and disease: towards clinical applications

The ability of chemical modifications of single nucleotides to alter the electrostatic charge, hydrophobic surface and base pairing of RNA molecules is exploited for the clinical use of stable artificial RNAs such as mRNA vaccines and synthetic small RNA molecules - to increase or decrease the expression of therapeutic proteins. Furthermore, naturally occurring biochemical modifications of nucleotides regulate RNA metabolism and function to modulate crucial cellular processes. Studies showing the mechanisms by which RNA modifications regulate basic cell functions in higher organisms have led to greater understanding of how aberrant RNA modification profiles can cause disease in humans. Together, these basic science discoveries have unravelled the molecular and cellular functions of RNA modifications, have provided new prospects for therapeutic manipulation and have led to a range of innovative clinical approaches.

Selective gene expression maintains human tRNA anticodon pools during differentiation

Transfer RNAs are essential for translating genetic information into proteins. The human genome contains hundreds of predicted tRNA genes, many in multiple copies. How their expression is regulated to control tRNA repertoires is unknown. Here we combined quantitative tRNA profiling and chromatin immunoprecipitation with sequencing to measure tRNA expression following the differentiation of human induced pluripotent stem cells into neuronal and cardiac cells. We find that tRNA transcript levels vary substantially, whereas tRNA anticodon pools, which govern decoding rates, are more stable among cell types. Mechanistically, RNA polymerase III transcribes a wide range of tRNA genes in human induced pluripotent stem cells but on differentiation becomes constrained to a subset we define as housekeeping tRNAs. This shift is mediated by decreased mTORC1 signalling, which activates the RNA polymerase III repressor MAF1. Our data explain how tRNA anticodon pools are buffered to maintain decoding speed across cell types and reveal that mTORC1 drives selective tRNA expression during differentiation.

Implications of tRNA abundance on translation elongation across bovine tissues

Introduction: Translation is a crucial stage of gene expression. It may also act as an additional layer of regulation that plays an important role in gene expression and function. Highly expressed genes are believed to be codon-biased to support increased protein production, in which quickly translated codons correspond to highly abundant tRNAs. Synonymous SNPs, considered to be silent due to the degeneracy of the genetic code, may shift protein abundance and function through alterations in translational efficiency and suboptimal pairing to lowly abundant tRNAs. Methods: Here, we applied Quantitative Mature tRNA sequencing (QuantM-tRNAseq) and ribosome profiling across bovine tissues in order to investigate the relationship between tRNA expression and slowed translation. Results: Moreover, we have identified genes modulated at transcriptional and/or translational levels underlying tissue-specific biological processes. We have also successfully defined pausing sites that depict the regulatory information encoded within the open reading frame of transcripts, which could be related to translation rate and facilitate proper protein folding. This work offers an atlas of distinctive pausing sites across three bovine tissues, which provides an opportunity to predict codon optimality and understand tissue-specific mechanisms of regulating protein synthesis.

Inhibition of lncRNA RPPH1 activity decreases tumor proliferation and metastasis through down-regulation of inflammation-related oncogenes

OBJECTIVE

Ribonuclease P RNA component H1 (RPPH1) is a long non-coding RNA (lncRNA) associated with cancer progression. Higher RPPH1 expression in breast and cervical cancer samples than that in normal tissues were observed through the lncRNASNP2 database; therefore, silencing RPPH1 expression might be a potential strategy for cancer treatments, even though RPPH1 is also an RNA subunit of ribonuclease P involved in processing transfer RNA (tRNA) precursors and the effect of RPPH1 knockdown is not yet fully understood.

METHODS

Differentially expressed genes (DEGs) were identified through RNA sequencing in each shRNA-transfected RPPH1 knockdown MDA-MB-231, RPPH1 knockdown HeLa cell, and respective control cells, then the gene ontology enrichment analysis was performed by IPA and MetaCore database according to these DEGs, with further in vitro experiments validating the effect of RPPH1 silencing in MDA-MB-231 and HeLa cells.

RESULTS

Hundreds of down-regulated DEGs were identified in RPPH1 knockdown MDA-MB-231 and HeLa cells while bioinformatics analysis revealed that these genes were involved in pathways related to immune response and cancerogenesis. Compared to mock- and vector-transfected cells, the production of mature tRNAs, cell proliferation and migration capacity were inhibited in RPPH1-silenced HeLa and MDA-MB-231 cells. Additionally, RPPH1 knockdown promoted G1 cell cycle arrest mainly through the down-regulation of cyclin D1, although glycolytic pathways were only affected in RPPH1 knockdown HeLa cells but not MDA-MB-231 cells.

CONCLUSION

This study demonstrated that knockdown RPPH1 affected tRNA production, cell proliferation and metabolism. Our findings might provide insight into the role of RPPH1 in tumor development.

Pan-cancer tRNA-derived fragment CAT1 coordinates RBPMS to stabilize NOTCH2 mRNA to promote tumorigenesis

Transfer RNA-derived fragments (tRFs) are a class of small non-coding regulatory RNAs that are involved in the pathophysiology of many diseases. However, the role of tRFs in cancer progression remains largely elusive. Here, we demonstrate that a pan-cancer 3'-tRF, CAT1 (cancer associated tRF 1), is ubiquitously upregulated in tumors and associated with poor prognosis of a variety of cancers, including lung cancer. The upregulated CAT1 in cancer cells binds to RNA-binding protein with multiple splicing (RBPMS) and displaces NOTCH2 association from RBPMS, thereby inhibiting the subsequent CCR4-NOT deadenylation-complex-mediated NOTCH2 mRNA decay. The CAT1-enhanced NOTCH2 expression promotes lung cancer cell proliferation and metastasis in vitro and in vivo. In addition, plasma CAT1 levels are substantially increased in patients with lung cancer compared to non-cancer control subjects. Our findings reveal an intrinsic connection between cancer-specific upregulation of CAT1 and cancer progression, show the regulation of NOTCH signaling in cancer by a 3'-tRF, and highlight its great clinical potential.

3'-tRF-CysGCA overexpression in HEK-293 cells alters the global expression profile and modulates cellular processes and pathways

tRNA fragments (tRFs) are small non-coding RNAs generated through specific cleavage of tRNAs and involved in various biological processes. Among the different types of tRFs, the 3'-tRFs have attracted scientific interest due to their regulatory role in gene expression. In this study, we investigated the role of 3'-tRF-CysGCA, a tRF deriving from cleavage in the T-loop of tRNACysGCA, in the regulation of gene expression in HEK-293 cells. Previous studies have shown that 3'-tRF-CysGCA is incorporated into the RISC complex and interacts with Argonaute proteins, suggesting its involvement in the regulation of gene expression. However, the general role and effect of the deregulation of 3'-tRF-CysGCA levels in human cells have not been investigated so far. To fill this gap, we stably overexpressed 3'-tRF-CysGCA in HEK-293 cells and performed transcriptomic and proteomic analyses. Moreover, we validated the interaction of this tRF with putative targets, the levels of which were found to be affected by 3'-tRF-CysGCA overexpression. Lastly, we investigated the implication of 3'-tRF-CysGCA in various pathways using extensive bioinformatics analysis. Our results indicate that 3'-tRF-CysGCA overexpression led to changes in the global gene expression profile of HEK-293 cells and that multiple cellular pathways were affected by the deregulation of the levels of this tRF. Additionally, we demonstrated that 3'-tRF-CysGCA directly interacts with thymopoietin (TMPO) transcript variant 1 (also known as LAP2α), leading to modulation of its levels. In conclusion, our findings suggest that 3'-tRF-CysGCA plays a significant role in gene expression regulation and highlight the importance of this tRF in cellular processes.

Cell type-specific translational regulation by human DUS enzymes

Dihydrouridine is an abundant and conserved modified nucleoside present on tRNA, but characterization and functional studies of modification sites and associated DUS writer enzymes in mammals is lacking. Here we use a chemical probing strategy, RNABPP-PS, to identify 5-chlorouridine as an activity-based probe for human DUS enzymes. We map D modifications using RNA-protein crosslinking and chemical transformation and mutational profiling to reveal D modification sites on human tRNAs. Further, we knock out individual DUS genes in two human cell lines to investigate regulation of tRNA expression levels and codon-specific translation. We show that whereas D modifications are present across most tRNA species, loss of D only perturbs the translational function of a subset of tRNAs in a cell type-specific manner. Our work provides powerful chemical strategies for investigating D and DUS enzymes in diverse biological systems and provides insight into the role of a ubiquitous tRNA modification in translational regulation.

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.

Emerging functional principles of tRNA-derived small RNAs and other regulatory small RNAs

Recent advancements in small RNA sequencing have unveiled a previously hidden world of regulatory small noncoding RNAs (sncRNAs) that extend beyond the well-studied small interfering RNAs, microRNAs, and piwi-interacting RNAs. This exploration, starting with tRNA-derived small RNAs, has led to the discovery of a diverse universe of sncRNAs derived from various longer structured RNAs such as rRNAs, small nucleolar RNAs, small nuclear RNAs, Y RNAs, and vault RNAs, with exciting uncharted functional possibilities. In this perspective, we discuss the emerging functional principles of sncRNAs beyond the well-known RNAi-like mechanisms, focusing on those that operate independent of linear sequence complementarity but rather function in an aptamer-like fashion. Aptamers use 3D structure for specific interactions with ligands and are modulated by RNA modifications and subcellular environments. Given that aptamer-like sncRNA functions are widespread and present in species lacking RNAi, they may represent an ancient functional principle that predates RNAi. We propose a rethinking of the origin of RNAi and its relationship with these aptamer-like functions in sncRNAs and how these complementary mechanisms shape biological processes. Lastly, the aptamer-like function of sncRNAs highlights the need for caution in using small RNA mimics in research and therapeutics, as their specificity is not restricted solely to linear sequence.

The functions of a 5' tRNA-Ala-derived fragment in gene expression

Transfer RNA (tRNA) can produce smaller RNA fragments called tRNA-derived fragments (tRFs). tRFs play critical roles in multiple cellular programs, although the functional mechanisms of tRFs remain largely unknown in plants. In this study, we examined the phenotype associated with 5' tRF-Ala (tRF-Ala, produced from tRNA-Ala) overexpression and knockdown lines (tDR-Ala-OE and tDR-Ala-kd, respectively) and the mechanisms by which tRF-Ala affects mRNA levels in Arabidopsis (Arabidopsis thaliana). We investigated the candidate proteins associated with tRF-Ala by quantitative proteomics and confirmed the direct interaction between tRF-Ala and the splicing factor SERINE-ARGININE RICH PROTEIN 34 (SR34). A transcriptome sequencing analysis showed that 318 genes among all the genes (786) with substantial alternative splicing (AS) variance in tDR-Ala-OE lines are targets of SR34. tRF-Ala diminished the binding affinity between SR34 and its targets by direct competition for interaction with SR34. These findings reveal the critical roles of tRF-Ala in regulating mRNA levels and splicing.

Cysteine tRNA acts as a stop codon readthrough-inducing tRNA in the human HEK293T cell line

Under certain circumstances, any of the three termination codons can be read through by a near-cognate tRNA; i.e., a tRNA whose two out of three anticodon nucleotides base pair with those of the stop codon. Unless programed to synthetize C-terminally extended protein variants with expanded physiological roles, readthrough represents an undesirable translational error. On the other side of a coin, a significant number of human genetic diseases is associated with the introduction of nonsense mutations (premature termination codons [PTCs]) into coding sequences, where stopping is not desirable. Here, the tRNA's ability to induce readthrough opens up the intriguing possibility of mitigating the deleterious effects of PTCs on human health. In yeast, the UGA and UAR stop codons were described to be read through by four readthrough-inducing rti-tRNAs-tRNATrp and tRNACys, and tRNATyr and tRNAGln, respectively. The readthrough-inducing potential of tRNATrp and tRNATyr was also observed in human cell lines. Here, we investigated the readthrough-inducing potential of human tRNACys in the HEK293T cell line. The tRNACys family consists of two isoacceptors, one with ACA and the other with GCA anticodons. We selected nine representative tRNACys isodecoders (differing in primary sequence and expression level) and tested them using dual luciferase reporter assays. We found that at least two tRNACys can significantly elevate UGA readthrough when overexpressed. This indicates a mechanistically conserved nature of rti-tRNAs between yeast and human, supporting the idea that they could be used in the PTC-associated RNA therapies.

A ligation-independent sequencing method reveals tRNA-derived RNAs with blocked 3' termini

Despite the numerous sequencing methods available, the vast diversity in size and chemical modifications of RNA molecules makes the capture of the full spectrum of cellular RNAs a difficult task. By combining quasi-random hexamer priming with a custom template switching strategy, we developed a method to construct sequencing libraries from RNA molecules of any length and with any type of 3' terminal modification, allowing the sequencing and analysis of virtually all RNA species. Ligation-independent detection of all types of RNA (LIDAR) is a simple, effective tool to comprehensively characterize changes in small non-coding RNAs and mRNAs simultaneously, with performance comparable to separate dedicated methods. With LIDAR, we comprehensively characterized the coding and non-coding transcriptome of mouse embryonic stem cells, neural progenitor cells, and sperm. LIDAR detected a much larger variety of tRNA-derived RNAs (tDRs) compared to traditional ligation-dependent sequencing methods, and uncovered the presence of tDRs with blocked 3' ends that had previously escaped detection. Our findings highlight the potential of LIDAR to systematically detect all RNAs in a sample and uncover new RNA species with potential regulatory functions.

An old friend with a new face: tRNA-derived small RNAs with big regulatory potential in cancer biology

Transfer RNAs (tRNAs) are small non-coding RNAs (sncRNAs) essential for protein translation. Emerging evidence suggests that tRNAs can also be processed into smaller fragments, tRNA-derived small RNAs (tsRNAs), a novel class of sncRNAs with powerful applications and high biological relevance to cancer. tsRNAs biogenesis is heterogeneous and involves different ribonucleases, such as Angiogenin and Dicer. For many years, tsRNAs were thought to be just degradation products. However, accumulating evidence shows their roles in gene expression: either directly via destabilising the mRNA or the ribosomal machinery, or indirectly via regulating the expression of ribosomal components. Furthermore, tsRNAs participate in various biological processes linked to cancer, including apoptosis, cell cycle, immune response, and retroviral insertion into the human genome. It is emerging that tsRNAs have significant therapeutic potential. Endogenous tsRNAs can be used as cancer biomarkers, while synthetic tsRNAs and antisense oligonucleotides can be employed to regulate gene expression. In this review, we are recapitulating the regulatory roles of tsRNAs, with a focus on cancer biology.

The role and mechanism of action of tRNA-derived fragments in the diagnosis and treatment of malignant tumors

Cancer is a leading cause of morbidity and death worldwide. While various factors are established as causing malignant tumors, the mechanisms underlying cancer development remain poorly understood. Early diagnosis and the development of effective treatments for cancer are important research topics. Transfer RNA (tRNA), the most abundant class of RNA molecules in the human transcriptome, participates in both protein synthesis and cellular metabolic processes. tRNA-derived fragments (tRFs) are produced by specific cleavage of pre-tRNA and mature tRNA molecules, which are highly conserved and occur widely in various organisms. tRFs were initially thought to be random products with no physiological function, but have been redefined as novel functional small non-coding RNA molecules that help to regulate RNA stability, modulate translation, and influence target gene expression, as well as other biological processes. There is increasing evidence supporting roles for tRFs in tumorigenesis and cancer development, including the regulation of tumor cell proliferation, invasion, migration, and drug resistance. Understanding the regulatory mechanisms by which tRFs impact these processes has potential to inform malignant tumor diagnosis and treatment. Further, tRFs are expected to become new biological markers for early diagnosis and prognosis prediction in patients with tumors, as well as a targets for precision cancer therapies. Video abstract.

Perseverance of protein homeostasis despite mistranslation of glycine codons with alanine

By linking amino acids to their codon assignments, transfer RNAs (tRNAs) are essential for protein synthesis and translation fidelity. Some human tRNA variants cause amino acid mis-incorporation at a codon or set of codons. We recently found that a naturally occurring tRNASer variant decodes phenylalanine codons with serine and inhibits protein synthesis. Here, we hypothesized that human tRNA variants that misread glycine (Gly) codons with alanine (Ala) will also disrupt protein homeostasis. The A3G mutation occurs naturally in tRNAGly variants (tRNAGlyCCC, tRNAGlyGCC) and creates an alanyl-tRNA synthetase (AlaRS) identity element (G3 : U70). Because AlaRS does not recognize the anticodon, the human tRNAAlaAGC G35C (tRNAAlaACC) variant may function similarly to mis-incorporate Ala at Gly codons. The tRNAGly and tRNAAla variants had no effect on protein synthesis in mammalian cells under normal growth conditions; however, tRNAGlyGCC A3G depressed protein synthesis in the context of proteasome inhibition. Mass spectrometry confirmed Ala mistranslation at multiple Gly codons caused by the tRNAGlyGCC A3G and tRNAAlaAGC G35C mutants, and in some cases, we observed multiple mistranslation events in the same peptide. The data reveal mistranslation of Ala at Gly codons and defects in protein homeostasis generated by natural human tRNA variants that are tolerated under normal conditions. This article is part of the theme issue 'Reactivity and mechanism in chemical and synthetic biology'.

Piwil2 (Mili) sustains neurogenesis and prevents cellular senescence in the postnatal hippocampus

Adult neural progenitor cells (aNPCs) ensure lifelong neurogenesis in the mammalian hippocampus. Proper regulation of aNPC fate has thus important implications for brain plasticity and healthy aging. Piwi proteins and the small noncoding RNAs interacting with them (piRNAs) have been proposed to control memory and anxiety, but the mechanism remains elusive. Here, we show that Piwil2 (Mili) is essential for proper neurogenesis in the postnatal mouse hippocampus. RNA sequencing of aNPCs and their differentiated progeny reveal that Mili and piRNAs are dynamically expressed in neurogenesis. Depletion of Mili and piRNAs in the adult hippocampus impairs aNPC differentiation toward a neural fate, induces senescence, and generates reactive glia. Transcripts modulated upon Mili depletion bear sequences complementary or homologous to piRNAs and include repetitive elements and mRNAs encoding essential proteins for proper neurogenesis. Our results provide evidence of a critical role for Mili in maintaining fitness and proper fate of aNPCs, underpinning a possible involvement of the piRNA pathway in brain plasticity and successful aging.

Identification of a tsRNA Contributor to Impaired Diabetic Wound Healing via High Glucose-Induced Endothelial Dysfunction

PURPOSE

Delayed skin healing in diabetic wounds is a major clinical problem. The tRNA-derived small RNAs (tsRNAs) were reported to be associated with diabetes. However, the role of tsRNAs in diabetic wound healing is unclear. Our study was designed to explore the tsRNA expression profile and mine key potential tsRNAs and their mechanism in diabetic wounds.

METHODS

Skin tissues of patients with diabetic foot ulcers and healthy controls were subjected to small RNA sequencing. The role of candidate tsRNA was explored by loss- and gain-of-function experiments in HUVECs.

RESULTS

A total of 55 differentially expressed tsRNAs were identified, including 12 upregulated and 43 downregulated in the diabetes group compared with the control group. These tsRNAs were mainly concentrated in intercellular interactions and neural function regulation in GO terms and enriched in MAPK, insulin, FoxO, calcium, Ras, ErbB, Wnt, T cell receptor, and cGMP-PKG signaling pathways. tRF-Gly-CCC-039 expression was upregulated in vivo and in vitro in the diabetic model. High glucose disturbed endothelial function in HUVECs, and tRF-Gly-CCC-039 mimics further harmed HUVECs function, characterized by the suppression of proliferation, migration, tube formation, and the expression of Coll1a1, Coll4a2, and MMP9. Conversely, the tRF-Gly-CCC-039 inhibitor could attenuate high-glucose-induced endothelial injury to HUVECs.

CONCLUSION

We investigated the tsRNAs expression profile in diabetic foot ulcers and defined the impairment role of tRF-Gly-CCC-039 in endothelial function in HUVECs. This study may provide novel insights into accelerating diabetic skin wound healing.

Absolute quantification of a plasma tRNA-derived fragment for the diagnosis and prognosis of gastric cancer

Background

The transition from a healthy gastric mucosa to gastric cancer is a multi-step process. Early screening can significantly improve the survival rate of gastric cancer patients. A reliable liquid biopsy for gastric cancer prediction is urgently needed and since tRNA-derived fragments (tRFs) are abundant in various body fluids, tRFs are possible new biomarkers for gastric cancer.

Methods

A total of 438 plasma samples from patients with different gastric mucosal lesions as well as healthy individuals were collected. A specific reverse transcription primer, a forward primer, a reverse primer, and a TaqMan probe were designed. A standard curve was constructed and an absolute quantitation method was devised for detection of tRF-33-P4R8YP9LON4VDP in plasma samples of individuals with differing gastric mucosa lesions. Receiver operating characteristic curves were constructed to evaluate the diagnostic values of tRF-33-P4R8YP9LON4VDP for individual with differing gastric mucosa. A Kaplan-Meier curve was established to calculate the prognostic value of tRF-33-P4R8YP9LON4VDP for advanced gastric cancer patients. Finally, a multivariate Cox regression analysis was performed to assess the independent prognostic value of tRF-33-P4R8YP9LON4VDP for advanced gastric cancer patients.

Results

A detection method for plasma tRF-33-P4R8YP9LON4VDP was successfully established. Levels of plasma tRF-33-P4R8YP9LON4VDP were shown to reflect a gradient change from healthy individuals to gastritis patients to early and advanced gastric cancer patients. Significant differences were found among individuals with differing gastric mucosa, with reduced levels of tRF-33-P4R8YP9LON4VDP significantly related to a poor prognosis. tRF-33-P4R8YP9LON4VDP was found to be an independent predictor of an unfavorable survival outcome.

Conclusions

In this study, we developed a quantitative detection method for plasma tRF-33-P4R8YP9LON4VDP that exhibited hypersensitivity, convenience, and specificity. Detection of tRF-33-P4R8YP9LON4VDP was found to be a valuable means by which to monitor different gastric mucosa and to predict patient prognosis.

tRFs and tRNA Halves: Novel Cellular Defenders in Multiple Biological Processes

tRNA fragments derived from angiogenin or Dicer cleavage are referred to as tRNA-derived fragments (tRFs) and tRNA halves. tRFs and tRNA halves have been identified in both eukaryotes and prokaryotes and are precisely cleaved at specific sites on either precursor or mature tRNA transcripts rather than via random degradation. tRFs and tRNA halves are highly involved in regulating transcription and translation in a canonical or non-canonical manner in response to cellular stress. In this review, we summarize the biogenesis and types of tRFs and tRNA halves, clarify the biological functions and molecular mechanisms of tRNA fragments in both physiological and pathological processes with a particular focus on their cytoprotective roles in defending against oxidation and apoptosis, and highlight their potential application as biomarkers in determining cell fate.

tModBase: deciphering the landscape of tRNA modifications and their dynamic changes from epitranscriptome data

tRNA molecules contain dense, abundant modifications that affect tRNA structure, stability, mRNA decoding and tsRNA formation. tRNA modifications and related enzymes are responsive to environmental cues and are associated with a range of physiological and pathological processes. However, there is a lack of resources that can be used to mine and analyse these dynamically changing tRNA modifications. In this study, we established tModBase (https://www.tmodbase.com/) for deciphering the landscape of tRNA modification profiles from epitranscriptome data. We analysed 103 datasets generated with second- and third-generation sequencing technologies and illustrated the misincorporation and termination signals of tRNA modification sites in ten species. We thus systematically demonstrate the modification profiles across different tissues/cell lines and summarize the characteristics of tRNA-associated human diseases. By integrating transcriptome data from 32 cancers, we developed novel tools for analysing the relationships between tRNA modifications and RNA modification enzymes, the expression of 1442 tRNA-derived small RNAs (tsRNAs), and 654 DNA variations. Our database will provide new insights into the features of tRNA modifications and the biological pathways in which they participate.

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.

Dysfunctional tRNA reprogramming and codon-biased translation in cancer

Many cancers hijack translation to increase the synthesis of tumor-driving proteins, the messenger mRNAs of which have specific codon usage patterns. Termed 'codon-biased translation' and originally identified in stress response regulation, this mechanism is supported by diverse studies demonstrating how the 50 RNA modifications of the epitranscriptome, specific tRNAs, and codon-biased mRNAs are used by oncogenic programs to promote proliferation and chemoresistance. The epitranscriptome writers METTL1-WDR4, Elongator complex protein (ELP)1-6, CTU1-2, and ALKBH8-TRM112 illustrate the principal mechanism of codon-biased translation, with gene amplifications, increased RNA modifications, and enhanced tRNA stability promoting cancer proliferation. Furthermore, systems-level analyses of 34 tRNA writers and 493 tRNA genes highlight the theme of tRNA epitranscriptome dysregulation in many cancers and identify candidate tRNA writers, tRNA modifications, and tRNA molecules as drivers of pathological codon-biased translation.

Signature of circulating small non-coding RNAs during early fracture healing in mice

Fracture healing is a complex process with multiple overlapping metabolic and differentiation phases. Small non-coding RNAs are involved in the regulation of fracture healing and their presence in circulation is under current interest due to their obvious value as potential biomarkers. Circulating microRNAs (miRNAs) have been characterized to some extent but the current knowledge on tRNA-derived small RNA fragments (tsRNAs) is relatively scarce, especially in circulation.

In this study, the spectrum of circulating miRNAs and tsRNAs was analysed by next generation sequencing to show their differential expression during fracture healing in vivo. Analysed tsRNA fragments included stress-induced translation interfering tRNA fragments (tiRNAs or tRNA halves) and internal tRNA fragments (i-tRF), within the size range of 28–36 bp. To unveil the expression of these non-coding RNAs, genome-wide analysis was performed on two months old C57BL/6 mice on days 1, 5, 7, 10, and 14 (D1, D5, D7, D10, and D14) after a closed tibial fracture.

Valine isoacceptor tRNA-derived Val-AAC 5′end and Val-CAC 5′end fragments were the major types of 5′end tiRNAs in circulation, comprising about 65 % of the total counts. Their expression was not affected by fracture. After a fracture, the levels of two 5′end tiRNAs Lys-TTT 5′ and Lys-CTT 5′ were decreased and His-GTG 5′ was increased through D1-D14. The level of miR-451a was decreased on the first post-fracture day (D1), whereas miR-328-3p, miR-133a-3p, miR-375-3p, miR-423-5p, and miR-150-5p were increased post-fracture. These data provide evidence on how fracture healing could provoke systemic metabolic effects and further pinpoint the potential of small non-coding RNAs as biomarkers for tissue regeneration.

Transfer RNA-derived small RNAs (tsRNAs): Versatile regulators in cancer

tRNA-derived small RNAs (tsRNAs) represent a novel class of regulatory small non-coding RNAs (sncRNAs), produced by the specific cleavage of transfer RNAs (tRNAs). In recent years, pilot studies one after the other have uncovered the critical roles of tsRNAs in various fundamental biological processes as well as in the development of human diseases including cancer. Based on the newly updated hallmarks of cancer, we provide a comprehensive review regarding the dysregulation, functional implications and complicated molecular mechanisms of tsRNAs in cancer. In addition, the potential technical challenges and future prospects in the fields of tsRNA research are discussed in this review.

tRNA derived small RNAs—Small players with big roles

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

Genetic background and mistranslation frequency determine the impact of mistranslating tRNASerUGG

Transfer RNA variants increase the frequency of mistranslation, the misincorporation of an amino acid not specified by the “standard” genetic code, to frequencies approaching 10% in yeast and bacteria. Cells cope with these variants by having multiple copies of each tRNA isodecoder and through pathways that deal with proteotoxic stress. In this study, we define the genetic interactions of the gene encoding tRNA^Ser_UGG,G26A, which mistranslates serine at proline codons. Using a collection of yeast temperature-sensitive alleles, we identify negative synthetic genetic interactions between the mistranslating tRNA and 109 alleles representing 91 genes, with nearly half of the genes having roles in RNA processing or protein folding and turnover. By regulating tRNA expression, we then compare the strength of the negative genetic interaction for a subset of identified alleles under differing amounts of mistranslation. The frequency of mistranslation correlated with the impact on cell growth for all strains analyzed; however, there were notable differences in the extent of the synthetic interaction at different frequencies of mistranslation depending on the genetic background. For many of the strains, the extent of the negative interaction with tRNA^Ser_UGG,G26A was proportional to the frequency of mistranslation or only observed at intermediate or high frequencies. For others, the synthetic interaction was approximately equivalent at all frequencies of mistranslation. As humans contain similar mistranslating tRNAs, these results are important when analyzing the impact of tRNA variants on disease, where both the individual’s genetic background and the expression of the mistranslating tRNA variant need to be considered.

Strong Selectional Forces Fine-Tune CpG Content in Genes Involved in Neurological Disorders as Revealed by Codon Usage Patterns

Neurodegenerative disorders cause irreversible damage to the neurons and adversely affect the quality of life. Protein misfolding and their aggregation in specific parts of the brain, mitochondrial dysfunction, calcium load, proteolytic stress, and oxidative stress are among the causes of neurodegenerative disorders. In addition, altered metabolism has been associated with neurodegeneration as evidenced by reductions in glutamine and alanine in transient global amnesia patients, higher homocysteine-cysteine disulfide, and lower methionine decline in serum urea have been observed in Alzheimer's disease patients. Neurodegeneration thus appears to be a culmination of altered metabolism. The study's objective is to analyze various attributes like composition, physical properties of the protein, and factors like selectional and mutational forces, influencing codon usage preferences in a panel of genes involved directly or indirectly in metabolism and contributing to neurodegeneration. Various parameters, including gene composition, dinucleotide analysis, Relative synonymous codon usage (RSCU), Codon adaptation index (CAI), neutrality and parity plots, and different protein indices, were computed and analyzed to determine the codon usage pattern and factors affecting it. The correlation of intrinsic protein properties such as the grand average of hydropathicity index (GRAVY), isoelectric point, hydrophobicity, and acidic, basic, and neutral amino acid content has been found to influence codon usage. In genes up to 800 amino acids long, the GC3 content was highly variable, while GC12 content was relatively constant. An optimum CpG content is present in genes to maintain a high expression level as required for genes involved in metabolism. Also observed was a low codon usage bias with a higher protein expression level. Compositional parameters and nucleotides at the second position of codons played essential roles in explaining the extent of bias. Overall analysis indicated that the dominance of selection pressure and compositional constraints and mutational forces shape codon usage.

Role of the Ribonuclease ONCONASE in miRNA Biogenesis and tRNA Processing: Focus on Cancer and Viral Infections

The majority of transcribed RNAs do not codify for proteins, nevertheless they display crucial regulatory functions by affecting the cellular protein expression profile. MicroRNAs (miRNAs) and transfer RNA-derived small RNAs (tsRNAs) are effectors of interfering mechanisms, so that their biogenesis is a tightly regulated process. Onconase (ONC) is an amphibian ribonuclease known for cytotoxicity against tumors and antiviral activity. Additionally, ONC administration in patients resulted in clinical effectiveness and in a well-tolerated feature, at least for lung carcinoma and malignant mesothelioma. Moreover, the ONC therapeutic effects are actually potentiated by cotreatment with many conventional antitumor drugs. This review not only aims to describe the ONC activity occurring either in different tumors or in viral infections but also to analyze the molecular mechanisms underlying ONC pleiotropic and cellular-specific effects. In cancer, data suggest that ONC affects malignant phenotypes by generating tRNA fragments and miRNAs able to downregulate oncogenes expression and upregulate tumor-suppressor proteins. In cells infected by viruses, ONC hampers viral spread by digesting the primer tRNAs necessary for viral DNA replication. In this scenario, new therapeutic tools might be developed by exploiting the action of ONC-elicited RNA derivatives.

Reaction hijacking of tyrosine tRNA synthetase as a new whole-of-life-cycle antimalarial strategy

Aminoacyl transfer RNA (tRNA) synthetases (aaRSs) are attractive drug targets, and we present class I and II aaRSs as previously unrecognized targets for adenosine 5'-monophosphate-mimicking nucleoside sulfamates. The target enzyme catalyzes the formation of an inhibitory amino acid-sulfamate conjugate through a reaction-hijacking mechanism. We identified adenosine 5'-sulfamate as a broad-specificity compound that hijacks a range of aaRSs and ML901 as a specific reagent a specific reagent that hijacks a single aaRS in the malaria parasite Plasmodium falciparum, namely tyrosine RS (PfYRS). ML901 exerts whole-life-cycle-killing activity with low nanomolar potency and single-dose efficacy in a mouse model of malaria. X-ray crystallographic studies of plasmodium and human YRSs reveal differential flexibility of a loop over the catalytic site that underpins differential susceptibility to reaction hijacking by ML901.

Structural variants shape the genomic landscape and clinical outcome of multiple myeloma

Deciphering genomic architecture is key to identifying novel disease drivers and understanding the mechanisms underlying myeloma initiation and progression. In this work, using the CoMMpass dataset, we show that structural variants (SV) occur in a nonrandom fashion throughout the genome with an increased frequency in the t(4;14), RB1, or TP53 mutated cases and reduced frequency in t(11;14) cases. By mapping sites of chromosomal rearrangements to topologically associated domains and identifying significantly upregulated genes by RNAseq we identify both predicted and novel putative driver genes. These data highlight the heterogeneity of transcriptional dysregulation occurring as a consequence of both the canonical and novel structural variants. Further, it shows that the complex rearrangements chromoplexy, chromothripsis and templated insertions are common in MM with each variant having its own distinct frequency and impact on clinical outcome. Chromothripsis is associated with a significant independent negative impact on clinical outcome in newly diagnosed cases consistent with its use alongside other clinical and genetic risk factors to identify prognosis.

tiRNAs: Insights into Their Biogenesis, Functions, and Future Applications in Livestock Research

Transfer RNA (tRNA)-derived small RNAs (tsRNAs) belong to a group of transfer ribonucleic acid (tRNA)-derived fragments that have recently gained interest as molecules with specific biological functions. Their involvement in the regulation of physiological processes and pathological phenotypes suggests molecular roles similar to those of miRNAs. tsRNA biogenesis under specific physiological conditions will offer new perspectives in understanding diseases, and may provide new sources for biological marker design to determine and monitor the health status of farm animals. In this review, we focus on the latest discoveries about tsRNAs and give special attention to molecules initially thought to be mainly associated with tRNA-derived stress-induced RNAs (tiRNAs). We present an outline of their biological functions, offer a collection of useful databases, and discuss future research perspectives and applications in livestock basic and applied research.

tRNAstudio: facilitating the study of human mature tRNAs from deep sequencing datasets

SUMMARY

High-throughput sequencing of transfer RNAs (tRNA-Seq) is a powerful approach to characterize the cellular tRNA pool. Currently, however, analyzing tRNA-Seq datasets requires strong bioinformatics and programming skills. tRNAstudio facilitates the analysis of tRNA-Seq datasets and extracts information on tRNA gene expression, post-transcriptional tRNA modification levels, and tRNA processing steps. Users need only running a few simple bash commands to activate a graphical user interface that allows the easy processing of tRNA-Seq datasets in local mode. Output files include extensive graphical representations and associated numerical tables, and an interactive html summary report to help interpret the data. We have validated tRNAstudio using datasets generated by different experimental methods and derived from human cell lines and tissues that present distinct patterns of tRNA expression, modification and processing.

AVAILABILITY AND IMPLEMENTATION

Freely available at https://github.com/GeneTranslationLab-IRB/tRNAstudio under an open-source GNU GPL v3.0 license.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Data-Independent Acquisition for the Detection of Mononucleoside RNA Modifications by Mass Spectrometry

RNA is dynamically modified in cells by a plethora of chemical moieties to modulate molecular functions and processes. Over 140 modifications have been identified across species and RNA types, with the highest density and diversity of modifications found in tRNA (tRNA). The methods used to identify and quantify these modifications have developed over recent years and continue to advance, primarily in the fields of next-generation sequencing (NGS) and mass spectrometry (MS). Most current NGS methods are limited to antibody-recognized or chemically derivatized modifications and have limitations in identifying multiple modifications simultaneously. Mass spectrometry can overcome both of these issues, accurately identifying a large number of modifications in a single run. Here, we present advances in MS data acquisition for the purpose of RNA modification identification and quantitation. Using this approach, we identified multiple tRNA wobble position modifications in Arabidopsis thaliana that are upregulated in salt-stressed growth conditions and may stabilize translation of salt stress induced proteins. This work presents improvements in methods for studying RNA modifications and introduces a possible regulatory role of wobble position modifications in A. thaliana translation.

tReasure: R-based GUI package analyzing tRNA expression profiles from small RNA sequencing data

BACKGROUND

Recent deep sequencing technologies have proven to be valuable resources to gain insights into the expression profiles of diverse tRNAs. However, despite these technologies, the association of tRNAs with diverse diseases has not been explored in depth because analytical tools are lacking.

RESULTS

We developed a user-friendly tool, tRNA Expression Analysis Software Utilizing R for Easy use (tReasure), to analyze differentially expressed tRNAs (DEtRNAs) from deep sequencing data of small RNAs using R packages. tReasure can quantify individual mature tRNAs, isodecoders, and isoacceptors. By adopting stringent mapping strategies, tReasure supports the precise measurement of mature tRNA read counts. The whole analysis workflow for determining DEtRNAs (uploading FASTQ files, removing adapter sequences and poor-quality reads, mapping and quantifying tRNAs, filtering out low count tRNAs, determining DEtRNAs, and visualizing statistical analysis) can be performed with the tReasure package.

CONCLUSIONS

tReasure is an open-source software available for download at https://treasure.pmrc.re.kr and will be indispensable for users who have little experience with command-line software to explore the biological implication of tRNA expression.

Characterization of tRNA expression profiles in large offspring syndrome

Background

Assisted Reproductive Technologies (ART) use can increase the risk of congenital overgrowth syndromes, such as large offspring syndrome (LOS) in ruminants. Epigenetic variations are known to influence gene expression and differentially methylated regions (DMRs) were previously determined to be associated with LOS in cattle. We observed DMRs overlapping tRNA clusters which could affect tRNA abundance and be associated with tissue specificity or overgrowth. Variations in tRNA expression have been identified in several disease pathways suggesting an important role in the regulation of biological processes. Understanding the role of tRNA expression in cattle offers an opportunity to reveal mechanisms of regulation at the translational level. We analyzed tRNA expression in the skeletal muscle and liver tissues of day 105 artificial insemination-conceived, ART-conceived with a normal body weight, and ART-conceived bovine fetuses with a body weight above the 97^th percentile compared to Control-AI.

Results

Despite the centrality of tRNAs to translation, in silico predictions have revealed dramatic differences in the number of tRNA genes between humans and cattle (597 vs 1,659). Consistent with reports in human, only a fraction of predicted tRNA genes are expressed. We detected the expression of 474 and 487 bovine tRNA genes in the muscle and liver with the remainder being unexpressed. 193 and 198 unique tRNA sequences were expressed in all treatment groups within muscle and liver respectively. In addition, an average of 193 tRNA sequences were expressed within the same treatment group in different tissues. Some tRNA isodecoders were differentially expressed between treatment groups. In the skeletal muscle and liver, we categorized 11 tRNA isoacceptors with undetected expression as well as an isodecoder that was unexpressed in the liver (Ser^GGA). Our results identified variation in the proportion of tRNA gene copies expressed between tissues and differences in the highest contributing tRNA anticodon within an amino acid family due to treatment and tissue type. Out of all amino acid families, roughly half of the most highly expressed tRNA isoacceptors correlated to their most frequent codon in the bovine genome.

Conclusion

Although the number of bovine tRNA genes is nearly triple of that of the tRNA genes in human, there is a shared occurrence of transcriptionally inactive tRNA genes in both species. We detected differential expression of tRNA genes as well as tissue- and treatment- specific tRNA transcripts with unique sequence variations that could modulate translation during protein homeostasis or cellular stress, and give rise to regulatory products targeting genes related to overgrowth in the skeletal muscle and/or tumor development in the liver of LOS individuals. While the absence of certain isodecoders may be relieved by wobble base pairing, missing tRNA species could increase the likelihood of mistranslation or mRNA degradation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08496-7.

Ionizing radiation and chemical oxidant exposure impacts on Cryptococcus neoformans transfer RNAs

Cryptococcus neoformans is a fungus that is able to survive abnormally high levels of ionizing radiation (IR). The radiolysis of water by IR generates reactive oxygen species (ROS) such as H₂O₂ and OH^-. C. neoformans withstands the damage caused by IR and ROS through antioxidant production and enzyme-catalyzed breakdown of ROS. Given these particular cellular protein needs, questions arise whether transfer ribonucleic acids molecules (tRNAs) undergo unique chemical modifications to maintain their structure, stability, and/or function under such environmental conditions. Here, we investigated the effects of IR and H₂O₂ exposure on tRNAs in C. neoformans. We experimentally identified the modified nucleosides present in C. neoformans tRNAs and quantified changes in those modifications upon exposure to oxidative conditions. To better understand these modified nucleoside results, we also evaluated tRNA pool composition in response to the oxidative conditions. We found that regardless of environmental conditions, tRNA modifications and transcripts were minimally affected. A rationale for the stability of the tRNA pool and its concomitant profile of modified nucleosides is proposed based on the lack of codon bias throughout the C. neoformans genome and in particular for oxidative response transcripts. Our findings suggest that C. neoformans can rapidly adapt to oxidative environments as mRNA translation/protein synthesis are minimally impacted by codon bias.

Selective cleavage of ncRNA and antiviral activity by RNase2/EDN in THP1-induced macrophages

RNase2 is the member of the RNaseA family most abundant in macrophages. Here, we knocked out RNase2 in THP-1 cells and analysed the response to Respiratory Syncytial Virus (RSV). RSV induced RNase2 expression, which significantly enhanced cell survival. Next, by cP-RNAseq sequencing, which amplifies the cyclic-phosphate endonuclease products, we analysed the ncRNA population. Among the ncRNAs accumulated in WT vs KO cells, we found mostly tRNA-derived fragments (tRFs) and second miRNAs. Differential sequence coverage identified tRFs from only few parental tRNAs, revealing a predominant cleavage at anticodon and d-loops at U/C (B1) and A (B2) sites. Selective tRNA cleavage was confirmed in vitro using the recombinant protein. Likewise, only few miRNAs were significantly more abundant in WT vs RNase2-KO cells. Complementarily, by screening of a tRF & tiRNA array, we identified an enriched population associated to RNase2 expression and RSV exposure. The results confirm the protein antiviral action and provide the first evidence of its cleavage selectivity on ncRNAs.

Graphical abstract

Structural basis for impaired 5' processing of a mutant tRNA associated with defects in neuronal homeostasis

Understanding and treating neurological disorders are global priorities. Some of these diseases are engendered by mutations that cause defects in the cellular synthesis of transfer RNAs (tRNAs), which function as adapter molecules that translate messenger RNAs into proteins. During tRNA biogenesis, ribonuclease P catalyzes removal of the transcribed sequence upstream of the mature tRNA. Here, we focus on a cytoplasmic tRNA^Arg_UCU that is expressed specifically in neurons and, when harboring a particular point mutation, contributes to neurodegeneration in mice. Our results suggest that this mutation favors stable alternative structures that are not cleaved by mouse ribonuclease P and motivate a paradigm that may help to understand the molecular basis for disease-associated mutations in other tRNAs.

n-Tr20 is a neuron-specific, cytoplasmic transfer RNA^Arg_UCU (tRNA^Arg_UCU) isodecoder that affects seizure susceptibility, neuronal excitability, and translation signaling in mice. In addition, the C50U substitution in n-Tr20 (n-Tr20^C50U), which is found in the widely used C57BL/6J (B6J) inbred mouse line, contributes to neurodegeneration by epistatic interactions with mutations in Gtpbp1 or Gtpbp2, two factors that rescue ribosomal stalling. The brains of B6J mice have high levels of immature and low levels of mature n-Tr20^C50U, implicating defective tRNA biogenesis as the basis for neuronal dysfunction, a hypothesis tested in this study. We demonstrate that partially purified mouse brain ribonuclease P, the endonuclease responsible for 5′ maturation of tRNAs, exhibits at 1 mM magnesium a 20-fold lower apparent cleavage rate for processing the n-Tr20^C50U precursor than the wild-type precursor. Thermal denaturation studies unexpectedly revealed that substituting the native C₅₀-G₆₄ base pair with the U₅₀-G₆₄ wobble pair in n-Tr20^C50U increased the T_m by as much as 8 °C, with the magnitude of the change dependent on [Mg²⁺]. Moreover, results from thermal denaturation, native gel electrophoresis, ¹H nuclear magnetic resonance spectroscopy, and kinetic studies collectively support the presence of stable alternative folds for n-Tr20^C50U that may also be sampled transiently and in low abundance by the wild-type tRNA. These findings suggest that mutation-driven restructuring could foster nonnative folds and that conformational toggling of tRNAs poses an intrinsic risk for dysfunction when point mutations selectively stabilize nonnative states. This model may be applicable for understanding the molecular basis of other diseases that are associated with tRNA mutations.

Disorders and roles of tsRNA, snoRNA, snRNA and piRNA in cancer

Most small non-coding RNAs (sncRNAs) with regulatory functions are encoded by majority sequences in the human genome, and the emergence of high-throughput sequencing technology has greatly expanded our understanding of sncRNAs. sncRNAs are composed of a variety of RNAs, including tRNA-derived small RNA (tsRNA), small nucleolar RNA (snoRNA), small nuclear RNA (snRNA), PIWI-interacting RNA (piRNA), etc. While for some, sncRNAs' implication in several pathologies is now well established, the potential involvement of tsRNA, snoRNA, snRNA and piRNA in human diseases is only beginning to emerge. Recently, accumulating pieces of evidence demonstrate that tsRNA, snoRNA, snRNA and piRNA play an important role in many biological processes, and their dysregulation is closely related to the progression of cancer. Abnormal expression of tsRNA, snoRNA, snRNA and piRNA participates in the occurrence and development of tumours through different mechanisms, such as transcriptional inhibition and post-transcriptional regulation. In this review, we describe the research progress in the classification, biogenesis and biological function of tsRNA, snoRNA, snRNA and piRNA. Moreover, we emphasised their dysregulation and mechanism of action in cancer and discussed their potential as diagnostic and prognostic biomarkers or therapeutic targets.

Plasma cell-free RNA characteristics in COVID-19 patients

The pathogenesis of COVID-19 is still elusive, which impedes disease progression prediction, differential diagnosis, and targeted therapy. Plasma cell-free RNAs (cfRNAs) carry unique information from human tissue and thus could point to resourceful solutions for pathogenesis and host-pathogen interactions. Here, we performed a comparative analysis of cfRNA profiles between COVID-19 patients and healthy donors using serial plasma. Analyses of the cfRNA landscape, potential gene regulatory mechanisms, dynamic changes in tRNA pools upon infection, and microbial communities were performed. A total of 380 cfRNA molecules were up-regulated in all COVID-19 patients, of which seven could serve as potential biomarkers (AUC > 0.85) with great sensitivity and specificity. Antiviral (NFKB1A, IFITM3, and IFI27) and neutrophil activation (S100A8, CD68, and CD63)–related genes exhibited decreased expression levels during treatment in COVID-19 patients, which is in accordance with the dynamically enhanced inflammatory response in COVID-19 patients. Noncoding RNAs, including some microRNAs (let 7 family) and long noncoding RNAs (GJA9-MYCBP) targeting interleukin (IL6/IL6R), were differentially expressed between COVID-19 patients and healthy donors, which accounts for the potential core mechanism of cytokine storm syndromes; the tRNA pools change significantly between the COVID-19 and healthy group, leading to the accumulation of SARS-CoV-2 biased codons, which facilitate SARS-CoV-2 replication. Finally, several pneumonia-related microorganisms were detected in the plasma of COVID-19 patients, raising the possibility of simultaneously monitoring immune response regulation and microbial communities using cfRNA analysis. This study fills the knowledge gap in the plasma cfRNA landscape of COVID-19 patients and offers insight into the potential mechanisms of cfRNAs to explain COVID-19 pathogenesis.

Adult Neural Stem Cell Regulation by Small Non-coding RNAs: Physiological Significance and Pathological Implications

The adult neurogenic niches are complex multicellular systems, receiving regulatory input from a multitude of intracellular, juxtacrine, and paracrine signals and biological pathways. Within the niches, adult neural stem cells (aNSCs) generate astrocytic and neuronal progeny, with the latter predominating in physiological conditions. The new neurons generated from this neurogenic process are functionally linked to memory, cognition, and mood regulation, while much less is known about the functional contribution of aNSC-derived newborn astrocytes and adult-born oligodendrocytes. Accumulating evidence suggests that the deregulation of aNSCs and their progeny can impact, or can be impacted by, aging and several brain pathologies, including neurodevelopmental and mood disorders, neurodegenerative diseases, and also by insults, such as epileptic seizures, stroke, or traumatic brain injury. Hence, understanding the regulatory underpinnings of aNSC activation, differentiation, and fate commitment could help identify novel therapeutic avenues for a series of pathological conditions. Over the last two decades, small non-coding RNAs (sncRNAs) have emerged as key regulators of NSC fate determination in the adult neurogenic niches. In this review, we synthesize prior knowledge on how sncRNAs, such as microRNAs (miRNAs) and piwi-interacting RNAs (piRNAs), may impact NSC fate determination in the adult brain and we critically assess the functional significance of these events. We discuss the concepts that emerge from these examples and how they could be used to provide a framework for considering aNSC (de)regulation in the pathogenesis and treatment of neurological diseases.

The nuclear and cytoplasmic activities of RNA polymerase III, and an evolving transcriptome for surveillance

A 1969 report that described biochemical and activity properties of the three eukaryotic RNA polymerases revealed Pol III as highly distinguishable, even before its transcripts were identified. Now known to be the most complex, Pol III contains several stably-associated subunits referred to as built-in transcription factors (BITFs) that enable highly efficient RNA synthesis by a unique termination-associated recycling process. In vertebrates, subunit RPC7(α/β) can be of two forms, encoded by POLR3G or POLR3GL, with differential activity. Here we review promoter-dependent transcription by Pol III as an evolutionary perspective of eukaryotic tRNA expression. Pol III also provides nonconventional functions reportedly by promoter-independent transcription, one of which is RNA synthesis from DNA 3'-ends during repair. Another is synthesis of 5'ppp-RNA signaling molecules from cytoplasmic viral DNA in a pathway of interferon activation that is dysfunctional in immunocompromised patients with mutations in Pol III subunits. These unconventional functions are also reviewed, including evidence that link them to the BITF subunits. We also review data on a fraction of the human Pol III transcriptome that evolved to include vault RNAs and snaRs with activities related to differentiation, and in innate immune and tumor surveillance. The Pol III of higher eukaryotes does considerably more than housekeeping.

Identification of conserved slow codons that are important for protein expression and function

ABSTRASTDue to the redundancy of the genetic code most amino acids are encoded by several 'synonymous' codons. These codons are used unevenly, and each organism demonstrates its own unique codon usage bias, where the 'preferred' codons are associated with tRNAs that are found in high concentrations. Therefore, for decades, the prevailing view had been that preferred and non-preferred codons are linked to high or slow translation rates, respectively.However, this simplified view is contrasted by the frequent failures of codon-optimization efforts and by evidence of non-preferred (i.e. 'slow') codons having specific roles important for efficient production of functional proteins. One such specific role of slower codons is the regulation of co-translational protein folding, a complex biophysical process that is very challenging to model or to measure.Here, we combined a genome-wide approach with experiments to investigate the role of slow codons in protein production and co-translational folding. We analysed homologous gene groups from divergent bacteria and identified positions of inter-species conservation of bias towards slow codons. We then generated mutants where the conserved slow codons are substituted with 'fast' ones, and experimentally studied the effects of these codon substitutions. Using cellular and biochemical approaches we find that at certain locations, slow-to-fast codon substitutions reduce protein expression, increase protein aggregation, and impair protein function.This report provides an approach for identifying functionally relevant regions with slower codons and demonstrates that such codons are important for protein expression and function.

THUMPD3-TRMT112 is a m2G methyltransferase working on a broad range of tRNA substrates

Post-transcriptional modifications affect tRNA biology and are closely associated with human diseases. However, progress on the functional analysis of tRNA modifications in metazoans has been slow because of the difficulty in identifying modifying enzymes. For example, the biogenesis and function of the prevalent N2-methylguanosine (m2G) at the sixth position of tRNAs in eukaryotes has long remained enigmatic. Herein, using a reverse genetics approach coupled with RNA-mass spectrometry, we identified that THUMP domain-containing protein 3 (THUMPD3) is responsible for tRNA: m2G6 formation in human cells. However, THUMPD3 alone could not modify tRNAs. Instead, multifunctional methyltransferase subunit TRM112-like protein (TRMT112) interacts with THUMPD3 to activate its methyltransferase activity. In the in vitro enzymatic assay system, THUMPD3-TRMT112 could methylate all the 26 tested G6-containing human cytoplasmic tRNAs by recognizing the characteristic 3'-CCA of mature tRNAs. We also showed that m2G7 of tRNATrp was introduced by THUMPD3-TRMT112. Furthermore, THUMPD3 is widely expressed in mouse tissues, with an extremely high level in the testis. THUMPD3-knockout cells exhibited impaired global protein synthesis and reduced growth. Our data highlight the significance of the tRNA: m2G6/7 modification and pave a way for further studies of the role of m2G in sperm tRNA derived fragments.