SCODiA Report by MLBI Lab

Single-Cell Atlas of Human Colon Cancer Reveals Aneuploidy-Driven Tumorigenesis and Reprogrammed Microenvironment

This report leverages single-cell RNA sequencing to map the cellular and molecular landscape of human colon cancer, comparing tumor tissue with adjacent normal tissue. We identify malignant Intestinal Epithelial cells through widespread aneuploidy and characterize their extensive cell cycle dysregulation and metabolic reprogramming. Furthermore, the tumor microenvironment exhibits significant remodeling, marked by distinct shifts in immune cell populations, particularly the expansion of immunosuppressive T regulatory cells and pro-tumorigenic M2B macrophages, alongside activated cancer-associated fibroblasts that remodel the extracellular matrix. These findings highlight key cellular and pathway alterations that drive colon cancer progression and immune evasion.

Contents

  1. Dataset overview
  2. UMAP Visualization of Single-Cell RNA-seq Data by Key Metadata
  3. Major Cell Type Score and Ploidy Distribution on UMAP
  4. Celltype Subset Marker Expression Validation
  5. Intestinal Epithelial Cell Copy Number Variation Analysis in Colon Tissue
  6. CNV-Based UMAPs for Cell Type, Ploidy, Condition, and Sample Characterization
  7. Minor Cell Type Population Analysis in Colon Tumor vs. Adjacent Normal Tissue
  8. T Cell and ILC Subset Population Analysis in Colon Tumor vs. Adjacent Normal Tissue
  9. T Cell Subset Population Dynamics in Colon Tumor vs. Adjacent Normal Tissue
  10. Macrophage Subset Reprogramming in Colon Cancer
  11. Macrophage Subset Population Shifts in Colon Tumor Microenvironment
  12. Intestinal Epithelial Cell Ploidy in Colon Cancer
  13. Colon Cancer Microenvironment Cell-Cell Interaction Analysis: Tumor vs. Adjacent Normal Tissue
  14. Immune Checkpoint and Cell Cycle Gene-Focused Cell-Cell Interaction Analysis in Colon Tissue
  15. Condition-Specific Cell-Cell Interaction Patterns in Colon Tissue
  16. Intestinal Epithelial Cell Condition-Specific Surfaceome Markers
  17. Macrophage Condition-Specific Surfaceome Markers in Colon Tissue
  18. Fibroblast Condition-Specific Surfaceome Markers in Colon Tissue
  19. T cell CD4+ Condition-Specific Surfaceome Markers in Colon Tissue
  20. Intestinal Epithelial Cell Cycle Genes Are Significantly Upregulated in Colon Tumor Microenvironment
  21. Gene Ontology (GSA) Analysis for Intestinal Epithelial Cells
  22. Gene Set Enrichment Analysis Reveals Condition- and Cell Type-Specific Pathway Alterations in Colon Tissue
  23. Discussion
  24. Query List

0. Dataset overview

데이터셋 요약

사전 계산된 분석 결과

1. UMAP Visualization of Single-Cell RNA-seq Data by Key Metadata

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents Uniform Manifold Approximation and Projection (UMAP) plots, visualizing the global transcriptional landscape of 63,689 single cells from human colon tissue. The cells are colored and grouped according to various metadata annotations: condition (Tumor vs. Adj_normal), sample, celltype_major, celltype_minor, ploidy_dec (Aneuploid vs. Diploid), and celltype_subset. These visualizations provide an essential overview of the dataset structure, the quality of cell type annotations, and the distribution of biological and technical factors within the embedding space.

Visual Summary

  1. Condition UMAP:
  1. Sample UMAP:
  1. Cell_type_major UMAP:
  1. Cell_type_minor UMAP:
  1. Ploidy_dec UMAP:
  1. Cell_type_subset UMAP:

Biological Interpretation

The UMAP visualizations collectively provide a powerful overview of the cellular heterogeneity in human colon tissue, particularly in the context of cancer.

Annotation Notes

The UMAP plots demonstrate excellent quality in cell type annotation at multiple levels of granularity. The distinct clustering of celltype_major, celltype_minor, and celltype_subset indicates that the cell identities are well-defined by their transcriptional profiles. Furthermore, the clear separation of aneuploid cells from diploid cells, particularly within the epithelial compartment of tumor samples, provides strong validation for the ploidy_dec annotation as a reliable indicator of malignant cells. The minimal presence of 'unassigned' or 'unclear' labels further supports the comprehensive nature of the annotations.

---

References:

  1. Aneuploidy as a Hallmark of Cancer:
  1. Chromosomal Instability and Cancer:

2. Major Cell Type Score and Ploidy Distribution on UMAP

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis visualizes the distribution of major cell type scores and inferred ploidy status across the UMAP embedding of single-cell RNA-seq data from human Colon tissue. The primary goal is to assess the spatial organization of different cell populations and to examine the ploidy status (Aneuploid vs. Diploid) in relation to these cell types, particularly the Intestinal Epithelial cells, which are identified as the tumor origin cell type.

Visual Summary

The UMAP plots effectively display the distinct clusters of major cell types based on their computed scores and ground truth annotations, along with the distribution of ploidy.

Biological Interpretation

The visualizations provide crucial insights into cell identity validation and tumor biology within the colon tissue sample:

  1. Validation of Cell Type Annotation and Scoring: The strong correspondence between the HiCAT_major_score heatmaps and the celltype_major categorical assignments confirms the accuracy and robustness of the cell type identification. Cells with high scores for a particular major cell type are indeed annotated as that cell type, indicating a reliable classification of cell populations. This is a critical step for downstream analyses, ensuring that differential gene expression or cell-cell interaction findings are attributed to correctly identified cell types.
  2. Malignant Epithelial Cell Identification via Ploidy: The striking co-localization of Aneuploid cells with the Intestinal Epithelial cell cluster is a key finding. Given that "Intestinal Epithelial cell" is specified as the "Tumor origin celltype" and aneuploidy is a well-established hallmark of cancer [1], this observation strongly suggests that the Aneuploid Intestinal Epithelial cells represent the malignant tumor cell population. Conversely, the Diploid Intestinal Epithelial cells likely represent normal epithelial cells, potentially from adjacent normal tissue or non-malignant epithelial components within the tumor.
  3. Tumor Microenvironment Composition: The presence of various immune cells (T cells, B cells, Myeloid cells, Mast cells) and stromal cells (Fibroblasts, Endothelial cells) that are largely Diploid suggests these populations constitute the tumor microenvironment (TME) or elements of the adjacent normal tissue. These non-malignant cells, while generally diploid, can play significant roles in tumor progression, immune evasion, and therapeutic response [2].

Annotation Notes

References

  1. Aneuploidy as a hallmark of cancer:
  1. Tumor microenvironment composition and function:

3. Celltype Subset Marker Expression Validation

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents a dot plot visualizing the expression of marker genes across various celltype_subset populations from single-cell RNA-seq data of human colon tissue. The primary goal is to validate the assigned cell identities by examining whether the identified marker genes exhibit specific and high expression patterns within their corresponding cell type subsets. The plot displays the mean expression level (color intensity) and the fraction of cells expressing the gene (dot size) for each gene-cell type pair. Markers were selected to prioritize surfaceome genes with high specificity, as indicated by the surfaceome_only: True and rem_mkrs_common_in_N_groups_or_more: 3 parameters.

Visual Summary

The dot plot demonstrates a clear and well-defined pattern of marker gene expression, with distinct clusters of genes highly specific to particular celltype_subset populations.

Biological Interpretation

The marker expression patterns strongly support the assigned celltype_subset annotations, indicating robust cell type identification within the colon tissue data.

Intestinal Epithelial Cells (IECs)

Immune Cells

Stromal and Endothelial Cells

Annotation Notes

The comprehensive and specific marker gene expression patterns observed across the celltype_subset populations provide strong evidence for the high quality and accuracy of the single-cell annotations. The selection of surfaceome markers further enhances the practical utility of these findings for future experimental validation, such as flow cytometry or immunohistochemistry. The distinct expression profiles confirm that each celltype_subset represents a biologically coherent population within the human colon.

---

References:

  1. SOX9: Furuyama K, Kawaguchi Y, Akiyama H, et al. Continuous cell supply from a Sox9-expressing progenitor zone in adult liver, exocrine pancreas and stomach. Nat Genet. 2011;43(11):1091-1100. PubMed search link for SOX9 and intestinal stem cell
  2. ASCL2: Sangiorgi E, Capecchi MR. Bmi1 is indispensable for intestinal stem cell self-renewal. Nat Genet. 2008;40(7):915-920. PubMed search link for ASCL2 intestinal crypt
  3. FABP1: GeneCards for FABP1. GeneCards link for FABP1
  4. CDX2, KRT20, VIL1: GeneCards for CDX2, KRT20, VIL1. GeneCards link for CDX2, GeneCards link for KRT20, GeneCards link for VIL1
  5. MUC2, TFF3: GeneCards for MUC2, TFF3. GeneCards link for MUC2, GeneCards link for TFF3
  6. LYZ: GeneCards for LYZ. GeneCards link for LYZ
  7. GP2: GeneCards for GP2. GeneCards link for GP2
  8. DCLK1: GeneCards for DCLK1. GeneCards link for DCLK1
  9. CHGA: GeneCards for CHGA. GeneCards link for CHGA
  10. POU2F2, CD79A/B: GeneCards for POU2F2, CD79A, CD79B. GeneCards link for POU2F2, GeneCards link for CD79A, GeneCards link for CD79B
  11. MZB1, XBP1, PRDM1: GeneCards for MZB1, XBP1, PRDM1. GeneCards link for MZB1, GeneCards link for XBP1, GeneCards link for PRDM1
  12. CD83, CD86, IRF7, IRF8: GeneCards for CD83, CD86, IRF7, IRF8. GeneCards link for CD83, GeneCards link for CD86, GeneCards link for IRF7, GeneCards link for IRF8
  13. MSR1: GeneCards for MSR1. GeneCards link for MSR1
  14. TPSAB1: GeneCards for TPSAB1. GeneCards link for TPSAB1
  15. GZMB: GeneCards for GZMB. GeneCards link for GZMB
  16. TBX21, STAT4: GeneCards for TBX21, STAT4. GeneCards link for TBX21, GeneCards link for STAT4
  17. RORC: GeneCards for RORC. GeneCards link for RORC
  18. GATA3: GeneCards for GATA3. GeneCards link for GATA3
  19. FOXP3: GeneCards for FOXP3. GeneCards link for FOXP3
  20. ILCs: Spits H, Cupedo T. Innate lymphoid cells: innate protectors and regulators of immunity. Cell. 2012;146(5):750-761. PubMed link for ILCs overview
  21. COL1A1, COL1A2, DCN, LUM: GeneCards for COL1A1, COL1A2, DCN, LUM. GeneCards link for COL1A1, GeneCards link for COL1A2, GeneCards link for DCN, GeneCards link for LUM
  22. ANGPT2, DLL4: GeneCards for ANGPT2, DLL4. GeneCards link for ANGPT2, GeneCards link for DLL4
  23. PROX1, PDPN: GeneCards for PROX1, PDPN. GeneCards link for PROX1, GeneCards link for PDPN
  24. S100B: GeneCards for S100B. GeneCards link for S100B
  25. ACTA2, TAGLN: GeneCards for ACTA2, TAGLN. GeneCards link for ACTA2, GeneCards link for TAGLN

4. Intestinal Epithelial Cell Copy Number Variation Analysis in Colon Tissue

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates copy number variations (CNVs) specifically within Intestinal Epithelial cells, which are identified as the tumor origin cell type, across various patient samples from colon tissue. The results include a heatmap displaying log2(CNR) (Copy Number Ratio) values across the genome for individual cells grouped by sample, alongside a summary heatmap detailing regions with significantly amplified copy numbers and associated genes. This helps to identify recurrent genomic alterations characteristic of the tumor state and understand cellular heterogeneity in terms of ploidy and CNV burden.

Visual Summary

CNV Heatmap (log2(CNR) by Genomic Spot and Cell Group)

The first heatmap visualizes the log2(CNR) values across all chromosomes for Intestinal Epithelial cells.

Summary of Significantly Amplified Copy Numbers

The second visualization provides a concise summary of the most significantly amplified cytogenetic bands across samples, with associated frequencies.

Key Amplified Regions and Genes:

Biological Interpretation

The observed CNVs in Intestinal Epithelial cells from colon tumors provide crucial insights into the genetic basis of colon cancer development and progression.

Clinical or Translational Implications

5. CNV-Based UMAPs for Cell Type, Ploidy, Condition, and Sample Characterization

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents UMAP (Uniform Manifold Approximation and Projection) visualizations based on estimated Copy Number Variation (CNV) patterns (obsm['X_cnv']) across 63,689 single cells from Colon tissue. The UMAPs are colored by various metadata annotations: major cell type (celltype_major), minor cell type (celltype_minor), ploidy status (ploidy_dec), tissue condition (condition), and individual patient sample (sample). The primary goal is to visualize how distinct cell populations, particularly malignant cells, cluster based on their genomic alteration profiles (CNVs and ploidy).

Visual Summary

  1. celltype_major and celltype_minor UMAPs:
  1. ploidy_dec UMAP:
  1. condition UMAP:
  1. sample UMAP:

Biological Interpretation

The CNV-based UMAPs provide robust evidence for distinguishing malignant cells from non-malignant cells in the Colon tissue samples.

Annotation Notes

The UMAPs generated using CNV data are highly effective for:

6. Minor Cell Type Population Analysis in Colon Tumor vs. Adjacent Normal Tissue

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents a population bar plot illustrating the proportional distribution of minor cell types across individual samples from both "Adj_normal" (adjacent normal) and "Tumor" conditions within the colon tissue. Each bar represents a single sample, with stacked segments indicating the relative abundance of different cell types, providing insight into the cellular heterogeneity and shifts associated with the tumor microenvironment.

Visual Summary

The visualization clearly displays the cellular composition for 9 adjacent normal samples and 20 tumor samples.

Biological Interpretation

The observed cellular landscape provides crucial insights into the profound biological changes occurring in the colon tumor microenvironment (TME) compared to adjacent normal tissue.

  1. Neoplastic Expansion: The most striking finding is the overwhelming dominance of Intestinal Epithelial cells in tumor samples. Given that Intestinal Epithelial cells are identified as the "Tumor origin celltype" in the data context, this expansion reflects the uncontrolled proliferation characteristic of colorectal carcinogenesis. This massive expansion of tumor cells dramatically alters the overall cellular proportions, effectively diluting the relative contribution of other stromal and immune components.
  2. Altered Tumor Microenvironment: The relative decrease in the proportional abundance of non-epithelial cells, such as Fibroblasts, Endothelial cells, and various immune cells, indicates a restructuring of the TME.
  1. Colon Tissue Specificity: The presence of cell types like Smooth muscle cells, Enteric glial cells, and Enteric neurons in the adjacent normal tissue is consistent with the complex multi-tissue structure of the colon. Their altered proportions in tumors reflect the invasive and disruptive nature of cancer on the surrounding tissue architecture.

Clinical or Translational Implications

The findings from this cell type population analysis have several important clinical and translational implications for colon cancer:

  1. Biomarker for Disease State: The significant shift in cellular composition, particularly the increased dominance of Intestinal Epithelial cells, serves as a clear biomarker distinguishing tumor tissue from normal tissue. This confirms the quality of cell type annotation and the ability to differentiate pathological states at a cellular level.
  2. Heterogeneity and Patient Stratification: The observed sample-to-sample heterogeneity in immune and stromal cell infiltration within the tumor condition suggests that colon cancers from different patients can have distinct tumor microenvironments. This heterogeneity is critical for patient stratification. For example, tumors with high immune infiltration might be more responsive to immunotherapies, while those with a "cold" phenotype might require different therapeutic strategies or pre-treatment to enhance immune cell recruitment [2].
  3. Therapeutic Targeting of the TME: Understanding the specific cellular composition of the TME is crucial for developing and applying targeted therapies. For instance, if a specific patient's tumor exhibits high CAF content, therapies targeting CAF-mediated desmoplasia or immunosuppression might be beneficial. Similarly, profiling immune cell populations can guide the selection of immunotherapeutic agents.
  4. Foundation for Deeper Analysis: This population analysis serves as a foundational step. Further analyses, such as differential gene expression (DEG), gene set enrichment analysis (GSEA), and cell-cell interaction (CCI) studies, within specific cell types and conditions will be essential to uncover the functional consequences of these population shifts and identify potential therapeutic targets or mechanisms of resistance.

---

References:

  1. Macrophages in TME: Orecchioni, S., et al. (2019). Macrophage Polarization: Different Molecular Mechanisms, Different Biological Functions. *Trends Immunol*. PubMed search: Macrophage polarization tumor microenvironment
  2. Tumor Immunophenotypes: Galon, J., et al. (2020). The Immunoscore (for colorectal cancer): From research to clinical practice. *J Immunother Cancer*. PubMed search: Immunoscore tumor stratification

7. T Cell and ILC Subset Population Analysis in Colon Tumor vs. Adjacent Normal Tissue

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents a population bar plot illustrating the relative proportions of various T cell and Innate Lymphoid Cell (ILC) subsets across individual samples from both "Adj_normal" (adjacent normal colon tissue) and "Tumor" conditions. The visualization provides insights into the immune cell composition changes occurring in the colorectal tumor microenvironment.

Visual Summary

The stacked bar plots visually represent the cellular composition of T cells and ILCs for each sample, grouped by condition (Adj_normal vs. Tumor).

Biological Interpretation

The observed shifts in T cell subset populations provide critical biological insights into the immune landscape of colorectal cancer.

Clinical or Translational Implications

These findings have several important clinical and translational implications for colorectal cancer:

8. T Cell Subset Population Dynamics in Colon Tumor vs. Adjacent Normal Tissue

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates the proportional changes of various T cell and innate lymphoid cell (ILC) subsets within colon tumor tissue compared to adjacent normal tissue. The box plots visually represent the celltype proportion for each subset across the two conditions, with statistical significance (p-values) indicating differences. This provides insights into the immune landscape shifts occurring in the tumor microenvironment.

Visual Summary

The box plots reveal statistically significant shifts in the proportions of several T cell and ILC subsets when comparing colon tumor tissue to adjacent normal tissue.

Biological Interpretation

The observed shifts in T cell and ILC populations are consistent with known immune evasion mechanisms and inflammatory responses within the tumor microenvironment of colorectal cancer.

Decreased LTI, ILC1, and ILC2 cells in Tumors:

Overall, the data paints a picture of immune dysregulation within the colon tumor microenvironment, characterized by an enrichment of immunosuppressive (Treg) and pro-inflammatory (Th17) T cell subsets, alongside a depletion of certain innate anti-tumor (ILC1) and regulatory (LTI, ILC2) populations.

Clinical or Translational Implications

These findings have significant clinical implications for understanding colon cancer progression and designing immunotherapeutic strategies:

9. Macrophage Subset Reprogramming in Colon Cancer

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents a population bar plot illustrating the relative proportions of different macrophage subsets (M1, M2A, M2B, M2C, M2D) within the overall macrophage population across individual samples from both "Adj_normal" (adjacent normal colon tissue) and "Tumor" (colon tumor tissue) conditions. The goal is to identify shifts in macrophage polarization states associated with the tumor microenvironment.

Visual Summary

The visualization displays stacked bar plots for macrophage subsets, stratified by individual samples and conditions (Adj_normal vs. Tumor).

Biological Interpretation

Macrophages are critical components of the tumor microenvironment (TME) and play diverse roles in cancer progression, largely dictated by their polarization state.

  1. Shift Towards M1 and M2B Dominance in Tumor: The most striking observation is the clear shift in macrophage composition in tumor samples.
  1. Reduced M2D in Tumor: While M2D macrophages are often associated with pro-tumorigenic functions, angiogenesis, and immune suppression in various cancers, this analysis shows that M2D is generally low in colon tumor samples, especially compared to the high proportions observed in certain adjacent normal samples (e.g., SMC03-N). This suggests that in this specific colon cancer cohort, M2D might not be the primary pro-tumorigenic macrophage subset or that its role is diminished, with other M2 types (like M2B) taking precedence in shaping the TME.
  2. Heterogeneity in Adjacent Normal Tissue: The varied macrophage profiles in adjacent normal samples highlight the baseline heterogeneity of resident macrophages in the healthy colon, which can be influenced by local immunological states, commensal microbiota, or individual differences. The distinct shift in the tumor context indicates significant reprogramming due to cancer-specific stimuli.

Overall, the macrophage landscape in colon cancer appears to be characterized by a significant presence of M1 macrophages, coupled with a prominent M2B population. This suggests a TME that is both pro-inflammatory (M1) and immunomodulatory/pro-tumorigenic (M2B), representing a complex interplay of immune responses.

Clinical or Translational Implications

  1. Biomarker Potential: The distinct shifts in macrophage subsets (especially the relative increase in M1 and M2B, and reduction in M2D compared to some normal tissues) could serve as potential diagnostic or prognostic biomarkers for colon cancer. The specific macrophage polarization signature might correlate with disease stage, aggressiveness, or patient outcomes.
  2. Immunotherapeutic Targets: Understanding the dominant macrophage subsets in the colon TME opens avenues for targeted immunotherapies.

The observed macrophage dynamics underscore the importance of precision immunomodulation in colon cancer, moving beyond a simplistic M1/M2 dichotomy to address the nuanced roles of individual M2 subtypes.

10. Macrophage Subset Population Shifts in Colon Tumor Microenvironment

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates the proportional changes of specific macrophage subsets (M1, M2A, M2B) in the colon tissue, comparing tumor samples against adjacent normal tissues. The goal is to identify macrophage populations that exhibit statistically significant differences in their representation within the tumor microenvironment, providing insights into their potential roles in tumor progression or suppression. The analysis used a p-value cutoff of 0.1 for statistical significance, and the proportions are calculated relative to the total number of cells in the respective samples.

Visual Summary

The box plots display the celltype proportion for three macrophage subsets: Mac (M2A), Mac (M2B), and Mac (M1), stratified by 'Adj_normal' and 'Tumor' conditions.

Biological Interpretation

Macrophages are highly plastic immune cells that polarize into distinct functional states, broadly categorized as M1 (classically activated) and M2 (alternatively activated) phenotypes, each with diverse roles in inflammation, tissue repair, and immune regulation. M2 macrophages are further sub-categorized (M2A, M2B, M2C, M2D) based on their activation pathways and effector functions. The observed shifts in macrophage subsets within the colon tumor microenvironment are biologically significant:

Collectively, these findings suggest a dominant shift towards M2B-like macrophage polarization within the colon tumor microenvironment, favoring an immunosuppressive and pro-tumorigenic milieu.

Clinical or Translational Implications

The distinct shifts in macrophage subsets, particularly the significant increase of M2B macrophages and decrease of M2A macrophages in colon tumors, have important clinical implications:

11. Intestinal Epithelial Cell Ploidy in Colon Cancer

Report figure

[Analysis Visualization Results]...

Analysis Overview

이 분석은 단일 세포 RNA 시퀀싱(scRNA-seq) 데이터를 사용하여 대장 조직 내 장 상피 세포(Intestinal Epithelial cells)의 ploidy 상태(배수성, 즉 염색체 수의 비정상)를 평가한 것입니다. Adj_normal(인접 정상 조직) 및 Tumor(종양 조직) 샘플에서 각 세포의 ploidy 상태(Aneuploid, Diploid, Unclear) 분포를 막대 그래프로 시각화하여, 종양 발생과 관련된 장 상피 세포의 유전체 변화를 탐색합니다. Intestinal Epithelial cell은 데이터 컨텍스트에 따라 종양의 기원 세포(Tumor origin celltype)로 지정되어 있습니다.

Visual Summary

제공된 막대 그래프는 Adj_normal 및 Tumor 조건에서 각 샘플별 장 상피 세포의 ploidy 분포를 보여줍니다.

Biological Interpretation

이 분석 결과는 대장암에서 장 상피 세포의 유전체 불안정성(genomic instability)과 관련된 중요한 생물학적 통찰을 제공합니다.

  1. 암의 특징으로서의 Aneuploidy: 정상 장 상피 세포가 거의 전적으로 Diploid인 반면, 많은 종양 조직의 장 상피 세포에서 Aneuploidy가 지배적으로 나타나는 것은 암 발생 및 진행의 핵심적인 특징인 염색체 이수성(aneuploidy)을 명확히 보여줍니다. Aneuploidy는 염색체 수의 비정상적인 변화로, 세포 증식 촉진, 세포사멸 회피, 전이 능력 강화 등 암세포 특성을 유발하는 데 기여합니다 PubMed search: aneuploidy cancer hallmark.
  2. 종양 내 이질성: 종양 샘플 내에서 Aneuploid 세포의 비율이 다양하다는 것은 종양 미세 환경의 복잡성 또는 종양 자체의 이질성을 반영할 수 있습니다.

Clinical or Translational Implications

이 분석은 대장암의 진단, 분류 및 치료 전략 수립에 다음과 같은 잠재적인 임상적 또는 번역적 함의를 가집니다.

  1. 종양 세포 식별 및 순도 평가: Aneuploidy는 종양 기원 세포인 장 상피 세포를 악성 세포로 식별하는 강력한 마커로 활용될 수 있습니다. 단일 세포 데이터셋에서 종양 세포를 정확히 구분하고, 샘플 내 종양 세포의 순도를 평가하는 데 기여합니다.
  2. 질병 진행 및 예후 마커: Aneuploidy의 정도는 암의 공격성, 진행 단계, 그리고 환자의 예후와 연관될 수 있습니다 GeneCards: Aneuploidy. 높은 Aneuploidy는 더 공격적인 종양을 나타낼 수 있으며, 이는 특정 치료 전략 선택에 영향을 미칠 수 있습니다.
  3. 치료 반응 예측: Aneuploidy는 일부 항암 치료에 대한 반응성을 예측하는 인자로 연구되기도 합니다. 예를 들어, 염색체 불안정성을 표적으로 하는 약물의 개발과 관련하여 Aneuploid 세포의 존재는 해당 치료의 잠재적 대상 환자군을 식별하는 데 도움이 될 수 있습니다.
  4. 저종양성(low tumor cellularity) 샘플 해석: 주로 Diploid 세포를 보이는 종양 샘플의 존재는 종양 세포 함량이 낮은 생검 샘플의 해석에 주의를 기울여야 함을 시사합니다. 이러한 경우, ploidy 상태 외에 다른 유전자 발현 패턴이나 유전적 변이 정보를 통합하여 종양 여부와 특성을 판단하는 것이 중요합니다.

12. Colon Cancer Microenvironment Cell-Cell Interaction Analysis: Tumor vs. Adjacent Normal Tissue

Report figure

[Analysis Visualization Results]...

Analysis Overview

이 분석은 인접 정상 조직(Adj_normal)과 종양(Tumor) 조건에서 Intestinal Epithelial cell, Fibroblast, Macrophage, T cell CD4+, T cell CD8+ 간의 세포-세포 상호작용(Cell-Cell Interaction, CCI) 패턴을 비교합니다. CellPhoneDB를 사용하여 리간드-수용체 상호작용을 예측하고, 그 결과를 닷 플롯으로 시각화하여 상호작용의 유의성(p-value)과 평균 발현 수준을 나타냅니다. 각 조건에서 최대 80개의 가장 유의미한 상호작용 쌍을 보여주며, 이를 통해 종양 미세환경에서 나타나는 세포 간 통신 변화를 이해하고자 합니다.

Visual Summary

두 개의 닷 플롯은 각각 'Adj_normal'과 'Tumor' 조건에서의 세포-세포 상호작용을 보여줍니다. X축은 리간드-수용체 쌍을, Y축은 상호작용하는 세포 쌍을 나타냅니다. 닷의 크기는 상호작용의 유의성(-log10(p))을 나타내며(큰 닷일수록 유의), 색상은 리간드와 수용체의 평균 발현 수준(-log2(m))을 나타냅니다(밝은 노란색에 가까울수록 발현 수준이 높음).

Adj_normal 조건 (상단 플롯):

Tumor 조건 (하단 플롯):

Biological Interpretation

이 분석 결과는 결장암 발병 및 진행 과정에서 세포-세포 상호작용 네트워크에 상당한 재편이 일어남을 명확히 보여줍니다.

  1. 미세환경의 재구성 (Remodeling of the TME):
  1. 염증 및 면역 회피 환경 조성:
  1. 혈관신생 및 세포 이동 촉진:
  1. 세포외 기질 상호작용의 변화:

Clinical or Translational Implications

이러한 세포-세포 상호작용 패턴의 변화는 결장암 치료를 위한 중요한 표적을 제시합니다.

이러한 결과는 특정 리간드-수용체 쌍에 대한 추가적인 실험적 검증(예: 특정 상호작용을 차단하는 항체 또는 소분자 억제제 사용)과 생체 내 모델에서의 효능 평가를 통해 새로운 결장암 치료제 개발로 이어질 수 있습니다.

13. Immune Checkpoint and Cell Cycle Gene-Focused Cell-Cell Interaction Analysis in Colon Tissue

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates cell-cell interactions (CCI) in human colon tissue, comparing "Adj_normal" (adjacent normal) and "Tumor" conditions. The focus is specifically on a curated list of genes associated with immune checkpoint and cell cycle pathways, identified from single-cell RNA sequencing data. CellPhoneDB was used to infer ligand-receptor interactions between different cell types (major and minor cell types, with ploidy information for epithelial cells). The results highlight significant interactions by p-value (dot size) and interaction strength (log2(mean expression) by dot color), offering insights into how these critical pathways mediate communication within the tissue microenvironment under normal and diseased states.

Visual Summary

CCI for Adj_normal

The "Adj_normal" dot plot displays a diverse array of cell-cell interactions, particularly involving T CD8+ cells, Fibroblasts, Endothelial cells, and Diploid Intestinal Epithelial cells. Key observations include:

CCI for Tumor

The "Tumor" dot plot reveals a distinct pattern of interactions, with a strong emphasis on Macrophages, T cells (CD8+ and CD4+), and Aneuploid Intestinal Epithelial cells (representing tumor cells).

Biological Interpretation

The differential cell-cell interactions observed between adjacent normal tissue and tumor tissue provide critical insights into the altered microenvironment in colon cancer, particularly concerning immune regulation and cellular growth pathways.

  1. EGFR Pathway Remodeling: In "Adj_normal" tissue, EGFR signaling (AREG-EGFR, HBEGF-EGFR) primarily supports homeostasis and repair, mediating communication between Diploid Intestinal Epithelial cells, Fibroblasts, and Endothelial cells. This is essential for tissue integrity. In "Tumor" tissue, however, EREG-EGFR interactions become prominent, notably between Aneuploid Intestinal Epithelial cells (tumor cells) and immune cells (T CD8+, Macrophages). Epiregulin (EREG) is known to be overexpressed in various cancers and can promote tumor cell proliferation and survival, and also influence the immune landscape by impacting immune cell function [1]. The shift from AREG/HBEGF to EREG as a dominant EGFR ligand in the tumor context suggests a reprogramming of growth factor signaling to support oncogenic processes.
  2. Immune Cell Activation and Interaction: The "Tumor" microenvironment exhibits strong CD86-CD28 interactions between T CD8+|Macrophage and Macrophage|T CD4+. CD86 (on antigen-presenting cells like macrophages) binding to CD28 (on T cells) is a critical co-stimulatory signal required for full T cell activation [2]. While this can indicate active anti-tumor immune responses, persistent co-stimulation in the tumor microenvironment can also lead to T cell exhaustion, especially in the context of other inhibitory signals not highlighted in this specific gene panel. The sustained IFNG-Type II IFN receptor interactions in both conditions, especially among T CD8+ cells, points to the presence of IFN-gamma, a key cytokine in anti-tumor immunity.
  3. Immune-Stromal Interactions in Homeostasis vs. Disease: In "Adj_normal," interactions involving Fibroblasts and T cells through TGFB1-TGFB_receptor suggest a role for TGF-beta in immune regulation and maintaining tissue architecture. TGF-beta is a potent immunosuppressive cytokine and a driver of fibrosis [3]. While not prominent in the top interactions shown for "Tumor" here, its absence from the top 80 pairs does not mean it's inactive, but rather that other pathways like CD86-CD28 and EREG-EGFR may be relatively more dominant or frequent in the tumor microenvironment based on the selection criteria.
  4. Ploidy-Specific Interactions: The distinction between "Diploid Intestinal Epi" in "Adj_normal" and "Aneuploid Intestinal Epi" in "Tumor" is critical. Aneuploidy is a hallmark of cancer. The interactions involving Aneuploid Intestinal Epithelial cells with immune cells (T CD8+, Macrophages) via EREG-EGFR axis directly implicate tumor cells in shaping the tumor microenvironment through growth factor signaling.
  5. Ambiguity in Gene Pair Interpretation: The presence of "LCK-CD8_receptor" and "CD93-IFNGR1" warrants cautious interpretation. LCK is an intracellular tyrosine kinase essential for T cell receptor signaling; it is not a secreted ligand or a transmembrane receptor. Its appearance as an interaction partner may reflect its critical involvement in CD8+ T cell signaling pathways rather than a direct ligand-receptor binding event, or it could be a tool artifact. Similarly, CD93 is a C-type lectin that is not known to directly bind IFNGR1 (which binds IFNG). These specific entries might represent indirect associations or require further validation beyond canonical ligand-receptor databases.

Clinical or Translational Implications

The identified cell-cell interactions, particularly those altered in the tumor microenvironment and involving immune checkpoint and cell cycle-related genes, offer several translational implications for colon cancer.

  1. EGFR Pathway as a Therapeutic Target: The prominent EREG-EGFR interactions between Aneuploid Intestinal Epithelial cells and immune cells in the tumor environment highlight the EGFR pathway as a potential therapeutic target. EGFR inhibitors are already used in metastatic colorectal cancer [4], and these findings suggest that targeting this specific ligand-receptor pair (EREG-EGFR) could be particularly relevant in preventing tumor proliferation and modulating immune evasion.
  2. Modulating Immune Co-stimulation: The strong CD86-CD28 interactions in the tumor microenvironment emphasize the activity of this T cell co-stimulatory axis. While CD28 activation is generally pro-inflammatory, the balance with inhibitory signals (like CTLA-4 or PD-1) is crucial. Understanding the precise context of these interactions (e.g., leading to activation vs. exhaustion) could inform strategies for immunotherapy, potentially combining CD28 agonists (if exhaustion is a major factor) or targeting upstream factors that regulate CD86 expression on macrophages.
  3. Context-Specific Immune Checkpoint Regulation: The persistent IFNG-Type II IFN receptor signaling suggests that IFN-gamma is present and active. While IFN-gamma can have anti-tumor effects, it also can induce PD-L1 expression on tumor cells, leading to immune evasion [5]. Further investigation into the balance of these signals, possibly through correlating with PDCD1 (PD-1) or CD274 (PD-L1) expression, would be crucial.
  4. Biomarker Discovery: The identified critical ligand-receptor pairs (e.g., EREG-EGFR, CD86-CD28) could serve as potential biomarkers for patient stratification or response prediction to targeted therapies or immunotherapies. Increased expression of EREG or specific immune cell populations highly engaging in CD86-CD28 interactions might indicate particular disease subtypes responsive to specific interventions.

---

References

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

[2] PubMed search for CD28-CD86 signaling T cell activation: https://pubmed.ncbi.nlm.nih.gov/?term=CD28+CD86+T+cell+activation

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

[4] PubMed search for EGFR inhibitors colorectal cancer: https://pubmed.ncbi.nlm.nih.gov/?term=EGFR+inhibitors+colorectal+cancer

[5] PubMed search for IFN-gamma PD-L1 upregulation: https://pubmed.ncbi.nlm.nih.gov/?term=IFN-gamma+PD-L1+upregulation

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

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis identifies cell-cell interactions (CCIs) involving major immune and stromal cells (B cell, Myeloid cell, Stromal cell, T cell, Mast cell) that differ significantly between 'Tumor' and 'Adj_normal' conditions in human colon tissue. The results are visualized as a dot plot, where dot size reflects the statistical significance (-log10(p-value)) of the interaction, and dot color indicates the scaled interaction strength. The analysis aimed to reveal how the communication landscape within the tumor microenvironment is altered compared to adjacent normal tissue.

Visual Summary

The dot plot effectively illustrates distinct patterns of cell-cell interactions between 'Adj_normal' and 'Tumor' conditions, across different samples.

Overall Condition-Specific Patterns:

Differential Interaction Landscape:

Cross-Condition Comparison for Samples:

Biological Interpretation

The observed shifts in cell-cell interactions underscore a significant remodeling of the tumor microenvironment (TME) in colon cancer, particularly involving stromal and immune cells.

Pro-tumorigenic Immunomodulation:

Clinical or Translational Implications

The distinctive and robust cell-cell interaction patterns identified in colon tumor samples hold significant clinical and translational potential.

15. Intestinal Epithelial Cell Condition-Specific Surfaceome Markers

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis identifies condition-specific surfaceome markers in Intestinal Epithelial cells, comparing Tumor and Adjacent Normal (Adj_normal) tissues from human colon single-cell RNA-seq data. The dot plot visualizes the expression levels and prevalence of the top 50 surface markers for each condition across individual samples. This approach prioritizes the discovery of cell-type-specific surface markers that could serve as potential therapeutic targets or biomarkers.

Visual Summary

The dot plot effectively displays differential gene expression patterns of surfaceome markers in Intestinal Epithelial cells across various samples categorized by their condition (Adj_normal vs. Tumor) and ploidy status (Diploid vs. Aneuploid within Tumor).

Biological Interpretation

The analysis reveals distinct molecular signatures on the surface of Intestinal Epithelial cells, differentiating healthy colon tissue from tumor tissue, and further highlighting differences related to cellular ploidy within the tumor.

Clinical or Translational Implications

The identified condition-specific surfaceome markers in Intestinal Epithelial cells offer several potential clinical and translational avenues:

Therapeutic Targets for Colorectal Cancer:

16. Macrophage Condition-Specific Surfaceome Markers in Colon Tissue

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis identifies surfaceome markers that are specifically expressed by Macrophage populations in either adjacent normal colon tissue or tumor colon tissue. By comparing gene expression patterns across different samples and conditions (Adj_normal vs. Tumor), the goal is to pinpoint surface proteins that can distinguish macrophages in these distinct microenvironments. This is crucial for understanding macrophage plasticity, their functional roles in the tumor microenvironment (TME), and for identifying potential diagnostic or therapeutic targets.

Visual Summary

The dot plot effectively visualizes the expression patterns of 30 selected surfaceome markers across various patient samples, grouped by condition (Adj_normal vs. Tumor). Each dot's size represents the fraction of cells expressing the gene within that sample group, while its color intensity reflects the mean expression level.

Markers Enriched in Adjacent Normal Macrophages:

A set of markers shows high expression (dark red, large dots) almost exclusively in "Adj_normal" samples (SMC01-N to SMC06-N). Key examples include:

Markers Enriched in Tumor Macrophages:

Conversely, a distinct set of markers is highly expressed (dark red, large dots) predominantly in "Tumor" samples (SMC10-T to SMC19-T). Notable examples include:

Biological Interpretation

The distinct sets of surfaceome markers highlight significant differences in the biological states and functions of macrophages residing in the adjacent normal colon tissue compared to those infiltrating colon tumors.

Macrophages in Adjacent Normal Tissue

The markers enriched in adjacent normal macrophages suggest a role in tissue homeostasis, immune regulation, and basal macrophage functions:

These markers collectively describe a macrophage population poised for tissue maintenance, iron metabolism, and potentially an anti-inflammatory or regulatory role characteristic of healthy tissue resident macrophages.

Macrophages in Tumor Tissue (Tumor-Associated Macrophages, TAMs)

The markers elevated in tumor macrophages indicate a distinct phenotype, often associated with pro-tumoral functions within the tumor microenvironment:

These markers collectively define a population of tumor-associated macrophages (TAMs) that are highly active in ECM remodeling, immune modulation, and likely contribute to the immunosuppressive and pro-metastatic environment characteristic of colon cancer.

Clinical or Translational Implications

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

17. Fibroblast Condition-Specific Surfaceome Markers in Colon Tissue

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis identifies and visualizes condition-specific surfaceome markers for Fibroblasts in Colon tissue, comparing "Tumor" and "Adj_normal" conditions. The dot plot displays the mean expression level (color intensity) and the fraction of cells expressing the marker (dot size) for each gene across individual samples, grouped by condition. This approach helps pinpoint cell-surface proteins that are differentially expressed, offering insights into condition-specific fibroblast phenotypes and potential therapeutic targets.

Visual Summary

The dot plot clearly differentiates two major groups of surfaceome markers corresponding to the "Adj_normal" and "Tumor" conditions in Fibroblast cells.

Biological Interpretation

The observed condition-specific surfaceome markers reflect distinct functional states of fibroblasts in normal colon tissue versus the tumor microenvironment (TME).

The marked differences highlight a profound phenotypic reprogramming of fibroblasts in the presence of tumor, indicative of their active participation in shaping the TME.

Clinical or Translational Implications

The identification of condition-specific surfaceome markers for fibroblasts holds significant clinical and translational potential, particularly for targeting CAFs in colorectal cancer.

Targeting these specific surface proteins on CAFs could offer strategies to modulate the tumor microenvironment, enhance immune responses, reduce desmoplasia, and improve the efficacy of conventional or immunotherapeutic treatments for colon cancer. Further validation in preclinical models and clinical studies would be crucial.

18. T cell CD4+ Condition-Specific Surfaceome Markers in Colon Tissue

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis identifies and visualizes condition-specific surfaceome markers for CD4+ T cells in human colon tissue, comparing Tumor (colorectal cancer) and Adjacent Normal conditions across multiple patient samples. The dot plot displays the mean expression level (color intensity) and the fraction of cells expressing each marker (dot size) for a selected panel of surface proteins. This provides insights into the distinct surface phenotypes of CD4+ T cells in different tissue microenvironments.

Visual Summary

The dot plot clearly delineates two major groups of surface markers, distinguishing CD4+ T cells from Adjacent Normal (Adj_normal) tissue from those in Tumor tissue.

Biological Interpretation

The distinct surfaceome profiles reveal significant phenotypic changes in CD4+ T cells residing within the colon tumor microenvironment compared to adjacent normal tissue.

Modulation of T cell Co-stimulation and Co-inhibition:

Clinical or Translational Implications

The identified condition-specific surfaceome markers have significant clinical and translational implications, particularly for biomarker discovery and therapeutic targeting in colorectal cancer.

19. Intestinal Epithelial Cell Cycle Genes Are Significantly Upregulated in Colon Tumor Microenvironment

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates the expression of a curated set of cell cycle pathway genes in Intestinal Epithelial cells, comparing colon tumor tissue (Tumor) with adjacent normal tissue (Adj_normal). The tool plot_box_for_gene_expression_with_signif_difference was used to identify and visualize genes with statistically significant differences in their "expressing cell fraction (sample)" between these two conditions. The "expressing cell fraction (sample)" represents the proportion of cells within each sample that express a given gene. The Intestinal Epithelial cell population is the designated cell type of tumor origin in this dataset, making this analysis highly relevant to understanding tumor biology.

Visual Summary

The visualization displays box plots for 24 distinct cell cycle-related genes. A striking and consistent pattern is observed across all plotted genes:

Biological Interpretation

The observed widespread and highly significant upregulation of cell cycle pathway genes in Intestinal Epithelial cells within colon tumors strongly points to uncontrolled proliferation as a key hallmark of these malignant cells. Given that Intestinal Epithelial cells are identified as the tumor origin cell type, these findings are directly reflective of the transformed state of these cells.

Overall, these findings provide strong evidence that Intestinal Epithelial cells in colon tumors are characterized by a highly active and dysregulated cell cycle machinery, driving tumor growth and expansion.

Clinical or Translational Implications

The pervasive upregulation of cell cycle genes in Intestinal Epithelial cells from colon tumors has significant clinical and translational implications:

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

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents Gene Ontology (GO) enrichment results (Gene Set Analysis, GSA_up) for Intestinal Epithelial cells, comparing three distinct cellular states:

  1. Adj_normal_vs_others: Pathways upregulated in Intestinal Epithelial cells from "Adjacent Normal" tissue compared to all other Intestinal Epithelial cells (primarily those from "Tumor" tissue).
  2. Diploid_vs_others: Pathways upregulated in Intestinal Epithelial cells inferred to be "Diploid" compared to those inferred as "Aneuploid."
  3. Tumor_vs_others: Pathways upregulated in Intestinal Epithelial cells from "Tumor" tissue compared to all other Intestinal Epithelial cells (primarily those from "Adjacent Normal" tissue).

The results are displayed as bar plots, where the length of the bar corresponds to the statistical significance (-log(p-val) and -log(q-val)) of the enrichment.

Visual Summary

The provided bar plots highlight distinct functional differences across the analyzed conditions and ploidy states within Intestinal Epithelial cells.

Biological Interpretation

These GSA results provide a comprehensive biological understanding of Intestinal Epithelial cells in the context of colorectal cancer, considering both tissue origin and ploidy status.

  1. Metabolic Reprogramming in Normal vs. Tumor Epithelium: Adjacent normal Intestinal Epithelial cells exhibit high metabolic activity, particularly in lipid and amino acid catabolism, essential for maintaining gut homeostasis and energy supply. The PPAR signaling pathway, crucial for lipid metabolism and anti-inflammatory responses, is highly active [1]. In contrast, tumor cells display a shift towards pathways related to increased protein synthesis and processing (ER protein processing, ribosome, spliceosome) and RNA metabolism (RNA transport, RNA degradation), indicative of high proliferative demand and stress within the malignant cells. This metabolic shift is a hallmark of cancer, where cells prioritize anabolism to support rapid growth [2].
  2. Immune Surveillance and Ploidy: Diploid Intestinal Epithelial cells show a robust upregulation of immune-related pathways, including Toll-like receptor signaling and NF-kappa B signaling, which are critical for recognizing pathogens and initiating inflammatory responses [3]. The enrichment of pathways related to various infections (viral, bacterial) suggests that diploid cells, potentially representing a healthier or more immune-competent subpopulation, are actively involved in defending against pathogens, a constant challenge in the gut environment. This contrasts with aneuploid cells (implied by 'others'), which might have compromised immune functions or altered cellular stress responses.
  3. Hallmarks of Colorectal Cancer: The enrichment in tumor Intestinal Epithelial cells directly reflects established hallmarks of cancer. The activation of "Cell cycle" pathways underscores uncontrolled proliferation. "Protein processing in endoplasmic reticulum" and "Ubiquitin mediated proteolysis" highlight significant cellular stress and increased protein turnover, often associated with rapid growth and altered protein quality control in cancer. The appearance of "Colorectal cancer" and "Pathways in cancer" terms explicitly confirms the cancerous state. The concurrent upregulation of "Autophagy" and "Cellular senescence" pathways suggests complex adaptive and maladaptive stress responses within the tumor cells, potentially promoting survival or contributing to resistance [4]. Changes in "Adherens junction" integrity are crucial for epithelial-mesenchymal transition (EMT) and metastasis [5].

Clinical or Translational Implications

---

References:

[1] PPAR signaling pathway. *KEGG Pathway Database*. https://www.genome.jp/pathway/hsa03320

[2] Hanahan, D., & Weinberg, R. A. (2011). Hallmarks of cancer: the next generation. *Cell*, 144(5), 646–674. https://pubmed.ncbi.nlm.nih.gov/21376720/

[3] Toll-like receptor signaling pathway. *KEGG Pathway Database*. https://www.genome.jp/pathway/hsa04620

[4] Levy, J. M., Towers, C. G., & Thorburn, A. (2017). Targeting autophagy in cancer. *Nature Reviews Cancer*, 17(9), 528–542. https://pubmed.ncbi.nlm.nih.gov/28775439/

[5] Adherens junction. *KEGG Pathway Database*. https://www.genome.jp/pathway/hsa04520

21. Gene Set Enrichment Analysis Reveals Condition- and Cell Type-Specific Pathway Alterations in Colon Tissue

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents Gene Set Enrichment Analysis (GSEA) results for key cell types found in colon tissue, comparing Tumor and Adjacent Normal (Adj_normal) conditions. For Intestinal Epithelial Cells (IECs), the analysis further distinguishes between diploid and aneuploid cell populations. The dot plot visualizes the Normalized Enrichment Score (NES) and statistical significance (-log(p-val)) for 120 selected gene sets, providing insights into the biological pathways differentially activated or suppressed in various cell types under different disease and ploidy states.

Visual Summary

The dot plot effectively summarizes the GSEA results, with each dot representing a gene set (y-axis) enriched in a specific cell type and comparison context (x-axis).

Overall, the plot reveals widespread and distinct pathway perturbations across different cell types and conditions, with particularly strong and numerous signals observed in Intestinal Epithelial Cells, immune cells (Macrophages, T cells, B cells, Plasma cells, ILCs), and stromal cells (Fibroblasts, Endothelial cells). The inclusion of ploidy-specific comparisons for IECs highlights specific metabolic and proliferative differences.

Biological Interpretation

  1. Intestinal Epithelial Cells (IECs) Undergoing Malignant Transformation:
  1. Tumor Microenvironment and Immune Cell Modulation:
  1. Stromal Contributions to Tumor Progression:

Clinical or Translational Implications

  1. Therapeutic Targeting of Wnt and HIF-1 Pathways: The consistent upregulation of Wnt signaling pathway in tumor IECs, fibroblasts, macrophages, and T cells, and HIF-1 signaling pathway across multiple cell types (IECs, macrophages, fibroblasts, endothelial cells), suggests these are central drivers of colorectal cancer pathobiology. Targeting these pathways could offer broad therapeutic benefits, impacting not only cancer cells but also key components of the tumor microenvironment.
  2. Metabolic Reprogramming as a Therapeutic Vulnerability: The pronounced metabolic shift towards glycolysis and increased anabolism in tumor IECs (Warburg effect) highlights the potential for anti-cancer therapies that target specific metabolic enzymes or pathways.
  3. Immuno-oncology Strategies: The widespread activation and differentiation of various immune cells (T cells, B cells, Plasma cells, Macrophages, ILCs) and the enrichment of inflammatory pathways underscore the dynamic immune landscape in colorectal cancer. Understanding the precise pro- or anti-tumorigenic roles of these activated immune cell subsets is crucial for developing effective immunotherapies, such as checkpoint inhibitors or adoptive cell therapies. The Inflammatory bowel disease pathway enrichment in T cells suggests a chronic inflammatory component that could be targeted to modulate immune responses.
  4. Aneuploidy-Specific Therapeutic Approaches: The distinct activation of DNA replication, cell cycle, and DNA repair pathways, coupled with p53 downregulation in aneuploid tumor IECs, suggests that these cells might be particularly susceptible to therapies that disrupt cell cycle progression or DNA integrity. Conversely, they might exhibit resistance to conventional treatments due to enhanced repair mechanisms. This highlights the importance of ploidy status in guiding precision medicine strategies.
  5. Targeting Stromal-Tumor Interactions: The activation of CAFs and endothelial cells and their involvement in ECM remodeling and angiogenesis (indicated by HIF-1, Wnt, and ECM-related pathways) emphasizes the importance of targeting the tumor microenvironment. Anti-angiogenic therapies or agents that modulate CAF function could impede tumor growth and metastasis.

22. Discussion

The comprehensive single-cell analysis of human colon tissue reveals a profoundly altered cellular and molecular landscape in colon cancer compared to adjacent normal tissue. A central finding is the clear identification of malignant Intestinal Epithelial cells, the tumor origin cell type, which are largely defined by pervasive aneuploidy and extensive genomic instability, as evidenced by widespread copy number variations (CNVs). Notably, amplification of the *EGFR*-containing region (7p14.1:7q11.23) is a recurrent event in these tumor cells, underscoring a key oncogenic driver. The distinct clustering of aneuploid, epithelial, and tumor-derived cells in CNV-based UMAPs further validates their malignant identity and highlights intratumoral heterogeneity, with some tumor samples containing significant diploid epithelial populations, possibly representing less transformed or reactive cells.

Beyond the malignant epithelium, the tumor microenvironment (TME) undergoes substantial reprogramming. Immune cell populations exhibit significant shifts, notably a pronounced increase in immunosuppressive T regulatory (Treg) cells and Th17 cells, alongside a depletion of anti-tumorigenic ILC1 and ILC2 populations. This suggests an immune landscape skewed towards immune evasion and chronic inflammation. Macrophages within the TME also display altered polarization, with a significant increase in pro-tumorigenic M2B macrophages and a decrease in M2A macrophages, indicating their active contribution to an immunosuppressive and angiogenesis-promoting environment. Cancer-associated fibroblasts (CAFs), identified by markers such as FAP and PDGFRB, are highly activated in tumor tissue, driving extensive extracellular matrix remodeling through upregulated collagen-integrin interactions.

Cell-cell interaction analysis further elucidates the complex communication network within the TME. A shift from epithelial-stromal homeostatic interactions in normal tissue to dominant interactions between aneuploid epithelial cells, macrophages, and activated fibroblasts is observed in tumors. Prominently, EREG-EGFR interactions mediate communication between tumor cells and immune cells, promoting growth. Pro-inflammatory and pro-angiogenic signals are amplified, including IL6-IL6R, VEGFA-VEGFR1/VEGFR2, and Ephrin-Eph receptor signaling, particularly involving macrophages and tumor epithelial cells. Crucial immune checkpoint and co-stimulatory axes, such as CD86-CD28, are highly active, suggesting ongoing but potentially exhausted immune responses. Upregulation of Prostaglandin E2, CXCL12-CXCR4, and TGFB1-TGFbeta_receptor1 interactions within the TME further points to key pathways driving immunosuppression and stromal remodeling.

Gene expression and pathway analyses underscore these cellular and microenvironmental changes. Tumor Intestinal Epithelial cells exhibit pervasive and highly significant upregulation of cell cycle genes (e.g., CCND1, MYC), confirming uncontrolled proliferation. GSEA results highlight metabolic reprogramming in tumor IECs towards glycolysis (Warburg effect) and away from oxidative phosphorylation, along with activated Wnt and HIF-1 signaling pathways. Aneuploid tumor IECs specifically show activation of DNA replication and repair pathways alongside downregulation of p53 signaling, indicating both genomic instability and compromised tumor suppression. Across the TME, Wnt and HIF-1 signaling are consistently upregulated in various cell types (fibroblasts, macrophages, endothelial cells), suggesting their central role in shaping the tumor. These findings collectively paint a picture of highly coordinated oncogenic processes impacting tumor cells and their supportive microenvironment in colon cancer.

Hypotheses:

  1. Aneuploid Intestinal Epithelial cells in colon tumors, characterized by recurrent *EGFR* amplification and pervasive cell cycle dysregulation, are the primary drivers of tumor growth and metastatic potential, distinct from diploid epithelial cells.
  2. The profound immune dysregulation in the colon TME, marked by increased Treg cells and M2B macrophages and decreased ILC1 cells, directly contributes to immune evasion by suppressing anti-tumor effector functions and promoting chronic inflammation.
  3. Reprogrammed cancer-associated fibroblasts (CAFs) actively remodel the extracellular matrix and engage in pro-tumorigenic cell-cell interactions (e.g., via CXCL12-CXCR4, TGFB1-TGFbeta_receptor1, FAP) that recruit immunosuppressive cells and foster tumor cell invasion.
  4. The shift in EGFR ligand usage from AREG/HBEGF in normal epithelium to EREG in tumor cells, coupled with its interaction with immune cells, represents a mechanism by which tumor cells actively manipulate the microenvironment for their survival and proliferation.
  5. Metabolic reprogramming, including heightened glycolysis (Warburg effect) and cholesterol/purine/pyrimidine metabolism, is a fundamental vulnerability of colon cancer cells, particularly aneuploid Intestinal Epithelial cells, supporting their rapid proliferation and biomass accumulation.

Potential therapeutic targets:

  1. EGFR (Epidermal Growth Factor Receptor): Recurrent amplification of the *EGFR*-containing region (7p14.1:7q11.23) in tumor Intestinal Epithelial cells, along with upregulated expression of its ligand EREG and activation of EREG-EGFR interactions between tumor and immune cells. EGFR signaling is a known driver of proliferation and survival in colorectal cancer. Evidence: CNV analysis (Section 4), condition-specific surfaceome markers for Intestinal Epithelial cells (MET, EREG are co-upregulated) (Section 15), and immune checkpoint/cell cycle focused CCI (EREG-EGFR interactions between Aneuploid Intestinal Epi and T CD8+/Mac) (Section 13). Validation: Test existing EGFR inhibitors (e.g., cetuximab, panitumumab) or novel EREG-specific blocking agents in preclinical models or in patient cohorts selected for EGFR amplification or high EREG expression; monitor tumor cell proliferation and immune modulation.
  2. TACSTD2 (TROP2): Highly specific and significantly upregulated surface marker on tumor Intestinal Epithelial cells, particularly in aneuploid populations. TROP2 is a well-established target for antibody-drug conjugates (ADCs) in epithelial cancers. Evidence: Condition-specific surfaceome markers for Intestinal Epithelial cells show strong, pervasive expression of TACSTD2 in tumor samples (Section 15). Validation: Evaluate the efficacy of TROP2-targeting ADCs (e.g., Sacituzumab govitecan) in colon cancer models or in clinical trials, especially in patients with high TROP2 expression in their tumors.
  3. FAP (Fibroblast Activation Protein): FAP is a canonical and highly specific surface marker for cancer-associated fibroblasts (CAFs) in colon tumors. CAFs drive ECM remodeling, immune suppression, and tumor growth, representing a critical component of the pro-tumorigenic microenvironment. Evidence: Condition-specific surfaceome markers for Fibroblasts show strong upregulation of FAP in tumor samples (Section 17). Increased COL-integrin interactions in tumor suggest active CAF remodeling (Section 14). Validation: Develop or test FAP-targeting agents (e.g., FAP-specific ADCs, FAP-CAR-T cells, or FAP inhibitors) to deplete or reprogram CAFs in colon cancer models, assessing impact on tumor growth, metastasis, and immune infiltration/function.
  4. Treg cells / CTLA4: Treg cells are significantly enriched in the tumor microenvironment and are potent immunosuppressors, contributing to immune evasion. CTLA4 is a key immune checkpoint receptor expressed on Tregs and activated T cells, mediating their suppressive function. Evidence: T cell subset analysis shows significantly higher Treg proportion in Tumor (Section 7, 8). Condition-specific surfaceome markers for T cell CD4+ show upregulation of CTLA4 in tumor (Section 18). Validation: Test CTLA4 blocking antibodies (e.g., ipilimumab) alone or in combination with other immunotherapies in colon cancer patients, evaluating reduction in Treg activity and enhancement of anti-tumor effector T cell responses.
  5. M2B Macrophages / MMP14 / TREM2: M2B macrophages are significantly increased in the tumor microenvironment and are implicated in immune suppression, angiogenesis, and tumor progression. MMP14 and TREM2 are highly expressed surface markers on tumor macrophages, mediating ECM remodeling and TAM survival/immunosuppression, respectively. Evidence: Macrophage subset analysis shows significantly higher M2B proportion in Tumor (Section 9, 10). Condition-specific surfaceome markers for Macrophages show high MMP14 and TREM2 expression in tumor macrophages (Section 16). Validation: Explore strategies to reprogram M2B macrophages or specifically target MMP14 or TREM2 activity (e.g., using small molecule inhibitors or blocking antibodies) in preclinical colon cancer models, assessing impact on tumor growth, angiogenesis, and immune context.

Follow-up validation ideas:

  1. Targeted qPCR/Immunostaining/Flow Cytometry: Validate the differential expression of key surface markers (e.g., TROP2, MET, CEACAM6 on tumor IECs; FAP, PDGFRB, B7-H3 on CAFs; MMP14, TREM2 on TAMs; TIGIT, CTLA4, OX40 on T cells) in larger patient cohorts using bulk tissue, flow cytometry of dissociated cells, or spatial transcriptomics/immunostaining to confirm their cell-type and condition specificity and explore prognostic value.
  2. In vitro Perturbation Assays: Use patient-derived organoids or cell lines to functionally investigate the impact of perturbing key pathways (e.g., EGFR signaling using EREG, Wnt signaling, HIF-1 signaling) on cell proliferation, survival, and differentiation in aneuploid Intestinal Epithelial cells, especially in 3D culture models that mimic tissue architecture.
  3. In vivo Efficacy Studies: Test the therapeutic potential of targeting identified pathways or surface markers (e.g., TROP2 ADCs, FAP inhibitors, B7-H3 blocking antibodies, or agents modulating M2B macrophage polarization or Treg function) in patient-derived xenograft (PDX) models or syngeneic mouse models of colon cancer.
  4. Functional Immune Assays: Conduct co-culture experiments with tumor cells, CAFs, and specific immune cell subsets from colon cancer patients to validate critical cell-cell interactions (e.g., CD86-CD28, CXCL12-CXCR4, TGFB1-TGFbeta_receptor1, IL6-IL6R, Ephrin-Eph) and their effects on immune cell activation, differentiation, and anti-tumor effector functions.
  5. CNV-FISH/ArrayCGH Validation: Perform Fluorescence In Situ Hybridization (FISH) or array Comparative Genomic Hybridization (aCGH) on macro-dissected tumor regions or single cells to confirm recurrent genomic amplifications (e.g., *EGFR*) and validate aneuploidy status in a larger cohort.
  6. Metabolic Flux Analysis: Use Seahorse Analyzer or stable isotope tracing in aneuploid Intestinal Epithelial cell lines or organoids to quantify metabolic shifts (e.g., glycolysis, oxidative phosphorylation rates) and confirm the Warburg effect and other metabolic vulnerabilities.

Limitations:

This single-cell RNA sequencing analysis provides a high-resolution view of cellular and molecular changes in colon cancer. However, it is based on a snapshot of gene expression and does not fully capture dynamic processes, protein-level modifications, or the spatial organization of cells, which are crucial aspects of tumor biology. While computational methods infer copy number variations and cell-cell interactions, these require orthogonal experimental validation. The observed associations between molecular profiles and disease state or ploidy status are correlative, and functional causality needs to be established through further experimental studies. The patient cohort size, while providing robust trends, might not capture the full spectrum of inter-patient heterogeneity in colorectal cancer.

23. Query List

  1. Show and save UMAPs for condition, sample, major cell type, minor cell type, ploidy_dec, and celltype_subset in two columns.
  2. Show major cell type scores on UMAP and save the result.
  3. Show and save a marker expression dot plot for celltype_subset. SET target_cell = None and var_group_rotation = 45. Leave all other arguments at their default values.
  4. Select Intestinal Epithelial cells, show a CNV heatmap grouped by sample together with a summary of regions with significantly amplified copy numbers, and save the results.
  5. Show and save UMAPs of CNV patterns colored by major cell type, minor cell type, ploidy results, condition, and sample in two columns.
  6. Show and save a population bar plot of minor cell types.
  7. Show and save a population bar plot of T cell subsets.
  8. Show and save box plots of T cell subset populations with statistically significant differences between conditions, if any. Choose ncols based on the total number of panels.
  9. Show and save a population bar plot of macrophage subsets.
  10. Show and save box plots of macrophage subset populations with statistically significant differences between conditions, if any. Choose ncols based on the total number of panels.
  11. Select Intestinal Epithelial cells, show their ploidy populations as a bar plot, and save the result.
  12. Show and save cell-cell interaction patterns involving Intestinal Epithelial cells, Fibroblast, Macrophage, T cell CD4+, and T cell CD8+. Select at most 80 cell-cell interactions per group.
  13. Select only genes associated with immune checkpoint and cell cycle pathways, show cell-cell interactions for these genes, and save the result.
  14. Find cell-cell interactions involving major immune and stromal cells that differ significantly between conditions, show them as a dot plot, and save the result. Set max_n_items_per_group to 60.
  15. Extract condition-specific markers for Intestinal Epithelial cells, show their expression as a dot plot, and save the result. Use only surfaceome markers, with at most 50 markers per condition.
  16. Extract condition-specific markers for Macrophage to show the markers expression as a dot plot, and save the result. Use only surfaceome markers, with at most 30 markers per condition.
  17. Extract condition-specific markers for Fibroblast to show the markers expression as a dot plot, and save the result. Use only surfaceome markers, with at most 30 markers per condition.
  18. Extract condition-specific markers for T cell CD4+ to show the markers expression as a dot plot, and save the result. Use only surfaceome markers, with at most 30 markers per condition.
  19. Among cell cycle pathway genes, select those with statistically significant expression differences between conditions in Intestinal Epithelial cells, show box plots, and save the result. Set max_n_items_to_plot to 24 and choose ncols based on the total number of panels for an overall width-to-height ratio of approximately 2:3.
  20. Show and save Gene Ontology (GSA) analysis results for Intestinal Epithelial cells as a bar plot.
  21. Show Gene Set Enrichment Analysis results as a dot plot, and save it. Set the color map to RdBu_r and n_pws_to_show to 120.
↑ Top