Transcriptome profiling and co-expression network analysis of lncRNAs and mRNAs in colorectal cancer by RNA sequencing

Abstract

Background: Long non-coding RNAs (lncRNAs) are widely involved in the pathogenesis of cancers. However, biological roles of lncRNAs in occurrence and progression of colorectal cancer (CRC) remain unclear. The current study aimed to evaluate the expression pattern of lncRNAs and messenger RNAs (mRNAs). Methods: RNA sequencing (RNA-Seq) in CRC tissues and adjacent normal tissues from 6 CRC patients was performed and functional lncRNA-mRNA co-expression network was constructed afterwards. Gene enrichment analysis was demonstrated using DAVID 6.8 tool. Reverse transcription quantitative polymerase chain reaction (RT-qPCR) was used to validate the expression pattern of differentially expressed lncRNAs. Pearson correlation analysis was applied to evaluate the relationships between selected lncRNAs and mRNAs. Results: One thousand seven hundred and sixteenth differentially expressed mRNAs and 311 differentially expressed lncRNAs were screened out. Among these, 568 mRNAs were up-regulated while 1148 mRNAs down-regulated, similarly 125 lncRNAs were up-regulated and 186 lncRNAs down-regulated. In addition, 1448 lncRNA–mRNA co-expression pairs were screened out from 940,905 candidate lncRNA-mRNA pairs. Gene enrichment analysis revealed that these lncRNA-related mRNAs are associated with cell adhesion, collagen adhesion, cell differentiation, and mainly enriched in ECM-receptor interaction and PI3K-Akt signaling pathways. Finally, RT-qPCR results verified the expression pattern of lncRNAs, as well as the relationships between lncRNAs and mRNAs in 60 pairs of CRC tissues. Conclusions: In conclusion, these results of the RNA-seq and bioinformatic analysis strongly suggested that the dysregulation of lncRNA is involved in the complicated process of CRC development, and providing important insight regarding the lncRNAs involved in CRC.