The impact of spread through air spaces on intrapulmonary metastasis in lung cancer: a retrospective study of multiple primary lung adenocarcinomas with confirmed common origin based on histological classification alone and combined histological-molecular classification
Original Article

The impact of spread through air spaces on intrapulmonary metastasis in lung cancer: a retrospective study of multiple primary lung adenocarcinomas with confirmed common origin based on histological classification alone and combined histological-molecular classification

Chenxi Yan1 ORCID logo, Shaojie Hu1, Yitao Tian1, Xiaoxue Zhang2, Wei Sun1

1Department of Thoracic Surgery, Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China; 2Department of Cardiac and Great Vascular Surgery, Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China

Contributions: (I) Conception and design: S Hu, W Sun; (II) Administrative support: W Sun; (III) Provision of study materials or patients: C Yan, Y Tian; (IV) Collection and assembly of data: C Yan, X Zhang; (V) Data analysis and interpretation: S Hu, Y Tian; (VI) Manuscript writing: All authors; (VII) Final approval of manuscript: All authors.

Correspondence to: Wei Sun, MD, PhD. Department of Thoracic Surgery, Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology, 1095 Jiefang Avenue, Wuhan 430030, China. Email: petersw@tjh.tjmu.edu.cn.

Background: Distinguishing multiple primary lung adenocarcinomas (MPLAs) from intrapulmonary metastases (IPMs) is essential for accurate staging and treatment planning. However, existing histologic and molecular criteria remain limited in diagnostic accuracy. Spread through air spaces (STAS) has been recognized as a marker of tumor aggressiveness, but its role in differentiating MPLA from IPM remains unclear. This study aimed to assess whether STAS can serve as a supportive marker to differentiate IPM from MPLA under both histologic and molecular frameworks.

Methods: We retrospectively enrolled 170 patients who were preoperatively diagnosed with MPLA and underwent surgical resection of multiple lung lesions. Clinical, histopathologic, and molecular data—including broad-panel next-generation sequencing (NGS) results [e.g., epidermal growth factor receptor (EGFR), Kirsten rat sarcoma viral oncogene homolog (KRAS), tumor protein p53 (TP53), among others]—were collected. Tumors were classified as MPLA or IPM using histologic criteria alone and a combined histologic-molecular algorithm. The presence and quantity of STAS were evaluated and analyzed as potential predictors of IPM.

Results: STAS was more frequent in the IPM group across both classification methods and was an independent predictor of IPM. Driver gene mutations (e.g., EGFR, KRAS) differed significantly between MPLA and IPM. The 83.5% concordance between histologic and combined classifications suggests potential bias in histology alone. In resource-limited settings, STAS may offer a cost-effective adjunct to support IPM diagnosis when molecular data are unavailable.

Conclusions: STAS is associated with IPM and may serve as a supportive diagnostic marker when distinguishing IPM from MPLA. In resource-limited settings where molecular testing is not always available, STAS may offer a practical, cost-effective histologic adjunct to guide clinical decision-making.

Keywords: Spread through air spaces (STAS); multiple primary lung adenocarcinomas (MPLAs); intrapulmonary metastasis (IPM); histological classification; molecular testing


Submitted Mar 15, 2025. Accepted for publication Jun 05, 2025. Published online Aug 26, 2025.

doi: 10.21037/tlcr-2025-304


Highlight box

Key findings

• This study demonstrated that the presence of spread through air spaces (STAS) was significantly more frequent in patients with intrapulmonary metastasis (IPM) than in those with multiple primary lung adenocarcinomas (MPLAs), and that STAS is an independent predictor of IPM across both histologic and combined histologic-molecular classification systems.

What is known and what is new?

• STAS is a recognized adverse prognostic factor in lung adenocarcinoma and is associated with recurrence. However, its role in distinguishing IPM from MPLA has not been well established.

• This study provides new evidence supporting the diagnostic value of STAS as a supportive histopathological marker to differentiate IPM from MPLA, particularly in the postoperative setting.

What is the implication, and what should change now?

• Accurate differentiation between MPLA and IPM is crucial for proper staging and treatment planning, especially to avoid overtreatment or undertreatment. In clinical settings where molecular testing is not always accessible, such as in resource-limited regions, STAS may serve as a cost-effective adjunct to guide further diagnostic workup or raise suspicion of IPM. STAS should be incorporated into multidisciplinary evaluation strategies to enhance diagnostic precision and optimize patient management.


Introduction

Lung cancer, particularly lung adenocarcinoma in non-small cell lung cancer (NSCLC), is a leading cause of cancer-related deaths worldwide (1), despite advances in treatment methods, the pathological diagnosis of multiple lung cancers, especially the differentiation between multiple primary lung adenocarcinomas (MPLAs) and intrapulmonary metastases (IPMs), remains a significant challenge. These two entities require distinctly different clinical management approaches. IPM is generally managed with systemic palliative chemotherapy, whereas MPLA are typically treated with curative surgical resection when feasible. As a result, accurate classification directly influences treatment decisions and long-term outcomes. Recent study has shown that overall survival and recurrence-free intervals are significantly longer in MPLA compared to IPM (2).

Spread through air spaces (STAS), a unique mode of metastasis, has recently gained attention among researchers. This metastatic pattern refers to tumor cells entering the alveolar spaces through alveoli and colonizing other regions of the lung (3,4), potentially having a close association with IPM. Although STAS differs from traditional IPM (such as hematogenous or lymphatic spread), whether it directly leads to IPM, particularly in multiple lung adenocarcinomas, remains unclear. The presence of STAS may indicate local extension of the primary tumor or more aggressive biological behavior, which directly affects tumor staging, surgical decision-making, and the need for systemic therapy. Misclassifying IPM as MPLA may lead to inappropriate curative surgery, while misdiagnosing MPLA as IPM could deny patients potentially curative treatment. Therefore, understanding the diagnostic relevance of STAS is critical for optimizing treatment planning and prognostication (5-7).

Accurately distinguishing between MPLAs and IPMs is crucial for clinical management. Traditional histological methods, such as the martini and melamed criteria (8), while helpful in judgment, have shown significant limitations in the context of rapid advances in molecular biology (9). They rely heavily on morphological similarities and lack molecular assessment, making them prone to misclassification when tumors share histologic features but differ at the genomic level. Most lung adenocarcinomas display mixed histological subtypes, and lymph node status—a key component of the martini and melamed criteria—has limited discriminatory value. Furthermore, study has shown that up to 44% of tumors classified as MPLA by histology alone may actually be clonally related based on next-generation sequencing (NGS) (10). This underscores the need to incorporate molecular data and complementary markers such as STAS to improve diagnostic precision. If STAS indeed promotes IPM, it may influence tumor biological behavior and alter clinical management strategies. Furthermore, the role of gene mutations and molecular pathways in lung cancer metastasis has gradually been recognized. Mutations in genes such as Kirsten rat sarcoma viral oncogene homolog (KRAS) and epidermal growth factor receptor (EGFR) may play important roles in IPM (11,12), but their role in STAS remains unclear (13,14). This is particularly problematic when tumors share common driver mutations—such as EGFR or KRAS—since these mutations may independently arise due to field cancerization or environmental exposure, especially in high-prevalence populations (15). Therefore, relying on shared mutations alone may lead to misclassification. Given these limitations, there is a growing need for additional markers to support the differentiation of IPM from MPLA. STAS, a histologic pattern associated with tumor aggressiveness and recurrence, has recently emerged as a potential supportive feature.

This study aimed to investigate the potential association between STAS and IPM in patients with multiple lung adenocarcinomas. By integrating histological and molecular features, we sought to explore whether STAS could serve as a complementary pathological marker to assist in distinguishing IPM from MPLA, thereby providing additional insight into the existing diagnostic framework. We present this article in accordance with the STROBE reporting checklist (available at https://tlcr.amegroups.com/article/view/10.21037/tlcr-2025-304/rc).


Methods

Patient and tumor selection

A total of 170 patients who were preoperatively diagnosed with MPLA were retrospectively enrolled. Clinical data (including gender, smoking history, tumor timing, location, among others), pathological specimens, and molecular testing results [including at least EGFR, KRAS, tumor protein p53 (TP53), among others] were collected for analysis.

Inclusion criteria were: (I) preoperative diagnosis of MPLA and postoperative pathological confirmation of lung adenocarcinoma; (II) presence of two or more lesions; (III) complete clinical, radiological, and pathological data; and (IV) molecular testing covering at least key genes such as EGFR, KRAS, TP53, ERBB2, ALK, among others.

Exclusion criteria included: (I) non-adenocarcinoma pathological types; (II) patients who had received radiotherapy, chemotherapy, or targeted therapy before surgery; (III) incomplete clinical data; and (IV) absence of any detectable mutations.

This study was approved by the hospital ethics committee, and all data were anonymized (see Figure 1).

Figure 1 The flowchart of patient selection and analysis. IPM, intrapulmonary metastasis; MPLA, multiple primary lung adenocarcinoma.

This study was conducted in accordance with the Declaration of Helsinki and its subsequent amendments. The study was approved by the Medical Ethics Committee of Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology (No. TJ-IRB202502150). As this was a retrospective study, the requirement for informed consent was waived by the ethics committee.

Power and sample size calculation

A priori power analysis was conducted using G*Power 3.1 software to estimate the required sample size for detecting a difference in STAS positivity between the IPM and MPLA groups. Assuming a two-sided test with α =0.05, power =0.80, and the observed proportions (IPM: 47.9%, MPLA: 23.8%), the minimum required sample size was calculated as 153 (43 in IPM, and 110 in MPLA). Our actual sample size of 170 exceeded this threshold, confirming that the study had sufficient power (actual power =0.80).

Clinical and imaging data collection

Demographic information (gender, age, smoking history, etc.), tumor timing (patients had synchronous tumors or metachronous, defined by a second tumor occurring at least 2 years after the first one), tumor location (same lobe or different lobes in the same lung), and lymph node status were collected. All clinical data were obtained from the hospital’s electronic medical record system and independently verified by two researchers before being entered into the database.

Histological evaluation

All surgical specimens were reviewed by two experienced lung cancer pathologists (blinded). Pathological evaluation included tumor size, major histological subtypes (acinar, papillary, lepidic, etc.), pleural invasion, vascular invasion, and the presence of STAS. Major histological subtypes were determined based on the 2015 World Health Organization (WHO) classification criteria for adenocarcinoma. The proportion of each histological component was recorded in 5% increments, with the top three components documented. STAS was assessed by examining all tumor borders and selecting one to three STAS-rich regions. Semi-quantitative evaluation was performed under a microscope at 20× objective and 10× eyepiece magnification. Tumor cells detached into alveolar spaces, forming scattered or clustered distributions, were observed without artificial manipulation. The STAS was classified into three categories: high STAS (≥5 clusters of STAS), low STAS (1–4 clusters of the micropapillary or solid nest-predominant STAS), and no STAS (16). And the only clusters showing a continuous trajectory from the main tumor to the alveolar spaces were considered true STAS, while clusters that appeared randomly scattered, lacked directional orientation, or were distant without continuity were excluded (4,17). Additional histological features recorded included: (I) histological structural patterns; (II) cytological features (cell size, cuboidal or columnar cells, nuclear heterogeneity, intranuclear inclusions, nucleolar size, mitotic rate, cytoplasmic content, and mucin cell count); (III) tumor stroma (cellular components, immune cell count, and type); and (IV) thyroid transcription factor 1 (TTF-1) expression. Based on these assessments, patients were preliminarily classified into MPLA and IPM groups.

Molecular analysis

For each patient’s tumor specimen, paraffin-embedded tissue blocks with high tumor cell content (≥30%) were prioritized for DNA extraction using standardized protocols. Molecular testing was performed using a NGS platform, covering key driver genes such as EGFR, KRAS, TP53, among others, as well as hotspot mutations in other relevant genes. Although the testing panels and mutation sites varied among patients, all cases included at least EGFR, KRAS, TP53, among others. Mutations were defined as having a variant allele frequency ≥2% and coverage depth ≥300×. All variant results were validated through bioinformatics analysis and manual review to ensure data accuracy.

Patient grouping and classification criteria

Two grouping methods were employed in this study:

Histological classification: the distinction between IPM and MPLA was based on histological criteria from previous studies (18,19), including: (I) histological structural patterns; (II) cytological features (cell size, cuboidal or columnar cells, nuclear heterogeneity, intranuclear inclusions, nucleolar size, mitotic rate, cytoplasmic content, and mucin cell count); (III) tumor stroma (cellular components, immune cell count, and type); and (IV) TTF-1 expression (SP24 clone, Novocastra, Newcastle, UK). Based on these criteria, patients were initially classified into MPLA and IPM groups (detailed in Figure 2A).

Figure 2 Histological and histomolecular classification strategies for differentiating MPLA and IPM. (A) A histology-based classification approach, aligned with the WHO guidelines, is used to distinguish multiple lung adenocarcinomas as either MPLA or IPM. This method involves a detailed assessment of the tumor’s architectural patterns, cytological and stromal characteristics, as well as expression of TTF-1. (B) A combined histological and molecular classification model has been proposed to improve diagnostic accuracy. In this approach, rare mutations are defined as those excluding classical EGFR mutations (such as exon 19 deletions and p.L858R) and KRAS p.G12X mutations. AIS, adenocarcinoma in situ; BRAF, B-Raf proto-oncogene; EGFR, epidermal growth factor receptor; ERBB2, erb-b2 receptor tyrosine kinase 2; IPM, intrapulmonary metastasis; KRAS, Kirsten rat sarcoma viral oncogene homolog; MET, MNNG HOS transforming gene; MIA, minimally invasive adenocarcinoma; MPLA, multiple primary lung adenocarcinoma; NRAS, neuroblastoma RAS viral oncogene homolog; TP53, tumor protein p53; TTF-1, thyroid transcription factor-1; WHO, World Health Organization.

Combined histological-molecular classification: a classification algorithm integrating clinical, histological, and molecular features from previous study was used (20). First, if one lesion was adenocarcinoma in situ (AIS) or minimally invasive adenocarcinoma (MIA), the tumor pair was classified as MPLAs. Second, for metachronous lesions, if the interval was ≥5 years, they were classified as MPLAs. Subsequently, molecular analysis results were used for further classification, except for cases where tumors shared a single common hotspot mutation—specifically EGFR exon 19 deletion, EGFR L858R, or KRAS G12X—and no additional shared alterations were detected, classification was based on histologic criteria. The final classification results were compared with the Martini and Melamed classification. Cases with shared hotspot mutations were classified as IPM based on molecular information, while others were classified as MPLA (detailed in Figure 2B).

Statistical analysis

All statistical analyses were performed using R software and RStudio on Windows. Continuous variables were expressed as mean ± standard deviation (SD) for normally distributed data, and as median with interquartile range (IQR) for non-normally distributed data. Normality was assessed using the Shapiro-Wilk test. Comparisons between groups were performed using independent t-tests for normally distributed variables or Mann-Whitney U tests for non-normally distributed variables. Categorical variables were expressed as percentages, and comparisons between groups were performed using chi-square tests or Fisher’s exact tests. Univariate and multivariate logistic regression analyses were used to assess the impact of STAS, pleural invasion, vascular invasion, and lymph node status on IPM, with results expressed as odds ratios (ORs) and 95% confidence intervals (95% CIs). A two-sided P<0.05 was considered statistically significant.


Results

Patient and histological characteristics

A total of 170 patients were included in this study, of whom 63 (37.1%) were male and 30 (17.6%) were smokers. Among the patients, 18 (10.6%) had metachronous tumors with an interval of more than 5 years, 65 (38.2%) had metachronous tumors with an interval of less than 5 years, and 87 (51.2%) had synchronous tumors. Lymph node-positive disease was observed in 13 patients (7.6%). In 59 patients (34.7%), both tumors were located in the same lobe, while in 97 patients (57.1%), the tumors were located in different lobes of the same lung. The three most common major histological subtypes were acinar (125 cases, 73.5%), papillary (80 cases, 47.1%), and lepidic (47 cases, 27.6%) (see Table 1). STAS was observed in 52 patients (30.6%), pleural invasion was detected in 23 patients (13.5%), and vascular invasion was present in 6 patients (3.5%) (see Table 1).

Table 1

Clinical, pathological, and molecular characteristics of 170 MPLA and IPM patients based on histological classification and combined histological-molecular classification

Variables Total
(n=170)
Histological classification Combined histological-molecular classification
MPLA (n=94) IPM (n=76) P MPLA (n=122) IPM (n=48) P
Sex 0.17 0.19
   Female 107 (62.9) 64 (68.1) 43 (56.6) 81 (66.4) 26 (54.2)
   Male 63 (37.1) 30 (31.9) 33 (43.4) 41 (33.6) 22 (45.8)
Age at first resection (years) 59.00
[52.25, 66.00]
59.00
[52.00, 64.00]
60.00
[53.00, 68.00]
0.17 58.50
[52.00, 65.00]
60.00
[54.50, 67.25]
0.17
Smoking history 0.10 0.36
   No 140 (82.4) 82 (87.2) 58 (76.3) 103 (84.4) 37 (77.1)
   Yes 30 (17.6) 12 (12.8) 18 (23.7) 19 (15.6) 11 (22.9)
Patients with different tumor chronology 0.56 0.005
   Synchronous tumors 87 (51.2) 46 (48.9) 41 (53.9) 55 (45.1) 32 (66.7)
   Metachronous tumors (<5 years) 65 (38.2) 36 (38.3) 29 (38.2) 49 (40.2) 16 (33.3)
   Metachronous tumors (>5 years) 18 (10.6) 12 (12.8) 6 (7.9) 18 (14.8) 0 (0.0)
Distribution of tumors 0.73 0.65
   Ipsilateral (same lobe) 59 (34.7) 31 (33.0) 28 (36.8) 40 (32.8) 19 (39.6)
   Ipsilateral (other lobe) 97 (57.1) 54 (57.4) 43 (56.6) 71 (58.2) 26 (54.2)
   Contralateral 14 (8.2) 9 (9.6) 5 (6.6) 11 (9.0) 3 (6.2)
Type of resection 0.18 0.66
   Nonanatomic resection (wedge resection) 9 (5.3) 4 (4.3) 5 (6.6) 6 (4.9) 3 (6.2)
   Segmentectomy 62 (36.5) 40 (42.6) 22 (28.9) 47 (38.5) 15 (31.2)
   Lobectomy 99 (58.2) 50 (53.2) 49 (64.5) 69 (56.6) 30 (62.5)
Side 0.45 0.84
   Left 62 (36.5) 37 (39.4) 25 (32.9) 44 (36.1) 18 (37.5)
   Right 94 (55.3) 48 (51.1) 46 (60.5) 67 (54.9) 27 (56.2)
   Bilateral 14 (8.2) 9 (9.6) 5 (6.6) 11 (9.0) 3 (6.2)
Tumor size (cm) 1.15±0.53 0.99±0.44 1.35±0.56 <0.001 1.06±0.49 1.36±0.55 0.001
Nodal status 0.03 0.24
   N0 or Nx 157 (92.4) 91 (96.8) 66 (86.8) 115 (94.3) 42 (87.5)
   N-positive 13 (7.6) 3 (3.2) 10 (13.2) 7 (5.7) 6 (12.5)
STAS 0.006 0.004
   No (both tumors) 118 (69.4) 74 (78.7) 44 (57.9) 93 (76.2) 25 (52.1)
   Yes (≥1 tumor) 52 (30.6) 20 (21.3) 32 (42.1) 29 (23.8) 23 (47.9)
STAS quantity 0.01 0.009
   No 118 (69.4) 74 (78.7) 44 (57.9) 93 (76.2) 25 (52.1)
   Low 34 (20.0) 13 (13.8) 21 (27.6) 19 (15.6) 15 (31.2)
   High 18 (10.6) 7 (7.4) 11 (14.5) 10 (8.2) 8 (16.7)
Pleural invasion 0.32 0.32
   No (both tumors) 147 (86.5) 84 (89.4) 63 (82.9) 108 (88.5) 39 (81.2)
   Yes (≥1 tumor) 23 (13.5) 10 (10.6) 13 (17.1) 14 (11.5) 9 (18.8)
Vascular involvement 0.49 0.10
   No (both tumors) 164 (96.5) 92 (97.9) 72 (94.7) 120 (98.4) 44 (91.7)
   Yes (≥1 tumor) 6 (3.5) 2 (2.1) 4 (5.3) 2 (1.6) 4 (8.3)
EGFR 0.23 0.02
   WT (both tumors) 73 (42.9) 36 (38.3) 37 (48.7) 45 (36.9) 28 (58.3)
   Mut (≥1 tumor) 97 (57.1) 58 (61.7) 39 (51.3) 77 (63.1) 20 (41.7)
KRAS 0.37 0.09
   WT (both tumors) 109 (64.1) 57 (60.6) 52 (68.4) 73 (59.8) 36 (75.0)
   Mut (≥1 tumor) 61 (35.9) 37 (39.4) 24 (31.6) 49 (40.2) 12 (25.0)
Histological subtype
   AIS <0.001 0.005
    No (both tumors) 149 (87.6) 73 (77.7) 76 (100.0) 101 (82.8) 48 (100.0)
    Yes (≥1 tumor) 21 (12.4) 21 (22.3) 0 (0.0) 21 (17.2) 0 (0.0)
   MIA <0.001 <0.001
    No (both tumors) 124 (72.9) 48 (51.1) 76 (100.0) 76 (62.3) 48 (100.0)
    Yes (≥1 tumor) 46 (27.1) 46 (48.9) 0 (0.0) 46 (37.7) 0 (0.0)
   Acinar 0.049 0.10
    No (both tumors) 45 (26.5) 31 (33.0) 14 (18.4) 37 (30.3) 8 (16.7)
    Yes (≥1 tumor) 125 (73.5) 63 (67.0) 62 (81.6) 85 (69.7) 40 (83.3)
   Lepidic 0.02 0.01
    No (both tumors) 123 (72.4) 61 (64.9) 62 (81.6) 81 (66.4) 42 (87.5)
    Yes (≥1 tumor) 47 (27.6) 33 (35.1) 14 (18.4) 41 (33.6) 6 (12.5)
   Tubular 0.50 >0.99
    No (both tumors) 165 (97.1) 90 (95.7) 75 (98.7) 118 (96.7) 47 (97.9)
    Yes (≥1 tumor) 5 (2.9) 4 (4.3) 1 (1.3) 4 (3.3) 1 (2.1)
   Papillary 0.59 0.98
    No (both tumors) 90 (52.9) 52 (55.3) 38 (50.0) 64 (52.5) 26 (54.2)
    Yes (≥1 tumor) 80 (47.1) 42 (44.7) 38 (50.0) 58 (47.5) 22 (45.8)
   Solid >0.99 0.32
    No (both tumors) 162 (95.3) 90 (95.7) 72 (94.7) 118 (96.7) 44 (91.7)
    Yes (≥1 tumor) 8 (4.7) 4 (4.3) 4 (5.3) 4 (3.3) 4 (8.3)
   Mucinous 0.47 0.68
    No (both tumors) 166 (97.6) 93 (98.9) 73 (96.1) 120 (98.4) 46 (95.8)
    Yes (≥1 tumor) 4 (2.4) 1 (1.1) 3 (3.9) 2 (1.6) 2 (4.2)
   Cribriform 0.77 0.48
    No (both tumors) 166 (97.6) 91 (96.8) 75 (98.7) 118 (96.7) 48 (100.0)
    Yes (≥1 tumor) 4 (2.4) 3 (3.2) 1 (1.3) 4 (3.3) 0 (0.0)
   Micropapillary >0.99 0.47
    No (both tumors) 161 (94.7) 89 (94.7) 72 (94.7) 117 (95.9) 44 (91.7)
    Yes (≥1 tumor) 9 (5.3) 5 (5.3) 4 (5.3) 5 (4.1) 4 (8.3)

Data are presented as n (%), median [IQR] or mean ± SD. AIS, adenocarcinoma in situ; EGFR, epidermal growth factor receptor; IQR, interquartile range; IPM, intrapulmonary metastasis; KRAS, Kirsten rat sarcoma viral oncogene homolog; MIA, minimally invasive adenocarcinoma; MPLA, multiple primary lung adenocarcinoma; Mut, mutated; SD, standard deviation; STAS, spread through air spaces; WT, wild type.

Molecular characteristics of tumors

Among the 340 tumors analyzed for gene mutations, 72 tumors (42.3%) harbored at least two mutations. Commonly mutated genes included TP53 (n=123), KRAS (n=93, of which 68 were KRAS p.G12C), EGFR (n=171, including 62 with exon 19 deletions, 79 with EGFR p.L858R, 13 with exon 20 insertions, and 17 with other mutations), STK11 (n=5), and BRAF (n=5) (see Figure 3).

Figure 3 Global distribution of gene mutations in the 170 studied lung adenocarcinomas. A total of 340 lung lesions (2 per patient) were analyzed, corresponding to 170 patients. EGFR, epidermal growth factor receptor; ERBB2, erb-b2 receptor tyrosine kinase 2; KRAS, Kirsten rat sarcoma viral oncogene homolog; TP53, tumor protein p53.

Classification of IPM and MPLA patients

A total of 170 patients were included in the study, classified into two groups based on histological classification: MPLA (n=94) and IPM (n=76). When integrating histological and molecular classifications, 122 patients were categorized as MPLA and 48 as IPM. No significant differences were observed in sex distribution, median age at first resection, smoking history, tumor laterality, or type of resection between the groups (all P>0.05).

There was no significant difference in the proportion of synchronous and metachronous tumors between the MPLA and IPM groups under histological classification (P=0.56). However, under the combined classification, synchronous tumors were more common in the IPM group (66.7%) compared to the MPLA group (45.1%) (P=0.005). The distribution of tumors (same lobe, ipsilateral different lobe, or contralateral) did not differ significantly between the groups (all P>0.05).

The median tumor size was significantly larger in the IPM group compared to the MPLA group in both classification systems (histological: 1.35±0.56 vs. 0.99±0.44 cm, P<0.001; combined: 1.36±0.55 vs. 1.06±0.49 cm, P=0.001). Nodal involvement (N-positive) was significantly more frequent in the IPM group than in the MPLA group under histological classification (13.2% vs. 3.2%, P=0.03), but this difference was not statistically significant in the combined classification (12.5% vs. 5.7%, P=0.24).

The presence of tumor STAS in at least one tumor was significantly higher in the IPM group than in the MPLA group under both classification systems (histological: 42.1% vs. 21.3%, P=0.006; combined: 47.9% vs. 23.8%, P=0.004). Additionally, STAS quantity was significantly higher in the IPM group, with a greater proportion of cases exhibiting low and high STAS (all P<0.05). No significant differences were observed in pleural invasion or vascular involvement between groups (all P>0.05).

EGFR mutations were more frequently detected in at least one tumor in the MPLA group than in the IPM group under the combined classification (63.1% vs. 41.7%, P=0.02), while no significant difference was observed under histological classification (P=0.23). KRAS mutation rates did not significantly differ between the groups in either classification system (all P>0.05).

Significant differences in histological subtypes were observed between groups. The MPLA group had a higher proportion of AIS and MIA compared to the IPM group (both P<0.001). The IPM group under histological classification exhibited a higher frequency of acinar subtype (P), while the MPLA group had a higher prevalence of lepidic growth patterns (all P<0.05). No significant differences were found in the distribution of papillary, tubular, solid, cribriform, micro papillary or mucinous subtypes (all P>0.05) (see Table 1).

Logistic regression analysis of STAS and IPM

We first classified patients based on histological methods and performed logistic regression analysis to evaluate the impact of STAS, pleural invasion, vascular involvement, and nodal status on IPM. In univariate analysis, STAS (OR =2.691; 95% CI: 1.385–5.335; P=0.004) and nodal status (OR =4.596; 95% CI: 1.346–21.089; P=0.02) were significantly associated with IPM. In contrast, pleural invasion (OR =1.733; 95% CI: 0.717–4.305; P=0.22) and vascular involvement (OR =2.556; 95% CI: 0.485–18.805; P=0.29) did not show statistical significance. In multivariate analysis, STAS remained an independent predictor of IPM (OR =2.267; 95% CI: 1.107–4.714; P=0.03), while the association of nodal status lost statistical significance after adjusting for other factors (OR =3.073; 95% CI: 0.750–16.730; P=0.14).

After further grouping patients using combined histological and molecular methods, we repeated the logistic regression analysis. The results showed that, in univariate analysis, STAS (OR =2.950; 95% CI: 1.462–5.998; P=0.003) remained significantly associated with IPM, while the OR for vascular involvement increased (OR =5.455; 95% CI: 1.028–40.347; P=0.055) but did not reach statistical significance. In multivariate analysis, STAS remained an independent predictor (OR =2.576; 95% CI: 1.200–5.526; P=0.01), while pleural invasion (OR =1.498; 95% CI: 0.551–3.864; P=0.41), vascular involvement (OR =2.753; 95% CI: 0.422–23.963; P=0.30), and nodal status (OR =0.952; 95% CI: 0.218–3.662; P=0.94) did not show statistical significance (see Table 2 and Figure 4).

Table 2

Results of logistic regression analysis based on histological classification and combined histological-molecular classification

Variable Histological classification Combined histological-molecular classification
Univariable Multivariable Univariable Multivariable
OR (95% CI) P OR (95% CI) P OR (95% CI) P OR (95% CI) P
STAS 2.691
(1.385–5.335)
0.004 2.267
(1.107–4.714)
0.03 2.950
(1.462–5.998)
0.003 2.576
(1.200–5.526)
0.01
Pleural invasion 1.733
(0.717–4.305)
0.22 1.451
(0.569–3.739)
0.43 1.780
(0.692–4.395)
0.22 1.498
(0.551–3.864)
0.41
Vascular involvement 2.556
(0.485–18.805)
0.29 0.701
(0.089–6.217)
0.73 5.455
(1.028–40.347)
0.055 2.753
(0.422–23.963)
0.30
Nodal status 4.596
(1.346–21.089)
0.02 3.073
(0.750–16.730)
0.14 2.347
(0.719–7.461)
0.14 0.952
(0.218–3.662)
0.94

CI, confidence interval; OR, odds ratio; STAS, spread through air spaces.

Figure 4 Forest plot of multivariate logistic regression analysis stratified by classification methods. (A) Histology-based classification. (B) Integrated histology and molecular classification. Odds ratios and 95% confidence intervals are displayed for each variable. STAS, spread through air spaces.

In summary, STAS was identified as an independent predictor of IPM in both histological and combined histological-molecular classification systems. However, STAS was also observed in a subset of patients classified as having MPLA, indicating that it is not exclusive to IPM. These findings suggest that while STAS can serve as a supportive indicator of intrapulmonary spread, it should be interpreted alongside histological and molecular features to improve diagnostic accuracy. In contrast, pleural invasion, vascular involvement, and nodal status did not demonstrate significant associations with IPM in multivariate analysis.


Discussion

This study systematically analyzed the clinicopathological features, molecular profiles, and STAS expression patterns in patients with multiple lung adenocarcinomas, aiming to explore the potential diagnostic value of STAS in distinguishing MPLA from IPM. While we observed that histological classification alone may result in ambiguity, especially in cases with overlapping morphology, the integration of molecular data helped clarify clonal relationships in several discordant cases. Specifically, 28 patients exhibited inconsistency between histological and combined classifications, mainly due to divergent driver gene mutations such as KRAS and EGFR, reinforcing prior evidence that histological similarity does not guarantee clonal identity.

Our findings demonstrated that STAS was significantly more prevalent in IPM and independently associated with intrapulmonary spread, suggesting that STAS may serve as a supportive marker for differentiating IPM from MPLA. Given the potential for histologic misclassification and the limitations of molecular testing in some settings (10), STAS could provide additional diagnostic value.

STAS, as an independent poor prognostic factor, was significantly more prevalent in the IPM group compared to the MPLA group (47.9% vs. 23.8%). The higher frequency of STAS in IPM suggests a potential role in intrapulmonary dissemination and distal seeding. Previous studies have shown that STAS-positive lung cancer patients are more likely to experience postoperative recurrence (21-24), supporting its association with aggressive tumor behavior. Our findings reinforce this observation and indicate that STAS may serve as a supportive pathological marker in differentiating IPM from MPLA. However, the presence of STAS in approximately one-quarter of patients classified as MPLA highlights that STAS is not specific to IPM. Therefore, while STAS can aid in diagnostic evaluation, it should be interpreted in the context of molecular alterations, histological features, and clinical presentation, rather than used as a standalone criterion. Additionally, the correlation between STAS and multifocal lesions may reflect underlying mechanisms of local invasion and spread, warranting further investigation into its molecular basis in lung adenocarcinoma progression. Furthermore, in clinical practice—particularly in developing countries—access to comprehensive molecular testing may be limited due to cost constraints and resource availability. In such settings, STAS assessment via routine pathological examination offers a practical, cost-effective adjunct to support the diagnostic differentiation between MPLA and IPM. Although it cannot replace molecular analysis, STAS may help to indicate a propensity for intrapulmonary spread when molecular data are unavailable, thus informing treatment decisions in resource-limited environments.

We observed significant differences in EGFR mutation patterns between MPLA and IPM. EGFR mutations were more frequent in the MPLA group (63.1% vs. 41.7%, P=0.02), which is consistent with previous studies reporting that EGFR-mutant lung cancers are more likely to present as multiple primary lesions (25,26). In contrast, KRAS mutations did not differ significantly between the two groups. Notably, prior research has suggested that KRAS mutations may confer greater invasiveness and metastatic potential (27,28), indicating that their role in distinguishing MPLA from IPM remains complex and may extend beyond simple prevalence differences. Although BRAF and STK11 mutations were rare in the overall cohort, they warrant further investigation to explore their roles in multifocal lung lesions.

However, this study has several limitations. First, due to its retrospective design, some patients lacked complete molecular testing data. Although current diagnostic practice primarily treats STAS as a binary variable (present or absent), previous study has suggested that the extent or density of STAS—including the number of detached tumor cell clusters or their distance from the main tumor—may be associated with recurrence risk and clinical outcomes (16). Therefore, future studies with larger sample sizes and objective grading systems are needed to further determine whether the “degree” of STAS can enhance its diagnostic or prognostic value. Additionally, the single-center nature and limited sample size may affect the generalizability of the results. Furthermore, our study did not incorporate emerging technologies such as spatial transcriptomics or single-cell sequencing, which could provide higher-resolution insights into clonal evolution and further elucidate the molecular mechanisms underlying MPLA and IPM.

In conclusion, our study demonstrates that STAS is significantly more prevalent in IPM than in MPLA and serves as an independent predictor of IPM across both histologic and molecular classification systems. While not specific to IPM, the presence of STAS may reflect aggressive tumor behavior and can support diagnostic interpretation when integrated with other clinicopathologic and genomic data. In clinical settings where access to broad molecular testing is limited—such as in many developing countries—STAS may serve as a practical and cost-effective histologic adjunct. However, given that our analysis was based on surgically resected specimens, the utility of STAS in preoperative biopsy samples remains limited and warrants further investigation.


Conclusions

STAS was significantly associated with IPM and identified as an independent predictor of IPM under both histologic and combined histologic–molecular classification systems. Although STAS is not exclusive to IPM, its presence may provide valuable supportive information for differentiating IPM from MPLA. Given the limitations of histologic and molecular criteria and the diagnostic challenges in resource-limited settings, STAS may serve as a practical and cost-effective histopathologic adjunct to guide postoperative clinical decision-making. Future studies are warranted to further validate the diagnostic value of STAS and explore its integration into multidisciplinary evaluation frameworks.


Acknowledgments

None.


Footnote

Reporting Checklist: The authors have completed the STROBE reporting checklist. Available at https://tlcr.amegroups.com/article/view/10.21037/tlcr-2025-304/rc

Data Sharing Statement: Available at https://tlcr.amegroups.com/article/view/10.21037/tlcr-2025-304/dss

Peer Review File: Available at https://tlcr.amegroups.com/article/view/10.21037/tlcr-2025-304/prf

Funding: This study was supported by the Chen Xiaoping Science and Technology Development Foundation (No. CXPJJH1180008 to W.S.).

Conflicts of Interest: All authors have completed the ICMJE uniform disclosure form (available at https://tlcr.amegroups.com/article/view/10.21037/tlcr-2025-304/coif). W.S. reports that the study was supported by the Chen Xiaoping Science and Technology Development Foundation (No. CXPJJH1180008). The other authors have no conflicts of interest to declare.

Ethical Statement: The authors are accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. This study was conducted in accordance with the Declaration of Helsinki and its subsequent amendments. The study was approved by the Medical Ethics Committee of Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology (No. TJ-IRB202502150). As this was a retrospective study, informed consent was waived by the ethics committee.

Open Access Statement: This is an Open Access article distributed in accordance with the Creative Commons Attribution-NonCommercial-NoDerivs 4.0 International License (CC BY-NC-ND 4.0), which permits the non-commercial replication and distribution of the article with the strict proviso that no changes or edits are made and the original work is properly cited (including links to both the formal publication through the relevant DOI and the license). See: https://creativecommons.org/licenses/by-nc-nd/4.0/.


References

  1. Thai AA, Solomon BJ, Sequist LV, et al. Lung cancer. Lancet 2021;398:535-54. [Crossref] [PubMed]
  2. Yang CY, Yeh YC, Wang LC, et al. Genomic Profiling With Large-Scale Next-Generation Sequencing Panels Distinguishes Separate Primary Lung Adenocarcinomas From Intrapulmonary Metastases. Mod Pathol 2023;36:100047. [Crossref] [PubMed]
  3. Giraud P, Antoine M, Larrouy A, et al. Evaluation of microscopic tumor extension in non-small-cell lung cancer for three-dimensional conformal radiotherapy planning. Int J Radiat Oncol Biol Phys 2000;48:1015-24. [Crossref] [PubMed]
  4. Kadota K, Nitadori JI, Sima CS, et al. Tumor Spread through Air Spaces is an Important Pattern of Invasion and Impacts the Frequency and Location of Recurrences after Limited Resection for Small Stage I Lung Adenocarcinomas. J Thorac Oncol 2015;10:806-14. [Crossref] [PubMed]
  5. Kozower BD, Larner JM, Detterbeck FC, et al. Special treatment issues in non-small cell lung cancer: Diagnosis and management of lung cancer, 3rd ed: American College of Chest Physicians evidence-based clinical practice guidelines. Chest 2013;143:e369S-99S.
  6. Felip E, Altorki N, Zhou C, et al. Adjuvant atezolizumab after adjuvant chemotherapy in resected stage IB-IIIA non-small-cell lung cancer (IMpower010): a randomised, multicentre, open-label, phase 3 trial. Lancet 2021;398:1344-57. [Crossref] [PubMed]
  7. Wu YL, Tsuboi M, He J, et al. Osimertinib in Resected EGFR-Mutated Non-Small-Cell Lung Cancer. N Engl J Med 2020;383:1711-23. [Crossref] [PubMed]
  8. Martini N, Melamed MR. Multiple primary lung cancers. J Thorac Cardiovasc Surg 1975;70:606-12.
  9. Girard N, Deshpande C, Lau C, et al. Comprehensive histologic assessment helps to differentiate multiple lung primary nonsmall cell carcinomas from metastases. Am J Surg Pathol 2009;33:1752-64. [Crossref] [PubMed]
  10. Chang JC, Alex D, Bott M, et al. Comprehensive Next-Generation Sequencing Unambiguously Distinguishes Separate Primary Lung Carcinomas From Intrapulmonary Metastases: Comparison with Standard Histopathologic Approach. Clin Cancer Res 2019;25:7113-25. [Crossref] [PubMed]
  11. Wang Z, Yuan X, Jiang G, et al. Towards the molecular era of discriminating multiple lung cancers. EBioMedicine 2023;90:104508. [Crossref] [PubMed]
  12. Iwata T, Sugio K, Uramoto H, et al. Detection of EGFR and K-ras mutations for diagnosis of multiple lung adenocarcinomas. Front Biosci (Landmark Ed) 2011;16:2961-9. [Crossref] [PubMed]
  13. Lee JS, Kim EK, Kim M, et al. Genetic and clinicopathologic characteristics of lung adenocarcinoma with tumor spread through air spaces. Lung Cancer 2018;123:121-6. [Crossref] [PubMed]
  14. Laville D, Désage AL, Fournel P, et al. Spread Through Air Spaces in Stage I to III Resected Lung Adenocarcinomas: Should the Presence of Spread Through Air Spaces Lead to an Upstaging? Am J Surg Pathol 2024;48:596-604. [Crossref] [PubMed]
  15. Chou TY, Dacic S, Wistuba I, et al. Differentiating Separate Primary Lung Adenocarcinomas From Intrapulmonary Metastases With Emphasis on Pathological and Molecular Considerations: Recommendations From the International Association for the Study of Lung Cancer Pathology Committee. J Thorac Oncol 2025;20:311-30. [Crossref] [PubMed]
  16. Uruga H, Fujii T, Fujimori S, et al. Semiquantitative Assessment of Tumor Spread through Air Spaces (STAS) in Early-Stage Lung Adenocarcinomas. J Thorac Oncol 2017;12:1046-51. [Crossref] [PubMed]
  17. Aly RG, Rekhtman N, Li X, et al. Spread Through Air Spaces (STAS) Is Prognostic in Atypical Carcinoid, Large Cell Neuroendocrine Carcinoma, and Small Cell Carcinoma of the Lung. J Thorac Oncol 2019;14:1583-93. [Crossref] [PubMed]
  18. Detterbeck FC, Franklin WA, Nicholson AG, et al. The IASLC Lung Cancer Staging Project: Background Data and Proposed Criteria to Distinguish Separate Primary Lung Cancers from Metastatic Foci in Patients with Two Lung Tumors in the Forthcoming Eighth Edition of the TNM Classification for Lung Cancer. J Thorac Oncol 2016;11:651-65.
  19. Nicholson AG, Chansky K, Crowley J, et al. The International Association for the Study of Lung Cancer Lung Cancer Staging Project: Proposals for the Revision of the Clinical and Pathologic Staging of Small Cell Lung Cancer in the Forthcoming Eighth Edition of the TNM Classification for Lung Cancer. J Thorac Oncol 2016;11:300-11.
  20. Mansuet-Lupo A, Barritault M, Alifano M, et al. Proposal for a Combined Histomolecular Algorithm to Distinguish Multiple Primary Adenocarcinomas from Intrapulmonary Metastasis in Patients with Multiple Lung Tumors. J Thorac Oncol 2019;14:844-56. [Crossref] [PubMed]
  21. Dagher S, Sulaiman A, Bayle-Bleuez S, et al. Spread Through Air Spaces (STAS) Is an Independent Prognostic Factor in Resected Lung Squamous Cell Carcinoma. Cancers (Basel) 2022;14:2281. [Crossref] [PubMed]
  22. Isgir BB, Kocaman G, Kahya Y, et al. Combination of grade and spread through air spaces (STAS) predicts recurrence in early stage lung adenocarcinoma: a retrospective cohort study. Updates Surg 2025;77:201-8. [Crossref] [PubMed]
  23. Liu A, Hou F, Qin Y, et al. Predictive value of a prognostic model based on pathologic features in lung invasive adenocarcinoma. Lung Cancer 2019;131:14-22. [Crossref] [PubMed]
  24. Kang YK, Song YS, Cho S, et al. Prognostic stratification model for patients with stage I non-small cell lung cancer adenocarcinoma treated with surgical resection without adjuvant therapies using metabolic features measured on F-18 FDG PET and postoperative pathologic factors. Lung Cancer 2018;119:1-6. [Crossref] [PubMed]
  25. Fan J, Dai X, Wang Z, et al. Concomitant EGFR Mutation and EML4-ALK Rearrangement in Lung Adenocarcinoma Is More Frequent in Multifocal Lesions. Clin Lung Cancer 2019;20:e517-30. [Crossref] [PubMed]
  26. Wu C, Zhao C, Yang Y, et al. High Discrepancy of Driver Mutations in Patients with NSCLC and Synchronous Multiple Lung Ground-Glass Nodules. J Thorac Oncol 2015;10:778-83. [Crossref] [PubMed]
  27. Jones GD, Caso R, Tan KS, et al. KRAS (G12C) Mutation Is Associated with Increased Risk of Recurrence in Surgically Resected Lung Adenocarcinoma. Clin Cancer Res 2021;27:2604-12. [Crossref] [PubMed]
  28. Calles A, Sholl LM, Rodig SJ, et al. Immunohistochemical Loss of LKB1 Is a Biomarker for More Aggressive Biology in KRAS-Mutant Lung Adenocarcinoma. Clin Cancer Res 2015;21:2851-60. [Crossref] [PubMed]
Cite this article as: Yan C, Hu S, Tian Y, Zhang X, Sun W. The impact of spread through air spaces on intrapulmonary metastasis in lung cancer: a retrospective study of multiple primary lung adenocarcinomas with confirmed common origin based on histological classification alone and combined histological-molecular classification. Transl Lung Cancer Res 2025;14(8):2969-2982. doi: 10.21037/tlcr-2025-304

Download Citation