DNA methylation profiling identifies two distinct subgroups in breast cancers with low hormone receptor expression, mainly associated with HER2 amplification status
Clinical Epigenetics volume 13, Article number: 184 (2021)
Current clinical guidelines suggest that breast cancers with low hormone receptor expression (LowHR) in 1–10% of tumor cells should be regarded as hormone receptor positive. However, clinical data show that these patients have worse outcome compared to patients with hormone receptor expression above 10%. We performed DNA methylation profiling on 23 LowHR breast cancer specimens, including 13 samples with HER2 amplification and compared our results with a reference breast cancer cohort from The Cancer Genome Atlas to clarify the status for this infrequent but important patient subgroup.
In unsupervised clustering and dimensionality reduction, breast cancers with low hormone receptor expression that lacked HER2 amplification usually clustered with triple negative breast cancer (TNBC) reference samples (8/10; “LowHR TNBC-like”). In contrast, most specimens with low hormone receptor expression and HER2 amplification grouped with hormone receptor positive cancers (11/13; “LowHR HRpos-like”). We observed highly similar DNA methylation patterns of LowHR TNBC-like samples and true TNBCs. Furthermore, the Ki67 proliferation index of LowHR TNBC-like samples and clinical outcome parameters were more similar to TNBCs and differed from LowHR HRpos-like cases.
We here demonstrate that LowHR breast cancer comprises two epigenetically distinct groups. Our data strongly suggest that LowHR TNBC-like samples are molecularly, histologically and clinically closely related to TNBC, while LowHR HRpos-like specimens are closely related to hormone receptor positive tumors.
The classification of breast cancer is based on the expression of the estrogen and progesterone receptor as well as the presence or absence of HER2 amplification . Patients with hormone receptor positive tumors usually benefit from endocrine therapy and have better disease specific survival compared to hormone receptor negative tumors . Additionally, tumors with HER2 amplification can be treated with different types of anti-HER2-therapy .
The evaluation of hormone receptor expression and the HER2 status is primarily performed by using immunohistochemistry (IHC) and in situ hybridization. In the 2010 revision of the guidelines of the American Society of Clinical Oncology and the College of American Pathologists (ASCO/CAP), the cut-off value for the definition of hormone receptor positive breast cancers was reduced from 10 to 1% . Only recently, the ASCO/CAP has updated this assessment, but also introduced a “ER low positive” category for samples with estrogen receptor expression between 1 and 10% . However, breast cancers with low hormone receptor expression are usually still considered as hormone receptor positive (= luminal) subtype . Various studies showed that therapy response and outcome of patients with tumors showing low hormone receptor expression between 1 and 10% (LowHR) are more similar to triple negative breast cancer (TNBC), the so-called basal subtype [6,7,8,9,10,11,12,13]. These results are supported by molecular studies, using RNA-based classification algorithms . This led to a dilemma in treatment decision, as it is unclear if patients with LowHR tumors should receive endocrine therapy or not or should be treated with more aggressive approaches like adjuvant chemotherapy, including platin compounds and even more, whether these tumors might be candidates for targeted therapeutic approaches like CDK4/6 inhibitors. More importantly, these patients are typically not included in clinical trials investigating new therapy strategies for TNBC, like immune checkpoint inhibitors. Furthermore, there are no data on the treatment of low hormone receptor expression and the additional HER2 amplification.
More recently, numerous studies showed that DNA methylation profiles are a reliable tool to classify different cancer types, reflecting the tumor’s cell-of-origin [15,16,17]. With regards to breast cancer, TNBCs exhibit global hypomethylation which is distinct from hormone receptor positive tumors . To find new molecular support for the currently proposed hypothesis that breast cancers with low hormone receptor expression resemble TNBCs, we subjected a cohort of these two subgroups to genome-wide DNA methylation profiling and compared the epigenetic profile of these specimens with a large reference cohort of hormone receptor positive tumors and TNBCs from The Cancer Genome Atlas (TCGA). Furthermore, we also correlated our results with clinical patient outcome.
Dimensionality reduction and cluster analysis
Unsupervised hierarchical clustering of LowHR and TCGA samples based on DNA methylation profiles revealed four distinct clusters (Fig. 1a), mainly representing normal breast tissue (“Normal Cluster”), hormone receptor positive breast cancers (“HR+ Cluster 1” and “HR+ Cluster 2”) as well as TNBCs (“TNBC Cluster”) with little overlap between the different groups. The two HR+ clusters were not associated with the luminal A or B subgroups as assessed using the RNA sequencing-based AIMS classification scheme. Similar results were seen using a t-SNE analysis (Fig. 1b–f), although there was no clear separation between the two groups of hormone receptor positive tumors that were observed in hierarchical clustering (Fig. 1f). Tumors in the periphery of their respective groups as well as tumor samples that fell in the normal breast tissue group tended to have relatively low tumor cell content (Additional file 1: Fig. S1).
Specimens with low hormone receptor expression either fell into one of the HR+ clusters or the TNBC cluster. This separation was mainly associated with the HER2 amplification status: samples that clustered with hormone receptor positive tumors (LowHR HRpos-like) usually showed HER2 amplification (11/13; 85%), while this was a rare event in LowHR TNBC-like samples (2/10; 20%). Of note, the HER2 status in copy number data derived from DNA methylation data was in line with the results from IHC and/or SISH. In particular, we confirmed the presence of HER2 amplification in the two HER2 positive LowHR TNBC-like samples and did not find any evidence for HER2 copy number changes in the two HER2 negative LowHR HRpos-like samples.
Summary copy number plots derived from DNA methylation data revealed additional similarities between LowHR TNBC-like samples and TNBCs, such as frequent losses of chromosome 4q and 5q as well as recurrent gains of chromosome 10p (Fig. 2).
Analysis of differential DNA methylation
To compare DNA methylation profiles between TNBCs and the two LowHR subtypes, we tested our dataset for differentially methylated positions (DMP) and differentially methylated regions (DMR).
A correlation analysis of the beta values of all CpGs revealed general hypermethylation in LowHR HRpos-like samples when compared to LowHR TNBC-like and TNBC specimens (Fig. 3a, b). There was a very strong, almost linear correlation between LowHR TNBC-like and TNBC samples (Fig. 3c). We identified 1446 DMPs between TNBC and LowHR HRpos-like samples. 1352 CpGs (93%) were hypomethylated in TNBCs (Fig. 3d). Comparing LowHR TNBC-like with LowHR HRpos-like samples, 1272 site were significantly differentially methylated. Again, the majority of CpGs (1151, 89%) were hypomethylated in LowHR TNBC-like samples (Fig. 3e). There was a considerable overlap between the identified DMPs, with 666 shared CpGs sites (Fig. 3g). Furthermore, we observed almost no difference between LowHR TNBC-like and TNBC samples, as our analysis only revealed a single differentially methylated probe (Fig. 3f). Similar results were obtained for DMRs. The DMRs identified when comparing LowHR HRpos-like samples and TNBCs included cancer-relevant genes such as EN1, TFF3 or IRX1 which have previously been described to be differentially methylated and expressed in TNBCs and hormone receptor positive breast cancers [26,27,28]. All DMRs are listed in Additional files 5–7: Tables S3–S5.
Correlation with proliferation index and clinical outcome
To further evaluate potential differences regarding the biological behavior and clinical outcome of the two identified LowHR subtypes and TNBCs, we compared the Ki67 tumor proliferation rate as well as pCR rate and OS between these groups using clinical trial information. DNA methylation classes did not correlate with local tumor stage or nodal status (Additional file 8: Table S6).
With 56% (p = 0.03) in TNBCs and 77% (p = 0.0002) in LowHR TNBC-like samples, the proliferation rate (Ki67) of these groups was significantly higher than in LowHR HRpos-like tumors (37%; Fig. 4a).
Furthermore, the pCR rate in TNBC (4 of 14 patients, 29%) and LowHR TNBC-like samples (3 of 10 patients, 30%) was lower than in LowHR HRpos-like cases (9 of 13 patients, 69%), although this difference had only borderline significance (p = 0.05 and p = 0.09; Fig. 4b). Notably, OS tended to be shorter in patients with LowHR TNBC-like tumors, although pairwise log rank test was not significant (p = 0.07) when compared with LowHR HRpos-like samples (Fig. 4c).
With this study, we demonstrate that breast cancers with low hormone receptor expression can be separated into two highly distinct groups by DNA methylation profiling. One group (“LowHR HRpos-like”) shows a DNA methylation profile similar to hormone receptor positive tumors, is strongly enriched with HER2 amplification, closely mirrors the chromosomal pattern of hormone receptor positive tumors and has a lower proliferation rate and a better clinical course. The other group (LowHR TNBC-like) shows a DNA methylation pattern of TNBC, has low rates of HER2 amplification, closely resembles the chromosomal pattern of TNBC and has a high proliferation rate and worse clinical outcome.
Since the lowering of the cut-off value for hormone receptor positive breast cancers by the ASCO/CAP in 2010, repeated concerns were raised that the outcome of patients with LowHR tumors is worse, as they do not benefit from endocrine therapy. However, despite from RNA expression profiles, there were no molecular data to support this clinical observation .
Using hierarchical clustering as well as t-SNE as a method of dimensionality reduction that has previously been proven to be a valuable tool to identify subgroups in DNA methylation data sets, we observed that normal breast tissue, hormone receptor positive breast cancers and TNBCs can be reliably distinguished based on their epigenetic profile. The small overlap between tumor samples and normal breast tissue was most likely due to low tumor cell content within these specimens, as we have recently demonstrated for other tumor entities . The overlap between the two hormone receptor positive and the TNBC cluster was also very low. This is in line with previous studies that demonstrated that TNBCs exhibit epigenetic signatures characterized by global DNA hypomethylation that are distinct from hormone receptor positive breast cancers . These findings provided the necessary rationale to compare the epigenetic signatures of these two subtypes with LowHR specimens. As described, LowHR tumors lacking HER2 amplification predominantly clustered with TNBCs, while specimens with additional HER2 amplification were typically very similar to hormone receptor positive tumors. Four specimens did not fit into this pattern. Differential DNA methylation analysis revealed almost no difference between LowHR TNBC-like and true TNBC samples, strongly supporting the assumption that these subtypes do not only show similar clinical outcome, but are also identical on the molecular level. Our investigation for differentially methylated genes between LowHR HRpos-like and TNBC specimens revealed a considerable overlap with cancer-relevant genes that have previously been described to be differentially methylated and expressed in hormone receptor positive tumors and TNBCs . This indicates that despite the low expression of hormone receptors, most HER2 amplified tumors are still very similar to hormone receptor positive breast cancers.
As DNA methylation signatures are considered to represent the cell-of-origin of a tumor, these findings could indicate that LowHR TNBC-like samples derive from the same cell type as TNBCs. On the other hand, the presence of HER2 amplification in LowHR HRpos-like specimens suggests a common origin with hormone receptor positive breast cancer, in line with previous cell-of-origin theories . This is also supported by the observation that all three hormone receptor negative samples with HER2 amplification that were included in the TCGA reference data set aggregated in one of the HR+ clusters.
Using genome wide copy number profiles derived from raw DNA methylation data, we were also able to study the differences in frequency of genome copy number alterations between hormone receptor positive tumors, TNBCs and LowHR specimens. In line with our results from DNA methylation profiling, LowHR TNBC-like tumors and TNBCs generally showed comparable numbers of recurrent genome copy abnormalities. On the other hand, copy number profiles of LowHR HRpos-like samples generally resembled hormone receptor positive breast cancers, with the exception of the previously described high rate of HER2 amplification. Previous studies already showed that there are differences in frequency of numeric alterations between TNBCs and hormone receptor positive tumors and the respective chromosomal locations we observed in our study were in line with these reports [30, 31]. These findings provide additional and independent proof for the correctness of the DNA methylation-based classification of LowHR tumors.
Based on our results, DNA methylation-based classification of LowHR tumors could be used to assess the individual prognosis and could guide treatment decisions in further studies. Based on the differences in their epigenetic profiles and clinical outcome data, one might suggest that LowHR TNBC-like samples require a similar treatment as TNBCs, whereas patients with LowHR HRpos-like tumors might still benefit from endocrine therapy. Strikingly, we observed high pCR rates for LowHR HRpos-like samples, which is unexpected due to their molecular similarity with hormone receptor positive breast cancer. However, it remains uncertain if this behavior is characteristic for this specific subgroup or if this observation is biased by other factors. Due to the retrospective design of our study and the small sample size, the clinical results should be interpreted with caution and need to be addressed in larger and prospective clinical trials.
In conclusion, using DNA methylation our study provides new molecular support that LowHR breast cancer comprises two molecularly distinct groups that can be separated by DNA methylation profiling, resembling TNBC and hormone receptor positive breast cancer with potential implications for individual patient prognosis and therapy selection.
Patients and samples
A flowchart which summarizes the composition of the study cohort is available as Additional file 2: Fig. S2.
For this study, we selected a total of 23 cases from the GeparSixto (NCT01426880) and GeparSepto (NCT01583426) trial of the German Breast Group (GBG) and the Arbeitsgemeinschaft Gynäkologische Onkologie – Breast Study Group (AGO-B) with estrogen and progesterone receptor expression between 1 and 10%, including 13 samples with HER2 amplification [19, 20]. Informed consent was obtained from all patients for study participation and translational research projects. Histopathological characteristics are summarized in Table 1. As a reference, we also included 14 randomly selected TNBC specimens from the GeparSepto study with estrogen and progesterone receptor expression < 1%. This study included only female early stage breast cancer patients without distant metastases. For all analyses, formalin-fixed and paraffin-embedded biopsy specimens from the initial diagnosis were used. Therefore, all patients were treatment-naïve. All samples were tested for hormone receptor expression and HER2 status as part of the central pathology assessment at the Institute of Pathology of the Charité – University Hospital Berlin.
Raw DNA methylation data from 709 samples as well as the corresponding pathology reports were obtained from the TCGA BRCA dataset available at the TCGA legacy archive (https://portal.gdc.cancer.gov/legacy-archive/). In this study, the definition of different breast cancer subtypes is based on the results from IHC. Therefore, we excluded patient samples without documented percentage of estrogen or progesterone receptor positive cells by IHC (n = 271) in the corresponding pathology reports as well specimens with bad quality in DNA methylation analysis (n = 16), as defined below. The resulting reference dataset consisted of 422 samples, including 334 breast cancer and 88 normal breast tissue specimens. The corresponding sample list is available as Additional file 3: Table S1.
IHC and silver in situ hybridization (SISH)
Antibodies, SISH reagents and scoring systems are listed in Additional file 4: Table S2. All antibodies for IHC were used in combination with the ultraView Universal DAB Detection Kit (VENTANA). Stainings were performed on the VENTANA BenchMark XT automated slide stainer using the “Cell Conditioning 1, Standard” setting for antigen retrieval and an incubation time of 16 min at 37 °C (estrogen and progesterone receptor) and 32 min (HER2) at 37 °C, respectively. HER2 SISH was also performed on the VENTANA BenchMark XT automated slide stainer with ISH-Protease 3 digestion for 8 min, followed by hybridization for 6 h and silver staining for 4 min. All slides were subsequently counterstained using the “Hematoxylin II” setting for 8 min and the “Bluing reagent” option for 4 min.
The results from central pathological investigation were used for further analyses.
We identified representative tumor areas using light microscopy of hematoxylin and eosin stained sections. If necessary, macrodissection was performed to reach a tumor cell content of at least 70%. Semi-automated DNA extraction was performed on the Maxwell RSC Instrument using the Maxwell RSC FFPE Plus DNA Purification Kit (Custom, AX4920; Promega), according to the manufacturer’s instructions. Extracted total DNA quantities were measured using the Qubit™ HS DNA Assay (Thermo Fisher Scientific).
DNA methylation analysis
We used the Illumina Infinium HD FFPE DNA Restore Kit for DNA restoration of FFPE samples. Following restoration, the EpiTect Bisulfite Kit (Qiagen) was used for bisulfite conversion. DNA methylation analysis was performed using the Illumina Infinium MethylationEPIC BeadChip, according to protocols supplied by the manufacturer.
Statistical analysis was performed using RStudio version 1.1.463 based on the statistical language R version 3.5.1. DNA methylation data were processed using the minfi package. The combineArrays function was used to merge data from the 450 K and 850 K array generation into a virtual 450 K dataset. The pfilter function of the wateRmelon package with the parameter perc = 5 was used to filter samples and probes with low quality . Samples were preprocessed using the normal-exponential out-of-band (Noob) normalization. Furthermore, we used the pwod function of the wateRmelon package for probe outlier detection. Testing for differentially methylation positions (DMP) and differentially methylation regions (DMR) was performed using the DMPfinder function of the minfi package and the bumphunter function of the bumphunter package, respectively . For both tests, hits with a false-discovery rate (FDR) below 0.05 and a median beta fold-change of at least 0.4 were considered to be significant. For the bumphunter function, we performed 1000 permutations.
T-distributed stochastic neighbor embedding (t-SNE) plots were generated based on the 5000 most variant CpG sites and using the Rtsne function of the Rtsne package with 2000 iterations while the perplexity was set to 25 . Heatmaps were generated using the ComplexHeatmap package with the “ward.D2” method and “euclidean” distance measuring . Tumor purity estimations for TCGA samples were derived from previously published datasets .
The means of Ki67 proliferation rates were tested for significance using the Wilcoxon–Mann–Whitney-Test. Pathological complete response (pCR) rates were compared using the Fisher’s test. Overall survival (OS) was visualized using Kaplan–Meier curves and tested for significant differences using pairwise logrank test.
Availability of data and materials
Raw DNA methylation data of the samples analyzed in this paper have been deposited in GEO (GSE163521). To review GEO accession GSE163521: Go to https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE163521. Enter token kvinqqgylnwdhyb into the box.
Goldhirsch A, Wood WC, Coates AS, Gelber RD, Thürlimann B, Senn H-J. Strategies for subtypes—dealing with the diversity of breast cancer: highlights of the St Gallen International Expert Consensus on the Primary Therapy of Early Breast Cancer 2011. Ann Oncol. 2011;22(8):1736–47.
Early Breast Cancer Trialists' Collaborative Group (EBCTCG). Effects of chemotherapy and hormonal therapy for early breast cancer on recurrence and 15-year survival: an overview of the randomised trials. Lancet. 2005;365(9472):1687–717.
Slamon DJ, Leyland-Jones B, Shak S, Fuchs H, Paton V, Bajamonde A, et al. Use of chemotherapy plus a monoclonal antibody against HER2 for metastatic breast cancer that overexpresses HER2. New Engl J Med. 2001;344(11):783–92.
Hammond MEH, Hayes DF, Dowsett M, Allred DC, Hagerty KL, Badve S, et al. American Society of Clinical Oncology/College of American Pathologists Guideline Recommendations for immunohistochemical testing of estrogen and progesterone receptors in breast cancer. J Clin Oncol. 2010;28(16):2784–95.
Allison KH, Hammond MEH, Dowsett M, McKernin SE, Carey LA, Fitzgibbons PL, et al. Estrogen and progesterone receptor testing in breast cancer: ASCO/CAP Guideline update. J Clin Oncol. 2020;38(12):1346–66.
Ding Y, Ding K, Yu K, Zou D, Yang H, He X, et al. Prognosis and endocrine therapy selection for patients with low hormone receptor-positive breast cancer following neoadjuvant chemotherapy: a retrospective study of 570 patients in China. Oncol Lett. 2019;18(6):6690–6.
Landmann A, Farrugia DJ, Zhu L, Diego EJ, Johnson RR, Soran A, et al. Low estrogen receptor (er)–positive breast cancer and neoadjuvant systemic chemotherapyis response similar to typical ER-positive or ER-negative disease? Am J Clin Pathol. 2018;150(1):34–42.
Raghav KPS, Hernandez-Aya LF, Lei X, Chavez-MacGregor M, Meric-Bernstam F, Buchholz TA, et al. Impact of low estrogen/progesterone receptor expression on survival outcomes in breast cancers previously classified as triple negative breast cancers. Cancer. 2012;118(6):1498–506.
Bouchard-Fortier A, Provencher L, Blanchette C, Diorio C. Prognostic and predictive value of low estrogen receptor expression in breast cancer. Curr Oncol. 2017;24(2):106–14.
Yoon MR, Rhu JY, Song BJ, Chae BJ, Yoo T-K. Treatment outcomes of weakly positive hormone receptor breast cancer and triple-negative breast cancer. J Breast Dis. 2019;7(1):1–8.
Fujii T, Kogawa T, Dong W, Sahin AA, Moulder S, Litton JK, et al. Revisiting the definition of estrogen receptor positivity in HER2-negative primary breast cancer. Ann Oncol. 2017;28(10):2420–8.
Prabhu JS, Korlimarla A, Desai K, Alexander A, Raghavan R, Coleegan EA, et al. A majority of low (1–10%) ER positive breast cancers behave like hormone receptor negative tumors. J Cancer. 2014;5(2):156–65.
Liao G-S, Dai M-S, Hsu H-M, Chu C-H, Hong Z-J, Fu C-Y, et al. Survival outcome of weak estrogen/progesterone receptor expression in HER2 negative breast cancer is similar to triple negative breast cancer. Eur J Surg Oncol Ejso. 2017;43(10):1855–61.
Iwamoto T, Booser D, Valero V, Murray JL, Koenig K, Esteva FJ, et al. Estrogen receptor (ER) mRNA and ER-related gene expression in breast cancers that are 1% to 10% ER-positive by immunohistochemistry. J Clin Oncol. 2012;30(7):729–34.
Bormann F, Rodríguez-Paredes M, Lasitschka F, Edelmann D, Musch T, Benner A, et al. Cell-of-origin DNA methylation signatures are maintained during colorectal carcinogenesis. Cell Rep. 2018;23(11):3407–18.
Capper D, Jones DTW, Sill M, Hovestadt V, Schrimpf D, Sturm D, et al. DNA methylation-based classification of central nervous system tumours. Nature. 2018;555(7697):469–74.
Jurmeister P, Bockmayr M, Seegerer P, Bockmayr T, Treue D, Montavon G, et al. Machine learning analysis of DNA methylation profiles distinguishes primary lung squamous cell carcinomas from head and neck metastases. Sci Transl Med. 2019;11(509):8513.
Koboldt DC, Fulton RS, McLellan MD, Schmidt H, Kalicki-Veizer J, McMichael JF, et al. Comprehensive molecular portraits of human breast tumours. Nature. 2012;490(7418):61–70.
von Minckwitz G, Schneeweiss A, Loibl S, Salat C, Denkert C, Rezai M, et al. Neoadjuvant carboplatin in patients with triple-negative and HER2-positive early breast cancer (GeparSixto; GBG 66): a randomised phase 2 trial. Lancet Oncol. 2014;15(7):747–56.
Untch M, Jackisch C, Schneeweiss A, Schmatloch S, Aktas B, Denkert C, et al. NAB-paclitaxel improves disease-free survival in early breast cancer: GBG 69–GeparSepto. J Clin Oncol. 2019;37(25):2226–34.
Pidsley R, Wong CCY, Volta M, Lunnon K, Mill J, Schalkwyk LC. A data-driven approach to preprocessing Illumina 450K methylation array data. BMC Genomics. 2013;14(1):293.
Jaffe AE, Murakami P, Lee H, Leek JT, Fallin MD, Feinberg AP, et al. Bump hunting to identify differentially methylated regions in epigenetic epidemiology studies. Int J Epidemiol. 2012;41(1):200–9.
Krijthe JH. Rtsne: T-distributed stochastic Neighbor embedding using a Barnes–Hut implementation. https://github.com/jkrijthe/Rtsne. undefined.
Gu Z, Eils R, Schlesner M. Complex heatmaps reveal patterns and correlations in multidimensional genomic data. Bioinform Oxf Engl. 2016;32(18):2847–9.
Aran D, Sirota M, Butte AJ. Systematic pan-cancer analysis of tumour purity. Nat Commun. 2015;6(1):8971.
Guo M, Sinha S, Wang SM. Coupled genome-wide DNA methylation and transcription analysis identified rich biomarkers and drug targets in triple-negative breast cancer. Cancers. 2019;11(11):1724.
Peluffo G, Subedee A, Harper NW, Kingston N, Jovanović B, Flores F, et al. EN1 is a transcriptional dependency in triple-negative breast cancer associated with brain metastasis. Cancer Res. 2019;79(16):4173–83.
Ahmed ARH, Griffiths AB, Tilby MT, Westley BR, May FEB. TFF3 Is a normal breast epithelial protein and is associated with differentiated phenotype in early breast cancer but predisposes to invasion and metastasis in advanced disease. Am J Pathol. 2012;180(3):904–16.
Prat A, Perou CM. Deconstructing the molecular portraits of breast cancer. Mol Oncol. 2011;5(1):5–23.
Neve RM, Chin K, Fridlyand J, Yeh J, Baehner FL, Fevr T, et al. A collection of breast cancer cell lines for the study of functionally distinct cancer subtypes. Cancer Cell. 2006;10(6):515–27.
Fridlyand J, Snijders AM, Ylstra B, Li H, Olshen A, Segraves R, et al. Breast tumor copy number aberration phenotypes and genomic instability. BMC Cancer. 2006;6(1):96.
We gratefully acknowledge the expert technical assistance by Peggy Wolkenstein, Daniel Teichmann, Alexandra Förster and Carola Geiler. We thank Damian Stichel for providing a modified version of the conumee package to generate frequency copy number profiles. The results shown here are, in part, based on data generated by the TCGA Research Network (http://cancergenome.nih.gov/).
This study was supported by the Deutsches Konsortium für Translationale Krebsforschung (DKTK partner Site Berlin and Joint Funding Upgrade “EPIC G8”) and the Deutsche Krebshilfe (“INTEGRATE-TN”). Philipp Jurmeister is participant in the BIH-Charité Digital Clinician Scientist Program funded by the Charité -Universitätsmedizin Berlin and the Berlin Institute of Health and the German Research Foundation (DFG).
Ethics approval and consent to participate
Ethics approval was granted by the local ethics committee of the Charité – Universitätsmedizin Berlin (EA1/139/05).
Ethics approval and consent to participate
Marion van Mackelenbergh reports personal fees from AstraZeneca, personal fees from Amgen, personal fees from Novartis, personal fees from Genomic Health, personal fees from Lilly, personal fees from Pfizer, personal fees from Roche, personal fees from Pierre Fabre, outside the submitted work. All other authors state that there are no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
: t-distributed stochastic neighbor embedding (t-SNE) plots showing the estimated tumor purity of The Cancer Genome Atlas samples based on manual estimation using light microscopy (a), RNAseq data (ESTIMATE; b), copy number variations (ABSOLUTE; c), DNA methylation (LUMP; d) and a combined score (CPE; e).
: Flowchart showing the composition of the study cohort with samples from the publicly available “The Cancer Genome Atlas” (TCGA) dataset as well as specimens from the clinical GeparSixto and GeparSepto trial. Abbreviations: BRCA = breast cancer; ER = estrogen receptor; PR = progesterone receptor; TNBC = triple negative breast cancer; t-SNE = t-distributed stochastic neighbor embedding.
: Sample annotation list for samples from the TCGA dataset.
: Table displaying the antibodies used for immunohistochemistry and the corresponding scoring systems.
: Results from DMR analysis comparing LowHR HRpos-like and TNBC samples
: Results from DMR analysis comparing LowHR TNBC-like and TNBC samples
: Results from DMR analysis comparing LowHR HRpos-like and LowHR TNBC-like samples
: Correlation of clinicopathological metadata with DNA methylation groups
About this article
Cite this article
Jurmeister, P., Weber, K., Villegas, S. et al. DNA methylation profiling identifies two distinct subgroups in breast cancers with low hormone receptor expression, mainly associated with HER2 amplification status. Clin Epigenet 13, 184 (2021). https://doi.org/10.1186/s13148-021-01176-5