Abstract
The bidirectional causal relationship between type 2 diabetes mellitus (T2DM) and breast cancer (BC) has been established by numerous epidemiological studies. However, the underlying molecular mechanisms are not yet fully understood. Identification of hub genes implicated in T2DM-BC molecular crosstalk may help elucidate on the causative mechanisms. For this, expression series GSE29231 (T2DM-adipose tissue), GSE70905 (BC- breast adenocarcinoma biopsies) and GSE150586 (diabetes and BC breast biopsies) were extracted from Gene Expression Omnibus (GEO) database, and analyzed to obtain differentially expressed genes (DEGs). The overlapping DEGs were determined using FunRich. Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) and Transcription Factor (TF) analyses were performed on EnrichR software and a protein-protein interaction (PPI) network was constructed using STRING software. The network was analyzed on Cytoscape to determine hub genes and Kaplan-Meier plots were obtained. A total of 94 overlapping DEGs were identified between T2DM and BC samples. These DEGs were mainly enriched for GO terms RNA polymerase II core promoter proximal region sequence and its DNA binding, and cAMP response element binding protein, and KEGG pathways including bladder cancer, thyroid cancer and PI3K-AKT signaling. Eight hub genes were identified: interleukin 6 (IL6), tumor protein 53 (TP53), interleukin 8 (CXCL8), MYC, matrix metalloproteinase 9 (MMP9), beta-catenin 1 (CTNNB1), nitric oxide synthase 3 (NOS3) and interleukin 1 beta (IL1β). MMP9 and MYC associated unfavorably with overall survival (OS) in breast cancer patients, IL6, TP53, IL1β and CTNNB1 associated favorably, whereas NOS3 did not show any correlation with OS. Salt inducible kinase 1 (SIK1) was identified as a significant key DEG for comorbid samples when compared with BC, also dysregulated in T2DM and BC samples (adjusted p <0.05). Furthermore, four of the significant hub genes identified, including IL6, CXCL8, IL1B and MYC were also differentially expressed for comorbid samples, however at p < 0.05. Our study identifies key genes including SIK1, for comorbid state and 8 hub genes that may be implicated in T2DM-BC crosstalk. However, limitations associated with the insilico nature of this study necessitates for subsequent validation in wet lab. Hence, further investigation is crucial to study the molecular mechanisms of action underlying these genes to fully explore their potential as diagnostic and prognostic biomarkers and therapeutic targets for T2DM-BC association.