Abstract
OBJECTIVE: TRNA-derived fragments (tRFs) and tRNA-derived stress-induced RNAs (tiRNAs) are recognized as novel and potential types of non-coding RNAs (ncRNAs), and several tRF/tiRNA signatures are closely associated with tumor diagnosis. This study aimed to analyze the expression profiles of plasma tRFs/tiRNAs and to clarify their diagnostic value in lung adenocarcinoma (LUAD). METHODS: The differential expression profiles of plasma tRFs/tiRNAs in patients with four patients with early LUAD, four patients with advanced LUAD, and four healthy controls were analyzed using high-throughput sequencing technology. Then, plasma tRFs/tiRNAs were validated by quantitative real-time polymerase chain reaction (qRT-PCR), and their diagnostic efficiency was appraised by receiver operating characteristic curve analysis. The correlation of candidate plasma tRFs/tiRNAs with clinicopathological features was also analyzed. Finally, bioinformatics analysis was performed to explore and identify the potential biological pathways induced by tRFs/tiRNAs. RESULTS: The sequencing results revealed that tRFs/tiRNAs from plasma samples in patients with LUAD were differently expressed, supporting the necessity of exploring their potential as biomarkers. The validation results of qRT-PCR demonstrated that the expression level of tRF-1:29-Pro-AGG-1-M6 was downregulated in LUAD, while that of tRF-55:76-Tyr-GTA-1-M2 was upregulated, which was consistent with the sequencing data. The areas under the receiver operating characteristic curve of tRF-1:29-Pro-AGG-1-M6 and tRF-55:76-Tyr-GTA-1-M2 were 0.882 and 0.896, respectively, which have significant values in the diagnosis of LUAD. The expressions of tRF-1:29-Pro-AGG-1-M6 and tRF-55:76-Tyr-GTA-1-M2 in LUAD were obviously correlated with various clinicopathological features such as tumor-node-metastasis stage, node stage, and the expression levels of carcinoembryonic antigen. In addition, their expression was significantly altered from before to after tumor resection in LUAD patients. The results of Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses further indicated that tRF-1:29-Pro-AGG-1-M6 and tRF-55:76-Tyr-GTA-1-M2 are widely distributed and apparently enriched in several tumor-related signaling pathways. CONCLUSIONS: Plasma tRF-1:29-Pro-AGG-1-M6 and tRF-55:76-Tyr-GTA-1-M2 may be promising components in the development of highly sensitive and non-invasive biomarkers for LUAD diagnosis.