Clinical Study of Artificial Intelligence-assisted Diagnosis System in Predicting the Invasive Subtypes of Early-stage Lung Adenocarcinoma Appearing as Pulmonary Nodules

Background: To test the value of Computed tomography (CT) features in predicting the infiltration degree and pathological subtype of ground glass lung adenocarcinoma (≤ 3 cm). Materials and Methods: Data from 412 lung adenocarcinoma patients with mixed ground glass nodules on CT from Jan. 2017 to Dec. 2021 were tested retrospectively. The patients were separated by the infiltrating degree into a minimally invasive adenocarcinoma (MIA) group and an invasive adenocarcinoma (IAC) group. Then the IAC group was subdivided into low-, medium-and high-risk groups by the prognosis differences among subtypes, which were of lepidic, papillary, and micropapillary predominance respectively. Results: Average diameter of nodules, average CT value, solid component ratio, lobe sign, and burr sign were independent risk factors of IAC. The average diameter of nodules ≥ 12.5 mm, solid component ratio ≥ 20.96%, average CT value ≥ -473.07 HU, positive lobe sign and positive burr sign indicated the nodules were more likely IAC. Average CT value, and solid component ratio were independent risk factors for the high-risk pathological type of lung adenocarcinoma. The average CT value ≥ -242.92 HU and solid component ratio ≥ 69.536% indicated nodules were more likely the high-risk pathological type of lung adenocarcinoma. Conclusion: CT imaging features improve the diagnostic efficacy of ground glass nodules, and have certain clinical value.

OA12.05 Noninvasive CT-Based Image Biopsy System (iBiopsy) for Early Stage Lung Adenocarcinoma

Objective: To explore the value of quantitative parameters of artificial intelligence (AI) and computed tomography (CT) signs in identifying pathological subtypes of lung adenocarcinoma appearing as ground-glass nodules (GGNs). Methods: CT images of 224 GGNs from 210 individuals were collected retrospectively and classified into atypical adenomatous hyperplasia (AAH)/adenocarcinoma in situ (AIS), minimally invasive adenocarcinoma (MIA), and invasive adenocarcinoma (IAC) groups. AI was used to identify GGNs and to obtain quantitative parameters, and CT signs were recognized manually. The mixed predictive model based on logistic multivariate regression was built and evaluated. Results: Of the 224 GGNs, 55, 93, and 76 were AAH/AIS, MIA, and IAC, respectively. In terms of AI parameters, from AAH/AIS to MIA, and IAC, there was a gradual increase in two-dimensional mean diameter, three-dimensional mean diameter, mean CT value, maximum CT value, and volume of GGNs (all P<0.0001). Except for the CT signs of the location, and the tumor–lung interface, there were significant differences among the three groups in the density, shape, vacuolar signs, air bronchogram, lobulation, spiculation, pleural indentation, and vascular convergence signs (all P<0.05). The areas under the curve (AUC) of predictive model 1 for identifying the AAH/AIS and MIA and model 2 for identifying MIA and IAC were 0.779 and 0.918, respectively, which were greater than the quantitative parameters independently (all P<0.05). Conclusion: AI parameters are valuable for identifying subtypes of early lung adenocarcinoma and have improved diagnostic efficacy when combined with CT signs.

Objective: To analyze the clinical features and prognostic factors of different histological subtypes of lung adenocarcinoma. Methods: Data from 370 lung adenocarcinoma patients who underwent surgical resection for pathologically supported adenocarcinoma in our hospital between 2000 and 2003 were retro- spectively reviewed. The Kaplan-Meier method was used to estimate patient survival, and Cox’s proportional hazards model was performed for multivariate analysis. Results: The 5-year overall survival rate was 25.26%, and the mean survival time was 3.89 years. In multivariate analysis, histological subtype, incised margin residual, TNM stage, tumor size, and adjuvant chemotherapy were identified as independent survival predictors. The 5-year survival rate in bronchioloalveolar adenocarcinoma (BAC) patients was 41.30%, higher than in patients with other subtypes of lung adenocarcinoma (P=0.002). No significant difference was found in the prognosis among patients with different subtypes of adenocarcinoma without a BAC component. Conclusion: Ade-nocarcinoma with a BAC component is an independent subtype of lung adenocarcinoma. Its prognosis lies between those of BAC and adenocarcinoma without BAC. Histological subtype, incised margin residual, TNM stage, tumor size, and adjuvant chemotherapy are independent survival predictors.

Objective: Solid and micropapillary pattern are highly invasive histologic subtypes in lung adenocarcinoma and are associated with poor prognosis while the biopsy sample is not enough for the accurate histological diagnosis. This study aims to assess the correlation and predictive efficacy between metabolic parameters in (18)F-fluorodeoxy glucose positron emission tomography/computed tomography ((18)F-FDG PET-CT), including the maximum SUV (SUV(max)), metabolic tumor volume (MTV), total lesion glycolysis (TLG) and solid and micropapillary histological subtypes in lung adenocarcinoma. Methods: A total of 145 resected lung adenocarcinomas were included. The clinical data and preoperative (18)F-FDG PET-CT data were retrospectively analyzed. Mann-Whitney U test was used for the comparison of the metabolic parameters between solid and micropapillary subtype group and other subtypes group. Receiver operating characteristic (ROC) curve and areas under curve (AUC) were used for evaluating the prediction efficacy of metabolic parameters for solid or micropapillary patterns. Univariate and multivariate analyses were conducted to determine the prediction factors of the presence of solid or micropapillary subtypes. Results: Median SUV(max) and TLG in solid and papillary predominant subtypes group (15.07 and 34.98, respectively) were significantly higher than those in other subtypes predominant group (6.03 and 10.16, respectively, P<0.05). ROC curve revealed that SUV(max) and TLG had good efficacy for prediction of solid and micropapillary predominant subtypes [AUC=0.811(95% CI: 0.715~0.907) and 0.725(95% CI: 0.610~0.840), P<0.05]. Median SUV(max) and TLG in lung adenocarcinoma with the solid or micropapillary patterns (11.58 and 22.81, respectively) were significantly higher than those in tumors without solid and micropapillary patterns (4.27 and 6.33, respectively, P<0.05). ROC curve revealed that SUV(max) and TLG had good efficacy for predicting the presence of solid or micropapillary patterns [AUC=0.757(95% CI: 0.679~0.834) and 0.681(95% CI: 0.595~0.768), P<0.005]. Multivariate logistic analysis showed that the clinical stage (Stage Ⅲ-Ⅳ), SUV(max) ≥10.27 and TLG≥7.12 were the independent predictive factors of the presence of solid or micropapillary patterns (P<0.05). Conclusions: Preoperative SUV(max) and TLG of lung adenocarcinoma have good prediction efficacy for the presence of solid or micropapillary patterns, especially for the solid and micropapillary predominant subtypes and are independent factors of the presence of solid or micropapillary patterns.

Background: Gene expression-based subtyping in lung adenocarcinoma (AD) and lung squamous cell carcinoma (SQ) classifies AD and SQ tumors into distinct subtypes with variable expression of underlying biology including DNA damage response genes. These subtypes are linked to differences in chemotherapy sensitivity, and may impact response to therapeutics like PARP inhibitors. Methods: Using The Cancer Genome Atlas (TCGA) lung cancer gene expression datasets (AD n=515 and SQ n=501), AD subtypes (Terminal Respiratory Unit (TRU), Proximal Proliferative (PP), and Proximal Inflammatory (PI)) and SQ subtypes (Primitive, Classical, Secretory, and Basal) were defined using gene expression based centroid predictors. Association between AD and SQ expression subtypes and 3 published BRCAness/PARP inhibitor response signatures developed in ovarian and/or breast cancer (Konstantinopoulos et al., PMID 20547991; Daemen et al., PMID 22875744; McGrail et al., PMID 28649435) was examined using linear regression. Association between subtypes and expression of 15 recognized homologous recombination (HR) related genes (ATM, ATR, BRCA1, BRCA2, BRIP1, CDK12, CHEK1, CHEK2, FANCA, FANCI, FANCD2, MRE11A, RAD51C, RAD51L1, PTEN) was also examined using linear regression, and association tests included adjustment for the 3 BRCAness/PARP inhibitor response signatures and proliferation score. Results: AD and SQ subtypes showed strong association with the 3 BRCAness/PARP inhibitor response signatures (F-test p-values 7.7e-05, 5.9e-13, 9.4e-33 in AD and 1.9e-05, 9.0e-13, 2.7e-19 in SQ). AD and SQ subtypes showed strong association with 15 HR genes (max and median F-test p-values were 8.5e-04 and 7.5e-25 in AD, and 7.3e-04 and 1.4e-12 in SQ). The TRU subtype in AD showed low expression relative to the other AD subtypes for a majority of the HR genes, including BRCA1. In SQ, the same was true for the basal and secretory subtypes. Simultaneous adjustment for the 3 BRCAness/PARP inhibitor response signatures, as well as for proliferation, reduced association strength between subtype and HR gene expression in AD and less so in SQ. In AD, association between subtype and gene expression remained significant for 4 HR genes (using Bonferroni correction for 15 tests), including CHECK2, FANCI, BRIP1, and RAD51L1. In SQ, association between subtype and gene expression remained significant for all HR genes except CHEK1 and FANCA, (median and min Bonferroni-adjusted p-value 2.9e-04 and 2.6e-21). Conclusions: Intrinsic biologic subtypes of lung AD and SQ are associated with published BRCAness/PARP inhibitor response signatures and reveal differential expression of several HR-related genes. Evaluation of these subtypes, particularly in SQ, as potential biomarkers of PARP inhibitor sensitivity should be investigated. Citation Format: Gregory Mayhew, Chuck Perou, D Neil Hayes, Myla Lai-Goldman, Hawazin Faruki. Differences in BRCAness/PARP inhibitor response signatures and homologous recombination gene expression across lung adenocarcinoma and squamous cell carcinoma gene expression subtypes [abstract]. In: Proceedings of the AACR-NCI-EORTC International Conference: Molecular Targets and Cancer Therapeutics; 2017 Oct 26-30; Philadelphia, PA. Philadelphia (PA): AACR; Mol Cancer Ther 2018;17(1 Suppl):Abstract nr A037.

To investigate the value of computed tomography (CT)-based radiomics signatures in combination with clinical and CT morphological features to identify epidermal growth factor receptor (EGFR)-mutation subtypes in lung adenocarcinoma (LADC). From February 2012 to October 2019, 608 patients were confirmed with LADC and underwent chest CT scans. Among them, 307 (50.5%) patients had a positive EGFR-mutation and 301 (49.5%) had a negative EGFR-mutation. Of the EGFR-mutant patients, 114 (37.1%) had a 19del -mutation, 155 (50.5%) had a L858R-mutation, and 38 (12.4%) had other rare mutations. Three combined models were generated by incorporating radiomics signatures, clinical, and CT morphological features to predict EGFR-mutation status. Patients were randomly split into training and testing cohorts, 80% and 20%, respectively. Model 1 was used to predict positive and negative EGFR-mutation, model 2 was used to predict 19del and non-19del mutations, and model 3 was used to predict L858R and non-L858R mutations. The receiver operating characteristic curve and the area under the curve (AUC) were used to evaluate their performance. For the three models, model 1 had AUC values of 0.969 and 0.886 in the training and validation cohorts, respectively. Model 2 had AUC values of 0.999 and 0.847 in the training and validation cohorts, respectively. Model 3 had AUC values of 0.984 and 0.806 in the training and validation cohorts, respectively. Combined models that incorporate radiomics signature, clinical, and CT morphological features may serve as an auxiliary tool to predict EGFR-mutation subtypes and contribute to individualized treatment for patients with LADC.

Hyaluronidase, hyaluronan synthase, E-cadherin and TGF-Β profile in lung adenocarcinoma subtypes and squamous cell carcinoma of smokers/nonsmokers

This study aims to develop a CT-based hybrid deep learning network to predict pathological subtypes of early-stage lung adenocarcinoma by integrating residual network (ResNet) with Vision Transformer (ViT).A total of 1411 pathologically confirmed ground-glass nodules (GGNs) retrospectively collected from two centers were used as internal and external validation sets for model development. 3D ResNet and ViT were applied to investigate two deep learning frameworks to classify three subtypes of lung adenocarcinoma namely invasive adenocarcinoma (IAC), minimally invasive adenocarcinoma and adenocarcinoma in situ, respectively. To further improve the model performance, four Res-TransNet based models were proposed by integrating ResNet and ViT with different ensemble learning strategies. Two classification tasks involving predicting IAC from Non-IAC (Task1) and classifying three subtypes (Task2) were designed and conducted in this study.For Task 1, the optimal Res-TransNet model yielded area under the receiver operating characteristic curve (AUC) values of 0.986 and 0.933 on internal and external validation sets, which were significantly higher than that of ResNet and ViT models (p < 0.05). For Task 2, the optimal fusion model generated the accuracy and weighted F1 score of 68.3% and 66.1% on the external validation set.The experimental results demonstrate that Res-TransNet can significantly increase the classification performance compared with the two basic models and have the potential to assist radiologists in precision diagnosis.

Purpose: The aims of this study were to combine CT images with Ki-67 expression to distinguish various subtypes of lung adenocarcinoma and to pre-operatively predict the Ki-67 expression level based on CT radiomic features.Methods: Data from 215 patients with 237 pathologically proven lung adenocarcinoma lesions who underwent CT and immunohistochemical Ki-67 from January 2019 to April 2021 were retrospectively analyzed. The receiver operating curve (ROC) identified the Ki-67 cut-off value for differentiating subtypes of lung adenocarcinoma. A chi-square test or t-test analyzed the differences in the CT images between the negative expression group (n = 132) and the positive expression group (n = 105), and then the risk factors affecting the expression level of Ki-67 were evaluated. Patients were randomly divided into a training dataset (n = 165) and a validation dataset (n = 72) in a ratio of 7:3. A total of 1,316 quantitative radiomic features were extracted from the Analysis Kinetics (A.K.) software. Radiomic feature selection and radiomic classifier were generated through a least absolute shrinkage and selection operator (LASSO) regression and logistic regression analysis model. The predictive capacity of the radiomic classifiers for the Ki-67 levels was investigated through the ROC curves in the training and testing groups.Results: The cut-off value of the Ki-67 to distinguish subtypes of lung adenocarcinoma was 5%. A comparison of clinical data and imaging features between the two groups showed that histopathological subtypes and air bronchograms could be used as risk factors to evaluate the expression of Ki-67 in lung adenocarcinoma (p = 0.005, p = 0.045, respectively). Through radiomic feature selection, eight top-class features constructed the radiomic model to pre-operatively predict the expression of Ki-67, and the area under the ROC curves of the training group and the testing group were 0.871 and 0.8, respectively.Conclusion: Ki-67 expression level with a cut-off value of 5% could be used to differentiate non-invasive lung adenocarcinomas from invasive lung adenocarcinomas. It is feasible and reliable to pre-operatively predict the expression level of Ki-67 in lung adenocarcinomas based on CT radiomic features, as a non-invasive biomarker to predict the degree of malignant invasion of lung adenocarcinoma, and to evaluate the prognosis of the tumor.

The histopathological subtype of lung adenocarcinoma (LUAD) is closely associated with prognosis. Micropapillary or solid predominant LUAD tends to relapse after surgery at an early stage, whereas lepidic pattern shows a favorable outcome. However, the molecular mechanism underlying this phenomenon remains unknown. Here, we recruited 31 lepidic predominant LUADs (LR: low-risk subtype group) and 28 micropapillary or solid predominant LUADs (HR: high-risk subtype group). Tissues of these cases were obtained and label-free quantitative proteomic and bioinformatic analyses were performed. Additionally, prognostic impact of targeted proteins was validated using The Cancer Genome Atlas databases (n = 492) and tissue microarrays composed of early-stage LUADs (n = 228). A total of 192 differentially expressed proteins were identified between tumor tissues of LR and HR and three clusters were identified via hierarchical clustering excluding eight proteins. Cluster 1 (65 proteins) showed a sequential decrease in expression from normal tissues to tumor tissues of LR and then to HR and was predominantly enriched in pathways such as tyrosine metabolism and ECM-receptor interaction, and increased matched mRNA expression of 18 proteins from this cluster predicted favorable prognosis. Cluster 2 (70 proteins) demonstrated a sequential increase in expression from normal tissues to tumor tissues of LR and then to HR and was mainly enriched in pathways such as extracellular organization, DNA replication and cell cycle, and high matched mRNA expression of 25 proteins indicated poor prognosis. Cluster 3 (49 proteins) showed high expression only in LR, with high matched mRNA expression of 20 proteins in this cluster indicating favorable prognosis. Furthermore, high expression of ERO1A and FEN1 at protein level predicted poor prognosis in early-stage LUAD, supporting the mRNA results. In conclusion, we discovered key differentially expressed proteins and pathways between low-risk and high-risk subtypes of early-stage LUAD. Some of these proteins could serve as potential biomarkers in prognostic evaluation.

To develop a model for differentiating the predominant subtype-based prognostic groups of lung adenocarcinoma using CT radiomic features, and to validate its performance in comparison with radiologists' assessments. A total of 993 patients presenting with invasive lung adenocarcinoma between March 2010 and June 2016 were identified. Predominant histologic subtypes were categorized into three groups according to their prognosis (group 0: lepidic; group 1: acinar/papillary; group 2: solid/micropapillary). Seven hundred eighteen radiomic features were extracted from segmented lung cancers on contrast-enhanced CT. A model-development set was formed from the images of 893 patients, while 100 image sets were reserved for testing. A least absolute shrinkage and selection operator method was used for feature selection. Performance of the radiomic model was evaluated using receiver operating characteristic curve analysis, and accuracy on the test set was compared with that of three radiologists with varying experiences (6, 7, and 19years in chest CT). Our model differentiated the three groups with areas under the curve (AUCs) of 0.892 and 0.895 on the development and test sets, respectively. In pairwise discrimination, the AUC was highest for group 0 vs. 2 (0.984). The accuracy of the model on the test set was higher than the averaged accuracy of the three radiologists without statistical significance (73.0% vs. 61.7%, p = 0.059). For group 2, the model achieved higher PPV than the observers (85.7% vs. 35.0-48.4%). Predominant subtype-based prognostic groups of lung adenocarcinoma were classified by a CT-based radiomic model with comparable performance to radiologists. • A CT-based radiomic model differentiated three prognosis-based subtype groups of lung adenocarcinoma with areas under the curve (AUCs) of 0.892 and 0.895 on development and test sets, respectively. • The CT-based radiomic model showed near perfect discrimination between group 0 and group 2 (AUCs, 0.984-1.000). • The accuracy of the CT-based radiomic model was comparable to the averaged accuracy of the three radiologists with 6, 7, and 19years of clinical experience in chest CT (73.0% vs. 61.7%, p = 0.059), achieving a higher positive predictive value for group 2 than the observers (85.7% vs. 35.0-48.4%).

Lung adenocarcinoma (LUAD) is one of the most widespread and fatal types of lung cancer. Oxidative stress, resulting from an imbalance in the production and accumulation of reactive oxygen species (ROS), is considered a promising therapeutic target for cancer treatment. Currently, immune checkpoint blockade (ICB) therapy is being explored as a potentially effective treatment for early-stage LUAD. In this research, we aim to identify distinct subtypes of LUAD patients by investigating genes associated with oxidative stress and immunotherapy. Additionally, we aim to propose subtype-specific therapeutic strategies. We conducted a thorough search of the Gene Expression Omnibus (GEO) datasets. From this search, we pinpointed datasets that contained both expression data and survival information. We selected genes associated with oxidative stress and immunotherapy using keyword searches on GeneCards. We then combined expression data of LUAD samples from both The Cancer Genome Atlas (TCGA) and 11 GEO datasets, forming a unified dataset. This dataset was subsequently divided into two subsets, Dataset_Training and Dataset_Testing, using a random bifurcation method, with each subset containing 50% of the data. We applied consensus clustering (CC) analysis to identify distinct LUAD subtypes within the Dataset_Training. Molecular variances associated with oxidative stress levels, the tumor microenvironment (TME), and immune checkpoint genes (ICGs) were then investigated among these subtypes. Employing feature selection combined with machine learning techniques, we constructed models that achieved the highest accuracy levels. We validated the identified subtypes and models from Dataset_Training using Dataset_Testing. A hub gene with the highest importance values in the machine learning model was identified. We then utilized virtual screening to discover potential compounds targeting this hub gene. In the unified dataset, we integrated 2,154 LUAD samples from TCGA-LUAD and 11 GEO datasets. We specifically selected 1,311 genes associated with immune and oxidative stress processes. The expression data of these genes were then employed for subtype identification through CC analysis. Within Dataset_Training, two distinct subtypes emerged, each marked by different levels of immune and oxidative stress pathway values. Consequently, we named these as the OX+ and IM+ subtypes. Notably, the OX+ subtype showed increased oxidative stress levels, correlating with a worse prognosis than the IM+ subtype. Conversely, the IM+ subtype demonstrated enhanced levels of immune pathways, immune cells, and ICGs compared to the OX+ subtype. We reconfirmed these findings in Dataset_Testing. Through gene selection, we identified an optimal combination of 12 genes for predicting LUAD subtypes: ACP1, AURKA, BIRC5, CYC1, GSTP1, HSPD1, HSPE1, MDH2, MRPL13, NDUFS1, SNRPD1, and SORD. Out of the four machine learning models we tested, the support vector machine (SVM) stood out, achieving the highest area under the curve (AUC) of 0.86 and an accuracy of 0.78 on Dataset_Testing. We focused on HSPE1, which was designated as the hub gene due to its paramount importance in the SVM model, and computed the docking structures for four compounds: ZINC3978005 (Dihydroergotamine), ZINC52955754 (Ergotamine), ZINC150588351 (Elbasvir), and ZINC242548690 (Digoxin). Our study identified two subtypes of LUAD patients based on oxidative stress and immunotherapy-related genes. Our findings provided subtype-specific therapeutic strategies.

The aim of this study was to investigate the differentiation of computed tomography (CT)-based entropy parameters between minimally invasive adenocarcinoma (MIA) and invasive adenocarcinoma (IAC) lesions appearing as pulmonary subsolid nodules (SSNs). This study was approved by the institutional review board in our hospital. From July 2015 to November 2018, 186 consecutive patients with solitary peripheral pulmonary SSNs that were pathologically confirmed as pulmonary adenocarcinomas (74 MIA and 112 IAC lesions) were included and subdivided into the training data set and the validation data set. Chest CT scans without contrast enhancement were performed in all patients preoperatively. The subjective CT features of the SSNs were reviewed and compared between the MIA and IAC groups. Each SSN was semisegmented with our in-house software, and entropy-related parameters were quantitatively extracted using another in-house software developed in the MATLAB platform. Logistic regression analysis and receiver operating characteristic analysis were performed to evaluate the diagnostic performances. Three diagnostic models including subjective model, entropy model, and combined model were built and analyzed using area under the curve (AUC) analysis. There were 119 nonsolid nodules and 67 part-solid nodules. Significant differences were found in the subjective CT features among nodule type, lesion size, lobulated shape, and irregular margin between the MIA and IAC groups. Multivariate analysis revealed that part-solid type and lobulated shape were significant independent factors for IAC (P < 0.0001 and P < 0.0001, respectively). Three entropy parameters including Entropy-0.8, Entropy-2.0-32, and Entropy-2.0-64 were identified as independent risk factors for the differentiation of MIA and IAC lesions. The median entropy model value of the MIA group was 0.266 (range, 0.174-0.590), which was significantly lower than the IAC group with value 0.815 (range, 0.623-0.901) (P < 0.0001). Multivariate analysis revealed that the combined model had an excellent diagnostic performance with sensitivity of 88.2%, specificity of 73.0%, and accuracy of 82.1%. The AUC value of the combined model was significantly higher (AUC, 0.869) than that of the subjective model (AUC, 0.809) or the entropy model alone (AUC, 0.836) (P < 0.0001). The CT-based entropy parameters could help assess the aggressiveness of pulmonary adenocarcinoma via quantitative analysis of intratumoral heterogeneity. The MIA can be differentiated from IAC accurately by using entropy-related parameters in peripheral pulmonary SSNs.

Lung adenocarcinoma is the most prevalent type of lung cancer, with invasive lung adenocarcinoma being the most common subtype. Screening and early treatment of high-risk individuals have improved survival; however, significant differences in prognosis still exist among patients at the same stage, especially in the early stages. Invasive lung adenocarcinoma has different histological morphologies and biological characteristics that can distinguish its prognosis. Notably, several studies have found that the pathological subtypes of invasive lung adenocarcinoma are closely associated with clinical treatment. This review summarised the distribution of various pathological subtypes of invasive lung adenocarcinoma in the population and their relationship with sex, smoking, imaging features, and other histological characteristics. We comprehensively analysed the genetic characteristics and biomarkers of the different pathological subtypes of invasive lung adenocarcinoma. Understanding the interaction between the pathological subtypes of invasive lung adenocarcinoma and the tumour microenvironment helps to reveal new therapeutic targets for lung adenocarcinoma. We also extensively reviewed the prognosis of various pathological subtypes and their effects on selecting surgical methods and adjuvant therapy and explored future treatment strategies.