Glioblastoma (GBM) is the most frequent and malignant form of primary brain tumor. GBM is essentially incurable and its resistance to therapy is attributed to a subpopulation of cells called glioma stem cells (GSCs). To meet the present shortage of relevant GBM cell (GC) lines we developed a library of annotated and validated cell lines derived from surgical samples of GBM patients, maintained under conditions to preserve GSC characteristics. This collection, which we call the Human Glioblastoma Cell Culture (HGCC) resource, consists of a biobank of 48 GC lines and an associated database containing high-resolution molecular data. We demonstrate that the HGCC lines are tumorigenic, harbor genomic lesions characteristic of GBMs, and represent all four transcriptional subtypes. The HGCC panel provides an open resource for in vitro and in vivo modeling of a large part of GBM diversity useful to both basic and translational GBM research.

Intratumoral heterogeneity is a hallmark of glioblastoma multiforme and thought to negatively affect treatment efficacy. Here, we establish libraries of glioma-initiating cell (GIC) clones from patient samples and find extensive molecular and phenotypic variability among clones, including a range of responses to radiation and drugs. This widespread variability was observed as a continuum of multitherapy resistance phenotypes linked to a proneural-mesenchymal shift in the transcriptome. Multitherapy resistance was associated with a semi-stable cell state that was characterized by an altered DNA methylation pattern at promoter regions of mesenchymal master regulators and enhancers. The gradient of cell states within the GIC compartment constitutes a distinct form of heterogeneity. Our findings may open an avenue toward the development of new therapeutic rationales designed to reverse resistant cell states.

SummarySignificant advances have been made in methods to analyze genomes and transcriptomes of single cells, but to fully define cell states, proteins must also be accessed as central actors defining a cell’s phenotype. Methods currently used to analyze endogenous protein expression in single cells are limited in specificity, throughput, or multiplex capability. Here, we present an approach to simultaneously and specifically interrogate large sets of protein and RNA targets in lysates from individual cells, enabling investigations of cell functions and responses. We applied our method to investigate the effects of BMP4, an experimental therapeutic agent, on early-passage glioblastoma cell cultures. We uncovered significant heterogeneity in responses to treatment at levels of RNA and protein, with a subset of cells reacting in a distinct manner to BMP4. Moreover, we found overall poor correlation between protein and RNA at the level of single cells, with proteins more accurately defining responses to treatment.