AUTHOR=Guo Jun-Nan , Xia Tian-Yi , Deng Shen-Hui , Xue Wei-Nan , Cui Bin-Bin , Liu Yan-Long 

TITLE=Prognostic Immunity and Therapeutic Sensitivity Analyses Based on Differential Genomic Instability-Associated LncRNAs in Left- and Right-Sided Colon Adenocarcinoma

JOURNAL=Frontiers in Molecular Biosciences

VOLUME=8

YEAR=2021

URL=https://www.frontiersin.org/journals/molecular-biosciences/articles/10.3389/fmolb.2021.668888

DOI=10.3389/fmolb.2021.668888

ISSN=2296-889X

ABSTRACT=<p><bold>Background:</bold> The purpose of our study was to develop a prognostic risk model based on differential genomic instability-associated (DGIA) long non-coding RNAs (lncRNAs) of left-sided and right-sided colon cancers (LCCs and RCCs); therefore, the prognostic key lncRNAs could be identified.</p><p><bold>Methods:</bold> We adopted two independent gene datasets, corresponding somatic mutation and clinical information from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. Identification of differential DGIA lncRNAs from LCCs and RCCs was conducted with the appliance of “Limma” analysis. Then, we screened out key lncRNAs based on univariate and multivariate Cox proportional hazard regression analysis. Meanwhile, DGIA lncRNAs related prognostic model (DRPM) was established. We employed the DRPM in the model group and internal verification group from TCGA for the purpose of risk grouping and accuracy verification of DRPM. We also verified the accuracy of key lncRNAs with GEO data. Finally, the differences of immune infiltration, functional pathways, and therapeutic sensitivities were analyzed within different risk groups.</p><p><bold>Results:</bold> A total of 123 DGIA lncRNAs were screened out by differential expression analysis. We obtained six DGIA lncRNAs by the construction of DRPM, including AC004009.1, AP003555.2, BOLA3-AS1, NKILA, LINC00543, and UCA1. After the risk grouping by these DGIA lncRNAs, we found the prognosis of the high-risk group (HRG) was significantly worse than that in the low-risk group (LRG) (all <italic>p</italic> &lt; 0.05). In all TCGA samples and model group, the expression of CD8<sup>+</sup> T cells in HRG was lower than that in LRG (all <italic>p</italic> &lt; 0.05). The functional analysis indicated that there was significant upregulation with regard to pathways related to both genetic instability and immunity in LRG, including cytosolic DNA sensing pathway, response to double-strand RNA, RIG-Ⅰ like receptor signaling pathway, and Toll-like receptor signaling pathway. Finally, we analyzed the difference and significance of key DGIA lncRNAs and risk groups in multiple therapeutic sensitivities.</p><p><bold>Conclusion:</bold> Through the analysis of the DGIA lncRNAs between LCCs and RCCs, we identified six key DGIA lncRNAs. They can not only predict the prognostic risk of patients but also serve as biomarkers for evaluating the differences of genetic instability, immune infiltration, and therapeutic sensitivity.</p>