Dataset Information

Tumor microenvironment-related multigene prognostic prediction model for breast cancer.

ABSTRACT:

Background

Breast cancer is an invasive disease with complex molecular mechanisms. Prognosis-related biomarkers are still urgently needed to predict outcomes of breast cancer patients.

Methods

Original data were download from The Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO). The analyses were performed using perl-5.32 and R-x64-4.1.1.

Results

In this study, 1086 differentially expressed genes (DEGs) were identified in the TCGA cohort; 523 shared DEGs were identified in the TCGA and GSE10886 cohorts. Eight subtypes were estimated using non-negative matrix factorization clustering with significant differences seen in overall survival (OS) and progression-free survival (PFS) (P < 0.01). Univariate Cox analysis and least absolute shrinkage and selection operator (LASSO) regression analysis were performed to develop a related risk score related to the 17 DEGs; this score separated breast cancer into low- and high-risk groups with significant differences in survival (P < 0.01) and showed powerful effectiveness (TCGA all group: 1-year area under the curve [AUC] = 0.729, 3-year AUC = 0.778, 5-year AUC = 0.781). A nomogram prediction model was constructed using non-negative matrix factorization clustering, the risk score, and clinical characteristics. Our model was confirmed to be related with tumor microenvironment. Furthermore, DEGs in high-risk breast cancer were enriched in histidine metabolism (normalized enrichment score [NES] = 1.49, P < 0.05), protein export (NES = 1.58, P < 0.05), and steroid hormone biosynthesis signaling pathways (NES = 1.56, P < 0.05).

Conclusions

We established a comprehensive model that can predict prognosis and guide treatment.

SUBMITTER: Hong K

PROVIDER: S-EPMC8833129 | biostudies-literature | 2022 Jan

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Tumor microenvironment-related multigene prognostic prediction model for breast cancer.

Hong Kai K Zhang Yingjue Y Yao Lingli L Zhang Jiabo J Sheng Xianneng X Guo Yu Y

Aging 20220120 2

<h4>Background</h4>Breast cancer is an invasive disease with complex molecular mechanisms. Prognosis-related biomarkers are still urgently needed to predict outcomes of breast cancer patients.<h4>Methods</h4>Original data were download from The Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO). The analyses were performed using perl-5.32 and R-x64-4.1.1.<h4>Results</h4>In this study, 1086 differentially expressed genes (DEGs) were identified in the TCGA cohort; 523 shared DEGs wer ...[more]

PMID: 35060926

Similar Datasets

Project description:BackgroundTumor microenvironment (TME) status is closely related to breast cancer (BC) prognosis and systemic therapeutic effects. However, to date studies have not considered the interactions of immune and stromal cells at the gene expression level in BC as a whole. Herein, we constructed a predictive model, for adjuvant decision-making, by mining TME molecular expression information related to BC patient prognosis and drug treatment sensitivity.MethodsClinical information and gene expression profiles were extracted from The Cancer Genome Atlas (TCGA), with patients divided into high- and low-score groups according to immune/stromal scores. TME-related prognostic genes were identified using Kaplan-Meier analysis, functional enrichment analysis, and protein-protein interaction (PPI) networks, and validated in the Gene Expression Omnibus (GEO) database. Least absolute shrinkage and selection operator (LASSO) Cox regression analysis was used to construct and verify a prognostic model based on TME-related genes. In addition, the patients' response to chemotherapy and immunotherapy was assessed by survival outcome and immunohistochemistry (IPS). Immunohistochemistry (IHC) staining laid a solid foundation for exploring the value of novel therapeutic target genes.ResultsBy dividing patients into low- and high-risk groups, a significant distinction in overall survival was found (p < 0.05). The risk model was independent of multiple clinicopathological parameters and accurately predicted prognosis in BC patients (p < 0.05). The nomogram-integrated risk score had high prediction accuracy and applicability, when compared with simple clinicopathological features. As predicted by the risk model, regardless of the chemotherapy regimen, the survival advantage of the low-risk group was evident in those patients receiving chemotherapy (p < 0.05). However, in patients receiving anthracycline (A) therapy, outcomes were not significantly different when compared with those receiving no-A therapy (p = 0.24), suggesting these patients may omit from A-containing adjuvant chemotherapy. Our risk model also effectively predicted tumor mutation burden (TMB) and immunotherapy efficacy in BC patients (p < 0.05).ConclusionThe prognostic score model based on TME-related genes effectively predicted prognosis and chemotherapy effects in BC patients. The model provides a theoretical basis for novel driver-gene discover in BC and guides the decision-making for the adjuvant treatment of early breast cancer (eBC).

Project description:Background: The tumor microenvironment (TME) has been reported to have significant value in the diagnosis and prognosis of cancers. This study aimed to identify key biomarkers in the TME of luminal breast cancer (BC). Methods: We obtained immune scores (ISs) and stromal scores (SSs) for The Cancer Genome Atlas (TCGA) luminal BC cohort from the online ESTIMATE (Estimation of STromal and Immune cells in MAlignant Tumor tissues using Expression data) portal. The relationships between ISs and SSs and the overall survival of luminal BC patients were assessed by the Kaplan-Meier method. The differentially expressed messenger RNAs (DEmRNAs) related to the ISs and SSs were subjected to functional enrichment analysis. Additionally, a competing endogenous RNA (ceRNA) network was constructed with differentially expressed microRNAs (DEmiRNAs) and long noncoding RNAs (DElncRNAs). Furthermore, a protein-protein interaction (PPI) network was established to analyze the DEmRNAs in the ceRNA network. Then, survival analysis of biomarkers involved in the ceRNA network was carried out to explore their prognostic value. Finally, these biomarkers were validated using the luminal BC dataset from the Gene Expression Omnibus (GEO) database. Results: The results showed that ISs were significantly associated with longer survival times of luminal BC patients. Functional enrichment analysis showed that the DEmRNAs were mainly associated with immune response, antigen binding, and the extracellular region. In the PPI network, the top 10 DEmRNAs were identified as hub genes that affected the TME of luminal BC. Finally, two DEmiRNAs, two DElncRNAs, and 17 DEmRNAs of the ceRNA network associated with the TME were shown to have prognostic value. Subsequently, the expression of 15 prognostic biomarkers was validated in one additional dataset (GSE81002). In particular, one lncRNA (GVINP1) and five mRNAs (CCDC69, DOCK2, IKZF1, JCHAIN, and NCKAP1L) were novel biomarkers. Conclusions: Our studies demonstrated that ISs were associated with the survival of luminal BC patients, and a set of novel biomarkers that might play a prognostic role in the TME of luminal BC was identified.

Project description:The outcomes of some breast cancer patients remain poor due to being susceptible to recurrence, metastasis and drug resistance, and lactate metabolism has been described as a hallmark of cancer and a contributor to cancer progression and immune escape. Hence, it is worthy of seeking potentially novel biomarkers from lactate metabolism relevant perspectives for this particular cohort of patients. In this context, 205 available lactate metabolism-related genes (LMGs) were obtained by a search of multiple genesets, and the landscape of somatic mutation, copy number variation, and mRNA expression levels was investigated among these genes. Crucially, 9 overall survival-related LMGs were identified through univariate Cox regression analysis in The Cancer Genome Atlas (TCGA) and Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) databases. Subsequently, a prognostic signature, defined as Lactate Metabolism Index (LMI), was established with 5 OS-related LMGs using Least Absolute Shrinkage and Selection Operator (LASSO) Cox hazard regression analysis in TCGA training set, and then validated in two external cohorts (METABRIC and GSE96058). From the comprehensive results, breast cancer patients with high LMI had considerably poorer survival probability across all cohorts, and the degree of clinical features tended to be more severe as the LMI value increased. Furthermore, a prognostic nomogram incorporating LMI, age, and AJCC stage was constructed and demonstrated great prediction performance for OS of breast cancer patients, which was evaluated by the calibration plot and the decision curve analysis. Moreover, the potential effect of different LMI values on levels of immune checkpoints, tumor-infiltrating immune cells, and cytokines were explored ultimately, and patients with higher LMI values might gain an immunosuppressive tumor microenvironment that contributed to immune escape of breast cancer and inferior prognosis. Collectively, all findings in the study indicated the potential prognostic value of LMI in breast cancer, providing further implications for the role of lactate metabolism in breast cancer prognosis, tumor immune microenvironment, and immunotherapy.

Project description:BackgroundBreast cancer is one of the most frequently diagnosed cancers among women worldwide. Alterations in the tumor microenvironment (TME) have been increasingly recognized as key in the development and progression of breast cancer in recent years. To deeply comprehend the gene expression profiling of the TME and identify immunological targets, as well as determine the relationship between gene expression and different prognoses is highly critical.MethodsThe stromal/immune scores of breast cancer patients from The Cancer Genome Atlas (TCGA) were employed to comprehensively evaluate the TME. Then, TME characteristics were assessed, overlapping genes of the top 3 Gene Ontology (GO) terms and upregulated differentially expressed genes (DEGs) were analyzed. Finally, through combined analyses of overall survival, time-dependent receiver operating characteristic (ROC), and protein-protein interaction (PPI) network, novel immune related genes with good prognosis were screened and validated in both TCGA and GEO database.ResultsAlthough the TME did not correlate with the stages of breast cancer, it was closely associated with the subtypes of breast cancer and gene mutations (CDH1, TP53 and PTEN), and had immunological characteristics. Based on GO functional enrichment analysis, the upregulated genes from the high vs low immune score groups were mainly involved in T cell activation, the external side of the plasma membrane, and receptor ligand activity. The top GO terms of the upregulated DEGs from the high vs low immune score groups exhibited better prognosis in breast cancer; 15 of them were related to good prognosis in breast cancer, especially CD226 and KLRC4-KLRK1.ConclusionsHigh CD226 and KLRC4-KLRK1 expression levels were identified and validated to correlate with better overall survival in specific stages or subtypes of breast cancer. CD226, KLRC4-KLRK1 and other new targets seem to be promising avenues for promoting antitumor targeted immunotherapy in breast cancer.

Project description:BRAF mutation is a representative oncogenic mutation, with a frequency of 60% in papillary thyroid carcinoma (PTC), but the reasons for the poor prognosis and more aggressive course of BRAF-mutated PTC are controversial. Tumor immune microenvironment (TIME) is an essential factor permitting the development and progression of malignancy, but whether TIME participates in the prognosis of BRAF-mutated PTC has not yet been reported. The primary goal of the present study was to provide a comprehensive TIME-related prognostic model to increase the predictive accuracy of progression-free survival (PFS) in patients with BRAF-mutated PTC. In this study, we analyzed the mRNA-seq data and corresponding clinical data of PTC patients obtained from the TCGA database. By calculating the TIME scores (immune score, stromal score and ESTIMATE score), the BRAF mutation group (n=237) was dichotomized into the high- and low-score groups. By functional analysis of differentially expressed genes (DEGs) in different high/low score groups, we identified 2 key TIME-related genes, HTR3A and NIPAL4, which affected PFS in BRAF-mutated PTC. A risk scoring system was developed by multivariate Cox analysis based on the abovementioned 2 TIME-related genes. Then, the BRAF-mutated cohort was divided into the high- and low-risk groups using the median risk score as a cutoff. A high risk score correlated positively with a higher HTR3A/NIPAL4 expression level but negatively with PFS in BRAF-mutated PTC. Ultimately, a nomogram was constructed by combining risk score with clinical parameter (Tumor stage), and the areas under the ROC curve (AUCs) of the nomogram for predicting 1-, 3- and 5-year PFS were then calculated and found to be 0.694, 0.707 and 0.738, respectively, indicating the improved accuracy and clinical utility of the nomogram versus the risk score model in the BRAF-mutated PTC cohort. Moreover, we determined the associations between prognostic genes or risk score and immune cell infiltration by two-way ANOVA. In the high-risk score, high HTR3A expression, and high NIPAL4 expression groups, higher infiltration of immune cells was found. Collectively, these findings confirm that the nomogram is effective in predicting the outcome of BRAF-mutated PTC and will add a spatial dimension to the developing risk stratification system.

Project description:Breast cancer is one of the most common female malignancies worldwide. Due to its early metastases formation and a high degree of malignancy, the 10 year-survival rate of metastatic breast cancer does not exceed 30%. Thus, more precise biomarkers are urgently needed. In our study, we first estimated the tumor microenvironment (TME) infiltration using the xCell algorithm. Based on TME infiltration, the three main TME clusters were identified using consensus clustering. Our results showed that the three main TME clusters cause significant differences in survival rates and TME infiltration patterns (log-rank test, p = 0.006). Then, multiple machine learning algorithms were used to develop a nine-pathway-based TME-related risk model to predict the prognosis of breast cancer (BRCA) patients (the immune-related pathway-based risk score, defined as IPRS). Based on the IPRS, BRCA patients were divided into two subgroups, and patients in the IPRS-low group presented significantly better overall survival (OS) rates than the IPRS-high group (log-rank test, p < 0.0001). Correlation analysis revealed that the IPRS-low group was characterized by increases in immune-related scores (cytolytic activity (CYT), major histocompatibility complex (MHC), T cell-inflamed immune gene expression profile (GEP), ESTIMATE, immune, and stromal scores) while exhibiting decreases in tumor purity, suggesting IPRS-low patients may have a strong immune response. Additionally, the gene-set enrichment analysis (GSEA) result confirmed that the IPRS-low patients were significantly enriched in several immune-associated signaling pathways. Furthermore, multivariate Cox analysis revealed that the IPRS was an independent prognostic biomarker after adjustment by clinicopathologic characteristics. The prognostic value of the IPRS model was further validated in three external validation cohorts. Altogether, our findings demonstrated that the IPRS was a powerful predictor to screen out certain populations with better prognosis in breast cancer and may serve as a potential biomarker guiding clinical treatment decisions.

Dataset Information

Tumor microenvironment-related multigene prognostic prediction model for breast cancer.

Background

Methods

Results

Conclusions

Publications

Tumor microenvironment-related multigene prognostic prediction model for breast cancer.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets