A Competing Endogenous RNA Network Reveals Novel lncRNA, miRNA and mRNA Biomarkers With Diagnostic and Prognostic Value for Early Breast Cancer

Background: This study aims to reveal early breast cancer (BC) specific competing endogenous RNA (ceRNA) network through the expression profiles of microRNAs (miRNAs), long non-coding RNAs (lncRNAs) and mRNAs. Methods: Based on The Cancer Genome Atlas (TCGA), we obtained the differentially expressed mRNAs, miRNAs, and lncRNAs (DEmRNAs, DEmiRNAs and DElncRNAs) between early BC and normal samples. The lncRNA-miRNA-mRNA interaction network was constructed using Cytoscape. Functional enrichment were performed using GeneCoDis3. The expression of selected genes were validated by qRT-PCR. Based on the published dataset, we validated the result of TCGA integration analysis. The diagnostic and prognostic value of candidate genes was evaluated by ROC curve analysis and survival analysis, respectively. Results: Totally, 1207 DEmRNAs, 194 DElncRNAs and 37 DEmiRNAs were obtained. Functional enrichment analysis results showed that all of DEmRNAs were enriched in pathway of cytokine-cytokine receptor interaction, PPAR signaling pathway and pathways in cancer. The DEmRNA-DEmiRNA-DElncRNA interaction network in early BC was consisted of 23 DEmiRNAs, 95 DElncRNAs and 309 DEmRNAs. Among ceRNA network, IL-6-hsa-miR-182-5p-ADAMTS9-AS1 interactions, LIFR-hsa-miR-21-5p-ADAMTS9-AS1 interactions and MMP1/MMP11-hsa-miR-145-5p-CDKN2B-AS1 interactions were speculated to involve in the development of early BC. The qRT-PCR results were consistent with our integrated analysis. Except for ADAMTS9-AS1 and CDKN2B-AS1, expression of the others results in the Gene Expression Omnibus (GEO) dataset were generally consistent with TCGA integrated analysis. The area under curve (AUC) of the ADAMTS9-AS1, CDKN2B-AS1, IL-6, MMP11, hsa-miR-145-5p and hsa-miR-182-5p were 0.947, 0.862, 0.842, 0.993, 0.960 and 0.944, and the specificity and sensitivity of the 6 biomarkers were 83.4% and 95.6%, 72.2% and 90.3%, 80.1% and 74.3%, 96.2% and 96.5%, 90.1% and 92.3%, and 88.7% and 90.4%, respectively. In addition, IL-6 had potential prognostic value for early BC. Conclusion: These findings may provide novel insights into the lncRNA-miRNA-mRNA network and uncover potential therapeutic targets in early BC.