Identification of key lncRNAs in the carcinogenesis and progression of colon adenocarcinoma by co-expression network analysis

Colon adenocarcinoma (COAD) is one of the most common cancers, and its carcinogenesis and progression is influenced by multiple long non-coding RNAs (lncRNA), especially through the miRNA sponge effect. In this study, more than 4000 lncRNAs were re-annotated from the microarray datasets through probe sequence mapping to obtain reliable lncRNA expression profiles. As a systems biology method for describing the correlation patterns among genes across microarray samples, weighted gene co-expression network analysis was conducted to identify lncRNA modules associated with the five stepwise stages from normal colonic samples to COAD (n = 94). In the most relevant module (R-2 = -0.78, P = 4E-20), four hub lncRNAs were identified (CTD-2396E7.11, PCGF5, RP11-33O4.1, and RP11-164P12.5). Then, these four hub lncRNAs were validated using two other independent datasets including GSE20916 (n = 145) and GSE39582 (n = 552). The results indicated that all hub lncRNAs were significantly negatively correlated with the three-stage colonic carcinogenesis, as well as TNM stages in COAD (one-way analysis of variance P < 0.05). Kaplan-Meier survival curve showed that patients with higher expression of each hub lncRNA had a significantly higher overall survival rate and lower relapse risk (log-rank P < 0.05). In conclusion, through co-expression analysis, we identified and validated four key lncRNAs in association with the carcinogenesis and progression of COAD, and these lncRNAs might have important clinical implications for improving the risk stratification, therapeutic decision and prognosis prediction in COAD patients.