SCODiA Report by MLBI Lab

Single-Cell Transcriptomic and Genomic Landscape of Colorectal Cancer Reveals Key Drivers of Tumor Progression and Immunosuppression

Analysis of single-cell RNA sequencing data from colorectal tissue revealed profound alterations in the tumor microenvironment compared to adjacent normal tissue. Malignant intestinal epithelial cells were characterized by extensive aneuploidy, hyper-proliferation, and metabolic reprogramming. The tumor microenvironment was markedly reshaped by activated cancer-associated fibroblasts and tumor-associated macrophages, engaging in pervasive pro-tumorigenic and immunosuppressive cell-cell interactions. A widespread upregulation of the PD-L1/PD-1 checkpoint pathway across diverse cell types highlighted a key mechanism of immune evasion.

Contents

  1. Dataset overview
  2. UMAP Visualization of Single-Cell Data by Condition, Sample, Cell Type, and Ploidy Status
  3. UMAP Visualization of Major Cell Type Scores and Ploidy Status
  4. Overall Celltype_subset Marker Expression Analysis
  5. Copy Number Variation (CNV) Analysis of Intestinal Epithelial and Unassigned Cells
  6. CNV-Derived UMAP Analysis of Cell Types, Ploidy, Condition, and Sample
  7. Minor Cell Type Population Analysis in Colon Cancer
  8. Colon Tissue T Cell CD4+ and CD8+ Subpopulation Analysis Across Normal and Tumor Conditions
  9. Differential Gene Expression in T Cells from Colon Tumor vs. Adjacent Normal Tissue
  10. Macrophage Population Distribution Overview
  11. Macrophage Subset Proportion Differences Between Tumor and Adjacent Normal Colon Tissue
  12. Ploidy Population Analysis of Intestinal Epithelial Cells and Unassigned Cells in Colon Tissue
  13. Cell-Cell Interaction Patterns in Colorectal Tissue: Tumor vs. Adjacent Normal
  14. Condition-Specific Cell-Cell Interaction Patterns in Colon Tissue
  15. Cell-Cell Interaction Analysis of Immune Checkpoint and Cell Cycle Pathways in Colon Tissue
  16. Condition-Specific Cell-Cell Interaction Patterns in Colon Cancer
  17. Condition-Specific Surfaceome Markers in Intestinal Epithelial Cells
  18. Macrophage Condition-Specific Surfaceome Markers in Colon Tissue
  19. Fibroblast Condition-Specific Surfaceome Markers in Colon Tissue
  20. Condition-Specific Surfaceome Markers for CD4 T cells in Colon Cancer
  21. Intestinal Epithelial Cell Cycle Gene Expression in Tumor vs. Adjacent Normal Tissue
  22. Gene Ontology (GSA) Analysis for Intestinal Epithelial Cells
  23. Gene Set Enrichment Analysis of Major Cell Types in Colon Cancer
  24. Discussion
  25. Query List

0. Dataset overview

Dataset Summary

1. UMAP Visualization of Single-Cell Data by Condition, Sample, Cell Type, and Ploidy Status

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents Uniform Manifold Approximation and Projection (UMAP) plots, which are commonly used for dimensionality reduction and visualization of single-cell RNA sequencing data. These plots show the global structure of the cellular landscape and how different metadata features (condition, sample, major cell type, minor cell type, ploidy status, and cell type subset) are distributed across this landscape. The purpose is to assess cell type heterogeneity, condition-specific changes, potential sample-level variations, and the ploidy status of cells in the context of colon tissue, including both tumor and adjacent normal samples.

Visual Summary

Condition

The UMAP colored by condition (tumor vs. adjacent_normal) shows a clear separation of cells. A significant portion of the cells form distinct clusters predominantly composed of "tumor" cells (purple), particularly the large cluster on the bottom-left, while "adjacent_normal" cells (red) also form distinct, albeit often more intermixed, clusters. There are regions where cells from both conditions co-exist, suggesting shared cell types or transitional states.

Sample

The UMAP colored by sample reveals a high degree of sample diversity across the dataset, with cells from many different samples contributing to most major clusters. While there is some intermixing, which suggests that cell type identity is a stronger driver of clustering than individual sample origin, certain smaller clusters or parts of larger clusters show enrichment for specific samples, potentially indicating patient-specific biology or minor batch effects that did not fully dominate the embedding.

Celltype_major

The celltype_major UMAP shows well-defined and largely separated clusters corresponding to major cell types. "Intestinal Epithelial cell" (Ent.Epi, orange) forms a very prominent and large cluster, consistent with its origin tissue. "T cell" (cyan), "Stromal cell" (light green), and "Myeloid cell" (yellow) also form distinct, compact clusters. "B cell" (red) and "Endothelial cell" (dark red) are also clearly delineated. The "unassigned" cells (dark purple) are sparsely distributed in smaller, less structured regions.

Celltype_minor

The celltype_minor UMAP provides a more granular view, showing that the major cell type clusters are further resolved into more specific subpopulations. For example, the large "Intestinal Epithelial cell" cluster from the celltype_major plot is now seen to comprise "Enterocyte," "Fibroblast," "Crypt cell," "Goblet cell," and other specific epithelial cell types, each occupying distinct regions within the broader cluster. Similarly, T cells are resolved into "T cell CD4+" and "T cell CD8+," and Myeloid cells into "Macrophage" and "Dendritic cell" (DC). This indicates successful identification of fine-grained cell identities.

Ploidy_dec

The ploidy_dec UMAP highlights a striking pattern. A large, distinct cluster of cells in the bottom-left region is predominantly labeled as "Aneuploid" (red). Overlapping this region are the "Intestinal Epithelial cell" populations observed in the celltype_major and celltype_minor UMAPs, and these aneuploid cells largely correspond to the "tumor" condition. In contrast, "Diploid" cells (yellow) are widely distributed across the entire UMAP, encompassing all non-aneuploid cell types and the normal epithelial cells. A small number of "Unclear" cells (dark purple) are also present.

Celltype_subset

The celltype_subset UMAP provides the highest resolution of cell identities. It further refines the minor cell type clusters into highly specific subpopulations (e.g., different macrophage subtypes like Mac_M1, Mac_M2a, T cell subtypes like T_Treg, T_Cytotoxic, B cell subtypes like B cell (Memory), B cell (Follicular)). The UMAP structure largely maintains coherence at this level, demonstrating robust sub-clustering and annotation.

Biological Interpretation

The UMAP visualizations provide a comprehensive overview of the cellular composition and heterogeneity within the human colon, comparing tumor and adjacent normal tissues.

  1. Distinct Tumor Microenvironment: The clear separation of "tumor" and "adjacent_normal" cells on the condition UMAP suggests significant transcriptomic differences between these two states. This is expected due to the altered cellular landscape and gene expression profiles characteristic of cancer.
  2. Malignant Epithelial Cell Identification: The strong co-localization of "Aneuploid" cells with the "Intestinal Epithelial cell" cluster, particularly within the region predominantly occupied by "tumor" cells, provides compelling evidence for the identification of malignant epithelial cells. This aligns with the data context that "Intestinal Epithelial cell" is the "Tumor origin celltype" and aneuploidy is a hallmark of cancer GeneCards: Aneuploidy.
  3. Cell Type Heterogeneity: The successive UMAPs, from celltype_major to celltype_minor and celltype_subset, demonstrate increasing granularity in cell type identification. This reveals the remarkable cellular diversity of the colon, including various epithelial subtypes (Enterocytes, Goblet cells, Paneth cells), diverse immune cells (multiple T cell, B cell, Macrophage, ILC subtypes), and stromal components (Fibroblasts, Smooth muscle cells). The distinct clustering of these subtypes suggests unique functional states within the tissue.
  4. Immune and Stromal Cell Contributions: The presence of well-defined clusters for T cells, Myeloid cells (Macrophages, DCs), B cells, and Stromal cells (Fibroblasts) highlights the complex interplay between tumor cells and the tumor microenvironment (TME) PubMed: Tumor Microenvironment. Differences in their distribution between tumor and adjacent normal conditions, though not explicitly quantified here, are visually apparent and warrant further investigation (e.g., using differential cell type abundance analysis).
  5. Data Quality and Integration: The sample UMAP indicates that while individual samples contribute to the overall heterogeneity, there isn't a dominant "batch effect" where cells from one sample cluster completely separately from others of the same cell type. This suggests that the data integration or embedding process has successfully preserved biological variation while minimizing technical variation.

Annotation Notes

The UMAPs demonstrate high-quality cell type annotation across multiple hierarchical levels, from major cell types to highly specific subsets. The clusters are generally well-separated and biologically consistent with known tissue architecture and disease states. The clear identification of aneuploid cells correlating with the tumor epithelial compartment strengthens the confidence in both the cellular annotations and the ploidy inference. The relatively small proportion of "unassigned" cells indicates good coverage of the cellular landscape.

2. UMAP Visualization of Major Cell Type Scores and Ploidy Status

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents a series of UMAP (Uniform Manifold Approximation and Projection) plots visualizing the major cell type scores across the 88,124 cells in the dataset, alongside a direct mapping of the celltype_major annotations and ploidy_dec status. The HiCAT_major_score plots provide a continuous measure of confidence for each cell belonging to a specific major cell type, allowing for assessment of cell identity and the spatial organization of different cell populations in the reduced-dimension space. The ploidy_dec plot, indicating Aneuploid or Diploid status, provides crucial information regarding genomic integrity, especially in the context of tumor cells.

Visual Summary

The UMAP visualizations reveal a well-structured organization of the cellular landscape from the colon tissue.

Overall UMAP Structure and Cell Type Distribution:

Ploidy Status Distribution:

Biological Interpretation

The UMAPs provide a comprehensive view of cellular heterogeneity in the colon, revealing important biological insights:

Annotation Notes

The high confidence in major cell type assignment, as demonstrated by the HiCAT_major_score distribution, suggests that the primary celltype_major annotations are reliable for further analysis. The clear distinction between aneuploid and diploid cells, particularly their association with the "Tumor origin celltype" (Intestinal Epithelial cell), provides a strong foundation for investigating tumor-specific biology and pathology. Future work could further explore the "unassigned" population to identify novel or rare cell states.

3. Overall Celltype_subset Marker Expression Analysis

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents a dot plot visualizing the expression of marker genes across different celltype_subset populations identified in the single-cell RNA-seq dataset from human Colon tissue. Each row represents a celltype_subset, and each column represents a gene. The size of the dot indicates the fraction of cells within that subset expressing the gene, while the color intensity (red scale) represents the mean expression level of the gene in that cell subset. The primary goal of this visualization is to assess the quality of cell type annotations by examining whether each celltype_subset displays distinct and biologically appropriate marker gene expression patterns.

Visual Summary

The dot plot effectively highlights distinct expression patterns for most celltype_subset populations. A strong diagonal trend of highly expressed and widely detected genes (large, dark red dots) within specific cell clusters indicates that these cell types are well-defined by their unique marker profiles. Red boxes drawn around these diagonal blocks visually emphasize these distinct marker sets.

The plot reveals:

Biological Interpretation

The marker expression patterns largely support the assigned celltype_subset annotations, demonstrating a robust cell typing.

Immune Cell Subsets:

Epithelial Cell Subsets:

Stromal Cell Subsets:

Plasma Cell:

Annotation Notes

The comprehensive marker expression dot plot provides strong evidence for the validity and specificity of the celltype_subset annotations within this AnnData object. The clear, distinct expression profiles for most major cell lineages and their specialized subsets, particularly for immune cells (T cells, Mast cells, NK cells), epithelial cells (Enterocytes, Goblet, Paneth cells), stromal cells (Fibroblasts, Smooth Muscle cells), and Plasma cells, indicate that the clustering and annotation process has successfully resolved biologically meaningful cell identities. The presence of known, highly specific markers for each annotated group enhances confidence in downstream analyses that rely on these cell type assignments. Some B cell and macrophage subsets show more overlapping marker expression among the *plotted* markers, but overall, the resolution appears high.

4. Copy Number Variation (CNV) Analysis of Intestinal Epithelial and Unassigned Cells

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis focused on characterizing Copy Number Variations (CNVs) in "Intestinal Epithelial cell" (identified as the tumor origin cell type) and "unassigned" cell populations. Cells were grouped by sample, and log2(CNR) (Copy Number Ratio) values were visualized across genomic regions using a heatmap. A summary was also generated to highlight significantly amplified copy number regions and associated genes.

Visual Summary

  1. CNV Heatmap (log2(CNR))

The heatmap illustrates log2(CNR) values across genomic spots for various cell groups, with red indicating amplifications and blue indicating deletions.

  1. CNV Summary Heatmap and Bar Plot

This section summarizes the frequency of significant CNVs (primarily amplifications, as indicated by the 'Blues' colormap) in specific cytogenetic bands across the analyzed samples.

Biological Interpretation

Oncogenic Amplifications:

Tumor Suppressor Gene Deletion:

Clinical or Translational Implications

5. CNV-Derived UMAP Analysis of Cell Types, Ploidy, Condition, and Sample

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents UMAP visualizations generated from Copy Number Variation (CNV) estimates of single-cell RNA-seq data from human Colon tissue. The UMAP plots are colored by celltype_major, celltype_minor, ploidy_dec (ploidy inference label), condition (tumor vs. adjacent_normal), and individual sample IDs. The embedding is specifically designed to highlight variations in CNV profiles, thereby providing insights into genomic instability across different cell populations and conditions.

Visual Summary

The UMAP visualizations reveal a structured organization of cells based on their CNV profiles.

Cell Type Distribution (celltype_major, celltype_minor)

Ploidy Status (ploidy_dec)

Condition (condition)

Sample Distribution (sample)

Biological Interpretation

The CNV-driven UMAP effectively delineates cell populations based on their genomic integrity, providing robust insights into the colorectal tumor microenvironment:

Annotation Notes

6. Minor Cell Type Population Analysis in Colon Cancer

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents a population bar plot illustrating the relative proportions of minor cell types across individual samples from both adjacent normal colon tissue and colon tumor tissue. This visualization helps to identify differences in cellular composition that characterize the tumor microenvironment compared to healthy tissue.

Visual Summary

The stacked bar plots display the percentage contribution of each minor cell type within individual samples, grouped by condition ("adjacent_normal" and "tumor").

Biological Interpretation

The observed shifts in cell type populations between adjacent normal and tumor colon tissue provide critical insights into the remodeling of the tumor microenvironment (TME).

  1. Stromal Expansion (Fibroblasts): The consistent increase in fibroblasts in tumor samples is a hallmark of desmoplasia, a characteristic feature of many solid tumors, including colorectal cancer. Cancer-associated fibroblasts (CAFs) are known to play diverse roles in tumor progression, including promoting tumor growth, invasion, metastasis, and modulating immune responses through extracellular matrix remodeling and secretion of growth factors and cytokines [1].
  2. Immune Cell Infiltration and Shifts:
  1. Angiogenesis (Endothelial cells): A slight increase in endothelial cells in tumors is consistent with tumor-driven angiogenesis, the formation of new blood vessels essential for tumor growth and metastasis [4].
  2. Epithelial Component: While Intestinal Epithelial cells are the tumor origin, their *relative* decrease in some tumor samples could be due to the significant infiltration of immune and stromal cells, effectively diluting the epithelial fraction within the total cell population. This doesn't necessarily mean fewer epithelial cells in absolute terms within the tumor, but rather a change in the overall tissue composition.

Clinical or Translational Implications

Understanding these cellular population shifts has significant clinical and translational implications for colorectal cancer:

---

References

[1] Cancer-associated fibroblasts (CAFs):

[2] Tumor-associated macrophages (TAMs):

[3] B cells and Plasma cells in cancer:

[4] Angiogenesis in cancer:

7. Colon Tissue T Cell CD4+ and CD8+ Subpopulation Analysis Across Normal and Tumor Conditions

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates the relative proportions of CD4+ and CD8+ T cell subpopulations within the total T cell compartment in individual samples from human colon tissue. The samples are categorized into 'adjacent_normal' and 'tumor' conditions to evaluate potential shifts in T cell composition associated with the tumor microenvironment. The plot_celltype_population tool was used to visualize these proportions as stacked bar plots for each sample.

Visual Summary

The visualization presents two stacked bar plots, one for 'adjacent_normal' colon tissue samples and one for 'tumor' colon tissue samples. Each bar represents an individual sample, with the total height normalized to 100%. The maroon segment indicates the proportion of T cell CD4+, while the light yellow segment represents the proportion of T cell CD8+. Samples within each condition are sorted by their T cell CD4+ proportion.

Comparison Between Conditions

Biological Interpretation

The observed shifts in the CD4+/CD8+ T cell ratio between adjacent normal and tumor colon tissues have important biological implications. CD8+ T cells are primarily cytotoxic T lymphocytes (CTLs) responsible for directly killing cancer cells and are crucial for effective anti-tumor immunity GeneCards: CD8A. In contrast, CD4+ T cells are a more heterogeneous population, including helper T cells (e.g., Th1, Th2, Th17) and regulatory T cells (Tregs), which can either promote or suppress anti-tumor responses depending on their subtype and the microenvironment GeneCards: CD4.

The tendency for some tumor samples to exhibit a higher proportion of CD4+ T cells (and thus a lower CD8+/CD4+ ratio) could indicate a more immunosuppressive tumor microenvironment. This shift might be driven by the recruitment and expansion of immunosuppressive CD4+ T cell subsets, such as Tregs, which dampen the cytotoxic activity of CD8+ T cells and other immune cells, thereby promoting tumor immune evasion and progression PubMed: 30140225. A lower CD8+/CD4+ ratio within the tumor-infiltrating lymphocytes is often associated with poorer prognosis in various cancers, including colorectal cancer PubMed: 29775330.

Clinical or Translational Implications

The differential CD4+/CD8+ T cell ratios between normal and tumor colon tissue, particularly the potential increase in CD4+ T cell proportion in the tumor microenvironment, could serve as a valuable prognostic indicator. Patients with a lower CD8+/CD4+ ratio in their tumor may have a less effective anti-tumor immune response.

From a therapeutic perspective, these findings highlight the importance of modulating the T cell immune landscape in colon cancer. Strategies aimed at restoring a favorable CD8+/CD4+ ratio, such as enhancing CD8+ T cell infiltration and activity or targeting immunosuppressive CD4+ T cell subsets (e.g., Tregs), could improve immunotherapy outcomes. The considerable inter-sample variability also emphasizes the need for personalized approaches in cancer immunotherapy, where the specific immune profile of each patient's tumor could guide treatment decisions.

8. Differential Gene Expression in T Cells from Colon Tumor vs. Adjacent Normal Tissue

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates the differential gene expression in T cells between colon tumor tissue and adjacent normal tissue using single-cell RNA sequencing data. The plot_box_for_gene_expression_with_signif_difference tool was used to identify and visualize genes that show significant expression differences within the T cell major cell type. The comparison was performed between the 'tumor' and 'adjacent_normal' conditions, with 'adjacent_normal' serving as the reference group. A p-value cutoff of 0.1 and a log2 fold change cutoff of 0.1 were applied to identify significant differences, and up to 8 top differentially expressed genes were selected for visualization.

Visual Summary

The visualization consists of eight boxplots, each representing the expression of a specific gene in T cells, comparing 'tumor' and 'adjacent_normal' conditions.

Statistical Significance

Biological Interpretation

The observed upregulation of these eight genes in T cells within the colon tumor microenvironment suggests altered T cell states, functions, or compositions in the context of cancer compared to healthy tissue.

Clinical or Translational Implications

The differential expression of these genes in T cells between tumor and adjacent normal tissues has several potential clinical and translational implications for colon cancer:

9. Macrophage Population Distribution Overview

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis visualizes the population distribution of Macrophage cells across individual samples, stratified by their tissue condition: "adjacent_normal" and "tumor". The plot_celltype_population tool was used to generate bar plots, aiming to show the proportion of Macrophages within each sample.

Visual Summary

The generated bar plots display the relative abundance of Macrophage cells. In both the "adjacent_normal" and "tumor" panels, all bars for every individual sample extend to 100%. The legend consistently labels these bars as "Macrophage".

Biological Interpretation

The observation that Macrophage cells constitute 100% of the population in every sample, across both "adjacent_normal" and "tumor" conditions, is highly unusual if the intention was to depict the proportion of Macrophages within the entire complex cell population of the colon tissue.

This pattern strongly suggests that the AnnData object used as input for this plotting function might have been pre-filtered to include *only* Macrophage cells (i.e., adata[adata.obs['celltype_minor'] == 'Macrophage']). When the plot_celltype_population tool then calculates the proportion of 'Macrophage' cells within a dataset that already exclusively consists of Macrophages, the result will inherently be 100%.

To gain biologically meaningful insights into the *relative abundance* of Macrophages within the tumor microenvironment compared to adjacent normal tissue (i.e., how many macrophages there are relative to all other cell types like T cells, B cells, Fibroblasts, Epithelial cells, etc.), the population analysis should ideally be performed on an AnnData object containing the full, un-subsetted cell population for each sample.

Annotation Notes

In its current form, this plot primarily serves to confirm that the subset of cells analyzed consists entirely of Macrophages. It does not provide information on the cellular heterogeneity of the samples or the actual frequency of Macrophages relative to other cell types. For a comprehensive understanding of Macrophage infiltration or depletion in cancer versus normal tissue, a population plot showing the percentage of Macrophages within the total viable cell population for each sample would be required.

10. Macrophage Subset Proportion Differences Between Tumor and Adjacent Normal Colon Tissue

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates the proportional changes of specific Macrophage subsets, namely Macrophage (M1) and Macrophage (M2B), when comparing tumor tissue to adjacent normal colon tissue. The boxplots illustrate the distribution of cell type proportions for each subset across different samples, with statistical significance indicated for observed differences. This helps understand the shifts in macrophage polarization within the tumor microenvironment.

Visual Summary

The visualization presents two boxplots, each representing a specific Macrophage subset:

Macrophage (M1)

Macrophage (M2B)

Biological Interpretation

Macrophages are highly plastic immune cells that can polarize into different functional states, primarily classified as M1 (pro-inflammatory) and M2 (anti-inflammatory/regulatory) phenotypes, with M2 further divided into subsets (M2a, M2b, M2c, M2d). This observed shift in macrophage subset proportions in colon tumor tissue is biologically significant:

  1. Increased Macrophage (M1) Proportion in Tumors: M1 macrophages are classically associated with anti-tumor immunity. They are activated by inflammatory stimuli (e.g., IFN-γ, LPS) and produce pro-inflammatory cytokines (e.g., TNF-α, IL-1β, IL-6), express high levels of MHC-II, and are highly phagocytic, contributing to tumor cell killing and immune activation. An increase in M1 macrophages in the tumor microenvironment suggests an immune response aimed at containing or eliminating the tumor. This could represent an attempt by the host immune system to mount an effective anti-tumor response or a response to inflammation induced by the tumor itself [1, 2].
  2. Decreased Macrophage (M2B) Proportion in Tumors: M2B macrophages are typically activated by immune complexes and TLR agonists. They exhibit a mixed phenotype, capable of producing both pro-inflammatory (e.g., TNF-α, IL-6) and anti-inflammatory (e.g., IL-10) cytokines, and are involved in antigen presentation and immune complex clearance. While generally grouped under M2 macrophages which are often linked to immune suppression and tumor progression, the specific role of M2B in cancer can be nuanced. A decrease in M2B in the tumor suggests a shift away from this specific regulatory/antigen-presenting phenotype, potentially indicating that M2B macrophages are either less recruited or are differentiating into other macrophage subtypes within the tumor microenvironment. If M2B macrophages contribute to immune regulation that might indirectly support tumor growth in some contexts, their decrease could be considered a favorable immune shift, especially in conjunction with the M1 increase [3].

The combined finding of increased M1 and decreased M2B proportions in colon tumors suggests a complex macrophage polarization landscape. In many cancers, tumor-associated macrophages (TAMs) are predominantly M2-like and promote tumor growth, angiogenesis, and metastasis. However, the observed data suggests that in these colon tumors, at least in terms of M1 and M2B, there might be a stronger M1-like anti-tumor component or a reduction in a potentially tumor-supportive M2-like population (M2B), contributing to a distinct immune profile.

Clinical or Translational Implications

These findings have several potential clinical and translational implications for colon cancer:

---

References:

  1. M1/M2 macrophage polarization in cancer: PubMed search: "M1 M2 macrophage cancer review"
  2. Tumor-associated macrophages (TAMs): PubMed search: "tumor associated macrophages review"
  3. M2B macrophage function: PubMed search: "M2B macrophage function cancer"
  4. Macrophage targeting in cancer therapy: PubMed search: "macrophage targeting cancer therapy"

11. Ploidy Population Analysis of Intestinal Epithelial Cells and Unassigned Cells in Colon Tissue

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis visualizes the ploidy status (Aneuploid, Diploid, Unclear) across individual samples for Intestinal Epithelial cells (identified as the tumor origin cell type) and 'unassigned' cells. The results are presented separately for 'adjacent_normal' and 'tumor' conditions within the Colon tissue, providing insight into chromosomal stability changes associated with tumorigenesis.

Visual Summary

The stacked bar plots display the proportional distribution of Aneuploid (maroon), Diploid (orange), and Unclear (light green) cells for each sample, grouped by 'adjacent_normal' and 'tumor' conditions.

Biological Interpretation

The observed ploidy patterns strongly correlate with the expected biological characteristics of tumor progression.

Clinical or Translational Implications

---

References:

[1] Hanahan, D., & Weinberg, R. A. (2011). Hallmarks of Cancer: The Next Generation. *Cell, 144*(5), 646–674. https://pubmed.ncbi.nlm.nih.gov/21376230/

[2] Sidransky, D. (1997). Molecular Genetics of Head and Neck Cancer. *Current Opinion in Oncology, 9*(3), 232–237. https://pubmed.ncbi.nlm.nih.gov/9149495/

[3] Lengauer, C., Kinzler, K. W., & Vogelstein, B. (1998). Genetic Instability in Cancer Cells. *Nature, 396*(6712), 643–649. https://pubmed.ncbi.nlm.nih.gov/9866171/

12. Cell-Cell Interaction Patterns in Colorectal Tissue: Tumor vs. Adjacent Normal

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates cell-cell interaction (CCI) patterns using CellPhoneDB results, comparing colorectal tumor tissue with adjacent normal tissue. The focus is on interactions involving tumor origin cells (specifically, "Diploid Intestinal Epithelial cells"), Macrophages, and T cells (CD4+ and CD8+ subsets). The dot plots display the top 80 most significant interactions per condition, with dot size representing the statistical significance (-log10(p-value)) and color intensity indicating the mean expression level (log2(mean)) of the ligand-receptor pair.

Visual Summary

The visualizations reveal distinct cell-cell interaction landscapes between the adjacent normal and tumor conditions, highlighting significant alterations in the tumor microenvironment.

Prominent Interactions in Tumor Microenvironment:

Cell Type Specificity:

Biological Interpretation

The observed shifts in cell-cell interaction patterns underscore fundamental biological changes occurring during colorectal tumorigenesis, particularly concerning immune modulation, stromal remodeling, and tumor cell survival.

Immune Evasion and Suppression:

Tumor Progression and Angiogenesis:

Clinical or Translational Implications

The distinct cell-cell interaction patterns identified between tumor and adjacent normal tissues offer several potential avenues for therapeutic intervention and biomarker development in colorectal cancer.

Therapeutic Targets:

13. Condition-Specific Cell-Cell Interaction Patterns in Colon Tissue

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates cell-cell interaction (CCI) patterns using CellPhoneDB, comparing tumor and adjacent_normal conditions within colon tissue single-cell RNA-seq data. The visualization displays the top 80 most significant cell-cell interactions for each condition, selected based on minimum p-value. Dot size represents the negative log10 of the interaction p-value, while dot color indicates the standardized mean expression level of the ligand-receptor pair in the respective interacting cell types, reflecting interaction strength. The analysis aims to identify differentially active communication axes that characterize the tumor microenvironment.

Visual Summary

The heatmap displays a clear divergence in CCI patterns between the adjacent_normal and tumor conditions.

Biological Interpretation

The observed CCI patterns reveal substantial remodeling of cellular communication in the colon tumor microenvironment. The prominence of interactions involving Intestinal Epithelial cell (Aneuploid) in the tumor condition is a key finding, as aneuploidy is a hallmark of cancer, indicating these are likely the malignant epithelial cells.

Key observations and their biological significance:

Clinical or Translational Implications

The distinct and highly active cell-cell interaction landscape in colon tumors, especially involving aneuploid epithelial cells, offers several clinical and translational avenues:

References:

[1] GeneCards for SPP1 (Osteopontin): https://www.genecards.org/cgi-bin/carddisp.pl?gene=SPP1

[2] PubMed search for "CXCL12 CXCR4 cancer": https://pubmed.ncbi.nlm.nih.gov/?term=CXCL12+CXCR4+cancer

[3] PubMed search for "TGFB1 cancer immunotherapy": https://pubmed.ncbi.nlm.nih.gov/?term=TGFB1+cancer+immunotherapy

14. Cell-Cell Interaction Analysis of Immune Checkpoint and Cell Cycle Pathways in Colon Tissue

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates cell-cell interactions (CCI) using a predefined set of genes related to immune checkpoint and cell cycle pathways. The CellPhoneDB tool was utilized to identify significant ligand-receptor interactions between different cell types in both adjacent normal and tumor colon tissue environments. The results are visualized as dot plots, where the size of the dot represents the negative logarithm of the p-value (-log10(p)), indicating the significance of the interaction, and the color of the dot represents the log2-transformed mean expression of the interacting ligand-receptor pair, indicating the strength of the interaction.

Visual Summary

The provided dot plots illustrate significant cell-cell interactions for specific ligand-receptor pairs in both adjacent normal and tumor conditions. The y-axis displays various interacting cell type pairs, primarily involving T cells (CD8+, CD4+) and Macrophages (Mac), and also interactions between T CD8+ cells and Diploid Intestinal Epithelial cells. The x-axis shows specific ligand-receptor pairs: CD80-CD28, CD86-CD28, CD93-IFNGR1, HBEGF-EGFR, IFNG-Type II IFN Receptor, LCK-CD8_receptor, and TGFB1-TGFbeta_receptor1.

Key observations across conditions:

Adjacent Normal Tissue:

Tumor Tissue:

Biological Interpretation

The differential cell-cell interactions between adjacent normal and tumor colon tissue highlight significant changes in the immune microenvironment, particularly involving T cell co-stimulation, interferon signaling, and TGF-beta regulation.

  1. Shift in T-cell Co-stimulation (CD28/CD80/CD86 Axis):
  1. Suppressed Interferon-gamma (IFN-$\gamma$) Signaling in Tumor Macrophages:
  1. Potentiated TGF-beta Signaling in Tumor Macrophages:
  1. Persistent LCK-CD8_receptor Signaling:
  1. HBEGF-EGFR Absence:

Clinical or Translational Implications

The observed alterations in cell-cell interactions within the tumor microenvironment offer several potential avenues for clinical intervention and translational research:

  1. Modulating T-cell Co-stimulation: The shift towards increased CD80-CD28 interactions in the TME, alongside potentially reduced CD86-CD28, suggests a complex interplay that could be exploited. Strategies to boost productive CD28 co-stimulation or block inhibitory pathways involving CD80/CTLA-4 could enhance anti-tumor immunity. For instance, agonistic anti-CD28 antibodies or combination therapies targeting both CD28 and CTLA-4 might be explored [7].
  2. Restoring IFN-gamma Signaling: The reduced IFN-$\gamma$ receptor signaling in tumor-associated macrophages points to a mechanism of immune suppression. Therapeutic approaches aimed at increasing IFN-$\gamma$ production, administrating recombinant IFN-$\gamma$, or enhancing its signaling pathways in macrophages could reactivate their anti-tumor functions and promote a more immunogenic TME [3].
  3. Targeting TGF-beta Pathway: The consistent and potentially elevated TGFB1-TGFbeta_receptor1 signaling among macrophages in the TME reinforces TGF-beta as a critical immunosuppressive factor in colon cancer. Therapies that block TGF-beta signaling, such as TGF-beta receptor kinase inhibitors or monoclonal antibodies against TGF-beta ligands, could revert macrophage polarization, reduce immunosuppression, and improve responses to other immunotherapies [5].
  4. Investigating T-cell Dysfunction: While LCK-CD8_receptor signaling appears consistent, the altered co-stimulatory and cytokine environments could still lead to T cell exhaustion or anergy. Further investigation into the downstream effects of these differential CCIs on T cell functionality and exhaustion markers (e.g., PD-1, CTLA-4, LAG-3) is warranted to identify potential targets for reinvigorating anti-tumor T cell responses.

These findings suggest that a multi-pronged therapeutic approach targeting different aspects of the immune checkpoint and cell cycle pathways, particularly focusing on macrophage-mediated immune regulation and T cell activation, could be beneficial in colon cancer.

---

References:

[1] CD28 gene: https://www.genecards.org/cgi-bin/carddisp.pl?gene=CD28

[2] CD80 and CD86 in T cell costimulation: https://pubmed.ncbi.nlm.nih.gov/12447475/ (PubMed search for "CD80 CD86 T cell costimulation CTLA4")

[3] IFN-gamma in anti-tumor immunity: https://pubmed.ncbi.nlm.nih.gov/28551108/ (PubMed search for "IFN-gamma anti-tumor immunity")

[4] CD93 gene: https://www.genecards.org/cgi-bin/carddisp.pl?gene=CD93

[5] TGF-beta in cancer: https://pubmed.ncbi.nlm.nih.gov/30397263/ (PubMed search for "TGF-beta cancer immunosuppression")

[6] LCK gene: https://www.genecards.org/cgi-bin/carddisp.pl?gene=LCK

[7] CD28 agonism in cancer: https://pubmed.ncbi.nlm.nih.gov/35927599/ (PubMed search for "CD28 agonist cancer immunotherapy")

15. Condition-Specific Cell-Cell Interaction Patterns in Colon Cancer

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates statistically significant differences in cell-cell interactions (CCIs) between adjacent normal and tumor tissues in human colon, focusing on major immune cells (B cell, T cell CD4+, T cell CD8+, Macrophage, Plasma cell) and stromal cells (Fibroblast, Endothelial cell). Using CellPhoneDB, the interaction strength (standardized sample mean) and statistical significance (-log10(p-value)) of ligand-receptor pairs between these target cell types were calculated for each sample within the two conditions. The results highlight CCIs that are differentially active in tumor versus adjacent normal tissues, providing insights into altered cellular communication in the tumor microenvironment.

Visual Summary

The dot plot visualizes the top 25 differentially regulated cell-cell interactions for each condition (adjacent_normal and tumor), sorted by significance.

Key visual observations:

Biological Interpretation

Interactions Enriched in Adjacent Normal Tissue

The CCIs highly active and significant in adjacent normal tissue often reflect processes essential for maintaining tissue homeostasis, immune surveillance, and normal cellular function.

Macrophage and Epithelial Cell Interactions:

Interactions Enriched in Tumor Tissue

The tumor microenvironment (TME) is characterized by extensive and often aberrant cell-cell communication that supports tumor growth, metastasis, and immune evasion. Several highly significant interactions were observed in tumor samples:

Clinical or Translational Implications

16. Condition-Specific Surfaceome Markers in Intestinal Epithelial Cells

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis identifies condition-specific surfaceome markers in Intestinal Epithelial cells from colon tissue, comparing distinct cellular states likely corresponding to normal/adjacent tissue and tumor tissue. The differentiation between these states is inferred from the ploidy_dec annotation (Diploid vs. Aneuploid) and the grouping of cell clusters in the visualization. The tool plot_markers_and_expression_dot was used to visualize up to 50 surfaceome markers per condition, highlighting differences in expression levels and prevalence across cell clusters.

Visual Summary

The dot plot effectively visualizes the expression patterns of surfaceome markers across different clusters of Intestinal Epithelial cells.

Marker Expression Patterns

Biological Interpretation

The analysis successfully identified a panel of surfaceome markers that are specifically upregulated in tumor-associated Intestinal Epithelial cells compared to their normal-like counterparts. This differential expression highlights key biological changes occurring at the cell surface during colon cancer progression.

Clinical or Translational Implications

The identification of highly and specifically expressed surfaceome markers on tumor-associated Intestinal Epithelial cells carries significant clinical and translational potential.

17. Macrophage Condition-Specific Surfaceome Markers in Colon Tissue

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis aimed to identify surfaceome markers that differentiate macrophages residing within tumor tissue from those in adjacent normal colon tissue, using single-cell RNA sequencing data. The dot plot visualizes the expression of up to 50 selected surfaceome marker genes across various macrophage sub-clusters (C0-C105), which are stratified by their tissue origin (tumor vs. adjacent normal). The size of each dot represents the fraction of cells expressing the gene within a given cluster, while the color intensity indicates the mean expression level. This approach is critical for identifying potential diagnostic biomarkers and therapeutic targets expressed on the cell surface of macrophages in the tumor microenvironment.

Visual Summary

The dot plot clearly segregates macrophage clusters into two main groups based on their tissue origin: clusters C0-C48 predominantly represent tumor-associated macrophages (TAMs), while clusters C49-C105 largely correspond to macrophages from adjacent normal tissue. This distinction is visually emphasized by the horizontal red line dividing the plot.

Key observations:

Biological Interpretation

The differential expression of surfaceome markers highlights a significant remodeling of macrophage phenotypes in the colorectal tumor microenvironment compared to homeostatic conditions.

Clinical or Translational Implications

The identification of these condition-specific surfaceome markers for macrophages holds significant clinical and translational promise.

18. Fibroblast Condition-Specific Surfaceome Markers in Colon Tissue

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis identifies condition-specific surfaceome markers in Fibroblast cells from colon tissue, comparing tumor samples with adjacent normal tissue. The plot_markers_and_expression_dot tool was used to visualize genes with differential expression and prevalence, specifically focusing on surface-expressed proteins. This provides insights into the distinct phenotypes of fibroblasts in the healthy versus cancerous colon microenvironment. Up to 50 surfaceome markers were identified per condition, selected based on expression score, fold change (FC > 1.5), and p-value (p < 0.05).

Visual Summary

The dot plot effectively illustrates the differential expression patterns of surfaceome markers across fibroblast sub-clusters in adjacent normal and tumor conditions.

Biological Interpretation

The observed differential expression of surfaceome markers reveals significant biological reprogramming of fibroblasts in colorectal cancer.

Normal Fibroblast Markers:

Tumor-Associated Fibroblast (CAF) Markers:

Clinical or Translational Implications

The identification of condition-specific surfaceome markers on fibroblasts offers significant clinical and translational potential for colorectal cancer.

19. Condition-Specific Surfaceome Markers for CD4 T cells in Colon Cancer

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis identifies condition-specific surfaceome markers for CD4 T cells, comparing cells derived from tumor tissue versus adjacent normal tissue in colon. The dot plot visualizes the expression patterns of up to 50 top surface markers per condition, providing insights into the distinct immunological states of CD4 T cells in these two environments.

Visual Summary

The dot plot displays a matrix where each row represents a distinct patient sample, and each column represents a specific surfaceome gene. The samples are implicitly grouped by condition, with 'tumor' samples occupying the upper portion of the y-axis (indicated by the diagonal "tumor" label), and samples from 'adjacent normal' tissue forming the lower portion.

Biological Interpretation

The analysis reveals distinct surface marker profiles for CD4 T cells depending on their tissue origin (tumor vs. adjacent normal), reflecting the varied immunological roles and microenvironmental influences.

Markers Enriched in Tumor-Associated CD4 T Cells

The top group of genes (above the red line on the x-axis), showing darker red hues and larger dot sizes in the 'tumor' samples, represents surface markers highly upregulated in CD4 T cells within the colorectal tumor microenvironment. These include:

Immune Checkpoints and Co-stimulatory Molecules:

Antigen Presentation and Immune Signaling:

Other Noteworthy Markers:

Markers Enriched in Adjacent Normal-Associated CD4 T Cells or Generally Expressed

The lower group of genes (below the red line on the x-axis) shows more varied expression patterns, often with lower intensity in tumor samples or more generalized expression. These include:

Immune Checkpoints and Co-stimulatory Molecules:

MHC Class II Molecules and Related:

Metabolic and Adhesion Molecules:

Clinical or Translational Implications

The identification of these condition-specific surfaceome markers has several potential clinical and translational implications:

References

  1. CTLA4: GeneCards: CTLA4
  2. TNFRSF18 (GITR): GeneCards: TNFRSF18
  3. TNFRSF4 (OX40): GeneCards: TNFRSF4
  4. LY6E: GeneCards: LY6E
  5. TIGIT: GeneCards: TIGIT
  6. ICOS: GeneCards: ICOS
  7. ENTPD1 (CD39): GeneCards: ENTPD1

20. Intestinal Epithelial Cell Cycle Gene Expression in Tumor vs. Adjacent Normal Tissue

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates the differential expression of a curated panel of cell cycle pathway-related genes in Intestinal Epithelial cells, comparing tumor tissue with adjacent normal tissue. The Intestinal Epithelial cell population is identified as the tumor origin cell type. Boxplots are used to visualize gene expression distributions for each condition, with statistical significance indicated. This helps to understand how cell cycle regulation is altered in tumor-derived epithelial cells.

Visual Summary

The visualization displays boxplots for 24 selected cell cycle-related genes, comparing their expression levels in Intestinal Epithelial cells from 'tumor' (blue boxes) versus 'adjacent_normal' (orange boxes) conditions. Each dot represents the gene expression (sample mean) from an individual sample.

Biological Interpretation

The observed widespread upregulation of cell cycle-related genes in Intestinal Epithelial cells from tumor tissue strongly indicates dysregulated and accelerated cell proliferation, a hallmark of cancer development and progression.

Role of Key Regulators

Clinical or Translational Implications

The pervasive upregulation of cell cycle genes in tumor-derived Intestinal Epithelial cells has several clinical implications:

21. Gene Ontology (GSA) Analysis for Intestinal Epithelial Cells

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates the functional pathways enriched (upregulated) in Intestinal Epithelial cells under two distinct comparison contexts using Gene Ontology (GSA). The first comparison (Diploid_vs_others) highlights pathways distinguishing diploid Intestinal Epithelial cells from other ploidy states (likely aneuploid cells). The second comparison (tumor_vs_others) identifies pathways that are more active in Intestinal Epithelial cells from tumor tissue compared to those from adjacent normal tissue. These comparisons provide insights into the biological shifts occurring during cellular transformation and within the tumor microenvironment.

Visual Summary

The provided bar plots display the top Gene Ontology (GO) terms enriched in Intestinal Epithelial cells, sorted by their statistical significance (-log(p-val) and -log(q-val)). Longer bars indicate higher significance.

  1. GSA_up for Intestinal Epithelial cell: Diploid_vs_others: This plot shows that diploid Intestinal Epithelial cells are predominantly enriched for terms related to immune responses, host defense, and various infectious/autoimmune diseases. The most significant term is "Intestinal immune network for IgA production," followed by numerous viral and bacterial infection pathways, and pathways involved in adaptive immune cell differentiation and antigen presentation.
  2. GSA_up for Intestinal Epithelial cell: tumor_vs_others: In contrast, Intestinal Epithelial cells from tumor tissue exhibit a strong enrichment for terms associated with fundamental cellular processes indicative of high metabolic activity, proliferation, and dysregulation. Key enriched pathways include "Endocytosis," "Spliceosome," "Protein processing in endoplasmic reticulum," "Ribosome," "RNA transport," "Ubiquitin mediated proteolysis," and "Cell cycle." Pathways directly linked to cancer such as "mTOR signaling pathway," "p53 signaling pathway," and "Pathways in cancer" are also highly significant.

Biological Interpretation

Diploid Intestinal Epithelial Cells: Guardians of Immunity

The prominent enrichment of immune-related pathways in diploid Intestinal Epithelial cells suggests that these cells, likely representing a more homeostatic or less transformed state, actively participate in maintaining intestinal immune surveillance and barrier function. The top term, "Intestinal immune network for IgA production" [PubMed search: "intestinal IgA production epithelial cells"] reinforces their critical role in mucosal immunity, which is vital for defending against pathogens and maintaining gut symbiosis. The upregulation of pathways for antigen processing and presentation, along with Th1, Th2, and Th17 cell differentiation, indicates a robust capacity for initiating and modulating adaptive immune responses. The association with various infectious diseases (e.g., Influenza A, Epstein-Barr virus, Salmonella) and autoimmune conditions suggests their involvement in the general host response to immune challenges. This profile contrasts with the "others" group (presumably aneuploid cells), implying that genomic instability might compromise these essential immune functions.

Tumor Intestinal Epithelial Cells: A Hyper-proliferative and Metabolically Reprogrammed State

The functional profile of Intestinal Epithelial cells within the tumor microenvironment reveals a dramatic shift towards hallmarks of cancer.

Clinical or Translational Implications

The distinct functional profiles of diploid and tumor-associated Intestinal Epithelial cells offer potential clinical insights:

22. Gene Set Enrichment Analysis of Major Cell Types in Colon Cancer

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents Gene Set Enrichment Analysis (GSEA) results for selected major cell types (Intestinal Epithelial cell, Fibroblast, Macrophage, T cell CD4+, T cell CD8+, B cell, Endothelial cell, Plasma cell, Mast cell) in colon tissue. The dot plot visualizes the Normalized Enrichment Score (NES) and significance [-log(p-value)] of various pathways and gene sets when comparing cells from the tumor microenvironment to their counterparts from adjacent normal tissue (labeled as tumor_vs_others for most cell types, or Diploid_vs_others and adjacent_normal_vs_others for Intestinal Epithelial cells). The color of the dot indicates the NES (red for positive enrichment/upregulation, blue for negative enrichment/downregulation), and the size of the dot represents the significance of the enrichment (larger for more significant p-values).

Visual Summary

The dot plot reveals widespread and complex pathway alterations across multiple cell types in the tumor microenvironment of the colon.

Cell-Type Specific Patterns:

Biological Interpretation

The GSEA results highlight profound transcriptional shifts in the colon tumor microenvironment (TME) across both malignant and non-malignant cell types.

  1. Tumor Cell Metabolism and Survival (Intestinal Epithelial cells): The Intestinal Epithelial cells, representing the tumor origin, exhibit strong enrichment in pathways supporting altered metabolism (e.g., amino sugar, glycosaminoglycan, sphingolipid metabolism), a hallmark of cancer. Upregulation of "HIF-1 signaling" points to adaptation to the hypoxic TME, while "Ras signaling" and "p53 signaling" indicate dysregulated proliferation and stress responses. The concurrent enrichment of "PD-L1 expression and PD-1 checkpoint pathway" suggests an intrinsic immune evasion mechanism by tumor cells. The upregulation of "Apoptosis" in tumor IECs might reflect cellular stress or an attempt at programmed cell death that is ultimately circumvented by anti-apoptotic mechanisms, or a complex regulation of the pathway.
  2. Remodeling of the Tumor Microenvironment by Stromal Cells:
  1. Immune Cell Dysfunction and Pro-tumorigenic Polarization:
  1. Shared Inflammatory and Immune Checkpoint Axis: The consistent upregulation of "PD-L1 expression and PD-1 checkpoint pathway" across nearly all major cell types (Intestinal Epithelial, Fibroblast, Macrophage, T cell, B cell, Endothelial, Plasma, Mast) in the tumor context is a striking finding. This indicates a pervasive mechanism of immune evasion orchestrated by both malignant and non-malignant cells in the colon TME.

Clinical or Translational Implications

The comprehensive GSEA findings offer several clinical and translational insights for colon cancer:

  1. Immune Checkpoint Blockade: The widespread upregulation of the PD-L1/PD-1 checkpoint pathway across diverse cell types in the colon tumor microenvironment strongly supports the rationale for immune checkpoint inhibitor therapies (e.g., anti-PD-1/PD-L1 antibodies). This is particularly relevant for colorectal cancers with high microsatellite instability (MSI-H) or deficient mismatch repair (dMMR), where such therapies have shown efficacy. The broad expression of PD-L1 suggests multiple cellular sources contribute to immune suppression. https://www.cancer.gov/about-cancer/treatment/types/immunotherapy/checkpoint-inhibitors
  2. Targeting the Tumor Microenvironment: The significant alterations in fibroblasts, macrophages, and endothelial cells highlight their critical roles in shaping the TME.
  1. Metabolic Reprogramming as a Therapeutic Vulnerability: The broad enrichment of metabolic pathways (e.g., amino sugar, glycosaminoglycan, sphingolipid metabolism) across multiple tumor-associated cell types suggests that metabolic inhibitors could represent a novel class of drugs for colon cancer, potentially targeting multiple cell populations simultaneously.
  2. Inflammation as a Driver: The pervasive activation of inflammatory signaling (AGE-RAGE, JAK-STAT, TNF) underscores the role of chronic inflammation in colon cancer progression. Modulating these inflammatory pathways, perhaps through anti-inflammatory agents or specific signaling inhibitors, could have therapeutic benefits.
  3. Differentiating Tumor Cell States: The distinct pathway enrichments observed in Intestinal Epithelial cell: Diploid_vs_others highlight potential differences in biological behavior between diploid and aneuploid/transformed cells, which might be leveraged for earlier detection or targeted intervention strategies based on ploidy status.

23. Discussion

The single-cell analysis of colorectal tissue provides a high-resolution view of the intricate changes occurring during tumorigenesis. A central finding is the pervasive genomic instability within the malignant Intestinal Epithelial cells, characterized by extensive aneuploidy and recurrent copy number variations (CNVs), including amplifications of oncogenes such as EGFR and ERBB2, and deletions of tumor suppressors like CDKN2A. The detection of aneuploid cells even in some histologically 'adjacent normal' samples suggests the presence of field cancerization or early, pre-malignant changes, underscoring the dynamic nature of tumor development.

The tumor microenvironment (TME) undergoes significant remodeling. Cellular composition shifts indicate an expansion of cancer-associated fibroblasts (CAFs) and plasma cells, along with altered macrophage subsets, notably an increase in M1-like macrophages and a decrease in M2B-like macrophages. CAFs are highly activated, expressing canonical markers such as FAP and PDGFRB, alongside immune-modulatory molecules like CD276 (B7-H3). Tumor-associated macrophages (TAMs) also exhibit a distinct surfaceome, including GPMNB, HAVCR2 (TIM-3), SIRPA, MMP14, NRP1, reflecting a phenotype that contributes to tumor growth, angiogenesis, and immune suppression.

Malignant Intestinal Epithelial cells themselves are transcriptionally reprogrammed. They display widespread upregulation of cell cycle-related genes, indicative of uncontrolled proliferation. Gene Ontology analysis reveals a shift from immune-active processes in diploid epithelial cells to hyper-proliferative, metabolically reprogrammed pathways (e.g., mTOR signaling, oxidative phosphorylation) and hypoxic adaptation (HIF-1 signaling) in tumor epithelial cells.

Cell-cell interaction analysis highlights a significantly intensified and aberrant communication network within the TME. Key pro-tumorigenic interactions include SPP1-integrin (involving macrophages and tumor epithelial cells), CXCL12-CXCR4 (fibroblast-tumor epithelial cell crosstalk, mediating migration and metastasis), HBEGF-EGFR, and VEGFA-NRP1 (tumor epithelial cell-endothelial cell interactions, promoting angiogenesis). Immune evasion is significantly driven by CD47-CD47R (SIRPA) signaling between macrophages and tumor cells, pervasive TGFB1-TGFB_receptor interactions across multiple cell types (macrophages, T cells, tumor epithelial cells), and a dampening of anti-tumor IFN-gamma signaling in TAMs. Furthermore, the GSEA results reveal a striking and widespread upregulation of the PD-L1/PD-1 checkpoint pathway across virtually all major cell types in the tumor, including malignant epithelial cells, stromal cells, and various immune cells, signifying a dominant mechanism of immune suppression.

T cells, while activated in the TME, appear functionally constrained. The increased proportion of CD4+ T cells in some tumors, coupled with the upregulation of inhibitory receptors like CTLA4 on tumor-infiltrating CD4 T cells and the broad activation of the PD-L1/PD-1 pathway, suggests a shift towards T cell exhaustion or a regulatory/immunosuppressive phenotype. This complex interplay between malignant cells, the reprogrammed stroma, and dysfunctional immune cells orchestrates an environment conducive to tumor progression and resistance to immune surveillance.

Hypotheses:

  1. The extensive aneuploidy and specific recurrent CNVs (e.g., EGFR, ERBB2 amplifications and CDKN2A deletion) in Intestinal Epithelial cells act as primary drivers of uncontrolled proliferation and malignant transformation in colorectal cancer.
  2. Cancer-associated fibroblasts (CAFs) and tumor-associated macrophages (TAMs) in the colorectal tumor microenvironment actively promote tumor growth, angiogenesis, and immune evasion through specific pro-tumorigenic and immunosuppressive ligand-receptor interactions (e.g., SPP1-integrin, CXCL12-CXCR4, CD47-SIRPA, TGFB1-TGFBR).
  3. The altered balance of CD4+ and CD8+ T cells, along with the widespread activation of the PD-L1/PD-1 checkpoint pathway and other inhibitory receptors (CTLA4, HAVCR2/TIM-3) on immune and tumor cells, contributes to T cell exhaustion and an overall immunosuppressive state in colorectal cancer.
  4. The observed metabolic reprogramming and hyper-proliferation in tumor intestinal epithelial cells, mediated by pathways like mTOR and HIF-1 signaling, represent critical vulnerabilities that can be therapeutically exploited to inhibit tumor growth.

Potential therapeutic targets:

  1. EGFR/ERBB2 signaling pathway: EGFR and ERBB2 amplifications are recurrent oncogenic events in tumor epithelial cells, driving proliferation and survival, and are well-established targets in cancer. Evidence: CNV analysis (Section 4) shows frequent amplifications of EGFR (74% frequency) and ERBB2 (37% frequency) in aneuploid Intestinal Epithelial cells. ERBB3 is also upregulated on tumor epithelial cells (Section 16). Validation: Test efficacy of anti-EGFR (e.g., Cetuximab, Panitumumab) or anti-HER2 (e.g., Trastuzumab) therapies in patient-derived tumor organoids or xenograft models stratified by these amplifications.
  2. CD47-SIRPA axis: CD47-SIRPA interactions constitute a 'don't eat me' signal, allowing tumor cells to evade macrophage phagocytosis. Blocking this interaction could re-enable anti-tumor immunity. Evidence: CCI analysis (Sections 12, 15) shows strong CD47-CD47R (SIRPA) interactions between Macrophages and Intestinal Epithelial cells (tumor cells) in the tumor microenvironment. SIRPA is highly expressed on TAMs (Section 17). Validation: Assess the impact of anti-CD47 or anti-SIRPA antibodies on macrophage phagocytosis of patient-derived tumor cells in vitro and on tumor growth in in vivo models.
  3. TGF-beta signaling pathway: TGF-beta is a potent immunosuppressive cytokine driving tumor growth, epithelial-mesenchymal transition (EMT), metastasis, and immune evasion within the TME. Inhibiting this pathway could reverse immunosuppression and inhibit tumor progression. Evidence: CCI analysis (Sections 12, 13, 14, 15) shows prominent and often elevated TGFB1-TGFbeta_receptor interactions involving Macrophages, T cells, and Intestinal Epithelial cells in tumor tissue. GSEA (Section 22) also shows enrichment of inflammatory pathways. Surfaceome analysis of normal fibroblasts also shows TGFBR3, suggesting a role in normal tissue that is dysregulated in tumor. Validation: Evaluate the effect of TGF-beta inhibitors (e.g., receptor kinase inhibitors or neutralizing antibodies) on T cell activity, macrophage polarization, EMT markers, and tumor progression in preclinical models, potentially in combination with other immunotherapies.
  4. FAP (Fibroblast Activation Protein alpha) on CAFs: FAP is a highly specific marker for activated cancer-associated fibroblasts (CAFs), which are critical for supporting tumor growth, ECM remodeling, and immune suppression. Targeting FAP can deplete pro-tumorigenic CAFs or deliver cytotoxic payloads specifically to the tumor stroma. Evidence: Surfaceome marker analysis (Section 18) identifies FAP as a strongly upregulated and prevalent marker on tumor-associated fibroblasts, contrasting with its minimal expression in adjacent normal fibroblasts. Validation: Test FAP-targeted therapies (e.g., FAP-targeting antibody-drug conjugates or CAR-T cells) in colorectal cancer models to assess their ability to deplete CAFs and inhibit tumor growth.
  5. PD-L1/PD-1 checkpoint pathway: Widespread upregulation of the PD-L1/PD-1 pathway across malignant and TME cell types indicates pervasive immune evasion, making it a prime target for immune checkpoint blockade to restore anti-tumor immunity. Evidence: GSEA (Section 22) shows strong enrichment of 'PD-L1 expression and PD-1 checkpoint pathway in cancer' across Intestinal Epithelial cells, Fibroblasts, Macrophages, T cells, B cells, Endothelial cells, Plasma cells, and Mast cells in the tumor condition. CTLA4 is also upregulated on tumor-infiltrating CD4 T cells (Section 19). Validation: Evaluate the efficacy of anti-PD-1 or anti-PD-L1 antibodies, alone or in combination with other immunotherapies or targeted agents, in preclinical colorectal cancer models, especially those demonstrating high PD-L1 expression.

Follow-up validation ideas:

  1. Spatial Transcriptomics/Proteomics: Use spatial multi-omics technologies (e.g., Visium, GeoMx, MIBI) to directly validate the spatial proximity and co-expression of key ligand-receptor pairs (e.g., SPP1-integrin, CD47-SIRPA, CXCL12-CXCR4, TGFB1-TGFBR) within the colorectal tumor microenvironment on patient tissue sections.
  2. Immunohistochemistry/Flow Cytometry: Validate the abundance and localization of specific surfaceome markers (e.g., CEACAM6, ERBB3, SDC1 on tumor epithelial cells; FAP, PDGFRB, CD276 on CAFs; GPMNB, HAVCR2, SIRPA on TAMs; CTLA4, GITR on CD4 T cells) in larger cohorts of colorectal cancer patients to assess their diagnostic, prognostic, or predictive value.
  3. In vitro/Ex vivo Perturbation Assays: Conduct co-culture experiments using patient-derived tumor organoids, CAFs, and TAMs to functionally test the impact of blocking identified therapeutic targets (e.g., anti-SPP1, anti-CD47, anti-TGFB, anti-FAP antibodies) on tumor cell proliferation, invasion, angiogenesis, and immune cell function.
  4. CRISPR/siRNA knockdown in cell lines/organoids: Investigate the functional role of highly upregulated genes (e.g., MCMs, CDKs, WEE1 in epithelial cells; MMP14, NRP1 in TAMs) in proliferation, survival, and migratory capacity of colorectal cancer cell lines or patient-derived organoids.
  5. FISH/Targeted qPCR for CNVs: Validate the presence and extent of recurrent CNVs (e.g., EGFR/ERBB2 amplification, CDKN2A deletion) in individual tumor cells and 'adjacent normal' aneuploid cells using Fluorescence In Situ Hybridization (FISH) or targeted quantitative PCR on sorted populations or tissue sections.
  6. Bulk/Single-cell Validation Cohorts: Apply the identified gene signatures and cell type proportion shifts to independent bulk or single-cell RNA-seq cohorts of colorectal cancer to confirm their robustness and association with clinical outcomes.

Limitations:

This single-cell analysis provides a comprehensive overview of cellular and molecular changes in colorectal cancer, but represents a snapshot in time and may not fully capture disease progression or dynamic responses to therapy. The spatial relationships of interacting cells, while inferred from ligand-receptor expression, require direct validation with spatial technologies. While robust, cell type assignments and ploidy inference rely on computational methods and gene expression profiles, which can have inherent limitations in resolving rare or ambiguous cell states. The specific roles of identified markers and pathways in driving tumor biology require further functional experimental investigation beyond correlative transcriptomic data. The macrophage population plot in Section 9 was uninformative and did not allow for assessment of overall macrophage abundance changes.

24. Query List

  1. Show UMAPs including condition, sample, major cell type, minor cell type, ploidy_dec, and celltype_subset in 2 columns and save it.
  2. Show major cell type scores on UMAP and save it.
  3. Show the overall Celltype_subset marker expression dot plot and save it. Keep target_cell as None and keep the other arguments at their default values.
  4. Select tumor origin cells (Intestinal Epithelial cell) and unassigned cells, group by sample, show a CNV heatmap, and include a summary of significantly amplified copy number regions. Save it.
  5. Show CNV patterns on UMAP, including major cell type, minor cell type, ploidy results, condition, and sample in 2 columns. Save it.
  6. Show a population bar plot for minor cell types and save it.
  7. Show a subset population bar plot for T cells and save it.
  8. For T cell subsets, show boxplots for those with significant differences between conditions and save it. Determine ncols appropriately based on the total number of panels.
  9. Show a subset population bar plot for Macrophages and save it.
  10. For Macrophage subsets, show boxplots for those with significant differences between conditions and save it. Determine ncols appropriately based on the total number of panels.
  11. Select tumor origin cells (Intestinal Epithelial cell) and unassigned cells, show their ploidy population as a bar plot, and save it.
  12. Show cell-cell interaction patterns by condition, including tumor origin cells (Intestinal Epithelial cell), Fibroblasts, Macrophages, and T cells. Select up to 80 interactions per condition and save it.
  13. Show cell-cell interactions by condition and save the result. Select at most 80 cell-cell interactions for each condition.
  14. Select only genes related to immune checkpoint and cell cycle pathways, and show cell-cell interactions involving these genes. Save it.
  15. For major immune cells (B cell, T cell CD4+, T cell CD8+, Macrophage, Plasma cell) and stromal cells (Fibroblast, Endothelial cell), find statistically significant differences in cell-cell interactions between conditions and show them as a dot plot. Set max_n_items_per_group = 25 and save it.
  16. Show the condition-specific markers for tumor-origin cells (Intestinal Epithelial cell) as a dot plot and save the result. Use only surfaceome markers, up to 50 markers per condition.
  17. Extract condition-specific markers for Macrophages and show them as a dot plot. Include only surfaceome markers, up to 50 per condition, and save it.
  18. Extract condition-specific markers for Fibroblasts and show them as a dot plot. Include only surfaceome markers, up to 50 per condition, and save it.
  19. Extract condition-specific markers for CD4 T cells and show them as a dot plot. Include only surfaceome markers, up to 50 per condition, and save it.
  20. For the tumor origin cell type (Intestinal Epithelial cell), select cell cycle pathway-related genes with statistically significant expression differences between conditions, and show boxplots. Set max_n_items_to_plot = 24 and determine ncols such that the panel ratio is approximately 2x3. Save it.
  21. Show Gene Ontology (GSA) analysis results for epithelial cells as bar plots and save the result.
  22. Show a dot plot of Gene Set Enrichment Analysis results for major cell types (Intestinal Epithelial cell, Fibroblast, Macrophage, T cell CD4+, T cell CD8+, B cell, Endothelial cell, Plasma cell, Mast cell). Use color map RdBu_r, set n_pws_to_show = 80, and save it.
Top