Exosomal tRNA-derived small RNA as a promising biomarker for cancer diagnosis

1tRNA-derived small RNAs (tsRNAs), usually 18~40 nucleotides (nt) in length, are novel small non-coding RNAs (sncRNAs) generated from precursor or mature tRNAs. tsRNAs can be grouped into three distinctive classes, including precursor tRNA-derived small RNAs with the characteristic of poly U residues at 3′ terminus (3′ U tRFs), and mature tRNA-derived fragments (tRFs) as well as halves (tRHs). And tRFs are further classified into 3 sub-groups: 5′ tRF, 3′ tRF and inter tRF. Biogenesis of different tsRNAs is regulated by distinct mechanisms. For example, 3′ U tRFs are generated by RNase Z during tRNA maturation, while tRHs are derived from angiogenin cleavage within anti-codon loop of mature tRNAs. Increasing evidences indicate that tsRNAs expression is spatially and temporally controlled under physiological conditions, thus playing an important role in many biological processes. For instance, sperm tsRNAs were discovered to be a paternal epigenetic factor mediating intergenerational inheritance of diet-induced metabolic disorders. Recently, Croce’s group found that tsRNA expression is dysregulated in lung cancer and chronic lymphocytic leukemia, suggesting that tsRNAs participate in tumor initiation and development.

2Exosomes are membrane-bound vehicles with 30–100 nm in diameter secreted by most cell types and present in various types of body fluids, including plasma/serum, urine and saliva. The exosome cargo consists of specific lipids, proteins and RNA molecules. Among them, microRNA (miRNA), circular RNA and long non-coding RNA have shown great potential as diagnostic or prognostic biomarkers. However, exosomal tsRNAs as diagnostic biomarkers have not been reported. In this study, we demonstrate the presence and expression pattern of tsRNAs in exosomes from cell culture media and patients’ plasma, highlighting their potential for cancer diagnosis.

3To comprehensively profile exosomal tsRNAs, we isolated exosomes from cultured medium of SK-Hep1 liver cancer cells and prepared small RNA-seq library for high-throughput sequencing (Fig. 1a). The prepared exosomes were confirmed by transmission electron microscopy (TEM), Western blot analysis (exosomal protein marker CD63 and non-exosomal marker Calnexin), and the particle size (30–100 nm) (Fig. 1b, Additional file 1: Figure S1, Additional file 2). The raw sequencing data were first cleaned through removing both 5′ and 3′ adaptors, and only high quality reads with 16–40 nt insertion were mapped to human genome and annotated (Additional file 2). We have performed exosomal small RNA sequencing twice with high reproducibility (R = 0.988; Fig. 1c). Our results showed that 5% of exosomal small RNAs were generated from tRNAs, suggesting that tsRNAs could be incorporated into exosomes and secreted into extracellular environment (Fig. 1d). Because tsRNAs generated from precursor tRNA have low levels in exosomes, so we focus on mature tRNA-derived tsRNAs that were further classified into tRNA-5 (5-termial of mature tRNA, 5′ tRF and 5′ tRH), tRF-i (internal of mature tRNA, i′ tRF) and tRNA-3 (3-termial of mature tRNA, 3′ tRF and 3′ tRH) according to their position (Additional file 1: Figure S2). As shown in Fig. 1e, tRNA-5 is the most abundant tsRNA in exosomes, accounting for 90%, whereas tRNA-3 and tRNA-i are 9 and 1%, respectively. We also analyzed the length distribution of tsRNAs and found that distributions of each type of tsRNA exhibited significant bias. Specifically, large majority of tRNA-5s were 32~33 nt in length (Fig. 1g), while tRNA-3s were mainly 16–18 nt in length and tRNA-i showed two peaks (16–17 nt and 28–32 nt) (Additional file 1: Figure S3).

4We analyzed the origin of these exosomal tsRNAs and found that 97% of tRNA-5s came from four tRNAs (Gly-, Glu-, Lys- and Val-tRNA) (Fig. 1f). Moreover, tRNA-5s from Gly-tRNA were mainly 32 nt in length, while the majority of tRNA-5s from both Lys- and Val-tRNA were 33 nt in length (Fig. 1h), indicating that tsRNAs are not random degradation products, but generated by unknown mechanism. In addition, 52% of tRNA-5 and 32% of RNA-is were mainly derived from Gly-tRNA (Fig. 1f, Additional file 1: Figure S3D), whereas tRNA-3s were almost equally contributed by four different tRNAs (Ala-, Asp-, Glu- and Val-tRNA) (Additional file 1: Figure S3C). To validate these results, we performed quantitative PCR after reverse transcription (RT) of cellular and exosomal RNAs, and confirmed the existence of tsRNAs in exosomes (Fig. 1i and Additional file 3: Table S1, Additional file 2), suggesting that cells could sort cellular tsRNAs into exosomes.

5To explore the potential value of exosomal tsRNA for cancer diagnosis, we extracted and sequenced small RNAs in plasma exosomes from liver cancer patients and healthy donors (Additional file 2). Data analysis showed that miRNA was the most abundant small RNA species within plasma exsomes (Additional file 1: FigureS4A), supporting the role of miRNA as a potential biomarker. Interestingly, tsRNAs still account for 0.2–2% of total small RNAs in plasma exosomes (Additional file 1: Figure S4A). Moreover, tsRNA level in plasma exosomes from liver cancer patients was significantly increased when compared with that from healthy donors (Fig. 2a). Consistent with previous results from cell culture medium, tRNA-5 was also the major type of tsRNA in plasma exosomes (Additional file 1: Figure S4B). Furthermore, significantly higher level of tRNA-5 was observed in plasma exosomes from liver cancer patients than that in healthy controls (Fig. 2b), suggesting abnormal secretion of specific tsRNA into plasma exosomes of cancer patients.

6Similar to the results from cell culture medium, the majority of tRNA-5s ranged from 31 nt to 34 nt in length (Fig. 2c) and tRNA-3s showed a peak at 18 nt (Additional file 1: Figure S5). Intriguingly, the level of tRNA-5s with 30~33 nt in length significantly increased in liver cancer patients (Fig. 2c), indicating that 30~33 nt tRNA-5s could be a potential biomarker for cancer diagnosis. In addition, most tRNA-5s in plasma exosomes were derived from Gly-tRNA, while the richest tRNA-3 was generated from Leu-tRNA (Additional file 1: Figure S6).

7To identify the potential of tsRNAs for biomarker, we compared tsRNA levels in plasma exosomes and found that 46 tsRNAs (35 up-regulated, 11 down-regulated) were differentially expressed between liver cancer patients and healthy donors (Fig. 2d, Additional file 3: Table S2). Moreover, the RNA-seq results were validated individually by RT-qPCR analyses. Compared with healthy donors, the liver cancer patients exhibited significantly higher levels of tRNA-ValTAC-3, tRNA-GlyTCC-5, tRNA-ValAAC-5 and tRNA-GluCTC-5 in the plasma exosomes (Fig. 3, Additional file 3: Table S3), demonstrating that plasma exosomal tsRNA could serve as a novel biomarker for cancer diagnosis. However, more independent patient samples are needed to further validate the relationship between tsRNA and liver cancer.

8In this study, we demonstrated the existence of abundant tsRNA in exosomes from cell culture medium and plasma, representing a novel small RNA species in exosomes. Moreover, the plasma exosomes in liver cancer patients have significantly higher tsRNA level than that in healthy control. Notably, four tsRNAs from plasma exosomes are differentially expressed between liver cancer patients and healthy donors, indicating their great potential as a novel “liquid biopsy” biomarker for cancer diagnosis. Taken together, our study not only expands non-coding RNA species in exosome, but also sheds light on the diagnostic value of tsRNA as a promising biomarker for cancer.