Document Type : Original Article(s)
Authors
- Mohadeseh Khoshandam 1
- Mohammad Rahmanian 2, 3
- Mohammad Taghi Hedayati Goudarzi 4
- Hossein Soltaninejad 5
- Sadegh Babashah 6
- Mahdiye Khoshandam 7
1 Department of Molecular Medicine, National Institute of Genetic Engineering and Biotechnology (NIGEB), Tehran, Iran
2 Gastroenterology and Liver Diseases Research Center, Research Institute for Gastroenterology and Liver Diseases, Shahid Beheshti University of Medical Sciences, Tehran, Iran
3 Student Research Committee, School of Medicine, Shahid Beheshti University of Medical Sciences, Tehran, Iran
4 Department of Cardiology, School of Medicine, Babol University of Medical Sciences, Babol, Iran
5 Department of Stem Cells Technology and Tissue Regeneration, Faculty of Interdisciplinary Science and Technologies, Tarbiat Modares University, Tehran, Iran
6 Department of Molecular Genetics, Faculty of Biological Sciences, Tarbiat Modares University, Tehran, Iran
7 Department of Cell Biology and Regenerative Medicine, Academic Center for Education, Culture and Research, Qom Branch, Qom, Iran
Abstract
Background: Breast cancer is the most common form of cancer among women worldwide, and the rates of both new cases and deaths have increased over the past two decades. The aim of the study was to identify and validate molecular pathways that could potentially be targeted for therapeutic interventions.
Methods: The bioinformatics resource WebGestalt was used to determine the functional annotation of the Gene Ontology, as well as enrichment analysis of Reactome and KEGG pathways in 2023-2024. GeneMANIA, a server for assessing protein-gene interactions, co-localization, pathways, co-expression, and protein-domain similarity of target genes and their interacting genes, was evaluated via this web tool. GEO was also used to determine mRNA expression levels in BRCA individuals. R packages were used to screen for differentially expressed genes for both datasets. On the other hand, the open cancer resources GENT2 TNMPlot, UCSCXena, ENCORI platform, BioXpress, OncoDB, OncoMX, and GEPIA2 were used to measure the differential expression of mRNAs in BRCA patients.
Results: Among the genes analyzed, matrix metalloproteinase-9 (MMP9) showed the greatest change. Similarly, matrix metallopeptidase 14 (MMP14) and Endogenous Vascular Endothelial Growth Factor-A (VEGFA) showed significant increases. Other up-regulated genes, including Apolipoprotein E (APOE), Hypoxia-Inducible Factor-1 Alpha (HIF1A), and Tumor Necrosis Factor (TNF) showed minimal expression changes with minor fluctuations. Finally, Interleukin-1 alpha precursor (IL1A) exhibited a slight increase in expression. Validation of gene expression changes through microarray studies on the GSE37751 and GSE42568 datasets provided consistent and significant results for several of the studied genes. GO analysis further revealed significant molecular functions, cellular components, KEGG pathways, and biological processes that were enriched among the differentially expressed genes. Among the top pathways identified based on FDR and P value were receptor binding signaling, regulation of cell migration, the extracellular matrix, and the AGE-RAGE signaling pathway.
Conclusion: The results predict that the hub genes correlated with angiogenesis may serve as potential therapeutic targets or could be biomarkers for breast cancer.
Highlights
Mohadeseh Khoshandam (Google Scholar)
Hossein Soltaninejad (Google Scholar)
Keywords
What’s Known
The literature highlights the central role of angiogenesis in tumor progression, metastasis, and therapeutic resistance.
What’s New
Interaction network analysis identified MMP9, VEGFA, HIF1A, APOE, TNF, IL1A, MMP14, and TNF as central hub genes with extensive interactions across co-expression, physical interactions, and shared pathways. The results predict that the hub genes correlated with angiogenesis may serve as potential therapeutic targets or could be biomarkers for breast cancer.
Introduction
Breast cancer is a global health concern with a significant effect on public health. 1 There were about 2.26 million new cases of breast cancer in 2020, accounting for about 11.6% of all new cancer incidents and surpassing lung cancer as the most common cancer globally. 2 Among individuals diagnosed with early-stage breast cancer, approximately 30% of node-negative breast cancers and up to 70% of node-positive cancers can progress to metastatic disease, depending on various risk factors. 3 Understanding the genetic background of this malignancy is important for early detection, accurate diagnosis, and the development of specific treatments. 4
The formation of new blood vessels from preexisting ones in postnatal life, known as angiogenesis, plays a vital role in both normal bodily functions and pathological conditions. 5 The angiogenesis process is tightly regulated and involves the interaction of various pro- and anti-angiogenic factors. 6 When the balance between these factors is disrupted, it can result in either excessive or inhibited angiogenesis, which plays an important role in the development of different diseases. 7 Tumors cannot expand more than a certain volume unless they establish a blood supply to provide them with the necessary nutrients and oxygen for their growth. 8 Angiogenic growth factors, including epidermal growth factor (EGF), vascular endothelial growth factor (VEGF), fibroblast growth factor (FGF), and transforming growth factor (TGF), are crucial in facilitating angiogenesis and tumor development. 9 Among these factors, VEGF is pivotal in controlling the growth of endothelial cells, which are essential for the formation of new blood vessels. 10 VEGF is expressed in response to some soluble mediators, including growth factors and cytokines. By enhancing endothelial permeability, VEGF stimulates endothelial cell proliferation, survival, migration, cell-cell contact, and lumen formation. VEGF carries out its functions by binding to cell surface receptors, which are a group of membrane tyrosine kinase receptors. 11 By interacting with cell surface receptors, especially at the site of inflammation, it initiates intracellular signaling pathways and activates multiple genes that regulate essential cellular processes involved in angiogenesis. 12 Because tumor angiogenesis is crucial for tumor growth, studying the pathways and hub genes involved in it is seen as a promising approach to developing therapies that can hinder tumor growth and development. 11
Gene network formalism is the primary method used to uncover the role of genes that are overexpressed in cancer samples compared to normal ones. Innovations in in-silico techniques have streamlined research into gene function, disease, and precision medicine at the molecular level by employing an integrated approach to identify a group of genes linked to angiogenesis, migration, inflammation, and apoptosis that have a significant effect on the survival of breast cancer patients. 13 The objective of this study was to conduct a deep bioinformatics analysis of publicly available datasets to pinpoint hub genes and key pathways involved in angiogenesis and migration, particularly in their correlation with the inflammation pathway in breast cancer.
Materials and Methods
This study was conducted at the National Institute of Genetic Engineering and Biotechnology (NIGEB) in collaboration with Tarbiat Modares University of Iran, between 2023 and 2025.
Data Collection and Processing
Statistical analysis of differential gene expression was performed using the Limma-Voom pipeline. The empirical Bayes method implemented in the Limma package was used to moderate standard errors across genes. P values were adjusted for multiple testing using the Benjamini–Hochberg procedure. Genes with an adjusted P<0.05 and absolute log2 fold change greater than 1 were considered significantly differentially expressed.
This dataset initially contained 114 normal samples; however, one sample was excluded due to incomplete metadata, resulting in 113 normal samples and 1109 tumor samples used in the analysis.
WebGestalt (WEB-based GEne SeT AnaLysis Toolkit)
The bioinformatics resource called WebGestalt (https://www.webgestalt.org/) was employed to determine gene ontology (GO) functional annotation (GO terms, including Biological Process and Molecular Function). Moreover, Kyoto Encyclopedia of Genes and Genomes (KEGG) and Reactome pathways enrichment analysis were performed. An adjusted P< 0.05 was set to show a statistically significant difference.
GeneMANIA
GeneMANIA (https://genemania.org/) is a server for evaluating protein and genetic interactions, co-localization, pathways, co-expression, and domain-protein similarity of target genes. In this research, we assessed the relationship between matrix metalloproteinase-9 (MMP9), matrix metallopeptidase 14 (MMP14), Endogenous Vascular Endothelial Growth Factor-A (VEGFA), Apolipoprotein E (APOE), Hypoxia-Inducible Factor-1 Alpha (HIF1A), Tumor necrosis factor (TNF), Interleukin-1 alpha precursor (IL1A), and their interactive genes via this web tool.
Verification of Expression
The GEO database (https://www.ncbi.nlm.nih.gov/geo/) was used to determine the mRNA expression levels in BRCA (BReast CAncer gene) subjects. Two datasets of microarray, namely GSE37751 and GSE42568, were downloaded. R packages (affy and Limma) were used to screen differentially expressed genes for both datasets. P value less than 0.05 and logFC higher and lower than 0 were the filter conditions to explore DEGs between BRCA and normal specimens. On the other hand, the open cancer sources GENT2 (http://gent2.appex.kr/gent2/), TNMplot (https://tnmplot.com/), UCSCXena (https://xenabrowser.net/), ENCORI platform (https://rnasysu.com/encori/index.php), BioXpress (https://hivelab.biochemistry.gwu.edu/bioxpress), OncoDB (https://oncodb.org/), OncoMX (https://www.oncomx.org/), and GEPIA2 (http://gepia2.cancer-pku.cn/#index) were employed to assay differential expression of mRNAs in patients with BRCA.
Results
We performed a comprehensive analysis of angiogenic hub genes (MMP9, MMP14, VEGFA, APOE, HIF1A, TNF, IL1A) to investigate their expression patterns and biological relevance in breast cancer, focusing on their roles in inflammation and cell migration
Differential Gene Expression
Using RNA-seq data from TCGA and subsequent normalization through Limma and Voom packages, we identified significant changes in the expression of angiogenic hub genes. Among the analyzed genes, MMP9 exhibited the highest positive fold change with a logFC of 2.75 (figure 1A and table 1). Similarly, MMP14 showed substantial upregulation with a logFC of 1.07 (figure 1B and table 1). Other notable upregulated genes included VEGFA with a logFC of 0.55 (figure 1C and table 1), and APOE with a logFC of 0.66 (figure 1D and table 1).
Figure 1. The expression of MMP9 in BRCA and normal tissue is shown in panel A. The expression of MMP14 in BRCA and normal tissues is shown in panel B. The expression of VEGFA in BRCA and normal tissues is shown in panel C. The expression of APOE in BRCA and normal tissues is shown in panel D. The expression of APOE in BRCA and normal tissues is shown in panel E. The expression of HIF1A in BRCA and normal tissue is shown in panel F. RCA: Breast invasive carcinoma; MMP9: Matrix metalloproteinase 9; MMP14: Matrix metalloproteinase 14; VEGFA: Vascular endothelial growth factor A; APOE: Apolipoprotein E; TCGA: The cancer genome atlas; HIF1A: Hypoxia-inducible factor 1 Alpha; TNF: Tumor necrosis factor; BRCA: Breast invasive carcinoma. RCA: Breast invasive carcinoma; MMP9: Matrix metalloproteinase 9; MMP14: Matrix metalloproteinase 14; VEGFA: Vascular endothelial growth factor A; APOE: Apolipoprotein E; TCGA: The cancer genome atlas; HIF1A: Hypoxia-inducible factor 1 Alpha; TNF: Tumor necrosis factor; BRCA: Breast invasive carcinoma
| Genes | LogFC | P value | Adjusted P value |
|---|---|---|---|
| MMP9 | 2.7495 | 4.84E-41 | 1.02E-39 |
| VEGFA | 0.5549 | 8.57E-10 | 3.30E-09 |
| APOE | 0.6576 | 2.04E-06 | 5.87E-06 |
| HIF1A | 0.0806 | 0.4244 | 0.4992 |
| TNF | 0.05270 | 0.7505 | 0.8001 |
| MMP14 | 1.07134 | 3.51E-23 | 3.11E-22 |
| IL1A | 0.4097 | 0.1063 | 0.1501 |
| The data includes statistical measures to assess the significance of these changes. MMP9: Matrix metalloproteinase 9; VEGFA: Vascular endothelial growth factor A; APOE: Apolipoprotein E; HIF1A: Hypoxia-inducible factor 1 Alpha; TNF: Tumor necrosis factor; MMP14: Matrix metalloproteinase 14; IL1A: Interleukin 1 Alpha; LOGFC: Log fold change | |||
In contrast, HIF1A displayed minimal expression changes (LogFC=0.08) with non-significant P values (figure 1E and table 1), indicating a limited role under the conditions analyzed. Similarly, TNF exhibited negligible changes in expression (figure 1F and table 1), suggesting no significant involvement in this dataset. Finally, IL1A demonstrated a slight increase in expression (LogFC=0.41), but this change was not statistically significant.
Validation Through Microarray Studies
The validation of gene expression changes through microarray studies in datasets GSE37751 and GSE42568 provided consistent and significant results for several genes. MMP9 showed strong upregulation in both datasets, with a LogFC of 1.41 in GSE37751 and 2.45 in GSE42568 (P=7.83E-05, adjusted P=5.40E-04), confirming its critical role in breast cancer progression. Similarly, MMP14 demonstrated significant upregulation in GSE37751, while in GSE42568, its expression was moderately elevated, providing partial validation of its role.
VEGFA showed moderate upregulation in GSE37751 and higher expression in GSE42568, suggesting its involvement in angiogenesis, although its statistical significance was weaker in GSE42568. APOE displayed slight upregulation in GSE37751, indicating its relevance to lipid metabolism in breast cancer, but its expression in GSE42568 was not significant, reflecting variability across datasets.
On the other hand, HIF1A showed limited but significant upregulation in GSE37751, while its expression in GSE42568 was not statistically significant, suggesting a minor or dataset-specific role. TNF exhibited minimal expression changes in both datasets, with a slight downregulation in GSE37751 and marginal upregulation in GSE42568, confirming no significant involvement. Similarly, IL1A showed negligible changes in expression, with LogFC values of 0.0373 in GSE37751 and 0.0722 in GSE42568, indicating no significant role in breast cancer (table 2).
| Genes | GSE37751 | GSE42568 | ||||
|---|---|---|---|---|---|---|
| LogFC | P value | Adjusted P value | Log FC | P value | Adjusted P value | |
| MMP9 | 1.41 | 4.12E-08 | 7.16E-07 | 2.45 | 7.83E-05 | 5.40E-04 |
| VEGFA | 0.28 | 0.00754 | 0.0237 | 0.418 | 0.0728 | 0.176 |
| APOE | 0.454 | 0.000067 | 0.000426 | 0.0497 | 0.8 | 0.8 |
| HIF1A | 0.264 | 0.0148 | 0.0412 | 0.255 | 0.269 | 0.46 |
| TNF | -0.0426 | 0.48 | 0.62 | 0.241 | 0.23 | 0.41 |
| IL1A | 0.0373 | 0.633 | 0.745 | 0.0722 | 0.688 | 0.81 |
| MMP14 | 0.558 | 1.59E-06 | 0.000017 | 0.244 | 3.56E-02 | 1.01E-01 |
| MMP9: Matrix metalloproteinase 9; VEGFA: Vascular endothelial growth factor A; APOE: Apolipoprotein E; HIF1A: Hypoxia-inducible factor 1 Alpha; TNF: Tumor Necrosis Factor; IL1A: Interleukin 1 Alpha; MMP14: Matrix metalloproteinase 14; GSE37751: Gene expression omnibus series 37751, a specific dataset identifier from microarray studies; GSE42568: Gene expression omnibus series 42568, another specific dataset identifier from microarray studies; logFC: Log fold change | ||||||
Functional Enrichment Analysis
GO analysis revealed significant molecular functions and biological processes enriched among the DEGs.
Molecular Functions
The top enriched activities include “chemoattractant activity” and “chemokine binding”, highlighting the involvement of these genes in immune cell recruitment and inflammatory responses (figure 2). Among the top pathways, signaling receptor binding, signaling receptor activator activity, signaling receptor regulator activity, receptor ligand activity, and molecular function activator activity were identified as superior based on FDR and P value.
Figure 2. Bar plot of the top 20 Molecular Function (MF) terms. Each row shows an enrichment function. The bar graph length shows the enrichment ratio. The enrichment ratio is determined by dividing the number of input genes by the number of background genes.
Biological Processes
The most enriched pathways include “extracellular matrix disassembly” and “cell chemotaxis”, indicating the involvement of these genes in cell migration and metastatic processes (figure 3). Among the top pathways, regulation of cell migration and cell motility, locomotion, regulation of locomotion, and cell migration were identified as superior based on FDR and P value.
Figure 3. The bar plot of the top 20 Biological Process (BP) terms is shown here. Each row shows an enrichment function. The bar graph length shows the enrichment ratio. The enrichment ratio is determined by dividing the number of input genes by the number of background genes.
Cellular Components
Genes were predominantly localized in “serine protease inhibitor complexes” and “collagen-containing extracellular matrix”, underscoring their roles in structural tissue remodeling (figure 4). Among the top pathways, extracellular matrix, external encapsulating structure, collagen-containing extracellular matrix, platelet alpha granule lumen, and cell surface were identified as superior based on FDR and P value.
Figure 4. The bar plot of the top 20 Cellular Component (CC) terms is shown here. Each row represents an enrichment function. Each row shows an enrichment function. The bar graph length shows the enrichment ratio. The enrichment ratio is calculated by dividing the count of input genes by the count of background genes.
KEGG Pathways
The genes were enriched in pathways, including “Bladder cancer”, “AGE-RAGE signaling in diabetic complications”, and “TNF signaling”, reflecting their diverse functional implications in cancer progression and angiogenesis (figure 5). Among the top pathways, the AGE-RAGE signaling pathway in diabetic complications, Rheumatoid arthritis, TNF signaling pathway, Fluid shear stress and atherosclerosis, and Pathways in cancer were identified as superior based on FDR and P value.
Figure 5. The bar chart of enrichment ratios of the top 20 significantly enriched KEGG pathways is shown in this figure. The enrichment ratio represents the proportion of observed genes to the expected genes within the gene list for the KEGG category.
Network Analysis
Interaction network analysis identified MMP9, VEGFA, HIF1A, APOE, TNF, IL1A, MMP14, and TNF as central hub genes with extensive interactions across co-expression, physical interactions, and shared pathways. VEGFA is primarily involved in the vascular endothelial growth factor receptor signaling pathway. HIF1A contributes to the positive regulation of lipid metabolic processes and small-molecule metabolic processes. APOE plays a role in the positive regulation of lipid metabolic processes and contributes to acylglycerol homeostasis. In addition, APOE is crucial in the particle remodeling of triglyceride-rich lipoprotein. TNF plays a role in regulating lipid metabolism, particularly in inflammatory responses, which can influence lipid storage and breakdown. Additionally, it is involved in regulating steroid metabolism, often in the context of stress and inflammatory signaling also participates in regulating various small-molecule metabolic pathways, such as glucose and amino acid metabolism, often during immune responses. These genes are pivotal in regulating angiogenesis, cell migration, and inflammatory responses in breast cancer (figure 6).
Figure 6. The MMP-9, MMP-14, VEGFA, APOE, HIF1A, TNF, and IL1A genes’ interaction network and related functions were obtained by GenMANIA analysis.
Cross-Database Validation
Validation of expression trends was performed using various databases, including UALCAN, GENT2, TNMPlot, and UCSC Xena. Genes such as MMP9, MMP14, VEGF, APOE, HIF1A, TNF, and IL1A were consistently upregulated, while prostaglandin-endoperoxide synthase 2 (PTGS2), B-cell leukemia/lymphoma 2 protein (BCL2), and interleukin 6 (IL6) were downregulated among the datasets that were not included in our study.
Additional Pathway Analysis
Further exploration of pathways emphasized the critical roles of angiogenic hub genes in breast cancer progression. The overlap between bladder cancer and breast cancer highlighted shared mechanisms in ECM remodeling and immune modulation, involving components such as MMP9, VEGFA, and TNF. The AGE-RAGE signaling pathway was relevant to chronic inflammation and angiogenesis, driven by VEGFA and TNF. The interleukin 17 (IL-17) signaling pathway involves genes such as IL1A, IL1B, and TNF, amplifying inflammatory responses and promoting tumor progression. The fluid shear stress and atherosclerosis pathway highlighted vascular adaptations in the tumor microenvironment, mediated by VEGFA and TNF. The HIF-1 signaling pathway played a key role in hypoxia responses through HIF1A, VEGFA, and MMP9.
Discussion
Interaction network analysis in our study identified MMP9, VEGFA, HIF1A, APOE, TNF, IL1A, MMP14, and TNF genes as core genes with extensive interactions in co-expression, physical interactions, and shared pathways. The results predicted that core genes related to angiogenesis may serve as potential therapeutic targets or biomarkers for breast cancer.
The breast cancer management continues to pose significant challenges. Understanding the protein families and pathways that contribute to the development and progression of cancer is a vital step in the development of targeted treatments.
VEGF-A is a well-known mediator of tumor angiogenesis in many solid tumors and has been used as an angiogenesis inhibitor in breast cancer. 14 Additionally, high levels of MMP-14 and MMP-9 expression are associated with increased TNF expression in the tumor stroma and different survival rates compared to each other. 15 - 17 Various studies have shown that increased levels of inflammatory factors, such as TNF, contribute to the promotion of breast cancer stem-like cells. 18 It was also demonstrated that TNF enhances breast cancer stem-like cells by upregulating TAZ expression through the non-canonical signaling (NCS) of the NF-κB pathway. Understanding the signaling cascade initiated by VEGF and the subsequent effect of the HIF gene on the genome of cancer cells may lead to breakthroughs in personalized treatment of solid cancers, including breast cancer. 19 , 20 APOE has been regarded as a potential biomarker that predicts the prognosis of cancer and has been associated with immune infiltration. 21 IL-1A produced by tumor cells influenced breast cancer progression and metastasis in mice by modulating the tumor immune microenvironment. 22
In this study, we investigated the relationship between VEGF, MMP-9, MMP-14, APOE, HIF1A, TNF, and IL-1A and breast cancer, and their upregulation in breast cancer cells was observed. Pathways involved in 1109 breast cancer tumor samples were compared with 113 normal samples, and the results were confirmed using online databases.
We also explored the MAPK signaling pathway, TNF signaling pathway, Rap1 signaling pathway, Ras signaling pathway, IL-17 signaling pathway, PI3K-Akt signaling pathway, lipid and atherosclerosis pathway, NF-kappa B signaling pathway, and HIF-1 signaling pathway.
MAPK signaling is associated with various types of cancer in humans, and its activation is influenced by a variety of external and internal signals and proteins. 23 Consequently, concentrating on these aspects of the signaling pathway often results in adverse effects, activation of alternative mechanisms, and reduced efficacy of the drug. This highlights the importance of a more comprehensive therapeutic approach. To mitigate these detrimental outcomes, it is advisable to explore alternative strategies that center on RAF modulators. 24 The spatiotemporal attributes of the MAPK pathway components could provide a novel perspective.
Structure, upstream activators, and downstream factors play a critical role in Erk (MAPK) kinase, a vital component of the Ras/Raf/MAPK pathway, which is responsible for regulating cellular metabolism in cancer cells. 25
The Ras/Raf/MAPK pathway can be activated through two mechanisms: a ligand-dependent mechanism, in which ligands such as hormones, cytokines, or growth factors bind physically to the receptors; and a ligand-independent mechanism, which can be initiated by physical stressors such as injury, radiation, or changes in osmotic pressure. This pathway plays a role in the growth of cancerous cells in humans. 25 , 26
Aside from mutations in the Ras gene, disruptions in components that regulate Ras can also affect Raf activation. During various biological processes, different growth factors, such as platelet-derived growth factor BB (PDGF-BB), EGF, TGF-α, and VEGF, can stimulate the canonical Ras-Raf-MEK-Erk pathway. 27 Further studies on Raf activation will help reveal the important molecular mechanisms driving cell proliferation, and metastasis and survival in cancer, especially through the roles of small Ras-GTPases and the EGF receptor (EGFR). As such, extensive research is currently underway to identify Raf kinase as a promising target for the development of anticancer drugs. 28 , 29
A small proportion of resistant melanomas depends on compensatory PI3K (phosphatidylinositol 3-kinase) and Akt (protein kinase B) signaling. 30 Early clinical trials involving PI3K inhibitors alongside vemurafenib demonstrated promising results. 31 EGFR has shown efficacy in BRAFV600E-mutated CRCs, by activating the PI3K/AKT pathway serving as a pathway for both acquired and innate resistance to BRAF inhibitors. 32
The resistance to chemotherapy caused by hypoxia poses a significant challenge in treating solid tumors. HIF1α, a gene that helps the body adapt to low oxygen levels, was found to have a profound impact on the efficacy of chemotherapy for various types of cancer and specific treatment plans. 33 The results of multivariate and univariate logistic regression analysis indicated that high levels of HIF1α were associated with lower pathologic complete response rates and a poorer prognosis. 34 The analysis of the geospatial dataset also revealed a negative correlation between HIF1α expression and the response to neoadjuvant drugs in breast cancer patients. 35
KEGG enrichment of DEGs of different HIF1α expression groups in the TCGA database revealed that HIF1α is involved in the IL-17 signaling pathway. 36 Additionally, correlation analysis showed a positive relationship between HIF1α and the IL-17 pathway. 37 Laboratory tests found that the HIF1α/IL-17 pathway affects the sensitivity of breast cancer cells to paclitaxel. 37 Furthermore, correlation analysis suggested that high expression of HIF1α/IL-17A/ C-X-C motif chemokine ligand 10 (CXCL10) genes is positively associated with neutrophil infiltration in breast cancer. 37
NF-κB serves as a key mediator of inflammation and cancer, with known implications in breast cancer and endocrine disorders. NF-κB is regulated by the canonical signaling (CS) pathway (activated by lipopolysaccharide [LPS], Reactive oxygen species [ROS], TNF-α, and IL-1) and the NCS pathways (triggered by inflammatory stimuli of inhibitor of nuclear factor kappa-B kinase subunit alpha [IKK-α]). 38 Dysregulation of NF-κB can lead to cancer, and its activation is observed in different cancers, such as breast cancer.
The NF-κB role in tumor development is intricate due to the multiple functions of the tumor microenvironment. 39 Upregulation of NF-κB in tumor cells can induce inflammation in the tumor microenvironment, leading to increased cytokine gene expression and cytokine release. 40 This results in activating the classical NF-κB signaling pathway and subsequent apoptosis of transformed cells. 41 Hence, NF-κB may serve as a potential marker and therapeutic target for breast cancer prediction and treatment. Despite the cytotoxic effects of drugs in tumor eradication, they are limited by potential side effects.
Although there is an imbalance between the number of tumors (n=1109) and normal (n=113) samples in the TCGA dataset, the Limma-Voom framework is well-suited for such scenarios. It uses precision weights to model the mean-variance relationship in RNA-seq data, which improves the reliability of differential expression estimates even when group sizes are unequal. Furthermore, the empirical Bayes moderation implemented in Limma stabilizes variance estimates across genes, mitigating the effects of sample size disparity. Therefore, while the imbalance is acknowledged as a potential limitation, the statistical framework employed here helps minimize its practical impact on the results.
Conclusion
The development of breast cancer is influenced by a complex interaction between environmental and genetic factors. Researchers found that specific genes, referred to as regulatory hub genes, could serve as potential targets for drug-based therapies. This discovery was made by analyzing extensive datasets of microarray and RNA sequencing data that are linked to pathways such as angiogenesis, migration, and inflammation. These identified hub genes are part of the protein family and have the potential to play a role in the progression of breast cancer. By recognizing these genes as potential targets for the treatment of breast cancer, new and more advanced therapeutic strategies could be developed.
Acknowledgment
This research was financially supported by the Iranian National Science Foundation (INSF) under grant no. 4039338.
Authors’ Contribution
Mo.Kh: Study design, visualization, validation, resources, data collection, drafting, and reviewing the manuscript; M.R: Study design, visualization, validation, resources, data collection, drafting, and reviewing the manuscript; MT.HGh: Conceptualization, formal analysis, and reviewing the manuscript, H.S: Formal analysis, study design, drafting, and reviewing the manuscript; S.B: Formal analysis and reviewing the manuscript; Ma.Kh: Study design, Formal analysis, and reviewing the manuscript. All authors have read and approved the final manuscript and agree to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved.
Declaration of AI
The authors declare that no AI tools were used in the preparation of this manuscript.
Conflict of Interest
None declared.
References
- Khoshandam M, Soheili ZS, Hosseinkhani S, Samiee S, Latifi-Navid H, Ahmadieh H, et al. In vivo inhibition of angiogenesis by htsFLT01/MiRGD nano complex. Transl Oncol. 2025; 56:102400. Publisher Full Text | DOI | PubMed [ PMC Free Article ]
- Wu B, Li Y, Shi B, Zhang X, Lai Y, Cui F, et al. Temporal trends of breast cancer burden in the Western Pacific Region from 1990 to 2044: Implications from the Global Burden of Disease Study 2019. J Adv Res. 2024; 59:189-99. Publisher Full Text | DOI | PubMed [ PMC Free Article ]
- Cherny NI, Paluch-Shimon S, Berner-Wygoda Y. Palliative care: needs of advanced breast cancer patients. Breast Cancer (Dove Med Press). 2018; 10:231-43. Publisher Full Text | DOI | PubMed [ PMC Free Article ]
- DeBerardinis RJ. Tumor Microenvironment, Metabolism, and Immunotherapy. N Engl J Med. 2020; 382:869-71. DOI | PubMed
- Manning D, Rivera EJ, Santana LF. The life cycle of a capillary: Mechanisms of angiogenesis and rarefaction in microvascular physiology and pathologies. Vascul Pharmacol. 2024; 156:107393. Publisher Full Text | DOI | PubMed [ PMC Free Article ]
- Shaw P, Dwivedi SKD, Bhattacharya R, Mukherjee P, Rao G. VEGF signaling: Role in angiogenesis and beyond. Biochim Biophys Acta Rev Cancer. 2024; 1879:189079. DOI | PubMed
- Akbarian M, Bertassoni LE, Tayebi L. Biological aspects in controlling angiogenesis: current progress. Cell Mol Life Sci. 2022; 79:349. Publisher Full Text | DOI | PubMed [ PMC Free Article ]
- Sriharikrishnaa S, Suresh PS, Prasada KS. An introduction to fundamentals of cancer biology. In: Optical Polarimetric Modalities for Biomedical Research. Cham: Springer; 2023. DOI
- Liu ZL, Chen HH, Zheng LL, Sun LP, Shi L. Angiogenic signaling pathways and anti-angiogenic therapy for cancer. Signal Transduct Target Ther. 2023; 8:198. Publisher Full Text | DOI | PubMed [ PMC Free Article ]
- Perez-Gutierrez L, Ferrara N. Biology and therapeutic targeting of vascular endothelial growth factor A. Nat Rev Mol Cell Biol. 2023; 24:816-34. DOI | PubMed
- Arya KR, Bharath Chand RP, Abhinand CS, Nair AS, Oommen OV, Sudhakaran PR. Identification of Hub Genes and Key Pathways Associated with Anti-VEGF Resistant Glioblastoma Using Gene Expression Data Analysis. Biomolecules. 2021; 11Publisher Full Text | DOI | PubMed [ PMC Free Article ]
- Mannan A, Dhiamn S, Garg N, Singh TG. Pharmacological modulation of Sonic Hedgehog signaling pathways in Angiogenesis: A mechanistic perspective. Dev Biol. 2023; 504:58-74. DOI | PubMed
- Kaddoura R, Alqutami F, Asbaita M, Hachim M. In Silico Analysis of Publicly Available Transcriptomic Data for the Identification of Triple-Negative Breast Cancer-Specific Biomarkers. Life (Basel). 2023; 13Publisher Full Text | DOI | PubMed [ PMC Free Article ]
- Khoshandam M, Soheili ZS, Hosseinkhani S, Samiee S, Latifi-Navid H, Kalhor N, et al. Leveraging the htsFLT01/MiRGD Complex to Enhance Apoptosis and Suppress Angiogenesis in MCF7 Breast Cancer Cells. Iran J Med Sci. 2025; 50:707-12. Publisher Full Text | DOI | PubMed [ PMC Free Article ]
- Kwon MJ. Matrix metalloproteinases as therapeutic targets in breast cancer. Front Oncol. 2022; 12:1108695. Publisher Full Text | DOI | PubMed [ PMC Free Article ]
- Luchian I, Goriuc A, Sandu D, Covasa M. The Role of Matrix Metalloproteinases (MMP-8, MMP-9, MMP-13) in Periodontal and Peri-Implant Pathological Processes. Int J Mol Sci. 2022; 23Publisher Full Text | DOI | PubMed [ PMC Free Article ]
- Eugeniu C, Eremei Z, Claudiu M, Radu N, Ruslan P. Immunoexpression of matrix metalloproteinases MMP-1, MMP-2, MMP-9 and MMP-14 in extragenital endometriosis and eutopic endometrium. Mold Med J. 2020; 63:6-11.
- Liu W, Lu X, Shi P, Yang G, Zhou Z, Li W, et al. TNF-alpha increases breast cancer stem-like cells through up-regulating TAZ expression via the non-canonical NF-kappaB pathway. Sci Rep. 2020; 10:1804. Publisher Full Text | DOI | PubMed [ PMC Free Article ]
- Jiang A, Li J, He Z, Liu Y, Qiao K, Fang Y, et al. Renal cancer: signaling pathways and advances in targeted therapies. MedComm (2020). 2024; 5:e676. Publisher Full Text | DOI | PubMed [ PMC Free Article ]
- Lopez de Andres J, Grinan-Lison C, Jimenez G, Marchal JA. Cancer stem cell secretome in the tumor microenvironment: a key point for an effective personalized cancer treatment. J Hematol Oncol. 2020; 13:136. Publisher Full Text | DOI | PubMed [ PMC Free Article ]
- Chen J, Zhu H, Chen S, Mi H. Apolipoprotein E is a Potential Biomarker for Predicting Cancer Prognosis and is Correlated with Immune Infiltration. Onco Targets Ther. 2024; 17:199-214. Publisher Full Text | DOI | PubMed [ PMC Free Article ]
- Krishnamohan M, Kaplanov I, Maudi-Boker S, Yousef M, Machluf-Katz N, Cohen I, et al. Tumor Cell-Associated IL-1alpha Affects Breast Cancer Progression and Metastasis in Mice through Manipulation of the Tumor Immune Microenvironment. Int J Mol Sci. 2024; 25Publisher Full Text | DOI | PubMed [ PMC Free Article ]
- Braicu C, Buse M, Busuioc C, Drula R, Gulei D, Raduly L, et al. A Comprehensive Review on MAPK: A Promising Therapeutic Target in Cancer. Cancers (Basel). 2019; 11Publisher Full Text | DOI | PubMed [ PMC Free Article ]
- Lee S, Rauch J, Kolch W. Targeting MAPK Signaling in Cancer: Mechanisms of Drug Resistance and Sensitivity. Int J Mol Sci. 2020; 21Publisher Full Text | DOI | PubMed [ PMC Free Article ]
- Bahar ME, Kim HJ, Kim DR. Targeting the RAS/RAF/MAPK pathway for cancer therapy: from mechanism to clinical studies. Signal Transduct Target Ther. 2023; 8:455. Publisher Full Text | DOI | PubMed [ PMC Free Article ]
- Acosta-Ramos E, Segovia-Mendoza M, Olivares-Reyes JA. Growth Factor Receptor Implications. Breast Cancer Treatment: An Interdisciplinary Approach. 2024; 7:237. DOI
- Baghel VS, Shinde S, Dixit V, et al. Dysregulated cell-signaling pathways in hepatocellular carcinoma: causes and therapeutic options. In: Theranostics and Precision Medicine for the Management of Hepatocellular Carcinoma, Volume 2. Amsterdam: Elsevier; 2022. DOI
- Karoulia Z, Gavathiotis E, Poulikakos PI. New perspectives for targeting RAF kinase in human cancer. Nat Rev Cancer. 2017; 17:676-91. Publisher Full Text | DOI | PubMed [ PMC Free Article ]
- Scardaci R, Berlinska E, Scaparone P, Vietti Michelina S, Garbo E, Novello S, et al. Novel RAF-directed approaches to overcome current clinical limits and block the RAS/RAF node. Mol Oncol. 2024; 18:1355-77. Publisher Full Text | DOI | PubMed [ PMC Free Article ]
- Tufail M, Wan WD, Jiang C, Li N. Targeting PI3K/AKT/mTOR signaling to overcome drug resistance in cancer. Chem Biol Interact. 2024; 396:111055. DOI | PubMed
- Lian S, Du Z, Chen Q, Xia Y, Miao X, Yu W, et al. From lab to clinic: The discovery and optimization journey of PI3K inhibitors. Eur J Med Chem. 2024; 277:116786. DOI | PubMed
- Wong AHN, Ma B, Lui RN. New developments in targeted therapy for metastatic colorectal cancer. Ther Adv Med Oncol. 2023; 15:17588359221148540. Publisher Full Text | DOI | PubMed [ PMC Free Article ]
- Qannita RA, Alalami AI, Harb AA, Aleidi SM, Taneera J, Abu-Gharbieh E, et al. Targeting Hypoxia-Inducible Factor-1 (HIF-1) in Cancer: Emerging Therapeutic Strategies and Pathway Regulation. Pharmaceuticals (Basel). 2024; 17Publisher Full Text | DOI | PubMed [ PMC Free Article ]
- Li Z, Wei R, Yao S, Meng F, Kong L. HIF-1A as a prognostic biomarker related to invasion, migration and immunosuppression of cervical cancer. Heliyon. 2024; 10:e24664. Publisher Full Text | DOI | PubMed [ PMC Free Article ]
- Sun H, Li Y, Zhang Y, Zhao X, Dong X, Guo Y, et al. The relevance between hypoxia-dependent spatial transcriptomics and the prognosis and efficacy of immunotherapy in claudin-low breast cancer. Front Immunol. 2022; 13:1042835. Publisher Full Text | DOI | PubMed [ PMC Free Article ]
- Dai H, Sheng X, Wang Y, Zhou L, Lin Y, Du Y, et al. HIF1alpha Regulates IL17 Signaling Pathway Influencing Sensitivity of Taxane-Based Chemotherapy for Breast Cancer. Front Cell Dev Biol. 2021; 9:729965. Publisher Full Text | DOI | PubMed [ PMC Free Article ]
- Chen X, Zhao J, Herjan T, Hong L, Liao Y, Liu C, et al. IL-17-induced HIF1alpha drives resistance to anti-PD-L1 via fibroblast-mediated immune exclusion. J Exp Med. 2022; 219Publisher Full Text | DOI | PubMed [ PMC Free Article ]
- Pavitra E, Kancharla J, Gupta VK, Prasad K, Sung JY, Kim J, et al. The role of NF-kappaB in breast cancer initiation, growth, metastasis, and resistance to chemotherapy. Biomed Pharmacother. 2023; 163:114822. DOI | PubMed
- Guo Q, Jin Y, Chen X, Ye X, Shen X, Lin M, et al. NF-kappaB in biology and targeted therapy: new insights and translational implications. Signal Transduct Target Ther. 2024; 9:53. Publisher Full Text | DOI | PubMed [ PMC Free Article ]
- Guo Q, Jin Y, Lin M, Zeng C, Zhang J. NF-kappaB signaling in therapy resistance of breast cancer: Mechanisms, approaches, and challenges. Life Sci. 2024; 348:122684. DOI | PubMed
- Cao Y, Yi Y, Han C, Shi B. NF-kappaB signaling pathway in tumor microenvironment. Front Immunol. 2024; 15:1476030. Publisher Full Text | DOI | PubMed [ PMC Free Article ]