A prognostic matrix code defines functional glioblastoma phenotypes and niches.

Interactions among tumor, immune and vascular niches play major roles in driving glioblastoma (GBM) malignancy and treatment responses. The composition, heterogeneity, and localization of extracellular core matrix proteins (CMPs) that mediate such interactions, however, are not well understood. Here, we characterize functional and clinical relevance of genes encoding CMPs in GBM at bulk, single cell, and spatial anatomical resolution. We identify a "matrix code" for genes encoding CMPs whose expression levels categorize GBM tumors into matrisome-high and matrisome-low groups that correlate with worse and better patient survival, respectively. The matrisome enrichment is associated with specific driver oncogenic alterations, mesenchymal state, infiltration of pro-tumor immune cells and immune checkpoint gene expression. Anatomical and single cell transcriptome analyses indicate that matrisome gene expression is enriched in vascular and leading edge/infiltrative anatomic structures that are known to harbor glioma stem cells driving GBM progression. Finally, we identified a 17-gene matrisome signature that retains and further refines the prognostic value of genes encoding CMPs and, importantly, potentially predicts responses to PD1 blockade in clinical trials for GBM. The matrisome gene expression profiles provide potential biomarkers of functionally relevant GBM niches that contribute to mesenchymal-immune cross talk and patient stratification which could be applied to optimize treatment responses.

A prognostic matrix code defines functional glioblastoma phenotypes and niches.

Interactions among tumor, immune and vascular niches play major roles in driving glioblastoma (GBM) malignancy and treatment responses. The composition, heterogeneity, and localization of extracellular core matrix proteins (CMPs) that mediate such interactions, however, are not well understood. Here, we characterize functional and clinical relevance of genes encoding CMPs in GBM at bulk, single cell, and spatial anatomical resolution. We identify a "matrix code" for genes encoding CMPs whose expression levels categorize GBM tumors into matrisome-high and matrisome-low groups that correlate with worse and better survival, respectively, of patients. The matrisome enrichment is associated with specific driver oncogenic alterations, mesenchymal state, infiltration of pro-tumor immune cells and immune checkpoint gene expression. Anatomical and single cell transcriptome analyses indicate that matrisome gene expression is enriched in vascular and leading edge/infiltrative anatomic structures that are known to harbor glioma stem cells driving GBM progression. Finally, we identified a 17-gene matrisome signature that retains and further refines the prognostic value of genes encoding CMPs and, importantly, potentially predicts responses to PD1 blockade in clinical trials for GBM. The matrisome gene expression profiles may provide biomarkers of functionally relevant GBM niches that contribute to mesenchymal-immune cross talk and patient stratification to optimize treatment responses.

The antidepressant imipramine inhibits breast cancer growth by targeting estrogen receptor signaling and DNA repair events.

Aberrant activities of various cell cycle and DNA repair proteins promote cancer growth and progression and render them resistant to therapies. Here, we demonstrate that the anti-depressant imipramine blocks growth of triple-negative (TNBC) and estrogen receptor-positive (ER+) breast cancers by inducing cell cycle arrest and by blocking heightened homologous recombination (HR) and non-homologous end joining-mediated (NHEJ) DNA repair activities. Our results reveal that imipramine inhibits the expression of several cell cycle- and DNA repair-associated proteins including E2F1, CDK1, Cyclin D1, and RAD51. In addition, we show that imipramine inhibits the growth of ER + breast cancers by inhibiting the estrogen receptor- α (ER-α) signaling. Our studies in preclinical mouse models and ex vivo explants from breast cancer patients show that imipramine sensitizes TNBC to the PARP inhibitor olaparib and endocrine resistant ER + breast cancer to anti-estrogens. Our studies suggest that repurposing imipramine could enhance routine care for breast cancer patients. Based on these results, we designed an ongoing clinical trial, where we are testing the efficacy of imipramine for treating patients with triple-negative and estrogen receptor-positive breast cancer. Since aberrant DNA repair activity is used by many cancers to survive and become resistant to therapy, imipramine could be used alone and/or with currently used drugs for treating many aggressive cancers.

Targeting aberrant replication and DNA repair events for treating breast cancers.

The major limitations of DNA-targeting chemotherapy drugs include life-threatening toxicity, acquired resistance and occurrence of secondary cancers. Here, we report a small molecule, Carbazole Blue (CB), that binds to DNA and inhibits cancer growth and metastasis by targeting DNA-related processes that tumor cells use but not the normal cells. We show that CB inhibits the expression of pro-tumorigenic genes that promote unchecked replication and aberrant DNA repair that cancer cells get addicted to survive. In contrast to chemotherapy drugs, systemic delivery of CB suppressed breast cancer growth and metastasis with no toxicity in pre-clinical mouse models. Using PDX and ex vivo explants from estrogen receptor (ER) positive, ER mutant and TNBC patients, we further demonstrated that CB effectively blocks therapy-sensitive and therapy-resistant breast cancer growth without affecting normal breast tissue. Our data provide a strong rationale to develop CB as a viable therapeutic for treating breast cancers.

M6A RNA Methylation Regulates Histone Ubiquitination to Support Cancer Growth and Progression.

Osteosarcoma is the most common malignancy of the bone, yet the survival for patients with osteosarcoma is virtually unchanged over the past 30 years. This is principally because development of new therapies is hampered by a lack of recurrent mutations that can be targeted in osteosarcoma. Here, we report that epigenetic changes via mRNA methylation holds great promise to better understand the mechanisms of osteosarcoma growth and to develop targeted therapeutics. In patients with osteosarcoma, the RNA demethylase ALKBH5 was amplified and higher expression correlated with copy-number changes. ALKBH5 was critical for promoting osteosarcoma growth and metastasis, yet it was dispensable for normal cell survival. Methyl RNA immunoprecipitation sequencing analysis and functional studies showed that ALKBH5 mediates its protumorigenic function by regulating m6A levels of histone deubiquitinase USP22 and the ubiquitin ligase RNF40. ALKBH5-mediated m6A deficiency in osteosarcoma led to increased expression of USP22 and RNF40 that resulted in inhibition of histone H2A monoubiquitination and induction of key protumorigenic genes, consequently driving unchecked cell-cycle progression, incessant replication, and DNA repair. RNF40, which is historically known to ubiquitinate H2B, inhibited H2A ubiquitination in cancer by interacting with and affecting the stability of DDB1-CUL4-based ubiquitin E3 ligase complex. Taken together, this study directly links increased activity of ALKBH5 with dysregulation of USP22/RNF40 and histone ubiquitination in cancers. More broadly, these results suggest that m6A RNA methylation works in concert with other epigenetic mechanisms to control cancer growth. SIGNIFICANCE: RNA demethylase ALKBH5 upregulates USP22 and RNF40 to inhibit histone H2A ubiquitination and induces expression of key replication and DNA repair-associated genes, driving osteosarcoma progression.

Single-cell analysis of human glioma and immune cells identifies S100A4 as an immunotherapy target.

A major rate-limiting step in developing more effective immunotherapies for GBM is our inadequate understanding of the cellular complexity and the molecular heterogeneity of immune infiltrates in gliomas. Here, we report an integrated analysis of 201,986 human glioma, immune, and other stromal cells at the single cell level. In doing so, we discover extensive spatial and molecular heterogeneity in immune infiltrates. We identify molecular signatures for nine distinct myeloid cell subtypes, of which five are independent prognostic indicators of glioma patient survival. Furthermore, we identify S100A4 as a regulator of immune suppressive T and myeloid cells in GBM and demonstrate that deleting S100a4 in non-cancer cells is sufficient to reprogram the immune landscape and significantly improve survival. This study provides insights into spatial, molecular, and functional heterogeneity of glioma and glioma-associated immune cells and demonstrates the utility of this dataset for discovering therapeutic targets for this poorly immunogenic cancer.

Cancer stem cells, not bulk tumor cells, determine mechanisms of resistance to SMO inhibitors.

The emergence of treatment resistance significantly reduces the clinical utility of many effective targeted therapies. Although both genetic and epigenetic mechanisms of drug resistance have been reported, whether these mechanisms are stochastically selected in individual tumors or governed by a predictable underlying principle is unknown. Here, we report that the dependence of cancer stem cells (CSCs), not bulk tumor cells, on the targeted pathway determines the molecular mechanism of resistance in individual tumors. Using both spontaneous and transplantable mouse models of sonic hedgehog (SHH) medulloblastoma (MB) treated with an SHH/Smoothened inhibitor, sonidegib/LDE225, we show that genetic-based resistance occurs only in tumors that contain SHH-dependent CSCs (SD-CSCs). In contrast, SHH MBs containing SHH-dependent bulk tumor cells but SHH-independent CSCs (SI-CSCs) acquire resistance through epigenetic reprogramming. Mechanistically, elevated proteasome activity in SMOi-resistant SI-CSC MBs alters the tumor cell maturation trajectory through enhanced degradation of specific epigenetic regulators, including histone acetylation machinery components, resulting in global reductions in H3K9Ac, H3K14Ac, H3K56Ac, H4K5Ac, and H4K8Ac marks and gene expression changes. These results provide new insights into how selective pressure on distinct tumor cell populations contributes to different mechanisms of resistance to targeted therapies. This insight provides a new conceptual framework to understand responses and resistance to SMOis and other targeted therapies.

MiR-584-5p potentiates vincristine and radiation response by inducing spindle defects and DNA damage in medulloblastoma.

Despite improvements in overall survival, only a modest percentage of patients survives high-risk medulloblastoma. The devastating side effects of radiation and chemotherapy substantially reduce quality of life for surviving patients. Here, using genomic screens, we identified miR-584-5p as a potent therapeutic adjuvant that potentiates medulloblastoma to radiation and vincristine. MiR-584-5p inhibited medulloblastoma growth and prolonged survival of mice in pre-clinical tumor models. MiR-584-5p overexpression caused cell cycle arrest, DNA damage, and spindle defects in medulloblastoma cells. MiR-584-5p mediated its tumor suppressor and therapy-sensitizing effects by targeting HDAC1 and eIF4E3. MiR-584-5p overexpression or HDAC1/eIF4E3 silencing inhibited medulloblastoma stem cell self-renewal without affecting neural stem cell growth. In medulloblastoma patients, reduced expression of miR-584-5p correlated with increased levels of HDAC1/eIF4E3. These findings identify a previously undefined role for miR-584-5p/HDAC1/eIF4E3 in regulating DNA repair, microtubule dynamics, and stemness in medulloblastoma and set the stage for a new way to treat medulloblastoma using miR-584-5p.

Cross-talk among writers, readers, and erasers of m6A regulates cancer growth and progression.

The importance of RNA methylation in biological processes is an emerging focus of investigation. We report that altering m6A levels by silencing either N 6-adenosine methyltransferase METTL14 (methyltransferase-like 14) or demethylase ALKBH5 (ALKB homolog 5) inhibits cancer growth and invasion. METTL14/ALKBH5 mediate their protumorigenic function by regulating m6A levels of key epithelial-mesenchymal transition and angiogenesis-associated transcripts, including transforming growth factor-ß signaling pathway genes. Using MeRIP-seq (methylated RNA immunoprecipitation sequencing) analysis and functional studies, we find that these target genes are particularly sensitive to changes in m6A modifications, as altered m6A status leads to aberrant expression of these genes, resulting in inappropriate cell cycle progression and evasion of apoptosis. Our results reveal that METTL14 and ALKBH5 determine the m6A status of target genes by controlling each other's expression and by inhibiting m6A reader YTHDF3 (YTH N 6-methyladenosine RNA binding protein 3), which blocks RNA demethylase activity. Furthermore, we show that ALKBH5/METTL14 constitute a positive feedback loop with RNA stability factor HuR to regulate the stability of target transcripts. We discover that hypoxia alters the level/activity of writers, erasers, and readers, leading to decreased m6A and consequently increased expression of target transcripts in cancer cells. This study unveils a previously undefined role for m6A in cancer and shows that the collaboration among writers-erasers-readers sets up the m6A threshold to ensure the stability of progrowth/proliferation-specific genes, and protumorigenic stimulus, such as hypoxia, perturbs that m6A threshold, leading to uncontrolled expression/activity of those genes, resulting in tumor growth, angiogenesis, and progression.

Novel post-transcriptional and post-translational regulation of pro-apoptotic protein BOK and anti-apoptotic protein Mcl-1 determine the fate of breast cancer cells to survive or die.

Deregulation of apoptosis is central to cancer progression and a major obstacle to effective treatment. The Bcl-2 gene family members play important roles in the regulation of apoptosis and are frequently altered in cancers. One such member is pro-apoptotic protein Bcl-2-related Ovarian Killer (BOK). Despite its critical role in apoptosis, the regulation of BOK expression is poorly understood in cancers. Here, we discovered that miR-296-5p regulates BOK expression by binding to its 3'-UTR in breast cancers. Interestingly, miR-296-5p also regulates the expression of anti-apoptotic protein myeloid cell leukemia 1 (Mcl-1), which is highly expressed in breast cancers. Our results reveal that Mcl-1 and BOK constitute a regulatory feedback loop as ectopic BOK expression induces Mcl-1, whereas silencing of Mcl-1 results in reduced BOK levels in breast cancer cells. In addition, we show that silencing of Mcl-1 but not BOK reduced the long-term growth of breast cancer cells. Silencing of both Mcl-1 and BOK rescued the effect of Mcl-1 silencing on breast cancer cell growth, suggesting that BOK is important for attenuating cell growth in the absence of Mcl-1. Depletion of BOK suppressed caspase-3 activation in the presence of paclitaxel and in turn protected cells from paclitaxel-induced apoptosis. Furthermore, we demonstrate that glycogen synthase kinase (GSK3) α/ß interacts with BOK and regulates its level post-translationally in breast cancer cells. Taken together, our results suggest that fine tuning of the levels of pro-apoptotic protein BOK and anti-apoptotic protein Mcl-1 may decide the fate of cancer cells to either undergo apoptosis or proliferation.