Computational repositioning and preclinical validation of mifepristone for human vestibular schwannoma

The computational repositioning of existing drugs represents an appealing avenue for identifying effective compounds to treat diseases with no FDA-approved pharmacotherapies. Here we present the largest meta-analysis to date of differential gene expression in human vestibular schwannoma (VS), a debilitating intracranial tumor, and use these data to inform the first application of algorithm-based drug repositioning for this tumor class. We apply an open-source computational drug repositioning platform to gene expression data from 80 patient tumors and identify eight promising FDA-approved drugs with potential for repurposing in VS. Of these eight, mifepristone, a progesterone and glucocorticoid receptor antagonist, consistently and adversely affects the morphology, metabolic activity, and proliferation of primary human VS cells and HEI-193 human schwannoma cells. Mifepristone treatment reduces VS cell viability more significantly than cells derived from patient meningiomas, while healthy human Schwann cells remain unaffected. Our data recommend a Phase II clinical trial of mifepristone in VS.

Vestibular schwannoma (VS) is the fourth most common intracranial tumor and the most common tumor of the cerebellopontine angle, arising from neoplastic Schwann cells of the vestibular nerve. No drug is FDA-approved to treat VS. In 95% of patients, these tumors cause debilitating sensorineural hearing loss (SNHL) and tinnitus, and can also lead to dizziness and facial paralysis. Bilateral VSs are the hallmark of neurofibromatosis type 2 (NF2), an autosomal dominant disorder caused by inactivation or loss of both alleles of the NF2 gene. If left untreated, growing VSs can compress the brainstem and lead to death. Mutations in the NF2 gene are identified in 100% of NF2-associated VSs and 66% of sporadically arising VSs 1 . Though mechanisms of VS-induced SNHL are multifactorial 2,3 , NF2-associated SNHL often correlates with VS size 2,4 . This observation suggests that slowing or inhibiting VS growth may not only prolong a patient's time to surgical intervention, but also minimize or prevent associated SNHL, substantially improving quality of life. Current treatment options for VS are limited to surgical resection and radiation therapy, both of which carry massive risks for patients, including facial nerve paralysis and loss of hearing. Identification of a drug with the potential to slow or halt VS growth, thereby preserving acoustic hearing and delaying life-threatening complications, represents a major unmet need.
The computational repositioning of FDA-approved drugs, in which data-driven analyses of gene-compound interactions catalyze the pursuit of new indications for approved drugs, provides a transformative avenue for therapeutic innovation 5 . Using novel computational strategies to capitalize on emerging high-throughput genomic SCiEntifiC REPORtS | (2018) 8:5437 | DOI: 10.1038/s41598-018-23609-7 information, it is possible to identify safe drugs with evidence-based potential for repurposing in other conditions, which can eliminate the need for a Phase I safety trial and expedite Phase II efficacy trials. Computational drug repositioning represents an appealing approach for diseases with no approved pharmacotherapies, narrowing down a potentially infinite biochemical playing field to a few high-potential candidates.
One approach for surveying drug-based perturbations relies on comparing the RNA-level gene expression profile specific to a given disease to large, pre-generated databases of multi-drug exposure profiles or known gene-drug interactions. The most commonly used tool to interrogate such data is the Broad Connectivity Map (CMAP) 6 , an online platform for matching differentially expressed gene sets from a disease of interest to a library of drug exposure profiles generated from human cell lines. CMAP relies on a modified Kolmogorov-Smirnov enrichment test to rank potentially effective compounds. The natural appeal of its resulting "connectivity score" has led to a large number of CMAP-relevant studies and expansion of these methodologies in targeted disease areas [7][8][9][10] . Such studies have led to several high-profile drug repositioning recommendations, such as the well-known case of the anticonvulsant topimarate for use in inflammatory bowel disease (IBD) 11 . Though gene expression data initially indicated that topimarate may be a sensible choice for repositioning in IBD, medical experts pushed back on this result because one of the most frequent severe adverse effects of topimarate is diarrhea, from which IBD patients already disproportionately suffer 12 . Therefore, a holistic evaluation of drug safety within the potential confounds of a new disease indication should never be neglected.
In order to identify an FDA-approved drug with potential for repositioning in VS, we conducted the largest meta-analysis of human VS transcriptomes to date and applied a computational drug repositioning algorithm to match differential gene expression patterns characteristic of VS with known interactions between genes and FDA-approved drugs. Then, after generating a shortlist of drugs with high potential for repositioning, we relied on the clinical expertise of neuro-otologists specializing in VS management to select mifepristone, a progesterone receptor antagonist approved for use in medical abortion, as a candidate worthy of further validation. We use primary human VS cells and immortalized human cell lines to validate mifepristone for repositioning in the treatment of this debilitating tumor. Mifepristone is known to cross the blood-brain barrier and carries a mild side effect profile even in long-term clinical trials. We recommend a Phase II clinical trial of mifepristone in VS.

Results
Computational drug repositioning based on differential gene expression reveals drugs with high potential for repurposing in VS. To identify FDA-approved drugs with potential for repositioning in VS, we conducted a computational screen using the open-source drug repositioning platform ksRepo, developed to screen expression profiles from any microarray or sequencing platform against any available database of gene-drug interactions 13 . ksRepo uses a modified Kolmogorov-Smirnov statistic to compare a ranked list of differentially expressed genes (DEGs) characteristic of a given disease with transcriptional signatures associated with drugs known to interact with those genes, as publicly stored in the Comparative Toxicogenomics Database (CTD) 14 . From that list of drugs, ksRepo selects for compounds with entries in DrugBank, a compendium of FDA-approved drugs 15 . The output is a list of FDA-approved drugs hypothesized to modulate genes with aberrant expression patterns in patients with disease (Fig. 1a). We selected ksRepo for this specific analysis because it does not require input data to be generated via a specific platform (as CMAP does); allows us to choose the database we use (the CTD is manually curated, allowing for high-fidelity associations); affords us flexibility in selecting which drugs to test (rather than relying on pre-compiled, pre-chosen gene signatures); and is compatible with gene-level meta-analysis. ksRepo was recently shown to be successful against a meta-analysis of DEGs from five independent prostate cancer datasets, from which this algorithm successfully predicted significance for five approved therapies in prostate cancer treatment 13 .
To provide robust input to ksRepo, we conducted the largest meta-analysis of primary human VS tissue to date, comprising genome-wide transcriptional microarray data from 80 tumors and 16 control nerves (Fig. 1b). Combined analysis of expression data from two large published datasets, one publicly available (NCBI GEO, GSE39645) 16 and one published and donated upon request 17 , yields 1,335 genes found to be commonly and concordantly dysregulated in VS, with 405 reaching significance after Bonferroni correction for multiple hypothesis testing (p < 0.05) (Supplementary Data Files S1, S2). ksRepo takes the entire meta-analytic expression profile as input in order to screen a comprehensive picture of tumor-related gene expression against the known interactions of 1,155 FDA-approved drugs. As 16 of 80 VSs in the meta-analysis were harvested from patients with NF2, we also conducted a separate, parallel analysis comprising only NF2-associated tumors. ksRepo returned 36 drugs with potential for repositioning from the complete VS meta-analysis and 68 drugs from the analysis specific to NF2-associated tumors (Supplementary Data Files S3, S4). Though there is no consensus frontline pharmacologic therapy for VS, ksRepo successfully identified multiple classes of drugs that have shown limited efficacy against this tumor, such as non-steroidal anti-inflammatory drugs (NSAIDs), glucocorticoids, tyrosine kinase inhibitors, and histone deacetylase inhibitors (Supplementary Data Files S2, S4).
Eight drugs appeared in both analyses and were prioritized for preclinical validation (Fig. 1c). Out of all FDA-approved drugs, this group of eight demonstrates significant enrichment for anti-inflammatory drugs, hormone-related compounds, and anti-neoplastic agents (Fig. 1d). Specifically, though anti-inflammatory drugs represent just 1.36% of all FDA-approved drugs, four of our eight drugs were classified as anti-inflammatory (50%). Similarly, hormone-related compounds demonstrated a 13.9-fold enrichment (37.5% of recommended drugs from 2.70% of all hormone-related drugs,), and anti-neoplastic agents a 10.2-fold increase (37.5% of recommended drugs from 3.68% of all anti-neoplastic drugs). After thorough theoretical and medical investigation of each of these eight drugs in relation to 1) its medical potential to slow or halt VS growth and 2) its potential to translate safely and effectively to VS patients, mifepristone emerged as the most promising lead candidate. Mifepristone reduces metabolic activity and proliferation of HEI-193 human schwannoma cells in culture. Mifepristone (RU486) is a progesterone and glucocorticoid receptor antagonist approved by the FDA for use in medical abortion. This steroid analog is able to cross the blood-brain barrier 18 and has been shown to provide palliative benefits to patients with other intracranial tumors, such as glioblastoma multiforme 18 and meningioma 19 . In vitro, mifepristone produces antiproliferative effects on cervical 20 , breast 21 , endometrial 22 , ovarian 23 , and prostate cancer cells 24 , regardless of progesterone receptor expression 25 . In human trials, mifepristone administration has been documented to improve quality of life for patients suffering from advanced thymic, renal, colon, leukemic, and pancreatic cancers 26,27 . Long-term administration of oral mifepristone is well tolerated by adults and carries only a mild toxicity profile 28 . Independently of our computational repositioning analysis, when gene expression data from our meta-analysis was input to Ingenuity Pathway Analysis (QIAGEN Inc., https://www.qiagenbioinformatics.com/products/ingenuitypathway-analysis), mifepristone was predicted as a significant upstream regulator of the resulting networks (p = 4.26*10 −5 ) and theorized to act upstream of inflammatory markers characteristic of VS, such as TNF and NFkB 29,30 (Supplementary Fig. S1). We did not perform further experiments to interrogate this hypothesis.
Administration of mifepristone to HEI-193 immortalized human schwannoma cells in culture produces a significant, dose-dependent response in metabolic activity (Fig. 2a); a significant reduction in cell confluence (Fig. 2b,c); and a dramatic decline in cellular proliferation ( Fig. 2d-f) without producing a significant increase in apoptosis ( Fig. 2g; Supplementary Fig. S2). Under mifepristone treatment, HEI-193 cells assume a long, thin, spindle-like shape ( Fig. 2c; Supplementary Videos S1, S2). This observation is reported in previous studies of mifepristone treatment of ovarian, breast, prostate, and nerve cells, where such morphological changes are attributed to dysregulated distribution of f-actin and tubulin proteins in the cytoskeleton 31 . Continuous, dynamic actin remodeling is characteristic of NF2-deficient schwannoma cells 32 , as the NF2 protein product, merlin, is known to selectively bind f-actin 33 . Phalloidin staining confirms the shrunken, crumpled appearance of f-actin in mifepristone-treated cells (Fig. 2h,i). Enzyme-linked immunosorbent assay for progesterone in conditioned cell culture medium reveals higher levels of progesterone in medium collected from mifepristone-treated cells than from untreated cells, suggesting that mifepristone may compete with progesterone for receptor binding ( Supplementary Fig. S3).

Mifepristone reduces metabolic activity and proliferation of primary human VS cells in cul-
ture, regardless of NF2 mutation. Current, large-scale meta-analyses of drug toxicology, bioavailability, and efficacy in animal models reveal a shocking lack of predictive power when compared to human data 34,35 . Accordingly, the U.S. National Research Council has recommended the substitution of model animal testing with in vitro human cell-based assays and in silico modeling of diseases and networks 36 . We evaluated the effect of mifepristone applied directly to primary human VS cells. Fresh VS tissue samples from ten human patients undergoing tumor resection surgery were collected and schwannoma cells grown in the laboratory according to published protocols 37,38 . When we received research tissue in excess of the amount necessary for cell culture, we also performed single-gene sequencing of the NF2 gene (Fig. 3a). Four of six sequenced VSs exhibited novel mutations in the NF2 gene, a fraction consistent with published literature 1 ( Fig. 3b; Supplementary Fig. S4).
When applied to primary human VS cells, mifepristone produced a dose-dependent response in metabolic activity and a dramatic reduction in cellular proliferation (Fig. 4a,b). Live-cell fluorescence imaging revealed a marked increase in cytotoxicity in primary cultures ( Fig. 4c-e). No correlation in drug response with tumor size or NF2 mutation type was observed.

Mifepristone produces a more significant effect on VS cells than arachnoid cells, but does not adversely affect healthy Schwann cells.
In a long-term clinical trial of mifepristone for unresectable meningioma, minor responses resulting in clinical benefit were noted in eight of 28 patients 28 , though a subsequent Phase III trial reported no difference between treatment and placebo 39 . To evaluate the effect of this drug on schwannoma cells in comparison to meningioma cells, we compared mifepristone-treated primary human VS and HEI-193 cells to immortalized human arachnoid cells in which the NF2 gene has been excised by CRISPR 40 . Primary VS cells responded more dramatically to mifepristone than arachnoid cells with or without the NF2 gene ( Fig. 4f), suggesting that schwannoma cells are more responsive than meningioma cells to this drug.
Additionally, to ensure that mifepristone administration did not lead to adverse effects among healthy human Schwann cells, primary human Schwann cells were cultured from eight great auricular nerves (GANs) harvested from patients undergoing benign parotidectomy or neck dissection 37,38 . Treatment of these cells with mifepristone did not cause appreciable changes in metabolic activity (Fig. 4g). Preliminary testing of clinically reasonable concentrations of other drugs recommended by ksRepo, including adenosine monophosphate, gold sodium thiomalate, succimer, and methylprednisolone, showed no effect on the metabolic activity of HEI-193 cells, though prednisolone produced modest effects when applied to primary human VS cells (Fig. 5a-h).

Discussion
The in silico repositioning of mifepristone for human VS using pooled human transcriptomic data, as well as the preclinical validation of this drug on primary human cells, constitutes a powerful case for the computational identification of novel indications for FDA-approved drugs. Future work will include more detailed pathway analyses of the shared transcriptional responses elicited by all top compounds selected by ksRepo in order to unveil new druggable targets and elucidate additional mechanisms that may drive VS growth. Additionally, as computational drug screening matures, integrating gene expression data with gene mutation data, comprehensive proteomics, pathway analyses, and human exposure data will allow us to refine the quality of algorithm-recommended results, decreasing the odds of false positives and false negatives.
Despite the fact that there is no approved pharmacotherapy for VS, ksRepo successfully identified a variety of drug classes that have shown variable degrees of promise when used against this tumor. Specifically, ksRepo's recommendation of the NSAID drugs celecoxib, nimeluside, piroxicam, and indomethacin highlights the finding that NSAIDs are cytostatic against VS in vitro 30 , and in vivo, retrospective studies show that aspirin intake correlates with decreased tumor volume 41,42 . Accordingly, a prospective, randomized clinical trial for aspirin (a drug that did not survive Bonferroni correction for multiple hypothesis testing after ksRepo identification) in VS is set to begin this year. Similarly, the recommendation of prednisolone and methylprednisolone supports the fact that glucocorticoids are the only treatment consistently used with any measure of success against VS-associated sudden SNHL, by mechanisms that remain unknown 43 .
Given the success of tyrosine kinase and histone deacetylase inhibitors in the treatment of other solid neoplasms 44 , it is not surprising that ksRepo spotlights four of these drugs for potential efficacy in VS. Gefitinib is an EGFR/tyrosine kinase inhibitor that showed early promise in the treatment of VS in vitro 45 ; a retrospective human study of its analog, erlotinib, in NF2 patients did not result in a measurable reduction in tumor volume overall, though three patients reported hearing improvement 46 . Additionally, ksRepo's recommendation of histone deacetylase inhibitors belinostat, vorinostat, and panobinostat reflects the preclinical success of AR-42, a novel histone deacetylase inhibitor that is the subject of an active clinical trial 47,48 .
Mifepristone belongs to a surprising class of drugs consistently recommended by ksRepo, comprising compounds specific to receptors for the sex hormones estrogen and progesterone, including norgestimate, levonorgestrel, raloxifene, and tamoxifen (Supplementary Data File S2). The existing VS literature contains conflicting reports regarding estrogen and progesterone receptor expression in these tumors, even when measured using nearly identical methods [49][50][51][52] . It is unclear what role sex hormone-related signaling might play in VS development or progression. However, a robust body of clinical literature dating back to Harvey Cushing documents the fact that VSs often present during pregnancy or accelerate growth during this period, indicating a potential correlation between sex hormone receptor-mediated signaling and rate of VS growth 53,54 .
Our approach is limited by the fact that ksRepo does not take into account bevacizumab, a monoclonal antibody targeting vascular endothelial growth factor A (VEGF-A). As a biologic drug designed to target a single protein, bevacizumab may not conform well to a drug repositioning approach focusing on genome-wide expression differences, and our 80-tumor meta-analysis reveals that VEGF-A is not significantly upregulated in VSs when compared to GAN controls (Bonferroni-corrected meta-analytic p = 1; Supplementary Data File S1). Regardless, though no human trial of any drug has shown outstanding results against VS, the most promising effects have been associated with bevacizumab. In six retrospective studies conducted among NF2 patients, bevacizumab decreased tumor volume by at least 20% and was associated with a measurable hearing response in most patients 55 . A recent open-label Phase II trial reported radiological responses in six of 14 NF2-associated tumors, with an accompanying hearing response in five 56 . However, the effect of this antibody is often transient, and when generalized across all studies, seems to produce a response in only approximately 50% of patients 55 . Interestingly, bevacizumab does not seem to be clinically effective against NF2-associated meningiomas, in which VEGF pathway upregulation may not be the driving force behind angiogenesis 57 , supporting the hypothesis that this drug may act on tumor vasculature but not VS cells. Uncomfortable administration and serious adverse events observed among long-term users, such as hypertension, renal failure, and hemorrhage, make this drug untenable for multi-year administration, and patients generally resort to surgery 58 .
A drug like mifepristone, which successfully crosses the blood-brain barrier 18 , drastically inhibits human VS cell proliferation, and is documented to have a mild adverse effect profile even when taken daily for up to 13 years 28 could potentially be administered throughout a patient's lifetime, prolonging time to surgery or eliminating the need for resection. Few prospective clinical trials are conducted in the primarily surgical management of this debilitating tumor. However, mifepristone is safe, approved for human use, associated with no serious adverse effects known to be compounded by VS, and deserving of further attention in a devastating disease with no FDA-approved therapy. We recommend a Phase II clinical trial of mifepristone in VS.

GEO Dataset
Processing. The NCBI GEO dataset used in this study is GSE39645, an Affymetrix Human Gene 1.0 ST chip-based gene expression study of VS which contained data from 28 patients with sporadic VS, 3 patients with NF2-associated VS, and 8 control nerve samples 16 . Data for GSE39645 was accessed through the NCBI GEO portal and analyzed using the integrated GEO2R tool 59 . As input for GEO2R, we classified each sample within a GEO series as either normal tissue or VS tissue. The GEO2R analysis was performed on both the full dataset (sporadic and NF2 combined), and a subset of samples containing only NF2-syndromic schwannomas. GEO2R provides a list of probes and corresponding gene symbols ranked according to their degree of differential expression (as calculated using the limma package in R 60 ), and includes p-values and t-statistics for differential expression. Following GEO2R analysis, all results were imported into R 61 and probe-level differential expression was consolidated to gene-level differential expression using a custom pipeline: t-statistic values were converted to Cohen's d statistic values and standard error values 62 . Resulting values were combined by gene using a fixed effects meta-analysis (as implemented in the meta.summaries function from the rmeta package in R 63 . Probes without gene annotations were removed from gene-level consolidation. Following consolidation, significantly differentially expressed genes were taken to be those with a Bonferroni-corrected significance of less than 0.05. Additional VS Dataset Processing. Raw Affymetrix Human Genome U219 gene expression data (.CEL files) for 36 patients with sporadic VS, 16 patients with NF2 syndrome-associated VS, and 7 control nerves were generously donated by Agnihotri et al. 17 . CEL files were loaded into R using the justRMA function from the affy package in R 64 . justRMA is an automated tool that both performs normalization using the Robust Multi-Array Average method 65 and also automatically annotates all probes in the normalized dataset using the Org.Hs.eg.db annotation database package 66 . Normalization was performed on the full dataset and the NF2-associated schwannomas, as above. Mirroring the GEO2R analysis, each normalized dataset was analyzed using limma and consolidated to gene-level differential expression using the custom pipeline described above. As above, significantly differentially expressed genes were taken to be those with a Bonferroni-corrected significance of less than 0.05.

Meta-Analysis of 80 VS Samples and ksRepo Prediction.
To robustly determine differential expression between VS and normal tissues, gene-level data from GSE39645 16 and Agnihotri et al. 17 were meta-analyzed by first removing genes that were not measured in both the Affymetrix Human Gene 1.0 ST chip and the Affymetrix Human Genome U219 chip, and subsequently combining Cohen's d and standard error values using a fixed-effects meta-analysis (again using meta.summaries). Meta-analysis was performed for the full GSE39645 and Agnihotri datasets, as well as for NF2-associated tumors exclusively. Following meta-analysis, the remaining genes were ranked according to their meta-analytic p-values to generate a gene list for further analysis using ksRepo (package available for download at https://github.com/adam-sam-brown/ksRepo, and described in Brown et al. 13 . ksRepo is a gene-based drug repositioning method that uses a modified Kolmogorov-Smirnov (KS) statistic to identify promising drug repositioning opportunities. ksRepo requires a database of compound-gene interactions, which are compared with the ranked meta-analytic gene lists from above. For this analysis, the ksRepo built-in Comparative Toxicogenomics Database (CTD) dataset was selected. The CTD provides a curated resource that links small chemical entities to genes (e.g., gene or protein expression influences) from the scientific literature on numerous model organisms and humans 14 . ksRepo contains a subset of the CTD, containing human-derived interactions between 1,268 unique drugs and 18,041 unique human genes. Drugs in the CTD subset were chosen based on case-insensitive matches between CTD names and names/synonyms for FDA-approved drugs downloaded from DrugBank 15 . The ksRepo output provides both the resampled p-value and FDR value. For the full dataset ksRepo analysis and the NF2-only ksRepo analysis, significant compounds were those for which the FDR was less than 0.05.
Human Specimen Collection and Primary Cell Culture. Surgical VS and GAN specimens were collected and processed according to protocols approved by the Human Studies Committee of Massachusetts General Hospital and Massachusetts Eye and Ear (Board Reference #14-148 H). Written informed consent was obtained from all subjects prior to inclusion in this study and all procedures were conducted in accordance with the Helsinki Declaration of 1975. Detailed methods for human surgical specimen collection, processing, and culture are previously published 37,38 . VS specimens were harvested from patients undergoing surgical tumor resection, and GAN specimens from healthy patients undergoing benign parotidectomy or neck dissection surgery, during which the GAN is routinely sacrificed. Patients who had received radiation therapy prior to surgery were excluded. After surgical resection, VS or GAN tissue was immediately placed in saline solution, transported to the laboratory, and cultured according to published protocols (Supplementary Methods). All downstream procedures were performed on primary cell cultures or collected culture medium up to but not exceeding two weeks of age in culture to ensure maximal Schwann or schwannoma cell purity 37 .
HEI-193 and Arachnoid Cell Culture. HEI-193 cells are derived from an NF2 patient with spontaneous bilateral vestibular schwannomas and a history of meningioma; these cells express a splice variant of the merlin protein (encoded by the NF2 gene), but neither typical isoform 67 . HEI-193 cells were cultured in DMEM/ F12-containing medium with 10% fetal bovine serum and 1% penicillin and streptomycin mix. Immortalized NF2-null and NF2-expressing arachnoid AC-CRISPR cell lines derived from primary human autopsy specimens were obtained via generous gift from Dr. Vijaya Ramesh at Massachusetts General Hospital 40 . NF2-null and NF2-expressing arachnoid cells were cultured in DMEM with 15% fetal bovine serum and 1% penicillin and streptomycin mix. All cell lines were maintained in an incubator at 37 degrees Celsius with 5% carbon dioxide and treated with drugs 24-36 hours after seeding at between 15,000-25,000 cells per well in 24-well plates. Phase contrast photos of healthy and drug-treated cultures were taken at 10X magnification on an IncuCyte S3 instrument (Essen Bio).
Drug Preparation and Treatment. Primary VS and GAN cultures were treated with mifepristone (Sigma Aldrich, lot #WXBC0031V). Fifteen, 25, 35, and 70 μM mifepristone were prepared by suspending the appropriate amount of drug (in powder form) in dimethyl sulfoxide (DMSO). The resulting drug suspension was diluted in culture medium to the concentration of interest, and drug-containing medium was applied to primary VS, GAN, and HEI-193 cells such that the amount of DMSO applied to cells in culture did not exceed 0.1% (24-well plate, 1 mL medium per well). Cultures were incubated with drug-containing medium or 0.1% DMSO vehicle for 72 hours and then processed for downstream applications. Proliferation Assay. 5′bromo-2′-deoxyuridine (BrdU) was added to label proliferating cells in culture 2 hours before fixation in 4% formalin (paraformaldehyde). Cell membranes were permeabilized with 10 minutes of incubation in 1% Triton X-100 and nuclear membranes with 20 minutes in 2N hydrochloric acid (HCl). Cells were blocked in 5% normal horse serum (NHS) and 1% Triton X-100 and incubated with a primary antibody against BrdU (#OBT0030G, AbD Serotec) overnight, followed by incubation with fluorescent anti-rat immunoglobulin G (AlexaFluor, Life Technologies). Cells were stained with Hoechst 33342 (Invitrogen) and phalloidin/f-actin (ThermoFisher Scientific) and coverslips mounted on slides with VectaShield (Vector Laboratories). The ratio of BrdU-positive to Hoechst-positive nuclei was determined by sampling three random fields of view using a Leica epifluorescence microscope. Manual counts were performed by J.E.S., who was blinded to treatment conditions. Metabolic Activity Assay. The metabolic activity of primary VS and HEI-193 cells was assessed using the colorimetric 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) assay (Life Technologies) according to manufacturer's instructions. Optical density (OD) of each well at 570 nm was read using a spectrophotometer. The average OD value of cells exposed to vehicle (0.1% DMSO) was set to 100% and used to normalize OD values of cells treated with drugs; metabolic activity was then reported as percent change. All statistical testing was performed on raw OD values (see Statistical Analysis).
Flow Cytometry. Apoptotic cell death was assessed using an Annexin V/propidium iodide (PI) staining kit (Miltenyi Biotec). Briefly, HEI-193 cells were seeded into T25 flasks and treated with mifepristone or DMSO vehicle in culture medium for 24 hours as described above. Adherent cells were collected by trypsinization, and non-adherent (floating) cells were collected from culture medium. Cells were centrifuged, washed in PBS, and incubated in 1X annexin binding buffer, annexin V-fluorescein isothiocyanate (FITC), and propidium iodide (PI) according to the manufacturer's recommendations. Stained cells were immediately analyzed using a Cytomics FC500 flow cytometer. Data were analyzed using CXP Analysis software (version 2.2, Coulter).

Cell Cycle Analysis.
Harvested HEI-193 cells were washed in PBS and fixed in cold 70% ethanol at −20 °C for 18-72 h. Before staining with propidium iodide (PI), cells were centrifuged and washed again in cold PBS. 2 × 10 6 or fewer cells were incubated with 500 µl staining solution (0.1% Triton X-100 (Sigma), 2 mg/ml RNase A (Qiagen), and 1 µg/ml PI (Miltenyi Biotec) in PBS) for 15 minutes at 37 °C. In order to exclude DNAse activity, RNase A was boiled for 5 minutes and cooled down before its addition to staining solution. Cells were analyzed on a Cytomics FC500 flow cytometer using CXP Analysis software (version 2.2, Coulter).
Enzyme-linked Immunosorbent Assay. Cell-conditioned medium was collected from mifepristone-treated and untreated HEI-193 cells after 72-hour incubation with the drug. Enzyme-linked immunosorbent assay (ELISA) was performed on each sample in triplicate to assess the quantity of progesterone in cell-conditioned medium, according to the manufacturer's protocol (Enzo Life Science, #ADI-900-011).
Cytotoxicity and Cell Confluence Assays. Cytotoxicity and cell confluence were measured using live-cell, time-lapse phase contrast and fluorescence imaging acquired at 10X by an IncuCyte S3 instrument (Essen Bioscience). Cytotoxicity was measured by incorporation of IncuCyte Cytotox Reagent (Essen Bioscience), applied according to manufacturer's instructions; the reagent fluoresces when it binds DNA after compromise of membrane integrity. Nine images were acquired per well every 2 hours for 72 hours, and cytotoxicity at each time point is reported as number of fluorescent objects per well, after thresholding to avoid the inclusion of small cellular debris. Phase object confluence was measured using time-lapse phase contrast imaging acquired at 10X by the same instrument, analyzing 9 images per well every 2 hours for 72 hours, and is reported as percent confluence per square millimeter of well space. Ingenuity Pathway Analysis. Ingenuity Pathway Analysis software (Qiagen) was used to perform standard Core Analysis on all genes in our meta-analysis that reached significance (p < 0.05) after Bonferroni correction for multiple hypothesis testing. Relevant upstream regulators for the resulting networks were identified and analyzed using publicly available Ingenuity Pathway Analysis documentation (Qiagen). gDNA Extraction. Genomic DNA (gDNA) was extracted from six vestibular schwannoma tissue samples using the DNeasy Blood and Tissue Kit (Qiagen) following manufacturer's specifications. The concentration of double-stranded DNA in each sample was evaluated using a Qubit dsDNA BR Assay Kit. A minimum measurement of 50 ng/μl was required for each sample to be included with HaloPlex target enrichment.

Library Preparation and Targeted Capture.
A library of DNA restriction fragments from all coding exons, introns, and UTRs (5′ and 3′) of the NF2 gene was prepared using a HaloPlex HS target enrichment kit (Agilent Technologies), following the manufacturer's instructions. The total region size was 95.045 kbp with an actual analyzed target of 89.408 kbp bases, which required 2581 amplicons to achieve this 94.07% coverage with maximum validation stringency. MiSeq instrument by using a MiSeq Reagent Kit v2 flow cell (Illumina). Raw data were processed and variants prioritized according to best practices of the core facility (Supplementary Methods).
Statistical analysis. Throughout this paper, though figures present metabolic activity and cellular proliferation data as percentage of vehicle-treated control, all statistical analyses were performed on raw data, in accordance with good statistical practice in pharmacology 68 . Specifically, in Figs 2a and f, 4a,b and g, raw vehicle-treated and mifepristone-treated cell data are compared using two-way ("randomized block") ANOVA, selected to minimize within-experiment variation by "blocking" treatment data with control data while meeting the equal-variance assumptions required by ANOVA 68 . In Fig. 2a, comparison between control group and mifepristone-treated group (35 μM) was conducted using randomized block ANOVA on mean optical density values per treatment condition measured in eight independent experiments (p = 0.006, F = 42.46, DF = 1). For Fig. 2f, randomized block ANOVA was performed between the ratio of BrdU+ to Hoescht+ cells per treatment condition in five independent experiments (p = 0.0007, F = 88.25, DF = 1). For Fig. 4a, comparison between control group and mifepristone-treated group (35 μM) was conducted using randomized block ANOVA on mean optical density values per treatment condition measured in ten independent experiments (p = 0.002, F = 23.02, DF = 1). In Fig. 4b, randomized block ANOVA was performed between the ratio of BrdU+ to Hoescht+ cells per treatment condition in seven independent experiments (p = 0.0002, F = 68.47, DF = 1). In Fig. 4g, randomized block ANOVA was performed on mean optical density values per treatment condition measured in eight independent experiments (p = 0.230, F = 1.255, DF = 1).

Data availability
The results of the 80-sample combined VS meta-analysis, 16