Dataset Information

A novel immune subtype classification of ER-positive, PR-negative and HER2-negative breast cancer based on the genomic and transcriptomic landscape.

ABSTRACT:

Background

The diversity and plasticity behind ER+/PR-/HER2- breast cancer have not been widely explored. It is essential to identify heterogeneous microenvironment phenotypes and investigate specific genomic events driving the formation of these phenotypes.

Methods

Based on the immune-related gene expression profiles of 411 ER+/PR-/HER2- breast cancers in the METABRIC cohort, we used consensus clustering to identify heterogeneous immune subtypes and assessed their reproducibility in an independent meta-cohort including 135 patients collected from GEO database. We further analyzed the differences of cellular and molecular characteristics, and potential immune escape mechanism among immune subtypes. In addition, we constructed a transcriptional trajectory to visualize the distribution of individual patient.

Results

Our analysis identified and validated five reproducible immune subtypes with distinct cellular and molecular characteristics, potential immune escape mechanisms, genomic drivers, as well as clinical outcomes. An immune-cold subtype, with the least amount of lymphocyte infiltration, had a poorer prognosis. By contrast, an immune-hot subtype, which demonstrated the highest infiltration of CD8+ T cells, DCs and NK cells, and elevated IFN-γ response, had a comparatively favorable prognosis. Other subtypes showed more diverse gene expression and immune infiltration patterns with distinct clinical outcomes. Finally, our analysis revealed a complex immune landscape consisting of both discrete cluster and continuous spectrum.

Conclusion

Overall, this study revealed five heterogeneous immune subtypes among ER+/PR-/HER2- breast cancer, also provided important implications for clinical translations.

SUBMITTER: Xie P

PROVIDER: S-EPMC8454077 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:While recent integrative analyses of copy number and gene expression data in breast cancer have revealed a complex molecular landscape with multiple subtypes and many oncogenic/tumour suppressor driver events, much less is known about the role of DNA methylation in shaping breast cancer taxonomy and defining driver events.Here, we applied a powerful integrative network algorithm to matched DNA methylation and RNA-Seq data for 724 estrogen receptor (ER)-positive (ER+) breast cancers and 111 normal adjacent tissue specimens from The Cancer Genome Atlas (TCGA) project, in order to identify putative epigenetic driver events and to explore the resulting molecular taxonomy. This revealed the existence of nine functionally deregulated epigenetic hotspots encompassing a total of 146 genes, which we were able to validate in independent data sets encompassing over 1000 ER+ breast cancers. Integrative clustering of the matched messenger RNA (mRNA) and DNA methylation data over these genes resulted in only two clusters, which correlated very strongly with the luminal-A and luminal B subtypes. Overall, luminal-A and luminal-B breast cancers shared the same epigenetically deregulated hotspots but with luminal-B cancers exhibiting increased aberrant DNA methylation patterns relative to normal tissue. We show that increased levels of DNA methylation and mRNA expression deviation from the normal state define a marker of poor prognosis. Our data further implicates epigenetic silencing of WNT signalling antagonists and bone morphogenetic proteins (BMP) as key events underlying both luminal subtypes but specially of luminal-B breast cancer. Finally, we show that DNA methylation changes within the identified epigenetic interactome hotspots do not exhibit mutually exclusive patterns within the same cancer sample, instead exhibiting coordinated changes within the sample.Our results indicate that the integrative DNA methylation and transcriptomic landscape of ER+ breast cancer is surprisingly homogeneous, defining two main subtypes which strongly correlate with luminal-A/B subtype status. In particular, we identify WNT and BMP signalling as key epigenetically deregulated tumour suppressor pathways in luminal ER+ breast cancer, with increased deregulation seen in luminal-B breast cancer.

Project description:Estrogen receptor-positive, progesterone receptor-negative, and human epidermal growth factor receptor 2 (HER2)-negative (ER+PR-HER2-) breast cancer comprise a special type of breast cancer that constitutes ~10% of all breast cancer patients. ER+PR-HER2- tumor benefits less from endocrine therapy, while its genomic features remain elusive. In this study, we systematically assessed the multiomic landscape and endocrine responsiveness of ER+PR-HER2- breast cancer. Methods: This study incorporated five cohorts. The first and second cohorts were from the Surveillance, Epidemiology, and End Results database (n=130,856) and Molecular Taxonomy of Breast Cancer International Consortium (n=1,055) for analyzing survival outcomes and endocrine responsiveness. The third cohort was from The Cancer Genome Atlas (n=630) for multiomic analysis and endocrine-resistant subgroup exploration. The fourth cohort, from the MD Anderson database (n=92), was employed to assist gene selection. The fifth cohort was a prospective observational cohort from Fudan University Shanghai Cancer Center (n=245) that was utilized to validate the gene-defined subgroup by immunohistochemistry (IHC). Results: Clinically, ER+PR-HER2- tumors showed lower endocrine responsiveness than did ER+PR+HER2- tumors. Genomically, copy number loss or promoter methylation of PR genes occurred in 75% of ER+PR-HER2- tumors, collectively explaining PR loss. ER+PR-HER2- tumors had higher TP53 (30.3% vs. 17.0%) and lower PIK3CA mutation rates (25.8% vs. 42.7%) and exhibited more ZNF703 (21.5% vs. 13.6%) and RPS6KB1 (18.5% vs. 7.8%) amplification events than ER+PR+HER2- tumors. Among ER+PR-HER2- tumors, nearly 20% were of the PAM50-defined non-luminal-like subgroup and manifested lower endocrine sensitivity scores and enriched biosynthesis, metabolism and DNA replication pathways. We further identified the non-luminal-like subgroup using three IHC markers, GATA3, CK5, and EGFR. These IHC-defined non-luminal-like (GATA3-negative, CK5-positive and/or EGFR-positive) tumors received limited benefit from adjuvant endocrine therapy. Conclusion: ER+PR-HER2- breast cancer consists of clinically and genomically distinct groups that may require different treatment strategies. The non-luminal-like subgroup was associated with reduced benefit from endocrine therapy.

Project description:Background:Treatment strategies for various subtypes of breast cancer (BC) are different based on their distinct molecular characteristics. Therefore, it is very important to identify key differentially expressed genes (DEGs) between ER-positive/HER2-negative BC and ER-negative/HER2-negative BC. Methods:Gene expression profiles of GSE22093 and GSE23988 were obtained from the Gene Expression Omnibus database. There were 74 ER-positive/HER2-negative BC tissues and 85 ER-negative/HER2-negative BC tissues in the two profile datasets. DEGs between ER-positive/HER2-negative tissues and ER-negative/HER2-negative BC tissues were identified by the GEO2R tool. The common DEGs among the two datasets were detected with Venn software online. Next, we made use of the Database for Annotation, Visualization and Integrated Discovery to analyze enriched Kyoto Encyclopedia of Gene and Genome (KEGG) pathways and gene ontology terms. Then, the protein-protein interactions (PPIs) of these DEGs were visualized by Cytoscape with the Search Tool for the Retrieval of Interacting Genes. Of the proteins in the PPI network, Molecular Complex Detection plug-in analysis identified nine core upregulated genes and one core downregulated gene. UALCAN and Gene Expression Profiling Interactive Analysis were applied to determine the expression of these 10 genes in BC. Furthermore, for the analysis of overall survival among those genes, the Kaplan-Meier method was implemented. Results:Ninety-three common DEGs (63 upregulated and 30 downregulated) were identified. KEGG pathway enrichment analysis showed that upregulated DEGs were particularly enriched in the progesterone-mediated oocyte maturation pathway. In addition, PGR might be a prognostic biomarker for ER-positive/HER2-negative BC. CCND1 is a poor prognostic biomarker for ER-positive/HER2-negative BC and ER-negative/HER2-negative BC. Moreover, TFF1, AGR2 and EGFR might be predictive biomarkers of node metastasis in ER-positive/HER2-negative BC and ER-negative/HER2-negative BC. Conclusions:CCND1, AGR2, PGR, TFF1 and EGFR are the key DEGs between ER-positive/HER2-negative BC and ER-negative/HER2-negative BC. Further studies are required to confirm the functions of the identified genes.