SCODiA Report by MLBI Lab

Single-Cell RNA-Seq Analysis Reveals the Heterogeneous Tumor Microenvironment and Molecular Landscape of Kidney Tumors

This single-cell RNA sequencing analysis of kidney tissue elucidates the complex tumor microenvironment. Aneuploid renal epithelial cells are confirmed as the malignant population, exhibiting genomic instability and reprogrammed oncogenic pathways. The tumor microenvironment features significant immune cell infiltration, characterized by T cell exhaustion and altered macrophage polarization, alongside extensive cell-cell interactions driving angiogenesis and immune evasion. Key molecular features and signaling pathways, including EGFR, MET, and HIF-1, emerge as potential diagnostic biomarkers and therapeutic targets for kidney tumors.

Contents

  1. Dataset overview
  2. UMAP Visualization of Single-Cell Data by Condition, Sample, Cell Type, and Ploidy
  3. Major Cell Type Score and Ploidy Visualization on UMAP
  4. Marker Expression Profile Across Kidney Cell Subsets
  5. Genomic Copy Number Variation (CNV) Analysis in Renal Epithelial Cells from Kidney Tumors
  6. UMAP Visualization of CNV Patterns Colored by Cell Type, Ploidy, Condition, and Sample
  7. Minor Cell Type Population Analysis in Kidney Tumor
  8. T Cell Subset Population Analysis in Kidney Tumor vs. Normal Tissue
  9. Analysis of T cell (CD8 Exhausted) Population Proportions in Kidney Tumor
  10. Macrophage Subset Composition in Normal and Kidney Tumor Samples
  11. Macrophage M2C Subpopulation Expansion in Kidney Tumor Microenvironment
  12. Ploidy Profile of Renal Epithelial Cells in Normal and Tumor Kidney Samples
  13. Cell-Cell Interaction Patterns in Kidney Tumor Microenvironment
  14. Condition-Specific Cell-Cell Interactions in Kidney Tumor Microenvironment
  15. Condition-Specific Cell-Cell Interaction Patterns Involving Immune Checkpoint and Cell Cycle Pathways in Kidney Tumor
  16. Renal Epithelial Cell Condition-Specific Surfaceome Markers in Kidney Tumor
  17. Macrophage Condition-Specific Surfaceome Markers in Kidney Tumor
  18. Surfaceome Marker Expression in T cell CD4+ Cells Across Normal and Kidney Tumor Tissues
  19. Differential Expression of Cell Cycle Pathway Genes in Renal Epithelial Cells of Kidney Tumors
  20. Gene Ontology (GSA) Analysis of Renal Epithelial Cells in Kidney Tumor
  21. Gene Set Enrichment Analysis (GSEA) of Kidney Tumor Microenvironment
  22. Discussion
  23. Query List

0. Dataset overview

Dataset Summary

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

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents a series of Uniform Manifold Approximation and Projection (UMAP) plots, visualizing the global transcriptional landscape of single-cell RNA sequencing data from kidney tissue. Each UMAP is colored by a different metadata feature: 'condition' (normal vs. tumor), 'sample', 'celltype_major', 'celltype_minor', 'ploidy_dec' (diploid vs. aneuploid), and 'celltype_subset'. These plots serve to assess the overall structure of the dataset, evaluate the quality of clustering and annotation, and reveal the distribution of various biological and technical factors across the cellular populations. The primary goal is to understand how these features segregate or co-localize within the dimensionally reduced space, offering insights into cell identity, disease state, and genomic stability.

Visual Summary

Condition and Ploidy Distribution

Cell Type Clustering

Sample Distribution

Biological Interpretation

The UMAP visualizations provide compelling evidence for a kidney tumor microenvironment, consistent with the inferred biological context.

  1. Tumor Cell Identity and Genomic Instability: The strong co-localization of 'tumor' cells and 'Aneuploid' cells, particularly within the large central clusters, strongly suggests that these regions represent the malignant cell populations. Given that 'Renal Epithelial cell' is identified as the resolved tumor-origin cell type, the significant overlap of the 'Renal Epithelial cell' major cluster with 'tumor' and 'Aneuploid' regions indicates these are indeed the cancerous cells undergoing genomic alterations characteristic of tumorigenesis. Aneuploidy, or an abnormal number of chromosomes, is a hallmark of many cancers, driving tumor evolution and heterogeneity GeneCards: Aneuploidy.
  2. Tumor Microenvironment Heterogeneity: While 'Renal Epithelial cells' define the tumor core, other cell types, such as 'Myeloid cells', 'T cells', and 'Endothelial cells', are observed in both normal and tumor-enriched regions. This suggests their active participation in the tumor microenvironment (TME). The granular resolution of 'celltype_minor' and 'celltype_subset' allows for identification of specific immune cell states (e.g., different macrophage polarizations, T cell exhaustion states) that are likely shaped by or contributing to the tumor context. For instance, specific macrophage subsets (e.g., M1 vs M2) are known to play distinct roles in tumor progression or suppression PubMed: Tumor-associated macrophages.
  3. Annotation Quality and Consistency: The clear separation of cell types at major, minor, and subset levels demonstrates the high quality of cell clustering and annotation. The hierarchical nature of these annotations is visually consistent, with subsets refining broader categories. This robust annotation foundation is critical for downstream differential expression, pathway, and cell-cell interaction analyses.
  4. Sample Integration: The general mixing of cells from different samples across the UMAP indicates successful data integration, which minimizes patient-specific batch effects and allows for more confident biological comparisons across the cohort. However, the observation of some sample-enriched clusters might warrant further investigation to determine if these represent genuine inter-patient biological variability (e.g., different tumor stages or patient-specific immune responses) or minor residual technical variations.

Annotation Notes

2. Major Cell Type Score and Ploidy Visualization on UMAP

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis visualizes the confidence scores for major cell type assignments across the single-cell UMAP embedding. It also presents the final major cell type annotations and the estimated ploidy status (diploid vs. aneuploid) for each cell. This allows for validation of cell type assignments and provides critical insights into the genomic stability of different cell populations, especially in the context of a kidney tumor.

Visual Summary

The UMAP plots display the distribution of 26,502 cells in a 2-dimensional space, where proximity generally indicates higher similarity in gene expression profiles.

Major Cell Type Scores:

Ploidy Status (ploidy_dec):

Major Cell Type Annotation (celltype_major):

Biological Interpretation

  1. Tumor Cell Identification and Ploidy: The most striking observation is the strong co-localization of aneuploid cells (maroon in ploidy_dec) with the Renal Epithelial cell population (light green in celltype_major, and high scores in HiCAT_major_score: Renal Epithelial cell) on the right side of the UMAP. This finding is highly significant as Renal Epithelial cells have been identified as the tumor-origin cell type. The presence of a large aneuploid population within this cell type strongly supports the presence of a kidney tumor, a condition often characterized by genomic instability and chromosomal abnormalities in malignant epithelial cells PubMed search: kidney cancer aneuploidy. This effectively demarcates the tumor cell compartment from non-malignant cells in the tissue.
  2. Immune Cell Infiltration:
  1. Stromal and Endothelial Components:

Annotation Notes

The high congruence between the individual HiCAT_major_score UMAPs and the final celltype_major UMAP demonstrates robust and reliable cell type annotation. The distinct clustering and strong scoring for each major cell type provide confidence in the cell identity assignments within this single-cell RNA-seq dataset. The presence of a small 'unassigned' population is common in single-cell datasets and does not detract from the overall quality of the major cell type annotations.

3. Marker Expression Profile Across Kidney Cell Subsets

Report figure

[Analysis Visualization Results]...

Analysis Overview

This dot plot visualizes the expression of key marker genes across various celltype_subset annotations derived from kidney single-cell RNA-seq data. The plot aims to confirm the identity of these cell subsets by showing the mean expression level (color intensity) and the fraction of cells expressing each gene within a given cell group (dot size). Specific, highly expressed markers are highlighted with red boxes to indicate distinct molecular signatures for different cell populations. The analysis specifically focuses on surfaceome markers, which are valuable for cell type identification and potential downstream applications such as cell sorting.

Visual Summary

The dot plot effectively showcases distinct and overlapping gene expression patterns across 20 identified celltype_subset populations.

Biological Interpretation

The observed marker expression patterns largely confirm the distinct identities of the annotated celltype_subset populations, providing strong support for the granularity of the single-cell annotations.

Annotation Notes

The marker expression patterns observed in this dot plot provide strong evidence supporting the quality and biological consistency of the celltype_subset annotations within the dataset. The clear, specific gene expression profiles for most cell types suggest robust clustering and accurate annotation assignment. The distinction between closely related cell types, such as different macrophage polarization states and T cell subsets, further indicates a high resolution of cellular identity.

However, it is noted that some celltype_subset categories from the provided metadata (e.g., Lymphatic Endothelial cell, DC (Plasmacytoid), T cell (Th1), T cell (Treg), T cell (Th9), various B cell subsets, and some renal tubule segments) are not shown as rows in this specific plot. This could be due to filters applied during marker selection or plotting (e.g., minimum cell count per group, N_cells_per_group_min=30, or marker specificity thresholds), rather than an absence of these cell types in the dataset. Further investigation of these unrepresented subsets with less stringent criteria may be warranted to fully validate all annotations. Overall, the available data confidently supports the current celltype_subset assignments for the displayed populations.

4. Genomic Copy Number Variation (CNV) Analysis in Renal Epithelial Cells from Kidney Tumors

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates genomic copy number variations (CNVs) in Renal Epithelial cells, identified as the tumor-origin cell type, across different patient samples from kidney tumor tissue. The primary goal is to visualize CNV patterns and identify regions with significantly amplified copy numbers. The heatmap displays log2(CNR) values, indicating amplifications (red) and deletions (blue) across genomic spots, grouped by sample. A supplementary bar plot summarizes the most frequently amplified regions and associated genes within a specific sample (SI_21561), highlighting potential drivers of tumor development.

Visual Summary

CNV Heatmap

The heatmap illustrates distinct and heterogeneous CNV profiles across the four analyzed samples of Renal Epithelial cells.

Summary of Significantly Amplified Regions (Sample SI_21561)

The bar plot for sample SI_21561 details specific cytogenetic bands and associated genes that show significant amplifications. The numerical labels on the left represent the mean log2(CNR) for each region, while the bar length indicates the frequency or significance of the amplification event within this sample.

Biological Interpretation

The CNV profiles of Renal Epithelial cells, the identified tumor-origin cells, provide critical insights into the genomic alterations driving kidney tumor development.

Key Oncogenic Drivers and Tumor Suppressors:

Clinical or Translational Implications

The identification of recurrent CNVs, particularly in a tumor-origin cell type like Renal Epithelial cells, has several potential clinical and translational implications for kidney tumor management:

These findings highlight the importance of high-resolution genomic profiling, such as single-cell CNV analysis, in understanding the molecular landscape of kidney tumors and guiding the development of more effective, personalized cancer therapies.

5. UMAP Visualization of CNV Patterns Colored by Cell Type, Ploidy, Condition, and Sample

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents UMAP (Uniform Manifold Approximation and Projection) plots derived from single-cell RNA-seq data, specifically incorporating CNV (Copy Number Variation) estimates into the embedding calculation (cnv=True). The cells are colored by celltype_major, celltype_minor, ploidy_dec (ploidy decision), condition (normal vs. tumor), and sample to visualize their distribution and relationships in this CNV-informed transcriptional space. The goal is to assess cell type composition, identify aneuploid cell populations, discern tumor-specific cell clusters, and evaluate potential batch effects or biological heterogeneity across samples.

Visual Summary

The UMAP plots reveal distinct clusters of cells, suggesting robust cell type identification and a structured underlying biology.

Biological Interpretation

The UMAP visualizations provide strong evidence for the presence of a kidney tumor and illuminate its cellular composition and genomic alterations:

Annotation Notes

The visualizations confirm that the major and minor cell type annotations are generally well-defined, leading to distinct clusters in the UMAP space. The ploidy_dec and condition metadata align strongly with the distribution of the "Renal Epithelial cell" population, reinforcing the inferred tumor origin. While some sample-level heterogeneity is visible, it does not appear to obscure the primary biological distinctions between cell types, ploidy states, and conditions. This indicates a good quality dataset and effective integration of CNV information.

6. Minor Cell Type Population Analysis in Kidney Tumor

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents a bar plot illustrating the proportional distribution of minor cell types across individual samples, distinguishing between 'normal' and 'tumor' kidney tissues. The goal is to identify shifts in cellular composition associated with the tumor microenvironment and disease state, using single-cell RNA sequencing data. The dataset includes 26,502 cells from human kidney tissue, with a substantial portion of cells identified as aneuploid (34.5%) and Renal Epithelial cells resolved as the tumor-origin cell type.

Visual Summary

The visualization displays two panels, one for 'normal' samples (n=6) and one for 'tumor' samples (n=8), each showing stacked bar plots representing the percentage contribution of different minor cell types within each sample.

Immune Cell Infiltration in Tumor:

Biological Interpretation

The observed shifts in cell type proportions provide insights into the cellular remodeling occurring during kidney tumor development.

This fraction could contain crucial information regarding tumor heterogeneity or novel cellular players.

Clinical or Translational Implications

7. T Cell Subset Population Analysis in Kidney Tumor vs. Normal Tissue

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents a population bar plot illustrating the relative proportions of various T cell subsets, including NK cells, across individual samples from both normal kidney tissue and kidney tumor tissue. The goal is to identify differences in the immune cell composition within the T cell compartment between healthy and cancerous conditions, providing insights into the immune microenvironment of kidney tumors.

Visual Summary

The visualization displays stacked bar plots, with each bar representing a single sample and segments indicating the percentage of different T cell subsets. Samples are distinctly grouped by their 'normal' or 'tumor' condition.

Biological Interpretation

The observed shifts in T cell subset populations in kidney tumor samples strongly suggest an altered and potentially immunosuppressive tumor microenvironment (TME), consistent with the presence of a kidney tumor inferred from the data context.

Clinical or Translational Implications

The distinct immunological landscape observed in kidney tumor samples, characterized by an increase in exhausted CD8 T cells and Tregs, alongside a potential decrease in NK cells, has significant implications for cancer diagnosis, prognosis, and therapeutic strategies.

8. Analysis of T cell (CD8 Exhausted) Population Proportions in Kidney Tumor

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates the proportional representation of T cell (CD8 Exhausted) subsets within the total cell population, comparing normal kidney tissue samples against kidney tumor samples. The objective is to identify any statistically significant differences in the abundance of this specific T cell subset between the two conditions, which can provide insights into the immune landscape of kidney tumors.

Visual Summary

The box plot illustrates the celltype proportion of "T cell (CD8 Exhausted)" cells in 'normal' and 'tumor' conditions.

Biological Interpretation

The observed increase in T cell (CD8 Exhausted) proportions in kidney tumor samples, albeit trending towards statistical significance, is biologically relevant in the context of cancer immunology.

Although the p-value of 0.0701 does not strictly reach the typical 0.05 threshold for significance, it indicates a strong trend. With more samples or further validation, this difference could become statistically significant. The substantial magnitude of the difference (median ~0% vs. ~2% and outlier at 11%) strongly supports a biological distinction.

Clinical or Translational Implications

The enrichment of T cell (CD8 Exhausted) populations in kidney tumors has significant clinical implications:

9. Macrophage Subset Composition in Normal and Kidney Tumor Samples

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents a population bar plot illustrating the relative proportions of different macrophage subsets (Macrophage (M1), Macrophage (M2A), Macrophage (M2B), Macrophage (M2C), and Macrophage (M2D)) within individual normal and kidney tumor samples. The data provides insights into the immune cell landscape, specifically focusing on macrophage polarization states, in the context of kidney tumor. This is crucial given the inferred disease context of Kidney Tumor and the known role of macrophages in cancer immunology.

Visual Summary

The stacked bar plots display the relative proportions of five distinct macrophage subsets across several normal and tumor kidney samples.

M2 Subsets Contribution:

Biological Interpretation

Macrophages are critical components of the tumor microenvironment, and their polarization into M1 (pro-inflammatory, anti-tumor) and M2 (pro-tumor, immune-suppressive, tissue repair) phenotypes significantly influences cancer progression. The observed data reveals an interesting pattern in kidney tumors:

Implications of M1 Dominance:

Clinical or Translational Implications

The distinct macrophage polarization profile observed in kidney tumors, particularly the apparent dominance of M1 macrophages, carries several clinical and translational implications:

References:

[1] Macrophage Polarization in Cancer: A Therapeutic Target. PubMed search: https://pubmed.ncbi.nlm.nih.gov/?term=macrophage+polarization+cancer+M1+M2

10. Macrophage M2C Subpopulation Expansion in Kidney Tumor Microenvironment

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents a box plot visualizing the proportion of Macrophage (M2C) cells within the overall cellular landscape, comparing kidney tumor samples against normal kidney tissue samples. The goal is to identify statistically significant differences in the abundance of this specific macrophage subset between the two conditions.

Visual Summary

The box plot illustrates a notable difference in the proportion of Macrophage (M2C) cells between normal and tumor conditions.

Biological Interpretation

The observed trend of an increased proportion of Macrophage (M2C) cells in kidney tumor tissue compared to normal tissue has significant biological implications in the context of tumor immunology and progression.

Clinical or Translational Implications

The enrichment of Macrophage (M2C) cells in kidney tumors holds several potential clinical and translational implications:

11. Ploidy Profile of Renal Epithelial Cells in Normal and Tumor Kidney Samples

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates the ploidy status (diploid vs. aneuploid) of Renal Epithelial cells, which have been identified as the tumor-origin cell type within the dataset. The ploidy populations are visualized for individual samples, stratified by their associated condition: "normal" or "tumor." This helps to assess the genomic stability of these crucial cells in both healthy and diseased kidney tissue contexts, providing insights into the molecular characteristics of the presumed kidney tumor.

Visual Summary

The bar plot effectively illustrates the distribution of diploid and aneuploid Renal Epithelial cells across various samples.

Biological Interpretation

The observed ploidy patterns strongly corroborate the inference of kidney tumor and highlight a fundamental aspect of cancer biology.

Clinical or Translational Implications

The distinct ploidy profiles observed between normal and tumor Renal Epithelial cells carry significant clinical and translational relevance in the context of kidney cancer.

12. Cell-Cell Interaction Patterns in Kidney Tumor Microenvironment

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates cell-cell interaction (CCI) patterns within the kidney tissue, comparing normal and tumor conditions. The focus is on key cell types identified in the tumor microenvironment (TME), including renal epithelial cells (distinguished by ploidy status as 'Diploid Renal Epi' and 'Aneuploid Renal Epi' given the tumor context), fibroblasts, macrophages, CD4+ T cells, and CD8+ T cells. The CellPhoneDB method was utilized to identify significant ligand-receptor interactions, aggregated by condition (normal vs. tumor). The expand_ploidy_from_tumor_origin parameter was enabled, allowing for a specific examination of aneuploid renal epithelial cells, which are identified as the tumor-origin cells. The results highlight critical changes in intercellular communication in the tumor setting, particularly involving aneuploid renal epithelial cells and immune cells like macrophages.

Visual Summary

The visualization presents two dot plots, one for the 'normal' and one for the 'tumor' condition, showcasing significant cell-cell interactions. The size of each dot corresponds to the -log10(p-value) of the interaction, indicating statistical significance (larger dots represent greater significance), while the color intensity reflects the log2(mean) expression of the interacting ligand-receptor pair (brighter colors, such as yellow/green, denote stronger interaction strength).

In the normal condition plot, a very limited number of interactions are observed. Specifically, only auto-crine interactions within 'Diploid Renal Epi' cells are evident, involving pairs such as APLP2-PIGR, PPIA-BSG, and ProstaglandinE2_byPTGES2-PTGER3. These interactions show high significance and moderate-to-strong expression. The scarcity of other interactions suggests a relatively quiescent or less complex communication landscape in healthy kidney tissue among the selected cell types.

In stark contrast, the tumor condition plot reveals a highly complex and active network of cell-cell interactions. A substantially larger number of significant ligand-receptor pairs are identified, particularly involving 'Aneuploid Renal Epi' cells (the inferred tumor-origin cells) and 'Mac' (Macrophages). Interactions are also prominent within 'Mac | Mac' (macrophage self-interactions) and 'Aneuploid Renal Epi | Aneuploid Renal Epi' (tumor cell auto-crine signaling).

Key observations from the tumor plot include:

Biological Interpretation

The stark contrast between the normal and tumor conditions underscores a profound rewiring of intercellular communication in the kidney tumor microenvironment. The minimal interactions in normal tissue emphasize the highly regulated and homeostatic nature of healthy renal epithelium.

The extensive interactions observed in the tumor condition, especially those involving 'Aneuploid Renal Epi' cells, provide crucial insights into tumor progression and immune evasion.

  1. Tumor-Immune Cell Crosstalk (Aneuploid Renal Epi ↔ Macrophage):
  1. Tumor Cell Autocrine and Homotypic Signaling (Aneuploid Renal Epi ↔ Aneuploid Renal Epi):
  1. Macrophage Self-Regulation (Mac ↔ Mac):

The overall increase in interaction complexity and strength in the tumor context, particularly involving the aneuploid renal epithelial cells, strongly supports their role as the driving force behind the altered cellular communication network. The identified interactions underscore mechanisms of immune evasion, angiogenesis, and sustained proliferation, which are hallmarks of cancer.

Clinical or Translational Implications

The identified cell-cell interaction patterns offer several clinically relevant insights for kidney tumor management:

  1. Therapeutic Target Prioritization:
  1. Biomarker Discovery: The specific patterns of cell-cell communication, or the expression levels of key interacting partners (e.g., CD47, SIRPA, HBEGF, EREG, NAMPT) in tumor tissue, could serve as prognostic or predictive biomarkers for patient stratification and response to targeted therapies.
  2. Experimental Validation: The numerous high-significance interactions identified warrant further experimental validation. *In vitro* co-culture models of aneuploid renal epithelial cells with macrophages, or *in vivo* xenograft models, could be used to confirm functional relevance. Immunohistochemistry or spatial transcriptomics on patient kidney tumor samples could validate the physical proximity and expression patterns of these ligand-receptor pairs. These validations are critical for translating these findings into novel therapeutic strategies for kidney cancer.

13. Condition-Specific Cell-Cell Interactions in Kidney Tumor Microenvironment

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis utilizes single-cell RNA sequencing data to identify cell-cell interactions (CCI) that significantly differ between normal and kidney tumor conditions. The dot plot visualizes the strength (scaled log expression mean, color intensity) and statistical significance (-log10(p-value), dot size) of various ligand-receptor interactions across different samples and conditions. The focus is on interactions involving renal epithelial cells (identified as the tumor-origin cell type), major immune cells (Myeloid, T cell, B cell, Mast cell, Plasma cell), and stromal/endothelial cells, highlighting the dynamic changes in intercellular communication within the tumor microenvironment. The analysis identifies interactions that are more prominent or stronger in the tumor context compared to normal tissue.

Visual Summary

The dot plot effectively illustrates condition-specific patterns of cell-cell interactions.

Biological Interpretation

The observed shifts in cell-cell interactions provide critical insights into kidney tumor biology:

  1. Enhanced Angiogenesis and Vascular Remodeling in Tumor:
  1. Tumor-Immune Cell Crosstalk and Immune Evasion:
  1. Pro-tumorigenic Signaling:

Clinical or Translational Implications

The profound changes in cell-cell interactions in kidney tumors offer multiple avenues for therapeutic intervention:

  1. Targeting Angiogenesis and Vascular Support:

Therapeutic Targets:

Therapeutic Agents:

  1. Modulating Tumor-Associated Macrophages (TAMs):

Therapeutic Targets:

Therapeutic Agents:

  1. Disrupting Extracellular Matrix Remodeling:

Overall, the increased and diversified cell-cell interactions in kidney tumors highlight the complexity of the tumor microenvironment. Targeting these specific interactions, particularly those driving angiogenesis, immune modulation, and tumor cell survival/migration, represents a promising strategy for developing novel and effective therapeutic approaches for kidney cancer.

14. Condition-Specific Cell-Cell Interaction Patterns Involving Immune Checkpoint and Cell Cycle Pathways in Kidney Tumor

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates condition-specific cell-cell interactions (CCI) within kidney tissue, focusing on genes related to immune checkpoints and cell cycle pathways. The study examines interactions involving renal epithelial cells (identified as the tumor-origin cell type), major immune cells (Myeloid, T cell, Mast cell, Plasma cell, B cell), and stromal cells. The interactions are compared between 'normal' and 'tumor' conditions across various patient samples. The dot plot visualizes the interaction strength (color intensity, scaled log strength) and statistical significance (dot size, -log10(p-value)) for identified ligand-receptor pairs and cell type combinations. Importantly, interactions are distinguished based on the ploidy status of renal epithelial cells (Diploid vs. Aneuploid), which is crucial for identifying tumor-specific processes.

Visual Summary

The dot plot displays a clear divergence in cell-cell interaction patterns between normal and tumor conditions.

Normal Condition (Left Panel):

Tumor Condition (Right Panel):

Biological Interpretation

The analysis reveals several critical cell-cell interaction axes that are distinctly altered in kidney tumor compared to normal tissue, predominantly involving aneuploid renal epithelial cells.

1. EGFR Signaling Axis in Kidney Tumor

Interactions Identified (Tumor Condition, predominantly involving aneuploid Renal Epithelial cells):

AREG_EGFR -- Mac | Renal.Epi (Aneup)

EREG_EGFR -- Mac | Renal.Epi (Aneup)

HBEGF_EGFR -- Mac | Renal.Epi (Aneup)

TGFA_EGFR -- Renal.Epi (Aneup) | Renal.Epi (Aneup)

These interactions show significantly increased strength and prevalence in tumor samples compared to their comparatively weaker or absent counterparts in normal tissue involving Renal.Epi (Dip).

Signaling Pathways and Mechanisms: The Epidermal Growth Factor Receptor (EGFR) is a receptor tyrosine kinase that plays a pivotal role in cell growth, proliferation, survival, and differentiation. Upon ligand binding (e.g., AREG, EREG, HBEGF, TGFA, EGF, BTC), EGFR dimerizes and activates intracellular signaling cascades, notably the RAS/MAPK, PI3K/AKT/mTOR, and STAT pathways. In the tumor microenvironment, tumor-associated macrophages (Mac) are a significant source of EGFR ligands (AREG, EREG, HBEGF), driving proliferation and survival of aneuploid renal epithelial cells. Furthermore, TGFA_EGFR interactions within aneuploid renal epithelial cells suggest an autocrine signaling loop, fostering self-sufficiency in growth signals, a hallmark of cancer. This widespread activation of EGFR signaling is a major contributor to uncontrolled cell growth, angiogenesis, and resistance to apoptosis in kidney tumors.

2. TGF-beta/Integrin Signaling in the Tumor Microenvironment

Interactions Identified (Tumor Condition, predominantly involving aneuploid Renal Epithelial cells):

TGFB3_integrin_aVb8_complex -- SMC | Renal.Epi (Aneup)

TGFB3_TGFbeta_receptor1 -- SMC | Renal.Epi (Aneup)

These interactions are highly active and robust in tumor samples, showing a marked increase compared to normal conditions, where TGFB1 related interactions involving Renal.Epi (Dip) and Mac were present but less pervasive.

Signaling Pathways and Mechanisms: Transforming Growth Factor beta (TGF-β) is a pleiotropic cytokine with context-dependent roles in cancer. While it can suppress early tumor development, it often promotes tumor progression, metastasis, angiogenesis, epithelial-mesenchymal transition (EMT), and immune evasion in advanced cancers. Integrins, particularly the αVβ8 complex, are known receptors for TGF-β, facilitating its activation and signaling. The strong interaction between TGFB3 from Smooth Muscle Cells (SMC, representing stromal components) and aneuploid Renal Epithelial cells via integrin_aVb8_complex and TGFbeta_receptor1 highlights a critical crosstalk between tumor cells and the surrounding stroma. This interaction likely contributes to the remodeling of the extracellular matrix, creation of a desmoplastic reaction, and further driving tumor progression and immune suppression within the kidney tumor microenvironment.

3. Immune Cell Co-stimulation and Communication in Tumor Infiltration

Interactions Identified (Tumor Condition):

CD86_CD28 -- Mac | T CD8+

IFNG_Type_II_IFNR -- T CD8+ | Mac

These interactions are visible with moderate to strong significance in several tumor samples, whereas they are largely absent or negligible in normal tissue samples. The LCK_CD8_receptor interaction involving T CD8+ cells is also observed in tumor samples, likely indicating intracellular signaling associated with T cell activation.

Signaling Pathways and Mechanisms:

Clinical or Translational Implications

The identified cell-cell interactions offer significant insights into kidney tumor biology and potential therapeutic strategies.

1. EGFR Signaling Axis as a Therapeutic Target

Potential as a Therapeutic Target: The widespread activation of EGFR signaling, driven by both macrophages and autocrine loops in aneuploid renal epithelial cells, strongly suggests that EGFR remains a viable and potent therapeutic target in kidney tumors. Inhibiting this pathway can directly curb tumor cell proliferation, survival, and angiogenesis.

Relevant Therapeutic Agents:

These agents are widely used in other cancers, and their efficacy in specific kidney tumor subtypes with activated EGFR signaling could be explored GeneCards: EGFR; PubMed search: "EGFR inhibitors renal cell carcinoma".

2. TGF-beta/Integrin Signaling for Modulating Tumor Microenvironment

Potential as a Therapeutic Target: The strong stromal-tumor cell interaction mediated by TGFB3 and integrin αVβ8 in tumor samples indicates that this axis is a critical driver of tumor progression and immune evasion. Targeting TGF-β signaling or specific integrins could attenuate tumor growth, reduce metastasis, and potentially reprogram the immunosuppressive tumor microenvironment.

Relevant Therapeutic Agents:

Modulating this pathway could also enhance the efficacy of other immunotherapies by reducing immune suppression UniProt: TGFB3; PubMed search: "TGF-beta integrin inhibitors kidney cancer".

3. Immune Cell Interactions for Immunomodulation

Potential as a Therapeutic Target: The presence of CD86_CD28 and IFNG_Type_II_IFNR interactions in tumor samples suggests an attempt by the immune system to respond, or at least active communication, but it might be insufficient or dysregulated. Modulating these interactions could enhance anti-tumor immunity.

Relevant Therapeutic Agents:

Understanding the balance of these co-stimulatory and co-inhibitory signals is crucial for designing effective immunotherapies for kidney cancer GeneCards: CD28; GeneCards: IFNG.

In conclusion, the analysis highlights crucial dysregulations in growth factor, stromal-tumor, and immune cell interactions within kidney tumors. Targeting these specific pathways, particularly EGFR and TGF-β signaling, and judiciously modulating immune cell crosstalk, hold promise for novel therapeutic interventions.

15. Renal Epithelial Cell Condition-Specific Surfaceome Markers in Kidney Tumor

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis identifies surfaceome markers differentially expressed in Renal Epithelial cells, comparing normal and tumor conditions using single-cell RNA-seq data. The goal is to highlight potential diagnostic and therapeutic targets specific to kidney tumor cells of epithelial origin. Markers were identified using a one-vs-rest comparison strategy with specific scoring and fold change cutoffs, focusing exclusively on surface-expressed proteins. The results are visualized using a dot plot, illustrating both the mean expression level and the fraction of cells expressing each marker across different samples grouped by condition.

Visual Summary

The dot plot displays a clear differential expression pattern of surfaceome markers in Renal Epithelial cells between normal and tumor conditions.

Biological Interpretation

The differential surfaceome expression in Renal Epithelial cells provides critical insights into the biological changes associated with kidney tumorigenesis.

Markers Downregulated in Tumor (Normal-Specific)

EPCAM (Epithelial Cell Adhesion Molecule)

DPEP1 (Dipeptidase 1)

Markers Upregulated in Tumor (Tumor-Specific)

EGFR (Epidermal Growth Factor Receptor)

MET (MET proto-oncogene, receptor tyrosine kinase)

CD47 (CD47 molecule)

VCAM1 (Vascular Cell Adhesion Molecule 1)

CDH6 (Cadherin 6)

CXCR4 (C-X-C Motif Chemokine Receptor 4)

CA12 (Carbonic Anhydrase 12)

  1. Aberrant Antigen Presentation: Tumor cells might aberrantly express MHC Class II to present tumor antigens, but with potential defects in processing or co-stimulatory signals that lead to T cell anergy or exhaustion.
  2. Immune Context: Reflecting an inflammatory tumor microenvironment where epithelial cells are induced to express MHC Class II, perhaps as an attempt to present antigens, which might be subverted by the tumor.
  3. Immune Evasion: In some contexts, MHC Class II expression on tumor cells can lead to T cell exclusion or a tolerogenic response.

Clinical or Translational Implications

The identified condition-specific surfaceome markers in Renal Epithelial cells present several clinical and translational opportunities, particularly for kidney tumor diagnosis, prognostication, and therapeutic targeting.

Established Therapeutic Targets

EGFR:

MET:

CD47:

Potential Therapeutic Targets / Biomarkers

VCAM1:

CDH6:

CXCR4:

CA12:

MHC Class II Associated Proteins (e.g., HLA-DRB1, HLA-DRA):

Overall, this analysis identifies a rich set of surfaceome markers in Renal Epithelial cells that are differentially regulated in kidney tumor. These findings provide a valuable foundation for further investigation into their roles in tumor biology and their potential as novel diagnostic biomarkers and therapeutic targets for kidney cancer.

16. Macrophage Condition-Specific Surfaceome Markers in Kidney Tumor

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis identifies surfaceome markers differentially expressed in macrophages isolated from kidney tumor tissue compared to macrophages from normal kidney tissue. The goal is to highlight potential biomarkers that characterize tumor-associated macrophages (TAMs) and evaluate their suitability as therapeutic targets. The dot plot visualizes the mean expression level (color intensity) and the fraction of cells expressing each marker (dot size) across individual samples grouped by 'normal' and 'tumor' conditions.

Visual Summary

The dot plot clearly segregates macrophage surfaceome markers into two main groups based on their expression patterns across conditions:

The stringent filtering criteria (nz_pct_score and nz_pct_FC) applied for both normal_vs_others and tumor_vs_others comparisons ensure that the displayed markers are robustly specific to their respective conditions within the macrophage population.

Biological Interpretation

The observed surfaceome changes in macrophages reflect their adaptation and functional specialization within the kidney tumor microenvironment (TME). Tumor-associated macrophages (TAMs) are crucial orchestrators of tumor progression, influencing angiogenesis, immune suppression, metastasis, and therapeutic resistance. The identified tumor-specific surface markers provide insights into these roles and represent promising targets for therapeutic intervention.

Key Tumor Macrophage Markers and Their Implications:

  1. CSF1R (Colony Stimulating Factor 1 Receptor)
  1. AXL (AXL Receptor Tyrosine Kinase)
  1. GPNMB (Glycoprotein Non-Metastatic Melanoma Protein B / Osteoactivin)
  1. LAIR1 (Leukocyte-associated immunoglobulin-like receptor 1)
  1. FOLR2 (Folate Receptor Beta)

Clinical or Translational Implications

The identification of specific surfaceome markers on macrophages within kidney tumors provides critical insights for the development of novel diagnostics and therapeutics.

Further validation of these markers in a larger cohort and functional studies are warranted to confirm their precise roles in kidney tumor progression and fully realize their therapeutic potential.

17. Surfaceome Marker Expression in T cell CD4+ Cells Across Normal and Kidney Tumor Tissues

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis identifies and visualizes surfaceome markers for T cell CD4+ cells, comparing their expression patterns between normal kidney tissue and kidney tumor tissue. The goal is to highlight condition-specific markers that may offer insights into the functional states of these crucial immune cells within the tumor microenvironment (TME) and potentially serve as therapeutic targets. Markers were selected based on their differential expression using a one-versus-rest comparison strategy for each condition, focusing on surfaceome-localized proteins.

Visual Summary

The dot plot displays the mean expression and fraction of cells expressing specific surfaceome markers in T cell CD4+ populations from normal and tumor conditions.

Biological Interpretation

The differential expression of surfaceome markers between normal and tumor-associated T cell CD4+ populations suggests distinct functional adaptations and roles of these cells within the kidney tumor microenvironment.

Markers Enriched in Normal T cell CD4+ (SI_22369)

PTPRC (CD45):

AREG (Amphiregulin):

CD52:

Markers Enriched in Tumor T cell CD4+ (SI_22368)

CD3D (CD3 delta chain):

KLRB1 (CD161):

CD74 (Invariant Chain):

CD4 (CD4 molecule):

HLA-DRA (MHC Class II DR Alpha chain):

LY6E (Lymphocyte Antigen 6 Family Member E):

CD63:

Other Markers Enriched in Tumor T cell CD4+ (SI_22368)

Clinical or Translational Implications

The identification of condition-specific surfaceome markers on T cell CD4+ cells in kidney tumors provides valuable insights for potential diagnostic, prognostic, and therapeutic strategies:

Therapeutic Targets:

18. Differential Expression of Cell Cycle Pathway Genes in Renal Epithelial Cells of Kidney Tumors

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates the "expressing cell fraction (sample)" of a selected set of cell cycle pathway genes in 'Renal Epithelial cells', which have been identified as the resolved tumor-origin cell type in this kidney dataset. The comparison is made between 'tumor' and 'normal' tissue conditions. The goal is to identify genes within this critical pathway that show statistically significant differences in the proportion of cells expressing them across conditions, providing insights into cell cycle dysregulation in kidney tumors.

Visual Summary

The box plots display the distribution of the "expressing cell fraction (sample)" for nine cell cycle-related genes across 'tumor' and 'normal' conditions in Renal Epithelial cells. Each point represents an individual sample. Statistical p-values comparing the two conditions are annotated above each pair of box plots.

Key observations from the plots include:

Biological Interpretation

The consistent observation of a significantly lower fraction of renal epithelial cells expressing key cell cycle pathway genes in the tumor condition compared to normal tissue is highly relevant to kidney tumor biology. This indicates that within the tumor microenvironment, a smaller proportion of the tumor-origin renal epithelial cells maintain expression of these critical regulators.

Collectively, these findings suggest that in kidney tumors, the population of renal epithelial cells shows a pervasive deregulation of fundamental processes controlling cell growth, division, and genomic integrity. The reduction in the *fraction of cells expressing* these genes, rather than just their average expression level, points towards a selection or shift in subpopulations of tumor cells that have lost or downregulated these crucial regulatory functions, thereby contributing to the malignant phenotype.

Clinical or Translational Implications

The observed downregulation of critical cell cycle regulators and DNA repair genes in a significant fraction of renal epithelial cells within kidney tumors holds several potential clinical implications:

19. Gene Ontology (GSA) Analysis of Renal Epithelial Cells in Kidney Tumor

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates the biological processes and pathways enriched in Renal Epithelial cells, identified as the tumor-origin cell type in this kidney dataset. Gene Ontology (GSA) analysis was performed on genes significantly up-regulated when comparing Renal Epithelial cells in the 'normal' condition against other cell types, and similarly, when comparing 'tumor' Renal Epithelial cells against other cell types. This provides insights into the characteristic functions of healthy renal epithelial cells versus the transformed pathways active in tumor-derived cells.

Visual Summary

The provided bar plots illustrate the top enriched Gene Ontology (GO) terms based on up-regulated differentially expressed genes (DEGs) for Renal Epithelial cells in two comparisons: 'normal_vs_others' and 'tumor_vs_others'. The x-axis represents the statistical significance as -log(p-val) and -log(q-val), indicating the strength of enrichment for each term.

Normal Renal Epithelial Cells (normal_vs_others):

Tumor Renal Epithelial Cells (tumor_vs_others):

Biological Interpretation

Normal Renal Epithelial Cells: Core Physiological Functions

The GSA results for normal Renal Epithelial cells highlight their specialized and highly active metabolic roles essential for kidney function. The strong enrichment of oxidative phosphorylation, glycolysis, pyruvate metabolism, and various amino acid and fatty acid degradation pathways underscores their high energy demand for processes like filtration, reabsorption, and secretion. For example, proximal tubule cells, a subtype of renal epithelial cells, are highly energetic and rely heavily on oxidative phosphorylation to drive active transport processes PubMed search: kidney proximal tubule oxidative phosphorylation. The presence of terms like "Collecting duct secretion" and "Proximal tubule bicarbonate reabsorption" directly confirms the specific physiological functions of these healthy kidney cells.

Tumor Renal Epithelial Cells: Malignant Transformation and Immune Dysregulation

The GO terms enriched in tumor Renal Epithelial cells reflect a profound shift in cellular biology, consistent with oncogenic transformation:

Comparative Insights

The comparison clearly delineates the functional divergence between healthy and malignant Renal Epithelial cells. Normal cells maintain a robust metabolic machinery geared towards kidney-specific physiological tasks. In contrast, tumor cells undergo a comprehensive reprogramming, activating pathways associated with uncontrolled growth, altered protein processing, and a complex interplay with inflammatory and immune-related processes, including those typically triggered by viral challenges. This suggests a highly adaptable and aggressive phenotype in the tumor cells.

Clinical or Translational Implications

  1. Biomarker Discovery and Diagnostics: The distinct sets of highly enriched pathways in tumor Renal Epithelial cells could yield novel diagnostic or prognostic biomarkers for kidney cancer. For instance, specific genes within the HIF-1, PI3K-Akt, or NF-kappa B pathways, or components of the ER stress response, could be explored as indicators of tumor presence or aggressiveness.
  2. Therapeutic Targets: The identification of activated oncogenic pathways (HIF-1, PI3K-Akt, NF-kappa B, MAPK) provides clear targets for therapeutic intervention. Inhibitors against these pathways are already available or in development for various cancers, suggesting potential for repurposing or developing new targeted therapies for kidney cancer PubMed search: targeted therapy renal cell carcinoma pathways. Targeting protein processing or ER stress could also be a viable strategy.
  3. Understanding Tumor Microenvironment and Immune Evasion: The strong enrichment of viral infection-associated and inflammatory pathways warrants further investigation. This could signify mechanisms by which tumor cells evade immune surveillance, promote a pro-tumorigenic inflammatory environment, or mimic infectious states. A deeper understanding could lead to novel immunotherapeutic strategies.
  4. Drug Resistance Mechanisms: Over-activation of some pathways might contribute to drug resistance. Monitoring the activity of these pathways could help predict treatment response and guide personalized therapy.

20. Gene Set Enrichment Analysis (GSEA) of Kidney Tumor Microenvironment

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents Gene Set Enrichment Analysis (GSEA) results across various cell types identified in single-cell RNA-seq data from kidney tissue, comparing gene expression in cells from the tumor condition against other cells (predominantly normal cells), and vice versa. The dot plot visualizes the Normalized Enrichment Score (NES) and statistical significance (-log(P)) for 120 selected pathways. The primary goal is to identify biological processes and pathways that are significantly altered in different cell populations within the kidney tumor microenvironment.

Visual Summary

The dot plot displays a matrix where rows represent specific gene sets or pathways, and columns represent comparisons for different cell types: Endothelial cell, Macrophage, Renal Epithelial cell, Smooth muscle cell, and T cell CD4+. For each cell type, two comparisons are shown: 'normal_vs_others' (genes enriched in normal cells compared to the tumor condition) and 'tumor_vs_others' (genes enriched in tumor cells compared to the normal condition).

Biological Interpretation

The GSEA results highlight significant shifts in cellular processes and pathways across different cell types within the kidney tumor microenvironment, consistent with cancer progression and host-tumor interactions.

Renal Epithelial Cells (Tumor Origin)

As the resolved tumor-origin cell type, Renal Epithelial cells exhibit hallmarks of malignant transformation:

Upregulation in Tumor Cells (Red, positive NES):

Downregulation in Tumor Cells (Blue, negative NES):

Immune Cells (Macrophages and CD4+ T Cells)

The immune microenvironment shows active but potentially dysfunctional responses:

Macrophages (Tumor vs. others):

T Cell CD4+ (Tumor vs. others):

Stromal and Vascular Cells (Endothelial Cells and Smooth Muscle Cells)

These cells play crucial roles in shaping the tumor microenvironment:

Endothelial Cells (Tumor vs. others):

Smooth Muscle Cells (Tumor vs. others):

Clinical or Translational Implications

The GSEA results provide critical insights into the underlying biology of kidney tumors, with several potential clinical implications:

21. Discussion

This comprehensive single-cell analysis of kidney tissue provides profound insights into the cellular composition and molecular reprogramming within the kidney tumor microenvironment (TME). The foundational observation is the unequivocal identification of renal epithelial cells as the tumor-origin population, strongly corroborated by their pervasive aneuploidy—a hallmark of malignancy. CNV analysis further dissects this genomic instability, revealing widespread alterations, including recurrent *KRAS* amplification, which is a known oncogenic driver, alongside unexpected amplifications of the tumor suppressors *CDKN2A* and *SMAD4* that warrant further investigation into co-amplified passenger genes or context-dependent roles.

The malignant renal epithelial cells exhibit substantial molecular reprogramming, as evidenced by GSA and GSEA. These analyses highlight a shift from the robust metabolic functions characteristic of normal kidney epithelial cells (e.g., oxidative phosphorylation, amino acid metabolism) to pathways indicative of uncontrolled growth, including hyperactivated oncogenic signaling cascades (PI3K-Akt, HIF-1, mTOR, NF-κB, MAPK), altered protein processing in the ER, and a Warburg-like metabolic reprogramming towards glycolysis. Furthermore, a significant fraction of tumor renal epithelial cells show reduced expression of key cell cycle regulators and DNA repair genes like *RB1*, *TGFB1*, *PRKDC*, *HDAC1*, and *CDC27*, underscoring a pervasive deregulation of genomic integrity and cell cycle control.

The TME itself is dynamically restructured. Population analyses reveal a reduction in normal renal epithelial cell types and a notable increase in macrophages and fibroblasts, indicative of an active stromal and immune infiltrate. The immune landscape is particularly complex: while a strong trend of increased T cell (CD8 Exhausted) populations and regulatory T cells (Tregs) suggests an immunosuppressive environment, Macrophage M1 (pro-inflammatory) phenotypes surprisingly dominate in tumor samples, coexisting with a significant increase in pro-tumorigenic M2C macrophages. This nuanced macrophage polarization implies either a dynamic balance of M1- and M2-like functions or a specific immune response profile in kidney tumors that deviates from other cancer types, demanding further functional characterization. GSEA for immune cells reinforces this complexity, showing T cell activation alongside PD-1 checkpoint pathway engagement and potential T cell receptor signaling downregulation, suggestive of an exhausted state.

Cell-cell interaction analyses reveal a dramatic increase in complexity and strength of communication networks within the tumor. Critical interactions include the 'don't eat me' signal via SIRPA-CD47 between aneuploid renal epithelial cells and macrophages, which is a primary mechanism of immune evasion. Robust EGFR-ligand interactions (e.g., AREG, EREG, HBEGF, TGFA) drive tumor cell proliferation and survival in an autocrine manner and through crosstalk with macrophages. Angiogenesis and stromal remodeling are highly active, characterized by elevated Notch, PDGF, VEGF, and Collagen-Integrin interactions involving endothelial and smooth muscle cells. Furthermore, TGF-beta/integrin signaling between stromal cells and tumor cells promotes tumor progression and immune suppression. The presence of CD86-CD28 and IFNG-Type_II_IFNR interactions in tumor-infiltrating immune cells suggests ongoing, albeit potentially dysregulated, immune responses.

Finally, the identification of condition-specific surfaceome markers across tumor-origin cells, macrophages, and CD4+ T cells offers a rich reservoir for clinical translation. Upregulated markers on tumor renal epithelial cells include established targets like EGFR, MET, and CD47, alongside VCAM1, CDH6, CXCR4, and unusual MHC Class II proteins. Tumor-associated macrophages uniquely express CSF1R, AXL, GPNMB, LAIR1, and FOLR2, reflecting their pro-tumorigenic roles. These markers, together with the insights from pathway and interaction analyses, converge to highlight critical nodes that can be leveraged for developing novel diagnostic tools, prognostic biomarkers, and targeted therapeutic strategies for kidney cancer.

Hypotheses:

  1. Genomic instability (aneuploidy, specific CNVs like KRAS amplification) in renal epithelial cells is a primary driver of their malignant transformation, leading to dysregulation of cell cycle control and DNA repair mechanisms and contributing to the aggressive phenotype of kidney tumors.
  2. The kidney tumor microenvironment actively promotes T cell exhaustion and modulates macrophage polarization towards a pro-tumorigenic state through specific cell-cell interactions (e.g., CD47-SIRPA, EGFR-ligand axes), thereby enabling immune evasion and sustained tumor progression.
  3. The diverse and condition-specific surfaceome profiles of tumor-associated renal epithelial cells, macrophages, and CD4+ T cells reflect their functional adaptations to the tumor microenvironment and serve as a rich source of novel diagnostic biomarkers and precision therapeutic targets.
  4. Metabolic reprogramming (e.g., increased glycolysis and protein synthesis) and the activation of specific oncogenic signaling cascades (e.g., PI3K-Akt, HIF-1, mTOR) are indispensable for the survival, proliferation, and adaptive capabilities of tumor-origin cells and supporting stromal components within the kidney tumor microenvironment.

Potential therapeutic targets:

  1. EGFR (Epidermal Growth Factor Receptor): Upregulated and activated in aneuploid renal epithelial cells, driven by both autocrine signaling and ligands from macrophages. This pathway is a major driver of cell proliferation, survival, and angiogenesis in kidney tumors. Evidence: Upregulated as a tumor-specific surfaceome marker in renal epithelial cells (Section 15). Strong cell-cell interactions involving EGFR and its ligands (AREG, EREG, HBEGF, TGFA) identified in tumor (Sections 12, 14). GSEA indicates enrichment of PI3K-Akt and MAPK signaling pathways (Section 20). Validation: Evaluate EGFR Tyrosine Kinase Inhibitors (e.g., erlotinib, gefitinib) or monoclonal antibodies (e.g., cetuximab) in kidney cancer cell lines and *in vivo* models to assess tumor growth inhibition and modulation of downstream signaling.
  2. MET (MET proto-oncogene, receptor tyrosine kinase): Overexpressed in renal epithelial tumor cells and associated with aggressive disease, invasion, and metastasis, promoting cell scattering and angiogenesis. Evidence: Identified as a significantly upregulated tumor-specific surfaceome marker in renal epithelial cells (Section 15). Validation: Test MET inhibitors (e.g., cabozantinib, crizotinib) in MET-overexpressing kidney cancer models. Cabozantinib's existing approval for advanced RCC provides strong clinical relevance for this target.
  3. CD47 (CD47 molecule): Highly upregulated on tumor renal epithelial cells, where it interacts with SIRPα on macrophages, providing a 'don't eat me' signal to evade phagocytosis and contribute to immune escape. Evidence: Identified as a significantly upregulated tumor-specific surfaceome marker in renal epithelial cells (Section 15). Strong and prevalent SIRPA-CD47 cell-cell interactions observed in tumor samples (Section 12). Validation: Evaluate anti-CD47 monoclonal antibodies (e.g., magrolimab) in preclinical kidney tumor models. Perform *in vitro* co-culture assays to demonstrate enhanced macrophage-mediated phagocytosis of tumor cells upon CD47 blockade.
  4. CSF1R (Colony Stimulating Factor 1 Receptor): Crucial for the recruitment, survival, and pro-tumorigenic M2-like polarization of tumor-associated macrophages (TAMs), which promote tumor growth, angiogenesis, and immune evasion within the TME. Evidence: Identified as a highly upregulated tumor-specific surfaceome marker in macrophages (Section 16). GSEA shows active recruitment, phagocytosis, and metabolic reprogramming of macrophages (Section 20). Validation: Test CSF1R small molecule inhibitors (e.g., Pexidartinib) or monoclonal antibodies (e.g., Cabiralizumab) in kidney tumor models to assess TAM depletion or reprogramming, and subsequent effects on tumor growth and immune responses.
  5. AXL (AXL Receptor Tyrosine Kinase): Highly expressed on TAMs; AXL signaling contributes to their immunosuppressive functions, promotes M2 polarization, efferocytosis, angiogenesis, metastasis, and confers resistance to therapies. Evidence: Identified as a highly upregulated tumor-specific surfaceome marker in macrophages (Section 16). Validation: Test AXL inhibitors (e.g., bemcentinib, gilteritinib) in kidney cancer models. Assess their impact on tumor growth, TAM phenotype, and the reversal of therapy resistance, particularly in combination with other agents.
  6. GPNMB (Glycoprotein Non-Metastatic Melanoma Protein B): Significantly expressed in TAMs and various tumor cells; promotes tumor cell proliferation, survival, invasion, and angiogenesis. Its restricted expression in normal tissues makes it an attractive target for antibody-drug conjugates (ADCs). Evidence: Identified as a highly upregulated tumor-specific surfaceome marker in macrophages (Section 16). Validation: Evaluate GPNMB-targeted antibody-drug conjugates (e.g., Glembatumumab vedotin) in preclinical kidney tumor models to demonstrate selective delivery and cytotoxicity against GPNMB-expressing tumor-promoting cells.
  7. FOLR2 (Folate Receptor Beta): Selectively expressed on activated TAMs but generally absent on quiescent macrophages and most normal cells. Mediates folate internalization, supporting tumor growth by providing necessary metabolites. Ideal for precision drug delivery. Evidence: Identified as a highly upregulated tumor-specific surfaceome marker in macrophages (Section 16). Validation: Develop and test folate-conjugated cytotoxic agents or anti-FOLR2 antibodies for targeted delivery to TAMs in kidney tumor models, aiming to enhance efficacy and minimize off-target effects.
  8. TGF-beta Signaling Pathway (via TGFB3 and integrin αVβ8): Strong stromal-tumor cell interaction, mediated by TGFB3 from smooth muscle cells and its receptors (integrin_aVb8_complex, TGFbeta_receptor1) on aneuploid renal epithelial cells. This pathway is a critical driver of tumor progression, metastasis, epithelial-mesenchymal transition, and immune evasion. Evidence: Highly active and robust TGFB3_integrin_aVb8_complex and TGFB3_TGFbeta_receptor1 interactions between Smooth Muscle Cells and aneuploid Renal Epithelial cells in tumor (Section 14). GSEA indicates ECM remodeling (Section 20). Validation: Test TGF-β pathway inhibitors (e.g., galunisertib) or integrin αVβ8 blocking antibodies in *in vitro* co-culture models of stromal and tumor cells and *in vivo* kidney tumor models to assess effects on tumor growth, metastasis, and TME immunosuppression.
  9. CXCR4 (C-X-C Motif Chemokine Receptor 4): Overexpressed in tumor renal epithelial cells. Its interaction with CXCL12 promotes tumor cell migration, invasion, and metastasis by guiding cells towards CXCL12-rich metastatic sites, and recruits macrophages to the tumor. Evidence: Identified as an upregulated tumor-specific surfaceome marker in renal epithelial cells (Section 15). CXCL14-CXCR4 interaction identified in tumor (Section 13). Validation: Evaluate CXCR4 antagonists (e.g., plerixafor, motixafortide) in preclinical models of kidney cancer metastasis, assessing reduction in metastatic burden and potential enhancement of immune cell infiltration.
  10. NOTCH Signaling Pathway: Upregulated Notch ligand-receptor interactions (e.g., JAG1/2/3/4-NOTCH1/2/3/4, DLL1-NOTCH4) are prominent in tumor endothelial and smooth muscle cells. Notch signaling is a critical regulator of endothelial cell differentiation, sprouting angiogenesis, and vascular remodeling, all vital for pathological neovascularization in tumors. Evidence: Strong cell-cell interactions involving Notch ligands and receptors detected in tumor endothelial and smooth muscle cells (Section 13). GSEA indicates active angiogenesis with enrichment of related pathways (Section 20). Validation: Test gamma-secretase inhibitors (GSIs), which block Notch activation, in angiogenesis assays and *in vivo* kidney tumor models to assess effects on tumor vasculature formation and overall tumor growth inhibition.
  11. PDGFRB (Platelet-Derived Growth Factor Receptor Beta): The PDGFB-PDGFRB interactions between endothelial and smooth muscle cells are highly active in tumors. This pathway is essential for recruiting pericytes to stabilize newly formed blood vessels, thereby supporting sustained tumor growth and resistance to anti-angiogenic therapies. Evidence: Prominent cell-cell interactions involving PDGFB and PDGFRB in tumor endothelial and smooth muscle cells (Section 13). GSEA shows active angiogenesis with enrichment of related pathways (Section 20). Validation: Evaluate PDGFR inhibitors (e.g., sunitinib, sorafenib, already cornerstones in RCC treatment) or specific PDGFRB blocking agents in kidney tumor models to further characterize their anti-angiogenic and anti-tumor effects, especially in contexts of resistance.

Follow-up validation ideas:

  1. Functional Validation of CNVs: Use FISH or array CGH on sorted aneuploid renal epithelial cells from kidney tumor samples to orthogonally confirm specific CNV regions identified (e.g., KRAS amplification). Perform CRISPR/shRNA-mediated perturbation of candidate genes (e.g., KRAS) in kidney cancer cell lines to assess their functional impact on proliferation, survival, and drug sensitivity *in vitro*.
  2. T Cell Exhaustion & Function: Conduct multiplexed flow cytometry or mass cytometry on tumor-infiltrating lymphocytes (TILs) from a larger patient cohort to precisely quantify the proportion and co-expression of exhaustion markers (e.g., PD-1, TIM-3, LAG-3) on CD8+ T cells, comparing tumor and adjacent normal tissue. Functional assays (e.g., cytokine production upon antigen restimulation) can assess the actual immune effector capacity of these subsets.
  3. Macrophage Polarization & Role: Employ spatial transcriptomics or multiplex immunostaining (e.g., IMC, CyTOF imaging) on kidney tumor tissue sections to validate the spatial distribution, density, and M1/M2 polarization state (e.g., CD68, CD80, CD163, FOLR2) of macrophages. *In vitro* co-culture experiments with kidney cancer cell lines and polarized macrophage subsets could investigate their reciprocal effects on tumor cell proliferation, migration, and immune suppression.
  4. Cell-Cell Interaction Validation: Utilize spatial transcriptomics or high-resolution imaging techniques (e.g., immunohistochemistry with proximity ligation assay) on patient kidney tumor samples to confirm the co-localization and protein expression of key ligand-receptor pairs (e.g., CD47-SIRPA, EGFR-ligands, TGFB3-Integrin αVβ8). *In vitro* blocking experiments using antibodies against these interactions in co-culture models could establish their functional relevance in tumor growth and immune evasion.
  5. Therapeutic Target Validation: Conduct *in vivo* studies using patient-derived xenografts (PDX) or syngeneic kidney tumor models, testing candidate therapeutic agents (e.g., anti-CD47 antibodies, MET inhibitors, AXL inhibitors, CSF1R inhibitors, CXCR4 antagonists). Evaluate their efficacy in reducing tumor growth, modulating the TME (e.g., TAM depletion/reprogramming), and improving survival outcomes.
  6. Cell Cycle Dysregulation: Perform multiplex RNA-FISH or immunohistochemistry on paired normal and tumor kidney tissues to quantify the fraction of renal epithelial cells expressing key cell cycle regulators (e.g., RB1, PRKDC, HDAC1) at the protein and RNA level, validating the single-cell RNA-seq findings within their spatial context.

Limitations:

The inference of tumor-origin cells and the disease context relies on metadata and pre-computed genomic features, rather than confirmed clinical diagnosis. Single-cell RNA-seq captures gene expression but does not directly quantify protein levels or post-translational modifications, which are crucial for full functional activity. While CNV estimates are available, they are based on RNA-seq and warrant orthogonal genomic validation (e.g., WGS, FISH). Computational predictions of cell-cell interactions require rigorous experimental validation through spatial analyses (e.g., spatial transcriptomics, multiplexed imaging) or functional co-culture assays. The substantial 'unassigned' cell population in the dataset might obscure novel or highly aberrant cell states critical to tumor biology. Furthermore, some observed trends, such as the increase in exhausted T cells (p=0.0701) and M2C macrophages (p=0.0577), suggest biological significance but would benefit from larger cohort sizes for definitive statistical conclusions. Inter-sample heterogeneity, as observed in CNV profiles and cell population compositions, necessitates cautious generalization of findings, and the precise functional consequences of some observed changes (e.g., MHC Class II expression on tumor cells, the nuanced M1/M2 macrophage balance) are complex and context-dependent, requiring further investigation beyond correlative data.

22. Query List

  1. Show UMAPs for condition, sample, major cell type, minor cell type, ploidy_dec, and celltype_subset in two columns, and save the results.
  2. Show major cell type scores on UMAP and save the result.
  3. Show a marker expression dot plot for celltype_subset and save the result. SET target_cell = None and var_group_rotation = 45. Leave all other arguments at their default values.
  4. Select Renal Epithelial cells (tumor-origin 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 UMAPs of CNV patterns colored by major cell type, minor cell type, ploidy results, condition, and sample in two columns, and save the results.
  6. Show a population bar plot of minor cell types and save the result.
  7. Show a population bar plot of T cell subsets and save the result.
  8. Show box plots of T cell subset populations with statistically significant differences between conditions, if any, and save the result. Choose ncols based on the total number of panels.
  9. Show a population bar plot of macrophage subsets and save the result.
  10. Show box plots of macrophage subset populations with statistically significant differences between conditions, if any, and save the result. Choose ncols based on the total number of panels.
  11. Select Renal Epithelial cells (tumor-origin cells), show their ploidy populations as a bar plot, and save the result.
  12. Show cell-cell interaction patterns involving Renal Epithelial cells (tumor-origin cells), fibroblasts, macrophages, T cells, and other relevant cell types, and save the results. Select at most 80 cell-cell interactions per group.
  13. Find cell-cell interactions involving tumor-origin cells, 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. For each interaction identified under each condition, using the provided dataset context, describe the associated signaling pathways and mechanisms, its potential as a therapeutic target, and relevant therapeutic agents.
  14. Find cell-cell interactions containing immune checkpoint and cell cycle pathway genes and involving tumor-origin cells, 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. For each interaction identified under each condition, using the provided dataset context, describe the associated signaling pathways and mechanisms, its potential as a therapeutic target, and relevant therapeutic agents.
  15. Show marker expressions for Renal Epithelial cells (tumor-origin cells) in a dot plot, and save the result. Use only surfaceome markers, with at most 60 markers per condition. For the markers identified under each condition, distinguish between established therapeutic targets and less well-characterized markers with potential as therapeutic targets. For each marker, using the provided dataset context, describe the associated signaling pathways and mechanisms, relevant therapeutic agents, and supporting evidence.
  16. Show marker expressions for Macrophage in a dot plot, and save the result. Use only surfaceome markers, with at most 40 markers per condition. For each marker identified under each condition, using the provided dataset context, describe the associated signaling pathways, the role of the corresponding cell type in the tumor microenvironment, the marker’s potential as a therapeutic target, relevant therapeutic agents, and other pertinent information.
  17. Show marker expressions for T cell CD4+ in a dot plot, and save the result. Use only surfaceome markers, with at most 40 markers per condition. For each marker identified under each condition, using the provided dataset context, describe the associated signaling pathways, the role of the corresponding cell type in the tumor microenvironment, the marker’s potential as a therapeutic target, relevant therapeutic agents, and other pertinent information.
  18. Among cell cycle pathway genes, select those with statistically significant expression differences between conditions in major disease-related 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.
  19. Show Gene Ontology (GSA) analysis results for Renal Epithelial cells (tumor-origin cells) as a bar plot, and save the result.
  20. Show Gene Set Enrichment Analysis results as a dot plot, and save it. Set the color map to RdBu_r, n_pws_to_show to 120, and target_cells to None.
↑ Top