Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2019 Aug 8;178(4):835-849.e21.
doi: 10.1016/j.cell.2019.06.024. Epub 2019 Jul 18.

An Integrative Model of Cellular States, Plasticity, and Genetics for Glioblastoma

Affiliations

An Integrative Model of Cellular States, Plasticity, and Genetics for Glioblastoma

Cyril Neftel et al. Cell. .

Abstract

Diverse genetic, epigenetic, and developmental programs drive glioblastoma, an incurable and poorly understood tumor, but their precise characterization remains challenging. Here, we use an integrative approach spanning single-cell RNA-sequencing of 28 tumors, bulk genetic and expression analysis of 401 specimens from the The Cancer Genome Atlas (TCGA), functional approaches, and single-cell lineage tracing to derive a unified model of cellular states and genetic diversity in glioblastoma. We find that malignant cells in glioblastoma exist in four main cellular states that recapitulate distinct neural cell types, are influenced by the tumor microenvironment, and exhibit plasticity. The relative frequency of cells in each state varies between glioblastoma samples and is influenced by copy number amplifications of the CDK4, EGFR, and PDGFRA loci and by mutations in the NF1 locus, which each favor a defined state. Our work provides a blueprint for glioblastoma, integrating the malignant cell programs, their plasticity, and their modulation by genetic drivers.

Keywords: CDK4; EGFR; NF1; PDGFRA; glioblastoma IDH-wildtype; glioblastoma stem cells; glioblastoma subtypes; lineage tracing; single-cell RNA-sequencing.

PubMed Disclaimer

Conflict of interest statement

DECLARATION OF INTERESTS

The authors declare no competing interests.

Figures

Figure 1.
Figure 1.. Classification of single cells from 28 glioblastomas.
(A) Inference of chromosomal CNAs based on average relative expression in windows of 100 analyzed genes. Row correspond to cells, with non-malignant (NM) cells that lack CNAs at the top followed by malignant cells (with CNAs, as defined in Figure S1) ordered by tumor and within a tumor clustered by overall CNA patterns. (B) tSNE plot of all single cells. Cells are colored based on presence of CNAs (blue) or high expression of sets of marker genes for macrophages (cyan), oligodendrocytes (magenta) or T-cells (green). (C) tSNE plot of all malignant cells, colored by tumor.
Figure 2.
Figure 2.. Expression signatures of intra-tumoral heterogeneity among malignant cells.
(A) Top: cell-to-cell correlation matrix of malignant cells from MGH105, with cells ordered by hierarchical clustering. Bottom: assignment of cells to potential overlapping clusters. (B) Hierarchical clustering of signatures for 269 potential clusters defined from 27 tumors. Groups of potential clusters are highlighted at the top and were used to define meta-modules. (C) Meta-modules, composed of genes consistently upregulated in potential clusters of the same group. Selected genes are indicated (see Table S2 for a full list). (D) Relative expression of meta-modules across neurodevelopment-related cell types as measured by scRNA-seq (Darmanis et al., 2017; Darmanis et al., 2015; Tirosh et al., 2016b). Error bars correspond to standard error.
Figure 3.
Figure 3.. Assignment of malignant cells to cellular states and their hybrids.
(A) Heatmap showing the meta-module scores of all non-cycling cells (left) and cycling cells (right). Within each group, the cells are ordered by their maximal score, for cells mapping to one meta-module, followed by cells mapping to two meta-modules (hybrid states, denoted as “H”). (B) Bar plot showing the percentage of cells with highest score for each meta-module. Adult and pediatric tumors are separated in order to demonstrate their distinct distributions. Error bars correspond to standard error, calculated by bootstrapping. (C) Bar plot showing the percentage of cycling cells among cells with highest score for each of the meta-modules. Error bars correspond to standard error, calculated by bootstrapping. (D) Bar plot showing the observed and expected percentages of hybrid cells (co-expressing two distinct meta-modules) out of all malignant cells. Expected percentages and their standard errors were calculated by shuffling the cell scores (STAR Methods). (E) In situ RNA hybridization of glioblastoma for NPC-like (CD24), MES-like (CD44) and proliferation (Ki67) markers. Arrows highlight representative cells positive for CD24 (blue) or CD44 (red). Arrowhead highlights a cell coexpressing CD24 and Ki67. (F) Two-dimensional representation of cellular states. Each quadrant corresponds to one cellular state, the exact position of malignant cells (dots) reflect their relative scores for the meta-modules, and their colors reflect the density of cycling cells (STAR Methods).
Figure 4.
Figure 4.. Intra-tumoral heterogeneity at the genetic and expression levels.
(A-B) Identification of genetic subclones by CNAs. Shown are the inferred CNAs of malignant cells in MGH125 (A) and MGH102 (B), separated into genetic subclones based on amplifications/deletions of specific chromosomes (STAR Methods). (C) Cell-state plots (as in Figure 3F) for six tumors with CNA-based subclones. Cells are colored by their subclone.
Figure 5.
Figure 5.. The distribution of glioblastoma cellular states is associated with chromosomal amplifications across the TCGA glioblastoma cohort.
(A) Pie charts displaying the fraction of cells in four cellular states in each glioblastoma from our cohort. Tumor indices are above each pie chart, with pediatric tumors indicated in red and recurrent tumors with “R”. Tumors are grouped by bulk TCGA subtype as labelled. (B) Analysis of the TCGA glioblastoma cohort shows that high-level amplifications of EGFR, PDGFRA and CDK4 are associated with high bulk scores for the AC-like, OPC-like, and NPC-like cellular states, respectively. Shown are significance values (-log10(P-value)) for the association of chromosomal amplifications with high bulk scores (shown above the zero line, indicating enrichment of the cellular state) or with low bulk scores (shown below the zero line, indicating depletion of the cellular state). Single chromosome gains are distinguished from high-level amplifications (Brennan et al., 2013) which are found to have significant associations with three cellular states (AC-like, OPC-like, NPC-like) while no associations of chromosomal amplifications were found with the MES-like state.
Figure 6.
Figure 6.. Glioblastoma oncogenes drive defined cellular states.
(A) Micrographs of immunofluorescence of mouse NPC over-expressing EGFR, CDK4 or eGFP immunostained for the astrocytic marker GFAP (red). (B) Quantification of GFAP+ cells shown in panel (A) (STAR Methods). (C) scRNA-seq scores for the AC-like signature (Y axis) of ranked cells (X axis) overexpressing EGFR (red) or GFP (black) (STAR Methods). (D) scRNA-seq scores for the NPC-like signature (Y axis) of ranked cells (X axis) overexpressing CDK4 (blue) or GFP (black). (E) Growth curve using NPCs over-expressing eGFP, EGFR or CDK4 shows increased proliferation (p<0.0001) in CDK4 expressing cells. RLU: Relative Light Units (arbitrary value). (F) Growth curve of astrocytes derived from the engineered NPCs (STAR Methods) shows significant (p<0.002, ANOVA) increase in growth of astrocytes overexpressing EGFR.
Figure 7.
Figure 7.. Cellular transitions in glioblastoma.
(A) Experimental workflow. Different fractions of cells were sorted from patient sample MGH143 and injected orthotopically into immunocompromised mice to generate PDXs. The patient sample and the PDXs subpopulations were subjected to scRNA-seq. (B) Samples described in (A) are each represented by a pie chart depicting the fraction of cells in four states. Pie charts are positioned on the X axis based on their sorted fraction and whether they represent injected or PDX sample, and on the Y axis based on their compositional similarity to the original patient sample (one minus the Manhattan distance over the fractions of four states). (C) Experimental workflow. Lentiviruses harboring oncogenes and unique barcodes were injected into the mouse hippocampus (STAR Methods) and the resulting tumors were analyzed by scRNA-seq. (D) Barcodes which were identified in multiple cells are each represented by a pie chart depicting the fraction of cells in each state. Pie charts are positioned based on the number of cells with the respective barcode (X axis), and the number of cellular states observed among these cells (Y axis). Pie chart sizes are proportional to log2 of the number of cells. (E) Experimental workflow. Primary cultures are established from glioblastoma samples (MGH143 and MGG23) and infected with lentiviruses harboring unique barcodes, xenografted into mouse brains and tumor formed are analyzed by scRNA-seq. (F) Unique barcodes from (E) are displayed as shown in (D). (G) Model for the cellular states of glioblastoma and their genetic and micro-environmental determinants. Mitotic spindles indicate cycling cells. Lighter/darker tones indicate strength of each program. Intermediate states are shown in between the four states and indicate transitions.

Comment in

Similar articles

Cited by

References

    1. Bao S, Wu Q, McLendon RE, Hao Y, Shi Q, Hjelmeland AB, Dewhirst MW, Bigner DD, and Rich JN (2006). Glioma stem cells promote radioresistance by preferential activation of the DNA damage response. Nature 444, 756–760. - PubMed
    1. Brennan CW, Verhaak RG, McKenna A, Campos B, Noushmehr H, Salama SR, Zheng S, Chakravarty D, Sanborn JZ, Berman SH, et al. (2013). The somatic genomic landscape of glioblastoma. Cell 155, 462–477. - PMC - PubMed
    1. Chen J, Li Y, Yu TS, McKay RM, Burns DK, Kernie SG, and Parada LF (2012). A restricted cell population propagates glioblastoma growth after chemotherapy. Nature 488, 522–526. - PMC - PubMed
    1. Darmanis S, Sloan SA, Croote D, Mignardi M, Chernikova S, Samghababi P, Zhang Y, Neff N, Kowarsky M, Caneda C, et al. (2017). Single-Cell RNA-Seq Analysis of Infiltrating Neoplastic Cells at the Migrating Front of Human Glioblastoma. Cell reports 21, 1399–1410. - PMC - PubMed
    1. Darmanis S, Sloan SA, Zhang Y, Enge M, Caneda C, Shuer LM, Gephart M.G. Hayden, Barres BA, and Quake SR (2015). A survey of human brain transcriptome diversity at the single cell level. Proceedings of the National Academy of Sciences of the United States of America - PMC - PubMed

Publication types

MeSH terms