Comprehensive analysis of the lncRNA-associated competing endogenous RNA network in breast cancer

Long noncoding RNAs (lncRNAs) have been confirmed to be potential prognostic markers in a variety of cancers and to interact with microRNAs (miRNAs) as competing endogenous RNAs (ceRNAs) to regulate target gene expression. However, the role of lncRNA-mediated ceRNAs in breast cancer (BC) remains unclear. In the present study, a ceRNA network was generated to explore their role in BC. The expression profiles of mRNAs, miRNAs and lncRNAs in 1,109 BC tissues and 113 normal breast tissues were obtained from The Cancer Genome Atlas database (TCGA). A total of 3,198 differentially expressed (DE) mRNAs, 150 differentially DEmiRNAs and 1,043 DElncRNAs were identified between BC and normal tissues. A lncRNA-miRNA-mRNA network associated with BC was successfully constructed based on the combined data obtained from RNA databases, and comprised 97 lncRNA nodes, 24 miRNA nodes and 74 mRNA nodes. The biological functions of the 74 DEmRNAs were further investigated by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. The results demonstrated that the DEmRNAs were significantly enriched in two GO biological process categories; the main biological process enriched term was ‘positive regulation of GTPase activity’. By KEGG analysis, four key enriched pathways were obtained, including the ‘MAPK signaling pathway’, the ‘Ras signaling pathway’, ‘prostate cancer’, and the ‘FoxO signaling pathway’. Kaplan-Meier survival analysis revealed that six DElncRNAs (INC AC112721.1, LINC00536, MIR7-3HG, ADAMTS9-AS1, AL356479.1 and LINC00466), nine DEmRNAs (KPNA2, RACGAP1, SHCBP1, ZNF367, NTRK2, ORS1, PTGS2, RASGRP1 and SFRP1) and two DEmiRNAs (hsa-miR-301b and hsa-miR-204) had significant effects on overall survival in BC. The present results demonstrated the aberrant expression of INC AC112721.1, AL356479.1, LINC00466 and MIR7-3HG in BC, indicating their potential prognostic role in patients with BC.