Combined MEK and JAK/STAT3 pathway inhibition effectively decreases SHH medulloblastoma tumor progression

Zagozewski, Jamie; Borlase, Stephanie; Guppy, Brent J.; Coudière-Morrison, Ludivine; Shahriary, Ghazaleh M.; Gordon, Victor; Liang, Lisa; Cheng, Stephen; Porter, Christopher J.; Kelley, Rhonda; Hawkins, Cynthia; Chan, Jennifer A.; Liang, Yan; Gong, Jingjing; Nör, Carolina; Saulnier, Olivier; Wechsler-Reya, Robert J.; Ramaswamy, Vijay; Werbowetski-Ogilvie, Tamra E.

doi:10.1038/s42003-022-03654-9

Download PDF

Article
Open access
Published: 14 July 2022

Combined MEK and JAK/STAT3 pathway inhibition effectively decreases SHH medulloblastoma tumor progression

Communications Biology volume 5, Article number: 697 (2022) Cite this article

2876 Accesses
6 Citations
15 Altmetric
Metrics details

Subjects

Abstract

Medulloblastoma (MB) is the most common primary malignant pediatric brain cancer. We recently identified novel roles for the MEK/MAPK pathway in regulating human Sonic Hedgehog (SHH) MB tumorigenesis. The MEK inhibitor, selumetinib, decreased SHH MB growth while extending survival in mouse models. However, the treated mice ultimately succumbed to disease progression. Here, we perform RNA sequencing on selumetinib-treated orthotopic xenografts to identify molecular pathways that compensate for MEK inhibition specifically in vivo. Notably, the JAK/STAT3 pathway exhibits increased activation in selumetinib-treated tumors. The combination of selumetinib and the JAK/STAT3 pathway inhibitor, pacritinib, further reduces growth in two xenograft models and also enhances survival. Multiplex spatial profiling of proteins in drug-treated xenografts reveals shifted molecular dependencies and compensatory changes following combination drug treatment. Our study warrants further investigation into MEK and JAK/STAT3 inhibition as a novel combinatory therapeutic strategy for SHH MB.

Cross-species analysis of SHH medulloblastoma models reveals significant inhibitory effects of trametinib on tumor progression

Article Open access 19 September 2023

Wnt activation as a therapeutic strategy in medulloblastoma

Article Open access 28 August 2020

EZH2 targeting reduces medulloblastoma growth through epigenetic reactivation of the BAI1/p53 tumor suppressor pathway

Article 03 October 2019

Introduction

Brain tumors are the deadliest form of childhood cancer and account for 20% of all new pediatric cancer cases. Medulloblastoma (MB) is the most common malignant primary pediatric brain tumor. Treatment currently consists of aggressive surgery, high doses of cytotoxic chemotherapy, and radiation to the whole brain and spinal cord. Unfortunately, 40% of patients still succumb to their disease, while survivors must deal with the long-term toxicities associated with chemotherapy and radiation^1,2,3.

MB consists of at least 4 highly distinct molecular subgroups: WNT, Sonic Hedgehog (SHH), Group 3, and Group 4^1,4. SHH MBs are characterized by activation of the SHH pathway, and include very high-risk groups of both children (>3 years old) and infants (<3 years old) that exhibit significant intertumoral heterogeneity and account for the majority of treatment failures despite aggressive therapy^{5,6,7,8,9,10,11}. For example, TP53 mutations confer a near universally fatal prognosis in older children with SHH MB tumors and constitute a very high-risk form of the disease⁷. Personalized therapies for SHH MB are lacking, as SHH pathway antagonists are not predicted to work in younger patients with tumors harboring mutations in downstream SHH pathway genes^12,13,14,15. Moreover, premature osseous fusion¹⁶ has limited the clinical utility of SHH pathway antagonists. Thus, novel targeted therapies are urgently needed.

A major contributing factor to MB cellular heterogeneity is the presence of putative brain tumor-propagating cells (TPC) or cancer stem/progenitor cells^17,18,19,20. These primitive cells are believed to underlie progression and drug resistance in several cancers, including MB^19,21,22. Therefore, the ability to identify robust TPC biomarkers and characterize their role in regulating MB tumorigenesis may reveal novel subgroup-specific therapies capable of targeting TPC subpopulations. We recently identified novel roles for the CD271/p75 neurotrophin receptor (CD271/p75NTR) and the MEK/ERK signaling pathway in contributing to SHH MB TPC growth and tumor progression^23,24. Bioinformatics analyses of large patient sample datasets and tumorspheres from SHH MB cultures²⁴ demonstrated that CD271 is a potential diagnostic marker for these tumors. CD271+ cells exhibit upregulated MEK/ERK signaling and inhibiting this pathway reduced endogenous CD271 levels, stem cell proliferation, survival, and migration in vitro. In addition, CD271 knockdown decreased cell proliferation and pERK activity^23,24. Importantly, treatment with the brain penetrant MEK inhibitor selumetinib extended survival and decreased CD271 levels in a biologically relevant preclinical animal model²⁴. Combined with pERK staining in primary SHH MB tumor samples, our previous in vitro and in vivo findings demonstrated that the MEK/ERK pathway is a therapeutic target in human SHH MB. While these data were promising, selumetinib-treated mice ultimately succumbed to disease progression.

Thus, we sought to enhance the efficacy of selumetinib by targeting pathways that compensate for MEK inhibition in vivo with the overarching goal of further enhancing survival in SHH MB tumors. Here, using a combination of transcriptomics, proteomics, multiplex spatial profiling as well as in vitro and in vivo functional assays, we show that JAK/STAT3 activity is upregulated following selumetinib treatment in SHH MB tumor xenografts and that targeting this pathway, in combination with MEK inhibition, significantly attenuates tumor growth. Dual inhibition of MEK and JAK/STAT3 pathway activity abrogates tumorigenic properties in three different MB cell models in vitro as well as in 2 in vivo models of SHH MB. Thus, combinatory JAK/STAT3 and MEK pathway inhibition may serve as a novel therapeutic strategy for future SHH MB patients.

Results

JAK/STAT3 pathway activity is upregulated following MEK inhibition in vivo

Our previous work demonstrated that selumetinib increased survival of NOD SCID mice bearing orthotopic SHH MB xenografts compared to vehicle controls²⁴. However, the animals eventually succumbed to tumor progression²⁴; thus, we sought to identify the altered molecular pathways that may compensate for MEK inhibition in vivo. We posit that these pathways can be therapeutically exploited to identify an optimized combinatorial therapeutic strategy to further decrease tumor growth.

To identify aberrant molecular changes following selumetinib treatment, we performed RNA sequencing (RNA-seq) on sorted control (N = 3 independent endpoint tumors at days 43 and 46 (2 mice) post-tumor cell injection) and selumetinib-treated (N = 3 independent endpoint tumors at days 56, 58 and 59 post-tumor cell injection) UI226 SHH MB cells (recently derived from a primary SHH MB tumor^23,25) from tumor xenografts previously characterized²⁴ (Fig. 1a, b). Endpoint tumors from both the control and selumetinib groups were histologically similar displaying high nuclear-to-cytoplasmic ratios and well-encapsulated margins (Fig. 1b). We chose in vivo characterization over in vitro screen-based analyses, as this will identify pathways that may compensate for MEK inhibition in the more complex milieu of the cerebellar tumor microenvironment. RNA-seq revealed 576 significantly (p < 0.05) differentially expressed genes, the majority (79%) of which were upregulated (Fig. 1c, d, Supplementary Fig. 1a, and Supplementary Data 1) in selumetinib-treated cells. Genes associated with JAK/STAT3, TNFα /NFκB, as well as apoptotic signaling pathways were enriched in gene sets that were upregulated in selumetinib-treated xenografts (Fig. 1e and Supplementary Data 2–7). Importantly, genes that show decreased expression upon MEK activation were also enriched in upregulated gene sets in selumetinib-treated tumors, consistent with sustained decreased pathway activity. A total of 119 genes exhibited decreased expression in selumetinib-treated tumors (Fig. 1f). Of note, CD271 was the 3^rd most significantly downregulated gene (5.2-fold) providing additional validation of our previous findings²⁴.

**Fig. 1: Pathways and genes differentially expressed in selumetinib-treated xenografts.**

Very little is known about the role of the JAK/STAT3 pathway in human SHH MB tumor progression and drug resistance. Several drugs targeting JAK/STAT3 signaling have been identified and are known to be brain penetrant in xenograft models including pacritinib²⁶ and AZD1480²⁷. Recent studies have also shown that STAT signaling activation accompanies MAPK pathway inhibition in melanoma^28,29, thyroid cancer³⁰, and breast cancer³¹, thus implicating upregulated JAK/STAT3 activity as a general mechanism associated with drug resistance or long-term adaptive responses in multiple cancers. To further explore this pathway in our model system, we performed pSTAT3 (Tyr705) immunohistochemistry (IHC) on formalin-fixed paraffin-embedded (FFPE) tissue sections derived from control and selumetinib treated UI226 SHH MB xenografts. STAT3 is phosphorylated upon upstream JAK activation and is an appropriate readout of activated JAK/STAT3 signaling^32,33. In line with the RNA-seq results, IHC revealed a significant 2.5-fold increase in the number of pSTAT3+ cells in selumetinib-treated xenografts relative to controls (Fig. 1g, h).

While we focused on in vivo molecular compensations following MEK inhibition, nevertheless, we examined pSTAT3 expression in SHH MB tumorspheres treated with selumetinib for 3 days in culture. pSTAT3 levels decrease with increasing selumetinib concentration, concomitant with an upregulation in apoptosis (Supplementary Fig. 1b). Given that our in vivo experiments were carried out long-term (30–40+ days treatment) in the more complex cerebellar tumor microenvironment, it is unsurprising that our short-term in vitro analysis did not recapitulate these findings. Therefore, we next cultured SHH MB tumorspheres over multiple passages using a lower, less toxic concentration of selumetinib to examine pSTAT3 levels following extended treatment. While Daoy tumorspheres are amenable to this long-term exposure and maintained over multiple passages, UI226 tumorspheres are not (Supplementary Fig. 1 c). Indeed, under these conditions, Daoy tumorspheres also exhibit an increase in pSTAT3 levels following MEK inhibition (Supplementary Fig. 1 d).

As lower pSTAT3 levels were observed in the UI226 xenograft vehicle controls in vivo, we next evaluated JAK/STAT pathway gene expression in primary MB patient samples. Analyses of Hallmark JAK/STAT signature genes³⁴ in a dataset consisting of 218 tumors across 7 types of pediatric brain cancers³⁵ demonstrate that the majority of JAK/STAT-related transcripts and proteins are lower across the 22 MB samples (Fig. 2a). However, the known negative regulator of JAK/STAT signaling, protein tyrosine phosphatase non-receptor type 2 (PTPN2), is strikingly higher across the MB subset (Fig. 2a). More specific hierarchical clustering selecting genes restricted to JAK/STAT pathway gene sets across an independent cohort of 763 MB patient samples stratified MB tumors into the SHH subgroup (Fig. 2b and Supplementary Fig. 2). These transcriptome signatures are significantly lower in SHH MBs relative to the other subgroups (Fig. 2c). Collectively, these data reveal low basal JAK/STAT pathway activity in both vehicle control xenografts and primary patient MB samples with the lowest levels in SHH MB tumors. Therefore, JAK/STAT3 signaling activation is a selumetinib-induced response in our UI226 SHH MB model.

**Fig. 2: The JAK/STAT3 signaling signature is significantly lower in SHH MB primary patient tumors.**

Treatment with JAK/STAT3 pathway inhibitors significantly inhibits tumor properties in vitro

To initially test the effect of inhibiting the JAK/STAT3 pathway on tumor properties, we developed a high-content 3D tumorsphere drug-screening assay to construct high-resolution concentration-response curves following treatment with the individual JAK/STAT3 inhibitors pacritinib and AZD1480 in vitro (Fig. 3a). Serum-free tumorsphere assays were specifically employed, as pSTAT3 is expressed in these stem cell-enriched, defined culture conditions. pSTAT3 is markedly decreased following treatment with both JAK/STAT3 inhibitors (Fig. 3b–d and Supplementary Fig. 3 a–c). Similar decreases in pSTAT3 levels were observed in tumorspheres from two different human cell lines, UI226 and Daoy³⁶, as well as Ptch+/−:p53+/− SHH mouse MB cultures derived from primary tumors. These mouse SHH MB cells are amenable to tumorsphere culture and retain SHH signaling in either stem-cell enriched or modified stem-cell conditions^37,38. We therefore utilized the tumorsphere assay to systematically score thousands of tumorspheres at various drug concentrations over 5 days using sphere size as a readout. To avoid inclusion of dead cells and aggregates, only true tumorspheres with diameters between 45 and 360 μm and circularities greater than 0.3 were scored. The results (Fig. 3e–f) were then used to determine the optimal μM range of concentrations for use in additional cellular assays. For comparison, we also evaluated the chemotherapeutic agent vincristine, which is quite potent in the nM range, in our tumorsphere model (Supplementary Fig. 3 d).

**Fig. 3: JAK/STAT3 inhibitors decrease tumorsphere size and cell migration in vitro.**

The inhibitory effects of pacritinib and AZD1480 on tumorsphere size were independently validated by cumulative frequency distribution analyses using Fiji/ImageJ for both UI226 and Daoy (Fig. 3g, h and Supplementary Fig. 3 e, f). At the highest doses, we observed reductions in tumorsphere number for all 3 cell models when treating with pacritinib and/or AZD1480 (Supplementary Fig. 3 g–l). However, at these concentrations, there were no statistically significant reductions in cell viability (Supplementary Fig. 3 m–r). When all 3 cell models were cultured as aggregates in 3D collagen to assess migration, significant decreases in cell motility were observed following treatment with both drugs (Fig. 3i–n). Collectively, these results show an overall decrease in tumorigenic properties of MB cells including migration, tumorsphere size, and number following treatment with single JAK/STAT3 inhibitors in vitro.

Dual MEK and JAK/STAT3 pathway inhibition elicits a further reduction in SHH MB tumorigenic properties in vitro

Next, we tested the hypothesis that simultaneous JAK/STAT3 + MEK inhibition would yield a further reduction in tumorsphere size in our high-content 3D tumorsphere assay. To evaluate in vitro drug synergy, we generated drug combination matrices in 96-well plates such that each well contained a unique combination of AZD1480 + selumetinib or pacritinib + selumetinib. Both JAK/STAT3 inhibitors produced a synergistic reduction in tumorsphere size when combined with selumetinib. Specifically, a concentration of 2.5 μM selumetinib in combination with 39 nM or 156 nM pacritinib elicited a synergistic reduction in both UI226 and Daoy tumorsphere size (Fig. 4a, b). When combining selumetinib with AZD1480, drug synergies were also observed at several concentrations (Fig. 4a, b). Additional reductions in tumorsphere size were independently validated by cumulative frequency distribution analyses using Fiji/ImageJ for both UI226 and Daoy cells (Fig. 4c, d).

**Fig. 4: Combination of MEK + JAK/STAT3 inhibition significantly inhibits SHH MB tumorsphere size, and cell migration.**

Having established a decrease in tumorsphere size following MEK + JAK/STAT3 inhibitor treatment, we sought to determine if loss of cell viability, as assessed by Trypan blue exclusion assays, could account for the smaller MB tumorspheres (Supplementary Fig. 4 a–c). A significant reduction in cell viability was observed following combination pacritinib + selumetinib treatment for both UI226 and Daoy tumorspheres (Supplementary Fig. 4 a, b). AZD1480 + selumetinib also significantly reduced cell viability for UI226 but not for Daoy tumorspheres (Supplementary Fig. 4 a, b). Similarly, Annexin V staining of both pacritinib + selumetinib and AZD1480 + selumetinib treated UI226 tumorspheres revealed significant increases in the proportion of dead cells (Annexin V+/7AAD+) compared to single-agent treatment (Supplementary Fig. 4 d, e). However, there was no change in Ptch+/−:p53+/− SHH MB cell viability following combination drug treatment (Supplementary Fig. 4 c) suggesting that MEK + JAK/STAT3 inhibition results in either cytostatic or a combination of cytostatic/cytotoxic effects on MB tumorspheres.

Finally, we evaluated the effect of dual MEK + JAK/STAT3 inhibitor treatment on migration. Statistically significant further reductions in migration compared to single agents and/or DMSO controls were observed for UI226, Daoy and Ptch+/−:p53+/− MB cells (Fig. 4e–h). Taken together, our in vitro assays demonstrate that dual MEK + JAK/STAT3 inhibition results in further reductions in SHH MB cell tumorigenic properties in vitro.

RNA sequencing reveals significant molecular changes in SHH MB tumorspheres following MEK inhibition or dual MEK and JAK/STAT3 pathway inhibition

As dual JAK/STAT3 + MEK inhibition reduced tumorsphere size and cell migration, we next performed RNA-seq on selumetinib and/or pacritinib treated UI226 tumorspheres to characterize the sustained transcriptomic changes 5 days following treatment in vitro. Interestingly, the pacritinib-treated tumorspheres at day 5 resembled the vehicle controls with only 9 sustained significantly differentially expressed genes (Fig. 5a). This suggests that the molecular changes associated with pacritinib alone are short-lived. In contrast, one-time selumetinib treatment was sufficient to drive the major transcriptomic changes at day 5 with 6628 genes exhibiting significant differential expression (FDR < 0.05) (Fig. 5b). In line with these findings, both selumetinib and combination drug treatment resulted in similar changes to hallmark gene sets, canonical pathways and oncogenic signatures (Fig. 5b–g). Genes/pathways associated with cell cycle and DNA replication/repair were the most significantly enriched downregulated gene sets following selumetinib or combination drug treatment (Fig. 5d). Strikingly, these pathways have also been shown to be strongly enriched specifically in high-risk SHH MB patient tumors exhibiting TP53 mutations^5,7. In addition, genes/pathways associated with neurogenesis and epithelial to mesenchymal transition (EMT) were significantly enriched among genes that were downregulated following selumetinib or combination drug treatment (Fig. 5d, e).

**Fig. 5: MEK or MEK + JAK/STAT3 inhibition significantly alters the transcriptome of UI226 SHH MB-treated tumorspheres.**

Further exploring the transcriptomic changes related to cell fate, we observed that gene sets associated with granule neuron progenitor cell proliferation after SHH stimulation were enriched for genes that were downregulated following selumetinib or combination drug treatment (Fig. 5f, upper). While most genes were not specifically associated with SHH signaling, significant changes in PTCH1, SUFU and GLI3 were observed (Fig. 5f, lower). In contrast, gene sets associated with NOTCH signaling, astrocytes, and neural crest differentiation (Fig. 5d, g) were significantly enriched among genes that were upregulated following combination drug treatment and/or selumetinib alone, likely reflecting the sustained, but indirect, effects after 5 days. Taken together, these results demonstrate that the transcriptome changes at day 5, in vitro, are mainly driven by the sustained effects of selumetinib treatment.

Dual MEK + JAK/STAT3 pathway inhibition significantly improves survival and reduces tumor growth in vivo

Next, we examined whether combinatorial drug treatment could result in a further reduction of tumor growth in vivo. As pacritinib has been shown to be effective at much lower concentrations, we focused on this drug in combination with selumetinib. Initial pilot studies with UI226 were completed to determine the maximum tolerated dose and revealed that 25, 50, and 100 mg/kg of pacritinib were all well-tolerated with no toxicities when administered via oral gavage twice daily, 5 days a week (Fig. 6a). Pacritinib treatment alone did not extend survival in vivo (Fig. 6b). Given the overall lower levels of pSTAT3 specifically in the UI226 vehicle controls (Fig. 1g, h), this was not surprising. Based on these findings, 50 mg/kg pacritinib was chosen for combination therapy in vivo along with a reduced dose of 37.5 mg/kg selumetinib (compared to 75 mg/kg selumetinib in Fig. 1 and our previous work²⁴). Following tumor engraftment, mice were administered vehicle, 50 mg/kg pacritinib, 37.5 mg/kg selumetinib, or a combination of 50 mg/kg pacritinib and 37.5 mg/kg selumetinib twice daily, 5 days a week until endpoint was reached. This treatment schedule resulted in highly favorable toxicity profiles, as weights steadily increased, even with pacritinib + selumetinib treatment (Fig. 6c). Importantly, the combination of pacritinib + selumetinib significantly improved survival in our UI226 xenograft model relative to vehicle control and each drug alone (Fig. 6d). In contrast, no significant improvements in survival were observed when HDMB03 Group 3 MB xenografts were treated with selumetinib and/or pacritinib relative to the vehicle (Fig. 6e). To compare differences in tumor growth across treatment groups, a subset of UI226 vehicle and treated xenografts was also extracted prior to endpoint, all at day 41, and samples were stained for the human specific marker STEM121 (Fig. 7a, b). While a small decrease in tumor size was observed for single agent therapies, combination drug treatment resulted in a substantial and significant reduction in tumor growth (Fig. 7a–c). Similarly, combination drug treatment of RCMB18 SHH MB patient-derived xenografts (PDX) in a small cohort of NOD SCID IL2Rg null (NSG) mice significantly reduced tumor growth, as single well-encapsulated nodules were observed relative to vehicle controls (Fig. 7d, e). Interestingly, this growth pattern was observed in the presence of pacritinib, with or without selumetinib. The tumors typically appeared as single well-encapsulated nodules as opposed to the vehicles and selumetinib-only samples which consisted of multiple infiltrating masses (Fig. 7d). Collectively, these data demonstrate that pacritinib + selumetinib treatment attenuates tumor progression in two biologically relevant SHH MB in vivo models.

**Fig. 6: MEK + JAK/STAT3 inhibition increases survival of NOD SCID mice transplanted with UI226 SHH MB cells.**

**Fig. 7: MEK + JAK/STAT3 inhibition decreases tumor growth in 2 SHH MB in vivo xenograft models.**

Digital spatial profiling reveals unique protein changes following in vivo treatment

To more comprehensively and quantitatively explore the molecular changes following combination drug treatment, we chose to employ multiplex digital spatial profiling, an approach that combines protein abundance and spatial distribution data using antibodies coupled to photocleavable oligonucleotide tags³⁹. Once unstained FFPE slides are labeled, the oligonucleotides are released from select regions of interest (ROIs) by UV exposure and quantified³⁹. Identification of protein/phosho-protein differences between samples provided immediate insight into the long-term effects of our inhibitors. Using this method, we profiled 56 proteins, including a neural cells core panel, a MAPK, a PI3K as well as a cell death protein panel, across 12 ROIs for representative vehicle, pacritinib, selumetinib, and combination therapy (pacritinib + selumetinib) UI226 FFPE sections collected prior to endpoint. For each sample, 6 ROIs in the tumor core and 6 ROIs along the tumor/mouse brain boundary were interrogated (Fig. 8a–d). CD271, MAP2, and Ki67 were used to visualize the tumor regions by immunofluorescence (Fig. 8a–d). Quantification of CD271 levels across all 12 ROIs/sample showed no significant changes following pacritinib treatment (Fig. 8e, f). However, a significant downregulation was observed after selumetinib treatment and a further reduction was seen with combination therapy (Fig. 8e, f). These results support our current (Fig. 1) and previous findings demonstrating decreased CD271 levels following selumetinib treatment²⁴ across platforms.

**Fig. 8: Multiplex digital spatial profiling of 56 proteins on MEK and/or JAK/STAT3 inhibitor-treated UI226 tumors in vivo.**

While our tumor samples were positive for the human cell marker STEM121 (Fig. 7a), we cannot exclude the possibility that there are mouse cells infiltrating the tumor, particularly along the border. As a significant portion of the 56 antibodies cross-react with mouse antigens, we next compared the 56 chosen proteins in the tumor core vs. the border across all 4 samples. No significant differences were observed for the vehicle, pacritinib, and combination-treated samples, while only 3 proteins showed a significant difference in the selumetinib sample (Supplementary Fig. 5 a–d). This was in stark contrast to our recently published studies for Group 3 MB in which mTOR inhibitor treatment induced neural differentiation regionally along the tumor border (Supplementary Fig. 5 e–h)⁴⁰. As we did not observe these regional differences in the current study, all 12 ROIs/sample were combined for additional analyses.

Multiplex spatial profiling revealed significant differences between the vehicle and inhibitor-treated samples (Fig. 8g–s and Supplementary Data 8–10). As expected, the selumetinib-treated tumor exhibited significant changes in proteins associated with the MAPK pathway (Fig. 8h). This includes downregulated pMEK1, total ERK1/2, and BRAF (Fig. 8j–k, m) but not pERK1/2 (Fig. 8l). There were no consistent and significant increases in differentiation proteins following treatment (Fig. 8g–i, p, q), but a reduction in OLIG2 was observed following selumetinib treatment alone (Fig. 8p). OLIG2 + progenitor cells have recently been shown to represent a putative tumor-initiating cell population in mouse models of SHH MB and are also enriched in human SHH MB stem-like cells⁴¹. Interestingly, sustained decreases in pMEK, total ERK, and BRAF were not observed following combination drug treatment; however, there was a significant reduction in pERK as well as significant upregulation of neurofibromin 1 (NF1), a tumor suppressor that reduces Ras activity and thus MAPK signaling (Fig. 8j–n). Of note, p53 was most significantly upregulated following combination treatment (Fig. 8o). In the pacritinib-treated tumor, no changes in pMEK, total ERK, and BRAF were observed; however, pERK was downregulated and NF1 was significantly upregulated (Fig. 8j–n), likely reflecting pathway crosstalk. Importantly, across all treatment groups, there were significant decreases in proteins associated with PI3K/AKT signaling such as PLCG1 and pAKT1 (Fig. 8r–s), again suggesting suppression through pathway crosstalk. As CD271 levels were sustained across several ROIs following pacritinib treatment (Fig. 8e, f), the data for this sample were also re-analyzed to compare protein levels between CD271^high and CD271^low ROIs; however, there were no significant differences between these groups for the proteins evaluated (Supplementary Fig. 6). Collectively, these results demonstrate the value of multiplex protein profiling in rapidly characterizing drug responses in xenograft models of SHH MB. While all 3 drug-treated samples exhibited significant alterations in MAPK and PI3K/AKT-associated proteins, differences at the individual protein level between samples and changes in molecular dependencies following treatment reveal compensatory changes in residual cells to sustain tumor growth suppression.

Discussion

Current large-scale sequencing efforts rely on sampling from bulk SHH MB tumors. While these technologies have revolutionized our understanding of intertumoral heterogeneity between the subgroups and have improved patient stratification, they do not take into account the extensive intratumoral heterogeneity, namely the putative stem cells, progenitors, and more differentiated phenotypes observed in these tumors. By deconstructing cellular heterogeneity and the mechanisms regulating self-renewal and differentiation, we recently revealed novel molecular vulnerabilities in specific subpopulations (ie. CD271+ cells) that could be exploited using brain penetrant small molecule inhibitors²⁴. Combined with pERK staining in primary patient samples, this led to the demonstration that the MEK/ERK pathway is a novel therapeutic target in human SHH MB models. Indeed, previous studies have shown that MEK/ERK signaling drives SHH pathway inhibitor resistance and promotes metastasis specifically in SHH MB mouse models³⁸. However, single-agent therapies typically fail due to the upregulation of compensatory pathways. Thus, we performed RNA-seq on selumetinib-treated xenograft tumors in lieu of in vitro drug screens to identify adaptive responses in a more biologically relevant model system. We discovered that JAK/STAT3 pathway activity is elevated following single agent selumetinib therapy leading to a comprehensive examination of combinatorial MEK + JAK/STAT3 inhibitor treatment. Our in vitro findings demonstrate significant reductions in tumorigenic properties of SHH MB, including tumorsphere size and migration, with combination MEK + JAK/STAT3 inhibitor treatment compared to controls and monotherapies. Importantly, pacritinib + selumetinib treatment resulted in a significant reduction in tumor burden in 2 in vivo models suggesting that combination treatment could potentially be used as a tumor debulking strategy for SHH MB patients.

Recently, pacritinib has been examined for its efficacy in the treatment of glioblastoma (GBM)²⁶. While the reduction in tumorigenic properties was observed in vitro, pacritinib alone was insufficient to increase the survival of orthotopic xenograft models in vivo, further underscoring the need for combination-based therapeutic approaches. Indeed, combination treatments with pacritinib and the tyrosine kinase inhibitor, afatinib, have also recently been evaluated in GBM⁴². Toxicity issues inhibited evaluation of survival in the combination treatment group in vivo; however, a significant reduction in tumor burden was observed following pacritinib + afatinib treatment compared to single-agent therapies and controls⁴². Of note, for the combination therapy studies, we chose 50 mg/kg pacritinib as opposed to the 100 mg/kg utilized in previous brain tumor studies^26,42. We believe that this lower dose, combined with the 5-day on, 2-day off drug holiday treatment schedule, helped to mitigate long-term toxicities in our model. However, similar to the previous studies^26,42, pacritinib alone did not improve survival in our animal models. Interestingly, Jensen et al.²⁶ demonstrated that pacritinib was rapidly metabolized by mouse but not human liver microsomes. This rapid clearance, combined with lower overall JAK/STAT3 pathway activity in vehicle controls, may account for the lack of efficacy as a single agent in the xenografts.

Resistance to MEK inhibition due to the upregulation of STAT3 signaling has been well documented in a number of cancers including melanoma^28,29, low-grade gliomas⁴³, colon⁴⁴, pancreatic⁴⁵, thyroid³⁰, lung⁴⁶, and breast cancer³¹. Consequently, combination-based therapeutic approaches utilizing both MEK and STAT3 inhibitors have been tested and led to reduced tumor growth and increased survival in preclinical animal models^44,45,46. Selumetinib has undergone extensive testing for the treatment of refractory low-grade gliomas^47,48 as well as plexiform neurofibromas^49,50, and has shown very good efficacy as well as excellent tolerability in both conditions. The MEK1/2 inhibitor trametinib is approved for adult BRAF V600E mutant cancers (ie. NCT01682083) and is also currently undergoing testing in both refractory low-grade gliomas and plexiform neurofibromas (NCT03363217, NCT02124772)^51,52. Both MEK inhibitors have excellent toxicity profiles in children. However, our results as well as previous studies in other cancers do raise the possibility that upregulation of other signaling pathways, particularly JAK/STAT3 activation, could potentially compensate for MEK inhibition following treatment of both diseases. Results of a recent retrospective centrally reviewed multi-center study underscore the need to further understand the molecular mechanisms leading to rapid tumor progression in some patients following MEK inhibitor treatment⁵³. Thus, more pre-clinical studies with additional brain tumor models are necessary to determine the appropriate timing for drug delivery and to more thoroughly assess sustained responses long-term. Interestingly, recent work in Group 3 MB demonstrated increased JAK/STAT3 pathway activity in chemoresistant cell lines⁵⁴. While pacritinib treatment did not improve survival in our Group 3 MB xenograft model, this does not exclude the possibility that treatment of chemoresistant Group 3 MB cells with JAK/STAT3 inhibitors could further improve survival in vivo.

In addition to JAK/STAT3 signaling, our RNA-seq data also suggested that the TNFα /NFκB pathway could also compensate for selumetinib treatment in SHH MB xenografts. While less widely studied than STAT3 compensation, elevated expression of TNFα has also been shown to confer resistance in response to MEK inhibition. For example, MEK inhibitor resistance was observed in melanoma due to the accumulation of macrophage-derived TNFα expression and consequent downstream upregulation of the transcription factor MITF, which plays an oncogenic role in melanoma⁵⁵. Elevation of TNFα expression in response to MEK inhibition has also been observed in lung cancer⁵⁶. Finally, NFκB inhibition was shown to reduce MB cell growth in multiple cell lines⁵⁷. These studies, along with our RNA-seq data, indicate that MEK inhibition in combination with TNFα /NFκB pathway inhibition should also be tested in SHH MB models.

The utilization of multiplex spatial profiling to assess the impact of treatment on xenograft models opens up new opportunities to comprehensively assess adaptive mechanisms and sustained drug responsiveness in situ. However, one limitation is the restricted number of relevant and validated antibody panels available for the spatial proteomics analyses. We were able to evaluate changes in proteins associated with neural differentiation, cell death, MAPK activity, and PI3K activity following long-term combination drug treatment in vivo. While the data enabled an initial assessment of some molecular changes, a direct comparison to the RNA-seq dataset from 5-day treated tumorspheres was difficult. Discrepancies between the short-term in vitro and long-term in vivo datasets are expected; however, future studies will expand the spatial profiling beyond proteomics to include comprehensive transcriptomic analyses to assess overlap.

We have shown that the combined inhibition of MEK and JAK/STAT3 pathway activity abrogates tumor properties in 3 different MB cell models in vitro as well as in 2 in vivo models of SHH MB. Our study reveals new insight into the altered molecular mechanisms associated with selumetinib treatment and warrants further investigation into dual MEK and JAK/STAT3 inhibition as a novel combinatory therapeutic strategy for SHH MB.

Methods

Cell culture and drug treatment

UI226 (SHH MB as analyzed by NanoString⁵⁸) cells were a kind gift from Dr. Timothy Ryken (Dartmouth-Hitchcock Medical Center, New Hampshire, USA) and have been adapted to cell culture in StemPro media as described^23,24,25. Daoy cells were purchased from ATCC and cultured as described previously^23,24. The HDMB03 Group 3 MB cell line was obtained from Dr. Till Milde⁵⁹. All cell lines have been authenticated by STR profiling (ATCC). Tumorspheres were also propagated directly from Ptch+/−:p53+/− mouse SHH MB tumors⁶⁰. Tumors were dissociated and Ptch+/−:p53+/− cells were grown and expanded in a NeuroCult Basal Medium (Mouse & Rat) (Stem Cell Technologies) supplemented with 7.5% BSA (1000X), 100X GlutaMAX (ThermoFisher, Scientific) in a final working solution consisting of 100X N2, 50X B27 (ThermoFisher Scientific), 10 ng/ml EGF (Sigma), 10 ng/ml bFGF (Sigma), and 2 μg/ml heparin (Sigma) as previously described³⁷. Tumorsphere assays were carried out in stem cell enriched conditions as described previously^23,24. Briefly, UI226, Daoy and Ptch+/−:p53+/− mouse SHH MB cells were plated in 24-well ultra-low attachment plates at optimal densities of 2500 cells/well (UI226), 10,000 cells/well (Ptch+/−:p53+/− cells) and 10,000 cells/well (Daoy) and then grown for 5 days at 37 °C/ 5% CO₂. Pacritinib was added at 10, 25, 50, 100, 250 nM to cells while AZD1480 was added at 0.1, 0.25, 0.5, 1.0 and 2.0 μM. Tumorsphere numbers were counted and compared to DMSO controls. For extended selumetinib treatment and subsequent evaluation of pSTAT3 levels in vitro, Daoy cells were cultured in 1.0 μM selumetinib in stem cell enriched conditions over 3 passages.

Cell sorting and RNA-sequencing

UI226 tumorspheres were grown in stem cell-enriched conditions as previously described²⁴. Tumorspheres were dissociated and 2.5 × 10⁵ cells were orthotopically injected into the cerebellum of NOD SCID mice. Following 14 days of tumor growth, animals were randomly divided into two groups with one receiving selumetinib (75 mg/kg) and the other receiving the 0.5% hydroxypropyl methyl cellulose, 0.1% polysorbate 80 vehicle control via oral gavage on a 5-day on, 2-day off cycle. Upon reaching endpoint, tumors were resected and human cells were isolated by fluorescent activated cell sorting (FACS) using a human leukocyte antigen (HLA) antibody on a MoFlo XDP Cell Sorter (Beckman Coulter). RNA was extracted (N = 3 tumors for control and selumetinib-treated xenografts) using the Norgen RNA extraction kit (Norgen Biotek) according to manufacturer’s instructions. Samples were sent for sequencing and analysis to StemCore labs and the Ottawa Bioinformatics Core Facility (OHRI) as previously described²⁴. RNA-seq was performed comparing sorted control and selumetinib-treated tumor cells from primary xenografts.

RNA-seq reads were mapped independently to the human reference genome (GRCh38 assembly) guided by transcripts from GENCODE release 23, and to the mouse reference genome (GRCm38 assembly) guided by transcripts from GENCODE release M12, using hisat2 2.1.0⁶¹. XenofilteR⁶² was run to deconvolute reads of human origin from those of mouse origin, and mapped human reads were assigned to genes from GENCODE release 23 using featureCounts⁶³. Data were loaded into R and the gene/count matrix was filtered to retain only genes with 5 or more mapped reads in two or more samples. Differential expression was assessed using DESeq2 (v1.20.0)⁶⁴, comparing the two conditions (vehicle control and selumetinib-treated). Expression differences were calculated using the DESeq results function. Lists of significantly differentially expressed genes were identified using a q-value (Benjamini–Hochberg corrected p-value) cut-off of 0.05, which coincides with a p-value of 0.00168.

For analyses of drug-treated tumorspheres, transcripts from tumorsphere samples were quantified with salmon v.1.4.0⁶⁵ against an index built from the GENCODE v35 assembly with inclusion of genomic decoy sequences. Data were loaded into R using the tximport library⁶⁶ and the gene/count matrix was filtered to retain only genes with five or more mapped reads in two or more samples. Differential expression was assessed using DESeq2 v1.30.1; expression differences were calculated using the DESeq2 lfcShrink() function, applying the apeglm method (v 1.12.0)⁶⁷. Multiple testing correction was performed using the Benjamini–Hochberg method, and lists of significantly DE genes were identified using a q-value (i.e. a corrected p value) cut-off of 0.05.

Gene set enrichment analysis

GSEA⁶⁸ was run using pre-ranked gene lists (the GSEA_preranked method) to explore enrichment of pathways and functional classes in gene expression changes following selumetinib or selumetinib/pacritinib treatment of xenografts and tumorspheres. The DE gene lists were filtered to retain protein-coding genes, and ranked by “-log10(p-value) * sign(fold change)”. GSEA was run using the MSigDb (v6.1)⁶⁹ collections H (hallmark gene sets), C4:CP (canonical pathways), C6 (oncogenic signatures), and C8 (cell type signatures). Enrichment maps⁷⁰ were generated in Cytoscape (v3.8.2) from GSEA results, following the protocol described in Reimand et al., 2019⁷¹. Enriched gene sets with a q value <0.005 were retained and clustered using the AutoAnnotate tool; figures were edited to remove unclustered gene sets, and small clusters unrelated to the cell types being studied.

Medulloblastoma patient sample transcriptome and proteomics data

Hierarchical clustering was performed based on the JAK/STAT signaling cascade or the entire JAK/STAT gene set from KEGG and used to determine subgroup affiliation of 763 subtyped primary MB samples previously profiled on the Affymetrix Gene 1.1ST arrays⁵ and normalized using the RMA method (GSE85217). Z-scores were used to generate a signature median across subgroups. The clustering was performed in the R2: Genomics Platform (r2.amc.nl). For proteomic analyses, ProTrack⁷² was used to retrieve and visualize protein data from 218 tumors across several histological classes of childhood brain cancers³⁵. JAK/STAT3 signaling signature genes were chosen from our Hallmark IL6 JAK/STAT3 signaling signature gene set in Supplementary Data 5.

High-Content Tumorsphere Assay

Approximately 4000 UI226 cells were seeded into each well of a 96-well ultra-low attachment plate. All plates contained 3 technical replicates per condition and final results are based on averages from 3 biological replicates. Cells were seeded in various drug concentrations or in the presence of vehicle (DMSO at 0.0025%). For AZD1480 or pacritinib, 15 serially diluted concentrations ranging from 10.0 μM to 0.60 nM were dispensed in triplicate across the plate. Fifteen serially diluted concentrations ranging from 0.002–40 nM were dispensed in triplicate for the chemotherapeutic agent vincristine. Cells were incubated at 37 °C for 5 days. Following incubation, tumorspheres were fixed and stained with a 4% formaldehyde + 300 ng/mL Hoechst 33342 (Thermo Scientific) solution at 4 °C overnight (O/N).

The following day, high-content imaging microscopy was performed using a Cytation 3 imager (BioTek, Winooski, Vermont) equipped with a 16-bit gray scale 1.25 megapixel Sony charge-coupled device camera and a motorized stage. Using a 4× objective, 30 partially overlapping images were acquired per well. Images were digitally aligned to produce a single composite image for analysis.

Automated software (Gen5) was employed to identify and score tumorsphere sizes. Sizes were determined using the ellipsoid equation and measured in μm. Perimeter wells were excluded from analyses to mitigate edge-associated artifacts. To identify and digitally mask tumorspheres, Hoechst detection threshold sensitivity was set to 10,000, and only spheres containing diameters between 45 and 360 μm with circularities greater than 0.3 were scored. Tumorsphere size data were imported into Prism and optimal concentration ranges were identified for further functional studies. To validate selected concentration ranges, Fiji/ImageJ⁷³ was used to evaluate individual tumorsphere size in multiple wells per treatment group. Results were displayed as cumulative frequency distribution of tumorsphere area. Only tumorspheres greater than 25 μm were analyzed.

Immunohistochemistry

Xenograft tumors were sectioned following formalin-fixation and paraffin-embedding. Slides were de-parafinized and re-hydrated in xylene and ethanol gradients followed by antigen retrieval in sodium citrate at 95–100 °C. Slides were washed and then treated for endogenous peroxidases for 10 min. Slides were blocked in 3% lamb serum for 30 min at room temperature (RT) followed by O/N incubation at 4 °C with primary antibody (Supplementary Table 1) prepared in 1% lamb serum/PBS. The following day, slides were first incubated in biotinylated secondary antibodies for 2 h, followed by streptavidin/HRP for 30 min at RT. Slides were then developed with DAB substrate (Sigma), counterstained with hematoxylin, and mounted in Permount (ThermoFisher Scientific).

Immunoblotting

UI226, Daoy, and Ptch+/−:p53+/− mouse SHH MB tumorspheres were grown in culture and treated with pacritinib or AZD1480 as described above. Proteins were loaded onto 10–12% Tris-glycine gels and resolved by SDS-PAGE. Protein was transferred on to nitrocellulose membranes, blocked with 5% milk, and incubated O/N in primary antibodies described in Supplementary Table 1. Secondary antibodies conjugated to HRP were added the following day for 1 h at RT. Signal detection was performed using SuperSignal West Pico. Images were captured with a Fusion FX Vilber Lourmat chemiluminescent imaging system (Marne-la-Vallee cedex 3, France).

High-content tumorsphere drug synergy assay

To evaluate in vitro drug synergy, 4000 UI226 cells/well were seeded into 25 wells (a five-by-five matrix) of a 96-well ultra-low attachment plate. Serially diluted concentrations of selumetinib were dispensed in the y-direction while serially diluted concentrations of AZD1480 or pacritinib were dispensed in the X-direction such that the uppermost left well and bottom most right well contained the most dilute and concentrated drug combinations, respectively. Cells were incubated for 5 days in the presence of drug combinations, fixed, counterstained, imaged, and analyzed in Prism. Each biological replicate employed 3 technical replicates and results were based on a total of three biological replicates. Data were averaged, normalized, and imported into Combenefit⁷⁴ software for synergy calculations. Data were analyzed using the classical Loewe synergy model and values presented in matrices represent the percentage of tumorsphere size reduction that is greater than predicted by an additive model.

Migration assays

Migration assays were performed as previously described²⁴. Briefly, Daoy, UI226, and Ptch+/−:p53+/− MB cells were seeded in round bottom 96-well ultra-low attachment plates at 1 × 10⁴ cells/well, 4 × 10³ cells, and 5 × 10³ cells, respectively. Individual aggregates or spheroids were formed after 4 days, and overlain with a Type I collagen mixture (Collagen type I, DMEM, and NaOH). Once the collagen solidified, DMSO, single agent or combination drug therapies were added. Migration was assessed after 3 days in culture by subtracting the diameter measured at day 0 from the diameter measured at day 3.

Annexin V staining

The Annexin V Apoptosis Detection Kit (BD Biosciences) was used to evaluate cell death as previously described²⁴. Dissociated UI226 SHH MB tumorspheres grown in stem cell conditions and treated with pacritinib + selumetinib or AZD1480+ selumetinib were stained with Annexin V, followed by 7AAD. Samples were run on a MoFlo XDP cell sorter (Beckman Coulter, Mississauga, ON, Canada) and analysis was performed using Kaluza software (Beckman Coulter).

Intracerebellar transplantation and drug treatment

The University of Manitoba Animal Care Committee approved all procedures and protocols. UI226 SHH or HDMB03 Group 3 MB cells were injected into NOD SCID mice as previously described^{20,23,24,40,75,76}. Briefly, UI226 or HDMB03 tumorspheres were dissociated and either 2.5 × 10⁵ cells (UI226) or 1.0 × 10⁵ cells (HDMB03) were injected into the cerebellum of 7–9 week old NOD SCID mice. RCMB18 is a MYCN amplified, TP53 mutated SHH MB PDX model¹⁵. RCMB18 cells were dissociated and 2.5 × 10⁵ cells were injected into the cerebellum of NOD SCID IL2Rg null (NSG) mice. Tumors were propagated and then transplanted into secondary recipients for drug treatment. Following tumor engraftment, animals were randomly separated to either receive vehicle control, selumetinib (37.5 mg/kg), pacritinib (50 mg/kg), or pacritinib + selumetinib (50 mg/kg pacritinib and 37.5 mg/kg selumetinib) twice daily 5 days a week by oral gavage, with a 2 day drug holiday on weekends. For UI226 and RCMB18, treatment was initiated 14 days following tumor cell injection. For HDMB03, treatment was initiated 5 days following tumor cell injection. Animals were treated until they reached endpoint. All drugs were resuspended in sterile 0.5% hydroxypropyl methyl cellulose + 0.1% polysorbate 80. RCMB18 tumor MRI was performed on a MR Solutions cryogen free FlexiScan 7T system (MR Solutions, Guildford, Surrey, UK). Endpoint was determined as a 20% weight loss from peak, lethargy, and ruffled fur. Following endpoint, animals were perfused and brains were extracted for histological analysis as previously described^23,40.

For tumor size calculations in both IHC (N = 4 samples per treatment group) and MRI images (N = 3 mice per treatment group), the ImageJ freehand tool was used to determine area and volume respectively. STEM121 staining was used to demarcate the human tumor cell margins in IHC images. Tumor tissue, consisting of 1 or multiple lobes, in the cerebellum was visually identified on MRI images and used to delineate tumor size using the freehand tool. For each tumor sample, 18 serial sections (300 µm thickness) were obtained by MRI and volume calculated for each slice containing visible tumor tissue. The sum from all the slices containing tumor tissue was used to calculate an overall tumor volume.

Digital spatial profiling on drug-treated xenografted medulloblastoma samples

Four representative slides from FFPE vehicle, pacritinib, selumetinib and pacritinib + selumetinib drug-treated xenografted tumor samples were processed following the GeoMx™ DSP slide preparation user manual (MAN-10087-04). Slides were baked at 60 °C for at least 1 h, and then rehydrated and blocked using NanoString blocking buffer for 1 h. CD271/NGFR-594 (Abcam, ab221212), Ki67-647 (Cell Signaling Technology, 9027S), and MAP2-532 (Novus Bio, NBP1-92711AF532) were added to the sections along with the NanoString protein antibody cocktail that included the available Neuro core panel, a MAPK panel, a PI3K panel and a cell death panel (NanoString Technologies), for O/N incubation at 4 °C. The slides were washed and stained with Syto13 (ThermoFisher, S7575) for 5 min before loading onto the GeoMx™ DSP microscope (NanoString Technologies). Fluorescent images were scanned at 20X and ROIs in the tumor core and along the tumor boarder were selected. Oligos from antibodies were cleaved and collected into 96-well plates, hybridized with NanoString barcodes O/N and then read with an nCounter machine the following day. The digital counts of each antibody for each ROI were generated for data analyses. Per manufacturer’s recommendations, the data were normalized by scaling to the negative control IgG probes, which reflect the rate at which antibodies bind to each region, called signal to noise ratio (SNR). The SNR of each protein target from vehicle, pacritinib, selumetinib, and pacritinib + selumetinib samples was calculated. Significance was determined using paired t-tests (volcano plots) or a one-way ANOVA and a Dunnett’s test for multiple comparisons (box plots). The HDMB03 PQR620-treated samples⁴⁰ underwent the same procedures and workflow; however, CD271 was not included as a morphology marker.

Fluorescence intensity measurements of CD271

CD271/NGFR-594 (Abcam, ab221212) was used as a morphology marker in guiding ROI selections as stated above. TIFF images from all ROIs were downloaded from the GeoMx™ DSP instrument and used to measure the mean fluorescence intensity of CD271 using Image J. CD271 levels in all samples were analyzed by one-way ANOVA and Tukey’s test for multiple comparisons.

Statistics and reproducibility

Prism 8.0 software (GraphPad) was used for all statistical analyses. The Log-rank (Mantel–Cox) test was utilized to evaluate xenograft survival data. ANOVA followed by a Tukey’s test for multiple comparisons were used to assess in vitro assays involving combination drug treatment. Single drug experiments were examined using one-way ANOVA and a Dunnett’s test for multiple comparisons. Brown–Forsythe tests were performed to assess homogeneity of variances. Additional tumorsphere size analysis was performed using two-sample Kolmogorov–Smirnov tests. Significance was considered at P values ≤0.05. All data reported as ± SEM.

Reporting summary

Further information on research design is available in the Nature Research Reporting Summary linked to this article.

Data availability

RNA-seq data have all been deposited into GEO under the accession code GSE165960. All supportive data corresponding to the GEO files are displayed in Figs. 1 and 5, as well as Supplementary Fig. 1 and Supplementary Data 1–7. Additional publicly available datasets utilized include GSE85217 and the Petralia et al.³⁵ proteomics data on ProTrack data portal (http://pbt.cptac-data-view.org/). Original, unprocessed immunoblots associated with Fig. 3c–e and Supplementary Figs. 1b, d and 3a–c can be found in Supplementary Fig. 7a, b online. All source data underlying the graphs and charts presented in the main figures can be found in Supplementary Data 11.

References

Taylor, M. D. et al. Molecular subgroups of medulloblastoma: The current consensus. Acta Neuropathologica 123, 465–472 (2012).
Article CAS PubMed Google Scholar
Mehta, M., Chang, S., Newton, H., Guha, A. & Vogelbaum, M. Principles and Practice of Neuro-oncology: A Multidisciplinary Approach 1st edn, 951 (Demos Medical Publishing, 2011).
Louis, D. N. et al. The 2021 WHO classification of tumors of the central nervous system: A summary. Neuro Oncol. 23, 1231–1251 (2021).
Article CAS PubMed PubMed Central Google Scholar
Northcott, P. A. et al. Medulloblastoma comprises four distinct molecular variants. J. Clin. Oncol. 29, 1408–1414 (2011).
Article PubMed Google Scholar
Cavalli, F. M. G. et al. Intertumoral heterogeneity within medulloblastoma subgroups. Cancer Cell 31, 737–754 e6 (2017).
Article CAS PubMed PubMed Central Google Scholar
Schwalbe, E. C. et al. Novel molecular subgroups for clinical classification and outcome prediction in childhood medulloblastoma: A cohort study. Lancet Oncol. 18, 958–971 (2017).
Article CAS PubMed PubMed Central Google Scholar
Zhukova, N. et al. Subgroup-specific prognostic implications of TP53 mutation in medulloblastoma. J. Clin. Oncol. 31, 2927–2935 (2013).
Article PubMed PubMed Central Google Scholar
Ramaswamy, V. et al. Medulloblastoma subgroup-specific outcomes in irradiated children: Who are the true high-risk patients? Neuro Oncol. 18, 291–297 (2016).
Article PubMed Google Scholar
Ramaswamy, V. et al. Risk stratification of childhood medulloblastoma in the molecular era: The current consensus. Acta Neuropathol. 131, 821–831 (2016).
Article CAS PubMed PubMed Central Google Scholar
Robinson, G. W. et al. Risk-adapted therapy for young children with medulloblastoma (SJYC07): Therapeutic and molecular outcomes from a multicentre, phase 2 trial. Lancet Oncol. 19, 768–784 (2018).
Article CAS PubMed PubMed Central Google Scholar
Garcia-Lopez, J., Kumar, R., Smith, K. S. & Northcott, P. A. Deconstructing Sonic Hedgehog medulloblastoma: Molecular subtypes, drivers, and beyond. Trends Genet. 37, 235–250 (2021).
Article CAS PubMed Google Scholar
Yauch, R. L. et al. Smoothened mutation confers resistance to a Hedgehog pathway inhibitor in medulloblastoma. Science 326, 572–574 (2009).
Article CAS PubMed PubMed Central Google Scholar
Dijkgraaf, G. J. et al. Small molecule inhibition of GDC-0449 refractory smoothened mutants and downstream mechanisms of drug resistance. Cancer Res. 71, 435–444 (2011).
Article CAS PubMed Google Scholar
Robinson, G. W. et al. Medulloblastoma: Results from phase II pediatric brain tumor consortium studies PBTC-025B and PBTC-032. J. Clin. Oncol. 33, 2646–2654 (2015).
Article CAS PubMed PubMed Central Google Scholar
Kool, M. et al. Genome sequencing of SHH medulloblastoma predicts genotype-related response to smoothened inhibition. Cancer Cell 25, 393–405 (2014).
Article CAS PubMed PubMed Central Google Scholar
Kimura, H., Ng, J. M. & Curran, T. Transient inhibition of the Hedgehog pathway in young mice causes permanent defects in bone structure. Cancer Cell 13, 249–260 (2008).
Article CAS PubMed Google Scholar
Read, T. A. et al. Identification of CD15 as a marker for tumor-propagating cells in a mouse model of medulloblastoma. Cancer Cell 15, 135–147 (2009).
Article CAS PubMed PubMed Central Google Scholar
de Antonellis, P. et al. MiR-34a targeting of Notch ligand delta-like 1 impairs CD15+/CD133+ tumor-propagating cells and supports neural differentiation in medulloblastoma. PLoS One 6, e24584 (2011).
Article PubMed CAS PubMed Central Google Scholar
Liu, J. et al. Isolation and characterization of cancer stem cells from medulloblastoma. Genet Mol. Res. 14, 3355–3361 (2015).
Article CAS PubMed Google Scholar
Singh, S. K. et al. Identification of human brain tumour initiating cells. Nature 432, 396–401 (2004).
Article CAS PubMed Google Scholar
Gate, D., Danielpour, M., Bannykh, S. & Town, T. Characterization of cancer stem cells and primary cilia in medulloblastoma. CNS Neurol. Disord. Drug Targets 14, 600–611 (2015).
Article CAS PubMed Google Scholar
Tang, X. et al. Differential proliferative index of cancer stem-like cells in primary and recurrent medulloblastoma in human. Childs Nerv. Syst. 28, 1869–1877 (2012).
Article PubMed Google Scholar
Liang, L. et al. Characterization of novel biomarkers in selecting for subtype specific medulloblastoma phenotypes. Oncotarget 6, 38881–38900 (2015).
Article PubMed PubMed Central Google Scholar
Liang, L. et al. CD271(+) cells are diagnostic and prognostic and exhibit elevated MAPK activity in SHH medulloblastoma. Cancer Res. 78, 4745–4759 (2018).
Article CAS PubMed Google Scholar
Markowitz, D. et al. Pharmacological inhibition of the protein kinase MRK/ZAK radiosensitizes medulloblastoma. Mol. Cancer Ther. 15, 1799–1808 (2016).
Article CAS PubMed Google Scholar
Jensen, K. V., Cseh, O., Aman, A., Weiss, S. & Luchman, H. A. The JAK2/STAT3 inhibitor pacritinib effectively inhibits patient-derived GBM brain tumor initiating cells in vitro and when used in combination with temozolomide increases survival in an orthotopic xenograft model. PLoS One 12, e0189670 (2017).
Article PubMed CAS PubMed Central Google Scholar
McFarland, B. C. et al. Therapeutic potential of AZD1480 for the treatment of human glioblastoma. Mol. Cancer Ther. 10, 2384–2393 (2011).
Article CAS PubMed PubMed Central Google Scholar
Huser, L. et al. SOX2-mediated upregulation of CD24 promotes adaptive resistance towards targeted therapy in melanoma. Int. J. Cancer 143, 3131–3142 (2018).
Article PubMed CAS Google Scholar
Liu, F. et al. Stat3-targeted therapies overcome the acquired resistance to vemurafenib in melanomas. J. Invest. Dermatol. 133, 2041–2049 (2013).
Article CAS PubMed Google Scholar
Notarangelo, T. et al. IL6/STAT3 axis mediates resistance to BRAF inhibitors in thyroid carcinoma cells. Cancer Lett. 433, 147–155 (2018).
Article CAS PubMed Google Scholar
Franklin, D. A., James, J. L., Axelrod, M. L. & Balko, J. M. MEK inhibition activates STAT signaling to increase breast cancer immunogenicity via MHC-I expression. Cancer Drug Resist. 3, 603–612 (2020).
CAS PubMed PubMed Central Google Scholar
Darnell, J. E. Jr. STATs and gene regulation. Science 277, 1630–1635 (1997).
Article CAS PubMed Google Scholar
Murakami, M. et al. IL-6-induced homodimerization of gp130 and associated activation of a tyrosine kinase. Science 260, 1808–1810 (1993).
Article CAS PubMed Google Scholar
Liberzon, A. et al. The Molecular Signatures Database (MSigDB) hallmark gene set collection. Cell Syst. 1, 417–425 (2015).
Article CAS PubMed PubMed Central Google Scholar
Petralia, F. et al. Integrated proteogenomic characterization across major histological types of pediatric brain cancer. Cell 183, 1962–1985 e31 (2020).
Article CAS PubMed PubMed Central Google Scholar
Jacobsen, P. F., Jenkyn, D. J. & Papadimitriou, J. M. Establishment of a human medulloblastoma cell line and its heterotransplantation into nude mice. J. Neuropathol. Exp. Neurol. 44, 472–485 (1985).
Article CAS PubMed Google Scholar
Ward, R. J. et al. Multipotent CD15+ cancer stem cells in patched-1-deficient mouse medulloblastoma. Cancer Res. 69, 4682–4690 (2009).
Article CAS PubMed Google Scholar
Zhao, X. et al. RAS/MAPK activation drives resistance to Smo inhibition, metastasis, and tumor evolution in Shh pathway-dependent tumors. Cancer Res. 75, 3623–3635 (2015).
Article CAS PubMed PubMed Central Google Scholar
Merritt, C. R. et al. Multiplex digital spatial profiling of proteins and RNA in fixed tissue. Nat. Biotechnol. 38, 586–599 (2020).
Article CAS PubMed Google Scholar
Zagozewski, J. et al. An OTX2-PAX3 signaling axis regulates Group 3 medulloblastoma cell fate. Nat. Commun. 11, 3627 (2020).
Article CAS PubMed PubMed Central Google Scholar
Zhang, L. et al. Single-cell transcriptomics in medulloblastoma reveals tumor-initiating progenitors and oncogenic cascades during tumorigenesis and relapse. Cancer Cell 36, 302–318.e7 (2019).
Article CAS PubMed PubMed Central Google Scholar
Jensen, K. V., Hao, X., Aman, A., Luchman, H. A. & Weiss, S. EGFR blockade in GBM brain tumor stem cells synergizes with JAK2/STAT3 pathway inhibition to abrogate compensatory mechanisms in vitro and in vivo. Neurooncol Adv. 2, vdaa020 (2020).
PubMed PubMed Central Google Scholar
Bid, H. K. et al. Development, characterization, and reversal of acquired resistance to the MEK1 inhibitor selumetinib (AZD6244) in an in vivo model of childhood astrocytoma. Clin. Cancer Res. 19, 6716–6729 (2013).
Article CAS PubMed PubMed Central Google Scholar
Jin, J., Guo, Q., Xie, J., Jin, D. & Zhu, Y. Combination of MEK inhibitor and the JAK2-STAT3 pathway inhibition for the therapy of colon cancer. Pathol. Oncol. Res. 25, 769–775 (2019).
Article CAS PubMed Google Scholar
Nagathihalli, N. S. et al. Inverse correlation of STAT3 and MEK signaling mediates resistance to RAS pathway inhibition in pancreatic cancer. Cancer Res. 78, 6235–6246 (2018).
Article CAS PubMed PubMed Central Google Scholar
Dai, B. et al. STAT3 mediates resistance to MEK inhibitor through microRNA miR-17. Cancer Res. 71, 3658–3668 (2011).
Article CAS PubMed PubMed Central Google Scholar
Banerjee, A. et al. A phase I trial of the MEK inhibitor selumetinib (AZD6244) in pediatric patients with recurrent or refractory low-grade glioma: A Pediatric Brain Tumor Consortium (PBTC) study. Neuro Oncol. 19, 1135–1144 (2017).
Article CAS PubMed PubMed Central Google Scholar
Fangusaro, J. et al. Selumetinib in paediatric patients with BRAF-aberrant or neurofibromatosis type 1-associated recurrent, refractory, or progressive low-grade glioma: a multicentre, phase 2 trial. Lancet Oncol. 20, 1011–1022 (2019).
Article CAS PubMed PubMed Central Google Scholar
Dombi, E. et al. Activity of selumetinib in neurofibromatosis type 1-related plexiform neurofibromas. N. Engl. J. Med. 375, 2550–2560 (2016).
Article CAS PubMed PubMed Central Google Scholar
Gross, A., Bishop, R. & Widemann, B. C. Selumetinib in plexiform neurofibromas. N. Engl. J. Med. 376, 1195 (2017).
Article PubMed Google Scholar
Perreault, S. et al. A phase 2 study of trametinib for patients with pediatric glioma or plexiform neurofibroma with refractory tumor and activation of the MAPK/ERK pathway: TRAM-01. BMC Cancer 19, 1250 (2019).
Article CAS PubMed PubMed Central Google Scholar
Bouffet, E. et al. Trametinib therapy in pediatric patients with low-grade gliomas (LGG) with BRAF gene fusion; a disease specific cohort in the first pediatric testing of trametinib. Neuro-Oncol. 20, i114 (2018).
Article PubMed Central Google Scholar
Selt, F. et al. Response to trametinib treatment in progressive pediatric low-grade glioma patients. J. Neurooncol. 149, 499–510 (2020).
Sreenivasan, L. et al. Autocrine IL-6/STAT3 signaling aids development of acquired drug resistance in Group 3 medulloblastoma. Cell Death Dis. 11, 1035 (2020).
Article CAS PubMed PubMed Central Google Scholar
Smith, M. P. et al. The immune microenvironment confers resistance to MAPK pathway inhibitors through macrophage-derived TNFalpha. Cancer Discov. 4, 1214–1229 (2014).
Article CAS PubMed PubMed Central Google Scholar
Xie, M. et al. MEK inhibition modulates cytokine response to mediate therapeutic efficacy in lung cancer. Cancer Res. 79, 5812–5825 (2019).
Article CAS PubMed PubMed Central Google Scholar
Spiller, S. E., Logsdon, N. J., Deckard, L. A. & Sontheimer, H. Inhibition of nuclear factor kappa-B signaling reduces growth in medulloblastoma in vivo. BMC Cancer 11, 136 (2011).
Article CAS PubMed PubMed Central Google Scholar
Northcott, P. A. et al. Rapid, reliable, and reproducible molecular sub-grouping of clinical medulloblastoma samples. Acta Neuropathol. 123, 615–626 (2012).
Article CAS PubMed Google Scholar
Milde, T. et al. HD-MB03 is a novel Group 3 medulloblastoma model demonstrating sensitivity to histone deacetylase inhibitor treatment. J. Neurooncol. 110, 335–348 (2012).
Article CAS PubMed Google Scholar
Wetmore, C., Eberhart, D. E. & Curran, T. Loss of p53 but not ARF accelerates medulloblastoma in mice heterozygous for patched. Cancer Res. 61, 513–516 (2001).
CAS PubMed Google Scholar
Pertea, M., Kim, D., Pertea, G. M., Leek, J. T. & Salzberg, S. L. Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown. Nat. Protoc. 11, 1650–1667 (2016).
Article CAS PubMed PubMed Central Google Scholar
Kluin, R. J. C. et al. XenofilteR: Computational deconvolution of mouse and human reads in tumor xenograft sequence data. BMC Bioinforma. 19, 366 (2018).
Article CAS Google Scholar
Liao, Y., Smyth, G. K. & Shi, W. featureCounts: An efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 30, 923–930 (2014).
Article CAS PubMed Google Scholar
Love, M. I., Huber, W. & Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 550 (2014).
Article PubMed CAS PubMed Central Google Scholar
Patro, R., Duggal, G., Love, M. I., Irizarry, R. A. & Kingsford, C. Salmon provides fast and bias-aware quantification of transcript expression. Nat. Methods 14, 417–419 (2017).
Article CAS PubMed PubMed Central Google Scholar
Soneson, C., Love, M. I. & Robinson, M. D. Differential analyses for RNA-seq: Transcript-level estimates improve gene-level inferences. F1000Res 4, 1521 (2015).
Article PubMed CAS Google Scholar
Zhu, A., Ibrahim, J. G. & Love, M. I. Heavy-tailed prior distributions for sequence count data: removing the noise and preserving large differences. Bioinformatics 35, 2084–2092 (2019).
Article CAS PubMed Google Scholar
Subramanian, A. et al. Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl Acad. Sci. USA 102, 15545–15550 (2005).
Article CAS PubMed PubMed Central Google Scholar
Liberzon, A. et al. Molecular signatures database (MSigDB) 3.0. Bioinformatics 27, 1739–1740 (2011).
Article CAS PubMed PubMed Central Google Scholar
Merico, D., Isserlin, R., Stueker, O., Emili, A. & Bader, G. D. Enrichment map: A network-based method for gene-set enrichment visualization and interpretation. PLoS One 5, e13984 (2010).
Article PubMed CAS PubMed Central Google Scholar
Reimand, J. et al. Pathway enrichment analysis and visualization of omics data using g:Profiler, GSEA, Cytoscape and EnrichmentMap. Nat. Protoc. 14, 482–517 (2019).
Article CAS PubMed PubMed Central Google Scholar
Calinawan, A. P. et al. ProTrack: An interactive multi-omics data browser for proteogenomic studies. Proteomics 20, e1900359 (2020).
Article PubMed CAS PubMed Central Google Scholar
Schindelin, J. et al. Fiji: An open-source platform for biological-image analysis. Nat. Methods 9, 676–682 (2012).
Article CAS PubMed Google Scholar
Di Veroli, G. Y. et al. Combenefit: An interactive platform for the analysis and visualization of drug combinations. Bioinformatics 32, 2866–2868 (2016).
Article PubMed CAS PubMed Central Google Scholar
Morrison, L. C. et al. Deconstruction of medulloblastoma cellular heterogeneity reveals differences between the most highly invasive and self-renewing phenotypes. Neoplasia 15, 384–398 (2013).
Article PubMed CAS PubMed Central Google Scholar
Werbowetski-Ogilvie, T. E., Morrison, L. C., Fiebig-Comyn, A. & Bhatia, M. In vivo generation of neural tumors from neoplastic pluripotent stem cells models early human pediatric brain tumor formation. Stem Cells 30, 392–404 (2012).
Article CAS PubMed Google Scholar

Download references

Acknowledgements

We thank Agnes Fresnoza and Shawn Blum from Vet Services at the University of Manitoba for outstanding technical support during the course of our animal studies. We thank Dr. Mike Jackson from the Small Animal and Materials Imaging Core Facility, University of Manitoba, for performing MRI. This work was funded by the Canada Research Chairs Tier II Program, and operating grants from the Canadian Institutes of Health Research, The Rally Foundation for Childhood Cancer Research, and the CancerCare Manitoba Foundation (all to TEWO). J.Z. was supported by a joint Postdoctoral Fellowship from the CancerCare Manitoba Foundation, the Children’s Hospital Research Institute of Manitoba and Research Manitoba. S.B. was supported by the Research Manitoba-University of Manitoba WM Ross Scholarship-CancerCare Manitoba Foundation Master’s Studentship Award. B.J.G. was supported by a William Donald Nash Brain Tumor Foundation of Canada Fellowship.

Author information

These authors contributed equally: Jamie Zagozewski, Stephanie Borlase.

Authors and Affiliations

Department of Biochemistry and Medical Genetics, Rady Faculty of Health Sciences, University of Manitoba, Winnipeg, MB, Canada
Jamie Zagozewski, Stephanie Borlase, Brent J. Guppy, Ludivine Coudière-Morrison, Ghazaleh M. Shahriary, Victor Gordon, Lisa Liang, Stephen Cheng & Tamra E. Werbowetski-Ogilvie
Ottawa Bioinformatics Core Facility, Ottawa Hospital Research Institute, Ottawa, ON, Canada
Christopher J. Porter
Central Animal Care Services, University of Manitoba, Winnipeg, MB, Canada
Rhonda Kelley
The Arthur and Sonia Labatt Brain Tumour Research Centre, The Hospital for Sick Children, Toronto, ON, Canada
Cynthia Hawkins & Vijay Ramaswamy
Department of Laboratory Medicine and Pathobiology, University of Toronto, Toronto, ON, Canada
Cynthia Hawkins
Department of Pathology & Laboratory Medicine, University of Calgary, Calgary, AB, Canada
Jennifer A. Chan
Arnie Charbonneau Cancer Institute, University of Calgary, Calgary, AB, Canada
Jennifer A. Chan
Pathology Department, NanoString Inc, Seattle, WA, USA
Yan Liang & Jingjing Gong
Developmental and Stem Cell Biology Program, The Hospital for Sick Children, Toronto, ON, Canada
Carolina Nör, Olivier Saulnier & Vijay Ramaswamy
Sanford Burnham Prebys Medical Discovery Institute, La Jolla, CA, USA
Robert J. Wechsler-Reya
Division of Haematology/Oncology, The Hospital for Sick Children, Toronto, ON, Canada
Vijay Ramaswamy
Department of Medical Biophysics, University of Toronto, Toronto, ON, Canada
Vijay Ramaswamy
CancerCare Manitoba Research Institute, Winnipeg, MB, Canada
Tamra E. Werbowetski-Ogilvie

Authors

Jamie Zagozewski
View author publications
You can also search for this author in PubMed Google Scholar
Stephanie Borlase
View author publications
You can also search for this author in PubMed Google Scholar
Brent J. Guppy
View author publications
You can also search for this author in PubMed Google Scholar
Ludivine Coudière-Morrison
View author publications
You can also search for this author in PubMed Google Scholar
Ghazaleh M. Shahriary
View author publications
You can also search for this author in PubMed Google Scholar
Victor Gordon
View author publications
You can also search for this author in PubMed Google Scholar
Lisa Liang
View author publications
You can also search for this author in PubMed Google Scholar
Stephen Cheng
View author publications
You can also search for this author in PubMed Google Scholar
Christopher J. Porter
View author publications
You can also search for this author in PubMed Google Scholar
Rhonda Kelley
View author publications
You can also search for this author in PubMed Google Scholar
Cynthia Hawkins
View author publications
You can also search for this author in PubMed Google Scholar
Jennifer A. Chan
View author publications
You can also search for this author in PubMed Google Scholar
Yan Liang
View author publications
You can also search for this author in PubMed Google Scholar
Jingjing Gong
View author publications
You can also search for this author in PubMed Google Scholar
Carolina Nör
View author publications
You can also search for this author in PubMed Google Scholar
Olivier Saulnier
View author publications
You can also search for this author in PubMed Google Scholar
Robert J. Wechsler-Reya
View author publications
You can also search for this author in PubMed Google Scholar
Vijay Ramaswamy
View author publications
You can also search for this author in PubMed Google Scholar
Tamra E. Werbowetski-Ogilvie
View author publications
You can also search for this author in PubMed Google Scholar

Contributions

J.Z., S.B., and T.E.W.O. were primarily responsible for conception and design of experiments, data collection, assembly and interpretation of data, manuscript writing, and final approval. Acquisition of data, data assembly, and interpretation: B.J.G., G.M.S., L.C.M., V.G., L.L., S.C., R.K., C.H., J.A.C., C.J.P., Y.L., J.G., C.N., O.S., R.W.R., and V.R. Additional manuscript writing and review: J.Z., S.B., L.C.M., and T.E.W.O.

Corresponding author

Correspondence to Tamra E. Werbowetski-Ogilvie.

Ethics declarations

Competing interests

The authors declare no competing interests.

Peer review

Peer review information

Communications Biology thanks Qing Liu and the other anonymous reviewer(s) for their contribution to the peer review of this work. Primary Handling Editors: Eve Rogers. Peer reviewer reports are available.

Additional information

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Peer Review File

Supplementary Information

Description of Additional Supplementary Files

Supplementary Data 1-10

Supplementary Data 11

Reporting Summary

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Cite this article

Zagozewski, J., Borlase, S., Guppy, B.J. et al. Combined MEK and JAK/STAT3 pathway inhibition effectively decreases SHH medulloblastoma tumor progression. Commun Biol 5, 697 (2022). https://doi.org/10.1038/s42003-022-03654-9

Download citation

Received: 15 October 2021
Accepted: 30 June 2022
Published: 14 July 2022
DOI: https://doi.org/10.1038/s42003-022-03654-9

This article is cited by

Mapping cancer biology in space: applications and perspectives on spatial omics for oncology
- Sumin Lee
- Gyeongjun Kim
- Sunghoon Kwon
Molecular Cancer (2024)
The novel oncogenic factor TET3 combines with AHR to promote thyroid cancer lymphangiogenesis via the HIF-1α/VEGF signaling pathway
- Liyun Yang
- Runyu Zhao
- Shuixian Huang
Cancer Cell International (2023)
Cross-species analysis of SHH medulloblastoma models reveals significant inhibitory effects of trametinib on tumor progression
- Stephanie Borlase
- Alexandria DeCarlo
- Tamra E. Werbowetski-Ogilvie
Cell Death Discovery (2023)

Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.