Construction of ceRNA network and identification of hub differentially expressed genes associated with breast cancer using reanalysis of microarray dataset: A systems biology approach

The interaction between long non-coding RNAs (lncRNAs), miRNAs and mRNAs has implications in the patho-genesis of different cancer, including breast cancer. In the current study, we developed an in-silico approach to ascertain the competing endogenous RNA (ceRNA) network in breast cancer. Our approach led to identification of 1816 differentially expressed (DE) mRNAs, including 1039 downregulated DEmRNAs (such as LEP and ADIPOQ) and 777 upregulated DEmRNAs (such as COL11A1 and COL10A1), 19 DElncRNAs, including 15 downregulated DElncRNAs (such as CARMN and COPG2IT1) and 4 upregulated DElncRNAs (such as MALAT1 and NRAV) and 27 DEmiRNAs, including 15 downregulated DEmiRNAs (such as MIR452 and MIR224) and 12 upregulated DEmiRNAs (such as MIR6787 and MIR21). Pathway analysis revealed down-regulation of PPAR, Fatty acid metabolism, Adipocytokine, Vascular smooth muscle contraction and Metabolism of xenobiotics by cytochrome P450, while up-regulation of Pyrimidine metabolism, p53 signaling pathway, Cell cycle, Oocyte meiosis and RNA transport pathways in breast cancer. Finally, we constructed an lncRNA/miRNA/mRNA ceRNA network consisted of 2 lncRNAs, 15 mRNAs, and 4 miRNAs. This network represents an appropriate target for design of anti-cancer modalities and documentation of novel markers for breast cancer.