Introduction: EGFR mutations holds the major targets for drug in lung adenocarcinoma (LUAD). Despite the tremendous study of EGFR mutant (MT) LUAD, the comprehensive interpretation of the heterogeneous character of LUAD harboring EGFR MT remains a key challenge. Here, we investigated the heterogeneity of EGFR MT LUAD and explored the tumor microenvironment (TME) in EGFR MT LUAD.
Method: We performed single-cell RNA sequencing (scRNA-seq) from 135 LUAD patients which consist of normal(n=24), EGFR wild (WT)(n=18), and MT(n=93). Also, we used whole genome sequencing and bulk-RNA sequencing to validate with scRNA-seq results. From 898,648 cells, main cell types were classified. To explore the various characteristics of MT LUAD tumor cells, we used two ways: i) We re-clustered epithelial cells populating the normal, WT, and MT. ii) We re-clustered only MT epithelial cells. In each analysis, we identified the tumor character in the clusters using differential expressed genes analysis, lineage tracing, clinical information, mutation, and trajectory analysis. Also, we extracted each main cell type except epithelial cells, and identified subtypes of main cell types. Finally, we revealed the interaction of cellular components in TME.
Results: In the analysis of epithelial cells, we identified characteristics of specific EGFR MT by comparing of EGFR WT and MT tumors in clusters with similar biological features. The cluster represented by alveolar type 2 (AT2) known as initiation of LUAD was populating normal, WT, and MT. In this cluster, MT- and WT-associated pathway shared but differently significant between MT and WT in the pathway analysis. The cluster represented by proliferative is mostly comprised tumor cells and we found significantly increased the expression of MDK, CD24 in the MT of the cluster. In the analysis of only MT epithelial cells, 2 of clusters were stage-specific cluster: i) The cluster annotated as early stage cluster, ii) The cluster annotated as advanced stage cluster. Trajectory showed that there is a pseudotemporal continuum, following the stage from early stage cluster to advanced stage cluster. Also, based on the lineage tracing, 2 of clusters revealed lineage-specific clusters: i) The cluster annotated as AT2 was enriched from early stage cells, ii) The cluster annotated as basal cell known as origin of lung squamous cell carcinoma(LUSC) was enriched from advanced stage cells. Psedotemporal ordering of these cluster revealed AT2 cluster transdifferentiate into basal cell cluster which implied the possibility of LUAD to LUSC transition by drug resistance. In the interaction of MT and WT TME, the number of signaling received epithelial cells from myeloid cells, endothelial cells, and fibroblasts as sender increased compared with the interaction of normal.
Conclusion: We shed light on the ecosystem of TME according to clinical and biological feature of tumor in EGFR mutant LUAD.
Citation Format: You Won Lee, Eun Ji Lee, Seung Yeon Oh, Kyoung-Ho Pyo, Seong Gu Heo, YoungJoon Park, Su-Jin Choi, Kyumin Lim, Ju-hyeon Lee, Jae Hwan Kim, Jii Bum Lee, Ji Yoon Lee, Sun Min Lim, Chang Gon Kim, Min Hee Hong, Mi Ran Yun, Byoung Chul Cho. Phenotype profiling of tumor microenvironment in EGFR mutant lung adenocarcinoma with multi-omics data. [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 5935.
