Abstract
Long non‑coding RNAs (lncRNAs) have been shown to contribute to progression and prognosis of hepatocellular carcinoma (HCC). However, expression profiling and interaction of lncRNAs with messenger RNAs (mRNAs) and microRNAs (miRNAs) remain largely unknown in HCC. The expression profiling of lncRNAs, mRNA and miRNAs was obtained using microarray. The Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes pathway analysis were used to characterize potential functions of differentially expressed mRNAs. Cytoscape was applied to construct an lncRNA‑miRNA‑mRNA co‑expression network and candidate lncRNAs were validated via quantitative PCR in 30 pairs of HCC and adjacent tumor‑free tissues. In this study, 1,056 upregulated and 1,288 downregulated lncRNAs were identified, while 2,687 mRNAs and 6 miRNAs were aberrantly expressed in HCC compared with adjacent tumor‑free tissues. Potential functions of differentially expressed mRNAs were demonstrated to significantly participate in modulating critical genes in the cell cycle, such as cyclin E1 and cyclin B2. After screening, 95 lncRNAs, 5 miRNAs and 36 mRNAs were recruited for construction of lncRNA‑mRNA‑miRNA co‑expression network in the cell cycle pathway. Subsequently, the top 5 lncRNAs that potentially modulate critical genes in the cell cycle were selected as the candidates for further verification. Kaplan‑Meier curves using the Cancer Genome Atlas database showed that 13 targeted mRNAs were associated with overall survival of HCC patients. Finally, three lncRNAs, including ENST00000522221, lnc‑HACE1‑6:1 and lnc‑ICOSLG‑11:1, are significantly upregulated in HCC tissues compared with adjacent tumor‑free tissues. These findings suggest that lncRNAs play essential roles in the pathogenesis of HCC via regulating coding genes and miRNAs, and may be important targets for diagnosis and treatment of this disease.