An atlas of causal and mechanistic drivers of interpatient heterogeneity in glioma.

Grade IV glioma, formerly known as glioblastoma multiforme (GBM) is the most aggressive and lethal type of brain tumor, and its treatment remains challenging in part due to extensive interpatient heterogeneity in disease driving mechanisms and lack of prognostic and predictive biomarkers. Using mechanistic inference of node-edge relationship (MINER), we have analyzed multiomics profiles from 516 patients and constructed an atlas of causal and mechanistic drivers of interpatient heterogeneity in GBM (gbmMINER). The atlas has delineated how 30 driver mutations act in a combinatorial scheme to causally influence a network of regulators (306 transcription factors and 73 miRNAs) of 179 transcriptional "programs", influencing disease progression in patients across 23 disease states. Through extensive testing on independent patient cohorts, we share evidence that a machine learning model trained on activity profiles of programs within gbmMINER significantly augments risk stratification, identifying patients who are super-responders to standard of care and those that would benefit from 2 nd line treatments. In addition to providing mechanistic hypotheses regarding disease prognosis, the activity of programs containing targets of 2 nd line treatments accurately predicted efficacy of 28 drugs in killing glioma stem-like cells from 43 patients. Our findings demonstrate that interpatient heterogeneity manifests from differential activities of transcriptional programs, providing actionable strategies for mechanistically characterizing GBM from a systems perspective and developing better prognostic and predictive biomarkers for personalized medicine.

Gene regulatory network topology governs resistance and treatment escape in glioma stem-like cells.

Poor prognosis and drug resistance in glioblastoma (GBM) can result from cellular heterogeneity and treatment-induced shifts in phenotypic states of tumor cells, including dedifferentiation into glioma stem-like cells (GSCs). This rare tumorigenic cell subpopulation resists temozolomide, undergoes proneural-to-mesenchymal transition (PMT) to evade therapy, and drives recurrence. Through inference of transcriptional regulatory networks (TRNs) of patient-derived GSCs (PD-GSCs) at single-cell resolution, we demonstrate how the topology of transcription factor interaction networks drives distinct trajectories of cell state transitions in PD-GSCs resistant or susceptible to cytotoxic drug treatment. By experimentally testing predictions based on TRN simulations, we show that drug treatment drives surviving PD-GSCs along a trajectory of intermediate states, exposing vulnerability to potentiated killing by siRNA or a second drug targeting treatment-induced transcriptional programs governing non-genetic cell plasticity. Our findings demonstrate an approach to uncover TRN topology and use it to rationally predict combinatorial treatments that disrupts acquired resistance in GBM.

Interindividual variation in human cortical cell type abundance and expression.

Single-cell transcriptomic studies have identified a conserved set of neocortical cell types from small postmortem cohorts. We extended these efforts by assessing cell type variation across 75 adult individuals undergoing epilepsy and tumor surgeries. Nearly all nuclei map to one of 125 robust cell types identified in the middle temporal gyrus. However, we found interindividual variance in abundances and gene expression signatures, particularly in deep-layer glutamatergic neurons and microglia. A minority of donor variance is explainable by age, sex, ancestry, disease state, and cell state. Genomic variation was associated with expression of 150 to 250 genes for most cell types. This characterization of cellular variation provides a baseline for cell typing in health and disease.

An expanded safety/feasibility study of the EMulate Therapeutics Voyager™ System in patients with recurrent glioblastoma.

Aim: The EMulate Therapeutics Voyager™ is a simple, wearable, home-use device that uses an alternating electromagnetic field to alter biologic signaling within cells. Objective: To assess the safety/feasibility of the Voyager in the treatment of recurrent glioblastoma (rGBM). Methods: In this study, patients with rGBM were treated with Voyager as monotherapy or in combination with standard chemotherapy at the Investigator's discretion. Safety was assessed by incidence of adverse events associated with the Voyager. Patients were followed until death. Results: A total of 75 patients were enrolled and treated for at least one day with the Voyager (safety population). Device-related adverse events were uncommon and generally did not result in interruption or withdrawal from treatment. There were no serious adverse events associated with Voyager. A total of 60 patients were treated for at least one month (clinical utility population). The median progression-free survival (PFS) was 17 weeks (4.3 months) in the Voyager only group (n = 24) and 21 weeks (5.3 months) in the Voyager + concurrent therapy group (n = 36). The median overall survival (OS) was 7 months in the Voyager only group and 9 months in the Voyager + concurrent therapy group. In patients treated with Voyager + concurrent therapy, the median OS for patients enrolled with their 1st or 2nd recurrence (n = 26) was 10 months, while in patients enrolled with their 3rd or 4th recurrence (n = 10) OS was 7 months. Conclusion: The data support the safety and feasibility of the Voyager for the treatment of rGBM. Further prospective study of the device is warranted. Trial Registration Number: NCT02296580 (ClinicalTrials.gov).

The potential role of human endogenous retrovirus K in glioblastoma.

The most active human endogenous retrovirus K (HERV-K) subtype, HML-2, has been implicated as a driver of oncogenesis in several cancers. However, the presence and function of HML-2 in malignant gliomas has remained unclear. In this issue of the JCI, Shah and colleagues demonstrate HML-2 overexpression in glioblastoma (GBM) and its role in maintaining the cancer stem cell phenotype. Given that stem-like cells are considered responsible for GBM heterogeneity and treatment resistance, targeting the stem cell niche may reduce tumor recurrence and improve clinical outcomes. The findings provide a foundation for future studies to determine whether antiretroviral and/or immunotherapy approaches targeting HML-2 could be used as therapeutics for GBM.

Target cell-specific synaptic dynamics of excitatory to inhibitory neuron connections in supragranular layers of human neocortex.

Rodent studies have demonstrated that synaptic dynamics from excitatory to inhibitory neuron types are often dependent on the target cell type. However, these target cell-specific properties have not been well investigated in human cortex, where there are major technical challenges in reliably obtaining healthy tissue, conducting multiple patch-clamp recordings on inhibitory cell types, and identifying those cell types. Here, we take advantage of newly developed methods for human neurosurgical tissue analysis with multiple patch-clamp recordings, post-hoc fluorescent in situ hybridization (FISH), machine learning-based cell type classification and prospective GABAergic AAV-based labeling to investigate synaptic properties between pyramidal neurons and PVALB- vs. SST-positive interneurons. We find that there are robust molecular differences in synapse-associated genes between these neuron types, and that individual presynaptic pyramidal neurons evoke postsynaptic responses with heterogeneous synaptic dynamics in different postsynaptic cell types. Using molecular identification with FISH and classifiers based on transcriptomically identified PVALB neurons analyzed by Patch-seq, we find that PVALB neurons typically show depressing synaptic characteristics, whereas other interneuron types including SST-positive neurons show facilitating characteristics. Together, these data support the existence of target cell-specific synaptic properties in human cortex that are similar to rodent, thereby indicating evolutionary conservation of local circuit connectivity motifs from excitatory to inhibitory neurons and their synaptic dynamics.

Association of Autologous Tumor Lysate-Loaded Dendritic Cell Vaccination With Extension of Survival Among Patients With Newly Diagnosed and Recurrent Glioblastoma: A Phase 3 Prospective Externally Controlled Cohort Trial.

Importance: Glioblastoma is the most lethal primary brain cancer. Clinical outcomes for glioblastoma remain poor, and new treatments are needed. Objective: To investigate whether adding autologous tumor lysate-loaded dendritic cell vaccine (DCVax-L) to standard of care (SOC) extends survival among patients with glioblastoma. Design, Setting, and Participants: This phase 3, prospective, externally controlled nonrandomized trial compared overall survival (OS) in patients with newly diagnosed glioblastoma (nGBM) and recurrent glioblastoma (rGBM) treated with DCVax-L plus SOC vs contemporaneous matched external control patients treated with SOC. This international, multicenter trial was conducted at 94 sites in 4 countries from August 2007 to November 2015. Data analysis was conducted from October 2020 to September 2021. Interventions: The active treatment was DCVax-L plus SOC temozolomide. The nGBM external control patients received SOC temozolomide and placebo; the rGBM external controls received approved rGBM therapies. Main Outcomes and Measures: The primary and secondary end points compared overall survival (OS) in nGBM and rGBM, respectively, with contemporaneous matched external control populations from the control groups of other formal randomized clinical trials. Results: A total of 331 patients were enrolled in the trial, with 232 randomized to the DCVax-L group and 99 to the placebo group. Median OS (mOS) for the 232 patients with nGBM receiving DCVax-L was 19.3 (95% CI, 17.5-21.3) months from randomization (22.4 months from surgery) vs 16.5 (95% CI, 16.0-17.5) months from randomization in control patients (HR = 0.80; 98% CI, 0.00-0.94; P = .002). Survival at 48 months from randomization was 15.7% vs 9.9%, and at 60 months, it was 13.0% vs 5.7%. For 64 patients with rGBM receiving DCVax-L, mOS was 13.2 (95% CI, 9.7-16.8) months from relapse vs 7.8 (95% CI, 7.2-8.2) months among control patients (HR, 0.58; 98% CI, 0.00-0.76; P < .001). Survival at 24 and 30 months after recurrence was 20.7% vs 9.6% and 11.1% vs 5.1%, respectively. Survival was improved in patients with nGBM with methylated MGMT receiving DCVax-L compared with external control patients (HR, 0.74; 98% CI, 0.55-1.00; P = .03). Conclusions and Relevance: In this study, adding DCVax-L to SOC resulted in clinically meaningful and statistically significant extension of survival for patients with both nGBM and rGBM compared with contemporaneous, matched external controls who received SOC alone. Trial Registration: ClinicalTrials.gov Identifier: NCT00045968.

Agent-based computational modeling of glioblastoma predicts that stromal density is central to oncolytic virus efficacy.

Oncolytic viruses (OVs) are emerging cancer immunotherapy. Despite notable successes in the treatment of some tumors, OV therapy for central nervous system cancers has failed to show efficacy. We used an ex vivo tumor model developed from human glioblastoma tissue to evaluate the infiltration of herpes simplex OV rQNestin (oHSV-1) into glioblastoma tumors. We next leveraged our data to develop a computational, model of glioblastoma dynamics that accounts for cellular interactions within the tumor. Using our computational model, we found that low stromal density was highly predictive of oHSV-1 therapeutic success, suggesting that the efficacy of oHSV-1 in glioblastoma may be determined by stromal-to-tumor cell regional density. We validated these findings in heterogenous patient samples from brain metastatic adenocarcinoma. Our integrated modeling strategy can be applied to suggest mechanisms of therapeutic responses for central nervous system cancers and to facilitate the successful translation of OVs into the clinic.

Cell-free plasma microRNAs that identify patients with glioblastoma.

Glioblastoma (GBM) is still one of the most commonly diagnosed advanced stage primary brain tumors. Current treatments for patients with primary GBM (pGBM) are often not effective and a significant proportion of the patients with pGBM recur. The effective treatment options for recurrent GBM (rGBM) are limited and survival outcomes are poor. This retrospective multicenter pilot study aims to determine potential cell-free microRNAs (cfmiRs) that identify patients with pGBM and rGBM tumors. 2,083 miRs were assessed using the HTG miRNA whole transcriptome assay (WTA). CfmiRs detection was compared in pre-operative plasma samples from patients with pGBM (n = 32) and rGBM (n = 13) to control plasma samples from normal healthy donors (n = 73). 265 cfmiRs were found differentially expressed in plasma samples from pGBM patients compared to normal healthy donors (FDR < 0.05). Of those 193 miRs were also detected in pGBM tumor tissues (n = 15). Additionally, we found 179 cfmiRs differentially expressed in rGBM, of which 68 cfmiRs were commonly differentially expressed in pGBM. Using Random Forest algorithm, specific cfmiR classifiers were found in the plasma of pGBM, rGBM, and both pGBM and rGBM combined. Two common cfmiR classifiers, miR-3180-3p and miR-5739, were found in all the comparisons. In receiving operating characteristic (ROC) curves analysis for rGBM miR-3180-3p showed a specificity of 87.7% and a sensitivity of 100% (AUC = 98.5%); while miR-5739 had a specificity of 79.5% and sensitivity of 92.3% (AUC = 90.2%). This study demonstrated that plasma samples from pGBM and rGBM patients have specific miR signatures. CfmiR-3180-3p and cfmiR-5739 have potential utility in diagnosing patients with pGBM and rGBM tumors using a minimally invasive blood assay.

Signature morpho-electric, transcriptomic, and dendritic properties of human layer 5 neocortical pyramidal neurons.

In the neocortex, subcerebral axonal projections originate largely from layer 5 (L5) extratelencephalic-projecting (ET) neurons. The unique morpho-electric properties of these neurons have been mainly described in rodents, where retrograde tracers or transgenic lines can label them. Similar labeling strategies are infeasible in the human neocortex, rendering the translational relevance of findings in rodents unclear. We leveraged the recent discovery of a transcriptomically defined L5 ET neuron type to study the properties of human L5 ET neurons in neocortical brain slices derived from neurosurgeries. Patch-seq recordings, where transcriptome, physiology, and morphology were assayed from the same cell, revealed many conserved morpho-electric properties of human and rodent L5 ET neurons. Divergent properties were often subtler than differences between L5 cell types within these two species. These data suggest a conserved function of L5 ET neurons in the neocortical hierarchy but also highlight phenotypic divergence possibly related to functional specialization of human neocortex.

Detection of Human Papillomavirus Integration in Brain Metastases from Oropharyngeal Tumors by Targeted Sequencing.

Human papillomavirus (HPV) positive and negative head and neck squamous cell carcinoma (HNSCC) are known to have differential phenotypes, including the incidence and location of metastases. HPV positive (HPV+) HNSCC are more likely to metastasize to distant sites, such as the lung, brain, and skin. Among these locations, metastasis to the brain is a rare event, and little is known about specific risk factors for this phenotype. In this report, we describe two patients who developed brain metastases from HNSCC. Both patient tumors had p16INK4a overexpression, suggesting these tumors were HPV+. This was confirmed after PCR, in situ hybridization, and mass spectrometry detected the presence of HPV type 16 (HPV16) DNA, RNA and protein. To further characterize the presence of HPV16, we used a target enrichment strategy on tumor DNA and RNA to isolate the viral sequences from the brain metastases. Analysis by targeted next generation sequencing revealed that both tumors had the HPV genome integrated into the host genome at known hotspots, 8q24.21 and 14q24.1. Applying a similar target enrichment strategy to a larger cohort of HPV+ HNSCC brain metastases could help to identify biomarkers that can predict metastasis and/or identify novel therapeutic options.

A Systems Approach to Brain Tumor Treatment.

Brain tumors are among the most lethal tumors. Glioblastoma, the most frequent primary brain tumor in adults, has a median survival time of approximately 15 months after diagnosis or a five-year survival rate of 10%; the recurrence rate is nearly 90%. Unfortunately, this prognosis has not improved for several decades. The lack of progress in the treatment of brain tumors has been attributed to their high rate of primary therapy resistance. Challenges such as pronounced inter-patient variability, intratumoral heterogeneity, and drug delivery across the blood-brain barrier hinder progress. A comprehensive, multiscale understanding of the disease, from the molecular to the whole tumor level, is needed to address the intratumor heterogeneity resulting from the coexistence of a diversity of neoplastic and non-neoplastic cell types in the tumor tissue. By contrast, inter-patient variability must be addressed by subtyping brain tumors to stratify patients and identify the best-matched drug(s) and therapies for a particular patient or cohort of patients. Accomplishing these diverse tasks will require a new framework, one involving a systems perspective in assessing the immense complexity of brain tumors. This would in turn entail a shift in how clinical medicine interfaces with the rapidly advancing high-throughput (HTP) technologies that have enabled the omics-scale profiling of molecular features of brain tumors from the single-cell to the tissue level. However, several gaps must be closed before such a framework can fulfill the promise of precision and personalized medicine for brain tumors. Ultimately, the goal is to integrate seamlessly multiscale systems analyses of patient tumors and clinical medicine. Accomplishing this goal would facilitate the rational design of therapeutic strategies matched to the characteristics of patients and their tumors. Here, we discuss some of the technologies, methodologies, and computational tools that will facilitate the realization of this vision to practice.

Functional enhancer elements drive subclass-selective expression from mouse to primate neocortex.

Viral genetic tools that target specific brain cell types could transform basic neuroscience and targeted gene therapy. Here, we use comparative open chromatin analysis to identify thousands of human-neocortical-subclass-specific putative enhancers from across the genome to control gene expression in adeno-associated virus (AAV) vectors. The cellular specificity of reporter expression from enhancer-AAVs is established by molecular profiling after systemic AAV delivery in mouse. Over 30% of enhancer-AAVs produce specific expression in the targeted subclass, including both excitatory and inhibitory subclasses. We present a collection of Parvalbumin (PVALB) enhancer-AAVs that show highly enriched expression not only in cortical PVALB cells but also in some subcortical PVALB populations. Five vectors maintain PVALB-enriched expression in primate neocortex. These results demonstrate how genome-wide open chromatin data mining and cross-species AAV validation can be used to create the next generation of non-species-restricted viral genetic tools.

Glioma Stem Cells as Immunotherapeutic Targets: Advancements and Challenges.

Glioblastoma is the most common and lethal primary brain malignancy. Despite major investments in research into glioblastoma biology and drug development, treatment remains limited and survival has not substantially improved beyond 1-2 years. Cancer stem cells (CSC) or glioma stem cells (GSC) refer to a population of tumor originating cells capable of self-renewal and differentiation. While controversial and challenging to study, evidence suggests that GCSs may result in glioblastoma tumor recurrence and resistance to treatment. Multiple treatment strategies have been suggested at targeting GCSs, including immunotherapy, posttranscriptional regulation, modulation of the tumor microenvironment, and epigenetic modulation. In this review, we discuss recent advances in glioblastoma treatment specifically focused on targeting of GCSs as well as their potential integration into current clinical pathways and trials.

Low-Grade Gemistocytic Morphology in H3 G34R-Mutant Gliomas and Concurrent K27M Mutation: Clinicopathologic Findings.

Mutations in histone H3 are key molecular drivers of pediatric and young adult high-grade gliomas. Histone H3 G34R mutations occur in hemispheric high-grade gliomas and H3 K27M mutations occur in aggressive, though histologically diverse, midline gliomas. Here, we report 2 rare cases of histologically low-grade gliomas with gemistocytic morphology and sequencing-confirmed histone H3 G34R mutations. One case is a histologically low-grade gemistocytic astrocytoma with a G34R-mutation in H3F3A. The second case is a histologically low-grade gemistocytic astrocytoma with co-occurring K27M and G34R mutations in HIST1H3B. Review of prior histone H3-mutant gliomas sequenced at our institution shows a divergent clinical and immunohistochemical pattern in the 2 cases. The first case is similar to prior histone H3 G34R-mutant tumors, while the second case most closely resembles prior histone H3 K27M-mutant gliomas. These represent novel cases of sequencing-confirmed histone H3 G34R-mutant gliomas with low-grade histology and add to the known rare cases of G34R-mutant tumors with gemistocytic morphology. Although K27M and G34R mutations are thought to be mutually exclusive, we document combined K27M and G34R mutations in HIST1H3B and present evidence suggesting the K27M-mutation drove tumor phenotype in this dual mutant glioma.

PHIP drives glioblastoma motility and invasion by regulating the focal adhesion complex.

The invasive behavior of glioblastoma is essential to its aggressive potential. Here, we show that pleckstrin homology domain interacting protein (PHIP), acting through effects on the force transduction layer of the focal adhesion complex, drives glioblastoma motility and invasion. Immunofluorescence analysis localized PHIP to the leading edge of glioblastoma cells, together with several focal adhesion proteins: vinculin (VCL), talin 1 (TLN1), integrin beta 1 (ITGB1), as well as phosphorylated forms of paxillin (pPXN) and focal adhesion kinase (pFAK). Confocal microscopy specifically localized PHIP to the force transduction layer, together with TLN1 and VCL. Immunoprecipitation revealed a physical interaction between PHIP and VCL. Targeted suppression of PHIP resulted in significant down-regulation of these focal adhesion proteins, along with zyxin (ZYX), and produced profoundly disorganized stress fibers. Live-cell imaging of glioblastoma cells overexpressing a ZYX-GFP construct demonstrated a role for PHIP in regulating focal adhesion dynamics. PHIP silencing significantly suppressed the migratory and invasive capacity of glioblastoma cells, partially restored following TLN1 or ZYX cDNA overexpression. PHIP knockdown produced substantial suppression of tumor growth upon intracranial implantation, as well as significantly reduced microvessel density and secreted VEGF levels. PHIP copy number was elevated in the classical glioblastoma subtype and correlated with elevated EGFR levels. These results demonstrate PHIP's role in regulating the actin cytoskeleton, focal adhesion dynamics, and tumor cell motility, and identify PHIP as a key driver of glioblastoma migration and invasion.

Transcriptomic evidence that von Economo neurons are regionally specialized extratelencephalic-projecting excitatory neurons.

von Economo neurons (VENs) are bipolar, spindle-shaped neurons restricted to layer 5 of human frontoinsula and anterior cingulate cortex that appear to be selectively vulnerable to neuropsychiatric and neurodegenerative diseases, although little is known about other VEN cellular phenotypes. Single nucleus RNA-sequencing of frontoinsula layer 5 identifies a transcriptomically-defined cell cluster that contained VENs, but also fork cells and a subset of pyramidal neurons. Cross-species alignment of this cell cluster with a well-annotated mouse classification shows strong homology to extratelencephalic (ET) excitatory neurons that project to subcerebral targets. This cluster also shows strong homology to a putative ET cluster in human temporal cortex, but with a strikingly specific regional signature. Together these results suggest that VENs are a regionally distinctive type of ET neuron. Additionally, we describe the first patch clamp recordings of VENs from neurosurgically-resected tissue that show distinctive intrinsic membrane properties relative to neighboring pyramidal neurons.

Cytomegalovirus is a tumor-associated virus: armed and dangerous.

Human cytomegalovirus (HCMV) gene products are present in multiple human malignancies, often in specific association with tumor cells and tumor vasculature. Emerging evidence from human and mouse models of CMV infection in cancer indicate that CMV can transform epithelial cells, promote epithelial to mesenchymal transition (EMT) and mesenchymal to epithelial (MET) in tumor cells, promote tumor angiogenesis and proliferation and incapacitate the host anti-CMV immune response. This review will discuss the increasing role of HCMV in human cancer by demonstrating how HCMV is well suited for impacting major themes in oncogenesis including initiation, promotion, progression, metastasis and immune evasion. What emerges is a picture of an extremely versatile pathogen that may play a significant role in human cancer progression and death.

Conserved cell types with divergent features in human versus mouse cortex.

Elucidating the cellular architecture of the human cerebral cortex is central to understanding our cognitive abilities and susceptibility to disease. Here we used single-nucleus RNA-sequencing analysis to perform a comprehensive study of cell types in the middle temporal gyrus of human cortex. We identified a highly diverse set of excitatory and inhibitory neuron types that are mostly sparse, with excitatory types being less layer-restricted than expected. Comparison to similar mouse cortex single-cell RNA-sequencing datasets revealed a surprisingly well-conserved cellular architecture that enables matching of homologous types and predictions of properties of human cell types. Despite this general conservation, we also found extensive differences between homologous human and mouse cell types, including marked alterations in proportions, laminar distributions, gene expression and morphology. These species-specific features emphasize the importance of directly studying human brain.

Publisher Correction: Toxoplasma Modulates Signature Pathways of Human Epilepsy, Neurodegeneration & Cancer.

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has been fixed in the paper.