Abstract:Introduction: Colorectal cancer (CRC) is a highly heterogeneous cancer. The molecular and cellular characteristics differ between the colon and rectal cancer type due to the differences in their anatomical location and pathological properties. With the advent of single-cell sequencing, it has become possible to analyze inter- and intra-tumoral tissue heterogeneities. Methods: A comprehensive CRC immune atlas, comprising 62,398 immune cells, was re-structured into 33 immune cell clusters at the single-cell level. Further, the immune cell lineage heterogeneity of colon, rectal, and paracancerous tissues was explored. Simultaneously, we characterized the TAM phenotypes and analyzed the transcriptomic factor regulatory network of each macrophage subset using SCENIC. In addition, monocle2 was used to elucidate the B cell developmental trajectory. The crosstalk between immune cells was explored using CellChat and the patterns of incoming and outgoing signals within the overall immune cell population were identified. Afterwards, the bulk RNA-sequencing data from The Cancer Genome Atlas (TCGA) were combined and the relative infiltration abundance of the identified subpopulations was analyzed using CIBERSORT. Moreover, cell composition patterns could be classified into five tumor microenvironment (TME) subtypes by employing a consistent non-negative matrix algorithm. Finally, the co-expression and interaction between SPP1+TAMs and Treg cells in the tumor microenvironment were analyzed by multiplex immunohistochemistry. Results: In the T cell lineage, we found that CXCL13+T cells were more widely distributed in colorectal cancer tissues, and the proportion of infiltration was increased. In addition, Th17 was found accounted for the highest proportion in CD39+CD101+PD1+T cells. Mover, Ma1-SPP1 showed the characteristics of M2 phenotypes and displayed an increased proportion in tumor tissues, which may promote angiogenesis. Plasma cells (PCs) displayed a significantly heterogeneous distribution in tumor as well as normal tissues. Specifically, the IgA+ PC population could be shown to be decreased in colorectal tumor tissues whereas the IgG+ PC one was enriched. In addition, information flow mediated by SPP1 and CD44, regulate signaling pathways of tumor progression. Among the five TME subtypes, the TME-1 subtype displayed a markedly reduced proportion of T-cell infiltration with the highest proportion of macrophages which was correlated to the worst prognosis. Finally, the co-expression and interaction between SPP1+TAMs and Treg cells were observed in the CD44 enriched region. Discussion: The heterogeneity distribution and phenotype of immune cells were analyzed in colon cancer and rectal cancer at the single-cell level. Further, the prognostic role of major tumor-infiltrating lymphocytes and TME subtypes in CRC was evaluated by integrating bulk RNA. These findings provide novel insight into the immunotherapy of CRC.

Cross‐modal integration of bulk RNA‐seq and single‐cell RNA sequencing data to reveal T‐cell exhaustion in colorectal cancer

Elucidating the role of T-cell exhaustion-related genes in colorectal cancer: a single-cell bioinformatics perspective

A Multiomics Analysis-Assisted Deep Learning Model Identifies a Macrophage-Oriented Module As a Potential Therapeutic Target in Colorectal Cancer.

Deep learning using bulk RNA-seq data expands cell landscape identification in tumor microenvironment

Single-cell transcriptome analysis reveals heterogeneity and convergence of the tumor microenvironment in colorectal cancer

Deep immunophenotyping reveals clinically distinct cellular states and ecosystems in large-scale colorectal cancer

Integration of single-cell RNA sequencing and bulk RNA transcriptome sequencing reveals a heterogeneous immune landscape and pivotal cell subpopulations associated with colorectal cancer prognosis

Disease associated human TCR characterization by deep-learning framework TCR-DeepInsight

Inference of Developmental Hierarchy and Functional States of Exhausted T Cells from Epigenetic Profiles with Deep Learning

Deep Single-Cell RNA Sequencing Data of Individual T Cells from Treatment-Naïve Colorectal Cancer Patients

Single‐cell analysis revealed a potential role of T‐cell exhaustion in colorectal cancer with liver metastasis

Single cell transcriptome of CD8+ T cells in multiple cancers reveals comprehensive exhaustion associated mechanisms

Integrating machine learning and single-cell trajectories to analyze T-cell exhaustion to predict prognosis and immunotherapy in colon cancer patients

Navigating colorectal cancer prognosis: A Treg‐related signature discovered through single‐cell and bulk transcriptomic approaches

Digitaldlsorter: Deep-Learning on scRNA-Seq to Deconvolute Gene Expression Data

Integrating spatial and single-cell transcriptomics reveals tumor heterogeneity and intercellular networks in colorectal cancer

Multi-omics Blueprint of Cellular Senescence in Deciphering Immune Characteristics and Prognosis Stratification of CRC

Integrated single-cell and bulk RNA-seq analysis identifies a prognostic T-cell signature in colorectal cancer

Integrated single-cell and bulk RNA-seq for constructing a prognostic risk model of immune-related genes in colorectal cancer

Integration of single-cell sequencing and bulk RNA-seq to identify and develop a prognostic signature related to colorectal cancer stem cells

Integrative analysis of single-cell and bulk RNA seq to reveal the prognostic model and tumor microenvironment remodeling mechanisms of cuproptosis-related genes in colorectal cancer