Meningiomas are the most common primary intracranial tumor with current classification offering limited therapeutic guidance. Here, we interrogated meningioma enhancer landscapes from 33 tumors to stratify patients based upon prognosis and identify novel meningioma-specific dependencies. Enhancers robustly stratified meningiomas into three biologically distinct groups (adipogenesis/cholesterol, mesodermal, and neural crest) distinguished by distinct hormonal lineage transcriptional regulators. Meningioma landscapes clustered with intrinsic brain tumors and hormonally responsive systemic cancers with meningioma subgroups, reflecting progesterone or androgen hormonal signaling. Enhancer classification identified a subset of tumors with poor prognosis, irrespective of histologic grading. Superenhancer signatures predicted drug dependencies with superior in vitro efficacy to treatment based upon the NF2 genomic profile. Inhibition of DUSP1, a novel and druggable meningioma target, impaired tumor growth in vivo. Collectively, epigenetic landscapes empower meningioma classification and identification of novel therapies. SIGNIFICANCE: Enhancer landscapes inform prognostic classification of aggressive meningiomas, identifying tumors at high risk of recurrence, and reveal previously unknown therapeutic targets. Druggable dependencies discovered through epigenetic profiling potentially guide treatment of intractable meningiomas.This article is highlighted in the In This Issue feature, p. 1611.

Brain tumors are dynamic complex ecosystems with multiple cell types. To model the brain tumor microenvironment in a reproducible and scalable system, we developed a rapid three-dimensional (3D) bioprinting method to construct clinically relevant biomimetic tissue models. In recurrent glioblastoma, macrophages/microglia prominently contribute to the tumor mass. To parse the function of macrophages in 3D, we compared the growth of glioblastoma stem cells (GSCs) alone or with astrocytes and neural precursor cells in a hyaluronic acid-rich hydrogel, with or without macrophage. Bioprinted constructs integrating macrophage recapitulate patient-derived transcriptional profiles predictive of patient survival, maintenance of stemness, invasion, and drug resistance. Whole-genome CRISPR screening with bioprinted complex systems identified unique molecular dependencies in GSCs, relative to sphere culture. Multicellular bioprinted models serve as a scalable and physiologic platform to interrogate drug sensitivity, cellular crosstalk, invasion, context-specific functional dependencies, as well as immunologic interactions in a species-matched neural environment.

Use of menopausal hormone therapy (MHT) is associated with increased risk for breast cancer. However, the relevant mechanisms and its interaction with genetic variants are not fully understood. We conducted a genome-wide interaction analysis between MHT use and genetic variants for breast cancer risk in 27,585 cases and 34,785 controls from 26 observational studies. All women were post-menopausal and of European ancestry. Multivariable logistic regression models were used to test for multiplicative interactions between genetic variants and current MHT use. We considered interaction p-values < 5 × 10^-8 as genome-wide significant, and p-values < 1 × 10^-5 as suggestive. Linkage disequilibrium (LD)-based clumping was performed to identify independent candidate variants. None of the 9.7 million genetic variants tested for interactions with MHT use reached genome-wide significance. Only 213 variants, representing 18 independent loci, had p-values < 1 × 10⁵. The strongest evidence was found for rs4674019 (p-value = 2.27 × 10^-7), which showed genome-wide significant interaction (p-value = 3.8 × 10^-8) with current MHT use when analysis was restricted to population-based studies only. Limiting the analyses to combined estrogen-progesterone MHT use only or to estrogen receptor (ER) positive cases did not identify any genome-wide significant evidence of interactions. In this large genome-wide SNP-MHT interaction study of breast cancer, we found no strong support for common genetic variants modifying the effect of MHT on breast cancer risk. These results suggest that common genetic variation has limited impact on the observed MHT-breast cancer risk association.

Abstract Background: Use of menopausal hormone therapy (MHT) is associated with increased risk for breast cancer. However, the relevant mechanisms and its interaction with genetic variants are not fully understood. Methods: We conducted a genome-wide interaction analysis between MHT use and genetic variants for breast cancer risk in 27,585 cases and 34,785 controls from 26 observational studies. All women were post-menopausal and of European ancestry. Multivariable logistic regression models were used to test for multiplicative interactions between genetic variants and current MHT use. We considered interaction p-values<5x10-8 as genome-wide significant, and p-values<1x10-5 as suggestive. Linkage disequilibrium (LD)-based clumping was performed to identify independent candidate variants.Results: None of the 9.7 million genetic variants tested for interactions with MHT use reached genome-wide significance. Only 213 variants, representing 18 independent loci, had p-values<1x105. The strongest evidence was found for rs4674019 (p-value=2.27x10-7), which showed genome-wide significant interaction (p-value=3.8x10-8) with current MHT use when analysis was restricted to population-based studies only. Limiting the analyses to combined estrogen-progesterone MHT use only or to estrogen receptor (ER) positive cases did not identify any genome-wide significant evidence of interactions. Conclusions: In this large genome-wide SNP-MHT interaction study of breast cancer, we found no strong support for common genetic variants modifying the effect of MHT on breast cancer risk. These results suggest that common genetic variation has limited impact on the observed MHT–breast cancer risk association.

Background

Genome-wide association studies (GWAS) have identified multiple common breast cancer susceptibility variants. Many of these variants have differential associations by estrogen receptor (ER) status, but how these variants relate with other tumor features and intrinsic molecular subtypes is unclear.

Methods

Among 106,571 invasive breast cancer cases and 95,762 controls of European ancestry with data on 173 breast cancer variants identified in previous GWAS, we used novel two-stage polytomous logistic regression models to evaluate variants in relation to multiple tumor features (ER, progesterone receptor (PR), human epidermal growth factor receptor 2 (HER2) and grade) adjusting for each other, and to intrinsic-like subtypes.

Results

Eighty-five of 173 variants were associated with at least one tumor feature (false discovery rate < 5%), most commonly ER and grade, followed by PR and HER2. Models for intrinsic-like subtypes found nearly all of these variants (83 of 85) associated at p < 0.05 with risk for at least one luminal-like subtype, and approximately half (41 of 85) of the variants were associated with risk of at least one non-luminal subtype, including 32 variants associated with triple-negative (TN) disease. Ten variants were associated with risk of all subtypes in different magnitude. Five variants were associated with risk of luminal A-like and TN subtypes in opposite directions.

Conclusion

This report demonstrates a high level of complexity in the etiology heterogeneity of breast cancer susceptibility variants and can inform investigations of subtype-specific risk prediction.