Construction of a circRNA-Related ceRNA Prognostic Regulatory Network in Breast Cancer
Huan Song,Jian Sun,Weimin Kong,Ye Ji,Dian Xu,Jianming Wang
DOI: https://doi.org/10.2147/OTT.S266507
IF: 4
2020-08-20
OncoTargets and Therapy
Abstract:Huan Song, 1, * Jian Sun, 2, * Weimin Kong, 2, * Ye Ji, 1 Dian Xu, 1 Jianming Wang 1 1 Department of Epidemiology, Center for Global Health, School of Public Health, Nanjing Medical University, Nanjing 211166, People's Republic of China; 2 Department of Thoracic Surgery, The First People's Hospital of Yancheng City, Yancheng 224006, People's Republic of China *These authors contributed equally to this work Correspondence: Jianming Wang Department of Epidemiology Center for Global Health, School of Public Health, Nanjing Medical University, Nanjing 211166, People's Republic of China Tel +86-25-86868438 Email jmwang@njmu.edu.cn Purpose: Accumulating evidence has indicated that circRNAs are closely involved in tumorigenesis and progression of human cancers. However, the molecular mechanism underlying function of circRNAs in breast cancer has not been thoroughly elucidated. Currently, we aimed to characterize the circRNA-related competing endogenous RNA (ceRNA) regulatory network in breast cancer and to construct prognostic model. Materials and Methods: First, we constructed circRNA expression profiles for paired breast cancer in a Chinese population using a human circRNA microarray. Expression profiles of circRNAs, miRNAs, and mRNAs were retrieved from our circRNA dataset, the Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) databases. We applied the limma and edgeR packages to identify differentially expressed RNAs. Weighted gene correlation network analysis (WGCNA) was used to identify critical modules of mRNAs. Next, a ceRNA network was established based on circRNA–miRNA and miRNA–mRNA intersections. Both Cox regression analysis and ROC curve analysis were performed to generate prognostic model. Additionally, we performed Gene Set Enrichment Analysis (GSEA) on prognostic signatures. Results: Total of 59 circRNAs, 98 miRNAs and 3966 mRNAs were identified as differentially expressed RNAs. We first identified 38 miRNA-mRNA pairs and 38 circRNA-miRNA pairs to construct the circRNA–miRNA-mRNA regulatory network and then generated a prognostic model based on 7 signatures (MMD, SLC29A4, CREB5, FOS, ANKRD29, MYOCD, and PIGR), and patients with high-risk scores presented poor prognosis. Several cancer-related pathways were enriched, including the TGF-β pathway, the focal adhesion pathway, and the JAK-STAT signaling pathway, and 20 prognostic ceRNA regulatory networks were subsequently identified. Conclusion: In all, we screened a series of dysregulated circRNAs, miRNAs, and mRNAs, and constructed circRNA-related ceRNA network in breast cancer. Our findings may help to deepen the understanding of circRNA-related regulatory mechanisms. Moreover, we generated a prognostic model that provided new insight into postoperative management for breast cancer. Keywords: ncRNA, ceRNA, circRNA, miRNA, mRNA, breast cancer, TCGA
oncology,biotechnology & applied microbiology