SCODiA Report by MLBI Lab

Single-Cell Resolution of Cellular Heterogeneity, Genomic Instability, and Immunosuppressive Microenvironment Remodeling in Human Lung Cancer Progression

This comprehensive single-cell analysis of human lung tissue reveals the intricate landscape of lung cancer progression. Malignant lung epithelial cells, characterized by pervasive aneuploidy and activated oncogenic pathways, drive tumor growth and extensively remodel the tumor microenvironment. Key findings highlight dynamic shifts in immune cell populations, particularly the early and sustained enrichment of immunosuppressive regulatory T cells and pro-tumorigenic M2-like macrophages. Altered cell-cell communication networks, involving prominent EGFR, TGF-beta, and immune checkpoint signaling, contribute to immune evasion, angiogenesis, and stromal support, offering critical insights for therapeutic intervention.

Contents

  1. Dataset overview
  2. UMAP Visualization of Lung Cancer Single-Cell Data
  3. Major Cell Type Score and Ploidy Analysis on UMAP
  4. Marker Expression Dot Plot for Cell Subtypes in Lung Tissue
  5. Copy Number Variation (CNV) Analysis in Tumor-Origin Lung Epithelial Cells
  6. UMAP Visualization of CNV Patterns by Cell Type, Ploidy, Condition, and Sample
  7. Cell Type Population Analysis by Disease Condition in Lung Tissue
  8. T Cell Subset Population Dynamics in Lung Cancer Progression
  9. Differential Abundance of T Cell Subsets and ILC2s in Lung Cancer Progression
  10. Macrophage Subpopulation Shifts Across Lung Cancer Stages
  11. Macrophage Subset Population Shifts Across Lung Cancer Progression
  12. Ploidy Analysis of Tumor-Origin Lung Epithelial Cells
  13. Cell-Cell Interaction Patterns in Lung Cancer Microenvironment
  14. Condition-Specific Cell-Cell Interaction Patterns in Lung Cancer
  15. Condition-Specific Cell-Cell Interaction Patterns in Lung Cancer
  16. Surfaceome Marker Expression in Lung Epithelial Cells of Lung Cancer
  17. Macrophage Condition-Specific Surfaceome Markers in Lung Cancer
  18. Fibroblast Condition-Specific Surfaceome Markers in Lung Cancer
  19. Condition-Specific Surfaceome Markers in CD4+ T Cells in Lung Cancer
  20. Differential Expression of Cell Cycle Genes in Lung Epithelial Cells Across Lung Cancer Stages
  21. Gene Ontology Analysis of Lung Epithelial Cells Across Conditions
  22. Gene Set Enrichment Analysis of Lung Cancer Microenvironment
  23. Discussion
  24. Query List

0. Dataset overview

Dataset Summary

Inferred Disease/Biological Context

1. UMAP Visualization of Lung Cancer Single-Cell Data

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents Uniform Manifold Approximation and Projection (UMAP) plots, an unsupervised dimensionality reduction technique, to visualize the global cellular heterogeneity of single-cell RNA-sequencing data from human lung tissue. The UMAPs are colored by various metadata attributes including condition (Normal, Tumor(early), Tumor(adv)), sample identifier, celltype_major, celltype_minor, ploidy_dec (Aneuploid, Diploid, Unclear), and celltype_subset. The purpose of these visualizations is to assess the quality of the cell embeddings, evaluate the consistency of cell type annotations, investigate the distribution of different conditions and samples, and identify populations based on their ploidy status within the complex tissue microenvironment of Lung Cancer.

Visual Summary

  1. Condition UMAP

The UMAP colored by condition displays a noticeable separation between Normal cells (dark red) and cells from Tumor(adv) (yellow) and Tumor(early) (purple) samples. Normal cells predominantly occupy distinct clusters, suggesting unique transcriptional profiles. Tumor(adv) and Tumor(early) cells show areas of intermingling, indicating some shared cellular states, but also form distinct clusters, hinting at specific gene expression programs associated with disease progression stages.

  1. Sample UMAP

The sample UMAP reveals a high degree of sample-specific clustering, with many individual samples forming discrete groupings. This indicates either strong patient-specific biological variation or potential technical batch effects. However, several clusters exhibit a mix of multiple samples, suggesting the presence of shared cell types or conditions across different patient samples.

  1. Celltype Major UMAP

Clear segregation of major cell types is observed. T cell (cyan/blue) populations form a large, distinct cluster, while Myeloid cell (light green), Lung Epithelial cell (orange), Stromal cell (teal), B cell (dark red), and Mast cell (yellow) each occupy well-defined regions. Lung Epithelial cell populations appear in several distinct clusters, including regions that overlap with cells identified as Tumor(adv) and Aneuploid in other plots.

  1. Celltype Minor UMAP

At a finer resolution, the celltype_minor UMAP further refines the major cell types. T cell CD8+ (blue) and T cell CD4+ (dark blue) are distinct but remain interconnected. Other populations like Macrophage (yellow), Dendritic cell (orange), Alveolar Epithelial cell (dark red), and Airway Epithelial cell (red) also form well-demarcated clusters, demonstrating the robust resolution of these populations.

  1. Ploidy_dec UMAP

The ploidy_dec UMAP is particularly informative, showing that Aneuploid cells (dark red) are largely confined to specific, distinct clusters. These clusters visually correspond to regions enriched for Tumor(adv) and Lung Epithelial cell populations. Diploid cells (yellow) are broadly distributed across most of the UMAP landscape, as expected for healthy and non-malignant cells. The Unclear cells (purple) are very few and sparsely distributed, indicating high confidence in the ploidy assignments for the vast majority of cells.

  1. Celltype Subset UMAP

The celltype_subset UMAP provides the most granular view of cellular heterogeneity. It resolves numerous fine-grained populations such as Alveolar type 2 (AT2, dark red), Alveolar type 1 (AT1, red), diverse Macrophage subtypes (e.g., Mac_M1, Mac_M2a, Mac_M2b), and various T cell subsets (e.g., T_Cyto, T_Naive, Tfh, Treg). These highly resolved populations underscore the extensive cellular complexity within the lung tissue and tumor microenvironment.

Biological Interpretation

Annotation Notes

2. Major Cell Type Score and Ploidy Analysis on UMAP

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis visualizes the distribution of major cell types and ploidy status across the single-cell RNA-seq dataset using Uniform Manifold Approximation and Projection (UMAP). The primary objective is to assess the quality of cell type annotations by examining cell type-specific scores on the UMAP embedding and to determine the spatial relationship between annotated cell types and chromosomal aneuploidy.

Visual Summary

The provided UMAP plots display the major cell type scores (HiCAT_major_score) for T cells, B cells, Myeloid cells, Mast cells, Endothelial cells, Stromal cells, and Lung Epithelial cells. These are shown alongside plots colored by ploidy status (ploidy_dec) and the assigned major cell type (celltype_major).

  1. Cell Type Score Distribution: Each HiCAT_major_score plot clearly delineates regions of high score (indicated by yellow/green color) corresponding to distinct clusters on the UMAP. For example, high T cell scores are concentrated in the large left-central and upper-left clusters, while high Myeloid cell scores are seen in the right-central region. Lung Epithelial cell scores are particularly high in the upper-right quadrant, indicating a distinct population.
  2. Ploidy Status Distribution: The ploidy_dec plot highlights a prominent population of aneuploid cells (red) predominantly located in the upper-right region of the UMAP. The vast majority of other cells are diploid (light yellow), with a negligible number classified as "Unclear."
  3. Major Cell Type Annotation: The celltype_major plot shows the manually assigned major cell types, with distinct colors for each type. These clusters largely align with the regions of high scores observed in the individual HiCAT_major_score plots. For instance, the T cell cluster in celltype_major corresponds directly to the high T cell score region. Crucially, the Lung Epithelial cell cluster (light orange) in the upper-right region directly overlaps with the aneuploid cell population shown in the ploidy_dec plot.

Biological Interpretation

The UMAP visualizations provide compelling insights into the cellular composition and disease-relevant characteristics of the lung tissue sample:

Annotation Notes

3. Marker Expression Dot Plot for Cell Subtypes in Lung Tissue

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents a marker expression dot plot for various celltype_subset populations identified in the single-cell RNA-seq dataset from human lung tissue. The plot visualizes both the mean expression level (color intensity, ranging from white to dark red) and the fraction of cells expressing a given gene (dot size) for a curated set of marker genes. The purpose of this plot is to assess the specificity and robustness of the celltype_subset annotations by examining the expression patterns of known and discovered marker genes, ensuring that each cluster exhibits distinct and appropriate transcriptional signatures. The analysis applies a var_group_rotation of 45 degrees for readability, and target_cell was set to None to display markers for all identified celltype_subset populations.

Visual Summary

The dot plot effectively visualizes marker gene expression across 37 distinct celltype_subset populations. Each row represents a celltype_subset, and each column represents a marker gene. Genes are grouped along the x-axis, with red rectangular boxes highlighting sets of genes highly specific to certain cell type groups.

Biological Interpretation

The marker expression patterns generally provide strong support for the assigned celltype_subset annotations, consistent with known cell biology in lung tissue and the context of Lung Cancer.

Lung Epithelial Cells

Immune Cells

Stromal and Endothelial Cells

Annotation Notes and Potential Refinements

The dot plot generally confirms the high quality of the celltype_subset annotations, with most cell types exhibiting distinct and biologically plausible marker gene expression profiles. The use of multiple specific markers for each subtype enhances confidence in the assignments.

However, a notable observation that warrants further investigation is the expression of *KRT17* in B cell (MZ). Keratin 17 is typically an intermediate filament protein expressed in epithelial cells, particularly in basal cells or certain proliferative epithelial states. Its expression in a B cell population could indicate a potential misannotation, a unique biological state (e.g., interaction with transformed epithelial cells leading to aberrant gene expression, or a rare cell state), or a technical artifact. This specific annotation could benefit from cross-referencing with additional B cell-specific markers not shown, or by re-evaluating the clustering of these specific B cells. Otherwise, the identified markers robustly support the existing cell type subset annotations.

4. Copy Number Variation (CNV) Analysis in Tumor-Origin Lung Epithelial Cells

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates copy number variations (CNVs) specifically within Lung Epithelial cells identified as tumor-origin, derived from single-cell RNA-seq data. The primary goal is to visualize the genomic landscape of CNVs across different samples and to identify recurrent regions of amplification. The results are presented as a heatmap illustrating log2(Copy Number Ratio) across genomic spots, grouped by sample, alongside a summary table and bar plot detailing significantly amplified cytogenetic bands and their frequencies. This analysis is crucial for understanding the genomic instability characteristic of lung cancer and for identifying potential oncogenic drivers.

Visual Summary

Figure 1: CNV Heatmap of Lung Epithelial Cells by Sample

The heatmap displays the log2(Copy Number Ratio) for Lung Epithelial cells, organized by genomic spot (chromosomes 1-22) on the x-axis and individual samples on the y-axis.

Figure 2: Summary of Significantly Amplified Cytogenetic Bands

This figure provides a detailed summary of recurrent amplification events.

1q21.3:1q23.2

20p11.21:20q11.23

5q31.1:5q31.3

8q24.3:8q24.3 (GSDMD)

These plots confirm the presence of significant and recurrent genomic amplifications in the tumor-origin Lung Epithelial cells.

Biological Interpretation

The CNV analysis provides strong evidence of genomic instability in the identified tumor-origin Lung Epithelial cells, consistent with the inferred Lung Cancer context. The clear distinction between diploid samples and tumor-associated samples validates the cell type and ploidy annotations.

Other Recurrent Amplifications:

Clinical or Translational Implications

The identification of recurrent and significant amplifications of oncogenes like *EGFR* and *ERBB2* has direct clinical implications for lung cancer management:

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

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents UMAP (Uniform Manifold Approximation and Projection) plots generated using CNV (Copy Number Variation) estimates. These visualizations aim to represent the cellular landscape based on genomic alterations and are colored by various cellular and sample-level metadata, including major cell types (celltype_major), minor cell types (celltype_minor), ploidy status (ploidy_dec), experimental condition (condition), and individual samples (sample). The CNV-based embedding allows for the exploration of how genomic instability structures the cell populations and relates to cell identity, disease state, and sample origin in the context of Lung Cancer.

Visual Summary

The UMAP plots reveal distinct clustering patterns driven by CNV profiles, effectively segregating different cellular populations.

Biological Interpretation

The UMAP visualizations provide compelling evidence regarding the cellular and genomic landscape of the lung tissue samples, particularly in the context of Lung Cancer.

  1. Tumor Cell Identification and Origin: The CNV-based embedding effectively separates malignant cells from non-malignant cells. The ploidy_dec plot clearly shows a distinct cluster of aneuploid cells, which, when cross-referenced with celltype_major and celltype_minor, is almost entirely comprised of Lung Epithelial cells. This observation strongly aligns with the pre-defined Resolved tumor-origin cell type: Lung Epithelial cell and the inferred disease context of Lung Cancer, where malignancy typically arises from epithelial cells.
  2. Genomic Instability as a Primary Driver of Cell Clustering: The profound separation of cells based on ploidy_dec (aneuploid vs. diploid) confirms that CNV status is a major factor shaping the UMAP embedding. This highlights the extensive genomic instability present in tumor cells and underscores its utility in distinguishing malignant from non-malignant populations in single-cell data.
  3. Tumor Microenvironment (TME) Composition: The large, diffuse diploid cluster represents the non-malignant cells of the TME, including diverse immune populations (T cells, Myeloid cells) and stromal components (Fibroblasts). These cells are found in both "Normal" and "Tumor" conditions, as expected for cells that infiltrate or reside within the tumor site but do not harbor the extensive aneuploidy characteristic of the tumor clone. The presence of non-malignant cells in tumor conditions reflects the complex ecosystem of the TME.
  4. Condition-Specific Cell Distributions: The condition UMAP shows a clear enrichment of "Normal" cells in the diploid clusters, while "Tumor(early)" and "Tumor(adv)" cells are enriched in both the aneuploid (tumor epithelial) and diploid (TME) compartments. This further supports the distinct genomic profiles of tumor cells and the heterogeneous cellular composition of the tumor environment at different disease stages. The mixing of early and advanced tumor stages within the aneuploid cluster suggests that the primary CNV features captured by this embedding might be broadly shared across tumor progression, or that more subtle CNV changes between stages are not fully resolved at this level.

Annotation Notes

6. Cell Type Population Analysis by Disease Condition in Lung Tissue

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents a population bar plot illustrating the relative proportions of different celltype_minor categories across individual samples, stratified by disease condition: Normal, Tumor (advanced), and Tumor (early). Each bar represents a single sample, with stacked segments showing the percentage contribution of each identified minor cell type. This visualization provides insights into the cellular composition shifts associated with lung cancer progression.

Visual Summary

The stacked bar plots reveal distinct differences in cell type proportions between normal lung tissue and tumor samples, and also between early and advanced tumor stages.

Biological Interpretation

The observed shifts in cell populations provide strong biological insights into lung cancer development and its impact on the tissue microenvironment.

Clinical or Translational Implications

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

7. T Cell Subset Population Dynamics in Lung Cancer Progression

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents a stacked bar plot visualizing the relative proportions of T cell subsets and related innate lymphoid cells (ILCs and NK cells) across different samples categorized into Normal, Tumor (early), and Tumor (advanced) conditions. The data is derived from single-cell RNA-seq of human lung tissue, focusing on the 'T cell' major population and its further refinement into 'subset' level categories. The primary goal is to identify shifts in the immune cell composition within the T cell compartment during lung cancer progression.

Visual Summary

The bar plot displays the percentage contribution of various T cell and innate lymphoid cell subsets for each individual sample, grouped by condition: Normal (9 samples), Tumor(adv) (3 samples), and Tumor(early) (8 samples).

Biological Interpretation

The observed changes in T cell subset populations provide critical insights into the immune microenvironment of lung cancer.

  1. Enrichment of Regulatory T cells (Tregs) in Tumor Conditions: A prominent finding is the consistent increase in T cell (Treg) populations in both early and advanced tumor samples compared to normal lung tissue. Tregs are potent immunosuppressive cells that play a crucial role in maintaining immune tolerance and preventing autoimmunity. Their enrichment within the tumor microenvironment is a well-established mechanism by which tumors evade anti-tumor immune responses, promoting tumor growth and progression PubMed search: regulatory T cells lung cancer immune evasion. This suggests that immunosuppression is an early event in lung cancer development and is sustained throughout progression.
  2. Implication for Anti-tumor Immunity: The relative increase in Tregs, alongside the likely unchanged or relatively decreased proportions of potentially anti-tumor effector cells like T cell (Cytotoxic) cells, indicates a shift towards an immunosuppressive tumor microenvironment. This altered balance can hinder effective anti-tumor immune responses, allowing cancer cells to proliferate unchecked.
  3. "Unassigned" Population in Advanced Tumors: The emergence of a more substantial "unassigned" population in some advanced tumor samples warrants further investigation. This could represent highly perturbed or novel T cell states induced by the advanced tumor microenvironment, or potentially cells with ambiguous transcriptomic profiles that do not clearly fit into predefined subsets. Characterizing these cells could reveal new insights into immune dysregulation in advanced disease.
  4. Minor Role of ILCs and NK cells in T cell Compartment Shifts: While important immune components, the proportions of ILCs (ILC1, ILC2, ILC3, ILCreg, LTI) and NK cells remain consistently low and do not show significant proportional shifts within the T cell compartment across disease stages in this analysis. This suggests that changes in T cell subsets, particularly Tregs, are the primary drivers of immune compositional changes observed within this specific analysis of the broad T cell population.

Clinical or Translational Implications

The findings have several potential clinical and translational implications for lung cancer management:

  1. Therapeutic Targeting of Tregs: The consistent and early enrichment of T cell (Treg) cells highlights them as attractive targets for immunotherapeutic interventions. Strategies aimed at depleting Tregs or inhibiting their suppressive function could help restore anti-tumor immunity and improve the efficacy of existing treatments, such as immune checkpoint inhibitors GeneCards: FOXP3, UniProt: FOXP3.
  2. Biomarker Potential: The proportion of T cell (Treg) cells in the tumor microenvironment could serve as a prognostic biomarker for lung cancer progression or as a predictive marker for response to immunotherapies. Monitoring these cellular shifts could inform personalized treatment strategies.
  3. Understanding Advanced Disease Heterogeneity: Further characterization of the "unassigned" cell population in advanced tumors could uncover novel immune cell states or therapeutic resistance mechanisms, potentially leading to the identification of new targets for advanced lung cancer.

8. Differential Abundance of T Cell Subsets and ILC2s in Lung Cancer Progression

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents box plots visualizing the proportional abundance of various T cell subsets and Group 2 innate lymphoid cells (ILC2s) across three conditions: "Tumor(early)", "Normal", and "Tumor(adv)" (advanced tumor). The goal is to identify and quantify statistically significant shifts in these immune cell populations during the progression of Lung Cancer. The tool plot_box_for_celltype_population_with_signif_difference was used, applying a p-value cutoff of 0.1 for statistical significance between conditions.

Visual Summary

The box plots illustrate the distribution of celltype proportions for eight distinct lymphoid populations (Tfh, Th22, NK, Th17, Th2, T_Naive, Treg, ILC2) across the "Tumor(early)", "Normal", and "Tumor(adv)" conditions. Statistically significant differences (p < 0.1) are highlighted with connecting lines and p-values.

Biological Interpretation

The observed shifts in T cell subsets and ILC2s reveal a dynamically changing immune microenvironment throughout lung cancer progression.

  1. Early-Stage Immune Dysregulation: The "Tumor(early)" condition is characterized by a significant increase in several helper T cell subsets, including Tfh, Th22, Th17, Th2, and Treg cells, as well as T_Naive cells, compared to "Normal" tissue. This suggests an initial host immune response, possibly attempting to combat the tumor or contributing to a pro-inflammatory, pro-tumorigenic environment. The simultaneous significant decrease in NK cells in early tumor is critical, indicating an early immune evasion mechanism where the innate anti-tumor surveillance is suppressed or diminished. The rise in Tregs is a classic hallmark of immune tolerance induced by developing tumors, allowing them to evade destruction PubMed: Treg cells lung cancer immune suppression.
  2. Progression to Advanced Disease and Immune Shifts: As the tumor progresses to "Tumor(adv)", distinct patterns emerge:
    • Persistent NK Cell Suppression: NK cell proportions remain significantly low in advanced tumors, reinforcing their role in immune evasion throughout disease progression PubMed: NK cell function lung cancer immune evasion.
    • Reduction in Helper T Cells: The proportions of Tfh, Th22, and Th2 cells significantly decrease from "Tumor(early)" to "Tumor(adv)", approaching or even falling below (Th2) "Normal" levels. This could indicate a shift from an initial, albeit dysregulated, immune activation towards immune exhaustion or a more immunosuppressive state in advanced cancer.
    • Sustained T_Naive Expansion: T_Naive cells remain significantly elevated in advanced tumors compared to normal, similar to early tumors. This might suggest a failure of effective T cell priming and activation, leading to an accumulation of uncommitted T cells in the tumor microenvironment.
    • ILC2 Decline in Advanced Disease: The significant decrease of ILC2s in advanced tumor stages compared to both normal and early tumor suggests a specific suppression of these innate lymphoid cells as the disease progresses. ILC2s are involved in type 2 immunity and tissue repair, and their reduction might reflect a changing inflammatory milieu or increased immunosuppression in the advanced tumor.

Clinical or Translational Implications

The differential abundance of these immune cell populations provides potential insights into lung cancer diagnosis, prognosis, and therapeutic strategies:

Immunotherapeutic Targets

9. Macrophage Subpopulation Shifts Across Lung Cancer Stages

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents a population bar plot visualizing the relative proportions of various macrophage subsets (Macrophage (M1), Macrophage (M2A), Macrophage (M2B), Macrophage (M2C), Macrophage (M2D)) across individual samples, grouped by physiological conditions: Normal, Tumor (early), and Tumor (adv) in the lung tissue. Each bar represents 100% of the macrophage population within a given sample, illustrating the internal composition of macrophages. This helps to identify shifts in macrophage polarization that may be associated with lung cancer progression.

Visual Summary

The stacked bar plots display the relative abundances of five distinct macrophage subsets for each sample across the three conditions.

Overall, the visual data indicates a progressive shift in macrophage polarization from a more balanced M1/M2A profile in normal tissue towards an enrichment of M2B and M2C macrophage subsets as lung cancer progresses from early to advanced stages.

Biological Interpretation

Macrophages are highly plastic immune cells that can adopt diverse phenotypes, broadly categorized as M1 (classically activated) or M2 (alternatively activated), with various M2 subtypes (M2A, M2B, M2C, M2D) reflecting specific activation states and functions.

The observed increase in Macrophage (M2B) and Macrophage (M2C) populations in both early and advanced lung tumor microenvironments strongly suggests a shift towards a pro-tumorigenic, immunosuppressive phenotype. This polarization of macrophages into tumor-associated macrophages (TAMs) is a well-established hallmark of cancer, where TAMs foster an environment conducive to tumor growth, invasion, and metastasis, and can inhibit anti-tumor immune responses GeneCards: Macrophage differentiation associated with tumor. Given that 'Lung Epithelial cell' is the resolved tumor-origin cell type, these M2-polarized macrophages likely interact with the tumor cells to facilitate disease progression.

The persistent presence of M1 macrophages even in advanced tumors indicates the complexity and heterogeneity of the tumor microenvironment, where both anti-tumor and pro-tumor forces may coexist. However, the proportional increase of M2-like subsets indicates a shift in the overall balance favoring tumor progression.

Clinical or Translational Implications

The progressive enrichment of M2B and M2C macrophage subsets in lung cancer suggests that targeting these pro-tumorigenic macrophages could be a viable therapeutic strategy.

10. Macrophage Subset Population Shifts Across Lung Cancer Progression

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates the proportional changes of specific macrophage subsets (Mac (M2B), Mac (M2A), and Mac (M2C)) across different lung conditions: "Normal," "Tumor(early)," and "Tumor(adv)" (advanced tumor). The goal is to identify statistically significant shifts in these immune cell populations during lung cancer development and progression, which can provide insights into their potential roles in the tumor microenvironment. Statistical significance was determined using a p-value cutoff of 0.1.

Visual Summary

The box plots illustrate the celltype proportion for three macrophage subsets across the specified conditions, with black dots representing individual sample data points.

Mac (M2B) Proportions:

Mac (M2A) Proportions:

Mac (M2C) Proportions:

Biological Interpretation

Macrophages are critical components of the tumor microenvironment (TME) and are often broadly classified into M1 (pro-inflammatory, anti-tumor) and M2 (anti-inflammatory, pro-tumor) phenotypes, though their polarization is more complex and exists on a spectrum. The observed shifts in macrophage subsets provide crucial insights into their evolving roles during lung cancer progression.

The most striking finding is the significant and progressive increase in Mac (M2B) proportions from normal to early tumor and further to advanced tumor stages. M2 macrophages, including subtypes like M2B, are generally associated with tumor progression, angiogenesis, immune suppression, and tissue remodeling, all of which support cancer growth and metastasis [1]. The rising proportion of Mac (M2B) as the tumor advances suggests their increasing recruitment and/or polarization within the TME, contributing to an immunosuppressive and pro-tumorigenic milieu in lung cancer.

Conversely, the significant decrease in Mac (M2A) proportions across tumor stages (Tumor(early) < Normal, Tumor(adv) < Normal, and Tumor(adv) < Tumor(early)) indicates a diminishing presence of this specific M2 subset as the disease progresses. While M2 macrophages are broadly pro-tumorigenic, different M2 subsets can have distinct functional roles. A decrease in M2A might reflect a shift in the specific types of pro-tumorigenic macrophages or a relative dominance of other subsets like M2B. Some studies suggest M2A can be associated with wound healing and fibrosis, and its reduction might imply a shift away from these specific activities or a selection for other M2 subtypes more directly involved in advanced tumor pathology.

For Mac (M2C), the significant decrease only in advanced tumors compared to normal tissue (p=0.0941) suggests a more specific alteration at later stages of the disease. M2C macrophages are known for their immunosuppressive functions, promoting immune tolerance, tissue repair, and phagocytosis of apoptotic cells [2]. Their reduction in advanced tumors could imply a complex rearrangement of the TME, where other immunosuppressive mechanisms might become dominant, or their specific roles are diminished, perhaps replaced by the more rapidly expanding M2B population.

The overall picture in lung cancer progression appears to be characterized by an accumulation of Mac (M2B) cells, which are strongly implicated in tumor promotion, while Mac (M2A) and to some extent Mac (M2C) populations decline. This suggests a dynamic remodeling of the macrophage compartment within the TME, favoring specific pro-tumorigenic M2 subsets as lung cancer advances.

Clinical or Translational Implications

The distinct and statistically significant changes in macrophage subset proportions across lung cancer stages have several potential clinical and translational implications:

  1. Biomarkers for Disease Progression: The increasing proportion of Mac (M2B) and decreasing proportion of Mac (M2A) could serve as potential biomarkers for lung cancer progression. Monitoring these ratios in patient samples might offer prognostic value, indicating the stage or aggressiveness of the tumor.
  2. Therapeutic Targets: Given the pro-tumorigenic role of M2B macrophages, targeting these cells or their polarization pathways could be a viable therapeutic strategy in lung cancer [3]. Interventions aimed at depleting M2B cells, reprogramming them towards an anti-tumor M1 phenotype, or inhibiting their recruitment could potentially slow tumor growth and enhance anti-tumor immunity.
  3. Immunotherapy Response Prediction: The specific balance of macrophage subsets in the TME is known to influence the response to various immunotherapies. Understanding these shifts could help predict which patients might benefit most from certain treatments or guide the development of combination therapies that modulate macrophage polarization alongside other immunomodulatory agents.
  4. Novel Drug Development: The differential changes in M2A, M2B, and M2C populations highlight the complexity of macrophage phenotypes in cancer. Developing drugs that specifically target the pathways driving the accumulation of pro-tumorigenic M2B macrophages, or prevent the decrease of potentially beneficial macrophage subtypes, could offer more precise and effective cancer treatments.

References

  1. M2 Macrophages and Tumor Progression:

PubMed Search: "M2 macrophages tumor progression"

  1. M2C Macrophage Functions:

PubMed Search: "M2c macrophages immune suppression"

  1. Targeting Macrophages in Cancer Therapy:

PubMed Search: "targeting M2 macrophages cancer therapy"

11. Ploidy Analysis of Tumor-Origin Lung Epithelial Cells

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis examines the ploidy status (Diploid, Aneuploid, Unclear) of Lung Epithelial cells, identified as the tumor-origin cell type, across samples from three conditions: Normal, Tumor (advanced), and Tumor (early). The objective is to understand how chromosomal instability, reflected by aneuploidy, correlates with disease state in lung cancer. The results are presented as stacked bar plots, where each bar represents a single sample, showing the proportional distribution of ploidy populations within its Lung Epithelial cells.

Visual Summary

The stacked bar plots effectively demonstrate distinct ploidy patterns in Lung Epithelial cells across the different conditions and individual samples.

Biological Interpretation

The observed ploidy changes in Lung Epithelial cells provide strong biological insights into lung cancer development and progression, consistent with their identification as the tumor-origin cell type.

Clinical or Translational Implications

12. Cell-Cell Interaction Patterns in Lung Cancer Microenvironment

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates cell-cell interaction (CCI) patterns involving key cell types within the lung tissue microenvironment across normal, early-stage tumor, and advanced-stage tumor conditions. The primary focus is on interactions concerning the identified tumor-origin cells (Lung Epithelial cells, specifically aneuploid), alongside immune cells (T cells, B cells, Macrophages, Dendritic cells) and stromal components (Fibroblasts, Endothelial cells). The goal is to identify significant ligand-receptor mediated communications that characterize tumor development and progression, highlighting potential therapeutic targets.

Visual Summary

The provided dot plots illustrate the strength (color intensity, log2(mean expression)) and statistical significance (dot size, -log10(p-value)) of various cell-cell interactions. Each plot represents a distinct condition: Normal, Tumor (early), and Tumor (adv).

Condition-Specific Patterns:

Biological Interpretation

1. Tumor-Origin Cell Interactions (Aneuploid Lung Epithelial Cells)

Tumor-Stromal Crosstalk:

Tumor-Immune Cell Interactions:

2. Macrophage-Specific Interactions

3. T Cell and Other Immune Cell Interactions

Clinical or Translational Implications

The differential cell-cell interaction patterns identified in this analysis offer several promising avenues for therapeutic targeting and deeper investigation in Lung Cancer:

Immunomodulation:

13. Condition-Specific Cell-Cell Interaction Patterns in Lung Cancer

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates condition-specific cell-cell interactions (CCIs) involving Lung Epithelial cells (identified as the tumor-origin cell type) and major immune and stromal cells across Normal, Tumor(early), and Tumor(adv) conditions in lung tissue samples. The goal is to identify significantly different interaction patterns that may underpin disease progression and reveal potential therapeutic targets. The CellPhoneDB tool was used to predict ligand-receptor interactions, and the results are visualized as a dot plot, showing interaction strength (color intensity) and significance (dot size) for each sample and condition.

Visual Summary

The dot plot visualizes a comprehensive landscape of cell-cell interactions, faceted by condition: Normal, Tumor(adv), and Tumor(early).

Biological Interpretation

The differential cell-cell interaction patterns highlight key biological processes disrupted or activated in lung cancer, providing insights into tumor-microenvironment co-evolution.

  1. Tumor Cell-Mediated Immunosuppression:
    • Prostaglandin E2 (PGE2) Signaling: The strong interaction of ProstaglandinE2_byPTGES3_PTGER2 from Lung.Epi (Aneupl) to T CD8+ cells (and other immune cells via PTGER4 on B cells and T CD4+ cells) is a critical finding. PGE2, often produced by tumor cells (indicated by PTGES3), is a potent immunosuppressant in the tumor microenvironment. It can inhibit the activation and proliferation of CD8+ T cells, promote the differentiation of regulatory T cells (Tregs), and influence myeloid-derived suppressor cells (MDSCs), ultimately contributing to immune evasion by the tumor. PubMed search for "PGE2 lung cancer immune suppression"
    • Retinoic Acid Signaling: The observed atRetinoicAcid_byALDH1A1_RAReceptor_RXRB -- Lung.Epi (Aneupl) | B cell and atRetinoicAcid_byALDH1A1_RAReceptor_RXRA -- Lung.Epi (Aneupl) | Mac interactions suggest that aneuploid lung epithelial cells (tumor cells) expressing ALDH1A1 might be producing retinoic acid (RA). RA can profoundly influence immune cell differentiation and function, often promoting tolerogenic immune responses by inducing Tregs or impairing effector T cell function, thus contributing to an immunosuppressive environment. PubMed search for "retinoic acid lung cancer immune microenvironment ALDH1A1"
  2. Angiogenesis and Extracellular Matrix (ECM) Remodeling:
    • VEGFB-NRP1 Axis: The prominent VEGFB_NRP1 -- Lung.Epi (Aneupl) | Endo interaction indicates active communication between tumor cells and endothelial cells. While VEGFA is the primary driver of pathological angiogenesis, VEGFB, often with Neuropilin-1 (NRP1), also contributes to vascular development and tumor growth. NRP1, expressed on endothelial cells, acts as a co-receptor for various growth factors, modulating angiogenesis and also influencing immune cell function, including Treg activity. GeneCards for NRP1
    • Collagen-Integrin Interactions: Numerous COL1A1_integrin_a1b1_complex -- Lung.Epi (Aneupl) | Fib and other collagen-integrin interactions highlight the extensive remodeling of the extracellular matrix (ECM) in the tumor microenvironment. Cancer-associated fibroblasts (CAFs) are key players in producing ECM components like Collagen I (COL1A1). Integrins mediate cell-ECM adhesion and signal transduction, which are crucial for tumor cell proliferation, survival, migration, and invasion. This intricate interplay creates a pro-tumorigenic niche. GeneCards for COL1A1
  3. Autocrine/Paracrine Tumor Growth Signaling:
    • TGFB1 Signaling: The TGFB1_TGFbeta_receptor1 -- Lung.Epi (Aneupl) | Lung.Epi (Aneupl) interaction points to autocrine or paracrine TGF-beta signaling within the tumor cell population. Transforming Growth Factor-beta (TGF-beta) is a pleiotropic cytokine that can promote tumor cell growth, epithelial-mesenchymal transition (EMT), immune suppression, and metastasis in advanced cancers. This self-sustaining signaling loop is a hallmark of many solid tumors. PubMed search for "TGF-beta lung cancer progression"
  4. Altered Adhesion and Immune Evasion:
    • ICAM1-Integrin: The reduced ICAM1_integrin_aLb2_complex -- Lung.Epi (Dip) | Mac in tumor conditions compared to normal, alongside the emergence of ICAM1_integrin_aLb2_complex -- Lung.Epi (Aneupl) | T CD8+ in tumor, suggests a shift in cell adhesion dynamics. Altered ICAM1 expression on tumor cells can impact immune cell recognition and cytotoxic activity, contributing to immune evasion.

Clinical or Translational Implications

The identified condition-specific CCIs represent promising avenues for therapeutic intervention in lung cancer.

Targeting Immunosuppressive Pathways:

Inhibiting Angiogenesis and ECM Remodeling:

Blocking Autocrine/Paracrine Tumor Growth:

These interactions provide a molecular blueprint of the lung cancer microenvironment, offering testable hypotheses for combination therapies and precision medicine approaches aimed at disrupting tumor-supporting communication networks.

14. Condition-Specific Cell-Cell Interaction Patterns in Lung Cancer

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates condition-specific cell-cell interaction (CCI) patterns involving Lung Epithelial cells (identified as tumor-origin), major immune cells, and stromal cells in lung tissue from Normal, Tumor (early), and Tumor (advanced) conditions. The focus is on interactions mediated by genes associated with immune checkpoints and cell cycle pathways. The results are visualized as a dot plot, where the size of the dot represents the significance (-log10(p-value)) and the color intensity represents the scaled interaction strength. The plot displays interactions that differ significantly between conditions, highlighting key communication pathways within the tumor microenvironment (TME).

Visual Summary

The dot plot effectively illustrates distinct patterns of cell-cell communication across Normal, Tumor (advanced), and Tumor (early) conditions.

  1. Condition-Specific Clustering: Samples largely cluster by condition, with Normal samples exhibiting a unique set of strong interactions (top-left blue box), while Tumor (advanced) and Tumor (early) samples share a highly overlapping and distinct pattern of interactions (middle-right blue box). This suggests a fundamental shift in cell communication in the cancerous state, which is largely established by the early tumor stage.
  2. Dominance of Aneuploid Lung Epithelial Cells in Tumor: Interactions prominent in tumor conditions (Tumor(adv) and Tumor(early)) overwhelmingly involve Lung.Epi (Aneuploid) cells, confirming their central role as the tumor-origin cell type and highlighting the involvement of malignant cells in these dysregulated interactions. In contrast, normal-specific interactions primarily involve Lung.Epi (Dip) (Diploid Lung Epithelial cells).
  3. Key Pathways and Cell Types in Tumor:
    • EGFR Signaling: Interactions involving EGFR ligands (AREG, EREG, HBEGF, BTC, TGFA) and the EGFR receptor are highly prevalent and strong in both tumor conditions. These involve Lung.Epi (Aneuploid) interacting with itself (autocrine/paracrine stimulation) as well as with Macrophages (Mac) and Fibroblasts (Fib).
    • TGF-beta Signaling: TGFB1 and TGFB3 interactions with TGF-beta receptors (TGFBeta_receptor1, TGFBeta_receptor2, TGFBR3) and integrin complexes (integrin_aVb6_complex) are also highly active in tumor samples. These primarily involve Lung.Epi (Aneuploid) interacting with Macrophages, Fibroblasts, and various immune cells (T CD8+, B cell, ILC).
    • Interferon (IFN) Signaling: While some IFN-gamma (IFNG_Type_II_IFNR) interactions are present in normal samples (e.g., with NK cells and T CD8+), they persist and often strengthen in tumor conditions with Lung.Epi (Aneuploid). Notably, CD93_IFNGR1 interactions involving Endothelial cells (Endo) and Macrophages with Lung.Epi (Aneuploid) also become prominent in tumor.
    • Immune Checkpoint/Co-stimulatory Interactions: Interactions like CD80 CD28 -- Lung.Epi (Aneuploid) | B cell and LCK_CD8_receptor -- Lung.Epi (Aneuploid) | ILC emerge as significant in tumor samples, suggesting altered immune regulation.
  4. No Major Differences Between Early and Advanced Tumor: For the selected top significant interactions, there doesn't appear to be a drastic qualitative or quantitative difference in the interaction patterns between Tumor (early) and Tumor (advanced) samples. The core dysregulated communication network seems to be largely established early in tumor progression.

Biological Interpretation

The observed CCI patterns reveal critical shifts in the lung tumor microenvironment, emphasizing the complex interplay between malignant epithelial cells, immune cells, and stromal cells.

EGFR Signaling in Tumor Progression

Key Tumor Interactions:

TGF-beta Pathway Remodeling the Tumor Microenvironment

Key Tumor Interactions:

Interferon Signaling and Immune Context

Key Tumor Interactions:

Immune Checkpoint and Co-stimulation Anomalies

Key Tumor Interactions:

Clinical or Translational Implications

The identified condition-specific cell-cell interactions, particularly those involving aneuploid lung epithelial cells in tumor samples, highlight several promising avenues for therapeutic intervention in Lung Cancer.

Targeting EGFR Signaling

Modulating TGF-beta Pathway

Leveraging or Counteracting Interferon Signaling

Novel Immune Checkpoint Strategies

15. Surfaceome Marker Expression in Lung Epithelial Cells of Lung Cancer

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis focuses on identifying and characterizing surfaceome markers expressed by Lung Epithelial cells from single-cell RNA-seq data. The goal is to delineate expression patterns across different conditions (Normal, Tumor(early), Tumor(adv)) and identify potential therapeutic targets. Markers were selected based on differential expression against a 'Normal' reference condition, specifically considering only cell surface proteins (surfaceome markers) with stringent fold change and non-zero percentage cutoffs. The interpretation distinguishes between established therapeutic targets and less well-characterized markers with potential clinical utility, providing biological context and relevant therapeutic information.

Visual Summary

The dot plot provides a comprehensive visualization of surfaceome gene expression and prevalence across various lung epithelial cell samples, grouped by condition and patient origin.

Biological Interpretation

The identification of these surfaceome markers offers critical insights into the molecular mechanisms driving lung cancer in epithelial cells. These proteins are often involved in cell communication, adhesion, and signaling, making them attractive candidates for therapeutic intervention.

Established Therapeutic Targets

  1. EGFR (Epidermal Growth Factor Receptor)
    • Expression Pattern: EGFR is broadly and highly expressed in many tumor samples (both early and advanced), with lower levels in normal lung epithelial cells. This is consistent with its known role in lung cancer.
    • Signaling Pathways & Mechanisms: EGFR is a receptor tyrosine kinase crucial for cell proliferation, survival, migration, and differentiation. Upon ligand binding, it activates multiple downstream pathways, including the RAS/RAF/MEK/ERK cascade (mitogenesis) and the PI3K/AKT/mTOR pathway (survival and growth). Aberrant activation, often due to activating mutations or gene amplification, is a primary oncogenic driver in non-small cell lung cancer (NSCLC).
    • Relevant Therapeutic Agents: Multiple EGFR tyrosine kinase inhibitors (TKIs) are approved for EGFR-mutant NSCLC, such as gefitinib, erlotinib, afatinib, osimertinib, and dacomitinib. Monoclonal antibodies targeting EGFR, like cetuximab and necitumumab, are also used in certain cancer contexts.
    • Supporting Evidence: EGFR mutations are a well-established and actionable biomarker in NSCLC, guiding the use of EGFR-TKIs as first-line therapy for eligible patients PubMed Search: EGFR lung cancer targeted therapy.
    • GeneCards: https://www.genecards.org/cgi-bin/carddisp.pl?gene=EGFR
  2. MET (Mesenchymal-Epithelial Transition factor)
    • Expression Pattern: Upregulated in a subset of tumor samples, notably in both early and advanced stages, showing minimal expression in most normal cells.
    • Signaling Pathways & Mechanisms: MET is a receptor tyrosine kinase that, when activated by its ligand hepatocyte growth factor (HGF), initiates intracellular signaling cascades (e.g., PI3K/AKT, RAS/MAPK, STAT3) that promote cell proliferation, survival, motility, invasion, angiogenesis, and resistance to apoptosis. Dysregulation of MET, through gene amplification, exon 14 skipping mutations, or protein overexpression, can drive tumorigenesis and metastasis in NSCLC.
    • Relevant Therapeutic Agents: Several MET-targeting TKIs, including capmatinib, tepotinib, and savolitinib, are approved or in late-stage clinical development for NSCLC patients with MET exon 14 skipping mutations or MET amplification.
    • Supporting Evidence: MET alterations are recognized oncogenic drivers in NSCLC, with specific inhibitors demonstrating clinical efficacy in patients harboring these genetic changes PubMed Search: MET inhibition lung cancer.
    • GeneCards: https://www.genecards.org/cgi-bin/carddisp.pl?gene=MET
  3. ERBB2 (HER2 - Erb-B2 Receptor Tyrosine Kinase 2)
    • Expression Pattern: Shows increased expression in a distinct subset of tumor samples, particularly in advanced disease, with low expression in normal cells.
    • Signaling Pathways & Mechanisms: ERBB2, also known as HER2, is a member of the EGFR family that acts as a co-receptor, enhancing the signaling of other ERBB family members. While lacking a direct ligand, it forms heterodimers that activate key pro-survival and pro-proliferative pathways such as PI3K/AKT/mTOR and RAS/MAPK. ERBB2 amplification or activating mutations are oncogenic drivers in various cancers, including a small but significant fraction of NSCLC cases.
    • Relevant Therapeutic Agents: Antibody-drug conjugates (ADCs) like trastuzumab deruxtecan (T-DXd) are effective in HER2-mutant NSCLC. Other agents, including HER2-targeted TKIs and monoclonal antibodies, are also under investigation.
    • Supporting Evidence: HER2 alterations, though less common than EGFR mutations, represent an important therapeutic target in NSCLC, with targeted therapies showing promising clinical responses PubMed Search: HER2 mutated lung cancer treatment.
    • GeneCards: https://www.genecards.org/cgi-bin/carddisp.pl?gene=ERBB2

Less Well-Characterized Markers with Potential as Therapeutic Targets

  1. CEACAM5 (Carcinoembryonic Antigen Related Cell Adhesion Molecule 5)
    • Expression Pattern: Exhibits strong upregulation in a large fraction of tumor cells in both early and advanced stages, largely absent in normal cells.
    • Signaling Pathways & Mechanisms: CEACAM5 is a cell adhesion glycoprotein involved in cell-cell interactions. In cancer, its overexpression can promote cell proliferation, survival, and immune evasion by modulating various signaling pathways. It can also contribute to resistance to chemotherapy.
    • Relevant Therapeutic Agents: CEACAM5 is being actively pursued as a target for antibody-drug conjugates (ADCs) due to its high and selective expression on cancer cells. Tusamitamab ravtansine, an anti-CEACAM5 ADC, is in clinical trials for NSCLC.
    • Supporting Evidence: CEACAM5 is a well-known tumor-associated antigen frequently overexpressed in various adenocarcinomas, including NSCLC, and its role as a therapeutic target is supported by preclinical and clinical studies with ADCs PubMed Search: CEACAM5 lung cancer ADC clinical trial.
    • GeneCards: https://www.genecards.org/cgi-bin/carddisp.pl?gene=CEACAM5
  2. EPHA2 (Ephrin Receptor A2)
    • Expression Pattern: Shows increased expression in tumor epithelial cells, particularly in advanced stages, with minimal to no expression in normal cells.
    • Signaling Pathways & Mechanisms: EPHA2 is a receptor tyrosine kinase that interacts with ephrin-A ligands. Its role in cancer is complex, but often, its overexpression and aberrant activation promote cell proliferation, survival, migration, and angiogenesis. It can activate downstream signaling pathways like Rho GTPases and MAPK, influencing cytoskeletal dynamics and cell motility, thereby facilitating metastasis.
    • Relevant Therapeutic Agents: Monoclonal antibodies, ADCs, and small molecule inhibitors targeting EPHA2 are under preclinical and early clinical investigation for various cancers, including lung cancer.
    • Supporting Evidence: EPHA2 overexpression is frequently observed in lung cancer and correlates with disease progression, metastasis, and poor patient prognosis, highlighting its potential as a therapeutic target PubMed Search: EPHA2 lung cancer therapeutic.
    • GeneCards: https://www.genecards.org/cgi-bin/carddisp.pl?gene=EPHA2
  3. PLAUR (Plasminogen Activator, Urokinase Receptor)
    • Expression Pattern: Consistently upregulated in tumor cells from both early and advanced stages, showing low expression in normal epithelial cells.
    • Signaling Pathways & Mechanisms: PLAUR is a GPI-anchored cell surface receptor that binds urokinase plasminogen activator (uPA), initiating a proteolytic cascade that degrades the extracellular matrix. This process is critical for cell migration, invasion, and metastasis. PLAUR also signals directly or indirectly through integrins and other receptors, affecting cell adhesion, proliferation, and survival. Its overexpression promotes an aggressive tumor phenotype.
    • Relevant Therapeutic Agents: Inhibitors of the uPA/uPAR system and monoclonal antibodies targeting PLAUR are being investigated in preclinical and early clinical studies to block tumor invasion and metastasis.
    • Supporting Evidence: High PLAUR expression is strongly associated with increased tumor invasiveness, metastatic potential, and poor prognosis across a wide range of cancers, including NSCLC PubMed Search: PLAUR lung cancer invasion.
    • GeneCards: https://www.genecards.org/cgi-bin/carddisp.pl?gene=PLAUR
  4. MUC4 (Mucin 4, Cell Surface Associated)
    • Expression Pattern: Clearly overexpressed in tumor epithelial cells across both early and advanced conditions, with minimal to absent expression in normal cells.
    • Signaling Pathways & Mechanisms: MUC4 is a large transmembrane mucin glycoprotein that forms a protective layer on epithelial surfaces. In cancer, its overexpression contributes to tumor progression by masking tumor antigens, promoting cell proliferation and survival (often through interaction with ERBB2), enhancing cell motility, and conferring resistance to apoptosis and chemotherapy. It activates downstream signaling pathways like MAPK and PI3K/AKT.
    • Relevant Therapeutic Agents: Therapeutic strategies include monoclonal antibodies targeting MUC4 or ADCs. Disruption of MUC4-ERBB2 interactions is also a focus of drug development.
    • Supporting Evidence: MUC4 overexpression is a common feature in lung adenocarcinoma and is often associated with poor prognosis, metastasis, and resistance to therapy, making it a promising candidate for targeted therapeutic development PubMed Search: MUC4 lung adenocarcinoma prognosis.
    • GeneCards: https://www.genecards.org/cgi-bin/carddisp.pl?gene=MUC4

Clinical or Translational Implications

This analysis provides a valuable roadmap for advancing lung cancer diagnostics and therapeutics:

In conclusion, this analysis provides compelling evidence for several surfaceome proteins as actionable targets in lung cancer, ranging from well-validated drivers to promising novel candidates, offering significant potential for future clinical translation.

16. Macrophage Condition-Specific Surfaceome Markers in Lung Cancer

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis identifies condition-specific surfaceome markers for Macrophages across Normal, advanced Tumor (Tumor(adv)), and early Tumor (Tumor(early)) lung tissue samples. The dot plot visualizes the mean expression level (color intensity) and the fraction of cells expressing (dot size) each marker within different samples. The goal is to highlight macrophage markers that are differentially expressed in the tumor microenvironment (TME) compared to normal tissue, providing insights into their biological roles and potential as therapeutic targets in Lung Cancer.

Visual Summary

The dot plot displays the expression of 74 macrophage surfaceome markers across 32 individual samples, categorized by their condition (Normal, Tumor(adv), Tumor(early)). Each row represents a sample, and each column represents a marker gene. The total number of cells per sample is indicated on the right-hand side.

Key Observations:

Biological Interpretation

Macrophages are key components of the tumor microenvironment (TME) in Lung Cancer, with diverse functions that can either promote or suppress tumor growth. The observed differential expression of surface markers indicates distinct macrophage phenotypes in normal versus tumor conditions, reflecting their adaptation and functional specialization within the TME.

Markers Enriched in Normal Macrophages

The high expression of specific surface markers in normal lung macrophages, and their downregulation in tumor-associated macrophages (TAMs), points to a shift away from a homeostatic or healthy inflammatory phenotype.

TREM1 (Triggering Receptor Expressed on Myeloid cells 1):

CD46 (Membrane Cofactor Protein):

Markers Enriched in Tumor-Associated Macrophages (TAMs)

The upregulation of these markers in TAMs indicates their active participation in tumor progression, often through immune suppression, angiogenesis, and extracellular matrix remodeling.

CSF1R (Colony Stimulating Factor 1 Receptor):

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

MMP14 (Matrix Metallopeptidase 14 / MT1-MMP):

TNFRSF1B (Tumor Necrosis Factor Receptor Superfamily Member 1B / TNFR2):

CD14 (Lipopolysaccharide Receptor):

Clinical or Translational Implications

The identification of condition-specific surfaceome markers on macrophages provides valuable insights with significant clinical and translational implications for Lung Cancer.

Diagnostic and Prognostic Biomarkers:

Therapeutic Targets for Macrophage Reprogramming/Depletion:

Personalized Medicine Approaches:

Monitoring Treatment Response:

In summary, the distinct surfaceome profiles of macrophages in normal versus tumor conditions highlight the dynamic role of these cells in Lung Cancer progression. Targeting these specific surface markers represents a promising avenue for developing novel immunotherapeutic strategies to reprogram the tumor microenvironment and improve patient outcomes.

17. Fibroblast Condition-Specific Surfaceome Markers in Lung Cancer

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis identifies condition-specific surfaceome markers for Fibroblasts in human lung tissue, comparing Normal samples to early-stage Tumor (Tumor(early)) samples. The markers were selected using differential gene expression analysis against a reference (Normal), with a focus on surface-expressed proteins. 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 condition. This approach aims to reveal phenotypic shifts in Fibroblast populations within the tumor microenvironment (TME) and identify potential therapeutic targets.

Visual Summary

The dot plot clearly delineates two distinct populations of Fibroblast markers based on their expression profiles across Normal and Tumor(early) conditions.

  1. Normal Condition-Specific Markers (e.g., GPNMB, TSPAN8, SCARA5, GAS1, LEPR, PTPRS, CD34, CDH11): These markers show high mean expression (darker red dots) and a high fraction of expressing cells (larger dot sizes) in the Fibroblast populations from Normal lung samples (LUNG_N31, LUNG_N18, LUNG_N30, LUNG_N34, LUNG_N09, LUNG_N19). Conversely, their expression is largely absent or significantly reduced in Tumor(early) samples, indicated by smaller, lighter-colored, or absent dots. This suggests a unique signature for homeostatic or non-tumor-activated fibroblasts.
  2. Tumor(early) Condition-Specific Markers (e.g., PDGFRB, ANTXR1, MMP14, FAP, ITGAV, NOTCH2, VCAM1, FZD1): A broad array of surface markers demonstrates strong and prevalent expression in Fibroblasts from early tumor samples (LUNG_T28, LUNG_T06, LUNG_T08, LUNG_T19, LUNG_T18, LUNG_T31, LUNG_T30) but are largely absent or minimally expressed in Normal samples. These markers include classic indicators of activated cancer-associated fibroblasts (CAFs), such as FAP and PDGFRB, alongside various matrix metalloproteinases, integrins, and signaling receptors. Their robust and consistent upregulation in early tumor stages highlights a significant phenotypic reprogramming of fibroblasts in the presence of malignancy.

Overall, the plot effectively illustrates a dramatic shift in fibroblast surfaceome composition as these cells transition from a normal tissue state to an early tumor-associated state, underscoring their active role in shaping the evolving tumor microenvironment.

Biological Interpretation

The observed differential expression of surfaceome markers in Fibroblasts between normal and early tumor conditions provides critical insights into their activation state and function within the lung cancer microenvironment. This shift indicates the emergence of Cancer-Associated Fibroblasts (CAFs), which are known to significantly contribute to tumor progression.

Below, we detail key markers prominently expressed in Tumor(early) Fibroblasts:

FAP (Fibroblast Activation Protein Alpha)

PDGFRB (Platelet-Derived Growth Factor Receptor Beta)

MMP14 (Matrix Metalloproteinase 14 / MT1-MMP)

Integrins (e.g., ITGAV, ITGB5, ITGA1, ITGA8)

VCAM1 (Vascular Cell Adhesion Molecule 1)

NOTCH2

Clinical or Translational Implications

The identification of distinct surfaceome marker profiles in Fibroblasts across normal and early tumor stages holds significant clinical and translational value for Lung Cancer.

  1. Biomarker Discovery for Early Detection and Prognosis: The specific upregulation of markers like FAP, PDGFRB, and MMP14 in early tumor-associated fibroblasts could serve as early diagnostic biomarkers for lung cancer progression or as prognostic indicators for patient outcomes. Detection of these markers, perhaps via liquid biopsy or specific imaging agents, could enable earlier intervention.
  2. Therapeutic Targeting of CAFs: The prominent and specific expression of these surface markers on CAFs, with limited expression in normal cells, presents opportunities for targeted therapies.
    • Direct Targeting: Agents like FAP-targeted antibodies, CAR T-cells, or radioconjugates could specifically eliminate or inactivate CAFs, thereby disrupting the pro-tumorigenic microenvironment.
    • Pathway Inhibition: Targeting key signaling pathways active in CAFs (e.g., PDGFRB with TKIs, Notch with GSIs) could inhibit CAF activation, proliferation, and their support for tumor growth.
    • ECM Remodeling: Inhibiting enzymes like MMP14 could reduce matrix degradation and hinder tumor invasion and metastasis, potentially making the tumor more accessible to other therapies.
  3. Enhancing Immunotherapy Efficacy: CAFs are known to contribute to immunosuppression in the TME. By targeting CAF-specific markers (e.g., FAP, VCAM1), it may be possible to reprogram the TME, reduce immunosuppression, and enhance the efficacy of immune checkpoint inhibitors or other immunotherapies in lung cancer patients.
  4. Monitoring Treatment Response: Changes in the expression of these CAF markers could potentially be used to monitor treatment response and detect residual disease or recurrence, providing a dynamic assessment of therapeutic interventions.

Overall, these findings underscore the critical role of Fibroblasts in lung cancer development and highlight a rich landscape of surfaceome targets for novel diagnostic and therapeutic strategies, particularly in the context of early disease intervention.

18. Condition-Specific Surfaceome Markers in CD4+ T Cells in Lung Cancer

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis identifies and visualizes surfaceome markers in CD4+ T cells, comparing expression patterns across Normal, Tumor (advanced), and Tumor (early) conditions in Lung tissue. The dot plot displays the mean expression level (color intensity) and the fraction of cells expressing each marker (dot size) for individual samples grouped by their condition. Markers were selected based on differential expression against a 'Normal' reference condition, focusing on surfaceome genes and applying specific cutoffs for non-zero percentage score and fold change, with at most 40 markers per condition considered. The analysis aims to uncover distinct immunological states of CD4+ T cells in the tumor microenvironment (TME) of Lung Cancer.

Visual Summary

The dot plot, titled "T cell CD4+ condition-specific markers," effectively categorizes samples into Normal, Tumor(adv), and Tumor(early) conditions, showcasing distinct marker expression profiles within each group.

Biological Interpretation

The identified surfaceome markers provide insights into the functional states and interactions of CD4+ T cells in the context of Lung Cancer. The distinctions between normal and tumor-associated CD4+ T cells highlight potential mechanisms of immune regulation and dysregulation in the TME.

Markers Enriched in Normal CD4+ T Cells (Compared to Tumor Conditions)

These markers suggest an active role in immune surveillance and maintenance of homeostasis in healthy lung tissue. Their downregulation in tumor conditions may indicate immune evasion or T cell dysfunction.

HLA-DRA, HLA-DRB1, HLA-DQA1 (MHC Class II molecules)

ADGRE5 (CD97)

MYADM (Myeloid-associated differentiation marker)

CD27

PIK3IP1 (FOAP1)

SELL (CD62L)

CD164

Markers Enriched in Tumor-Associated CD4+ T Cells (Compared to Normal Condition)

These markers indicate specific activation states, exhaustion, or regulatory functions acquired by CD4+ T cells within the immunosuppressive TME of Lung Cancer.

CD7

BST1 (CD157)

TNFRSF14 (HVEM)

EVI2A

ICOS (CD278)

Clinical or Translational Implications

The distinct surfaceome marker profiles observed in CD4+ T cells from Normal versus Tumor (early and advanced) Lung Cancer patients offer valuable clinical and translational insights:

  1. Biomarker Development: The identified markers can serve as potential diagnostic or prognostic biomarkers for Lung Cancer. For instance, a decrease in MHC Class II molecules (HLA-DRA/DRB1/DQA1) or CD27 on circulating or tumor-infiltrating CD4+ T cells could indicate tumor presence or disease progression. Conversely, elevated expression of *ICOS* or *TNFRSF14* might identify distinct populations of tumor-infiltrating lymphocytes with specific functional states.
  2. Therapeutic Target Identification: Markers predominantly expressed in tumor-associated CD4+ T cells, such as *ICOS* and *TNFRSF14*, represent promising immunotherapeutic targets.
    • ICOS agonists: Given ICOS's role in T cell activation and proliferation, agonistic antibodies could be explored to enhance the anti-tumor response of CD4+ T cells, potentially synergizing with existing checkpoint inhibitors.
    • TNFRSF14 (HVEM) modulators: The complex role of TNFRSF14 in immune co-stimulation and co-inhibition makes it an intriguing target. Therapies could aim to either activate pro-inflammatory pathways (e.g., via LIGHT-HVEM) or block inhibitory interactions (e.g., via BTLA-HVEM) to reshape the TME.
  3. Understanding Immune Evasion: The downregulation of MHC Class II molecules in tumor CD4+ T cells suggests a mechanism of immune evasion or altered antigen presentation capacity in the TME. Strategies to restore or enhance antigen presentation could improve the efficacy of immunotherapies.
  4. T Cell Phenotyping: These markers provide a refined understanding of CD4+ T cell heterogeneity within the Lung Cancer TME. This can guide the development of adoptive T cell therapies by selecting or engineering specific T cell subsets with desired functional properties. For example, identifying the ratio of CD27-high (less exhausted) vs. ICOS-high (activated/regulatory) T cells could inform treatment strategies.
  5. Monitoring Treatment Response: Changes in the expression of these markers on CD4+ T cells could potentially be used to monitor patient responses to various cancer treatments, including chemotherapy, radiation, and immunotherapy.

19. Differential Expression of Cell Cycle Genes in Lung Epithelial Cells Across Lung Cancer Stages

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates the expression of a curated panel of cell cycle pathway genes within Lung Epithelial cells across different conditions: 'Normal', 'Tumor(early)', and 'Tumor(adv)' (advanced tumor). The expression is quantified as the "expressing cell fraction (sample)", which represents the proportion of cells within each sample that express a given gene. The objective is to identify statistically significant differences in the expression of these genes across these conditions, providing insight into cell cycle dysregulation in lung cancer.

Visual Summary

The box plots display the distribution of expressing cell fractions for 24 selected cell cycle-related genes across the three conditions. Each dot represents a sample. Key observations include:

Biological Interpretation

The observed differential expression patterns of cell cycle genes in Lung Epithelial cells provide strong biological insights into lung cancer pathogenesis:

Compromised Cell Cycle Checkpoints and Tumor Suppression:

Dysregulation of Ubiquitin-Proteasome System and Signaling:

Clinical or Translational Implications

Therapeutic Targets:

References

  1. Hallmarks of Cancer: https://pubmed.ncbi.nlm.nih.gov/?term=Hanahan+Weinberg+hallmarks+of+cancer
  2. GADD45B Function: https://www.genecards.org/cgi-bin/carddisp.pl?gene=GADD45B
  3. TP53 in Cancer: https://www.genecards.org/cgi-bin/carddisp.pl?gene=TP53
  4. SMAD4 in Cancer: https://www.genecards.org/cgi-bin/carddisp.pl?gene=SMAD4
  5. CUL1 Function: https://www.genecards.org/cgi-bin/carddisp.pl?gene=CUL1
  6. MDM2 Function: https://www.genecards.org/cgi-bin/carddisp.pl?gene=MDM2
  7. TGFB1 in Cancer: https://pubmed.ncbi.nlm.nih.gov/?term=TGFB1+lung+cancer+dual+role

20. Gene Ontology Analysis of Lung Epithelial Cells Across Conditions

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents Gene Ontology (GO) enrichment results, specifically Gene Set Analysis (GSA), for Lung Epithelial cells. The analysis was performed on up-regulated differentially expressed genes (DEGs) in various conditions compared to a reference 'Normal' state. Four distinct comparisons are examined:

  1. Diploid vs. Normal: Comparing lung epithelial cells identified as diploid against normal lung epithelial cells.
  2. Normal vs. others: Comparing normal lung epithelial cells against all other cell conditions (Tumor(adv), Tumor(early), Diploid). This highlights pathways characteristic of healthy tissue.
  3. Tumor(adv) vs. Normal: Comparing advanced-stage tumor lung epithelial cells against normal lung epithelial cells.
  4. Tumor(early) vs. Normal: Comparing early-stage tumor lung epithelial cells against normal lung epithelial cells.

The results are displayed as bar plots, where terms are ranked by their statistical significance (-log(p-val) and -log(q-val)).

Visual Summary

The four bar plots illustrate distinct patterns of pathway enrichment across different conditions of Lung Epithelial cells. Each plot ranks GO terms based on the negative logarithm of their p-value and q-value (FDR-adjusted p-value) for up-regulated genes.

  1. GSA for Lung Epithelial cell: Diploid_vs_Normal: This plot shows enrichment for "Coronavirus disease," "Ribosome," "Glycolysis / Gluconeogenesis," "Glutathione metabolism," and "Protein processing in endoplasmic reticulum" among the top terms. The -log(p-val) values generally range up to 3.5, and -log(q-val) values are relatively low, indicating moderate significance.
  2. GSA for Lung Epithelial cell: Normal_vs_others: This plot, contrasting normal cells against all others, shows "Phagosome," "PPAR signaling pathway," "Human T-cell leukemia virus 1 infection," "Type I diabetes mellitus," and "Allograft rejection" as leading terms. The -log(p-val) values reach up to 4.5, with corresponding positive -log(q-val) values, suggesting pathways enriched in the normal state.
  3. GSA for Lung Epithelial cell: Tumor(adv)_vs_Normal: This comparison reveals highly significant enrichment for "Ribosome," "Spliceosome," "Pathways of neurodegeneration," "Coronavirus disease," "Huntington disease," "Parkinson disease," "RNA transport," "Alzheimer disease," and "Proteasome." The -log(p-val) and -log(q-val) scales are notably higher (up to 20 for p-val, 17 for q-val), indicating extremely strong statistical significance for these pathways in advanced tumor cells.
  4. GSA for Lung Epithelial cell: Tumor(early)_vs_Normal: Similar to the advanced tumor state, this plot also shows "Ribosome," "Spliceosome," "Salmonella infection," "Protein processing in endoplasmic reticulum," and several neurodegenerative disease-related pathways at the top. The significance levels (-log(p-val) up to 20, -log(q-val) up to 17) are comparable to the advanced tumor condition, suggesting these alterations emerge early in tumorigenesis.

Biological Interpretation

The Gene Ontology analysis provides critical insights into the functional shifts within Lung Epithelial cells during different stages of lung cancer.

Hallmarks of Cancer Metabolism and Proliferation

Immune Evasion and Inflammation

Unanticipated Overlap with Neurodegenerative Disease Pathways

Clinical or Translational Implications

The GSA results highlight several biologically significant pathways with potential clinical relevance for lung cancer:

  1. Therapeutic Targeting of Core Machinery: The consistent and highly significant up-regulation of "Ribosome," "Spliceosome," "Protein processing in endoplasmic reticulum," and "Proteasome" in both early and advanced tumor cells suggests that targeting these fundamental cellular processes could be effective therapeutic strategies. Inhibitors of proteasome (e.g., bortezomib in multiple myeloma) or spliceosome modulators are already in clinical use or under investigation, offering avenues for repurposing or developing novel agents for lung cancer.
  2. Metabolic Reprogramming as a Target: The enrichment of "Glycolysis / Gluconeogenesis" and "HIF-1 signaling pathway" in diploid cells compared to normal suggests that metabolic vulnerabilities emerge early. Targeting tumor metabolism, such as glycolysis inhibition, could be a strategy for early intervention or prevention in high-risk individuals.
  3. Understanding Immune Context: The distinct pathways enriched in normal cells (e.g., PPAR signaling, antigen processing) underscore the importance of maintaining or restoring anti-tumor immunity. Strategies that enhance antigen presentation or modulate inflammatory pathways could improve immunotherapy responses.
  4. Novel Insights from Neurodegenerative Pathways: The unexpected enrichment of neurodegenerative disease pathways warrants further investigation. This could potentially uncover shared molecular mechanisms that contribute to cancer cell resilience or identify existing drugs for neurodegenerative diseases that could be repurposed for lung cancer, based on shared targets related to protein quality control or cellular stress.

21. Gene Set Enrichment Analysis of Lung Cancer Microenvironment

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents Gene Set Enrichment Analysis (GSEA) results visualized as a dot plot. The analysis compares gene expression profiles across various cell types and conditions found in lung tissue, with all comparisons made against a 'Normal' reference condition, as indicated by the plot title and compute_cfg. The plot highlights biological pathways and gene sets that are significantly upregulated (positive Normalized Enrichment Score, NES, colored red) or downregulated (negative NES, colored blue). The size of each dot reflects the statistical significance of the enrichment, represented by -log(P). This comprehensive view allows for the identification of cell-type-specific and condition-dependent molecular alterations in the context of Lung Cancer.

Visual Summary

The dot plot reveals a complex landscape of pathway dysregulation across different cell types and tumor stages in lung cancer.

Biological Interpretation

The GSEA results provide a detailed molecular blueprint of the lung cancer microenvironment, highlighting crucial biological processes involved in tumor initiation, progression, and interaction with surrounding cells.

Clinical or Translational Implications

The detailed pathway analysis provides valuable insights with potential clinical and translational impact for Lung Cancer.

22. Discussion

The single-cell RNA sequencing analysis of human lung tissue provides an unprecedented resolution into the intricate cellular and molecular changes underlying lung cancer progression. A central finding is the robust identification of malignant lung epithelial cells, primarily originating from alveolar epithelial cells, characterized by extensive aneuploidy and pervasive genomic instability, including recurrent amplifications of key oncogenes such as EGFR and ERBB2 (HER2). These genomic alterations are evident from early disease stages and intensify with progression, confirming their foundational role in oncogenesis.

The tumor microenvironment (TME) undergoes dramatic remodeling, profoundly impacting both immune and stromal compartments. A critical shift towards immunosuppression is observed with the early and sustained enrichment of regulatory T cells (Tregs), which significantly increase in both early and advanced tumor stages. Concurrently, a progressive decrease in Natural Killer (NK) cells and ILC2s highlights a diminished innate anti-tumor response. Macrophages exhibit a striking phenotypic switch, with a significant and progressive increase in pro-tumorigenic M2B macrophage subsets and a decrease in M2A subsets, further contributing to an immunosuppressive and pro-growth milieu. These immune cell dynamics suggest a complex interplay where the host immune system attempts to respond initially but is progressively subverted and exhausted by the advancing tumor.

Cell-cell communication analysis reveals active and dysregulated crosstalk between tumor-origin lung epithelial cells and their TME constituents. Prominent interactions include autocrine/paracrine EGFR signaling among tumor cells and with macrophages and fibroblasts, which drives proliferation and TME conditioning. TGF-beta signaling also emerges as a crucial mediator of tumor-stromal-immune interactions, promoting epithelial-mesenchymal transition (EMT) and fostering an immunosuppressive and fibrotic niche. Notably, tumor cells actively engage in immune evasion by secreting immunosuppressive mediators like Prostaglandin E2 (PGE2) and Retinoic Acid, which modulate T cell and B cell functions. Angiogenesis is robustly supported through VEGF family interactions with endothelial cells, indicating active vascular remodeling.

Gene Ontology and Gene Set Enrichment Analysis further illuminate the molecular underpinnings of these cellular changes. Tumor-origin lung epithelial cells universally exhibit massive upregulation of machinery for protein synthesis, processing, and degradation (ribosomes, spliceosomes, proteasomes), indicative of their heightened proliferative state. Unexpectedly, these cells also show significant enrichment of pathways associated with neurodegenerative diseases, suggesting shared cellular stress responses or protein quality control disruptions critical for tumor cell adaptation. The observation that even diploid lung epithelial cells within the tumor context display significant enrichment of cancer-related pathways suggests that oncogenic transformation can precede overt aneuploidy, or these cells are rapidly conditioned by the evolving tumor microenvironment.

In contrast to commonly observed patterns in some 'hot' tumors, the TME here suggests a landscape trending towards immune exclusion or exhaustion in later stages, characterized by a relative reduction of certain helper T cell subsets from early to advanced disease alongside persistent NK cell suppression. The detailed identification of surfaceome markers across malignant epithelial cells, macrophages, fibroblasts, and CD4+ T cells provides a rich molecular landscape for identifying actionable therapeutic targets, ranging from well-established oncogenic drivers to novel modulators of tumor-microenvironment interactions.

Hypotheses:

  1. Aneuploidy drives early immune evasion: The early onset and progressive increase of aneuploidy in lung epithelial cells directly contribute to the observed shifts in immune cell populations (e.g., Treg enrichment, NK cell suppression) by altering antigen presentation or secreting immunosuppressive factors.
  2. Macrophage polarization is a primary determinant of tumor aggressiveness: The progressive increase in pro-tumorigenic M2B macrophages, coupled with the decline of M2A macrophages, directly correlates with tumor progression and predicts a more aggressive disease course in lung cancer.
  3. Co-targeting EGFR and TME-modulating pathways will overcome resistance: Given the prevalent EGFR amplifications and the active EGFR signaling in tumor cells and stromal cells, combining EGFR inhibitors with agents that target specific TME components (e.g., FAP, CSF1R, TGF-beta, or PGE2 pathways) will overcome adaptive resistance mechanisms and enhance anti-tumor responses.
  4. Neurodegenerative disease pathways represent novel tumor vulnerabilities: The strong and unexpected enrichment of neurodegenerative disease-related pathways in tumor-origin lung epithelial cells indicates shared mechanisms of cellular stress and protein dyshomeostasis that could be exploited therapeutically, potentially through drug repurposing.

Potential therapeutic targets:

  1. EGFR (Epidermal Growth Factor Receptor): Highly and broadly expressed in tumor-origin lung epithelial cells, often amplified (7p13:7q11.23 region), and involved in extensive autocrine/paracrine signaling with TME components (macrophages, fibroblasts), driving proliferation, survival, and TME remodeling. Evidence: Recurrent amplification in tumor-origin cells (Section 4, Image 5). Robustly upregulated in tumor epithelial cells (Section 15, Image 18). Involved in significant cell-cell interactions (EREG_EGFR, AREG_EGFR, HBEGF_EGFR) with tumor cells, macrophages, and fibroblasts (Section 14, Image 17). GSEA confirms activation of MAPK and PI3K-Akt pathways (Section 21, Image 27). The dataset context also notes 'EGFR_mutation' metadata. Validation: Clinical trials with EGFR TKIs (Gefitinib, Erlotinib, Osimertinib) are well-established. Further validation could involve combination therapies targeting EGFR and TME components, e.g., co-targeting EGFR with anti-fibrotic agents in models with high fibroblast activity.
  2. ERBB2 (HER2): Recurrently amplified (17q21.1:17q21.2 region) and overexpressed in a subset of tumor-origin lung epithelial cells, functioning as a co-receptor to enhance pro-survival and pro-proliferative signaling pathways. Evidence: Section 4 shows recurrent amplification (Image 5). Section 15 shows increased expression in a distinct subset of tumor samples (Image 18). GSEA indicates activated PI3K-Akt pathway (Section 21), consistent with ERBB2 signaling. Validation: Clinical trials are ongoing with HER2-targeted therapies like trastuzumab deruxtecan (T-DXd) for HER2-mutant NSCLC. Further validation could involve patient stratification based on ERBB2 amplification status for specific HER2-directed therapies and assessing clinical outcomes.
  3. CSF1R (Colony Stimulating Factor 1 Receptor): Highly expressed on tumor-associated macrophages (TAMs) in early and advanced tumor conditions, playing a critical role in their recruitment, survival, and pro-tumorigenic polarization within the TME. Targeting CSF1R can deplete or reprogram TAMs. Evidence: Section 16 (Image 19) clearly shows CSF1R as a marker with strong enrichment in tumor-associated macrophages. Section 9 and 10 confirm the shift towards M2-like pro-tumorigenic macrophages (M2B, M2C) in tumor conditions, which are often sustained by CSF1R signaling. Validation: Preclinical studies in lung cancer models with CSF1R inhibitors (e.g., Pexidartinib) to assess TAM depletion/reprogramming and anti-tumor efficacy, especially in combination with immune checkpoint inhibitors. Clinical trials are evaluating CSF1R inhibitors in various cancers.
  4. FAP (Fibroblast Activation Protein Alpha): A canonical marker of cancer-associated fibroblasts (CAFs), highly expressed in early tumor-associated fibroblasts but minimally in normal tissues. FAP+ CAFs promote tumor growth, invasion, and metastasis by remodeling the ECM and creating an immunosuppressive environment. Evidence: Section 17 (Image 20) demonstrates strong and prevalent expression of FAP in fibroblasts from early tumor samples, absent in normal. Section 13 and 14 highlight extensive collagen-integrin interactions involving fibroblasts and tumor cells, which FAP contributes to by ECM remodeling. Validation: Preclinical testing of FAP-targeted therapies (e.g., FAP-specific antibodies, CAR T-cells, radioconjugates) in lung cancer models to evaluate CAF depletion, ECM normalization, and enhanced anti-tumor immunity.
  5. TGF-beta Pathway (e.g., TGFB1, TGFbeta_receptor1/2): Active autocrine/paracrine signaling from aneuploid lung epithelial cells to themselves and with TME cells (macrophages, fibroblasts) drives tumor growth, EMT, fibrosis, and immune suppression. TGF-beta pathway components are upregulated in tumor epithelial cells. Evidence: Section 13 and 14 highlight multiple significant TGF-beta-mediated cell-cell interactions: TGFB1_TGFbeta_receptor1 (autocrine in tumor cells), TGFB1_TGFbeta_receptor2 (Mac to tumor cells), TGFB3_TGFBR3 (tumor cells to fibroblasts) (Images 16, 17). Section 19 shows upregulation of TGFB1 in tumor epithelial cells (Image 22). GSEA confirms activation of TGF-beta signaling pathways (Section 21). Validation: Clinical trials with TGF-beta inhibitors (e.g., Galunisertib, Fresolimumab) in lung cancer, often in combination with immune checkpoint inhibitors, to reverse immunosuppression and desmoplasia.
  6. ICOS (CD278): An inducible co-stimulatory molecule highly expressed on tumor-associated CD4+ T cells in both early and advanced stages. It plays a complex role, promoting T cell activation but also influencing immunosuppressive Tregs. Modulating ICOS can enhance anti-tumor T cell responses. Evidence: Section 18 (Image 21) shows ICOS as an upregulated marker in tumor-associated CD4+ T cells. Section 8 details the increase in several helper T cell subsets and Tregs in early tumor stages, which can express ICOS. Validation: Clinical trials with ICOS agonistic antibodies (e.g., JTX-2011) to enhance anti-tumor T cell immunity in lung cancer patients, potentially combined with other immunotherapies. Further studies to dissect its role in effector vs. regulatory T cells in the TME.

Follow-up validation ideas:

  1. Flow cytometry/Immunostaining: Quantify the proportions and activation states of T cell subsets (Tregs, NK cells, helper T cells) and macrophage subsets (M1/M2B/M2A) in a larger cohort of lung cancer patients at different stages and in normal tissues using multi-parameter flow cytometry or spatial transcriptomics combined with immunostaining to validate population shifts and their spatial localization.
  2. CRISPR/shRNA perturbation assays: In vitro and in vivo models of lung cancer (cell lines, patient-derived organoids, PDX models) to knock down or overexpress key genes implicated in specific cell-cell interactions (e.g., PTGES3, ALDH1A1, NRP1, COL1A1) in tumor cells or TME cells, followed by functional assays (proliferation, invasion, immune cell co-culture, angiogenesis).
  3. Drug sensitivity screens: Test candidate therapeutic agents (e.g., EGFR TKIs, CSF1R inhibitors, FAP inhibitors, TGF-beta pathway blockers, EP2/EP4 receptor antagonists) individually and in combinations on lung cancer cell lines and co-culture systems representing the identified TME interactions, followed by dose-response and synergy analyses.
  4. Spatial Transcriptomics/Proteomics: Apply spatial transcriptomics or high-plex spatial proteomics (e.g., CODEX, IMC) to lung cancer tissue sections to confirm the co-localization and spatial interactions of specific cell types (e.g., aneuploid epithelial cells, M2B macrophages, FAP+ fibroblasts) and their surface markers (e.g., CSF1R, FAP, EGFR, MUC4) within the tumor microenvironment.
  5. Genomic/Clinical Cohort Validation: Validate recurrent CNVs (EGFR, ERBB2) and their correlation with patient outcomes, treatment response, and specific histological subtypes using larger independent lung cancer patient cohorts with matched genomic and clinical data (e.g., TCGA, other single-cell or bulk-seq datasets).
  6. Mechanistic studies of neurodegenerative pathways: Conduct biochemical and cell biology experiments to investigate the functional consequences of activated neurodegenerative disease pathways (e.g., protein aggregation, ER stress, mitochondrial dysfunction) in lung cancer cells and how their modulation impacts tumor cell survival, proliferation, and drug sensitivity.

Limitations:

This single-cell analysis provides comprehensive insights into lung cancer; however, certain limitations should be noted. The dataset primarily relies on transcriptional measurements, which may not fully capture post-translational modifications, protein activity, or cell-state dynamics. The inference of CNVs from RNA-seq data, while robust, is an approximation and would benefit from orthogonal genomic validation (e.g., WGS, FISH). While numerous cell-cell interactions are predicted, these are computational inferences based on ligand-receptor expression and require experimental validation to confirm functional significance. The cross-sectional nature of the data does not allow for direct observation of dynamic processes of tumor evolution or immune cell fate decisions in individual patients over time. Finally, the observed patient-specific heterogeneity, while biologically informative, implies that generalized therapeutic strategies may require further stratification based on individual tumor microenvironment characteristics.

23. Query List

  1. Show and save UMAPs for condition, sample, celltype_major, celltype_minor, ploidy_dec, and celltype_subset in two columns.
  2. Show major cell type scores on UMAP and save the result.
  3. Show and save a marker expression dot plot for celltype_subset. SET target_cell = None and var_group_rotation = 45. Leave all other arguments at their default values.
  4. Select Lung Epithelial cells (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 celltype_major, celltype_minor, ploidy_dec, condition, and sample in two columns.
  6. Show and save a population bar plot of celltype_minor.
  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 Lung Epithelial cells (tumor-origin), show their ploidy populations as a bar plot, and save the result.
  12. Show and save cell-cell interaction patterns involving Lung Epithelial cells (tumor-origin), Fibroblast, Macrophage, T cell, B cell, Dendritic cell, and Endothelial cell. Select at most 80 cell-cell interactions per group.
  13. Find cell-cell interactions involving Lung Epithelial cells (tumor-origin), 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 Lung Epithelial cells (tumor-origin), 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 Lung Epithelial cells (tumor-origin) 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 Fibroblast 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. 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.
  19. Among cell cycle pathway genes, select those with statistically significant expression differences between conditions in Lung Epithelial cells (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.
  20. Show and save Gene Ontology (GSA) analysis results for Lung Epithelial cells (tumor-origin) as a bar plot.
  21. Show Gene Set Enrichment Analysis results as a dot plot, and save it. Set the color map to RdBu_r, n_pws_to_show to 120, and target_cells to None.
Top