{"doi":"10.3389/fimmu.2025.1637625","title":"Single-cell transcriptomic profiling reveals the heterogeneity of epithelial cells in lung adenocarcinoma lymph node metastasis and develops a prognostic signature","abstract":"<jats:sec><jats:title>Background</jats:title><jats:p>Lymph node metastasis markedly worsens prognosis in lung adenocarcinoma (LUAD); however, the evolutionary dynamics and regulatory mechanisms underlying the heterogeneity of malignant epithelial cells during this process remain poorly understood and warrant comprehensive investigation.</jats:p></jats:sec><jats:sec><jats:title>Methods</jats:title><jats:p>We performed a comprehensive single-cell transcriptomic analysis of epithelial cells from 18 samples comprising normal lung tissue and lymph node metastases. Malignant epithelial cells were identified via inferred copy number variation (CNV) profiles. Key malignant subpopulations were further characterized through trajectory inference, cell–cell communication mapping, gene set variation analysis (GSVA), and reconstruction of transcription factor regulatory networks. To assess clinical relevance, we developed and validated a prognostic model—termed the EAS score—based on the transcriptional signatures of malignant epithelial subsets, using integrated data from multiple TCGA and GEO cohorts. The functional role of the hub gene SELENBP1 was experimentally validated through quantitative PCR (qPCR), Western blotting, immunohistochemistry (IHC), Transwell migration assays, colony formation assays, flow cytometry, ROS quantification, and subcutaneous tumorigenesis assays <jats:italic>in vivo</jats:italic>.</jats:p></jats:sec><jats:sec><jats:title>Results</jats:title><jats:p>Single-cell transcriptomic analysis identified four distinct malignant epithelial subtypes (Clusters 0–3), each characterized by unique patterns of CNV. Leveraging these defined cellular subpopulations, we constructed a highly accurate model for prognostication in LUAD, enabling reliable classification of patients based on clinical outcomes. Through detailed comparisons between groups with divergent prognostic risks, the study revealed notable differences across the tumor microenvironment (TME), including alterations in pathway activity, gene enrichment distributions, mutation profiles, and anticipated responses to immune checkpoint blockade. In addition, functional validation experiments confirmed that SELENBP1 plays a tumor-suppressive role, further supporting its relevance as a potential intervention target in LUAD.</jats:p></jats:sec><jats:sec><jats:title>Conclusion</jats:title><jats:p>This research provides insights into the evolutionary complexity and heterogeneity of malignant epithelial populations in lymph node metastatic sites of LUAD. It also presents a scoring system based on prognostic indicators, which serves as a reliable tool for forecasting patient survival outcomes. Moreover, the discovery of SELENBP1 as a candidate tumor suppressor emphasizes its importance in guiding both clinical risk categorization and the design of personalized treatment strategies for individuals classified as high-risk LUAD cases.</jats:p></jats:sec>","journal":"Frontiers in Immunology","year":2025,"id":645034,"datarank":0.10397207708399181,"base_score":0.6931471805599453,"endowment":0.6931471805599453,"self_citation_contribution":0.10397207708399181,"citation_network_contribution":0.0,"self_endowment_contribution":0.10397207708399181,"citer_contribution":0.0,"corpus_percentile":null,"corpus_rank":null,"citation_count":1,"citer_count":0,"citers_with_citation_signal":0,"citers_with_endowment":0,"datacite_reuse_total":0,"is_dataset":false,"is_dataset_confidence":null,"is_data_producer":false,"deposit_databanks":null,"is_oa":false,"file_count":0,"downloads":0,"has_version_chain":false,"published_date":null,"fair_score":null,"fair_percentile":null,"algorithm_id":"datarank_citation_only_1hop_v6","ranking_scope":"data_only","authors":[{"id":6714,"name":"Han Zhang","orcid":"0000-0002-9991-9409","position":1,"is_corresponding":false},{"id":516511,"name":"Ying Shi","orcid":"0000-0002-8509-3503","position":2,"is_corresponding":false},{"id":1679350,"name":"Mengli Xue","orcid":null,"position":3,"is_corresponding":false},{"id":1679351,"name":"Jingxian Wang","orcid":null,"position":4,"is_corresponding":false},{"id":1250999,"name":"Yun Ding","orcid":"0000-0002-0354-2388","position":5,"is_corresponding":false},{"id":525840,"name":"Lin Tan","orcid":"0000-0002-4921-1293","position":6,"is_corresponding":false},{"id":495604,"name":"Hui Yuan","orcid":"0000-0002-9122-3407","position":7,"is_corresponding":false},{"id":801787,"name":"Xin Li","orcid":"0000-0002-8779-125X","position":8,"is_corresponding":false},{"id":1149868,"name":"Daqiang Sun","orcid":"0009-0000-6602-357X","position":9,"is_corresponding":false},{"id":1679349,"name":"Qiuqiao Mu","orcid":null,"position":0,"is_corresponding":false}],"reference_count":0,"raw_metadata":{"has_enrichment":true,"resolved":true,"title":"Single-cell transcriptomic profiling reveals the heterogeneity of epithelial cells in lung adenocarcinoma lymph node metastasis and develops a prognostic signature","abstract":"<jats:sec><jats:title>Background</jats:title><jats:p>Lymph node metastasis markedly worsens prognosis in lung adenocarcinoma (LUAD); however, the evolutionary dynamics and regulatory mechanisms underlying the heterogeneity of malignant epithelial cells during this process remain poorly understood and warrant comprehensive investigation.</jats:p></jats:sec><jats:sec><jats:title>Methods</jats:title><jats:p>We performed a comprehensive single-cell transcriptomic analysis of epithelial cells from 18 samples comprising normal lung tissue and lymph node metastases. Malignant epithelial cells were identified via inferred copy number variation (CNV) profiles. Key malignant subpopulations were further characterized through trajectory inference, cell–cell communication mapping, gene set variation analysis (GSVA), and reconstruction of transcription factor regulatory networks. To assess clinical relevance, we developed and validated a prognostic model—termed the EAS score—based on the transcriptional signatures of malignant epithelial subsets, using integrated data from multiple TCGA and GEO cohorts. The functional role of the hub gene SELENBP1 was experimentally validated through quantitative PCR (qPCR), Western blotting, immunohistochemistry (IHC), Transwell migration assays, colony formation assays, flow cytometry, ROS quantification, and subcutaneous tumorigenesis assays <jats:italic>in vivo</jats:italic>.</jats:p></jats:sec><jats:sec><jats:title>Results</jats:title><jats:p>Single-cell transcriptomic analysis identified four distinct malignant epithelial subtypes (Clusters 0–3), each characterized by unique patterns of CNV. Leveraging these defined cellular subpopulations, we constructed a highly accurate model for prognostication in LUAD, enabling reliable classification of patients based on clinical outcomes. Through detailed comparisons between groups with divergent prognostic risks, the study revealed notable differences across the tumor microenvironment (TME), including alterations in pathway activity, gene enrichment distributions, mutation profiles, and anticipated responses to immune checkpoint blockade. In addition, functional validation experiments confirmed that SELENBP1 plays a tumor-suppressive role, further supporting its relevance as a potential intervention target in LUAD.</jats:p></jats:sec><jats:sec><jats:title>Conclusion</jats:title><jats:p>This research provides insights into the evolutionary complexity and heterogeneity of malignant epithelial populations in lymph node metastatic sites of LUAD. It also presents a scoring system based on prognostic indicators, which serves as a reliable tool for forecasting patient survival outcomes. Moreover, the discovery of SELENBP1 as a candidate tumor suppressor emphasizes its importance in guiding both clinical risk categorization and the design of personalized treatment strategies for individuals classified as high-risk LUAD cases.</jats:p></jats:sec>","is_dataset_classified":null,"base_score":0.6931471805599453,"endowment":0.6931471805599453,"datacite_reuse_total":0,"file_count":0,"downloads":0,"views":0,"has_version_chain":false,"is_dataset":false,"is_oa":false,"pmid":"40787454","pmcid":"PMC12331613","openalex_id":"https://openalex.org/W4412924584","authors":[],"funders":[],"total_grants":0,"fwci":0.4217,"citation_percentile":0.60978391,"influential_citations":0,"citation_trend":[{"year":2026,"count":1}],"oa_status":"gold","license":"cc-by","oa_locations":[{"url":"https://doi.org/10.3389/fimmu.2025.1637625","host_type":"journal"},{"url":"https://doi.org/10.3389/fimmu.2025.1637625","host_type":"publisher"},{"url":"https://www.frontiersin.org/articles/10.3389/fimmu.2025.1637625/full","host_type":"publisher"},{"url":"https://pubmed.ncbi.nlm.nih.gov/40787454","host_type":"repository"},{"url":"https://doaj.org/article/5364a85ee6a54cbdb66f5772b9db269f","host_type":"repository"},{"url":"https://www.ncbi.nlm.nih.gov/pmc/articles/12331613","host_type":"repository"},{"url":"https://europepmc.org/articles/PMC12331613","host_type":"Europe_PMC"},{"url":"https://europepmc.org/articles/PMC12331613?pdf=render","host_type":"Europe_PMC"}],"fields_of_study":["Single-cell and spatial transcriptomics","Ferroptosis and cancer prognosis","Cancer Immunotherapy and Biomarkers","Humans","Lung Neoplasms","Single-Cell Analysis","Prognosis","Adenocarcinoma of Lung","Gene Expression Profiling","Lymphatic Metastasis","Epithelial Cells","Animals","Transcriptome","Gene Expression Regulation, Neoplastic","Female","DNA Copy Number Variations","Mice","Male","Biomarkers, Tumor","Middle Aged"],"mesh_terms":["Adenocarcinoma of Lung","Animals","Epithelial Cells","Female","Humans","Lung Neoplasms","Lymphatic Metastasis","Male","Middle Aged","Prognosis","Biomarkers, Tumor","Gene Expression Regulation, Neoplastic","Gene Expression Profiling","Mice","DNA Copy Number Variations","Single-Cell Analysis","Transcriptome"],"keywords":["Biology","Transcriptome","Carcinogenesis","Cancer research","Metastasis","Lymph node","Adenocarcinoma","Tumor microenvironment","Gene expression profiling","Gene signature","Flow cytometry","Gene expression","Cancer","Gene","Immunology","Genetics","Immunotherapy","Epithelial cell","lymph node metastasis","Selenbp1","Luad"],"sdg_mappings":[],"linked_datasets":[],"clinical_trials":[],"software_tools":[],"database_accessions":[],"source":"live","citation_network_status":"fetched"},"created_at":"2026-08-09T02:39:59.551911Z","pmid":null,"pmcid":null,"fwci":null,"citation_percentile":null,"influential_citations":0,"oa_status":null,"license":null,"views":0,"total_file_size_bytes":0,"version_count":0,"fair_f":null,"fair_a":null,"fair_i":null,"fair_r":null,"fair_zscore":null,"fair_rationale":null,"fair_model":null,"fair_agent_version":null,"fair_fulltext_source":null,"fair_has_llm":null,"fair_computed_at":null,"clinical_trials":[],"software_tools":[],"db_accessions":[],"linked_datasets":[],"topics":[]}