Construction of a Myc-associated ceRNA network reveals a prognostic signature in hepatocellular carcinoma

Abstract

Hepatocellular carcinoma (HCC) remains an extremely lethal disease worldwide. High-throughput methods have revealed global transcriptome dysregulation; however, a comprehensive investigation of the complexity and behavioral characteristics of the competing endogenous RNA (ceRNA) network in HCC is lacking. In this study, we extracted the transcriptome (RNA) sequencing data of 371 HCC patients from The Cancer Genome Atlas platform. With the comparison of the high Myc expression (Mychigh) tumor and low Myc expression (Myclow) tumor groups in HCC, we identified 1,125 differentially expressed (DE) mRNAs, 589 long non-coding RNAs (lncRNAs), and 93 microRNAs (miRNAs). DE RNAs predicted the interactions necessary to construct an associated Myc ceRNA network, including 19 DE lncRNAs, 5 miRNAs, and 72 mRNAs. We identified a significant signature (long intergenic non-protein-coding [LINC] RNA 2691 [LINC02691] and LINC02499) that effectively predicted overall survival and had protective effects. The target genes of microRNA (miR)-212-3p predicted to intersect with DE mRNAs included SEC14-like protein 2 (SEC14L2) and solute carrier family 6 member 1 (SLC6A1), which were strongly correlated with survival and prognosis. With the use of the lncRNA-miRNA-mRNA axis, we constructed a ceRNA network containing four lncRNAs (LINC02691, LINC02499, LINC01354, and NAV2 antisense RNA 4), one miRNA (miR-212-3p), and two mRNAs (SEC14L2 and SLC6A1). Overall, we successfully constructed a mutually regulated ceRNA network and identified potential precision-targeted therapies and prognostic biomarkers.

Keywords

• TCGA
• HCC
• ceRNA
• immune infiltration
• methylation expression