Integrated electrophysiological and genomic profiles of single cells reveal spiking tumor cells in human glioma.

Prior studies have described the complex interplay that exists between glioma cells and neurons; however, the electrophysiological properties endogenous to glioma cells remain obscure. To address this, we employed Patch-sequencing (Patch-seq) on human glioma specimens and found that one-third of patched cells in IDH mutant (IDHmut) tumors demonstrate properties of both neurons and glia. To define these hybrid cells (HCs), which fire single, short action potentials, and discern if they are of tumoral origin, we developed the single cell rule association mining (SCRAM) computational tool to annotate each cell individually. SCRAM revealed that HCs possess select features of GABAergic neurons and oligodendrocyte precursor cells, and include both tumor and non-tumor cells. These studies characterize the combined electrophysiological and molecular properties of human glioma cells and describe a cell type in human glioma with unique electrophysiological and transcriptomic properties that may also exist in the non-tumor brain.

Glioma epileptiform activity and progression are driven by IGSF3-mediated potassium dysregulation.

Seizures are a frequent pathophysiological feature of malignant glioma. Recent studies implicate peritumoral synaptic dysregulation as a driver of brain hyperactivity and tumor progression; however, the molecular mechanisms that govern these phenomena remain elusive. Using scRNA-seq and intraoperative patient ECoG recordings, we show that tumors from seizure patients are enriched for gene signatures regulating synapse formation. Employing a human-to-mouse in vivo functionalization pipeline to screen these genes, we identify IGSF3 as a mediator of glioma progression and dysregulated neural circuitry that manifests as spreading depolarization (SD). Mechanistically, we discover that IGSF3 interacts with Kir4.1 to suppress potassium buffering and found that seizure patients exhibit reduced expression of potassium handlers in proliferating tumor cells. In vivo imaging reveals that dysregulated synaptic activity emanates from the tumor-neuron interface, which we confirm in patients. Our studies reveal that tumor progression and seizures are enabled by ion dyshomeostasis and identify SD as a driver of disease.

Single-cell RNA sequencing reveals immunosuppressive myeloid cell diversity during malignant progression in a murine model of glioma.

Recent studies have shown the importance of the dynamic tumor microenvironment (TME) in high-grade gliomas (HGGs). In particular, myeloid cells are known to mediate immunosuppression in glioma; however, it is still unclear if myeloid cells play a role in low-grade glioma (LGG) malignant progression. Here, we investigate the cellular heterogeneity of the TME using single-cell RNA sequencing in a murine glioma model that recapitulates the malignant progression of LGG to HGG. LGGs show increased infiltrating CD4+ and CD8+ T cells and natural killer (NK) cells in the TME, whereas HGGs abrogate this infiltration. Our study identifies distinct macrophage clusters in the TME that show an immune-activated phenotype in LGG but then evolve to an immunosuppressive state in HGG. We identify CD74 and macrophage migration inhibition factor (MIF) as potential targets for these distinct macrophage populations. Targeting these intra-tumoral macrophages in the LGG stage may attenuate their immunosuppressive properties and impair malignant progression.

CaSpER identifies and visualizes CNV events by integrative analysis of single-cell or bulk RNA-sequencing data.

RNA sequencing experiments generate large amounts of information about expression levels of genes. Although they are mainly used for quantifying expression levels, they contain much more biologically important information such as copy number variants (CNVs). Here, we present CaSpER, a signal processing approach for identification, visualization, and integrative analysis of focal and large-scale CNV events in multiscale resolution using either bulk or single-cell RNA sequencing data. CaSpER integrates the multiscale smoothing of expression signal and allelic shift signals for CNV calling. The allelic shift signal measures the loss-of-heterozygosity (LOH) which is valuable for CNV identification. CaSpER employs an efficient methodology for the generation of a genome-wide B-allele frequency (BAF) signal profile from the reads and utilizes it for correction of CNVs calls. CaSpER increases the utility of RNA-sequencing datasets and complements other tools for complete characterization and visualization of the genomic and transcriptomic landscape of single cell and bulk RNA sequencing data.

Somatic POLE mutations cause an ultramutated giant cell high-grade glioma subtype with better prognosis.

BACKGROUND: Malignant high-grade gliomas (HGGs), including the most aggressive form, glioblastoma multiforme, show significant clinical and genomic heterogeneity. Despite recent advances, the overall survival of HGGs and their response to treatment remain poor. In order to gain further insight into disease pathophysiology by correlating genomic landscape with clinical behavior, thereby identifying distinct HGG molecular subgroups associated with improved prognosis, we performed a comprehensive genomic analysis. METHODS: We analyzed and compared 720 exome-sequenced gliomas (136 from Yale, 584 from The Cancer Genome Atlas) based on their genomic, histological, and clinical features. RESULTS: We identified a subgroup of HGGs (6 total, 4 adults and 2 children) that harbored a statistically significantly increased number of somatic mutations (mean = 9257.3 vs 76.2, P = .002). All of these "ultramutated" tumors harbored somatic mutations in the exonuclease domain of the polymerase epsilon gene (POLE), displaying a distinctive genetic profile, characterized by genomic stability and increased C-to-A transversions. Histologically, they all harbored multinucleated giant or bizarre cells, some with predominant infiltrating immune cells. One adult and both pediatric patients carried homozygous germline mutations in the mutS homolog 6 (MSH6) gene. In adults, POLE mutations were observed in patients younger than 40 years and were associated with a longer progression-free survival. CONCLUSIONS: We identified a genomically, histologically, and clinically distinct subgroup of HGGs that harbored somatic POLE mutations and carried an improved prognosis. Identification of distinctive molecular and pathological HGG phenotypes has implications not only for improved classification but also for potential targeted treatments.