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

The subcellular distribution of miRNA isoforms, tRNA-derived fragments, and rRNA-derived fragments depends on nucleotide sequence and cell type

BACKGROUND

MicroRNA isoforms (isomiRs), tRNA-derived fragments (tRFs), and rRNA-derived fragments (rRFs) represent most of the small non-coding RNAs (sncRNAs) found in cells. Members of these three classes modulate messenger RNA (mRNA) and protein abundance and are dysregulated in diseases. Experimental studies to date have assumed that the subcellular distribution of these molecules is well-understood, independent of cell type, and the same for all isoforms of a sncRNA.

RESULTS

We tested these assumptions by investigating the subcellular distribution of isomiRs, tRFs, and rRFs in biological replicates from three cell lines from the same tissue and same-sex donors that model the same cancer subtype. In each cell line, we profiled the isomiRs, tRFs, and rRFs in the nucleus, cytoplasm, whole mitochondrion (MT), mitoplast (MP), and whole cell. Using a rigorous mathematical model we developed, we accounted for cross-fraction contamination and technical errors and adjusted the measured abundances accordingly. Analyses of the adjusted abundances show that isomiRs, tRFs, and rRFs exhibit complex patterns of subcellular distributions. These patterns depend on each sncRNA's exact sequence and the cell type. Even in the same cell line, isoforms of the same sncRNA whose sequences differ by a few nucleotides (nts) can have different subcellular distributions.

CONCLUSIONS

SncRNAs with similar sequences have different subcellular distributions within and across cell lines, suggesting that each isoform could have a different function. Future computational and experimental studies of isomiRs, tRFs, and rRFs will need to distinguish among each molecule's various isoforms and account for differences in each isoform's subcellular distribution in the cell line at hand. While the findings add to a growing body of evidence that isomiRs, tRFs, rRFs, tRNAs, and rRNAs follow complex intracellular trafficking rules, further investigation is needed to exclude alternative explanations for the observed subcellular distribution of sncRNAs.

tsRNA modifications: An emerging layer of biological regulation in disease

BACKGROUND

Transfer RNA (tRNA)-derived small RNA (tsRNA) represents an important and increasingly valued type of small non-coding RNA (sncRNA). The investigation of tRNA and tsRNA modification crosswalks has not only provided novel insights into the information and functions of tsRNA, but has also expanded the diversity and complexity of the tsRNA biological regulation network.

AIM OF REVIEW

Comparing with other sncRNAs, tsRNA biogenesis show obvious correlation with RNA modifications from mature tRNA and harbor various tRNA modifications. In this review, we aim to present the current aspect of tsRNA modifications and that modified tsRNA shape different regulatory mechanisms in physiological and pathological processes.

KEY SCIENTIFIC CONCEPTS OF REVIEW

Strategies for studying tsRNA mechanisms include its specific generation and functional effects induced by sequence/RNA modification/secondary structure. tsRNAs could harbor more than one tRNA modifications such as 5-methylcytosine (m5C), N1-methyladenosine (m1A), pseudouridine (Ψ) and N7-methylguanosine (m7G). This review consolidates the current knowledge of tRNA modification regulating tsRNA biogenesis, outlines the functional roles of various modified tsRNA and highlights their specific contributions in various disease pathogenesis. Therefore, the improvement of tsRNA modification detection technology and the introduction of experimental methods of tsRNA modification are conducive to further broadening the understanding of tsRNA function at the level of RNA modification.

A novel tRNA-derived fragment tRF-17-18VBY9M works as a potential diagnostic biomarker for gastric cancer

BACKGROUND

Gastric cancer (GC) is one of the most prevalent malignant tumors worldwide. The low effectiveness of common biomarkers for the detection of early GC makes it essential to seek new biomarkers to improve diagnostic efficacy. tsRNAs (transfer RNA-derived small RNAs) are related to the growth of malignant tumors. In this article, we focused on whether tsRNAs may be employed as biomarkers for GC.

METHODS

tRF-17-18VBY9M was screened in the tsRFun database as a research object. The methodological efficacy of tRF-17-18VBY9M was evaluated using Sanger sequencing, agarose gel electrophoresis assays, and gradient dilution. The χ2 test was applied to assess the interaction between tRF-17-18VBY9M expression and clinicopathologic characteristics. The receiver operating characteristic (ROC) curve was utilized to investigate the clinical efficiency of tRF-17-18VBY9M in GC.

RESULTS

The Chi-square test demonstrated that high-expressed tRF-17-18VBY9M was closely associated with the T stage, tumor node metastasis stage (TNM), lymph node metastasis, and neurological/vascular invasion. ROC curve analysis revealed that the diagnostic value of tRF-17-18VBY9M in GC was superior to carcinoembryonic antigen (CEA), carbohydrate antigen 199 (CA199), and carbohydrate antigen 724 (CA724).

CONCLUSION

tRF-17-18VBY9M is up-regulated in both GC sera and tissues. Differential tRF-17-18VBY9M expression distinguishes GC patients from healthy donors and gastritis patients, which suggests tRF-17-18VBY9M could act as a diagnostic biomarker in GC.

Small RNA-seq and clinical evaluation of tRNA-derived fragments in multiple myeloma: Loss of mitochondrial i-tRFHisGTG results in patients' poor treatment outcome

Despite the substantial progress in multiple myeloma (MM) therapy nowadays, treatment resistance and disease relapse remain major clinical hindrances. Herein, we have investigated tRNA-derived fragment (tRF) profiles in MM and precursor stages (smoldering MM/sMM; monoclonal gammopathy of undetermined significance/MGUS), aiming to unveil potential MM-related tRFs in ameliorating MM prognosis and risk stratification. Small RNA-seq was performed to profile tRFs in bone marrow CD138+ plasma cells, revealing the significant deregulation of the mitochondrial internal tRFHisGTG (mt-i-tRFHisGTG) in MM versus sMM/MGUS. The screening cohort of the study consisted of 147 MM patients, and mt-i-tRFHisGTG levels were quantified by RT-qPCR. Disease progression was assessed as clinical end-point for survival analysis, while internal validation was performed by bootstrap and decision curve analyses. Screening cohort analysis highlighted the potent association of reduced mt-i-tRFHisGTG levels with patients' bone disease (p = 0.010), osteolysis (p = 0.023) and with significantly higher risk for short-term disease progression following first-line chemotherapy, independently of patients' clinical data (HR = 1.954; p = 0.036). Additionally, mt-i-tRFHisGTG-fitted multivariate models led to superior risk stratification of MM patients' treatment outcome and prognosis compared to disease-established markers. Notably, our study highlighted mt-i-tRFHisGTG loss as a powerful independent indicator of post-treatment progression of MM patients, leading to superior risk stratification of patients' treatment outcome.

Quantitative and qualitative detection of tRNAs, tRNA halves and tRFs in human cancer samples: Molecular grounds for biomarker development and clinical perspectives

Transfer RNAs (tRNAs) are small non-coding RNAs playing a central role during protein synthesis. Besides translation, growing evidence suggests that in many contexts, precursor or mature tRNAs can also be processed into smaller fragments playing many non-canonical regulatory roles in different biological pathways with oncogenic relevance. Depending on the source, these molecules can be classified as tRNA halves (also known as tiRNAs) or tRNA-derived fragments (tRFs), and furtherly divided into 5'-tRNA and 3'-tRNA halves, or tRF-1, tRF-2, tRF-3, tRF-5, and i-tRF, respectively. Unlike DNA and mRNA, high-throughput sequencing of tRNAs is challenging, because of technical limitations of currently developed sequencing methods. In recent years, different sequencing approaches have been proposed allowing the quantification and identification of an increasing number of tRNA fragments with critical functions in distinct physiological and pathophysiological processes. In the present review, we discussed pros and cons of recent advances in different sequencing methods, also introducing the expanding repertoire of bioinformatics tool and resources specifically focused on tRNA research and discussing current issues in the study of these small RNA molecules. Furthermore, we discussed the potential value of tRNA fragments as diagnostic and prognostic biomarkers for different types of cancers.

tRF-1:30-Gly-CCC-3 inhibits thyroid cancer via binding to PC and modulating metabolic reprogramming

tRFs and tiRNAs (tRNA-derived fragments) are an emerging class of small noncoding RNAs produced by the precise shearing of tRNAs in response to specific stimuli. They have been reported to regulate the pathological processes of numerous human cancers. However, the biofunction of tRFs and tiRNAs in the development and progression of papillary thyroid cancer (PTC) has not been reported yet. In this study, we aimed to explore the biological roles of tRFs and tiRNAs in PTC and discovered that a novel 5'tRNA-derived fragment called tRF-1:30-Gly-CCC-3 (tRF-30) was markedly down-regulated in PTC tissues and cell lines. Functionally, tRF-30 inhibited the proliferation and invasion of PTC cells. Mechanistically, tRF-30 directly bound to the biotin-dependent enzyme pyruvate carboxylase (PC), downregulated its protein level, interfered with the TCA cycle intermediate anaplerosis, and thus affected metabolic reprogramming and PTC progression. These findings revealed a novel regulatory mechanism for tRFs and a potential therapeutic target for PTC.

tRFUniverse: A comprehensive resource for the interactive analyses of tRNA-derived ncRNAs in human cancer

tRNA-derived ncRNAs are a heterogeneous class of non-coding RNAs recently proposed to be active regulators of gene expression and be involved in many diseases, including cancer. Consequently, several online resources on tRNA-derived ncRNAs have been released. Although interesting, such resources present only basic features and do not adequately exploit the wealth of knowledge available about tRNA-derived ncRNAs. Therefore, we introduce tRFUniverse, a novel online resource for the analysis of tRNA-derived ncRNAs in human cancer. tRFUniverse presents an extensive collection of classes of tRNA-derived ncRNAs analyzed across all the TCGA and TARGET tumor cohorts, NCI-60 cell lines, and biological fluids. Moreover, public AGO CLASH/CLIP-Seq data were analyzed to identify the molecular interactions between tRNA-derived ncRNAs and other transcripts. Importantly, tRFUniverse combines in a single resource a comprehensive set of features that we believe may be helpful to investigate the involvement of tRNA-derived ncRNAs in cancer biology.

MINRbase: a comprehensive database of nuclear- and mitochondrial-ribosomal-RNA-derived fragments (rRFs)

We describe the Mitochondrial and Nuclear rRNA fragment database (MINRbase), a knowledge repository aimed at facilitating the study of ribosomal RNA-derived fragments (rRFs). MINRbase provides interactive access to the profiles of 130 238 expressed rRFs arising from the four human nuclear rRNAs (18S, 5.8S, 28S, 5S), two mitochondrial rRNAs (12S, 16S) or four spacers of 45S pre-rRNA. We compiled these profiles by analyzing 11 632 datasets, including the GEUVADIS and The Cancer Genome Atlas (TCGA) repositories. MINRbase offers a user-friendly interface that lets researchers issue complex queries based on one or more criteria, such as parental rRNA identity, nucleotide sequence, rRF minimum abundance and metadata keywords (e.g. tissue type, disease). A 'summary' page for each rRF provides a granular breakdown of its expression by tissue type, disease, sex, ancestry and other variables; it also allows users to create publication-ready plots at the click of a button. MINRbase has already allowed us to generate support for three novel observations: the internal spacers of 45S are prolific producers of abundant rRFs; many abundant rRFs straddle the known boundaries of rRNAs; rRF production is regimented and depends on 'personal attributes' (sex, ancestry) and 'context' (tissue type, tissue state, disease). MINRbase is available at https://cm.jefferson.edu/MINRbase/.

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.

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.

Characterization of all small RNAs in and comparisons across cultured megakaryocytes and platelets of healthy individuals and COVID-19 patients

BACKGROUND

The small noncoding RNAs (sncRNAs) in megakaryocytes (MKs) and platelets are not well characterized. Neither is the impact of SARS-CoV-2 infection on the sncRNAs of platelets.

OBJECTIVES

To investigate the sorting of MK sncRNAs into platelets, and the differences in the platelet sncRNAomes of healthy donors (HDs) and COVID-19 patients.

METHODS

We comprehensively profiled sncRNAs from MKs cultured from cord blood-derived CD34+ cells, platelets from HDs, and platelets from patients with moderate and severe SARS-CoV-2 infection. We also comprehensively profiled Argonaute (AGO)-bound sncRNAs from the cultured MKs.

RESULTS

We characterized the sncRNAs in MKs and platelets and can account for ∼95% of all sequenced reads. We found that MKs primarily comprise microRNA isoforms (isomiRs), tRNA-derived fragments (tRFs), rRNA-derived fragments (rRFs), and Y RNA-derived fragments (yRFs) in comparable abundances. The platelets of HDs showed a skewed distribution by comparison: 56.7% of all sncRNAs are yRFs, 34.4% are isomiRs, and <2.0% are tRFs and rRFs. Most isomiRs in MKs and platelets are either noncanonical, nontemplated, or both. When comparing MKs and platelets from HDs, we found numerous isomiRs, tRFs, rRFs, and yRFs showing opposite enrichments or depletions, including molecules from the same parental miRNA arm, tRNA, rRNA, or Y RNA. The sncRNAome of platelets from patients with COVID-19 is skewed compared to that of HDs with only 19.8% of all sncRNAs now being yRFs, isomiRs increasing to 63.6%, and tRFs and rRFs more than tripling their presence to 6.1%.

CONCLUSION

The sncRNAomes of MKs and platelets are very rich and more complex than it has been believed. The evidence suggests complex mechanisms that sort MK sncRNAs into platelets. SARS-CoV-2 infection acutely alters the contents of platelets by changing the relative proportions of their sncRNAs.

tRF-21-LNK8KEP1B as a potential novel diagnostic biomarker for enthesitis-related arthritis

OBJECTIVE

tRNA-derived fragments (tRFs) play crucial roles in the progression of various diseases, and widely distribute in human tissues, including blood and urine. The diagnosis of enthesitis-related arthritis (ERA) is based on the observation of clinical manifestations. Therefore, we aimed to investigate whether serum tRFs can be used as diagnostic markers for ERA.

METHODS

Serum was collected from children admitted to the Children's Hospital Affiliated with Nanjing Medical University between February 2022 to October 2022. The expression profiles of tRFs in the serum of ERA patients (n = 5) and healthy controls (HCs; n = 5) were investigated using small RNA high-throughput sequencing. The level and diagnostic value of tRF-21-LNK8KEP1B were evaluated by real-time quantitative PCR in serum samples from 30 ERA patients and 31 HCs. The specificity and sensitivity of tRFs were determined using receiver operating characteristic analyses. Bioinformatics analysis was performed to explore and identify the potential biological pathways induced by tRFs.

RESULTS

Ninety-eight upregulated and 63 downregulated tRFs were identified in the serum. We selected tRF-21-LNK8KEP1B as a candidate marker using KEGG pathway enrichment and PCR validation. tRF-21-LNK8KEP1B was substantially increased in the serum of ERA patients compared with that in HCs. The area under the curve (AUC) for tRF-21-LNK8KEP1B in the ERA group was 0.7849.

CONCLUSIONS

Collectively, we demonstrated the promising role of serum tRF-21-LNK8KEP1B -levels as a diagnostic biomarker for ERA.

XRN2 suppresses aberrant entry of tRNA trailers into argonaute in humans and Arabidopsis

MicroRNAs (miRNAs) are a well-characterized class of small RNAs (sRNAs) that regulate gene expression post-transcriptionally. miRNAs function within a complex milieu of other sRNAs of similar size and abundance, with the best characterized being tRNA fragments or tRFs. The mechanism by which the RNA-induced silencing complex (RISC) selects for specific sRNAs over others is not entirely understood in human cells. Several highly expressed tRNA trailers (tRF-1s) are strikingly similar to microRNAs in length but are generally excluded from the microRNA effector pathway. This exclusion provides a paradigm for identifying mechanisms of RISC selectivity. Here, we show that 5' to 3' exoribonuclease XRN2 contributes to human RISC selectivity. Although highly abundant, tRF-1s are highly unstable and degraded by XRN2 which blocks tRF-1 accumulation in RISC. We also find that XRN mediated degradation of tRF-1s and subsequent exclusion from RISC is conserved in plants. Our findings reveal a conserved mechanism that prevents aberrant entry of a class of highly produced sRNAs into Ago2.

Stress Response in Entamoeba histolytica Is Associated with Robust Processing of tRNA to tRNA Halves

tRNA-derived fragments have been reported in many different organisms and have diverse cellular roles, such as regulating gene expression, inhibiting protein translation, silencing transposable elements, and modulating cell proliferation. In particular, tRNA halves, a class of tRNA fragments produced by the cleavage of tRNAs in the anti-codon loop, have been widely reported to accumulate under stress and regulate translation in cells. Here, we report the presence of tRNA-derived fragments in Entamoeba, with tRNA halves being the most abundant. We further established that tRNA halves accumulate in the parasites upon different stress stimuli such as oxidative stress, heat shock, and serum starvation. We also observed differential expression of tRNA halves during developmental changes of trophozoite-to-cyst conversion, with various tRNA halves accumulating during early encystation. In contrast to other systems, the stress response does not appear to be mediated by a few specific tRNA halves, as multiple tRNAs appear to be processed during the various stresses. Furthermore, we identified some tRNA-derived fragments associated with Entamoeba Argonaute proteins, EhAgo2-2 and EhAgo2-3, which have a preference for different tRNA-derived fragment species. Finally, we show that tRNA halves are packaged inside extracellular vesicles secreted by amoebas. The ubiquitous presence of tRNA-derived fragments, their association with the Argonaute proteins, and the accumulation of tRNA halves during multiple different stresses, including encystation, suggest a nuanced level of gene expression regulation mediated by different tRNA-derived fragments in Entamoeba. IMPORTANCE In the present study, we report for the first time the presence of tRNA-derived fragments in Entamoeba. tRNA-derived fragments were identified by bioinformatics analyses of small-RNA sequencing data sets from the parasites and also confirmed experimentally. We found that tRNA halves accumulated in parasites exposed to environmental stress or during the developmental process of encystation. We also found that shorter tRNA-derived fragments are bound to Entamoeba Argonaute proteins, indicating that they may have a potential role in the Argonaute-mediated RNA-interference pathway, which mediates robust gene silencing in Entamoeba. We noticed that in response to heat shock, the protein translation levels were elevated in the parasites. This effect was reversed in the presence of an analog of leucine, which also reduced the levels of the tRNA halves in the stressed cells. Our results suggest that tRNA-derived fragments in Entamoeba have a possible role in regulating gene expression during environmental stress.

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.

Cerebrospinal fluid and blood profiles of transfer RNA fragments show age, sex, and Parkinson's disease-related changes

Transfer RNA fragments (tRFs) have recently been shown to be an important family of small regulatory RNAs with diverse functions. Recent reports have revealed modified tRF blood levels in a number of nervous system conditions including epilepsy, ischemic stroke, and neurodegenerative diseases, but little is known about tRF levels in the cerebrospinal fluid (CSF). To address this issue, we studied age, sex, and Parkinson's disease (PD) effects on the distributions of tRFs in the CSF and blood data of healthy controls and PD patients from the NIH and the Parkinson's Progression Markers Initiative (PPMI) small RNA-seq datasets. We discovered that long tRFs are expressed in higher levels in the CSF than in the blood. Furthermore, the CSF showed a pronounced age-associated decline in the level of tRFs cleaved from the 3'-end and anti-codon loop of the parental tRNA (3'-tRFs, i-tRFs), and more pronounced profile differences than the blood profiles between the sexes. In comparison, we observed moderate age-related elevation of blood 3'-tRF levels. In addition, distinct sets of tRFs in the CSF and in the blood segregated PD patients from controls. Finally, we found enrichment of tRFs predicted to target cholinergic mRNAs (Cholino-tRFs) among mitochondrial-originated tRFs, raising the possibility that the neurodegeneration-related mitochondrial impairment in PD patients may lead to deregulation of their cholinergic tone. Our findings demonstrate that the CSF and blood tRF profiles are distinct and that the CSF tRF profiles are modified in a sex-, age-, and disease-related manner, suggesting that they reflect the inter-individual cerebral differences and calling for incorporating this important subset of small RNA regulators into future studies.

Insight on Non-Coding RNAs from Biofluids in Ovarian Tumors

Ovarian tumors are the most frequent adnexal mass, raising diagnostic and therapeutic issues linked to a large spectrum of tumors, with a continuum from benign to malignant. Thus far, none of the available diagnostic tools have proven efficient in deciding strategy, and no consensus exists on the best strategy between "single test", "dual testing", "sequential testing", "multiple testing options" and "no testing". In addition, there is a need for prognostic tools such as biological markers of recurrence and theragnostic tools to detect women not responding to chemotherapy in order to adapt therapies. Non-coding RNAs are classified as small or long based on their nucleotide count. Non-coding RNAs have multiple biological functions such as a role in tumorigenesis, gene regulation and genome protection. These ncRNAs emerge as new potential tools to differentiate benign from malignant tumors and to evaluate prognostic and theragnostic factors. In the specific setting of ovarian tumors, the goal of the present work is to offer an insight into the contribution of biofluid non-coding RNAs (ncRNA) expression.

mtR_find: A Parallel Processing Tool to Identify and Annotate RNAs Derived from the Mitochondrial Genome

RNAs originating from mitochondrial genomes are abundant in transcriptomic datasets produced by high-throughput sequencing technologies, primarily in short-read outputs. Specific features of mitochondrial small RNAs (mt-sRNAs), such as non-templated additions, presence of length variants, sequence variants, and other modifications, necessitate the need for the development of an appropriate tool for their effective identification and annotation. We have developed mtR_find, a tool to detect and annotate mitochondrial RNAs, including mt-sRNAs and mitochondria-derived long non-coding RNAs (mt-lncRNA). mtR_find uses a novel method to compute the count of RNA sequences from adapter-trimmed reads. When analyzing the published datasets with mtR_find, we identified mt-sRNAs significantly associated with the health conditions, such as hepatocellular carcinoma and obesity, and we discovered novel mt-sRNAs. Furthermore, we identified mt-lncRNAs in early development in mice. These examples show the immediate impact of miR_find in extracting a novel biological information from the existing sequencing datasets. For benchmarking, the tool has been tested on a simulated dataset and the results were concordant. For accurate annotation of mitochondria-derived RNA, particularly mt-sRNA, we developed an appropriate nomenclature. mtR_find encompasses the mt-ncRNA transcriptomes in unpreceded resolution and simplicity, allowing re-analysis of the existing transcriptomic databases and the use of mt-ncRNAs as diagnostic or prognostic markers in the field of medicine.

Emerging roles of tRNA-derived fragments in cancer

tRNA-derived fragments (tRFs) are an emerging category of small non-coding RNAs that are generated from cleavage of mature tRNAs or tRNA precursors. The advance in high-throughput sequencing has contributed to the identification of increasing number of tRFs with critical functions in distinct physiological and pathophysiological processes. tRFs can regulate cell viability, differentiation, and homeostasis through multiple mechanisms and are thus considered as critical regulators of human diseases including cancer. In addition, increasing evidence suggest the extracellular tRFs may be utilized as promising diagnostic and prognostic biomarkers for cancer liquid biopsy. In this review, we focus on the biogenesis, classification and modification of tRFs, and summarize the multifaceted functions of tRFs with an emphasis on the current research status and perspectives of tRFs in cancer.

Roles of mitochondrial genetics in cancer metastasis

The contributions of mitochondria to cancer have been recognized for decades. However, the focus on the metabolic role of mitochondria and the diminutive size of the mitochondrial genome compared to the nuclear genome have hindered discovery of the roles of mitochondrial genetics in cancer. This review summarizes recent data demonstrating the contributions of mitochondrial DNA (mtDNA) copy-number variants (CNVs), somatic mutations, and germline polymorphisms to cancer initiation, progression, and metastasis. The goal is to summarize accumulating data to establish a framework for exploring the contributions of mtDNA to neoplasia and metastasis.

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.

MIR retrotransposons link the epigenome and the transcriptome of coding genes in acute myeloid leukemia

DNMT3A and IDH1/2 mutations combinatorically regulate the transcriptome and the epigenome in acute myeloid leukemia; yet the mechanisms of this interplay are unknown. Using a systems approach within topologically associating domains, we find that genes with significant expression-methylation correlations are enriched in signaling and metabolic pathways. The common denominator across these methylation-regulated genes is the density in MIR retrotransposons of their introns. Moreover, a discrete number of CpGs overlapping enhancers are responsible for regulating most of these genes. Established mouse models recapitulate the dependency of MIR-rich genes on the balanced expression of epigenetic modifiers, while projection of leukemic profiles onto normal hematopoiesis ones further consolidates the dependencies of methylation-regulated genes on MIRs. Collectively, MIR elements on genes and enhancers are susceptible to changes in DNA methylation activity and explain the cooperativity of proteins in this pathway in normal and malignant hematopoiesis.

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.

Characterization of novel small non-coding RNAs and their modifications in bladder cancer using an updated small RNA-seq workflow

Background: Bladder cancer (BLCA) is one of the most common cancer types worldwide. The disease is responsible for about 200,000 deaths annually, thus improved diagnostics and therapy is needed. A large body of evidence reveal that small RNAs of less than 40 nucleotides may act as tumor suppressors, oncogenes, and disease biomarkers, with a major focus on microRNAs. However, the role of other families of small RNAs is not yet deciphered. Recent results suggest that small RNAs and their modification status, play a role in BLCA development and are promising biomarkers due to their high abundance in the exomes and body fluids (including urine). Moreover, free modified nucleosides have been detected at elevated levels from the urine of BLCA patients. A genome-wide view of small RNAs, and their modifications, will help pinpoint the molecules that could be used as biomarker or has important biology in BLCA development. Methods: BLCA tumor tissue specimens were obtained from 12 patients undergoing transurethral resection of non-muscle invasive papillary urothelial carcinomas. Genome-wide profiling of small RNAs less than 40 bases long was performed by a modified protocol with TGIRT (thermostable group II reverse transcriptase) to identify novel small RNAs and their modification status. Results: Comprehensive analysis identified not only microRNAs. Intriguingly, 57 ± 15% (mean ± S.D.) of sequencing reads mapped to non-microRNA-small RNAs including tRNA-derived fragments (tRFs), ribosomal RNA-derived fragments (rRFs) and YRNA-derived fragments (YRFs). Misincorporation (mismatch) sites identified potential base modification positions on the small RNAs, especially on tRFs, corresponding to m1A (N1-methyladenosine), m1G (N1-methylguanosine) and m2 2G (N2, N2-dimethylguanosine). We also detected mismatch sites on rRFs corresponding to known modifications on 28 and 18S rRNA. Conclusion: We found abundant non-microRNA-small RNAs in BLCA tumor samples. Small RNAs, especially tRFs and rRFs, contain modifications that can be captured as mismatch by TGIRT sequencing. Both the modifications and the non-microRNA-small RNAs should be explored as a biomarker for BLCA detection or follow-up.

tRForest: a novel random forest-based algorithm for tRNA-derived fragment target prediction

tRNA fragments (tRFs) are small RNAs comparable to the size and function of miRNAs. tRFs are generally Dicer independent, are found associated with Ago, and can repress expression of genes post-transcriptionally. Given that this expands the repertoire of small RNAs capable of post-transcriptional gene expression, it is important to predict tRF targets with confidence. Some attempts have been made to predict tRF targets, but are limited in the scope of tRF classes used in prediction or limited in feature selection. We hypothesized that established miRNA target prediction features applied to tRFs through a random forest machine learning algorithm will immensely improve tRF target prediction. Using this approach, we show significant improvements in tRF target prediction for all classes of tRFs and validate our predictions in two independent cell lines. Finally, Gene Ontology analysis suggests that among the tRFs conserved between mice and humans, the predicted targets are enriched significantly in neuronal function, and we show this specifically for tRF-3009a. These improvements to tRF target prediction further our understanding of tRF function broadly across species and provide avenues for testing novel roles for tRFs in biology. We have created a publicly available website for the targets of tRFs predicted by tRForest.

Serum transfer RNA-derived fragment tRF-31-79MP9P9NH57SD acts as a novel diagnostic biomarker for non-small cell lung cancer

Background

tRNA‐derived fragments (tRFs) have been found to have a crucial function in the pathophysiology of cancers. However, the function of tRFs in non‐small cell lung cancer (NSCLC) is yet unknown. The goal of this study was to assess the tRF‐31‐79MP9P9NH57SD serum expression from NSCLC patients and to determine its diagnostic usefulness.

Methods

By using stem‐loop quantitative real‐time PCR, we were able to detect various tRF‐31‐79MP9P9NH57SD expressions in 96 NSCLC serum samples, 96 healthy controls, and 20 pairs of NSCLC serum samples pre‐ and post‐surgery (qRT‐PCR). After that, we analyzed its diagnostic effectiveness using the receiver operating characteristic (ROC) curve.

Results

Serum tRF‐31‐79MP9P9NH57SD expression was higher in NSCLC patients, and levels of tRF‐31‐79MP9P9NH57SD were linked to the clinical stage (p = 0.002) and the malignancy of lymph node (p = 0.012). In addition, after the procedure, the serum tRF‐31‐79MP9P9NH57SD expression in NSCLC patients dropped. With 48.96 percent sensitivity and 90.62 percent specificity, the area under ROC curve (AUC) was 0.733.

Conclusion

serum tRF‐31‐79MP9P9NH57SD possibly is a new and groundbreaking biomarker for the NSCLC.

Identification of tRNA-Related Fragments and Their Potential Regulatory Effects in Thyroid-Associated Ophthalmopathy

Recently, the potential role of tRNA-related fragments (tRFs) in ophthalmic diseases has been extensively researched. However, systematic studies on the potential regulatory effects of tRFs in thyroid-associated ophthalmopathy (TAO) are lacking. We used high-throughput sequencing techniques to measure expression levels of mRNAs and tRFs in patients with TAO, and the results were verified by real-time quantitative reverse transcription polymerase chain reaction (q-PCR). Next, the potential biological regulatory effect of differentially expressed tRFs was analyzed, and potential downstream target RNAs of differentially expressed tRFs were predicted to explore the potential role of tRFs as therapeutic targets and biomarkers of TAO. A total of 50 tRFs and 361 mRNAs were dysregulated in the TAO group, and tRF5-GluCTC, PMAIP1, HSD17B2 and ATF3 were verified to be significantly differentially expressed in TAO. Our research reveals that several associated pathways likely play a role in the pathogenesis of TAO. By targeting ATF3, HSD17B2 and PMAIP1, tRF5-GluCTC may play a potential role in regulating the orbital fibroblast adipogenic response and fibrotic hyperplasia in patients with TAO.

Extracellular vesicles-associated tRNA-derived fragments (tRFs): biogenesis, biological functions, and their role as potential biomarkers in human diseases

Traditionally, transfer RNAs (tRNAs) specifically decoded messenger RNA (mRNA) and participated in protein translation. tRNA-derived fragments (tRFs), also known as tRNA-derived small RNAs (tsRNAs), are generated by the specific cleavage of pre- and mature tRNAs and are a class of newly defined functional small non-coding RNAs (sncRNAs). Following the different cleavage positions of precursor or mature tRNA, tRFs are classified into seven types, 5′-tRNA half, 3′-tRNA half, tRF-1, 5′U-tRF, 3′-tRF, 5′-tRF, and i-tRF. It has been demonstrated that tRFs have a diverse range of biological functions in cellular processes, which include inhibiting protein translation, modulating stress response, regulating gene expression, and involvement in cell cycles and epigenetic inheritance. Emerging evidences have indicated that tRFs in extracellular vesicles (EVs) seem to act as regulatory molecules in various cellular processes and play essential roles in cell-to-cell communication. Furthermore, the dysregulation of EV-associated tRFs has been associated with the occurrence and progression of a variety of cancers and they can serve as novel potential biomarkers for cancer diagnosis. In this review, the biogenesis and classification of tRFs are summarized, and the biological functions of EV-associated tRFs and their roles as potential biomarkers in human diseases are discussed.

Nuclear-Mitochondrial Interactions

Mitochondria, the cell’s major energy producers, also act as signaling hubs, interacting with other organelles both directly and indirectly. Despite having its own circular genome, the majority of mitochondrial proteins are encoded by nuclear DNA. To respond to changes in cell physiology, the mitochondria must send signals to the nucleus, which can, in turn, upregulate gene expression to alter metabolism or initiate a stress response. This is known as retrograde signaling. A variety of stimuli and pathways fall under the retrograde signaling umbrella. Mitochondrial dysfunction has already been shown to have severe implications for human health. Disruption of retrograde signaling, whether directly associated with mitochondrial dysfunction or cellular environmental changes, may also contribute to pathological deficits. In this review, we discuss known signaling pathways between the mitochondria and the nucleus, examine the possibility of direct contacts, and identify pathological consequences of an altered relationship.

An emerging role of the 5' termini of mature tRNAs in human diseases: Current situation and prospects

The fundamental biological roles of a class of small noncoding RNAs (sncRNAs), derived from mature tRNAs or pre-tRNAs, in human diseases have received increasing attention in recent years. These ncRNAs are called tRNA-derived fragments (tRFs) or tRNA-derived small RNAs (tsRNAs). tRFs mainly include tRF-1, tRF-5, tRF-3 and tRNA halves (tiRNAs or tRHs), which are produced by enzyme-specific cleavage of tRNAs. Here, we classify tRF-5 and 5' tiRNAs into the same category: 5'-tRFs and review the biological functions and regulatory mechanisms of 5'-tRFs in cancer and other diseases (metabolic diseases, neurodegenerative diseases, pathological stress injury and virus infection) to provide a new theoretical basis for the diagnosis and treatment of diseases.

Screening and Comprehensive Analysis of Cancer-Associated tRNA-Derived Fragments

tRNA-derived fragments (tRFs) constitute a novel class of small non-coding RNA cleaved from tRNAs. In recent years, researches have shown the regulatory roles of a few tRFs in cancers, illuminating a new direction for tRF-centric cancer researches. Nonetheless, more specific screening of tRFs related to oncogenesis pathways, cancer progression stages and cancer prognosis is continuously demanded to reveal the landscape of the cancer-associated tRFs. In this work, by combining the clinical information recorded in The Cancer Genome Atlas (TCGA) and the tRF expression profiles curated by MINTbase v2.0, we systematically screened 1,516 cancer-associated tRFs (ca-tRFs) across seven cancer types. The ca-tRF set collectively combined the differentially expressed tRFs between cancer samples and control samples, the tRFs significantly correlated with tumor stage and the tRFs significantly correlated with patient survival. By incorporating our previous tRF-target dataset, we found the ca-tRFs tend to target cancer-associated genes and onco-pathways like ATF6-mediated unfolded protein response, angiogenesis, cell cycle process regulation, focal adhesion, PI3K-Akt signaling pathway, cellular senescence and FoxO signaling pathway across multiple cancer types. And cell composition analysis implies that the expressions of ca-tRFs are more likely to be correlated with T-cell infiltration. We also found the ca-tRF expression pattern is informative to prognosis, suggesting plausible tRF-based cancer subtypes. Together, our systematic analysis demonstrates the potentially extensive involvements of tRFs in cancers, and provides a reasonable list of cancer-associated tRFs for further investigations.

Expression Profiles of tRNA-Derived Small RNAs and Their Potential Roles in Primary Nasopharyngeal Carcinoma

Introduction: tRNA-derived small RNAs (tsRNAs), a class of small non-coding RNAs, are divided into two categories: tRNA-related fragments (tRFs) and tRNA halves (tiRNAs). Abnormal expression of tsRNAs has been found in diverse cancers, which indicates that further understanding of the function of tsRNAs will help identify new biomarkers and potential therapeutic targets. Until now, the underlying roles of tsRNAs in primary nasopharyngeal carcinoma (NPC) are still unknown.

Methods: tRF and tiRNA sequencing was performed on four pairs of NPC tissues and healthy controls. Thirty pairs of NPC samples were used for quantitative real-time polymerase chain reaction (qRT-PCR) verification, and the ROC analysis was used to evaluate the diagnostic efficiency initially. Target prediction and bioinformatics analysis of validated tRFs and tiRNAs were conducted to explore the mechanisms of tsRNAs in NPC’s pathogenesis.

Results: A total of 158 differentially expressed tRFs and tiRNAs were identified, of which 88 are upregulated and 70 are downregulated in NPC. Three validated tRFs in the results of qRT-PCR were consistent with the sequencing data: two upregulations (tRF-1:28-Val-CAC-2 and tRF-1:24-Ser-CGA-1-M3) and one downregulation (tRF-55:76-Arg-ACG-1-M2). The GO and KEGG pathway enrichment analysis showed that the potential target genes of validated tRFs are widely enriched in cancer pathways. The related modules may play an essential role in the pathogenesis of NPC.

Conclusions: The tsRNAs may become a novel class of biological diagnostic indicators and possible targets for NPC.

tRNAGlyGCC-Derived Internal Fragment (i-tRF-GlyGCC) in Ovarian Cancer Treatment Outcome and Progression

Epithelial ovarian cancer (EOC) remains a highly-lethal gynecological malignancy, characterized by frequent recurrence, chemotherapy resistance and poor 5-year survival. Identifying novel predictive molecular markers remains an overdue challenge in the disease's clinical management. Herein, in silico analysis of TCGA-OV highlighted the tRNA-derived internal fragment (i-tRF-GlyGCC) among the most abundant tRFs in ovarian tumors, while target prediction and gene ontology (GO) enrichment analysis predicted its implication in key biological processes. Thereafter, i-tRF-GlyGCC levels were quantified in a screening EOC (n = 98) and an institutionally-independent serous ovarian cancer (SOC) validation cohort (n = 100, OVCAD multicenter study). Disease progression and patient death were used as clinical endpoints for the survival analysis. Internal validation was performed by bootstrap analysis and the clinical net benefit was estimated by decision curve analysis. The analysis highlighted the significant association of i-tRF-GlyGCC with advanced FIGO stages, suboptimal debulking and most importantly, with early progression and poor overall survival of EOC patients. The OVCAD validation cohort corroborated the unfavorable predictive value of i-tRF-GlyGCC in EOC. Ultimately, evaluation of i-tRF-GlyGCC with the established/clinically used prognostic markers offered superior patient risk-stratification and enhanced clinical benefit in EOC prognosis. In conclusion, i-tRF-GlyGCC assessment could aid towards personalized prognosis and support precision medicine decisions in EOC.

tsRFun: a comprehensive platform for decoding human tsRNA expression, functions and prognostic value by high-throughput small RNA-Seq and CLIP-Seq data

tRNA-derived small RNA (tsRNA), a novel type of regulatory small noncoding RNA, plays an important role in physiological and pathological processes. However, the understanding of the functional mechanism of tsRNAs in cells and their role in the occurrence and development of diseases is limited. Here, we integrated multiomics data such as transcriptome, epitranscriptome, and targetome data, and developed novel computer tools to establish tsRFun, a comprehensive platform to facilitate tsRNA research (http://rna.sysu.edu.cn/tsRFun/ or http://biomed.nscc-gz.cn/DB/tsRFun/). tsRFun evaluated tsRNA expression profiles and the prognostic value of tsRNAs across 32 types of cancers, identified tsRNA target molecules utilizing high-throughput CLASH/CLEAR or CLIP sequencing data, and constructed the interaction networks among tsRNAs, microRNAs, and mRNAs. In addition to its data presentation capabilities, tsRFun offers multiple real-time online tools for tsRNA identification, target prediction, and functional enrichment analysis. In summary, tsRFun provides a valuable data resource and multiple analysis tools for tsRNA investigation.

tRNA Metabolism and Lung Cancer: Beyond Translation

Lung cancer, one of the most malignant tumors, has extremely high morbidity and mortality, posing a serious threat to global health. It is an urgent need to fully understand the pathogenesis of lung cancer and provide new ideas for its treatment. Interestingly, accumulating evidence has identified that transfer RNAs (tRNAs) and tRNA metabolism–associated enzymes not only participate in the protein translation but also play an important role in the occurrence and development of lung cancer. In this review, we summarize the different aspects of tRNA metabolism in lung cancer, such as tRNA transcription and mutation, tRNA molecules and derivatives, tRNA-modifying enzymes, and aminoacyl-tRNA synthetases (ARSs), aiming at a better understanding of the pathogenesis of lung cancer and providing new therapeutic strategies for it.

tRNA-derived small RNAs: novel regulators of cancer hallmarks and targets of clinical application

tRNAs are a group of conventional noncoding RNAs (ncRNAs) with critical roles in the biological synthesis of proteins. Recently, tRNA-derived small RNAs (tsRNAs) were found to have important biological functions in the development of human diseases including carcinomas, rather than just being considered pure degradation material. tsRNAs not only are abnormally expressed in the cancer tissues and serum of cancer patients, but also have been suggested to regulate various vital cancer hallmarks. On the other hand, the application of tsRNAs as biomarkers and therapeutic targets is promising. In this review, we focused on the basic characteristics of tsRNAs, and their biological functions known thus far, and explored the regulatory roles of tsRNAs in cancer hallmarks including proliferation, apoptosis, metastasis, tumor microenvironment, drug resistance, cancer stem cell phenotype, and cancer cell metabolism. In addition, we also discussed the research progress on the application of tsRNAs as tumor biomarkers and therapeutic targets.

The tRNA-derived fragment 5026a inhibits the proliferation of gastric cancer cells by regulating the PTEN/PI3K/AKT signaling pathway

Background

Recently, tRNA-derived fragments (tRFs) have been shown to serve important biological functions. However, the role of tRFs in gastric cancer has not been fully elucidated. This study aimed to identify the tumor suppressor role of tRF-5026a (tRF-18-79MP9P04) in gastric cancer.

Methods

Quantitative reverse transcription-polymerase chain reaction (qRT-PCR) was first used to detect tRF-5026a expression levels in gastric cancer tissues and patient plasma. Next, the relationship between tRF-5026a levels and clinicopathological features in gastric cancer patients was assessed. Cell lines with varying tRF-5026a levels were assessed by measuring tRF-5026a using qRT-PCR. After transfecting cell lines with a tRF-5026a mimic or inhibitor, cell proliferation, colony formation, migration, apoptosis, and cell cycle were evaluated. The expression levels of related proteins in the PTEN/PI3K/AKT pathway were also analyzed by Western blotting. Finally, the effect of tRF-5026a on tumor growth was tested using subcutaneous tumor models in nude mice.

Results

tRF-5026a was downregulated in gastric cancer patient tissues and plasma samples. tRF-5026a levels were closely related to tumor size, had a certain diagnostic value, and could be used to predict overall survival. tRF-5026a was also downregulated in gastric cancer cell lines. tRF-5026a inhibited the proliferation, migration, and cell cycle progression of gastric cancer cells by regulating the PTEN/PI3K/AKT signaling pathway. Animal experiments showed that upregulation of tRF-5026a effectively inhibited tumor growth.

Conclusions

tRF-5026a (tRF-18-79MP9P04) is a promising biomarker for gastric cancer diagnostics and has tumor suppressor effects mediated through the PTEN/PI3K/AKT signaling pathway.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13287-021-02497-1.

The transcriptional trajectories of pluripotency and differentiation comprise genes with antithetical architecture and repetitive-element content

Background

Extensive molecular differences exist between proliferative and differentiated cells. Here, we conduct a meta-analysis of publicly available transcriptomic datasets from preimplantation and differentiation stages examining the architectural properties and content of genes whose abundance changes significantly across developmental time points.

Results

Analysis of preimplantation embryos from human and mouse showed that short genes whose introns are enriched in Alu (human) and B (mouse) elements, respectively, have higher abundance in the blastocyst compared to the zygote. These highly expressed genes encode ribosomal proteins or metabolic enzymes. On the other hand, long genes whose introns are depleted in repetitive elements have lower abundance in the blastocyst and include genes from signaling pathways. Additionally, the sequences of the genes that are differentially expressed between the blastocyst and the zygote contain distinct collections of pyknon motifs that differ between up- and down-regulated genes. Further examination of the genes that participate in the stem cell-specific protein interaction network shows that their introns are short and enriched in Alu (human) and B (mouse) elements. As organogenesis progresses, in both human and mouse, we find that the primarily short and repeat-rich expressed genes make way for primarily longer, repeat-poor genes. With that in mind, we used a machine learning-based approach to identify gene signatures able to classify human adult tissues: we find that the most discriminatory genes comprising these signatures have long introns that are repeat-poor and include transcription factors and signaling-cascade genes. The introns of widely expressed genes across human tissues, on the other hand, are short and repeat-rich, and coincide with those with the highest expression at the blastocyst stage.

Conclusions

Protein-coding genes that are characteristic of each trajectory, i.e., proliferation/pluripotency or differentiation, exhibit antithetical biases in their intronic and exonic lengths and in their repetitive-element content. While the respective human and mouse gene signatures are functionally and evolutionarily conserved, their introns and exons are enriched or depleted in organism-specific repetitive elements. We posit that these organism-specific repetitive sequences found in exons and introns are used to effect the corresponding genes’ regulation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-020-00928-8.

Small Non-Coding-RNA in Gynecological Malignancies

Gynecologic malignancies, which include cancers of the cervix, ovary, uterus, vulva, vagina, and fallopian tube, are among the leading causes of female mortality worldwide, with the most prevalent being endometrial, ovarian, and cervical cancer. Gynecologic malignancies are complex, heterogeneous diseases, and despite extensive research efforts, the molecular mechanisms underlying their development and pathology remain largely unclear. Currently, mechanistic and therapeutic research in cancer is largely focused on protein targets that are encoded by about 1% of the human genome. Our current understanding of 99% of the genome, which includes noncoding RNA, is limited. The discovery of tens of thousands of noncoding RNAs (ncRNAs), possessing either structural or regulatory functions, has fundamentally altered our understanding of genetics, physiology, pathophysiology, and disease treatment as they relate to gynecologic malignancies. In recent years, it has become clear that ncRNAs are relatively stable, and can serve as biomarkers for cancer diagnosis and prognosis, as well as guide therapy choices. Here we discuss the role of small non-coding RNAs, i.e., microRNAs (miRs), P-Element induced wimpy testis interacting (PIWI) RNAs (piRNAs), and tRNA-derived small RNAs in gynecological malignancies, specifically focusing on ovarian, endometrial, and cervical cancer.

tRFTars: predicting the targets of tRNA-derived fragments

Background

tRNA-derived fragments (tRFs) are 14–40-nucleotide-long, small non-coding RNAs derived from specific tRNA cleavage events with key regulatory functions in many biological processes. Many studies have shown that tRFs are associated with Argonaute (AGO) complexes and inhibit gene expression in the same manner as miRNAs. However, there are currently no tools for accurately predicting tRF target genes.

Methods

We used tRF-mRNA pairs identified by crosslinking, ligation, and sequencing of hybrids (CLASH) and covalent ligation of endogenous AGO-bound RNAs (CLEAR)-CLIP to assess features that may participate in tRF targeting, including the sequence context of each site and tRF-mRNA interactions. We applied genetic algorithm (GA) to select key features and support vector machine (SVM) to construct tRF prediction models.

Results

We first identified features that globally influenced tRF targeting. Among these features, the most significant were the minimum free folding energy (MFE), position 8 match, number of bases paired in the tRF-mRNA duplex, and length of the tRF, which were consistent with previous findings. Our constructed model yielded an area under the receiver operating characteristic (ROC) curve (AUC) = 0.980 (0.977–0.983) in the training process and an AUC = 0.847 (0.83–0.861) in the test process. The model was applied to all the sites with perfect Watson–Crick complementarity to the seed in the 3′ untranslated region (3′-UTR) of the human genome. Seven of nine target/nontarget genes of tRFs confirmed by reporter assay were predicted. We also validated the predictions via quantitative real-time PCR (qRT-PCR). Thirteen potential target genes from the top of the predictions were significantly down-regulated at the mRNA levels by overexpression of the tRFs (tRF-3001a, tRF-3003a or tRF-3009a).

Conclusions

Predictions can be obtained online, tRFTars, freely available at http://trftars.cmuzhenninglab.org:3838/tar/, which is the first tool to predict targets of tRFs in humans with a user-friendly interface.

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.

tRFtarget: a database for transfer RNA-derived fragment targets

Transfer RNA-derived fragments (tRFs) are a new class of small non-coding RNAs and play important roles in biological and physiological processes. Prediction of tRF target genes and binding sites is crucial in understanding the biological functions of tRFs in the molecular mechanisms of human diseases. We developed a publicly accessible web-based database, tRFtarget (http://trftarget.net), for tRF target prediction. It contains the computationally predicted interactions between tRFs and mRNA transcripts using the two state-of-the-art prediction tools RNAhybrid and IntaRNA, including location of the binding sites on the target, the binding region, and free energy of the binding stability with graphic illustration. tRFtarget covers 936 tRFs and 135 thousand predicted targets in eight species. It allows researchers to search either target genes by tRF IDs or tRFs by gene symbols/transcript names. We also integrated the manually curated experimental evidence of the predicted interactions into the database. Furthermore, we provided a convenient link to the DAVID® web server to perform downstream functional pathway analysis and gene ontology annotation on the predicted target genes. This database provides useful information for the scientific community to experimentally validate tRF target genes and facilitate the investigation of the molecular functions and mechanisms of tRFs.

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

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

Roles of mitochondria in the hallmarks of metastasis

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

Noncanonical Roles of tRNAs: tRNA Fragments and Beyond

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

Transfer-RNA-Derived Fragments Are Potential Prognostic Factors in Patients with Squamous Cell Carcinoma of the Head and Neck

Transfer-RNA-derived fragments (tRFs) are a class of small non-coding RNAs that are functionally different from their parental transfer RNAs (tRNAs). tRFs can regulate gene expression by several mechanisms, and are involved in a variety of pathological processes. Here, we aimed at understanding the composition and abundance of tRFs in squamous cell carcinoma of the head and neck (SCCHN), and evaluated the potential of tRFs as prognostic markers in this cancer type. We obtained tRF expression data from The Cancer Genome Atlas (TCGA) HNSC cohort (523 patients) using MINTbase v2.0, and correlated to available TCGA clinical data. RNA-binding proteins were predicted according to the calculated Position Weight Matrix (PWM) score from the RNA-Binding Protein DataBase (RBPDB). A total of 10,158 tRFs were retrieved and a high diversity in expression levels was seen. Fifteen tRFs were found to be significantly associated with overall survival (Kaplan-Meier survival analysis, log rank test p-value < 0.01). The top prognostic marker, tRF-20-S998LO9D (p < 0.001), was further measured in tumor and tumor-free samples from 16 patients with squamous cell carcinoma of the oral tongue and 12 healthy controls, and was significantly upregulated in tumor compared to matched tumor-free tongue (p < 0.001). Results suggest that tRFs are useful prognostic markers in SCCHN.

On the expanding roles of tRNA fragments in modulating cell behavior

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

Contribution of miRNAs, tRNAs and tRFs to Aberrant Signaling and Translation Deregulation in Lung Cancer

Simple Summary

The profiles of miRNAs, tRNA-derived fragments and tRNAs from lung cancer biopsy specimens indicate involvement of gene networks that modulate signaling and translation initiation. The current study highlights the important role of several regulatory small non-coding RNAs in aberrant signaling and translation deregulation in lung cancer.

Abstract

Transcriptomics profiles of miRNAs, tRNAs or tRFs are used as biomarkers, after separate examination of several cancer cell lines, blood samples or biopsies. However, the possible contribution of all three profiles on oncogenic signaling and translation as a net regulatory effect, is under investigation. The present analysis of miRNAs and tRFs from lung cancer biopsies indicated putative targets, which belong to gene networks involved in cell proliferation, transcription and translation regulation. In addition, we observed differential expression of specific tRNAs along with several tRNA-related genes with possible involvement in carcinogenesis. Transfection of lung adenocarcinoma cells with two identified tRFs and subsequent NGS analysis indicated gene targets that mediate signaling and translation regulation. Broader analysis of all major signaling and translation factors in several biopsy specimens revealed a crosstalk between the PI3K/AKT and MAPK pathways and downstream activation of eIF4E and eEF2. Subsequent polysome profile analysis and 48S pre-initiation reconstitution experiments showed increased global translation rates and indicated that aberrant expression patterns of translation initiation factors could contribute to elevated protein synthesis. Overall, our results outline the modulatory effects that possibly correlate the expression of important regulatory non-coding RNAs with aberrant signaling and translation deregulation in lung cancer.

Global identification and characterization of tRNA-derived RNA fragment landscapes across human cancers

Transfer RNA-derived RNA fragments (tRFs) are a class of small non-coding RNAs that are abundant in many organisms, but their role in cancer has not been fully explored. Here, we report a functional genomic landscape of tRFs in 8118 specimens across 15 cancer types from The Cancer Genome Atlas. These tRFs exhibited characteristics of widespread expression, high sequence conservation, cytoplasmic localization, specific patterns of tRNA cleavage and conserved cleavage in tissues. A cross-tumor analysis revealed significant commonality among tRF expression subtypes from distinct tissues of origins, characterized by upregulation of a group of tRFs with similar size and activation of cancer-associated signaling. One of the largest superclusters was composed of 22 nt 3'-tRFs upregulated in 13 cancer types, all of which share the activation of Ras/MAPK, RTK and TSC/mTOR signaling. tRF-based subgrouping provided clinically relevant stratifications and significantly improved outcome prediction by incorporating clinical variables. Additionally, we discovered 11 cancer driver tRFs using an effective approach for accurately exploring cross-tumor and platform trends. As a proof of concept, we performed comprehensive functional assays on a non-microRNA driver tRF, 5'-IleAAT-8-1-L20, and validated its oncogenic roles in lung cancer in vitro and in vivo. Our study also provides a valuable tRF resource for identifying diagnostic and prognostic biomarkers, developing cancer therapy and studying cancer pathogenesis.