Programs, Origins, and Niches of Immunomodulatory Myeloid Cells in Gliomas

Gliomas are incurable malignancies notable for an immunosuppressive microenvironment with abundant myeloid cells whose immunomodulatory properties remain poorly defined. Here, utilizing scRNA-seq data for 183,062 myeloid cells from 85 human tumors, we discover that nearly all glioma-associated myeloid cells express at least one of four immunomodulatory activity programs: Scavenger Immunosuppressive, C1Q Immunosuppressive, CXCR4 Inflammatory, and IL1B Inflammatory. All four programs are present in IDH1 mutant and wild-type gliomas and are expressed in macrophages, monocytes, and microglia whether of blood or resident myeloid cell origins. Integrating our scRNA-seq data with mitochondrial DNA-based lineage tracing, spatial transcriptomics, and organoid explant systems that model peripheral monocyte infiltration, we show that these programs are driven by microenvironmental cues and therapies rather than myeloid cell type, origin, or mutation status. The C1Q Immunosuppressive program is driven by routinely administered dexamethasone. The Scavenger Immunosuppressive program includes ligands with established roles in T-cell suppression, is induced in hypoxic regions, and is associated with immunotherapy resistance. Both immunosuppressive programs are less prevalent in lower-grade gliomas, which are instead enriched for the CXCR4 Inflammatory program. Our study provides a framework to understand immunomodulatory myeloid cells in glioma, and a foundation to develop more effective immunotherapies.

DrugMap: A quantitative pan-cancer analysis of cysteine ligandability

Cysteine-focused chemical proteomic platforms have accelerated the clinical development of covalent inhibitors of a wide-range of targets in cancer. However, how different oncogenic contexts influence cysteine targeting remains unknown. To address this question, we have developed DrugMap , an atlas of cysteine ligandability compiled across 416 cancer cell lines. We unexpectedly find that cysteine ligandability varies across cancer cell lines, and we attribute this to differences in cellular redox states, protein conformational changes, and genetic mutations. Leveraging these findings, we identify actionable cysteines in NFκB1 and SOX10 and develop corresponding covalent ligands that block the activity of these transcription factors. We demonstrate that the NFκB1 probe blocks DNA binding, whereas the SOX10 ligand increases SOX10-SOX10 interactions and disrupts melanoma transcriptional signaling. Our findings reveal heterogeneity in cysteine ligandability across cancers, pinpoint cell-intrinsic features driving cysteine targeting, and illustrate the use of covalent probes to disrupt oncogenic transcription factor activity.

Elucidating the diversity of malignant mesenchymal states in glioblastoma by integrative analysis

Background

Multiple glioblastoma studies have described a mesenchymal (MES) state, with each study defining the MES program by distinct sets of genes and highlighting distinct functional associations, including both immune activation and suppression. These variable descriptions complicate our understanding of the MES state and its implications. Here, we hypothesize that there is a range of glioma MES states, possibly reflecting distinct prior states in which a MES program can be induced, and/or distinct mechanisms that induce the MES states in those cells.

Methods

We integrated multiple published single-cell and bulk RNA sequencing datasets and MES signatures to define a core MES program that recurs across studies, as well as multiple function-specific MES signatures that vary across MES cells. We then examined the co-occurrence of these signatures and their associations with genetic and microenvironmental features.

Results

Based on co-occurrence of MES signatures, we found three main variants of MES states: hypoxia-related (MES-Hyp), astrocyte-related (MES-Ast), and an intermediate state. Notably, the MES states are differentially associated with genetic and microenvironmental features. MES-Hyp is preferentially associated with NF1 deletion, overall macrophage abundance, a high macrophage/microglia ratio, and M2-related macrophages, consistent with previous studies that associated MES with immune suppression. In contrast, MES-Ast is associated with T cell abundance and cytotoxicity, consistent with immune activation through expression of MHC-I/II.

Conclusions

Diverse MES states occur in glioblastoma. These states share a subset of core genes but differ primarily in their association with hypoxia vs. astrocytic expression programs, and with immune suppression vs. activation, respectively.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13073-022-01109-8.

GLI3 Is Associated With Neuronal Differentiation in SHH-Activated and WNT-Activated Medulloblastoma

Glioma-associated oncogene homolog 3 (GLI3), whose main function is to inhibit GLI1, has been associated with neuronal differentiation in medulloblastoma. However, it is not clear what molecular subtype(s) show increased GLI3 expression. GLI3 levels were assessed by immunohistochemistry in 2 independent cohorts, including a total of 88 cases, and found to be high in both WNT- and SHH-activated medulloblastoma. Analysis of bulk mRNA expression data and single cell RNA sequencing studies confirmed that GLI1 and GLI3 are highly expressed in SHH-activated medulloblastoma, whereas GLI3 but not GLI1 is highly expressed in WNT-activated medulloblastoma. Immunohistochemical analysis has shown that GLI3 is expressed inside the neuronal differentiated nodules of SHH-activated medulloblastoma, whereas GLI1/2 are expressed in desmoplastic areas. In contrast, GLI3 is diffusely expressed in WNT-activated medulloblastoma, whereas GLI1 is suppressed. Our data suggest that GLI3 may be a master regulator of neuronal differentiation and morphology in these subgroups.

Re-evaluating Biopsy for Recurrent Glioblastoma: A Position Statement by the Christopher Davidson Forum Investigators

Patients with glioblastoma (GBM) need bold new approaches to their treatment, yet progress has been hindered by a relative inability to dynamically track treatment response, mechanisms of resistance, evolution of targetable mutations, and changes in mutational burden. We are writing on behalf of a multidisciplinary group of academic neuro-oncology professionals who met at the collaborative Christopher Davidson Forum at Washington University in St Louis in the fall of 2019. We propose a dramatic but necessary change to the routine management of patients with GBM to advance the field: to routinely biopsy recurrent GBM at the time of presumed recurrence. Data derived from these samples will identify true recurrence vs treatment effect, avoid treatments with little chance of success, enable clinical trial access, and aid in the scientific advancement of our understanding of GBM.

EPEN-21. IMPAIRED NEURONAL-GLIAL FATE SPECIFICATION IN PEDIATRIC EPENDYMOMA REVEALED BY SINGLE-CELL RNA-SEQ

Ependymoma represents a heterogeneous disease affecting the entire neuraxis. Extensive molecular profiling efforts have identified molecular ependymoma subgroups based on DNA methylation. However, the intratumoral heterogeneity and developmental origins of these groups are only partially understood, and effective treatments are still lacking for about 50% of patients with high-risk tumors. We interrogated the cellular architecture of ependymoma using single cell/nucleus RNA-sequencing to analyze 24 tumor specimens across major molecular subgroups and anatomic locations. We additionally analyzed ten patient-derived ependymoma cell models and two patient-derived xenografts (PDXs). Interestingly, we identified an analogous cellular hierarchy across all ependymoma groups, originating from undifferentiated neural stem cell-like populations towards different degrees of impaired differentiation states comprising neuronal precursor-like, astro-glial-like, and ependymal-like tumor cells. While prognostically favorable ependymoma groups predominantly harbored differentiated cell populations, aggressive groups were enriched for undifferentiated subpopulations. Projection of transcriptomic signatures onto an independent bulk RNA-seq cohort stratified patient survival even within known molecular groups, thus refining the prognostic power of DNA methylation-based profiling. Furthermore, we identified novel potentially druggable targets including IGF- and FGF-signaling within poorly prognostic transcriptional programs. Ependymoma-derived cell models/PDXs widely recapitulated the transcriptional programs identified within fresh tumors and are leveraged to validate identified target genes in functional follow-up analyses. Taken together, our analyses reveal a developmental hierarchy and transcriptomic context underlying the biologically and clinically distinct behavior of ependymoma groups. The newly characterized cellular states and underlying regulatory networks could serve as basis for future therapeutic target identification and reveal biomarkers for clinical trials.

Abstract NG01: MADR: Rapid generation of somatic mosaics to model cancer in mice and human organoids

In situ transgenesis methods such as virus and electroporation can create somatic transgenic mice quickly, but they lack the exquisite control over copy number, zygosity, and locus specificity. We have recently established mosaic analysis by dual recombinase-mediated cassette exchange (MADR), which permits stable labeling of mutant cells expressing transgenic elements from precisely-defined chromosomal loci. MADR provides a toolkit of elements for combinatorial labeling, inducible/reversible transgene manipulation, VCre recombinase expression, and genetic manipulation of human cells. Further, we have demonstrated the versatility of MADR by creating glioma models with mixed, reporter-identified zygosity or with “personalized” driver mutations from pediatric glioma. For example, introducing H3f3a mutation variants with MADR regulates the spatiotemporal profile of glioma, and single-cell RNA and ATAC sequencing analysis demonstrates a recapitulation of developmental hierarchy seen in K27M-mutant human glioma. Moreover, we have generated novel models of supratentorial ependymoma using patient-derived oncofusion transgenes. These models display a high degree of fidelity and we now compare these models on a single-cell level with our previous models and human tumor cell transcriptomes. In addition, we now demonstrate the ability to generalize the MADR technology to other non-CNS tissues using local plasmid delivery. Finally, we have engineered human cells to allow for MADR transgenesis and somatic transgenic organoids. These combined approaches will enable researchers to discovery disease mechanisms and test therapeutics in more physiologically relevant cancer models.MADR is extensible to thousands of existing mouse lines and can be adapted to human cells, providing a flexible platform to democratize the generation of somatic transgenic disease models.

Citation Format: Gi Bum Kim, David Rincon Fernandez Pacheco, David Saxon, Amy Yang, Sara Sabat, Marina Dutra-Clarke, Rachelle Levy, Ashley Watkins, Hannah Park, Aslam Abbasi Akhtar, Paul W. Linesch, Naomi Kobritz, Swasty S. Shandra, Katie Grausam, Alberto Ayala-Sarmiento, Jessica Molina, Kristyna Sedivakova, Daniel S. Gareau, Mariella G. Filbin, Serguei Bannykh, Jie Tang, Mario L. Suva, Bin Chen, Moise Danielpour, Joshua J. Breunig. MADR: Rapid generation of somatic mosaics to model cancer in mice and human organoids [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr NG01.

Rapid Generation of Somatic Mouse Mosaics with Locus-Specific, Stably Integrated Transgenic Elements

In situ transgenesis methods such as viruses and electroporation can rapidly create somatic transgenic mice but lack control over copy number, zygosity, and locus specificity. Here we establish mosaic analysis by dual recombinase-mediated cassette exchange (MADR), which permits stable labeling of mutant cells expressing transgenic elements from precisely defined chromosomal loci. We provide a toolkit of MADR elements for combination labeling, inducible and reversible transgene manipulation, VCre recombinase expression, and transgenesis of human cells. Further, we demonstrate the versatility of MADR by creating glioma models with mixed reporter-identified zygosity or with "personalized" driver mutations from pediatric glioma. MADR is extensible to thousands of existing mouse lines, providing a flexible platform to democratize the generation of somatic mosaic mice. VIDEO ABSTRACT.

A Public BCR Present in a Unique Dual-Receptor-Expressing Lymphocyte from Type 1 Diabetes Patients Encodes a Potent T Cell Autoantigen

T and B cells are the two known lineages of adaptive immune cells. Here, we describe a previously unknown lymphocyte that is a dual expresser (DE) of TCR and BCR and key lineage markers of both B and T cells. In type 1 diabetes (T1D), DEs are predominated by one clonotype that encodes a potent CD4 T cell autoantigen in its antigen binding site. Molecular dynamics simulations revealed that this peptide has an optimal binding register for diabetogenic HLA-DQ8. In concordance, a synthetic version of the peptide forms stable DQ8 complexes and potently stimulates autoreactive CD4 T cells from T1D patients, but not healthy controls. Moreover, mAbs bearing this clonotype are autoreactive against CD4 T cells and inhibit insulin tetramer binding to CD4 T cells. Thus, compartmentalization of adaptive immune cells into T and B cells is not absolute, and violators of this paradigm are likely key drivers of autoimmune diseases.

Re-programing Chromatin with a Bifunctional LSD1/HDAC Inhibitor Induces Therapeutic Differentiation in DIPG

H3K27M mutations resulting in epigenetic dysfunction are frequently observed in diffuse intrinsic pontine glioma (DIPGs), an incurable pediatric cancer. We conduct a CRISPR screen revealing that knockout of KDM1A encoding lysine-specific demethylase 1 (LSD1) sensitizes DIPG cells to histone deacetylase (HDAC) inhibitors. Consistently, Corin, a bifunctional inhibitor of HDACs and LSD1, potently inhibits DIPG growth in vitro and in xenografts. Mechanistically, Corin increases H3K27me3 levels suppressed by H3K27M histones, and simultaneously increases HDAC-targeted H3K27ac and LSD1-targeted H3K4me1 at differentiation-associated genes. Corin treatment induces cell death, cell-cycle arrest, and a cellular differentiation phenotype and drives transcriptional changes correlating with increased survival time in DIPG patients. These data suggest a strategy for treating DIPG by simultaneously inhibiting LSD1 and HDACs.

Histone Variant and Cell Context Determine H3K27M Reprogramming of the Enhancer Landscape and Oncogenic State

Development of effective targeted cancer therapies is fundamentally limited by our molecular understanding of disease pathogenesis. Diffuse intrinsic pontine glioma (DIPG) is a fatal malignancy of the childhood pons characterized by a unique substitution to methionine in histone H3 at lysine 27 (H3K27M) that results in globally altered epigenetic marks and oncogenic transcription. Through primary DIPG tumor characterization and isogenic oncohistone expression, we show that the same H3K27M mutation displays distinct modes of oncogenic reprogramming and establishes distinct enhancer architecture depending upon both the variant of histone H3 and the cell context in which the mutation occurs. Compared with non-malignant pediatric pontine tissue, we identify and functionally validate both shared and variant-specific pathophysiology. Altogether, we provide a powerful resource of epigenomic data in 25 primary DIPG samples and 5 rare normal pediatric pontine tissue samples, revealing clinically relevant functional distinctions previously unidentified in DIPG.

Abstract 3647: Single cell RNA sequencing reveals mitogenic and progenitor gene programs inBRAF-rearranged pilocytic astrocytomas

Single-cell RNA sequencing (scRNAseq) has been performed across a range of intermediate to high-grade gliomas that each harbor multiple driver mutations. Pilocytic astrocytoma (PA), the most common childhood brain tumor, is a low-grade glioma. PAs frequently harbor oncogenic KIAA1549-BRAF fusions, and exhibit low rates of other driver mutations. While PAs exhibit a favorable prognosis compared to the higher-grade gliomas, treatment morbidity and tumor recurrence can represent major challenges for some PA patients. We performed scRNAseq of both tumor and non-tumor cells in newly-diagnosed PAs that contained KIAA1549-BRAF rearrangements. To confidently distinguish tumor cells from nontumor cells, we sorted cells by glial progenitor marker A2B5 status and profiled KIAA1549-BRAF fusion status of cells using a sensitive quantitative PCR-based assay. Results were validated using RNA in situ hybridization. When compared to higher-grade gliomas, a higher proportion of the PA tumor cells exhibited a differentiated, astrocyte-like phenotype. A smaller proportion of cells exhibited a progenitor-like phenotype with evidence of proliferation. These progenitor-like tumor cells expressed a mitogen-activated protein kinase (MAPK) program that was absent from higher-grade gliomas. Similar patterns of expression of genes associated with the astrocyte-like and MAPK gene programs were also seen in formalin fixed, paraffin embedded PA tissues using RNA in situ hybridization. Immune cells, especially microglia, comprised almost half of all cells in the PAs and accounted for differences in bulk expression profiles between tumor locations and subtypes. These single cell transcriptional data reveal a transcriptional developmental hierarchy in a pediatric low grade glioma that is skewed towards more mature brain cells compared to higher-grade gliomas. The results indicate that future analyses of bulk PA tissues should attempt to account for considerable infiltration by nontumor cells. Finally, the finding that a MAPK gene program is not uniformly expressed in PA tumor cells has implications for ongoing clinical investigations of therapies directed at the MAPK pathway for the treatment of PA.

Citation Format: Zachary J. Reitman, Brenton Paolella, Guillaume Bergthold, Kristine Pelton, Robert Jones, Sarah Becker, Claire E. Sinai, Haley Malkin, Ying Huang, Leslie Grimmett, Zachary T. Herbert, Yu Sun, Jessica Weatherbee, John Alberta, John Daley, Orit Rozenblatt-Rosen, Rosalind Segal, Daphne Haas-Kogan, Mariella G. Filbin, Mario L. Suva, Aviv Regev, Charles Stiles, Mark W. Kieran, Liliana Goumnerova, Keith L. Ligon, Alex K. Shalek, Pratiti Bandopadhayay, Rameen Beroukhim. Single cell RNA sequencing reveals mitogenic and progenitor gene programs inBRAF-rearranged pilocytic astrocytomas [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 3647.

Resolving the phylogenetic origin of glioblastoma via multifocal genomic analysis of pre-treatment and treatment-resistant autopsy specimens

Glioblastomas are malignant neoplasms composed of diverse cell populations. This intratumoral diversity has an underlying architecture, with a hierarchical relationship through clonal evolution from a common ancestor. Therapies are limited by emergence of resistant subclones from this phylogenetic reservoir. To characterize this clonal ancestral origin of recurrent tumors, we determined phylogenetic relationships using whole exome sequencing of pre-treatment IDH1/2 wild-type glioblastoma specimens, matched to post-treatment autopsy samples (n = 9) and metastatic extracranial post-treatment autopsy samples (n = 3). We identified “truncal” genetic events common to the evolutionary ancestry of the initial specimen and later recurrences, thereby inferring the identity of the precursor cell population. Mutations were identified in a subset of cases in known glioblastoma genes such as NF1(n = 3), TP53(n = 4) and EGFR(n = 5). However, by phylogenetic analysis, there were no protein-coding mutations as recurrent truncal events across the majority of cases. In contrast, whole copy-loss of chromosome 10 (12 of 12 cases), copy-loss of chromosome 9p21 (11 of 12 cases) and copy-gain in chromosome 7 (10 of 12 cases) were identified as shared events in the majority of cases. Strikingly, mutations in the TERT promoter were also identified as shared events in all evaluated pairs (9 of 9). Thus, we define four truncal non-coding genomic alterations that represent early genomic events in gliomagenesis, that identify the persistent cellular reservoir from which glioblastoma recurrences emerge. Therapies to target these key early genomic events are needed. These findings offer an evolutionary explanation for why precision therapies that target protein-coding mutations lack efficacy in GBM.

Non-coding and structural alterations may be early drivers of brain cancer development. A team led by Priscilla Brastianos and Tracy Batchelor from Massachusetts General Hospital, Boston, USA, analyzed the genetic landscape of glioblastoma by comparing pre-treatment and autopsy tumor specimens from 12 patients who died of the aggressive brain cancer. They identified a common set of four genetic events that occurred early in the evolution of nearly every patient’s cancer: three losses or gains of chromosome regions or entire chromosomes, and mutations in the gene-activating promoter of TERT, which encodes an enzyme implicated in the cancer’s growth. The findings help explain why therapies that target protein-coding mutations don’t work in brain cancer when they do in other tumor types. They also point to new drug targets.

Dissecting hematopoietic and renal cell heterogeneity in adult zebrafish at single-cell resolution using RNA sequencing

The work by Tang et al. provides a comprehensive, single-cell, transcriptomic analysis of kidney and blood cells from the adult zebrafish, identifying novel cell types, including two classes of NK immune cells, classically defined and erythroid-primed hematopoietic stem and progenitor cells, mucin-secreting kidney cells, and kidney stem/progenitor cells.

Recent advances in single-cell, transcriptomic profiling have provided unprecedented access to investigate cell heterogeneity during tissue and organ development. In this study, we used massively parallel, single-cell RNA sequencing to define cell heterogeneity within the zebrafish kidney marrow, constructing a comprehensive molecular atlas of definitive hematopoiesis and functionally distinct renal cells found in adult zebrafish. Because our method analyzed blood and kidney cells in an unbiased manner, our approach was useful in characterizing immune-cell deficiencies within DNA–protein kinase catalytic subunit (prkdc), interleukin-2 receptor γ a (il2rga), and double-homozygous–mutant fish, identifying blood cell losses in T, B, and natural killer cells within specific genetic mutants. Our analysis also uncovered novel cell types, including two classes of natural killer immune cells, classically defined and erythroid-primed hematopoietic stem and progenitor cells, mucin-secreting kidney cells, and kidney stem/progenitor cells. In total, our work provides the first, comprehensive, single-cell, transcriptomic analysis of kidney and marrow cells in the adult zebrafish.