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
. 2022 May 11;23(1):113.
doi: 10.1186/s13059-022-02677-z.

Refining colorectal cancer classification and clinical stratification through a single-cell atlas

Ateeq M Khaliq  1 , Cihat Erdogan  2 , Zeyneb Kurt  3 , Sultan Sevgi Turgut  4 , Miles W Grunvald  5 , Tim Rand  6 , Sonal Khare  6 , Jeffrey A Borgia  5 , Dana M Hayden  5 , Sam G Pappas  5 , Henry R Govekar  5 , Audrey E Kam  5 , Jochen Reiser  5 , Kiran Turaga  7 , Milan Radovich  1 , Yong Zang  1 , Yingjie Qiu  1 , Yunlong Liu  1 , Melissa L Fishel  1 , Anita Turk  1 , Vineet Gupta  5 , Ram Al-Sabti  5 , Janakiraman Subramanian  8 , Timothy M Kuzel  5 , Anguraj Sadanandam  9 , Levi Waldron  10 , Arif Hussain  11 , Mohammad Saleem  5 , Bassel El-Rayes  12 , Ameen A Salahudeen  6 , Ashiq Masood  13
Affiliations

Refining colorectal cancer classification and clinical stratification through a single-cell atlas

Ateeq M Khaliq et al. Genome Biol. .

Erratum in

  • Correction: Refining colorectal cancer classification and clinical stratification through a single-cell atlas.
    Khaliq AM, Erdogan C, Kurt Z, Turgut SS, Grunvald MW, Rand T, Khare S, Borgia JA, Hayden DM, Pappas SG, Govekar HR, Kam AE, Reiser J, Turaga K, Radovich M, Zang Y, Qiu Y, Liu Y, Fishel ML, Turk A, Gupta V, Al-Sabti R, Subramanian J, Kuzel TM, Sadanandam A, Waldron L, Hussain A, Saleem M, El-Rayes B, Salahudeen AA, Masood A. Khaliq AM, et al. Genome Biol. 2022 Jul 13;23(1):156. doi: 10.1186/s13059-022-02724-9. Genome Biol. 2022. PMID: 35831907 Free PMC article. No abstract available.

Abstract

Background: Colorectal cancer (CRC) consensus molecular subtypes (CMS) have different immunological, stromal cell, and clinicopathological characteristics. Single-cell characterization of CMS subtype tumor microenvironments is required to elucidate mechanisms of tumor and stroma cell contributions to pathogenesis which may advance subtype-specific therapeutic development. We interrogate racially diverse human CRC samples and analyze multiple independent external cohorts for a total of 487,829 single cells enabling high-resolution depiction of the cellular diversity and heterogeneity within the tumor and microenvironmental cells.

Results: Tumor cells recapitulate individual CMS subgroups yet exhibit significant intratumoral CMS heterogeneity. Both CMS1 microsatellite instability (MSI-H) CRCs and microsatellite stable (MSS) CRC demonstrate similar pathway activations at the tumor epithelial level. However, CD8+ cytotoxic T cell phenotype infiltration in MSI-H CRCs may explain why these tumors respond to immune checkpoint inhibitors. Cellular transcriptomic profiles in CRC exist in a tumor immune stromal continuum in contrast to discrete subtypes proposed by studies utilizing bulk transcriptomics. We note a dichotomy in tumor microenvironments across CMS subgroups exists by which patients with high cancer-associated fibroblasts (CAFs) and C1Q+TAM content exhibit poor outcomes, providing a higher level of personalization and precision than would distinct subtypes. Additionally, we discover CAF subtypes known to be associated with immunotherapy resistance.

Conclusions: Distinct CAFs and C1Q+ TAMs are sufficient to explain CMS predictive ability and a simpler signature based on these cellular phenotypes could stratify CRC patient prognosis with greater precision. Therapeutically targeting specific CAF subtypes and C1Q + TAMs may promote immunotherapy responses in CRC patients.

Keywords: CMS classification; Cancer-associated fibroblast; Colorectal cancer; Immunotherapy; Single-cell analysis; Stromal signatures.

PubMed Disclaimer

Conflict of interest statement

Ashiq Masood and Jeffrey A. Borgia received research funding from Tempus lab. Anjuraj Sadanandam receives research funding from Bristol-Myers Squibb; Merck KGaA, Pierre Fabre. Furthermore, Anjuraj Sadanandam. holds patent PCT/IB2013/060416, “Colorectal cancer classification with differential prognosis and personalized therapeutic responses,” and patent number 2011213.2 “Prognostic and Treatment Response Predictive Method.”

Figures

Fig. 1
Fig. 1
Identification and clustering of single cells. A Workflow of sample collection, sorting, and sequencing (methods contain full description for each step) and t-SNE characterization of the 49,859 cells profiled. B Identification of various cell types based on expression of specified marker genes. C Characterization of the proportion of cell types identified in each sample in tumor vs. normal colon tissue and Consensus Molecular Subtypes (CMS) of bulk RNA-seq data. D Characterization of the proportion of cell types identified in tumor vs. normal colon tissue, sidedness (right vs. left), microsatellite instability (MSI) status, single-cell Consensus Molecular Subtypes (scCMS) classification, Consensus Molecular Subtypes (CMS) of bulk RNA-seq data, and origin of sample. The graph represents total clusters and cell types identified after re-clustering of each cell compartment depicting global heterogeneous landscape of colorectal cancers
Fig. 2
Fig. 2
CMS Cell proportions, Gene set enrichment and trajectories of tumor cells. A Proportions of tumor cells classified as various CMS subtypes in various samples. Primary CRC datasets is labeled CAC and Lee et al. 2020 [18] dataset labeled as SMC. Annotation of tumor cells is based on bulk CMS classification. B Gene set variation expression analyses of combined primary CRC data and Lee et al. 2020 [18], in the tumor epithelial cell compartment. C Trajectory analysis of primary CRC dataset colored by bulk CMS status. D Trajectory analysis of Lee et al. 2020 [18] (Korean cohort) CRC dataset colored by bulk CMS status. E Trajectory analysis of Lee et al. 2020 [18] (Belgian cohort) colored by bulk CMS status
Fig. 3
Fig. 3
Fibroblast clusters in colon and colorectal tumors. A t-SNE of fibroblasts colored by Normal, CAF-S1 and CAF-S4 subtypes. B Dot plots showing the variable expression of fibroblast specific marker genes across CAF-S1 and CAF-S4. C Integration analysis of five CAFs subtypes from a breast cancer (BC) cohort to validate the existence of specific CAF subtypes in the CRC samples. D IHC representative images of CAF-S1 (FAP+, PDGFR-ß+) and CAF-S4 (RGS5+, MCAM+) in CRC sections from five independent patients. Asterisk (*) indicates tumor cells and arrows (>) indicate CAFs. Annotated by a board-certified GI pathologist. All images are × 20 magnification. High-resolution images are available on GitHub as source data [43]. E Boxplots show the distribution of cell types in two CRC bulk expression datasets, within tumors based on CMS status. The whiskers depict the 1.5 x IQR. The P-values for pairwise t-tests comparisons (with Benjamini-Hochberg correction) of cell abundance across CMS are shown in the figure. Note “NS”: P > 0.05, *P </=0.05, **P </=0.01, ***P </=0.001, ****P </=0.0001
Fig. 4
Fig. 4
Myeloid cell clusters in colon and colorectal tumors. A t-SNE of myeloid cells colored. by distinct cell types. B Heatmap showing the variable expression of myeloid specific marker genes across various myeloid cell types. C Identification of various myeloid cell subtypes based on expression of specified marker genes. D Boxplots show the distribution of cell types in two CRC bulk expression datasets, within tumors with varying CMS status [45, 48]. The whiskers depict the 1.5 x IQR. The P-values for pairwise t-tests comparisons (with Benjamini-Hochberg correction) of cell abundance across CMS are shown in the figure. E Box plot show TAMs infiltration and influence on CD8+ T cells in relation to CMS1 and CMS4. Note in GSE17536 Ythdf2 deficiency trended in CMS4 but did not reach the statistical significance. Note “NS”: P > 0.05, *P </=0.05, **P </=0.01, ***P </=0.001, ****P </=0.0001
Fig. 5
Fig. 5
T cell clusters in colon and colorectal tumors. A t-SNE of 22525 T cells colored by distinct clusters. B t-SNE plot showing the variable expression of T cell specific marker genes across various clusters. C Boxplots show the distribution of cell types in two CRC balk expression datasets, within tumors with varying CMS status. The whiskers depict the 1.5 x IQR. The P-values for pairwise t-tests comparisons (with Benjamini-Hochberg correction) of cell abundance across CMS are shown in the figure. Note “NS”: P > 0.05, *P </=0.05, **P </=0.01, ***P </=0.001, ****P </=0.0001
Fig. 6
Fig. 6
Continuous subtype scoring across cell type (GSE39582 [45], GSE1736 [48]). A, B Continuous scores reported by CMS classification across cell types show minimal separation in the top 2 principal components in GSE39582 and GSE17536 datasets respectively. All cell types are represented in CMS1-4 using PCSS1 and PCSS2 scores. Note that the cell types largely form a continuum along CMS status and are not clustered in discrete subtypes separate from one another. Cells and markers are colored by bulk CMS status accordingly to the tumor sample of origin. (PCSS1 = PC Cluster Subtype Scores1, PCSS2 = PC Cluster Subtype Scores1)
Fig. 7
Fig. 7
Survival analysis of two independent bulk dataset. The relationship between relative cell abundance and disease-free survival (DFS) in the GSE39582 [45] (A) and GSE17536 [48] (B) (COX regression analysis). Kaplan Meier curves depicting DFS in GSE39582 [45] (A) and GSE17536 [48] (B). Note, in addition to CMS4, CMS1-3 subgroups (good prognosis subtypes) with high CAF and C1Q+ TAMs signatures were associated with poor DFS. CAF's and C1Q+ TAMs stratified all CMS subgroups into high and low-poor survival subgroups beyond CMS categorization. HR and P values are indicated
Fig. 8
Fig. 8
CAF, TAM and CD8+ T-cell interactions in the CRC microenviroment. A Dotplot showing cell-cell interaction between CD8+ T cells and TAMs. B Circle plot directed ligand receptor interactions in CD8+ T cells and TAMs for better visualization. C Dotplot showing cell-cell interaction between CD8+ T cells and CAFs. D Circle plot directed ligand receptor interactions in CD8+ T cells and CAFs for better visualization. E Dotplot demonstrating cell-cell interaction between TAMs and CAFs. F Circle plot directed ligand receptor interactions in TAMs and CAFs for better visualization
See this image and copyright information in PMC

Similar articles

See all similar articles

Cited by

See all "Cited by" articles

References

    1. Arnold M, Sierra MS, Laversanne M, Soerjomataram I, Jemal A, Bray F. Global patterns and trends in colorectal cancer incidence and mortality. Gut. 2017;66:683–691. doi: 10.1136/gutjnl-2015-310912. - DOI - PubMed
    1. Molinari C, Marisi G, Passardi A, Matteucci L, De Maio G, Ulivi P. Heterogeneity in colorectal cancer: a challenge for personalized medicine? Int J Mol Sci. 2018;19 Available from: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC6321493/. Cited 2020 Dec 29. - PMC - PubMed
    1. Osterman E, Hammarström K, Imam I, Osterlund E, Sjöblom T, Glimelius B. Recurrence risk after radical colorectal cancer surgery—less than before, but how high is it? Cancers. 2020;12:3308. doi: 10.3390/cancers12113308. - DOI - PMC - PubMed
    1. André T, Shiu K-K, Kim TW, Jensen BV, Jensen LH, Punt C, et al. Pembrolizumab in Microsatellite-instability-high advanced colorectal cancer. N Engl J Med. 2020;383:2207–2218. doi: 10.1056/NEJMoa2017699. - DOI - PubMed
    1. Kopetz S, Grothey A, Yaeger R, Van Cutsem E, Desai J, Yoshino T, et al. Encorafenib, binimetinib, and cetuximab in BRAF V600E-mutated colorectal cancer. N Engl J Med. 2019;381:1632–1643. doi: 10.1056/NEJMoa1908075. - DOI - PubMed

Publication types

MeSH terms