The floral development of the allotetraploid Coffea arabica L. correlates with a small RNA dynamic reprogramming.

Noncoding and coding RNAs are key regulators of plant growth, development, and stress responses. To investigate the types of transcripts accumulated during the vegetative to reproductive transition and floral development in the Coffea arabica L., we sequenced small RNA libraries from eight developmental stages, up to anthesis. We combined these data with messenger RNA and PARE sequencing of two important development stages that marks the transition of an apparent latent to a rapid growth stage. In addition, we took advantage of multiple in silico tools to characterize genomic loci producing small RNAs such as phasiRNAs, miRNAs, and tRFs. Our differential and co-expression analysis showed that some types of small RNAs such as tRNAs, snoRNAs, snRNAs, and phasiRNAs preferentially accumulate in a stage-specific manner. Members of the miR482/miR2118 superfamily and their 21-nucleotide phasiRNAs originating from resistance genes show a robust co-expression pattern that is maintained across all the evaluated developmental stages. Finally, the majority of miRNAs accumulate in a family stage-specific manner, related to modulated hormonal responses and transcription factor expression.

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.

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.

miRPreM and tiRPreM: Improved methodologies for the prediction of miRNAs and tRNA-induced small non-coding RNAs for model and non-model organisms.

In recent years, microRNAs (miRNAs) and tRNA-derived RNA fragments (tRFs) have been reported extensively following different approaches of identification and analysis. Comprehensively analyzing the present approaches to overcome the existing variations, we developed a benchmarking methodology each for the identification of miRNAs and tRFs, termed as miRNA Prediction Methodology (miRPreM) and tRNA-induced small non-coding RNA Prediction Methodology (tiRPreM), respectively. We emphasized the use of respective genome of organism under study for mapping reads, sample data with at least two biological replicates, normalized read count support and novel miRNA prediction by two standard tools with multiple runs. The performance of these methodologies was evaluated by using Oryza coarctata, a wild rice species as a case study for model and non-model organisms. With organism-specific reference genome approach, 98 miRNAs and 60 tRFs were exclusively found. We observed high accuracy (13 out of 15) when tested these genome-specific miRNAs in support of analyzing the data with respective organism. Such a strong impact of miRPreM, we have predicted more than double number of miRNAs (186) as compared with the traditional approaches (79) and with tiRPreM, we have predicted all known classes of tRFs within the same small RNA data. Moreover, the methodologies presented here are in standard form in order to extend its applicability to different organisms rather than restricting to plants. Hence, miRPreM and tiRPreM can fulfill the need of a comprehensive methodology for miRNA prediction and tRF identification, respectively, for model and non-model organisms.

The function and therapeutic potential of transfer RNA-derived small RNAs in cardiovascular diseases: A review.

Transfer RNA-derived small RNAs (tsRNAs) are a class of small non-coding RNA (sncRNA) molecules derived from tRNA, including tRNA derived fragments (tRFs) and tRNA halfs (tiRNAs). tsRNAs can affect cell functions by participating in gene expression regulation, translation regulation, intercellular signal transduction, and immune response. They have been shown to play an important role in various human diseases, including cardiovascular diseases (CVDs). Targeted regulation of tsRNAs expression can affect the progression of CVDs. The tsRNAs induced by pathological conditions can be detected when released into the extracellular, giving them enormous potential as disease biomarkers. Here, we review the biogenesis, degradation process and related functional mechanisms of tsRNAs, and discuss the research progress and application prospects of tsRNAs in different CVDs, to provide a new perspective on the treatment of CVDs.

A versatile and efficient method for detecting tRNA-derived fragments.

Recently, it has been discovered surprisingly that tRNA can be cleaved into specific small fragments under certain conditions. Most importantly, these tRNA-derived fragments (tRFs) participate in the regulation of gene expression, playing pivotal roles in various physiological and pathological processes and thus attracting widespread attention. Detecting tRF expression in tissues and cells often involves using tRF-specific stem-loop primers for reverse transcription. However, the high specificity offered by this method limits it to transcribing only one specific tRF sequence per reaction, necessitating separate reverse transcription and qPCR steps for multiple tRFs, leading to substantially increased time and resource consumption. This becomes especially challenging in precious samples with limited RNA availability. To address these issues, there is an urgent need for a universal and cost-effective tRF identification method. This study introduces a versatile tRF detection approach based on the uniform polyadenylation of all tRFs, allowing reverse transcription with a universal oligo(dT) primer. This method enables simultaneous reverse transcription of all target tRFs in one reaction, greatly facilitating subsequent qPCR analysis. Furthermore, it demonstrates exceptional sensitivity and specificity, offering significant value in tRF-related research.

Gene expression differences consistent with water loss reduction underlie desiccation tolerance of natural Drosophila populations.

BACKGROUND: Climate change is one of the main factors shaping the distribution and biodiversity of organisms, among others by greatly altering water availability, thus exposing species and ecosystems to harsh desiccation conditions. However, most of the studies so far have focused on the effects of increased temperature. Integrating transcriptomics and physiology is key to advancing our knowledge on how species cope with desiccation stress, and these studies are still best accomplished in model organisms. RESULTS: Here, we characterized the natural variation of European D. melanogaster populations across climate zones and found that strains from arid regions were similar or more tolerant to desiccation compared with strains from temperate regions. Tolerant and sensitive strains differed not only in their transcriptomic response to stress but also in their basal expression levels. We further showed that gene expression changes in tolerant strains correlated with their physiological response to desiccation stress and with their cuticular hydrocarbon composition, and functionally validated three of the candidate genes identified. Transposable elements, which are known to influence stress response across organisms, were not found to be enriched nearby differentially expressed genes. Finally, we identified several tRNA-derived small RNA fragments that differentially targeted genes in response to desiccation stress. CONCLUSIONS: Overall, our results showed that basal gene expression differences across individuals should be analyzed if we are to understand the genetic basis of differential stress survival. Moreover, tRNA-derived small RNA fragments appear to be relevant across stress responses and allow for the identification of stress-response genes not detected at the transcriptional level.

Semen sEV tRF-Based Models Increase Non-Invasive Prediction Accuracy of Clinically Significant Prostate Cancer among Patients with Moderately Altered PSA Levels.

PSA screening has led to an over-diagnosis of prostate cancer (PCa) and unnecessary biopsies of benign conditions due to its low cancer specificity. Consequently, more accurate, preferentially non-invasive, tests are needed. We aim to evaluate the potential of semen sEV (small extracellular vesicles) tsRNAs (tRNA-derived small RNAs) as PCa indicators. Initially, following a literature review in the OncotRF database and high-throughput small RNA-sequencing studies in PCa tissue together with the sncRNA profile in semen sEVs, we selected four candidate 5'tRF tsRNAs for validation as PCa biomarkers. RT-qPCR analysis in semen sEVs from men with moderately elevated serum PSA levels successfully shows that the differential expression of the four tRFs between PCa and healthy control groups can be detected in a non-invasive manner. The combined model incorporating PSA and specific tRFs (5'-tRNA-Glu-TTC-9-1_L30 and 5'-tRNA-Val-CAC-3-1_L30) achieved high predictive accuracy in identifying samples with a Gleason score ≥ 7 and staging disease beyond IIA, supporting that the 5'tRF fingerprint in semen sEV can improve the PSA predictive value to discriminate between malignant and indolent prostate conditions. The in silico study allowed us to map target genes for the four 5'tRFs possibly involved in PCa. Our findings highlight the synergistic use of multiple biomarkers as an efficient approach to improve PCa screening and prognosis.

Transfer RNA-derived small RNAs in tumor microenvironment.

Transfer RNAs (tRNAs) are a class of non-coding RNAs responsible for amino acid translocation during protein synthesis and are ubiquitously found in organisms. With certain modifications and under specific conditions, tRNAs can be sheared and fragmented into small non-coding RNAs, also known as tRNA-derived small RNAs (tDRs). With the development of high-throughput sequencing technologies and bioinformatic strategies, more and more tDRs have been identified and their functions in organisms have been characterized. tRNA and it derived tDRs, have been shown to be essential not only for transcription and translation, but also for regulating cell proliferation, apoptosis, metastasis, and immunity. Aberrant expression of tDRs is associated with a wide range of human diseases, especially with tumorigenesis and tumor progression. The tumor microenvironment (TME) is a complex ecosystem consisting of various cellular and cell-free components that are mutually compatible with the tumor. It has been shown that tDRs regulate the TME by regulating cancer stem cells, immunity, energy metabolism, epithelial mesenchymal transition, and extracellular matrix remodeling, playing a pro-tumor or tumor suppressor role. In this review, the biogenesis, classification, and function of tDRs, as well as their effects on the TME and the clinical application prospects will be summarized and discussed based on up to date available knowledge.

Comprehensive annotation and characterization of planarian tRNA and tRNA-derived fragments (tRFs).

tRNA-derived fragments (tRFs) have recently gained a lot of scientific interest due to their diverse regulatory roles in several cellular processes. However, their function in dynamic biological processes such as development and regeneration remains unexplored. Here, we show that tRFs are dynamically expressed during planarian regeneration, suggesting a possible role for these small RNAs in the regulation of regeneration. In order to characterize planarian tRFs, we first annotated 457 tRNAs in S. mediterranea combining two tRNA prediction algorithms. Annotation of tRNAs facilitated the identification of three main species of tRFs in planarians-the shorter tRF-5s and itRFs, and the abundantly expressed 5'-tsRNAs. Spatial profiling of tRFs in sequential transverse sections of planarians revealed diverse expression patterns of these small RNAs, including those that are enriched in the head and pharyngeal regions. Expression analysis of these tRF species revealed dynamic expression of these small RNAs over the course of regeneration suggesting an important role in planarian anterior and posterior regeneration. Finally, we show that 5'-tsRNA in planaria interact with all three SMEDWI proteins and an involvement of AGO1 in the processing of itRFs. In summary, our findings implicate a novel role for tRFs in planarian regeneration, highlighting their importance in regulating complex systemic processes. Our study adds to the catalog of posttranscriptional regulatory systems in planaria, providing valuable insights on the biogenesis and the function of tRFs in neoblasts and planarian regeneration.

Elucidation of transfer RNAs as stress regulating agents and the experimental strategies to conceive the functional role of tRNA-derived fragments in plants.

In plants, the transfer RNAs (tRNAs) exhibit their profound influence in orchestrating diverse physiological activities like cell growth, development, and response to several surrounding stimuli. The tRNAs, which were known to restrict their function solely in deciphering the codons, are now emerging as frontline defenders in stress biology. The plants that are constantly confronted with a huge panoply of stresses rely on tRNA-mediated stress regulation by altering the tRNA abundance, curbing the transport of tRNAs, fragmenting the mature tRNAs during stress. Among them, the studies on the generation of transfer RNA-derived fragments (tRFs) and their biological implication in stress response have attained huge interest. In plants, the tRFs hold stable expression patterns and regulate biological functions under diverse environmental conditions. In this review, we discuss the fate of plant tRNAs upon stress and thereafter how the tRFs are metamorphosed into sharp ammunition to wrestle with stress. We also address the various methods developed to date for uncovering the role of tRFs and their function in plants.

Sex-specific declines in cholinergic-targeting tRNA fragments in the nucleus accumbens in Alzheimer's disease.

INTRODUCTION: Females with Alzheimer's disease (AD) suffer accelerated dementia and loss of cholinergic neurons compared to males, but the underlying mechanisms are unknown. Seeking causal contributors to both these phenomena, we pursued changes in transfer RNS (tRNA) fragments (tRFs) targeting cholinergic transcripts (CholinotRFs). METHODS: We analyzed small RNA-sequencing (RNA-Seq) data from the nucleus accumbens (NAc) brain region which is enriched in cholinergic neurons, compared to hypothalamic or cortical tissues from AD brains; and explored small RNA expression in neuronal cell lines undergoing cholinergic differentiation. RESULTS: NAc CholinotRFs of mitochondrial genome origin showed reduced levels that correlated with elevations in their predicted cholinergic-associated mRNA targets. Single-cell RNA seq from AD temporal cortices showed altered sex-specific levels of cholinergic transcripts in diverse cell types; inversely, human-originated neuroblastoma cells under cholinergic differentiation presented sex-specific CholinotRF elevations. DISCUSSION: Our findings support CholinotRFs contributions to cholinergic regulation, predicting their involvement in AD sex-specific cholinergic loss and dementia.

Characterization of tRNA-Derived Fragments in Lung Squamous Cell Carcinoma with Respect to Tobacco Smoke.

Lung squamous cell carcinoma (LUSC) is a highly heterogeneous cancer that is influenced by etiological agents such as tobacco smoke. Accordingly, transfer RNA-derived fragments (tRFs) are implicated in both cancer onset and development and demonstrate the potential to act as targets for cancer treatments and therapies. Therefore, we aimed to characterize tRF expression with respect to LUSC pathogenesis and clinical outcomes. Specifically, we analyzed the effect of tobacco smoke on tRF expression. In order to do so, we extracted tRF read counts from MINTbase v2.0 for 425 primary tumor samples and 36 adjacent normal samples. We analyzed the data in three primary cohorts: (1) all primary tumor samples (425 samples), (2) smoking-induced LUSC primary tumor samples (134 samples), and (3) non-smoking-induced LUSC primary tumor samples (18 samples). Differential expression analysis was performed to examine tRF expression in each of the three cohorts. tRF expression was correlated to clinical variables and patient survival outcomes. We identified unique tRFs in primary tumor samples, smoking-induced LUSC primary tumor samples, and non-smoking-induced LUSC primary tumor samples. In addition, many of these tRFs demonstrated correlations to worse patient survival outcomes. Notably, tRFs in the smoking-induced LUSC and non-smoking-induced LUSC primary tumor cohorts were significantly correlated to clinical variables pertaining to cancer stage and treatment efficacy. We hope that our results will be used to better inform future LUSC diagnostic and therapeutic modalities.

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.

The Fate and Functionality of Alien tRNA Fragments in Culturing Medium and Cells of Escherichia coli.

Numerous observations have supported the idea that various types of noncoding RNAs, including tRNA fragments (tRFs), are involved in communications between the host and its microbial community. The possibility of using their signaling function has stimulated the study of secreted RNAs, potentially involved in the interspecies interaction of bacteria. This work aimed at identifying such RNAs and characterizing their maturation during transport. We applied an approach that allowed us to detect oligoribonucleotides secreted by Prevotella copri (Segatella copri) or Rhodospirillum rubrum inside Escherichia coli cells. Four tRFs imported by E. coli cells co-cultured with these bacteria were obtained via chemical synthesis, and all of them affected the growth of E. coli. Their successive modifications in the culture medium and recipient cells were studied by high-throughput cDNA sequencing. Instead of the expected accidental exonucleolysis, in the milieu, we observed nonrandom cleavage by endonucleases continued in recipient cells. We also found intramolecular rearrangements of synthetic oligonucleotides, which may be considered traces of intermediate RNA circular isomerization. Using custom software, we estimated the frequency of such events in transcriptomes and secretomes of E. coli and observed surprising reproducibility in positions of such rare events, assuming the functionality of ring isoforms or their permuted derivatives in bacteria.

Spatiotemporal expression profile of novel and known small RNAs throughout rice plant development focussing on seed tissues.

BACKGROUND: Small RNAs (sRNAs) regulate numerous plant processes directly related to yield, such as disease resistance and plant growth. To exploit this yield-regulating potential of sRNAs, the sRNA profile of one of the world's most important staple crops - rice - was investigated throughout plant development using next-generation sequencing. RESULTS: Root and leaves were investigated at both the vegetative and generative phase, and early-life sRNA expression was characterized in the embryo and endosperm. This led to the identification of 49,505 novel sRNAs and 5581 tRNA-derived sRNAs (tsRNAs). In all tissues, 24 nt small interfering RNAs (siRNAs) were highly expressed and associated with euchromatic, but not heterochromatic transposable elements. Twenty-one nt siRNAs deriving from genic regions in the endosperm were exceptionally highly expressed, mimicking previously reported expression levels of 24 nt siRNAs in younger endosperm samples. In rice embryos, sRNA content was highly diverse while tsRNAs were underrepresented, possibly due to snoRNA activity. Publicly available mRNA expression and DNA methylation profiles were used to identify putative siRNA targets in embryo and endosperm. These include multiple genes related to the plant hormones gibberellic acid and ethylene, and to seed phytoalexin and iron content. CONCLUSIONS: This work introduces multiple sRNAs as potential regulators of rice yield and quality, identifying them as possible targets for the continuous search to optimize rice production.

Plasma tRNA-derived small RNAs signature as a predictive and prognostic biomarker in lung adenocarcinoma.

BACKGROUND: The prevalence of lung adenocarcinoma (LUAD) has increased, thus novel biomarkers for its early diagnosis is becoming more important than ever. tRNA-derived small RNA (tsRNA) is a new class of non-coding RNA which has important regulatory roles in cancer biology. This study was designed to identify novel predictive and prognostic tsRNA biomarkers. METHODS: tsRNAs were identified and performed differential expression analysis from 10 plasma samples (6 LUAD and 4 normal, SRP266333) and 96 tissue samples (48 LUAD and 48 normal, SRP133217). Then a tsRNA-mRNA regulatory network was constructed to find hub tsRNAs. Functional enrichment analysis was performed to infer the potential pathways associated with tsRNAs. Afterwards, a Support Vector Machine (SVM) algorithm was used to explore the potential biomarkers for diagnosing LUAD. Lastly, the function of tRF-21-RK9P4P9L0 was explored in A549 and H1299 cell lines. RESULTS: A significant difference of read distribution was observed between normal people and LUAD patients whether in plasma or tissue. A tsRNA-mRNA regulatory network consisting of 155 DEtsRNAs (differential expression tsRNAs) and 406 DEmRNAs (differential expression mRNAs) was established. Three tsRNAs (tRF-16-L85J3KE, tRF-21-RK9P4P9L0 and tRF-16-PSQP4PE) were identified as hub genes with degree > 100. We found Co-DEmRNAs (intersection of DEtsRNAs target mRNAs and differentially expressed mRNAs in LUAD) were engaged in a number of cancer pathways. The AUC of the three hub tsRNAs' expression for diagnosing LUAD reached 0.92. Furthermore, the qPCR validation of the three hub tsRNAs in 37 paired normal and LUAD tissues was consistent with the RNA-Seq results. In addition, tRF-21-RK9P4P9L0 was negatively associated with LUAD prognosis. Inhibition of tRF-21-RK9P4P9L0 expression reduced the proliferation, migration and invasion ability of A549 and H1299 cell lines. CONCLUSION: These findings will help us further understand the molecular mechanisms of LUAD and contribute to novel diagnostic biomarkers and therapeutic target discovery.

tRNA-derived fragments (tRFs): establishing their turf in post-transcriptional gene regulation.

Transfer RNA (tRNA)-derived fragments (tRFs) are an emerging class of conserved small non-coding RNAs that play important roles in post-transcriptional gene regulation. High-throughput sequencing of multiple biological samples have identified heterogeneous species of tRFs with distinct functionalities. These small RNAs have garnered a lot of scientific attention due to their ubiquitous expression and versatility in regulating various biological processes. In this review, we highlight our current understanding of tRF biogenesis and their regulatory functions. We summarize the diverse modes of biogenesis through which tRFs are generated and discuss the mechanism through which different tRF species regulate gene expression and the biological implications. Finally, we conceptualize research areas that require focus to strengthen our understanding of the biogenesis and function of tRFs.

RNA-sequencing reveals the expression profiles of tsRNAs and their potential carcinogenic role in cholangiocarcinoma.

BACKGROUND: Recently, the incidence of cholangiocarcinoma (CCA) has gradually increased. As CCA has a poor prognosis, the ideal survival rate is scarce for patients. The abnormal expressed tsRNAs may regulate the progression of a variety of tumors, and tsRNAs is expected to become a new diagnostic biomarker of cancer. However, the expression of tsRNAs is obscure and should be elucidated in CCA. METHODS: High-throughput RNA sequencing technology (RNA-seq) was utilized to determine the overall expression profiles of tsRNAs in three pairs CCA and adjacent normal tissues and to screen the tsRNAs that were differentially expressed. The target genes of dysregulated tsRNAs were predicted and the biological effects and potential signaling pathways of these target genes were explored by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses. Quantitative real-time polymerase chain reaction (qRT-PCR) was used to validate 11 differentially expressed tRFs with 12 pairs CCA and adjacent normal tissues. RESULTS: High-throughput RNA-seq totally demonstrated 535 dysregulated tsRNAs, of which 241 tsRNAs were upregulated, such as tRF-21-YLKZKWE5D，tRF-16-9NF5W8B，tRF-27-78YLKZKWE52，tRF-19-RLXN48KP，tRF-33-IK9NJ4S2I7L7DV，tRF-19-F8DHXYIV, and 294 tsRNAs were downregulated (tRF-20-739P8WQ0, tRF-34-JJ6RRNLIK898HR, tRF-17-VL8RPY5, tRF-23-YP9LON4VDP, tRF-39-EH623K76IR3DR2I2, tRF-17-18YKISM, tRF-19-Q1Q89PJZ, etc.) in CCA compared with adjacent normal tissues (|log2 [fold change] | ≥ 1 and p value <0.05). GO and KEGG enrichment analyses indicated that the target genes of dysregulated tRFs (tRF-34-JJ6RRNLIK898HR, tRF-38-0668K87SERM492V, and tRF-39-0668K87SERM492E2) were mainly enriched in the Notch signaling pathway, Hippo signaling pathway, cAMP signaling pathway and in growth hormone synthesis, secretion and action, etc. qRT-PCR result showed that tRF-34-JJ6RRNLIK898HR/tRF-38-0668K87SERM492V/tRF-39-0668K87SERM492E2 was downregulated (p = 0.021), and tRF-20-LE2WMK81 was upregulated in CCA (p = 0.033). CONCLUSION: Differentially expressed tRFs in CCA are enriched in many pathways associated with neoplasms, which may impact the tumor progression and have potential to be diagnostic biomarkers and therapeutic targets of CCA.

Dysregulation of different classes of tRNA fragments in chronic lymphocytic leukemia.

Chronic lymphocytic leukemia (CLL) is the most common human leukemia, and dysregulation of tRNA-derived short noncoding RNA (tsRNA) (tRF-1) expression is an accompanying event in the development of this disease. tsRNAs are fragments originating from the 3' end of tRNA precursors and do not contain mature tRNA sequences. In contrast to tsRNAs, mature tRFs (tRF-3s, tRF-5s, and internal tRFs) are produced from mature tRNA sequences and are redundant fragments. We investigated tsRNA expression in CLL and determined tsRNA signatures in indolent CLL and aggressive CLL vs. normal B cells. We noticed that both ts-43 and ts-44 are derived from distinct genes of pre-tRNAHis, and are down-regulated in CLL 3- to 5-fold vs. normal B cells. Thus, we investigated expression levels of tRF-5 fragments from tRNAHis in CLL samples and healthy controls, and determined that such fragments are down-regulated by 5-fold in CLLs vs. normal controls. Given these results, we investigated the expression of all mature tRFs in CLLs vs. normal controls. We found a drastic dysregulation of the expression of mature tRFs in CLL. In aggressive CLL, for the top 15 up-regulated fragments, linear fold change varied from 2,053- to 622-fold. For the top 15 down-regulated fragments in CLL, linear fold change varied from 314- to 52-fold. In addition, 964 mature tRFs were up-regulated at least 2-fold in CLL, while 701 fragments were down-regulated at least 2-fold. Similar results were obtained for indolent CLL. Our results suggest that mature tRFs may have oncogenic and/or tumor suppressor function in CLL.