SCODiA Report by MLBI Lab

Single-Cell Transcriptomics Reveals Colorectal Cancer Microenvironment Remodeling and Immunosuppression

This report details single-cell RNA-seq findings from human colon tissue, highlighting profound cellular and molecular changes in colorectal cancer. Key observations include the malignant transformation of Intestinal Epithelial cells, characterized by aneuploidy and dysregulated cell cycle, along with significant shifts in immune cell populations and extensive stromal reprogramming. The tumor microenvironment is marked by increased immunosuppressive T cell and macrophage subsets, altered cell-cell communication, and a strong pro-tumorigenic fibroblast phenotype. These insights uncover potential biomarkers and therapeutic targets for colorectal cancer.

Contents

  1. Dataset overview
  2. UMAP Visualization of Cell Populations by Condition, Sample, Cell Type, and Ploidy in Colon Tissue
  3. Major Cell Type Score and Ploidy Distribution on UMAP
  4. Celltype Subset Marker Expression Validation
  5. Intestinal Epithelial Cells (Tumor Origin) CNV Heatmap and Amplification Summary
  6. CNV-based UMAP Analysis of Colon Tissue Cells
  7. Minor Cell Type Population Analysis in Colon Normal vs. Tumor Tissue
  8. T cell Subsets Population Analysis in Normal and Colon Tumor Tissues
  9. T cell Subset Population Shifts in Colon Cancer
  10. Macrophage Subpopulation Shifts in Colon Cancer
  11. Macrophage Subset Population Shifts in Colon Tumor Microenvironment
  12. Intestinal Epithelial Cell Ploidy Distribution in Normal vs. Tumor Conditions
  13. Colon Cell-Cell Interaction Patterns in Normal and Tumor Conditions
  14. Immune Checkpoint and Cell Cycle Gene-Associated Cell-Cell Interactions in Colon Normal vs. Tumor Microenvironments
  15. Condition-Specific Cell-Cell Interaction Patterns in Normal and Tumor Colon Tissue
  16. Intestinal Epithelial Cell의 조건 특이적 표면 마커 발현 패턴 분석
  17. Surfaceome Marker Characterization of Colon Cell Subtypes
  18. Fibroblast Condition-Specific Surfaceome Markers in Colon Tissue
  19. T cell CD4+ Condition-Specific Surfaceome Markers in Colon Tissue
  20. Dysregulation of Cell Cycle Pathway Genes in Intestinal Epithelial Cells of Colon Tumors
  21. Gene Ontology (GSA) Analysis of Intestinal Epithelial Cells in Colon Tissue
  22. Gene Set Enrichment Analysis (GSEA) of Intestinal Epithelial Cells, CD4+ T Cells, and Fibroblasts in Colon Tissue
  23. Discussion
  24. Query List

0. Dataset overview

Dataset Summary

Precomputed Results

1. UMAP Visualization of Cell Populations by Condition, Sample, Cell Type, and Ploidy in Colon Tissue

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents UMAP visualizations of single-cell RNA-seq data from human colon tissue, dissecting cellular heterogeneity across various biological and technical annotations. The UMAP plots allow for an assessment of the overall structure of the dataset, the quality of cell type annotations, the presence of batch effects, and the distribution of cells from normal and tumor conditions, including those with aneuploid status.

Visual Summary

The UMAPs display the relationships between 48,033 cells based on their transcriptional profiles, colored by six different metadata categories: condition, sample, celltype_major, celltype_minor, ploidy_dec, and celltype_subset.

Biological Interpretation

  1. Robust Cell Type Identification: The clear and hierarchical segregation of cell populations from major to subset levels confirms a high-quality annotation strategy and robust clustering. The distinct clustering of celltype_major, celltype_minor, and celltype_subset demonstrates that the underlying gene expression patterns effectively distinguish diverse cell identities present in the colon tissue. The relatively low number of 'unassigned' cells further supports the comprehensiveness of the cell type annotations.
  2. Tumor-Specific Cell Compartments: The condition UMAP reveals populations predominantly associated with the 'tumor' state. The most striking observation is the large Intestinal Epithelial cell cluster (top right) being highly enriched for tumor cells. This is consistent with the Intestinal Epithelial cell being the designated 'Tumor origin celltype', implying these are the transformed malignant cells.
  3. Aneuploidy as a Hallmark of Cancer Cells: The ploidy_dec UMAP provides strong evidence for the identity of tumor cells. The 'Aneuploid' cells almost exclusively reside within the Intestinal Epithelial cell cluster, which itself is highly enriched in cells from 'tumor' conditions. Aneuploidy, or an abnormal number of chromosomes, is a well-established hallmark of cancer cells, confirming that these epithelial cells are indeed the malignant population within the tumor microenvironment [GeneCards - Aneuploidy: GeneCards].
  4. Heterogeneity within Tumor Epithelium and Microenvironment: Even within the putatively malignant Intestinal Epithelial cell cluster, the celltype_subset map indicates diverse epithelial subtypes (e.g., Goblet, Paneth, Enterocyte, Crypt cells). This suggests either the presence of distinct malignant epithelial subclones with varying differentiation states or the persistence of non-transformed epithelial cells within the tumor context. Furthermore, the diverse immune and stromal populations observed (T cells, B cells, Myeloid cells, Fibroblasts, Endothelial cells, etc.) indicate a complex and active tumor microenvironment, where different immune and stromal cell subsets are likely interacting with the tumor cells and influencing disease progression.

Annotation Notes

The comprehensive and well-segregated UMAPs for all levels of cell type annotation, coupled with the clear distinction of normal/tumor conditions and the biological validation of aneuploidy in tumor epithelial cells, indicate a high quality of the single-cell dataset and its annotations. The effective integration, as evidenced by sample mixing, ensures that downstream differential expression and cell-cell interaction analyses will likely reflect true biological variations rather than technical artifacts.

2. Major Cell Type Score and Ploidy Distribution on UMAP

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents Uniform Manifold Approximation and Projection (UMAP) visualizations, illustrating the distribution of various major cell type scores, inferred cellular ploidy status, and the final assigned major cell types across a single-cell RNA-seq dataset from human colon tissue. These plots collectively provide a comprehensive overview of the cellular landscape, the confidence of cell type annotations, and the genomic stability (ploidy) of different cell populations.

Visual Summary

The UMAP projections reveal a well-structured embedding with distinct clusters representing various cell populations.

Biological Interpretation

The UMAP visualizations provide critical biological insights into the cellular composition and disease state of the colon tissue:

Annotation Notes

The consistency observed between the major cell type scores and the final assigned celltype_major clusters indicates a high quality of cell identity annotation in this dataset. The distinct clustering and clear separation of cell types enhance confidence in further cell-type-specific analyses. The ploidy_dec overlay serves as an excellent internal validation and provides critical context for identifying the malignant epithelial cell population, which is essential for understanding tumor biology in this dataset. Further investigation into the 'unassigned' cluster and potentially low-scoring populations like Enteric neuron might reveal novel or very rare cell types, but for major cell populations, the current annotations appear robust.

3. Celltype Subset Marker Expression Validation

Report figure

[Analysis Visualization Results]...

Analysis Overview

This dot plot visualizes the expression of selected marker genes across different celltype_subset populations derived from single-cell RNA-seq data of human colon tissue. Each row represents a celltype_subset, and each column represents a marker gene. The size of the dot corresponds to the fraction of cells within that group expressing the gene, while the color intensity reflects the mean expression level of the gene in those cells. This analysis aims to confirm the identity and specificity of the assigned celltype_subset annotations based on their characteristic gene expression profiles.

Visual Summary

The dot plot displays a clear diagonal pattern of strong marker gene expression, indicating that most celltype_subset annotations are well-defined by distinct sets of genes. Red boxes highlight these cell-type-specific marker clusters.

Key observations include:

Biological Interpretation

The observed marker expression patterns largely align with known biological identities of these cell types in the human colon, providing strong evidence for the accuracy of the celltype_subset annotations.

  1. B Lymphocytes and Plasma Cells:
  1. Intestinal Epithelial Cells (Tumor Origin):
  1. Myeloid Cells:
  1. T Lymphocytes:
  1. Stromal and Endothelial Cells:

Annotation Notes

The marker expression dot plot strongly supports the current celltype_subset annotations. The distinct, specific expression patterns of known lineage markers for almost all cell types provide high confidence in the quality of the cell type assignments. The plot_markers_and_expression_dot tool successfully identified highly specific surfaceome markers (surfaceome_only: True) for each group, which is particularly valuable for cell identity validation. No major misannotations or ambiguous populations are immediately apparent from this visualization, suggesting robust and reliable cell type classification at the subset level. This foundational annotation work is critical for subsequent downstream analyses, such as differential gene expression or cell-cell interaction studies.

4. Intestinal Epithelial Cells (Tumor Origin) CNV Heatmap and Amplification Summary

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis aimed to visualize copy number variations (CNVs) in Intestinal Epithelial cells, designated as the tumor-origin cell type, grouped by sample. The visualization provides a heatmap of log2(Copy Number Ratio, CNR) across chromosomes for individual cell groups and a summary of frequently amplified cytogenetic regions.

Visual Summary

  1. CNV Heatmap (Image 1):
  1. CNV Amplification Summary (Image 2):

Biological Interpretation

The user query specifically targeted "Intestinal Epithelial cell" cells as the tumor origin. However, the provided heatmap and summary primarily display CNV patterns in cell groups labeled as "B cell" and "T cell" populations. This discrepancy is critical for interpretation:

Potential Explanations

  1. Cell Misclassification: The T_cacX groups might represent misclassified Intestinal Epithelial tumor cells that exhibit an altered gene expression profile, leading to their erroneous assignment as T cells. This is a common challenge in single-cell analysis of highly plastic tumor cells.
  2. Malignant Lymphoid Cells: Less likely in the context of a primary colon tumor analysis, these T-cell groups could represent a concurrent, potentially malignant, lymphoid population (e.g., lymphoma or leukemia within the tumor microenvironment).
  3. Technical Artifact: There could be technical issues in cell type annotation or CNV calling within these specific samples.

Clinical or Translational Implications

Given the interpretation that the observed CNVs, particularly the recurrent amplifications, are likely indicative of a malignant process (either within misclassified Intestinal Epithelial cells or an unexpected malignancy in the T-cell compartment):

Annotation Notes

The primary observation from this analysis is the presence of significant copy number amplifications in cell groups labeled as T cells (T_cac1, T_cac3, T_cac8, T_cac9), despite the user's query specifying "Intestinal Epithelial cell" as the tumor-origin cell type. The visualization provides a clear distinction between these aberrantly copy-numbered cell groups and other "Diploid B/T" cell groups that maintain genomic stability. This highlights a potential issue in cell type annotation for the CNV-positive cells or an unexpected biological finding that requires further investigation into cell identity and tumor heterogeneity.

5. CNV-based UMAP Analysis of Colon Tissue Cells

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents UMAP visualizations of single-cell RNA-seq data from human colon tissue, leveraging Copy Number Variation (CNV) estimates to define the embedding space. The UMAPs are colored by major cell type, minor cell type, ploidy status, condition (normal/tumor), and individual sample, providing an integrated view of cellular heterogeneity based on genomic alterations. The primary goal is to understand how CNV patterns differentiate cell populations and relate to disease status and cell identity.

Visual Summary

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

Biological Interpretation

The CNV-based UMAP effectively delineates cell populations based on their genomic integrity, providing strong biological insights into the dataset:

Annotation Notes

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

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents a stacked bar plot showing the proportional representation of various minor cell types across individual samples from both normal colon tissue and colon tumors. The purpose is to identify shifts in cellular composition associated with the disease state, providing an overview of the tumor microenvironment (TME) and its differences from healthy tissue.

Visual Summary

The visualization displays two main panels, one for "normal" samples and one for "tumor" samples, each comprising multiple individual sample bars. Each bar is stacked with different colored segments representing the proportion of various minor cell types.

Key observations from the plot include:

Biological Interpretation

The observed cellular landscape provides critical insights into the pathology of colon cancer:

Immune Cell Infiltration and Exclusion:

Clinical or Translational Implications

Understanding the cellular composition of the colon tumor microenvironment has direct clinical and translational relevance:

---

References

  1. Galon, J., et al. "Type, Density, and Location of Immune Cells within Human Colorectal Tumors Predict Clinical Outcome." *Science*, vol. 313, no. 5795, 2006, pp. 1960-64. PubMed Search: Colorectal cancer immune score prognosis
  2. Quail, D. F., and Joyce, J. A. "Microenvironmental Dynamics during Cancer Progression." *Nature Medicine*, vol. 19, no. 11, 2013, pp. 1423-37. PubMed Search: Tumor microenvironment stromal cells
  3. Mlecnik, B., et al. "The Immunoscore: A New Possible Approach to the Classification of Cancer." *OncoImmunology*, vol. 1, no. 6, 2012, pp. 780-87. PubMed Search: CD8+ T cells colorectal cancer prognosis
  4. Chen, D. S., and Mellman, I. "Elements of Cancer-Immunity and the Cancer-Immunity Cycle." *Immunity*, vol. 39, no. 1, 2013, pp. 1-10. PubMed Search: Immune cold tumors immunotherapy
  5. Ostman, A., and Augsten, M. "Cancer-Associated Fibroblasts and Tumor Growth: Cell-Cell Interactions and Signaling Networks." *Seminars in Cancer Biology*, vol. 22, no. 4, 2012, pp. 323-31. PubMed Search: Cancer associated fibroblasts therapy
  6. Hanahan, D., and Weinberg, R. A. "Hallmarks of Cancer: The Next Generation." *Cell*, vol. 144, no. 5, 2011, pp. 646-74. PubMed Search: Tumor angiogenesis therapy

7. T cell Subsets Population Analysis in Normal and Colon Tumor Tissues

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents a population bar plot of T cell subsets, along with other innate lymphoid cells (ILCs) and NK cells, across individual normal and colon tumor tissue samples. The proportions of these immune cell populations are normalized to 100% for each sample, providing a comparative view of their relative abundance in the tumor microenvironment versus healthy tissue.

Visual Summary

The stacked bar plots display the relative proportions of various T cell subsets, NK cells, and ILCs for each sample, grouped by 'normal' and 'tumor' conditions.

Biological Interpretation

The observed shifts in T cell subset populations, particularly the increased proportion of T cell (Treg) in tumor samples, provide critical biological insights into the immune landscape of colorectal cancer.

Clinical or Translational Implications

The findings from this T cell subset population analysis have several important clinical and translational implications for colorectal cancer:

8. T cell Subset Population Shifts in Colon Cancer

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates the proportional changes of specific T cell subsets (Treg, Th22, Th17, and Tfh) between normal and tumor colon tissue conditions using single-cell RNA sequencing data. Box plots were generated to visualize the celltype proportions, and statistical significance (p-values) was calculated for the observed differences. The goal is to identify T cell populations with statistically significant differences that may play a role in the colorectal tumor microenvironment.

Visual Summary

The box plots illustrate the distribution of celltype proportions for four T cell subsets (Treg, Th22, Th17, Tfh) across normal and tumor conditions. Each black dot represents the proportion from an individual sample.

Biological Interpretation

The observed shifts in T cell subset proportions between normal and tumor colon tissue highlight significant alterations in the immune landscape associated with colorectal cancer.

  1. Increased Immunosuppressive and Pro-inflammatory T cells in Tumors:
  1. Decreased Follicular Helper T cells (Tfh) in Tumors (p=0.0925):

Collectively, these findings suggest a profound reprogramming of the T cell compartment in the colorectal tumor microenvironment, favoring immunosuppression and pro-tumorigenic inflammation, while potentially compromising efficient humoral anti-tumor responses.

Clinical or Translational Implications

The observed shifts in T cell subset populations carry several potential clinical and translational implications for colorectal cancer:

Therapeutic Targets:

9. Macrophage Subpopulation Shifts in Colon Cancer

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates the relative proportions of various macrophage subsets (M1, M2A, M2B, M2C, M2D) within single-cell RNA-seq data derived from human colon tissue. The stacked bar plot specifically compares the macrophage composition between normal samples and tumor samples, providing insights into changes in the immune microenvironment associated with colon cancer.

Visual Summary

The stacked bar plot reveals substantial differences in macrophage subpopulation composition when comparing normal colon tissue to tumor tissue.

Biological Interpretation

Macrophages are highly adaptable immune cells that undergo diverse polarization states, playing crucial roles in both tissue homeostasis and disease, including cancer. The "M1" and "M2" classification represents functional extremes, with M1 macrophages generally pro-inflammatory and anti-tumorigenic, and M2 macrophages typically associated with immune suppression, tissue repair, and tumor promotion.

Clinical or Translational Implications

The observed shifts in macrophage subsets have important implications for understanding colon cancer biology and developing therapeutic strategies.

10. Macrophage Subset Population Shifts in Colon Tumor Microenvironment

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates the proportional representation of specific macrophage subsets (Macrophage M2A, M2B, and M2D) within the Colon tissue, comparing normal versus tumor conditions using single-cell RNA sequencing data. The goal is to identify statistically significant differences in the abundance of these subsets that may contribute to the distinct microenvironments found in normal and diseased states. The plot_box_for_celltype_population_with_signif_difference tool was used to visualize these population differences and their statistical significance.

Visual Summary

The box plots display the celltype proportion of three macrophage subsets – Mac (M2A), Mac (M2B), and Mac (M2D) – across normal and tumor conditions. Each plot includes individual data points (samples) overlaid on the box plot, along with p-values indicating statistical significance of differences between conditions.

Biological Interpretation

Macrophages are highly plastic immune cells that polarize into distinct functional phenotypes based on microenvironmental cues. M2-like macrophages, often termed tumor-associated macrophages (TAMs), are broadly associated with immunosuppression, angiogenesis, tissue remodeling, and tumor progression PubMed search: Tumor-associated macrophages M2 colorectal cancer. The observed shifts in macrophage subsets provide insights into the altered immune landscape in colorectal cancer.

Collectively, these findings indicate a significant re-programming of the macrophage compartment within the colon tumor microenvironment. While M2A macrophages decrease, the pro-tumorigenic M2B and M2D subsets become more abundant, underscoring a shift towards an immunosuppressive and pro-angiogenic phenotype that favors tumor development and progression.

Clinical or Translational Implications

The observed shifts in macrophage subsets highlight their potential as diagnostic biomarkers and therapeutic targets in colorectal cancer.

Further research into the specific molecular mechanisms driving these macrophage polarization shifts in colon cancer could uncover novel therapeutic avenues to modulate the tumor microenvironment and improve patient outcomes.

11. Intestinal Epithelial Cell Ploidy Distribution in Normal vs. Tumor Conditions

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates the ploidy status (Aneuploid, Diploid, or Unclear) of Intestinal Epithelial cells, which are identified as the tumor-origin cell type, across different individual samples under both normal and tumor conditions. The plot_celltype_population tool was used to visualize the proportion of each ploidy status as stacked bar plots for each sample, providing insights into genomic stability.

Visual Summary

The visualization presents the ploidy distribution for Intestinal Epithelial cells.

Biological Interpretation

The observed ploidy patterns provide critical biological insights into the genomic stability of Intestinal Epithelial cells in the context of colon cancer.

Clinical or Translational Implications

---

References:

  1. Aneuploidy as a Hallmark of Cancer:
  1. Aneuploidy and Prognosis in Cancer:

12. Colon Cell-Cell Interaction Patterns in Normal and Tumor Conditions

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates cell-cell interaction (CCI) patterns within the human colon, comparing normal tissue with tumor tissue using single-cell RNA sequencing data. Specifically, it focuses on interactions involving Intestinal Epithelial cells, Fibroblasts, Macrophages, and T cell subsets (CD4+ and CD8+ T cells). CellPhoneDB was used to identify ligand-receptor pairs, and the results are visualized as dot plots, highlighting the top 80 significant interactions based on p-value and mean expression.

Visual Summary

The analysis provides two dot plots: "CCI for normal" and "CCI for tumor".

CCI for Normal Tissue:

CCI for Tumor Tissue:

Biological Interpretation

The comparative analysis reveals profound shifts in cell-cell communication patterns between normal colon tissue and colorectal tumor.

  1. Loss of Epithelial-Immune Crosstalk in Tumor: A striking observation is the complete absence of interactions involving Intestinal Epithelial cells (the cell type of origin for the tumor) with T cells, or even among themselves, in the tumor microenvironment's top 80 interactions. In contrast, normal colon shows active communication between diploid intestinal epithelial cells and T cells, exemplified by interactions like HLA-E_NKG2A/E and prostaglandin pathways. This suggests a significant breakdown or suppression of direct communication between tumor cells and the immune infiltrate, potentially facilitating immune evasion. While the analysis parameter expand_ploidy_from_tumor_origin=True was set to include potentially aneuploid tumor cells as "Intestinal Epithelial cell," their interactions are not prominent in the tumor plot.
  2. Diminished Microenvironmental Diversity in Tumor: The tumor microenvironment (TME) displays a severely restricted repertoire of cell-cell interactions compared to normal tissue. The lack of Fibroblast and Macrophage involvement in the displayed top interactions, alongside the absence of Intestinal Epithelial cell interactions, indicates that the most significant communication events are concentrated within the T cell compartment itself, or that other interactions are effectively suppressed or outranked.
  3. Persistence of Intratumoral T Cell Interactions: Despite the overall reduction in diversity, T cell – T cell interactions, particularly involving CD58-CD2 and LCK-CD8-Treceptor, remain highly significant and active in the tumor.
  1. Disappearance of Immunoregulatory Pathways in Tumor: The absence of interactions such as HLA-E_CD94-NKG2A/E in the tumor is notable. In normal tissue, HLA-E presenting peptides to the inhibitory receptor NKG2A/E on T cells and NK cells can suppress immune responses [PubMed Search: HLA-E NKG2A immune evasion]. Its disappearance from the prominent interactions in tumor might suggest altered immune evasion strategies or a shift in the immune cell populations expressing these receptors. Similarly, the loss of Prostaglandin E2 related interactions, known for their immunomodulatory roles and frequent dysregulation in cancer, also points to altered immune regulation in the TME [PubMed Search: Prostaglandin E2 tumor microenvironment].

Clinical or Translational Implications

  1. Immune Evasion Mechanisms: The dramatic reduction in tumor cell-T cell interactions highlights a key mechanism of immune evasion in colorectal cancer. Therapeutic strategies could focus on re-establishing productive communication between tumor cells and immune cells, or on counteracting factors that lead to this disengagement.
  2. Targeting Intratumoral T Cell Dynamics: The persistent CD58-CD2 and LCK-CD8-Treceptor interactions within T cell infiltrates represent potential targets. While these are fundamental for T cell function, understanding their precise role (e.g., contributing to an effective vs. exhausted T cell state) in the TME is crucial. Modulating these interactions could enhance anti-tumor immunity.
  3. Biomarker Discovery: The distinct interaction profiles between normal and tumor tissue, particularly the disappearance of certain ligand-receptor pairs (e.g., HLA-E, prostaglandin-related interactions) in tumor, could serve as prognostic or predictive biomarkers for disease progression and response to immunotherapy.
  4. Re-engaging the Immune System: Identifying the reasons behind the loss of Intestinal Epithelial cell-T cell interactions in tumor could inform novel therapeutic approaches aimed at breaking immune tolerance or activating dormant anti-tumor responses. This might involve strategies to upregulate specific ligands on tumor cells or to enhance receptor expression/activity on immune cells.

13. Immune Checkpoint and Cell Cycle Gene-Associated Cell-Cell Interactions in Colon Normal vs. Tumor Microenvironments

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates cell-cell interactions (CCI) within single-cell RNA-seq data from human colon tissue, comparing normal and tumor conditions. The focus is specifically on genes associated with immune checkpoint and cell cycle pathways, as specified by the user's query. The CellPhoneDB method was utilized to identify significant ligand-receptor interactions between different cell types, and the plot_cci_dots tool visualized these interactions. The visualization presents -log10(p-value) as dot size (indicating significance) and log2(mean expression) as dot color (indicating interaction strength). The goal is to identify how these critical pathways mediate intercellular communication and if these patterns differ between healthy and cancerous states.

Visual Summary

The dot plots illustrate cell-cell interactions mediated by a predefined set of immune checkpoint and cell cycle-related genes, comparing normal and tumor conditions.

  1. Normal Condition (Left Plot):
  1. Tumor Condition (Right Plot):

In summary, the most striking difference is the disappearance of significant IFNG_Type_II_IFNR mediated interactions in the tumor microenvironment, while LCK_CD8_receptor interactions persist in both conditions but with consistently low mean expression.

Biological Interpretation

The observed differences in cell-cell interactions between normal and tumor colon tissue provide crucial insights into immune regulation and potential mechanisms of tumor immune evasion.

  1. Loss of IFNG_Type_II_IFNR Signaling in Tumor:
  1. Persistent LCK_CD8_receptor Interactions with Low Expression:

Clinical or Translational Implications

The findings have several important implications for understanding colon cancer immunology and developing therapeutic strategies:

  1. Immune Evasion Mechanism: The striking loss of significant IFNG_Type_II_IFNR interactions in the tumor microenvironment strongly suggests a mechanism of immune evasion. Tumors that can suppress or evade IFN-gamma signaling are more likely to escape immune surveillance and checkpoint inhibitor therapies that rely on robust T cell activation and IFN-gamma production.
  2. Therapeutic Target Prioritization:
  1. Experimental Validation: Further experimental validation is warranted to confirm the functional consequences of the observed CCI patterns. This could involve:
  1. Role of LCK/CD8 Interactions: While LCK_CD8_receptor interactions are present, their low mean expression warrants further investigation into the functional state of T cells in the tumor. Are these T cells anergic, exhausted, or simply less abundant in the tumor context? Modulating LCK activity or enhancing CD8-MHC-I engagement might still be relevant if T cell dysfunction is reversible.

14. Condition-Specific Cell-Cell Interaction Patterns in Normal and Tumor Colon Tissue

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates cell-cell interactions (CCIs) involving B cells, Myeloid cells, Stromal cells (Fibroblasts), T cells, Endothelial cells, and Mast cells that exhibit significant differences between normal and tumor colon tissue conditions. The results are visualized as a dot plot, where color intensity represents the scaled log strength of the interaction, and dot size indicates the statistical significance (-log10(p-value)). The plot highlights the top 60 most significantly different interactions per group, providing insights into the altered intercellular communication within the tumor microenvironment.

Visual Summary

The dot plot effectively delineates distinct patterns of cell-cell interactions that are enriched in either normal or tumor colon samples.

Biological Interpretation

The observed condition-specific CCI patterns provide critical insights into the distinct microenvironments of normal and tumor colon tissue.

Clinical or Translational Implications

The observed differences in CCI patterns between normal and tumor colon tissue hold significant clinical and translational potential.

Therapeutic Targeting:

15. Intestinal Epithelial Cell의 조건 특이적 표면 마커 발현 패턴 분석

Report figure

[Analysis Visualization Results]...

Analysis Overview

이 분석은 단일 세포 RNA 시퀀싱 데이터를 활용하여 대장(Colon) 조직의 Intestinal Epithelial cell에서 정상(normal) 및 종양(tumor) 조건에 따라 차등적으로 발현되는 표면 마커 유전자들을 식별하고 시각화합니다. 특히, surfaceome_only=True 파라미터를 사용하여 세포 표면에 존재하는 단백질만을 대상으로 하였으며, 각 조건에서 최대 50개의 유전자 마커를 선택하여 발현 패턴을 도트 플롯으로 나타냈습니다. 이는 Intestinal Epithelial cell이 종양의 기원 세포(Tumor origin celltype)임을 고려할 때, 종양 발생 및 진행과 관련된 핵심 표면 마커를 식별하는 데 중점을 둡니다.

Visual Summary

제공된 도트 플롯은 Intestinal Epithelial cell에서 정상(normal) 및 종양(tumor) 조건에 따른 유전자 발현 패턴을 시각적으로 보여줍니다.

Biological Interpretation

이 분석 결과는 대장암 발생 시 Intestinal Epithelial cell의 표면 단백질 구성에 상당한 변화가 있음을 명확히 보여줍니다.

정상 Intestinal Epithelial cell의 특징

종양 Intestinal Epithelial cell의 특징

전반적으로, 종양 특이적 표면 마커의 발현 증가는 Intestinal Epithelial cell이 종양 환경에서 증식, 침윤, 전이 및 면역 회피와 같은 악성 특징을 획득하는 분자적 메커니즘을 반영합니다.

Clinical or Translational Implications

이 분석 결과에서 식별된 조건 특이적 표면 마커들은 대장암의 진단, 예후 예측 및 치료에 중요한 임상적 의미를 가질 수 있습니다.

향후 이러한 마커들에 대한 기능적 검증 및 전임상/임상 연구를 통해 대장암 진단 및 치료 전략 개발에 기여할 수 있을 것입니다.

16. Surfaceome Marker Characterization of Colon Cell Subtypes

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis aimed to identify cell-type-specific surfaceome markers across various cell subsets found in human colon tissue, focusing particularly on Macrophage subtypes. The plot_markers_and_expression_dot tool was utilized, and parameters were set to specifically search for surfaceome markers (surfaceome_only: True). The resulting dot plot visualizes the expression fraction and mean expression levels of these identified markers across different celltype_subset populations. While the initial query requested "condition-specific markers for Macrophage," the generated plot displays markers that distinguish various cell subtypes from each other, which is crucial for detailed cell type characterization.

Visual Summary

The provided dot plot effectively displays the relative expression patterns of selected markers across a comprehensive list of colon cell subtypes (celltype_subset on the y-axis). Each row represents a distinct cell subtype, and each column corresponds to an identified marker gene.

Biological Interpretation

Macrophage Subtype Characterization

The analysis successfully identified distinct marker profiles for the various Macrophage subtypes (M1, M2A, M2B, M2C, M2D) present in the colon tissue:

Other Prominent Cell Type Markers

The plot also reveals well-known surfaceome markers for other key immune and stromal populations:

Marker Specificity

The parameter rem_mkrs_common_in_N_groups_or_more: 3 ensured that the plotted markers are relatively specific, distinguishing individual cell subtypes or small groups of related subtypes rather than being broadly expressed across many cell types. This enhances their utility for precise cell characterization.

Clinical or Translational Implications

The identification of cell-type-specific surfaceome markers, even with the inclusion of some non-surfaceome markers, has significant clinical and translational implications, particularly in the context of colon tissue and disease conditions like cancer:

Annotation Notes/Limitations

A notable observation is the presence of several intracellular proteins or transcription factors (e.g., *SOCS3, PTGS2, GATA2, LYZ, XBP1*) among the identified "surfaceome markers," despite the surfaceome_only: True parameter being applied. This suggests a potential limitation in the surfaceome annotation list utilized by the tool, or a broad definition of "surfaceome-associated" that includes genes whose expression is highly correlated with surface markers. While these non-surfaceome markers still contribute to the unique transcriptional profiles of their respective cell types, their non-surface localization should be considered if these markers are pursued for experimental validation involving surface staining or antibody-based targeting. Further validation would be required to confirm the surface localization of all genes presented as "surfaceome markers."

17. Fibroblast Condition-Specific Surfaceome Markers in Colon Tissue

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis aimed to identify surfaceome markers that distinguish fibroblasts in normal colon tissue from those in colorectal tumor tissue using single-cell RNA sequencing data. By focusing on surfaceome genes, the analysis prioritizes markers with potential for cell-surface targeting and clinical utility. The results are presented as a dot plot showing the mean expression and the fraction of cells expressing each marker across different fibroblast clusters under normal and tumor conditions.

Visual Summary

The dot plot effectively illustrates the differential expression of 30 selected surfaceome markers across various fibroblast clusters (represented as "B_cac" for normal and "T_cac" for tumor, likely indicating sample-derived clusters) under two conditions: "normal" and "tumor".

Biological Interpretation

The observed condition-specific surfaceome profiles highlight a fundamental phenotypic reprogramming of fibroblasts in the colorectal tumor microenvironment.

The diverse set of upregulated markers underscores the multifaceted roles of CAFs in promoting tumor progression in the colon, affecting aspects from ECM dynamics and angiogenesis to immune suppression.

Clinical or Translational Implications

The identification of condition-specific surfaceome markers for fibroblasts in colon cancer holds significant clinical and translational potential:

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

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis aimed to identify surfaceome markers that distinguish T cell CD4+ populations residing in normal colon tissue from those found in colon tumor tissue. The plot_markers_and_expression_dot tool was used to visualize the expression of the top condition-specific surface markers (up to 30 per condition) across different sub-clusters of T cell CD4+ cells. This provides insight into the phenotypic alterations of CD4+ T cells in the tumor microenvironment compared to normal homeostasis.

Visual Summary

The dot plot effectively illustrates distinct surface marker profiles for T cell CD4+ populations associated with normal versus tumor conditions.

Biological Interpretation

The identified surfaceome markers reveal significant biological shifts in CD4+ T cells within the colon tumor microenvironment.

Normal Homeostasis Markers

Immune Checkpoint and Dysregulation

Activation, Adhesion, and Trafficking

Other Notable Markers

The overall picture suggests that CD4+ T cells in colon tumors undergo substantial phenotypic remodeling, characterized by markers associated with chronic activation, exhaustion, immune suppression, altered migratory capacity, and potentially regulatory functions. The heterogeneity within tumor-associated clusters indicates diverse functional subsets contributing to the complex immune landscape of colorectal cancer.

Clinical or Translational Implications

The identification of condition-specific surfaceome markers for CD4+ T cells in colon cancer has several important clinical and translational implications:

19. Dysregulation of Cell Cycle Pathway Genes in Intestinal Epithelial Cells of Colon Tumors

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigated the expression levels of a predefined set of cell cycle pathway genes within Intestinal Epithelial cells, comparing normal colon tissue samples with tumor samples. The specific metric plotted is the "expressing cell fraction (sample)," which represents the proportion of cells within each sample that express a given gene. This allows for the identification of cell cycle genes whose expression patterns are significantly altered in the tumor microenvironment, specifically within the cell type identified as the tumor's origin.

Visual Summary

The box plots display the expressing cell fraction for 24 distinct cell cycle pathway genes, comparing normal (blue boxes) and tumor (orange boxes) conditions in Intestinal Epithelial cells. A consistent and statistically significant pattern emerges across all depicted genes:

Biological Interpretation

The Intestinal Epithelial cell is identified as the tumor origin celltype. The widespread and statistically significant upregulation of numerous cell cycle pathway genes in these cells within tumor samples strongly indicates a fundamental shift towards increased proliferation and uncontrolled cell division, which are hallmarks of cancer.

Overall, the elevated expressing cell fraction of these cell cycle genes in Intestinal Epithelial cells from tumors strongly supports a highly proliferative phenotype, which is a key characteristic of colon cancer.

Clinical or Translational Implications

The findings have several significant clinical and translational implications:

20. Gene Ontology (GSA) Analysis of Intestinal Epithelial Cells in Colon Tissue

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents Gene Ontology (GO) enrichment results for Intestinal Epithelial cells obtained from single-cell RNA-seq data of human colon tissue. The analysis identifies gene sets (pathways/terms) that are significantly *upregulated* in specific cellular states. Two main comparisons were performed:

  1. Diploid Intestinal Epithelial cells versus Aneuploid Intestinal Epithelial cells (Diploid_vs_others): This compares the transcriptional profiles of intestinal epithelial cells classified as diploid against those classified as aneuploid based on ploidy inference. The results show pathways enriched in diploid cells.
  2. Tumor Intestinal Epithelial cells versus other (likely Normal) Intestinal Epithelial cells (tumor_vs_others): This compares intestinal epithelial cells derived from tumor samples against those from non-tumor (normal) samples. The results highlight pathways enriched in tumor-derived epithelial cells.

Visual Summary

Intestinal Epithelial Cell: Diploid vs. Others

The bar plot for "Diploid_vs_others" shows GO terms that are significantly upregulated in Diploid Intestinal Epithelial cells compared to Aneuploid cells. The terms are sorted by their -log(p-val) (left panel) and also show -log(q-val) (right panel), indicating statistical significance.

Intestinal Epithelial Cell: Tumor vs. Others

The bar plot for "tumor_vs_others" displays GO terms that are significantly upregulated in Intestinal Epithelial cells from tumor samples compared to non-tumor cells.

Biological Interpretation

The Gene Ontology analysis of Intestinal Epithelial cells provides insight into the functional shifts associated with ploidy changes and disease state in the colon.

Insights from Diploid vs. Aneuploid Cells

The upregulation of terms like "Mineral absorption," "Aldosterone-regulated sodium reabsorption," and "Tight junction" in Diploid Intestinal Epithelial cells suggests that these cells maintain typical healthy epithelial functions, including ion transport and barrier integrity. The enrichment of key signaling pathways (FoxO, ErbB, p53, PI3K-Akt, MAPK) along with cellular quality control mechanisms such as "Cellular senescence," "Apoptosis," and "Mitophagy" indicates robust regulatory networks and active surveillance mechanisms in diploid cells. This profile is consistent with a healthy or less perturbed cellular state where growth, metabolism, and stress responses are tightly controlled, potentially preventing the accumulation of damaged or aberrant cells. The presence of "Colorectal cancer" and other cancer-related terms in this context could signify that genes involved in normal epithelial homeostasis, whose dysregulation can contribute to cancer, are appropriately regulated and active in diploid cells.

Insights from Tumor vs. Normal Cells

In contrast, Intestinal Epithelial cells from tumor samples exhibit a distinct biological signature. The striking upregulation of "Ribosome," "Protein processing in endoplasmic reticulum," "Proteasome," and "Cell cycle" pathways reflects the high metabolic demand, rapid protein synthesis, and increased proliferative activity characteristic of cancer cells. Alterations in energy metabolism, highlighted by "Oxidative phosphorylation" and "Citrate cycle (TCA cycle)," are consistent with the metabolic reprogramming observed in many cancers, often involving increased aerobic glycolysis and oxidative phosphorylation to support rapid growth [1].

The unexpected enrichment of numerous neurodegenerative disease pathways (e.g., Parkinson's, Huntington's, Alzheimer's) in tumor cells warrants careful consideration. While these diseases primarily affect the nervous system, they are often characterized by protein misfolding, aggregation, and cellular stress responses involving the ubiquitin-proteasome system, ER stress, and mitochondrial dysfunction [2]. Given the intense protein synthesis and turnover in rapidly proliferating cancer cells, these terms likely reflect a generalized cellular stress response, altered protein homeostasis, or activation of shared molecular pathways involved in managing proteotoxicity, rather than implying neurodegeneration in the colon.

Furthermore, the enrichment of viral and bacterial infection pathways suggests active host-pathogen interactions or inflammatory responses within the tumor microenvironment. Colon cancer progression is often linked to chronic inflammation and dysbiosis of the gut microbiome, which can influence epithelial cell behavior and contribute to carcinogenesis [3]. The direct upregulation of the "Colorectal cancer" pathway in tumor cells confirms the activation of specific disease-associated gene networks.

Clinical or Translational Implications

This analysis highlights key functional distinctions between healthy-like (Diploid) and disease-associated (Aneuploid/Tumor) Intestinal Epithelial cells in the colon.

References

  1. Metabolic Reprogramming in Cancer:
  1. Protein Misfolding and Stress Responses (shared mechanisms):
  1. Gut Microbiome and Colorectal Cancer:
  1. Targeting Cancer Metabolism:

21. Gene Set Enrichment Analysis (GSEA) of Intestinal Epithelial Cells, CD4+ T Cells, and Fibroblasts in Colon Tissue

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents Gene Set Enrichment Analysis (GSEA) results for key cell types in human colon tissue: Intestinal Epithelial Cells, CD4+ T cells, and Fibroblasts. The dot plot visualizes the enrichment of 120 selected pathways across different comparison groups. Each dot's color indicates the Normalized Enrichment Score (NES), where red signifies positive enrichment (upregulation of pathway genes) and blue signifies negative enrichment (downregulation). The size of each dot reflects the statistical significance, specifically -log(p-val), with larger dots indicating higher significance.

The comparisons performed are:

Intestinal Epithelial Cell:

T cell CD4+:

Fibroblast:

(Note: Although "Macrophage" was requested in the query, results for this cell type are not displayed in the provided plot based on the specified parameters.)

Visual Summary

The dot plot effectively highlights pathway enrichments and depletions across different cell types and conditions. A clear distinction is observed between normal/diploid states and tumor states, with tumor conditions generally showing stronger and broader pathway enrichments (larger, darker red dots) for proliferative, metabolic, and oncogenic signaling pathways. Conversely, pathways related to specific immune functions show differential enrichment between normal and tumor-infiltrating T cells. The use of the RdBu_r colormap clearly differentiates positively (red) and negatively (blue) enriched pathways, and varying dot sizes provide immediate insight into statistical significance.

Biological Interpretation

Intestinal Epithelial Cells (Tumor Origin)

As the tumor origin cell type, Intestinal Epithelial Cells (IECs) show profound changes in the tumor microenvironment:

CD4+ T Cells

CD4+ T cells, critical components of the adaptive immune response, display altered functional states in the tumor:

Fibroblasts

Fibroblasts in the colon show dramatic transformation in the tumor context:

Cross-Cell Type Observations

Clinical or Translational Implications

The GSEA results provide crucial insights into the biology of colon cancer and its microenvironment, with several potential clinical and translational implications:

22. Discussion

The comprehensive single-cell analysis of human colon tissue provides a detailed understanding of the cellular and molecular landscape differentiating normal physiology from colorectal cancer. A central finding is the clear identification of Intestinal Epithelial cells as the tumor origin, with a substantial proportion exhibiting aneuploidy in tumor samples, a hallmark of genomic instability and malignancy. These malignant epithelial cells demonstrate a widespread upregulation of cell cycle genes such as MYC, CCND1, CDK4, and CDK6, reflecting uncontrolled proliferation. Gene Ontology and GSEA further confirm intense metabolic reprogramming, protein synthesis, and activation of oncogenic pathways including Wnt, HIF-1, PI3K-Akt, MAPK, Ras, p53, and mTOR signaling, which collectively drive tumor growth and survival.

The tumor microenvironment undergoes extensive remodeling. While normal colon tissue maintains a balanced cellular composition with robust immune surveillance, the tumor context is characterized by a proportional reduction of overall immune populations and significant shifts within immune cell subsets. Specifically, CD4+ T cells in tumors show an increase in immunosuppressive (Treg), pro-inflammatory (Th22, Th17) populations, and a decrease in T follicular helper cells (Tfh). Their GSEA profile indicates T cell anergy or exhaustion, with depletion of TCR signaling and antigen processing pathways, alongside enrichment in metabolic adaptations like glycolysis and HIF-1 signaling. Similarly, macrophage populations shift from M2A dominance in normal tissue to an M1-dominant state in tumors, accompanied by increased pro-tumorigenic M2B and M2D subsets, suggesting a complex and often contradictory immune response within the TME.

Cell-cell interaction analysis reveals a dramatic loss of productive epithelial-immune crosstalk in tumor tissue, where the top interactions are strikingly confined to T cell-T cell communication or dominated by stromal interactions. Notably, the critical IFN-gamma Type II Receptor (IFNGR) signaling, active in normal immune cells, is significantly diminished in tumor contexts, pointing to a key mechanism of immune evasion. Furthermore, the presence of VSIR (VISTA)-HLA-E interactions on tumor-infiltrating T cells underscores chronic T cell suppression. Conversely, tumor-associated fibroblasts (CAFs) exhibit extensive activation, marked by strong collagen-integrin interactions indicative of desmoplasia. Their gene expression profile highlights enrichment in ECM remodeling, growth factor signaling (TGF-beta, Wnt), and proliferative pathways (Cell cycle, DNA replication), affirming their crucial role in fostering a pro-tumorigenic niche. Surfaceome analysis further identifies key markers such as CEACAM5/6, EPCAM, CD47, and ERBB2 on tumor epithelial cells, and PDGFRB, CD248, and MMP14 on CAFs, providing molecular signatures of these transformed cell states.

Collectively, these findings paint a picture of colorectal cancer as a disease driven by genomic instability and aberrant epithelial cell proliferation, fostered by a profoundly altered and immunosuppressive tumor microenvironment. The coordinated dysregulation across epithelial, immune, and stromal compartments highlights the multifaceted challenges in treating this complex disease.

Hypotheses:

  1. The observed aneuploidy in Intestinal Epithelial cells is a primary driver of oncogenic pathway activation (e.g., Wnt, PI3K-Akt) and uncontrolled proliferation in colorectal cancer.
  2. The shift towards increased Treg, Th22, and M2B/M2D macrophage populations in the tumor microenvironment is a key mechanism of immune evasion, actively suppressing anti-tumor immune responses and promoting tumor progression.
  3. The breakdown of IFN-gamma signaling and the upregulation of immune checkpoints like CTLA4 and VSIR (VISTA) in tumor-infiltrating T cells contribute to T cell exhaustion and limit the efficacy of anti-tumor immunity.
  4. Cancer-associated fibroblasts (CAFs) actively remodel the extracellular matrix and secrete pro-tumorigenic factors, creating a supportive niche that facilitates tumor growth, invasion, and immune evasion through pathways like TGF-beta and Wnt signaling.
  5. The metabolic reprogramming (e.g., increased glycolysis, altered oxidative phosphorylation) observed across tumor epithelial cells, T cells, and fibroblasts is an adaptive response to the hypoxic tumor microenvironment, supporting rapid proliferation and survival.

Potential therapeutic targets:

  1. ERBB2 (HER2): ERBB2 amplification is a well-established oncogenic driver in various cancers and was identified as a recurrent amplification in tumor-associated cells (potentially misclassified Intestinal Epithelial cells) and as an upregulated surface marker on tumor Intestinal Epithelial cells. Evidence: Section 4 highlights frequent ERBB2 amplification (17q12:17q21.31) in tumor-associated cell groups. Section 15 identifies ERBB2 as a highly expressed surface marker on tumor Intestinal Epithelial cells, the tumor origin celltype. Validation: Confirm ERBB2 amplification in Intestinal Epithelial tumor cells via FISH. Evaluate the efficacy of HER2-targeted therapies (e.g., trastuzumab, pertuzumab) in patient-derived organoids or xenograft models with high ERBB2 expression/amplification.
  2. CD47: CD47 acts as a 'don't eat me' signal, allowing cancer cells to evade macrophage-mediated phagocytosis. Its upregulation on tumor Intestinal Epithelial cells indicates a key immune evasion mechanism. Evidence: Section 15 identifies CD47 as a highly expressed surface marker on tumor Intestinal Epithelial cells, suggesting it actively contributes to immune evasion in the tumor microenvironment. Validation: Test anti-CD47 antibodies in vitro (e.g., co-culture with macrophages) and in vivo (e.g., syngeneic mouse models or PDX models) to assess enhanced phagocytosis of tumor cells and reduced tumor growth.
  3. CTLA4 and VSIR (VISTA): These are immune checkpoint molecules upregulated on tumor-infiltrating T cells, contributing to T cell exhaustion and immunosuppression, a key feature of the colorectal cancer microenvironment. Evidence: Section 18 identifies CTLA4 as a prominent upregulated surface marker on tumor-associated CD4+ T cells. Section 14 highlights VSIR-HLA-E interactions in tumor-infiltrating T cells, suggesting active immune suppression. Validation: Evaluate the efficacy of anti-CTLA4 or anti-VISTA antibodies (alone or in combination) in preclinical models. Monitor T cell activation, proliferation, and cytokine production in response to checkpoint blockade in patient-derived T cells.
  4. PDGFRB and Integrin pathways (e.g., ITGA1, ITGAV, ITGB1): PDGFRB is a key receptor for CAF activation and proliferation, while integrins mediate crucial CAF-ECM interactions, driving desmoplasia, tumor growth, and invasion in the tumor microenvironment. Evidence: Section 17 shows PDGFRB, ITGA1, and ITGAV are highly upregulated surface markers on tumor fibroblasts. Section 21's GSEA identifies significant enrichment in 'ECM-receptor interaction', 'Focal adhesion', and 'TGF-beta signaling pathway' in tumor fibroblasts. Section 14 highlights numerous collagen-integrin interactions in tumor fibroblasts. Validation: Investigate PDGFRB inhibitors or integrin-blocking antibodies in in vitro CAF-tumor co-culture models to assess effects on CAF activation, ECM remodeling, and tumor cell invasion. Test these agents in vivo in relevant animal models of colorectal cancer.
  5. Cell Cycle Regulators (e.g., CDK4/6, MYC, CCND1): These genes are consistently and significantly upregulated in tumor Intestinal Epithelial cells, indicating uncontrolled proliferation, a hallmark of cancer. Evidence: Section 19 shows a widespread and statistically significant increase in the expressing cell fraction of genes like MYC, CCND1, CDK4, and CDK6 in tumor Intestinal Epithelial cells. Validation: Assess the sensitivity of colorectal cancer cell lines and patient-derived organoids to CDK4/6 inhibitors or drugs targeting MYC pathways. Evaluate the anti-proliferative effects and tumor growth inhibition in xenograft models.

Follow-up validation ideas:

  1. Perform spatial transcriptomics or multiplexed immunohistochemistry (e.g., CODEX, IMC) to validate the spatial distribution and co-localization of aneuploid Intestinal Epithelial cells with high cell cycle gene expression, and to map the interaction patterns of immunosuppressive immune cells and activated fibroblasts within the tumor microenvironment.
  2. Utilize flow cytometry or mass cytometry (CyTOF) on dissociated tumor and normal samples to quantify the absolute numbers and confirm the phenotypic shifts of T cell subsets (Tregs, Th17, Th22) and macrophage subsets (M1, M2A, M2B, M2D) identified by surface markers (e.g., CTLA4, CD39 for T cells; CD36, FCGR3A for macrophages).
  3. Conduct in vitro co-culture experiments using patient-derived tumor organoids or cancer cell lines with immune cells (T cells, macrophages) and fibroblasts to functionally validate the identified cell-cell interaction pathways (e.g., IFN-gamma signaling, VSIR-HLA-E, collagen-integrin) and their impact on immune cell function or tumor growth.
  4. Employ CRISPR-Cas9 or shRNA perturbation assays in tumor epithelial cells or CAFs to knockdown/overexpress key oncogenic drivers (e.g., MYC, CCND1, PDGFRB, MMP14) or immune modulators (e.g., CD47, ERBB2) and assess their functional consequences on proliferation, invasion, and immune evasion in vitro and in vivo.
  5. Validate recurrent genomic amplifications (e.g., ERBB2) identified in the CNV analysis using FISH (Fluorescence In Situ Hybridization) on tumor tissue sections to confirm amplification status and cellular localization, especially in putative misclassified T cells vs. true epithelial cells.
  6. Analyze independent cohorts of colorectal cancer patients using bulk RNA-seq or proteomics to validate the differential expression of key gene sets (e.g., cell cycle genes, metabolic pathways) and surface markers in tumor tissue, correlating findings with clinical outcomes and response to therapy.

Limitations:

This report is based on single-cell RNA sequencing data, which provides a snapshot of cellular states and infers cell-cell interactions. While robust computational methods are used, inferences regarding cell-cell communication, lineage relationships, and functional consequences require experimental validation. The detection of CNVs relies on transcriptional data and may not capture all genomic alterations accurately, as highlighted by the potential cell type misclassification issue in one CNV analysis. Furthermore, the dataset represents a specific cohort and tissue type, and findings may not be universally generalizable to all colorectal cancer subtypes or patient populations. The interpretation of complex GSA/GSEA pathways, especially those with broader biological relevance (e.g., neurodegenerative disease pathways in tumor cells), should be approached with caution, recognizing that they may reflect generalized cellular stress rather than specific disease mechanisms.

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. save
  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 cell cells as tumor-origin, 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 cell cells as tumor-origin, show their ploidy populations as a bar plot, and save the result.
  12. Show and save cell-cell interaction patterns involving Intestinal Epithelial cell, Fibroblast, Macrophage, and T cell. 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 B cell, Myeloid cell, Stromal cell, T cell, Endothelial cell, and Mast cell 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 cell, 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, 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, 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+, 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 cell, 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 cell as a bar plot.
  21. Show Gene Set Enrichment Analysis results for Intestinal Epithelial cell, T cell CD4+, Macrophage, and Fibroblast as a dot plot, and save it. Set the color map to RdBu_r and n_pws_to_show to 120.
↑ Top