Epigenetic encoding, heritability and plasticity of glioma transcriptional cell states

Cell of origin epigenetic priming determines susceptibility to Tet2 mutation

Hematopoietic stem cell (HSC) mutations can result in clonal hematopoiesis (CH) with heterogeneous clinical outcomes. Here, we investigate how the cell state preceding Tet2 mutation impacts the pre-malignant phenotype. Using an inducible system for clonal analysis of myeloid progenitors, we find that the epigenetic features of clones at similar differentiation status are highly heterogeneous and functionally respond differently to Tet2 mutation. Cell differentiation stage also influences Tet2 mutation response indicating that the cell of origin's epigenome modulates clone-specific behaviors in CH. Molecular features associated with higher risk outcomes include Sox4 that sensitizes cells to Tet2 inactivation, inducing dedifferentiation, altered metabolism and increasing the in vivo clonal output of mutant cells, as confirmed in primary GMP and HSC models. Our findings validate the hypothesis that epigenetic features can predispose specific clones for dominance, explaining why identical genetic mutations can result in different phenotypes.

Glial and glioneuronal tumors: Navigating the complexity of evolving concepts and new classification

The World Health Organization (WHO) published the 5th edition classification of tumors of central nervous system in 2021, commonly abbreviated as WHO CNS5, which became the new standard for brain tumor diagnosis and therapy. This edition dramatically impacted tumor diagnostics. In short it introduced new tumors, changed the names of previously recognized tumors, and modified the working definition of previously known tumors. The new system appears complex due to the integration of morphological and multiple molecular criteria. The most radical changes occurred in the field of glial and glioneuronal tumors, which constitutes the lengthy first chapter of this new edition. Herein we present an illustrative outline of the evolving concepts of glial and glioneuronal tumors. We also attempt to explain the rationales behind this substantial change in tumor classification and the challenges to update and integrate it into clinical practice. We aim to present a concise and precise roadmap to aid navigation through the intricate conceptual framework of glial and glioneuronal tumors in the context of WHO CNS5.

Mutant IDH inhibitors induce lineage differentiation in IDH-mutant oligodendroglioma

A subset of patients with IDH-mutant glioma respond to inhibitors of mutant IDH (IDHi), yet the molecular underpinnings of such responses are not understood. Here, we profiled by single-cell or single-nucleus RNA-sequencing three IDH-mutant oligodendrogliomas from patients who derived clinical benefit from IDHi. Importantly, the tissues were sampled on-drug, four weeks from treatment initiation. We further integrate our findings with analysis of single-cell and bulk transcriptomes from independent cohorts and experimental models. We find that IDHi treatment induces a robust differentiation toward the astrocytic lineage, accompanied by a depletion of stem-like cells and a reduction of cell proliferation. Furthermore, mutations in NOTCH1 are associated with decreased astrocytic differentiation and may limit the response to IDHi. Our study highlights the differentiating potential of IDHi on the cellular hierarchies that drive oligodendrogliomas and suggests a genetic modifier that may improve patient stratification.

Cellular diversity through space and time: adding new dimensions to GBM therapeutic development

The current median survival for glioblastoma (GBM) patients is only about 16 months, with many patients succumbing to the disease in just a matter of months, making it the most common and aggressive primary brain cancer in adults. This poor outcome is, in part, due to the lack of new treatment options with only one FDA-approved treatment in the last decade. Advances in sequencing techniques and transcriptomic analyses have revealed a vast degree of heterogeneity in GBM, from inter-patient diversity to intra-tumoral cellular variability. These cutting-edge approaches are providing new molecular insights highlighting a critical role for the tumor microenvironment (TME) as a driver of cellular plasticity and phenotypic heterogeneity. With this expanded molecular toolbox, the influence of TME factors, including endogenous (e.g., oxygen and nutrient availability and interactions with non-malignant cells) and iatrogenically induced (e.g., post-therapeutic intervention) stimuli, on tumor cell states can be explored to a greater depth. There exists a critical need for interrogating the temporal and spatial aspects of patient tumors at a high, cell-level resolution to identify therapeutically targetable states, interactions and mechanisms. In this review, we discuss advancements in our understanding of spatiotemporal diversity in GBM with an emphasis on the influence of hypoxia and immune cell interactions on tumor cell heterogeneity. Additionally, we describe specific high-resolution spatially resolved methodologies and their potential to expand the impact of pre-clinical GBM studies. Finally, we highlight clinical attempts at targeting hypoxia- and immune-related mechanisms of malignancy and the potential therapeutic opportunities afforded by single-cell and spatial exploration of GBM patient specimens.

Cancer stem cell hypothesis 2.0 in glioblastoma: Where are we now and where are we going?

Over the past 2 decades, the cancer stem cell (CSC) hypothesis has provided insight into many malignant tumors, including glioblastoma (GBM). Cancer stem cells have been identified in patient-derived tumors and in some mouse models, allowing for a deeper understanding of cellular and molecular mechanisms underlying GBM growth and therapeutic resistance. The CSC hypothesis has been the cornerstone of cellular heterogeneity, providing a conceptual and technical framework to explain this longstanding phenotype in GBM. This hypothesis has evolved to fit recent insights into how cellular plasticity drives tumor growth to suggest that CSCs do not represent a distinct population but rather a cellular state with substantial plasticity that can be achieved by non-CSCs under specific conditions. This has further been reinforced by advances in genomics, including single-cell approaches, that have used the CSC hypothesis to identify multiple putative CSC states with unique properties, including specific developmental and metabolic programs. In this review, we provide a historical perspective on the CSC hypothesis and its recent evolution, with a focus on key functional phenotypes, and provide an update on the definition for its use in future genomic studies.

Long-term Multimodal Recording Reveals Epigenetic Adaptation Routes in Dormant Breast Cancer Cells

Patients with estrogen receptor-positive breast cancer receive adjuvant endocrine therapies (ET) that delay relapse by targeting clinically undetectable micrometastatic deposits. Yet, up to 50% of patients relapse even decades after surgery through unknown mechanisms likely involving dormancy. To investigate genetic and transcriptional changes underlying tumor awakening, we analyzed late relapse patients and longitudinally profiled a rare cohort treated with long-term neoadjuvant ETs until progression. Next, we developed an in vitro evolutionary study to record the adaptive strategies of individual lineages in unperturbed parallel experiments. Our data demonstrate that ETs induce nongenetic cell state transitions into dormancy in a stochastic subset of cells via epigenetic reprogramming. Single lineages with divergent phenotypes awaken unpredictably in the absence of recurrent genetic alterations. Targeting the dormant epigenome shows promising activity against adapting cancer cells. Overall, this study uncovers the contribution of epigenetic adaptation to the evolution of resistance to ETs.

SIGNIFICANCE

This study advances the understanding of therapy-induced dormancy with potential clinical implications for breast cancer. Estrogen receptor-positive breast cancer cells adapt to endocrine treatment by entering a dormant state characterized by strong heterochromatinization with no recurrent genetic changes. Targeting the epigenetic rewiring impairs the adaptation of cancer cells to ETs. See related commentary by Llinas-Bertran et al., p. 704. This article is featured in Selected Articles from This Issue, p. 695.

Epigenetic regulation of tumor-immune symbiosis in glioma

Glioma is a type of aggressive and incurable brain tumor. Patients with glioma are highly resistant to all types of therapies, including immunotherapies. Epigenetic reprogramming is a key molecular hallmark in tumors across cancer types, including glioma. Mounting evidence highlights a pivotal role of epigenetic regulation in shaping tumor biology and therapeutic responses through mechanisms involving both glioma cells and immune cells, as well as their symbiotic interactions in the tumor microenvironment (TME). In this review, we discuss the molecular mechanisms of epigenetic regulation that impacts glioma cell biology and tumor immunity in both a cell-autonomous and non-cell-autonomous manner. Moreover, we provide an overview of potential therapeutic approaches that can disrupt epigenetic-regulated tumor-immune symbiosis in the glioma TME.

VC-resist glioblastoma cell state: vessel co-option as a key driver of chemoradiation resistance

Glioblastoma (GBM) is a highly lethal type of cancer. GBM recurrence following chemoradiation is typically attributed to the regrowth of invasive and resistant cells. Therefore, there is a pressing need to gain a deeper understanding of the mechanisms underlying GBM resistance to chemoradiation and its ability to infiltrate. Using a combination of transcriptomic, proteomic, and phosphoproteomic analyses, longitudinal imaging, organotypic cultures, functional assays, animal studies, and clinical data analyses, we demonstrate that chemoradiation and brain vasculature induce cell transition to a functional state named VC-Resist (vessel co-opting and resistant cell state). This cell state is midway along the transcriptomic axis between proneural and mesenchymal GBM cells and is closer to the AC/MES1-like state. VC-Resist GBM cells are highly vessel co-opting, allowing significant infiltration into the surrounding brain tissue and homing to the perivascular niche, which in turn induces even more VC-Resist transition. The molecular and functional characteristics of this FGFR1-YAP1-dependent GBM cell state, including resistance to DNA damage, enrichment in the G2M phase, and induction of senescence/stemness pathways, contribute to its enhanced resistance to chemoradiation. These findings demonstrate how vessel co-option, perivascular niche, and GBM cell plasticity jointly drive resistance to therapy during GBM recurrence.

Somatic epimutations enable single-cell lineage tracing in native hematopoiesis across the murine and human lifespan

Current approaches to lineage tracing of stem cell clones require genetic engineering or rely on sparse somatic DNA variants, which are difficult to capture at single-cell resolution. Here, we show that targeted single-cell measurements of DNA methylation at single-CpG resolution deliver joint information about cellular differentiation state and clonal identities. We develop EPI-clone, a droplet-based method for transgene-free lineage tracing, and apply it to study hematopoiesis, capturing hundreds of clonal trajectories across almost 100,000 single-cells. Using ground-truth genetic barcodes, we demonstrate that EPI-clone accurately identifies clonal lineages throughout hematopoietic differentiation. Applied to unperturbed hematopoiesis, we describe an overall decline of clonal complexity during murine ageing and the expansion of rare low-output stem cell clones. In aged human donors, we identified expanded hematopoietic clones with and without genetic lesions, and various degrees of clonal complexity. Taken together, EPI-clone enables accurate and transgene-free single-cell lineage tracing at scale.

Single-cell nanobiopsy enables multigenerational longitudinal transcriptomics of cancer cells

Single-cell RNA sequencing has revolutionized our understanding of cellular heterogeneity, but routine methods require cell lysis and fail to probe the dynamic trajectories responsible for cellular state transitions, which can only be inferred. Here, we present a nanobiopsy platform that enables the injection of exogenous molecules and multigenerational longitudinal cytoplasmic sampling from a single cell and its progeny. The technique is based on scanning ion conductance microscopy (SICM) and, as a proof of concept, was applied to longitudinally profile the transcriptome of single glioblastoma (GBM) brain tumor cells in vitro over 72 hours. The GBM cells were biopsied before and after exposure to chemotherapy and radiotherapy, and our results suggest that treatment either induces or selects for more transcriptionally stable cells. We envision the nanobiopsy will contribute to transforming standard single-cell transcriptomics from a static analysis into a dynamic assay.

The Epigenetic Evolution of Glioma Is Determined by the IDH1 Mutation Status and Treatment Regimen

Tumor adaptation or selection is thought to underlie therapy resistance in glioma. To investigate longitudinal epigenetic evolution of gliomas in response to therapeutic pressure, we performed an epigenomic analysis of 132 matched initial and recurrent tumors from patients with IDH-wildtype (IDHwt) and IDH-mutant (IDHmut) glioma. IDHwt gliomas showed a stable epigenome over time with relatively low levels of global methylation. The epigenome of IDHmut gliomas showed initial high levels of genome-wide DNA methylation that was progressively reduced to levels similar to those of IDHwt tumors. Integration of epigenomics, gene expression, and functional genomics identified HOXD13 as a master regulator of IDHmut astrocytoma evolution. Furthermore, relapse of IDHmut tumors was accompanied by histologic progression that was associated with survival, as validated in an independent cohort. Finally, the initial cell composition of the tumor microenvironment varied between IDHwt and IDHmut tumors and changed differentially following treatment, suggesting increased neoangiogenesis and T-cell infiltration upon treatment of IDHmut gliomas. This study provides one of the largest cohorts of paired longitudinal glioma samples with epigenomic, transcriptomic, and genomic profiling and suggests that treatment of IDHmut glioma is associated with epigenomic evolution toward an IDHwt-like phenotype.

SIGNIFICANCE

Standard treatments are related to loss of DNA methylation in IDHmut glioma, resulting in epigenetic activation of genes associated with tumor progression and alterations in the microenvironment that resemble treatment-naïve IDHwt glioma.

Mechanistic insight into glioma through spatially multidimensional proteomics

Characterizing the tumor microenvironment at the molecular level is essential for understanding the mechanisms of tumorigenesis and evolution. However, the specificity of the blood proteome in localized region of the tumor and its linkages with other systems is difficult to investigate. Here, we propose a spatially multidimensional comparative proteomics strategy using glioma as an example. The blood proteome signature of tumor microenvironment was specifically identified by in situ collection of arterial and venous blood from the glioma region of the brain for comparison with peripheral blood. Also, by integrating with different dimensions of tissue and peripheral blood proteomics, the information on the genesis, migration, and exchange of glioma-associated proteins was revealed, which provided a powerful method for tumor mechanism research and biomarker discovery. The study recruited multidimensional clinical cohorts, allowing the proteomic results to corroborate each other, reliably revealing biological processes specific to gliomas, and identifying highly accurate biomarkers.

IDHwt glioblastomas can be stratified by their transcriptional response to standard treatment, with implications for targeted therapy

BACKGROUND

Glioblastoma (GBM) brain tumors lacking IDH1 mutations (IDHwt) have the worst prognosis of all brain neoplasms. Patients receive surgery and chemoradiotherapy but tumors almost always fatally recur.

RESULTS

Using RNA sequencing data from 107 pairs of pre- and post-standard treatment locally recurrent IDHwt GBM tumors, we identify two responder subtypes based on longitudinal changes in gene expression. In two thirds of patients, a specific subset of genes is upregulated from primary to recurrence (Up responders), and in one third, the same genes are downregulated (Down responders), specifically in neoplastic cells. Characterization of the responder subtypes indicates subtype-specific adaptive treatment resistance mechanisms that are associated with distinct changes in the tumor microenvironment. In Up responders, recurrent tumors are enriched in quiescent proneural GBM stem cells and differentiated neoplastic cells, with increased interaction with the surrounding normal brain and neurotransmitter signaling, whereas Down responders commonly undergo mesenchymal transition. ChIP-sequencing data from longitudinal GBM tumors suggests that the observed transcriptional reprogramming could be driven by Polycomb-based chromatin remodeling rather than DNA methylation.

CONCLUSIONS

We show that the responder subtype is cancer-cell intrinsic, recapitulated in in vitro GBM cell models, and influenced by the presence of the tumor microenvironment. Stratifying GBM tumors by responder subtype may lead to more effective treatment.

Glioblastoma evolution and heterogeneity from a 3D whole-tumor perspective

Treatment failure for the lethal brain tumor glioblastoma (GBM) is attributed to intratumoral heterogeneity and tumor evolution. We utilized 3D neuronavigation during surgical resection to acquire samples representing the whole tumor mapped by 3D spatial coordinates. Integrative tissue and single-cell analysis revealed sources of genomic, epigenomic, and microenvironmental intratumoral heterogeneity and their spatial patterning. By distinguishing tumor-wide molecular features from those with regional specificity, we inferred GBM evolutionary trajectories from neurodevelopmental lineage origins and initiating events such as chromothripsis to emergence of genetic subclones and spatially restricted activation of differential tumor and microenvironmental programs in the core, periphery, and contrast-enhancing regions. Our work depicts GBM evolution and heterogeneity from a 3D whole-tumor perspective, highlights potential therapeutic targets that might circumvent heterogeneity-related failures, and establishes an interactive platform enabling 360° visualization and analysis of 3D spatial patterns for user-selected genes, programs, and other features across whole GBM tumors.

Cancer Evolution: A Multifaceted Affair

Cancer cells adapt and survive through the acquisition and selection of molecular modifications. This process defines cancer evolution. Building on a theoretical framework based on heritable genetic changes has provided insights into the mechanisms supporting cancer evolution. However, cancer hallmarks also emerge via heritable nongenetic mechanisms, including epigenetic and chromatin topological changes, and interactions between tumor cells and the tumor microenvironment. Recent findings on tumor evolutionary mechanisms draw a multifaceted picture where heterogeneous forces interact and influence each other while shaping tumor progression. A comprehensive characterization of the cancer evolutionary toolkit is required to improve personalized medicine and biomarker discovery.

SIGNIFICANCE

Tumor evolution is fueled by multiple enabling mechanisms. Importantly, genetic instability, epigenetic reprogramming, and interactions with the tumor microenvironment are neither alternative nor independent evolutionary mechanisms. As demonstrated by findings highlighted in this perspective, experimental and theoretical approaches must account for multiple evolutionary mechanisms and their interactions to ultimately understand, predict, and steer tumor evolution.

Characterization of prevalent tyrosine kinase inhibitors and their challenges in glioblastoma treatment

Glioblastoma multiforme (GBM) is a highly aggressive malignant primary tumor in the central nervous system. Despite extensive efforts in radiotherapy, chemotherapy, and neurosurgery, there remains an inadequate level of improvement in treatment outcomes. The development of large-scale genomic and proteomic analysis suggests that GBMs are characterized by transcriptional heterogeneity, which is responsible for therapy resistance. Hence, knowledge about the genetic and epigenetic heterogeneity of GBM is crucial for developing effective treatments for this aggressive form of brain cancer. Tyrosine kinases (TKs) can act as signal transducers, regulate important cellular processes like differentiation, proliferation, apoptosis and metabolism. Therefore, TK inhibitors (TKIs) have been developed to specifically target these kinases. TKIs are categorized into allosteric and non-allosteric inhibitors. Irreversible inhibitors form covalent bonds, which can lead to longer-lasting effects. However, this can also increase the risk of off-target effects and toxicity. The development of TKIs as therapeutics through computer-aided drug design (CADD) and bioinformatic techniques enhance the potential to improve patients' survival rates. Therefore, the continued exploration of TKIs as drug targets is expected to lead to even more effective and specific therapeutics in the future.

Identification and validation of chromatin regulator-related signatures as a novel prognostic model for low-grade gliomas using translational bioinformatics

AIMS

The purpose of this study is to explore the potential biological role and prognostic significance of chromatin regulators (CRs) in low-grade gliomas (LGGs).

MAIN METHODS

CRs were obtained from the FACER database. Transcription profiles of LGG patients were collected from the TCGA and CGGA databases. Differentially expressed CRs (DECRs) between LGGs and normal controls were identified using DESeq2. The consensus clustering algorithm was employed to distinguish subtypes of LGGs based on prognosis-related DECRs. The differences in clinical and molecular characteristics between different subtypes were explored. R packages, GSVA, ssGSEA, and ESTIMATE were utilized to elucidate the tumor microenvironment and activated pathways in different subtypes. Subsequently, a CRs-related signature was developed using LASSO Cox regression. Its performance was evaluated by Kaplan-Meier curve and ROC curve analyses. In vitro experiments were performed to explore the function of JADE3 in LGGs, which predominantly expressed in glioma cells.

KEY FINDINGS

We identified 43 DECRs and two CRs-related subtypes of LGGs. The subtype characterized by shorter survival displayed significant enrichment for pathways associated with DNA damage response and repair, along with heightened immune cell infiltration. Furthermore, the CRs-based signature exhibited excellent prognostic performance in both the TCGA and CGGA databases. Knockdown of JADE3 significantly increased the invasion, migration, and proliferation abilities of Hs683.

SIGNIFICANCE

Our study reveals the aberrant expression and prognostic value of CRs in LGGs. It emphasizes the potential regulatory role of CRs in the microenvironment and DNA damage repair in LGGs. JADE3 could be a possible therapeutic target for LGGs.

Transcriptional state dynamics lead to heterogeneity and adaptive tumor evolution in urothelial bladder carcinoma

Intra-tumor heterogeneity contributes to treatment failure and poor survival in urothelial bladder carcinoma (UBC). Analyzing transcriptome from a UBC cohort, we report that intra-tumor transcriptomic heterogeneity indicates co-existence of tumor cells in epithelial and mesenchymal-like transcriptional states and bi-directional transition between them occurs within and between tumor subclones. We model spontaneous and reversible transition between these partially heritable states in cell lines and characterize their population dynamics. SMAD3, KLF4 and PPARG emerge as key regulatory markers of the transcriptional dynamics. Nutrient limitation, as in the core of large tumors, and radiation treatment perturb the dynamics, initially selecting for a transiently resistant phenotype and then reconstituting heterogeneity and growth potential, driving adaptive evolution. Dominance of transcriptional states with low PPARG expression indicates an aggressive phenotype in UBC patients. We propose that phenotypic plasticity and dynamic, non-genetic intra-tumor heterogeneity modulate both the trajectory of disease progression and adaptive treatment response in UBC.

Alterations in cellular metabolism under different grades of glioma staging identified based on a multi-omics analysis strategy

Glioma is a type of brain tumor closely related to abnormal cell metabolism. Firstly, multiple combinatorial sequencing studies have revealed this relationship. Genomic studies have identified gene mutations and gene expression disorders related to the development of gliomas, which affect cell metabolic pathways. In addition, transcriptome studies have revealed the genes and regulatory networks that regulate cell metabolism in glioma tissues. Metabonomics studies have shown that the metabolic pathway of glioma cells has changed, indicating their distinct energy and nutritional requirements. This paper focuses on the retrospective analysis of multiple groups combined with sequencing to analyze the changes in various metabolites during metabolism in patients with glioma. Finally, the changes in genes, regulatory networks, and metabolic pathways regulating cell metabolism in patients with glioma under different metabolic conditions were discussed. It is also proposed that multi-group metabolic analysis is expected to better understand the mechanism of abnormal metabolism of gliomas and provide more personalized methods and guidance for early diagnosis, treatment, and prognosis evaluation of gliomas.

Phenotypic noise and plasticity in cancer evolution

Non-genetic alterations can produce changes in a cell's phenotype. In cancer, these phenomena can influence a cell's fitness by conferring access to heritable, beneficial phenotypes. Herein, we argue that current discussions of 'phenotypic plasticity' in cancer evolution ignore a salient feature of the original definition: namely, that it occurs in response to an environmental change. We suggest 'phenotypic noise' be used to distinguish non-genetic changes in phenotype that occur independently from the environment. We discuss the conceptual and methodological techniques used to identify these phenomena during cancer evolution. We propose that the distinction will guide efforts to define mechanisms of phenotype change, accelerate translational work to manipulate phenotypes through treatment, and, ultimately, improve patient outcomes.

Disrupting cellular memory to overcome drug resistance

Gene expression states persist for varying lengths of time at the single-cell level, a phenomenon known as gene expression memory. When cells switch states, losing memory of their prior state, this transition can occur in the absence of genetic changes. However, we lack robust methods to find regulators of memory or track state switching. Here, we develop a lineage tracing-based technique to quantify memory and identify cells that switch states. Applied to melanoma cells without therapy, we quantify long-lived fluctuations in gene expression that are predictive of later resistance to targeted therapy. We also identify the PI3K and TGF-β pathways as state switching modulators. We propose a pretreatment model, first applying a PI3K inhibitor to modulate gene expression states, then applying targeted therapy, which leads to less resistance than targeted therapy alone. Together, we present a method for finding modulators of gene expression memory and their associated cell fates.

Glioblastoma revisited: from neuronal-like invasion to pacemaking

In recent years, two developments have helped us to better understand the fundamental biology of glioblastoma: the description of a striking intratumoral heterogeneity including gene expression-based cell states, and the discovery that neuro-cancer interactions and cancer-intrinsic neurodevelopmental mechanisms are fundamental features of glioblastoma. In this opinion article, we aim to integrate both developments. We explain how two key disease features are characterized by different neural mechanisms related to distinct but plastic cancer cell states: first, the single cell-dominated invasive parts and second, the more solid parts which are dominated by communicating cell networks constantly activated by pacemaker-like glioblastoma cells. The resulting integrative roadmap of molecular and functional heterogeneity contributes to the Cancer Neuroscience of glioblastoma and suggests novel therapeutic strategies.

Human cerebrospinal fluid affects chemoradiotherapy sensitivities in tumor cells from patients with glioblastoma

Cancers in the central nervous system resist therapies effective in other cancers, possibly due to the unique biochemistry of the human brain microenvironment composed of cerebrospinal fluid (CSF). However, the impact of CSF on cancer cells and therapeutic efficacy is unknown. Here, we examined the effect of human CSF on glioblastoma (GBM) tumors from 25 patients. We found that CSF induces tumor cell plasticity and resistance to standard GBM treatments (temozolomide and irradiation). We identified nuclear protein 1 (NUPR1), a transcription factor hampering ferroptosis, as a mediator of therapeutic resistance in CSF. NUPR1 inhibition with a repurposed antipsychotic, trifluoperazine, enhanced the killing of GBM cells resistant to chemoradiation in CSF. The same chemo-effective doses of trifluoperazine were safe for human neurons and astrocytes derived from pluripotent stem cells. These findings reveal that chemoradiation efficacy decreases in human CSF and suggest that combining trifluoperazine with standard care may improve the survival of patients with GBM.

Cancer vaccine strategies for the treatment of diffusely infiltrating gliomas

Diffusely infiltrating gliomas - including glioblastoma (GBM), isocitrate dehydrogenase (IDH) mutant gliomas, and histone 3 (H3) altered gliomas - are primary brain tumors with an invariably fatal outcome. Despite advances in the understanding of their biology, standard, targeted and immune checkpoint inhibitor immunotherapies have proven ineffective in arresting their inexorable progression and associated morbidity and mortality. Recognizing the unique aspects of the immunogenicity of cancer cells, the last decade has seen the development and evaluation of vaccine-based therapies for the treatment of solid tumors, including gliomas. Here we review the current vaccine strategies for the treatment of GBM, IDH-mutant gliomas and diffuse midline glioma H3 K27M-altered. We discuss potential benefits and challenges of vaccine therapies in these specific patient populations.

Decoding Cancer Evolution: Integrating Genetic and Non-Genetic Insights

The development of cancer begins with cells transitioning from their multicellular nature to a state akin to unicellular organisms. This shift leads to a breakdown in the crucial regulators inherent to multicellularity, resulting in the emergence of diverse cancer cell subpopulations that have enhanced adaptability. The presence of different cell subpopulations within a tumour, known as intratumoural heterogeneity (ITH), poses challenges for cancer treatment. In this review, we delve into the dynamics of the shift from multicellularity to unicellularity during cancer onset and progression. We highlight the role of genetic and non-genetic factors, as well as tumour microenvironment, in promoting ITH and cancer evolution. Additionally, we shed light on the latest advancements in omics technologies that allow for in-depth analysis of tumours at the single-cell level and their spatial organization within the tissue. Obtaining such detailed information is crucial for deepening our understanding of the diverse evolutionary paths of cancer, allowing for the development of effective therapies targeting the key drivers of cancer evolution.

Dysregulation of inflammasome activation in glioma

Gliomas are the most common brain tumors characterized by complicated heterogeneity. The genetic, molecular, and histological pathology of gliomas is characterized by high neuro-inflammation. The inflammatory microenvironment in the central nervous system (CNS) has been closely linked with inflammasomes that control the inflammatory response and coordinate innate host defenses. Dysregulation of the inflammasome causes an abnormal inflammatory response, leading to carcinogenesis in glioma. Because of the clinical importance of the various physiological properties of the inflammasome in glioma, the inflammasome has been suggested as a promising treatment target for glioma management. Here, we summarize the current knowledge on the contribution of the inflammasomes in glioma and therapeutic insights. Video Abstract.

Single-cell multi-omics defines the cell-type-specific impact of splicing aberrations in human hematopoietic clonal outgrowths

RNA splicing factors are recurrently mutated in clonal blood disorders, but the impact of dysregulated splicing in hematopoiesis remains unclear. To overcome technical limitations, we integrated genotyping of transcriptomes (GoT) with long-read single-cell transcriptomics and proteogenomics for single-cell profiling of transcriptomes, surface proteins, somatic mutations, and RNA splicing (GoT-Splice). We applied GoT-Splice to hematopoietic progenitors from myelodysplastic syndrome (MDS) patients with mutations in the core splicing factor SF3B1. SF3B1mut cells were enriched in the megakaryocytic-erythroid lineage, with expansion of SF3B1mut erythroid progenitor cells. We uncovered distinct cryptic 3' splice site usage in different progenitor populations and stage-specific aberrant splicing during erythroid differentiation. Profiling SF3B1-mutated clonal hematopoiesis samples revealed that erythroid bias and cell-type-specific cryptic 3' splice site usage in SF3B1mut cells precede overt MDS. Collectively, GoT-Splice defines the cell-type-specific impact of somatic mutations on RNA splicing, from early clonal outgrowths to overt neoplasia, directly in human samples.

Multi-omic integration of DNA methylation and gene expression data reveals molecular vulnerabilities in glioblastoma

Glioblastoma (GBM) is one of the most aggressive types of cancer and exhibits profound genetic and epigenetic heterogeneity, making the development of an effective treatment a major challenge. The recent incorporation of molecular features into the diagnosis of patients with GBM has led to an improved categorization into various tumour subtypes with different prognoses and disease management. In this work, we have exploited the benefits of genome-wide multi-omic approaches to identify potential molecular vulnerabilities existing in patients with GBM. Integration of gene expression and DNA methylation data from both bulk GBM and patient-derived GBM stem cell lines has revealed the presence of major sources of GBM variability, pinpointing subtype-specific tumour vulnerabilities amenable to pharmacological interventions. In this sense, inhibition of the AP-1, SMAD3 and RUNX1/RUNX2 pathways, in combination or not with the chemotherapeutic agent temozolomide, led to the subtype-specific impairment of tumour growth, particularly in the context of the aggressive, mesenchymal-like subtype. These results emphasize the involvement of these molecular pathways in the development of GBM and have potential implications for the development of personalized therapeutic approaches.

Cancer stem cells in brain tumors: From origin to clinical implications

Malignant brain tumors are highly heterogeneous tumors with a poor prognosis and a high morbidity and mortality rate in both children and adults. The cancer stem cell (CSC, also named tumor-initiating cell) model states that tumor growth is driven by a subset of CSCs. This model explains some of the clinical observations of brain tumors, including the almost unavoidable tumor recurrence after initial successful chemotherapy and/or radiotherapy and treatment resistance. Over the past two decades, strategies for the identification and characterization of brain CSCs have improved significantly, supporting the design of new diagnostic and therapeutic strategies for brain tumors. Relevant studies have unveiled novel characteristics of CSCs in the brain, including their heterogeneity and distinctive immunobiology, which have provided opportunities for new research directions and potential therapeutic approaches. In this review, we summarize the current knowledge of CSCs markers and stemness regulators in brain tumors. We also comprehensively describe the influence of the CSCs niche and tumor microenvironment on brain tumor stemness, including interactions between CSCs and the immune system, and discuss the potential application of CSCs in brain-based therapies for the treatment of brain tumors.

Glioblastoma modeling with 3D organoids: progress and challenges

Glioblastoma (GBM) is the most aggressive adult primary brain tumor with nearly universal treatment resistance and recurrence. The mainstay of therapy remains maximal safe surgical resection followed by concurrent radiation therapy and temozolomide chemotherapy. Despite intensive investigation, alternative treatment options, such as immunotherapy or targeted molecular therapy, have yielded limited success to achieve long-term remission. This difficulty is partly due to the lack of pre-clinical models that fully recapitulate the intratumoral and intertumoral heterogeneity of GBM and the complex tumor microenvironment. Recently, GBM 3D organoids originating from resected patient tumors, genetic manipulation of induced pluripotent stem cell (iPSC)-derived brain organoids and bio-printing or fusion with non-malignant tissues have emerged as novel culture systems to portray the biology of GBM. Here, we highlight several methodologies for generating GBM organoids and discuss insights gained using such organoid models compared to classic modeling approaches using cell lines and xenografts. We also outline limitations of current GBM 3D organoids, most notably the difficulty retaining the tumor microenvironment, and discuss current efforts for improvements. Finally, we propose potential applications of organoid models for a deeper mechanistic understanding of GBM and therapeutic development.

β-Catenin-Driven Differentiation Is a Tissue-Specific Epigenetic Vulnerability in Adrenal Cancer

Adrenocortical carcinoma (ACC) is a rare cancer in which tissue-specific differentiation is paradoxically associated with dismal outcomes. The differentiated ACC subtype CIMP-high is prevalent, incurable, and routinely fatal. CIMP-high ACC possess abnormal DNA methylation and frequent β-catenin-activating mutations. Here, we demonstrated that ACC differentiation is maintained by a balance between nuclear, tissue-specific β-catenin-containing complexes, and the epigenome. On chromatin, β-catenin bound master adrenal transcription factor SF1 and hijacked the adrenocortical super-enhancer landscape to maintain differentiation in CIMP-high ACC; off chromatin, β-catenin bound histone methyltransferase EZH2. SF1/β-catenin and EZH2/β-catenin complexes present in normal adrenals persisted through all phases of ACC evolution. Pharmacologic EZH2 inhibition in CIMP-high ACC expelled SF1/β-catenin from chromatin and favored EZH2/β-catenin assembly, erasing differentiation and restraining cancer growth in vitro and in vivo. These studies illustrate how tissue-specific programs shape oncogene selection, surreptitiously encoding targetable therapeutic vulnerabilities.

SIGNIFICANCE

Oncogenic β-catenin can use tissue-specific partners to regulate cellular differentiation programs that can be reversed by epigenetic therapies, identifying epigenetic control of differentiation as a viable target for β-catenin-driven cancers.

Chimeric Antigen Receptor T Cells in Glioblastoma-Current Concepts and Promising Future

Glioblastoma (GBM) is a highly aggressive primary brain tumor that is largely refractory to treatment and, therefore, invariably relapses. GBM patients have a median overall survival of 15 months and, given this devastating prognosis, there is a high need for therapy improvement. One of the therapeutic approaches currently tested in GBM is chimeric antigen receptor (CAR)-T cell therapy. CAR-T cells are genetically altered T cells that are redirected to eliminate tumor cells in a highly specific manner. There are several challenges to CAR-T cell therapy in solid tumors such as GBM, including restricted trafficking and penetration of tumor tissue, a highly immunosuppressive tumor microenvironment (TME), as well as heterogeneous antigen expression and antigen loss. In addition, CAR-T cells have limitations concerning safety, toxicity, and the manufacturing process. To date, CAR-T cells directed against several target antigens in GBM including interleukin-13 receptor alpha 2 (IL-13Rα2), epidermal growth factor receptor variant III (EGFRvIII), human epidermal growth factor receptor 2 (HER2), and ephrin type-A receptor 2 (EphA2) have been tested in preclinical and clinical studies. These studies demonstrated that CAR-T cell therapy is a feasible option in GBM with at least transient responses and acceptable adverse effects. Further improvements in CAR-T cells regarding their efficacy, flexibility, and safety could render them a promising therapy option in GBM.

Single-cell methylation sequencing data reveal succinct metastatic migration histories and tumor progression models

Recent studies exploring the impact of methylation in tumor evolution suggest that although the methylation status of many of the CpG sites are preserved across distinct lineages, others are altered as the cancer progresses. Because changes in methylation status of a CpG site may be retained in mitosis, they could be used to infer the progression history of a tumor via single-cell lineage tree reconstruction. In this work, we introduce the first principled distance-based computational method, Sgootr, for inferring a tumor's single-cell methylation lineage tree and for jointly identifying lineage-informative CpG sites that harbor changes in methylation status that are retained along the lineage. We apply Sgootr on single-cell bisulfite-treated whole-genome sequencing data of multiregionally sampled tumor cells from nine metastatic colorectal cancer patients, as well as multiregionally sampled single-cell reduced-representation bisulfite sequencing data from a glioblastoma patient. We show that the tumor lineages constructed reveal a simple model underlying tumor progression and metastatic seeding. A comparison of Sgootr against alternative approaches shows that Sgootr can construct lineage trees with fewer migration events and with more in concordance with the sequential-progression model of tumor evolution, with a running time a fraction of that used in prior studies. Lineage-informative CpG sites identified by Sgootr are in inter-CpG island (CGI) regions, as opposed to intra-CGIs, which have been the main regions of interest in genomic methylation-related analyses.

Haploinsufficiency of NFKBIA reshapes the epigenome antipodal to the IDH mutation and imparts disease fate in diffuse gliomas

Genetic alterations help predict the clinical behavior of diffuse gliomas, but some variability remains uncorrelated. Here, we demonstrate that haploinsufficient deletions of chromatin-bound tumor suppressor NFKB inhibitor alpha (NFKBIA) display distinct patterns of occurrence in relation to other genetic markers and are disproportionately present at recurrence. NFKBIA haploinsufficiency is associated with unfavorable patient outcomes, independent of genetic and clinicopathologic predictors. NFKBIA deletions reshape the DNA and histone methylome antipodal to the IDH mutation and induce a transcriptome landscape partly reminiscent of H3K27M mutant pediatric gliomas. In IDH mutant gliomas, NFKBIA deletions are common in tumors with a clinical course similar to that of IDH wild-type tumors. An externally validated nomogram model for estimating individual patient survival in IDH mutant gliomas confirms that NFKBIA deletions predict comparatively brief survival. Thus, NFKBIA haploinsufficiency aligns with distinct epigenome changes, portends a poor prognosis, and should be incorporated into models predicting the disease fate of diffuse gliomas.

Applications of single-cell RNA sequencing in drug discovery and development

Single-cell technologies, particularly single-cell RNA sequencing (scRNA-seq) methods, together with associated computational tools and the growing availability of public data resources, are transforming drug discovery and development. New opportunities are emerging in target identification owing to improved disease understanding through cell subtyping, and highly multiplexed functional genomics screens incorporating scRNA-seq are enhancing target credentialling and prioritization. ScRNA-seq is also aiding the selection of relevant preclinical disease models and providing new insights into drug mechanisms of action. In clinical development, scRNA-seq can inform decision-making via improved biomarker identification for patient stratification and more precise monitoring of drug response and disease progression. Here, we illustrate how scRNA-seq methods are being applied in key steps in drug discovery and development, and discuss ongoing challenges for their implementation in the pharmaceutical industry.

Considerations for modelling diffuse high-grade gliomas and developing clinically relevant therapies

Diffuse high-grade gliomas contain some of the most dangerous human cancers that lack curative treatment options. The recent molecular stratification of gliomas by the World Health Organisation in 2021 is expected to improve outcomes for patients in neuro-oncology through the development of treatments targeted to specific tumour types. Despite this promise, research is hindered by the lack of preclinical modelling platforms capable of recapitulating the heterogeneity and cellular phenotypes of tumours residing in their native human brain microenvironment. The microenvironment provides cues to subsets of glioma cells that influence proliferation, survival, and gene expression, thus altering susceptibility to therapeutic intervention. As such, conventional in vitro cellular models poorly reflect the varied responses to chemotherapy and radiotherapy seen in these diverse cellular states that differ in transcriptional profile and differentiation status. In an effort to improve the relevance of traditional modelling platforms, recent attention has focused on human pluripotent stem cell-based and tissue engineering techniques, such as three-dimensional (3D) bioprinting and microfluidic devices. The proper application of these exciting new technologies with consideration of tumour heterogeneity and microenvironmental interactions holds potential to develop more applicable models and clinically relevant therapies. In doing so, we will have a better chance of translating preclinical research findings to patient populations, thereby addressing the current derisory oncology clinical trial success rate.

macroH2A2 antagonizes epigenetic programs of stemness in glioblastoma

Self-renewal is a crucial property of glioblastoma cells that is enabled by the choreographed functions of chromatin regulators and transcription factors. Identifying targetable epigenetic mechanisms of self-renewal could therefore represent an important step toward developing effective treatments for this universally lethal cancer. Here we uncover an epigenetic axis of self-renewal mediated by the histone variant macroH2A2. With omics and functional assays deploying patient-derived in vitro and in vivo models, we show that macroH2A2 shapes chromatin accessibility at enhancer elements to antagonize transcriptional programs of self-renewal. macroH2A2 also sensitizes cells to small molecule-mediated cell death via activation of a viral mimicry response. Consistent with these results, our analyses of clinical cohorts indicate that high transcriptional levels of this histone variant are associated with better prognosis of high-grade glioma patients. Our results reveal a targetable epigenetic mechanism of self-renewal controlled by macroH2A2 and suggest additional treatment approaches for glioblastoma patients.

Effect of OR51E1 Single Nucleotide Polymorphisms on Glioma Susceptibility in the Chinese Han Population

BACKGROUND

Glioma is one of the common primary intracranial tumors, which is heterogeneous among individuals with a low cure rate. Our study aimed to investigate the association between single nucleotide polymorphisms (SNPs) of the OR51E1 gene and glioma susceptibility in the Chinese Han population.

METHODS

A total of six SNPs on OR51E1 in 1,026 subjects (526 cases and 500 controls) were genotyped by MassARRAY iPLEX GOLD assay. The association between these SNPs and glioma susceptibility was analyzed using logistic regression, and odds ratios (ORs) and 95% confidence intervals (CIs) were also calculated. The multifactor dimensionality reduction (MDR) method was applied to detect "SNP-SNP" interactions.

RESULTS

In the overall sample, polymorphisms rs10768148, rs7102992, and rs10500608 were identified to be associated with glioma risk. In the stratified analysis based on gender, only polymorphism rs10768148 was observed to be associated with the risk of glioma. In the age-stratified analysis, rs7102992, rs74052483, and rs10500609 contributed to the risk of glioma in subjects aged > 40 years. And polymorphisms rs10768148 and rs7102992 were associated with the risk of glioma in subjects aged ≤ 40 years and subjects with astrocytoma. In addition, a strong synergistic relationship between rs74052483 and rs10768148, and a strong redundant relationship between rs7102992 and rs10768148 were identified in the study.

CONCLUSIONS

This study demonstrated the association of OR51E1 polymorphisms with glioma susceptibility, providing a basis for assessing glioma risk-associated variants in the Chinese Han population.

Microenvironment-Driven Dynamic Chromatin Changes in Glioblastoma Recapitulate Early Neural Development at Single-Cell Resolution

The tumor microenvironment is necessary for recapitulating the intratumoral heterogeneity and cell state plasticity found in human primary glioblastoma (GBM). Conventional models do not accurately recapitulate the spectrum of GBM cellular states, hindering elucidation of the underlying transcriptional regulation of these states. Using our glioblastoma cerebral organoid model, we profiled the chromatin accessibility of 28,040 single cells in five patient-derived glioma stem cell lines. Integration of paired epigenomes and transcriptomes within the context of tumor-normal host cell interactions was used to probe the gene-regulatory networks underlying individual GBM cellular states in a way not readily possible in other in vitro models. These analyses identified the epigenetic underpinnings of GBM cellular states and characterized dynamic chromatin changes reminiscent of early neural development that underlie GBM cell state transitions. Despite large differences between tumors, a shared cellular compartment made up of neural progenitor-like cells and outer radial glia-like cells was observed. Together, these results shed light on the transcriptional regulation program in GBM and offer novel therapeutic targets across a broad range of genetically heterogenous GBMs.

SIGNIFICANCE

Single-cell analyses elucidate the chromatin landscape and transcriptional regulation of glioblastoma cellular states and identify a radial glia-like population, providing potential targets to disrupt cell states and improve therapeutic efficacy.

Applying single cell multi-omic analyses to understand treatment resistance in pediatric high grade glioma

Despite improvements in cancer patient outcomes seen in the past decade, tumor resistance to therapy remains a major impediment to achieving durable clinical responses. Intratumoral heterogeneity related to genetic, epigenetic, transcriptomic, proteomic, and metabolic differences between individual cancer cells has emerged as a driver of therapeutic resistance. This cell to cell heterogeneity can be assessed using single cell profiling technologies that enable the identification of tumor cell clones that exhibit similar defining features like specific mutations or patterns of DNA methylation. Single cell profiling of tumors before and after treatment can generate new insights into the cancer cell characteristics that confer therapeutic resistance by identifying intrinsically resistant sub-populations that survive treatment and by describing new cellular features that emerge post-treatment due to tumor cell evolution. Integrative, single cell analytical approaches have already proven advantageous in studies characterizing treatment-resistant clones in cancers where pre- and post-treatment patient samples are readily available, such as leukemia. In contrast, little is known about other cancer subtypes like pediatric high grade glioma, a class of heterogeneous, malignant brain tumors in children that rapidly develop resistance to multiple therapeutic modalities, including chemotherapy, immunotherapy, and radiation. Leveraging single cell multi-omic technologies to analyze naïve and therapy-resistant glioma may lead to the discovery of novel strategies to overcome treatment resistance in brain tumors with dismal clinical outcomes. In this review, we explore the potential for single cell multi-omic analyses to reveal mechanisms of glioma resistance to therapy and discuss opportunities to apply these approaches to improve long-term therapeutic response in pediatric high grade glioma and other brain tumors with limited treatment options.

Extensive intratumor regional epigenetic heterogeneity in clear cell renal cell carcinoma targets kidney enhancers and is associated with poor outcome

BACKGROUND

Clear cell renal cell cancer (ccRCC), the 8th leading cause of cancer-related death in the US, is challenging to treat due to high level intratumoral heterogeneity (ITH) and the paucity of druggable driver mutations. CcRCC is unusual for its high frequency of epigenetic regulator mutations, such as the SETD2 histone H3 lysine 36 trimethylase (H3K36me3), and low frequency of traditional cancer driver mutations. In this work, we examined epigenetic level ITH and defined its relationships with pathologic features, aspects of tumor biology, and SETD2 mutations.

RESULTS

A multi-region sampling approach coupled with EPIC DNA methylation arrays was conducted on a cohort of normal kidney and ccRCC. ITH was assessed using DNA methylation (5mC) and CNV-based entropy and Euclidian distances. We found elevated 5mC heterogeneity and entropy in ccRCC relative to normal kidney. Variable CpGs are highly enriched in enhancer regions. Using intra-class correlation coefficient analysis, we identified CpGs that segregate tumor regions according to clinical phenotypes related to tumor aggressiveness. SETD2 wild-type tumors overall possess greater 5mC and copy number ITH than SETD2 mutant tumor regions, suggesting SETD2 loss contributes to a distinct epigenome. Finally, coupling our regional data with TCGA, we identified a 5mC signature that links regions within a primary tumor with metastatic potential.

CONCLUSION

Taken together, our results reveal marked levels of epigenetic ITH in ccRCC that are linked to clinically relevant tumor phenotypes and could translate into novel epigenetic biomarkers.

LonP1 Drives Proneural Mesenchymal Transition in IDH1-R132H Diffuse Glioma

Malignant astroctyoma and glioblastoma are diffuse CNS tumors that have markedly similar features, including microvascular proliferation and necrosis, and the latter presents higher grade and poorer survival. The Isocitrate dehydrogenase 1/2 (IDH) mutation further predicts improved survival and is present in oligodendroglioma and astrocytoma. The latter are more prevalent in younger populations with a median age of 37 years at diagnosis as compared to glioblastoma with a median age of 641,2. These tumors frequently have co-occurring ATRX and/or TP53 mutations (Brat et al., 2021). The IDH mutation is known to cause dysregulation of the hypoxia response broadly in CNS tumors and subsequent reduction in both tumor growth and treatment resistance. The frequency of tumor recurrence is high for diffuse CNS tumors. Understanding the mechanism and potential molecular targets enhancing treatment resistance and local invasion in IDH mutant diffuse glioma is necessary for developing new treatment strategies for better tumor control and improving overall survival. Recent evidence highlights the importance of local foci in IDH mutant glioma with an accelerated stress response as responsible for recurrence in these tumors. Here, we demonstrate that LonP1 drives NRF2 and subsequent proneural mesenchymal transition interdependent with the IDH mutation in response to stress and other tumor microenvironment cues. Our findings provide further evidence that targeting LonP1 may be a crucial strategy for improving the standard-of-care treatment in IDH mutant diffuse astrocytoma.

Unraveling non-genetic heterogeneity in cancer with dynamical models and computational tools

Individual cells within an otherwise genetically homogenous population constantly undergo fluctuations in their molecular state, giving rise to non-genetic heterogeneity. Such diversity is being increasingly implicated in cancer therapy resistance and metastasis. Identifying the origins of non-genetic heterogeneity is therefore crucial for making clinical breakthroughs. We discuss with examples how dynamical models and computational tools have provided critical multiscale insights into the nature and consequences of non-genetic heterogeneity in cancer. We demonstrate how mechanistic modeling has been pivotal in establishing key concepts underlying non-genetic diversity at various biological scales, from population dynamics to gene regulatory networks. We discuss advances in single-cell longitudinal profiling techniques to reveal patterns of non-genetic heterogeneity, highlighting the ongoing efforts and challenges in statistical frameworks to robustly interpret such multimodal datasets. Moving forward, we stress the need for data-driven statistical and mechanistically motivated dynamical frameworks to come together to develop predictive cancer models and inform therapeutic strategies.

Single-cell heterogeneity of EGFR and CDK4 co-amplification is linked to immune infiltration in glioblastoma

Glioblastoma (GBM) is the most aggressive brain tumor, with a median survival of ∼15 months. Targeted approaches have not been successful in this tumor type due to the large extent of intratumor heterogeneity. Mosaic amplification of oncogenes suggests that multiple genetically distinct clones are present in each tumor. To uncover the relationships between genetically diverse subpopulations of GBM cells and their native tumor microenvironment, we employ highly multiplexed spatial protein profiling coupled with single-cell spatial mapping of fluorescence in situ hybridization (FISH) for EGFR, CDK4, and PDGFRA. Single-cell FISH analysis of a total of 35,843 single nuclei reveals that tumors in which amplifications of EGFR and CDK4 more frequently co-occur in the same cell exhibit higher infiltration of CD163+ immunosuppressive macrophages. Our results suggest that high-throughput assessment of genomic alterations at the single-cell level could provide a measure for predicting the immune state of GBM.

Neuronal and tumourigenic boundaries of glioblastoma plasticity

Glioblastoma (GBM) remains the most lethal primary brain cancer largely due to recurrence of treatment-resistant disease. Current therapies are ultimately ineffective as GBM tumour cells adapt their identity to escape treatment. Recent advances in single-cell epigenetics and transcriptomics highlight heterogeneous cell populations in GBM tumours originating from unique cancerous genetic aberrations. However, they also suggest that tumour cells conserve molecular properties of parent neuronal cells, with their permissive epigenetic profiles enabling them to morph along a finite number of reprogramming routes to evade treatment. Here, we review the known tumourigenic, neurodevelopmental and brain-injury boundaries of GBM plasticity, and propose that effective treatment of GBM requires the addition of therapeutics that restrain GBM plasticity.

Identification, validation and biological characterisation of novel glioblastoma tumour microenvironment subtypes: implications for precision immunotherapy

BACKGROUND

New precision medicine therapies are urgently required for glioblastoma (GBM). However, to date, efforts to subtype patients based on molecular profiles have failed to direct treatment strategies. We hypothesised that interrogation of the GBM tumour microenvironment (TME) and identification of novel TME-specific subtypes could inform new precision immunotherapy treatment strategies.

MATERIALS AND METHODS

A refined and validated microenvironment cell population (MCP) counter method was applied to >800 GBM patient tumours (GBM-MCP-counter). Specifically, partition around medoids (PAM) clustering of GBM-MCP-counter scores in the GLIOTRAIN discovery cohort identified three novel patient clusters, uniquely characterised by TME composition, functional orientation markers and immune checkpoint proteins. Validation was carried out in three independent GBM-RNA-seq datasets. Neoantigen, mutational and gene ontology analysis identified mutations and uniquely altered pathways across subtypes. The longitudinal Glioma Longitudinal AnalySiS (GLASS) cohort and three immunotherapy clinical trial cohorts [treatment with neoadjuvant/adjuvant anti-programmed cell death protein 1 (PD-1) or PSVRIPO] were further interrogated to assess subtype alterations between primary and recurrent tumours and to assess the utility of TME classifiers as immunotherapy biomarkers.

RESULTS

TME^High tumours (30%) displayed elevated lymphocyte, myeloid cell immune checkpoint, programmed cell death protein 1 (PD-1) and cytotoxic T-lymphocyte-associated protein 4 transcripts. TME^High/mesenchymal+ patients featured tertiary lymphoid structures. TME^Med (46%) tumours were enriched for endothelial cell gene expression profiles and displayed heterogeneous immune populations. TME^Low (24%) tumours were manifest as an 'immune-desert' group. TME subtype transitions upon recurrence were identified in the longitudinal GLASS cohort. Assessment of GBM immunotherapy trial datasets revealed that TME^High patients receiving neoadjuvant anti-PD-1 had significantly increased overall survival (P = 0.04). Moreover, TME^High patients treated with adjuvant anti-PD-1 or oncolytic virus (PVSRIPO) showed a trend towards improved survival.

CONCLUSIONS

We have established a novel TME-based classification system for application in intracranial malignancies. TME subtypes represent canonical 'termini a quo' (starting points) to support an improved precision immunotherapy treatment approach.

A variational algorithm to detect the clonal copy number substructure of tumors from scRNA-seq data

Single-cell RNA sequencing is the reference technology to characterize the composition of the tumor microenvironment and to study tumor heterogeneity at high resolution. Here we report Single CEll Variational ANeuploidy analysis (SCEVAN), a fast variational algorithm for the deconvolution of the clonal substructure of tumors from single-cell RNA-seq data. It uses a multichannel segmentation algorithm exploiting the assumption that all the cells in a given copy number clone share the same breakpoints. Thus, the smoothed expression profile of every individual cell constitutes part of the evidence of the copy number profile in each subclone. SCEVAN can automatically and accurately discriminate between malignant and non-malignant cells, resulting in a practical framework to analyze tumors and their microenvironment. We apply SCEVAN to datasets encompassing 106 samples and 93,322 cells from different tumor types and technologies. We demonstrate its application to characterize the intratumor heterogeneity and geographic evolution of malignant brain tumors.

Epigenetics as a mediator of plasticity in cancer

The concept of an epigenetic landscape describing potential cellular fates arising from pluripotent cells, first advanced by Conrad Waddington, has evolved in light of experiments showing nondeterministic outcomes of regulatory processes and mathematical methods for quantifying stochasticity. In this Review, we discuss modern approaches to epigenetic and gene regulation landscapes and the associated ideas of entropy and attractor states, illustrating how their definitions are both more precise and relevant to understanding cancer etiology and the plasticity of cancerous states. We address the interplay between different types of regulatory landscapes and how their changes underlie cancer progression. We also consider the roles of cellular aging and intrinsic and extrinsic stimuli in modulating cellular states and how landscape alterations can be quantitatively mapped onto phenotypic outcomes and thereby used in therapy development.

Scalable co-sequencing of RNA and DNA from individual nuclei

The ideal technology for directly investigating the relationship between genotype and phenotype would analyze both RNA and DNA genome-wide and with single-cell resolution. However, existing tools lack the throughput required for comprehensive analysis of complex tumors and tissues. We introduce a highly scalable method for jointly profiling DNA and expression following nucleosome depletion (DEFND-seq). In DEFND-seq, nuclei are nucleosome-depleted, tagmented, and separated into individual droplets for mRNA and genomic DNA barcoding. Once nuclei have been depleted of nucleosomes, subsequent steps can be performed using the widely available 10x Genomics droplet microfluidic technology and commercial kits without experimental modification. We demonstrate the production of high-complexity mRNA and gDNA sequencing libraries from thousands of individual nuclei from both cell lines and archived surgical specimens for associating gene expression phenotypes with both copy number and single nucleotide variants.