Comprehensive analysis of aberrantly expressed profiles of lncRNAs, miRNAs and mRNAs with associated ceRNA network in cholangiocarcinoma

Abstract

Cholangiocarcinoma (CCA) is a highly malignant and poorly differentiated bile duct cancer with an extremely poor prognosis, but the pathogenesis of CCA remains not well-known. Attention has been increasingly focused on long noncoding RNAs, which plays an important role in tumorigenesis. However, the roles of cancer specific lncRNA and its related competitive endogenous RNAs (ceRNA) network in CCA remain elusive. In this study, we comprehensively integrated expression profiles, including data on mRNAs, lncRNAs and miRNAs obtained from 36 CCA tissues and 9 normal tissues in The Cancer Genome Atlas. 1434 cancer specific lncRNAs, 68 miRNAs and 3538 mRNAs (|logFC|> 1, p< 0.05) were identified. Based on bioinformatics generated from miRcode, starBase, miRTarBase, TargetScan and miRDB, we constructed an lncRNA-miRNA-mRNA network (ceRNA network) in CCA. We constructed the lncRNA-miRNA-mRNA ceRNA network consisting of 206 molecules and 454 interactions. In addition, we used Cytoscape software to visualize the ceRNA network in WGCNA, 22 mRNA network modules were identified, five of which were significantly related to tumor grade and survival time. Moreover, three lncRNAs COL18A1-AS1, SLC6A1-AS1 and HULC were found to be significantly associated with overall survival. The present study provides novel insight for better understanding of lncRNA-related ceRNA network in CCA and useful resource for identifcation of novel biomarkers of CCA.

Keywords: Cholangiocarcinoma; TCGA; ceRNA; lncRNA; mRNA; miRNA.