.2014 Jun 20;344(6190):1396-401.

doi: 10.1126/science.1254257. Epub 2014 Jun 12.

Single-cell RNA-seq highlights intratumoral heterogeneity in primary glioblastoma

Anoop P Patel¹, Itay Tirosh², John J Trombetta², Alex K Shalek², Shawn M Gillespie³, Hiroaki Wakimoto⁴, Daniel P Cahill⁴, Brian V Nahed⁴, William T Curry⁴, Robert L Martuza⁴, David N Louis⁵, Orit Rozenblatt-Rosen², Mario L Suvà⁶, Aviv Regev⁷, Bradley E Bernstein⁸

Affiliations

¹ Department of Neurosurgery, Massachusetts General Hospital and Harvard Medical School, Boston, MA 02114, USA. Department of Pathology and Center for Cancer Research, Massachusetts General Hospital and Harvard Medical School, Boston, MA 02114, USA. Broad Institute of Harvard and Massachusetts Institute of Techonology (MIT), Cambridge, MA 02142, USA. Howard Hughes Medical Institute Chevy Chase, MD 20815, USA.
² Broad Institute of Harvard and Massachusetts Institute of Techonology (MIT), Cambridge, MA 02142, USA.
³ Department of Pathology and Center for Cancer Research, Massachusetts General Hospital and Harvard Medical School, Boston, MA 02114, USA. Broad Institute of Harvard and Massachusetts Institute of Techonology (MIT), Cambridge, MA 02142, USA. Howard Hughes Medical Institute Chevy Chase, MD 20815, USA.
⁴ Department of Neurosurgery, Massachusetts General Hospital and Harvard Medical School, Boston, MA 02114, USA.
⁵ Department of Pathology and Center for Cancer Research, Massachusetts General Hospital and Harvard Medical School, Boston, MA 02114, USA.
⁶ Department of Pathology and Center for Cancer Research, Massachusetts General Hospital and Harvard Medical School, Boston, MA 02114, USA. Broad Institute of Harvard and Massachusetts Institute of Techonology (MIT), Cambridge, MA 02142, USA. bernstein.bradley@mgh.harvard.edu aregev@broadinstitute.org suva.mario@mgh.harvard.edu.
⁷ Broad Institute of Harvard and Massachusetts Institute of Techonology (MIT), Cambridge, MA 02142, USA. Howard Hughes Medical Institute Chevy Chase, MD 20815, USA. Department of Biology, MIT, Cambridge, MA 02139, USA. bernstein.bradley@mgh.harvard.edu aregev@broadinstitute.org suva.mario@mgh.harvard.edu.
⁸ Department of Pathology and Center for Cancer Research, Massachusetts General Hospital and Harvard Medical School, Boston, MA 02114, USA. Broad Institute of Harvard and Massachusetts Institute of Techonology (MIT), Cambridge, MA 02142, USA. Howard Hughes Medical Institute Chevy Chase, MD 20815, USA. bernstein.bradley@mgh.harvard.edu aregev@broadinstitute.org suva.mario@mgh.harvard.edu.

PMID:24925914
PMCID: PMC4123637
DOI: 10.1126/science.1254257

Single-cell RNA-seq highlights intratumoral heterogeneity in primary glioblastoma

Anoop P Patel et al. Science.2014.

.2014 Jun 20;344(6190):1396-401.

doi: 10.1126/science.1254257. Epub 2014 Jun 12.

Authors

Affiliations

¹ Department of Neurosurgery, Massachusetts General Hospital and Harvard Medical School, Boston, MA 02114, USA. Department of Pathology and Center for Cancer Research, Massachusetts General Hospital and Harvard Medical School, Boston, MA 02114, USA. Broad Institute of Harvard and Massachusetts Institute of Techonology (MIT), Cambridge, MA 02142, USA. Howard Hughes Medical Institute Chevy Chase, MD 20815, USA.
² Broad Institute of Harvard and Massachusetts Institute of Techonology (MIT), Cambridge, MA 02142, USA.
³ Department of Pathology and Center for Cancer Research, Massachusetts General Hospital and Harvard Medical School, Boston, MA 02114, USA. Broad Institute of Harvard and Massachusetts Institute of Techonology (MIT), Cambridge, MA 02142, USA. Howard Hughes Medical Institute Chevy Chase, MD 20815, USA.
⁴ Department of Neurosurgery, Massachusetts General Hospital and Harvard Medical School, Boston, MA 02114, USA.
⁵ Department of Pathology and Center for Cancer Research, Massachusetts General Hospital and Harvard Medical School, Boston, MA 02114, USA.
⁶ Department of Pathology and Center for Cancer Research, Massachusetts General Hospital and Harvard Medical School, Boston, MA 02114, USA. Broad Institute of Harvard and Massachusetts Institute of Techonology (MIT), Cambridge, MA 02142, USA. bernstein.bradley@mgh.harvard.edu aregev@broadinstitute.org suva.mario@mgh.harvard.edu.
⁷ Broad Institute of Harvard and Massachusetts Institute of Techonology (MIT), Cambridge, MA 02142, USA. Howard Hughes Medical Institute Chevy Chase, MD 20815, USA. Department of Biology, MIT, Cambridge, MA 02139, USA. bernstein.bradley@mgh.harvard.edu aregev@broadinstitute.org suva.mario@mgh.harvard.edu.
⁸ Department of Pathology and Center for Cancer Research, Massachusetts General Hospital and Harvard Medical School, Boston, MA 02114, USA. Broad Institute of Harvard and Massachusetts Institute of Techonology (MIT), Cambridge, MA 02142, USA. Howard Hughes Medical Institute Chevy Chase, MD 20815, USA. bernstein.bradley@mgh.harvard.edu aregev@broadinstitute.org suva.mario@mgh.harvard.edu.

PMID:24925914
PMCID: PMC4123637
DOI: 10.1126/science.1254257

Abstract

Human cancers are complex ecosystems composed of cells with distinct phenotypes, genotypes, and epigenetic states, but current models do not adequately reflect tumor composition in patients. We used single-cell RNA sequencing (RNA-seq) to profile 430 cells from five primary glioblastomas, which we found to be inherently variable in their expression of diverse transcriptional programs related to oncogenic signaling, proliferation, complement/immune response, and hypoxia. We also observed a continuum of stemness-related expression states that enabled us to identify putative regulators of stemness in vivo. Finally, we show that established glioblastoma subtype classifiers are variably expressed across individual cells within a tumor and demonstrate the potential prognostic implications of such intratumoral heterogeneity. Thus, we reveal previously unappreciated heterogeneity in diverse regulatory programs central to glioblastoma biology, prognosis, and therapy.

PubMed Disclaimer

Figures

**Figure 1. Intratumoral glioblastoma heterogeneity quantified by single cell RNA-seq**
**(A)** Workflow depicts rapid dissociation and isolation of glioblastoma cells from primary tumors for generating single cell and bulk RNA-seq profiles and deriving glioblastoma culture models.(B) Clustering of copy number variation (CNV) profiles inferred from RNA-seq data for all single cells and a normal brain sample. Clusters (dendrogram) primarily reflect tumor-specific CNV (colored bar coded as in panel D). Topmost cluster (red, arrow) contains the normal brain sample and 10 single cells, nine of which correlate with normal oligodendrocyte expression profiles and one with normal monocytes (‘Oligo’ and ‘Mono’, black and white heatmap).(C) Heatmap of CNV signal normalized against the ‘normal’ cluster defined in (B) shows CNV changes by chromosome (columns) for individual cells (rows). All cells outside the normal cluster exhibit chromosome 7 gain (red) and chromosome 10 loss (blue), which are characteristic of glioblastoma.(D) Multidimensional scaling illustrates the relative similarity between all 430 single cells and population controls. The distance between any two cells reflects the similarity of their expression profiles. Cells group by tumor (color code), but each tumor also contains outliers that are more similar to cells in other tumors.(E) RNA-seq read densities (vertical scale of 10) over surface receptor genes are depicted for individual cells (rows) from MGH30. Cell-to-cell variability suggests a mosaic pattern of receptor expression, in contrast to constitutively expressedGAPDH.

**Figure 2. Unbiased analysis of intratumoral heterogeneity reveals coherent transcriptional modules**
**(A)** Gene sets that vary coherently between cells in specific tumors or across the global dataset (colored boxes) were identified by principal component analysis or clustering (24). Hierarchical clustering of these gene sets across all cells (tree) reveals four meta-signatures related to hypoxia, complement/immune response, oligodendrocytes and cell cycle.(B) Heatmap shows expression of the cell cycle meta-signature, selected cell cycle gene sets and representative genes from the signature (rows) in individual glioblastoma cells (columns). Cells were grouped by tumor and ordered by meta-signature score.(C) Heatmap depicts hypoxia meta-signature as in (B).

**Figure 3. Transcriptional signatures of a stem-like compartment in primary glioblastoma**
**(A)** Stem-like (GSC) and differentiated (DGC) culture models were derived from patient tumor MGH26. GSCs grow as spheres (left, top), initiate tumors in xenotransplantation (right, top), and express the stem cell marker CD133 (right, bottom).(B) Heatmap depicts expression of genes (rows) from a stemness signature in differentiated models (DGC, left columns), stem-like models (GSC, right columns) derived from 3 tumors, and in 70 individual cells from MGH31 (middle).(C) Bar plot depicts the Pearson correlation coefficient (y-axis) between the stemness signature and selected transcriptional modules in each tumor (colored bars): cell cycle, transcriptional targets of POU3F2, SOX2, SALL2, OLIG2 (core TF) (43), NFI transcriptional targets (NFI) (42), and the proneural (PN), classical (CL), mesenchymal (MES) and neural (N) subtypes defined by the Cancer Genome Atlas (21).(D) Plot depicts stemness score (y-axis) computed from stemness signature gene expression in individual cells from each tumor (xaxis) ordered by score. Bar plots depict the overall variance (y-axis, standard deviation) in the stemness score (red) and the average variance of simulated control gene sets (blue), confirming the significance of the gradient.

**Figure 4. Individual tumors contain a spectrum of glioblastoma subtypes and hybrid cellular states**
**(A)** Heatmap depicts average expression of classifier genes for each subtype (rows) across all classifiable cells grouped by tumor (columns). PN: proneural, CL: classical, MES: mesenchymal, N: neural. Each tumor contains a dominant subtype, but also has cells that conform to alternate subtypes.(B) Hexagonal plots depict bootstrapped classifier scores for all cells in each tumor. Each data point corresponds to a single cell and is positioned along three axes according to its relative scores for the indicated subtypes (Supplemental Methods). Cells corresponding to each subtype are indicated by solid color, while hybrid cells are depicted by two colors.(C) Kaplan-Meier survival curves are shown for proneural tumors from the Cancer Genome Atlas (21). Intratumoral heterogeneity was estimated based on detected signal for alternative subtypes, and used to partition the tumors into a pure proneural group and three groups with the indicated additional subtype. Tumors with mesenchymal signal had significantly worse outcome than pure proneural (p<0.05). (D) Kaplan-Meier survival curves shown for proneural tumors partitioned based on the relative strength of alternative subtype signatures in aggregate (24). Tumors with high signal for alternative subtypes had significantly worse outcome than pure proneural (p<0.05).

See this image and copyright information in PMC

Cited by

Single-cell and spatial transcriptome characterize coinhibitory cell-cell communications during histological progression of lung adenocarcinoma.
Luo J, Gao Q, Wang M, Liu H, Zhu H.Luo J, et al.Front Immunol. 2024 Oct 15;15:1430163. doi: 10.3389/fimmu.2024.1430163. eCollection 2024.Front Immunol. 2024.PMID:39474422Free PMC article.
Neuronal Activity Promotes Glioma Growth through Neuroligin-3 Secretion.
Venkatesh HS, Johung TB, Caretti V, Noll A, Tang Y, Nagaraja S, Gibson EM, Mount CW, Polepalli J, Mitra SS, Woo PJ, Malenka RC, Vogel H, Bredel M, Mallick P, Monje M.Venkatesh HS, et al.Cell. 2015 May 7;161(4):803-16. doi: 10.1016/j.cell.2015.04.012. Epub 2015 Apr 23.Cell. 2015.PMID:25913192Free PMC article.
Combined inhibition of Bcl-2/Bcl-xL and Usp9X/Bag3 overcomes apoptotic resistance in glioblastoma in vitro and in vivo.
Karpel-Massler G, Shu C, Chau L, Banu M, Halatsch ME, Westhoff MA, Ramirez Y, Ross AH, Bruce JN, Canoll P, Siegelin MD.Karpel-Massler G, et al.Oncotarget. 2015 Jun 10;6(16):14507-21. doi: 10.18632/oncotarget.3993.Oncotarget. 2015.PMID:26008975Free PMC article.
Localization of protoporphyrin IX during glioma-resection surgery via paired stimulated Raman histology and fluorescence microscopy.
Nasir-Moin M, Wadiura LI, Sacalean V, Juros D, Movahed-Ezazi M, Lock EK, Smith A, Lee M, Weiss H, Müther M, Alber D, Ratna S, Fang C, Suero-Molina E, Hellwig S, Stummer W, Rössler K, Hainfellner JA, Widhalm G, Kiesel B, Reichert D, Mischkulnig M, Jain R, Straehle J, Neidert N, Schnell O, Beck J, Trautman J, Pastore S, Pacione D, Placantonakis D, Oermann EK, Golfinos JG, Hollon TC, Snuderl M, Freudiger CW, Heiland DH, Orringer DA.Nasir-Moin M, et al.Nat Biomed Eng. 2024 Jun;8(6):672-688. doi: 10.1038/s41551-024-01217-3. Epub 2024 Jul 10.Nat Biomed Eng. 2024.PMID:38987630
Mosaicism for Receptor Tyrosine Kinase Activation in a Glioblastoma Involving Both PDGFRA Amplification and NTRK2 Fusion.
Shepherd DJ, Miller TE, Forst DA, Jones P, Nardi V, Martinez-Lage M, Stemmer-Rachamimov A, Gonzalez RG, Iafrate AJ, Ritterhouse LL.Shepherd DJ, et al.Oncologist. 2021 Nov;26(11):919-924. doi: 10.1002/onco.13835. Epub 2021 Jun 8.Oncologist. 2021.PMID:34041811Free PMC article.

See all "Cited by" articles

References

1. Yachida S, et al. Distant metastasis occurs late during the genetic evolution of pancreatic cancer. Nature. 2010;467:1114–1117. - PMC - PubMed
1. Eppert K, et al. Stem cell gene expression programs influence clinical outcome in human leukemia. Nat Med. 2011;17:1086–1093. - PubMed
1. Parsons DW, et al. An integrated genomic analysis of human glioblastoma multiforme. Science. 2008;321:1807–1812. - PMC - PubMed
1. Navin N, et al. Tumour evolution inferred by single-cell sequencing. Nature. 2011;472:90–94. - PMC - PubMed
1. Gerlinger M, et al. Intratumor heterogeneity and branched evolution revealed by multiregion sequencing. N Engl J Med. 2012;366:883–892. - PMC - PubMed

Publication types

Actions
Actions

MeSH terms

Actions
Actions
Actions
Actions
Actions
Actions
Actions
Actions
Actions
Actions
Actions
Actions
Actions

Substances

Actions

Associated data

Actions
- Search in PubMed
- Search in GEO

Related information

Grants and funding

LinkOut - more resources

Full Text Sources
Other Literature Sources
Medical
- MedlinePlus Health Information
Molecular Biology Databases
- NIAID Data Ecosystem - Find datasets on Infectious and Immune-mediated Diseases

Movatterモバイル変換

Account

Save citation to file

Email citation

Add to Collections

Add to My Bibliography

Your saved search

Create a file for external citation management software

Your RSS Feed

Full text links

Actions

Share

Single-cell RNA-seq highlights intratumoral heterogeneity in primary glioblastoma

Affiliations

Single-cell RNA-seq highlights intratumoral heterogeneity in primary glioblastoma

Authors

Affiliations

Abstract

Figures

Similar articles

Cited by

References

Publication types

MeSH terms

Substances

Associated data

Related information

Grants and funding

LinkOut - more resources

Full Text Sources

Other Literature Sources

Medical

Molecular Biology Databases