Prediction of tumor-infiltrating lymphocytes through habitat radiomics and exploration of response mechanisms in neoadjuvant immunochemotherapy-treated lung cancer.
Authors
Affiliations (10)
Affiliations (10)
- Department of Thoracic Surgery, Union Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, 430022, China.
- College of Medicine, Chongqing University, Chongqing, China.
- Radiation Oncology Centre, Chongqing University Cancer Hospital, Chongqing, China.
- Department of Pathology, Union Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China.
- Department of Radiology, Union Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China.
- School of Computer Science and Technology, Huazhong University of Science and Technology, Wuhan, China.
- Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China.
- Department of Thoracic Surgery, Union Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, 430022, China. [email protected].
- Department of Thoracic Surgery, The First Affiliated Hospital of Zhengzhou University, Zhengzhou, Henan, China. [email protected].
- Department of Thoracic Surgery, Union Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, 430022, China. [email protected].
Abstract
Neoadjuvant immunochemotherapy (NAIC) induces tumor microenvironment remodeling in non-small cell lung cancer (NSCLC), presenting challenges for treatment response assessment. This study developed and validated a habitat radiomics approach for non-invasive prediction of tumor-infiltrating lymphocyte (TIL) status to evaluate NAIC response in NSCLC. This retrospective study enrolled 238 NSCLC patients following NAIC for clinical analysis, of which 201 patients met criteria for radiomics analysis. Patients were classified into TIL-positive and TIL-negative groups based on pathological assessment. Post-treatment computed tomography (CT) images were analyzed using K-means clustering to identify tumor habitat sub-regions for radiomic feature extraction. Seven machine learning algorithms were evaluated for TIL status prediction. Model interpretability was assessed through SHapley Additive exPlanations (SHAP) analysis. Single-cell RNA sequencing (scRNA-seq) data were analyzed to compare major pathological response (MPR) and non-MPR tumor microenvironments through cell type annotation, differentiation trajectory analysis, and intercellular communication network analysis. Pre-treatment neutrophil-to-lymphocyte ratio (NLR) showed association with pathological response in multivariable analysis. The radiomics cohort was randomly divided 7:3 into training (n = 140) and test (n = 61) sets. The Random Forest model achieved an area under the receiver operating characteristic curve (AUC) of 0.823 (95% CI: 0.694-0.932) in the test set, and the habitat radiomics model stratified patients into high and low recurrence risk groups. Single-cell analysis identified immunosuppressive features in non-responding tumors, characterized by expansion of SERPINB9 + regulatory T cells (Tregs) that regulated suppressive intercellular communication networks. This study establishes a habitat radiomics model for non-invasive assessment of TIL status following neoadjuvant immunochemotherapy in NSCLC. The model shows reliable predictive performance and prognostic stratification capability, offering potential clinical utility for treatment response evaluation and patient selection.