Clinical associations of ESR2 (estrogen receptor beta) expression across thousands of primary breast tumors

Hina Dalal, Malin Dahlgren, Sergii Gladchuk, Christian Brueffer, Sofia K Gruvberger-Saal, Lao H Saal, Hina Dalal, Malin Dahlgren, Sergii Gladchuk, Christian Brueffer, Sofia K Gruvberger-Saal, Lao H Saal

Abstract

Estrogen receptor alpha (ERα, encoded by ESR1) is a well-characterized transcription factor expressed in more than 75% of breast tumors and is the key biomarker to direct endocrine therapies. On the other hand, much less is known about estrogen receptor beta (ERβ, encoded by ESR2) and its importance in cancer. Previous studies had some disagreement, however most reports suggested a more favorable prognosis for patients with high ESR2 expression. To add further clarity to ESR2 in breast cancer, we interrogated a large population-based cohort of primary breast tumors (n = 3207) from the SCAN-B study. RNA-seq shows ESR2 is expressed at low levels overall with a slight inverse correlation to ESR1 expression (Spearman R = -0.18, p = 2.2e-16), and highest ESR2 expression in the basal- and normal-like PAM50 subtypes. ESR2-high tumors had favorable overall survival (p = 0.006), particularly in subgroups receiving endocrine therapy (p = 0.03) and in triple-negative breast cancer (p = 0.01). These results were generally robust in multivariable analyses accounting for patient age, tumor size, node status, and grade. Gene modules consistent with immune response were associated to ESR2-high tumors. Taken together, our results indicate that ESR2 is generally expressed at low levels in breast cancer but associated with improved overall survival and may be related to immune response modulation.

Conflict of interest statement

The authors declare no competing interests.

© 2022. The Author(s).

Figures

Figure 1
Figure 1
ESR1 and ESR2 mRNA expression in SCAN-B and TCGA data sets. (A) Scatterplot of ESR1 and ESR2 mRNA expression (log2TPM). ESR1 and ESR2 data points are colored by density using 2D kernel density estimation function from the MASS R package. Adjacent to the scatterplot, histograms are shown indicating the frequency of expression values and color-coded according to ERα clinical status. (B) Expression of ESR1 and ESR2 is shown for the TCGA dataset as in panel A.
Figure 2
Figure 2
ESR2 mRNA by PAM50 molecular subtype in SCAN-B (A) and TCGA (B) cohorts, and by age groups (age at diagnosis) in SCAN-B (C) and TCGA (D).
Figure 3
Figure 3
ESR2 expression and association to overall survival (OS) and relapse-free interval (RFI) in the full SCAN-B cohort (A,B), the endocrine-treated sub-group (C,D) and the chemotherapy-treated subgroup (E,F).
Figure 4
Figure 4
ESR2 expression and association to overall survival (OS) and relapse-free interval (RFI) in the SCAN-B clinical groups. (A,B) Patients with ERα-positive, HER2-negative breast cancer; (C,D) ERα-positive, HER2-positive breast cancer; (E,F) Triple-negative breast cancer (TNBC); and (G,H) ERα-negative, HER2-positive breast cancer.
Figure 5
Figure 5
Multivariate analysis of high ESR2 expression in the full SCAN-B cohort, the endocrine-treated group, and the triple-negative breast cancer (TNBC) group.
Figure 6
Figure 6
Gene Set Enrichment Analysis (GSEA; GO category: Biological Process) of genes ranked by fold change (log2FC) and p-value < 0.05, associated with ESR2-high vs -low in SCAN-B. GSEA analysis based on ESR2-high vs -low was performed separately for the ERα-positive (A) and ERα-negative (B) subgroups. Categories found enriched in both subgroup analyses are indicated by red text.

References

    1. Bray F, et al. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin. 2018;68:394–424.
    1. Sung H, et al. Global Cancer Statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin. 2021;71:209–249.
    1. Green S, et al. Human oestrogen receptor cDNA: Sequence, expression and homology to v-erb-A. Nature. 1986;320:134–139.
    1. Kuiper GG, Enmark E, Pelto-Huikko M, Nilsson S, Gustafsson JA. Cloning of a novel receptor expressed in rat prostate and ovary. Proc. Natl. Acad. Sci. U. S. A. 1996;93:5925–5930.
    1. Jensen EV. On the mechanism of estrogen action. Perspect. Biol. Med. 1962;6:47–59.
    1. Jensen EV, Block GE, Smith S, Kyser K, DeSombre ER. Estrogen receptors and breast cancer response to adrenalectomy. Natl. Cancer Inst. Monogr. 1971;34:55–70.
    1. Bundred NJ. Prognostic and predictive factors in breast cancer. Cancer Treat. Rev. 2001;27:137–142.
    1. Cleator S, Ashworth A. Molecular profiling of breast cancer: Clinical implications. Br. J. Cancer. 2004;90:1120–1124.
    1. Park YR, et al. Absence of estrogen receptor is associated with worse oncologic outcome in patients who were received neoadjuvant chemotherapy for breast cancer. Asian J. Surg. 2020;43:467–475.
    1. Kumar M, Salem K, Tevaarwerk AJ, Strigel RM, Fowler AM. Recent advances in imaging steroid hormone receptors in breast cancer. J. Nucl. Med. 2020;61:172–176.
    1. van Kruchten M, et al. PET imaging of estrogen receptors as a diagnostic tool for breast cancer patients presenting with a clinical dilemma. J. Nucl. Med. 2012;53:182–190.
    1. Oueslati M, et al. ERα and ERβ co-expression: An indicator of aggressive tumors and hormonal sensitivity. Oncol. Lett. 2017;14:1675–1682.
    1. Mirza S, et al. Clinical significance of promoter hypermethylation of ERβ and RARβ2 in tumor and serum DNA in Indian breast cancer patients. Ann. Surg. Oncol. 2012;19:3107–3115.
    1. Elebro K, et al. High estrogen receptor β expression Is prognostic among adjuvant chemotherapy-treated patients-results from a population-based breast cancer cohort. Clin. Cancer Res. 2017;23:766–777.
    1. Girault I, et al. Altered expression pattern of alternatively spliced estrogen receptor beta transcripts in breast carcinoma. Cancer Lett. 2004;215:101–112.
    1. Zhou Y, Liu X. The role of estrogen receptor beta in breast cancer. Biomark. Res. 2020;8:39.
    1. Crandall DL, Busler DE, Novak TJ, Weber RV, Kral JG. Identification of estrogen receptor beta RNA in human breast and abdominal subcutaneous adipose tissue. Biochem. Biophys. Res. Commun. 1998;248:523–526.
    1. Kuiper GG, Shughrue PJ, Merchenthaler I, Gustafsson JA. The estrogen receptor beta subtype: A novel mediator of estrogen action in neuroendocrine systems. Front. Neuroendocrinol. 1998;19:253–286.
    1. Hawse JR, et al. Optimized immunohistochemical detection of estrogen receptor beta using two validated monoclonal antibodies confirms its expression in normal and malignant breast tissues. Breast Cancer Res. Treat. 2020;179:241–249.
    1. Omoto Y, Iwase H. Clinical significance of estrogen receptor β in breast and prostate cancer from biological aspects. Cancer Sci. 2015;106:337–343.
    1. Pons DG, et al. The presence of estrogen receptor β modulates the response of breast cancer cells to therapeutic agents. Int. J. Biochem. Cell Biol. 2015;66:85–94.
    1. Razandi M, Pedram A, Jordan VC, Fuqua S, Levin ER. Tamoxifen regulates cell fate through mitochondrial estrogen receptor beta in breast cancer. Oncogene. 2013;32:3274–3285.
    1. Ruddy SC, et al. Preferential estrogen receptor β ligands reduce Bcl-2 expression in hormone-resistant breast cancer cells to increase autophagy. Mol. Cancer Ther. 2014;13:1882–1893.
    1. Song P, et al. Estrogen receptor β inhibits breast cancer cells migration and invasion through CLDN6-mediated autophagy. J. Exp. Clin. Cancer Res. CR. 2019;38:354.
    1. Greish K, Nehoff H, Bahman F, Pritchard T, Taurin S. Raloxifene nano-micelles effect on triple-negative breast cancer is mediated through estrogen receptor-β and epidermal growth factor receptor. J. Drug Target. 2019;27:903–916.
    1. Tu Z, et al. The enhanced antiproliferative response to combined treatment of trichostatin A with raloxifene in MCF-7 breast cancer cells and its relevance to estrogen receptor β expression. Mol. Cell. Biochem. 2012;366:111–122.
    1. Anestis A, et al. Estrogen receptor beta increases sensitivity to enzalutamide in androgen receptor-positive triple-negative breast cancer. J. Cancer Res. Clin. Oncol. 2019;145:1221–1233.
    1. Chang J, et al. Expression of ERβ gene in breast carcinoma and the relevance in neoadjuvant therapy. Oncol. Lett. 2017;13:1641–1646.
    1. Yan, M., Rayoo, M., Takano, E. A., kConFab Investigators & Fox, S. B. Nuclear and cytoplasmic expressions of ERβ1 and ERβ2 are predictive of response to therapy and alters prognosis in familial breast cancers. Breast Cancer Res. Treat.126, 395–405 (2011).
    1. Gruvberger-Saal SK, et al. Estrogen receptor beta expression is associated with tamoxifen response in ERalpha-negative breast carcinoma. Clin. Cancer Res. 2007;13:1987–1994.
    1. Hopp TA, et al. Low levels of estrogen receptor beta protein predict resistance to tamoxifen therapy in breast cancer. Clin. Cancer Res. Off. J. Am. Assoc. Cancer Res. 2004;10:7490–7499.
    1. Hamilton-Burke W, et al. Phosphorylation of estrogen receptor beta at serine 105 is associated with good prognosis in breast cancer. Am. J. Pathol. 2010;177:1079–1086.
    1. Guo L, Zhang Y, Zhang W, Yilamu D. Correlation between estrogen receptor β expression and the curative effect of endocrine therapy in breast cancer patients. Exp. Ther. Med. 2014;7:1568–1572.
    1. Guo L, Zhang YU, Yilamu D, Liu S, Guo C. ERβ overexpression results in endocrine therapy resistance and poor prognosis in postmenopausal ERα-positive breast cancer patients. Oncol. Lett. 2016;11:1531–1536.
    1. Austin D, et al. Estrogen receptor-beta is a potential target for triple negative breast cancer treatment. Oncotarget. 2018;9:33912–33930.
    1. Hamilton N, et al. Biologic roles of estrogen receptor-β and insulin-like growth factor-2 in triple-negative breast cancer. BioMed Res. Int. 2015;2015:925703.
    1. Guo L, et al. Expression and prognostic value of estrogen receptor β in patients with triple-negative and triple-positive breast cancer. Exp. Ther. Med. 2015;9:2147–2150.
    1. Wang J, et al. ERβ1 inversely correlates with PTEN/PI3K/AKT pathway and predicts a favorable prognosis in triple-negative breast cancer. Breast Cancer Res. Treat. 2015;152:255–269.
    1. Honma N, et al. Clinical importance of estrogen receptor-beta evaluation in breast cancer patients treated with adjuvant tamoxifen therapy. J. Clin. Oncol. 2008;26:3727–3734.
    1. de Cremoux P, et al. Quantitation of estradiol receptors alpha and beta and progesterone receptors in human breast tumors by real-time reverse transcription-polymerase chain reaction. Correlation with protein assays. Biochem. Pharmacol. 2002;64:507–515.
    1. Stefanou D, Batistatou A, Briasoulis E, Arkoumani E, Agnantis NJ. Estrogen receptor beta (ERbeta) expression in breast carcinomas is not correlated with estrogen receptor alpha (ERalpha) and prognosis: The Greek experience. Eur. J. Gynaecol. Oncol. 2004;25:457–461.
    1. Andersson S, et al. Insufficient antibody validation challenges oestrogen receptor beta research. Nat. Commun. 2017;8:15840.
    1. Nelson AW, et al. Comprehensive assessment of estrogen receptor beta antibodies in cancer cell line models and tissue reveals critical limitations in reagent specificity. Mol. Cell. Endocrinol. 2017;440:138–150.
    1. Saal LH, et al. The Sweden Cancerome Analysis Network - Breast (SCAN-B) Initiative: A large-scale multicenter infrastructure towards implementation of breast cancer genomic analyses in the clinical routine. Genome Med. 2015;7:20.
    1. Rydén L, et al. Minimizing inequality in access to precision medicine in breast cancer by real-time population-based molecular analysis in the SCAN-B initiative. Br. J. Surg. 2018;105:e158–e168.
    1. Brueffer C, et al. Clinical value of RNA sequencing-based classifiers for prediction of the five conventional breast cancer biomarkers: A report from the population-based multicenter Sweden cancerome analysis network-breast initiative. JCO Precis. Oncol. 2018;2:1–18.
    1. Brueffer C, et al. The mutational landscape of the SCAN-B real-world primary breast cancer transcriptome. EMBO Mol. Med. 2020;12:12118.
    1. Dahlgren M, et al. Preexisting somatic mutations of estrogen receptor alpha (esr1) in early-stage primary breast cancer. JNCI Cancer Spectr. 2021;5:028.
    1. Bièche I, et al. Quantification of estrogen receptor α and β expression in sporadic breast cancer. Oncogene. 2001;20:8109–8115.
    1. Collins L, Botero M, Schnitt S. Bimodal frequency distribution of estrogen receptor immunohistochemical staining results in breast cancer. Am. J. Clin. Pathol. 2005;123:16–20.
    1. Thomas C, et al. Estrogen receptor β-mediated inhibition of actin-based cell migration suppresses metastasis of inflammatory breast cancer. Cancer Res. 2021;81:2399–2414.
    1. Koboldt DC, et al. Comprehensive molecular portraits of human breast tumours. Nature. 2012;490:61–70.
    1. Xie Y, et al. Breast cancer gene expression datasets do not reflect the disease at the population level. NPJ Breast Cancer. 2020;6:39.
    1. Benzina S, et al. Pax-5 is a potent regulator of E-cadherin and breast cancer malignant processes. Oncotarget. 2017;8:12052–12066.
    1. Leblanc N, Harquail J, Crapoulet N, Ouellette RJ, Robichaud GA. Pax-5 inhibits breast cancer proliferation through miR-215 up-regulation. Anticancer Res. 2018;38:5013–5026.
    1. Garaud S, et al. Tumor-infiltrating B cells signal functional humoral immune responses in breast cancer. JCI Insight. 2019;4:129641.
    1. Shen M, Wang J, Ren X. New insights into tumor-infiltrating B lymphocytes in breast cancer: Clinical impacts and regulatory mechanisms. Front. Immunol. 2018;9:470.
    1. Larsson O, Tian B, Sonenberg N. Toward a genome-wide landscape of translational control. Cold Spring Harb. Perspect. Biol. 2013;5:a012302.
    1. Yakimchuk K, Jondal M, Okret S. Estrogen receptor α and β in the normal immune system and in lymphoid malignancies. Mol. Cell. Endocrinol. 2013;375:121–129.
    1. Kovats S. Estrogen receptors regulate innate immune cells and signaling pathways. Cell. Immunol. 2015;294:63–69.
    1. Khan D, Ansar Ahmed S. The immune system is a natural target for estrogen action: Opposing effects of estrogen in two prototypical autoimmune diseases. Front. Immunol. 2016;6:635.
    1. National Cancer Institute. Surveillance, Epidemiology, and End Results (SEER) program website. . Published 2018. Accessed August 10, 2021.
    1. Nielsen TO, et al. Immunohistochemical and clinical characterization of the basal-like subtype of invasive breast carcinoma. Clin. Cancer Res. 2004;10:5367–5374.
    1. Dent R, et al. Triple-negative breast cancer: Clinical features and patterns of recurrence. Clin. Cancer Res. 2007;13:4429–4434.
    1. Bianchini G, Balko JM, Mayer IA, Sanders ME, Gianni L. Triple-negative breast cancer: Challenges and opportunities of a heterogeneous disease. Nat. Rev. Clin. Oncol. 2016;13:674–690.
    1. Wolff AC, et al. Recommendations for human epidermal growth factor receptor 2 testing in breast cancer: American Society of Clinical Oncology/College of American Pathologists clinical practice guideline update. J. Clin. Oncol. 2013;31:3997–4013.
    1. Saal LH, et al. BioArray Software Environment (BASE): A platform for comprehensive management and analysis of microarray data. Genome Biol. 2002;3:SOFTWARE0003.
    1. Vallon-Christersson J, Nordborg N, Svensson M, Häkkinen J. BASE - 2nd generation software for microarray data management and analysis. BMC Bioinformatics. 2009;10:330.
    1. Häkkinen, J., Nordborg, N., Månsson, O. & Vallon-Christersson, J. Implementation of an open source software solution for laboratory information management and automated rnaseq data analysis in a large-scale cancer genomics initiative using base with extension package reggie. bioRxiv 038976 (2016) 10.1101/038976.
    1. Picard toolkit. Broad Institute, GitHub repository (2019).
    1. Bolger AM, Lohse M, Usadel B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinforma. Oxf. Engl. 2014;30:2114–2120.
    1. Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat. Methods. 2012;9:357–359.
    1. Kim D, Paggi JM, Park C, Bennett C, Salzberg SL. Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat. Biotechnol. 2019;37:907–915.
    1. Pertea M, et al. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat. Biotechnol. 2015;33:290–295.
    1. Soneson C, Love M, Robinson M. Differential analyses for RNA-seq: Transcript-level estimates improve gene-level inferences. F1000Research. 2016;4:1521.
    1. Colaprico A, et al. TCGAbiolinks: An R/Bioconductor package for integrative analysis of TCGA data. Nucleic Acids Res. 2016;44:e71–e71.
    1. Ritchie ME, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43:e47.
    1. Sergushichev, A. A. An algorithm for fast preranked gene set enrichment analysis using cumulative statistic calculation. bioRxiv 060012 (2016). 10.1101/060012.
    1. Korotkevich, G. et al. Fast gene set enrichment analysis. bioRxiv 060012 (2021). 10.1101/060012.
    1. Wang J, Duncan D, Shi Z, Zhang B. WEB-based GEne SeT AnaLysis Toolkit (WebGestalt): Update 2013. Nucleic Acids Res. 2013;41:W77–W83.
    1. Zhang B, Kirov S, Snoddy J. WebGestalt: An integrated system for exploring gene sets in various biological contexts. Nucleic Acids Res. 2005;33:W741–748.
    1. Wang J, Vasaikar S, Shi Z, Greer M, Zhang B. WebGestalt 2017: A more comprehensive, powerful, flexible and interactive gene set enrichment analysis toolkit. Nucleic Acids Res. 2017;45:W130–W137.
    1. Liu J, et al. An integrated tcga pan-cancer clinical data resource to drive high-quality survival outcome analytics. Cell. 2018;173:400–416.e11.

Source: PubMed

3
Abonnieren