Identification of differentially expressed genes regulated by molecular signature in breast cancer-associated fibroblasts by bioinformatics analysis

Objective Breast cancer is a severe risk to public health and has adequately convoluted pathogenesis. Therefore, the description of key molecular markers and pathways is of much importance for clarifying the molecular mechanism of breast cancer-associated fibroblasts initiation and progression. Breast cancer-associated fibroblasts gene expression dataset was downloaded from Gene Expression Omnibus database. Methods A total of nine samples, including three normal fibroblasts, three granulin-stimulated fibroblasts and three cancer-associated fibroblasts samples, were used to identify differentially expressed genes (DEGs) between normal fibroblasts, granulin-stimulated fibroblasts and cancer-associated fibroblasts samples. The gene ontology (GO) and pathway enrichment analysis was performed, and protein-protein interaction (PPI) network of the DEGs was constructed by NetworkAnalyst software. Results Totally, 190 DEGs were identified, including 66 up-regulated and 124 down-regulated genes. GO analysis results showed that up-regulated DEGs were significantly enriched in biological processes (BP), including cell-cell signalling and negative regulation of cell proliferation; molecular function (MF), including insulin-like growth factor II binding and insulin-like growth factor I binding; cellular component (CC), including insulin-like growth factor binding protein complex and integral component of plasma membrane; the down-regulated DEGs were significantly enriched in BP, including cell adhesion and extracellular matrix organization; MF, including N -acetylgalactosamine 4-sulfate 6- O -sulfotransferase activity and calcium ion binding; CC, including extracellular space and extracellular matrix. WIKIPATHWAYS analysis showed the up-regulated DEGs were enriched in myometrial relaxation and contraction pathways. WIKIPATHWAYS, REACTOME, PID_NCI and KEGG pathway analysis showed the down-regulated DEGs were enriched endochondral ossification, TGF beta signalling pathway, integrin cell surface interactions, beta1 integrin cell surface interactions, malaria and glycosaminoglycan biosynthesis—chondroitin sulfate/dermatan sulphate. The top 5 up-regulated hub genes, CDKN2A, MME, PBX1, IGFBP3, and TFAP2C and top 5 down-regulated hub genes VCAM1, KRT18, TGM2, ACTA2, and STAMBP were identified from the PPI network, and subnetworks revealed these genes were involved in significant pathways, including myometrial relaxation and contraction pathways, integrin cell surface interacti