Identification of immune-related lncRNAs to improve the prognosis prediction for patients with papillary thyroid cancer.

Objective: To identify immune-related long non-coding RNAs (lncRNAs) in papillary thyroid cancer (PTC).

Methods: The Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) databases were used to obtain the gene expression profile. Immune-related lncRNAs were screened from the Molecular Signatures Database v4.0 (MsigDB). We performed a survival analysis of critical lncRNAs. Further, the function of prognostic lncRNAs was inferred using the Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) to clarify the possible mechanisms underlying their predictive ability. The assessment was performed in clinical samples and PTC cells.

Results: We obtained 4 immune-related lncRNAs, 15 microRNAs (miRNAs), and 375 mRNAs as the key mediators in the pathophysiological processes of PTC from the GEO database. Further, Lasso regression analysis identified seven prognostic markers (LINC02550, SLC26A4-AS1, ACVR2B-AS1, AC005479.2, LINC02454, and AL136366.1), most of which were related to tumor development. The KEGG pathway enrichment analysis showed different, changed genes mainly enriched in the cancer-related pathways, PI3K-Akt signaling pathway, and focal adhesion. Only SLC26A4-AS1 had an intersection in the results of the two databases.

Conclusion: LncRNA SLC26A4-AS1, which is the most associated with prognosis, may play an oncogenic role in the development of PTC.