Introduction

Anatomical pulmonary resection in combination with mediastinal lymph node dissection (MLND) or sampling remains the standard surgical approach even for early-stage nonā€“small-cell lung cancer (NSCLC) patients1,2. However, the nodal metastatic rate among early NSCLC patients is not high. The reported nodal metastatic rate is 16% to 29% in early-stage NSCLC patients with a peripheral tumor 3ā€‰cm or less in diameter, 15% to 20% in patients with a peripheral tumor 2ā€‰cm or less3,4,5, and as low as 6.6% in clinical T1aN0M0 lung adenocarcinoma patients6. Considering the surgical injury, MLND is unnecessary and harmful to patients without mediastinal node involvement2,7,8. Correctly predicting pathologic node-negative (pN0) disease in patients with NSCLC before operation may avoid unnecessary MLND, although it is an extremely difficult problem in clinical practice. Because the relative low diagnostic sensitivity, these tumors with no nodal involvement on computed tomography (CT) or even positron emission tomography (PET)/CT image, which are staged as clinical N0 (cN0) disease, may also have pathological regional lymph node metastases9,10. Consequently, MLND is still adopted as a gold standard in almost NSCLC surgical treatments.

Recent advancements in the quality of high-resolution CT (HRCT) and the widespread practice of low-dose CT (LDCT) for lung cancer screening have resulted in a greater increase in the early detection of smaller-sized NSCLC, especially lung adenocarcinoma appearing as subsolid nodule11,12. As a result, limited resection like wedge resection or segmenectomy is more favored and the necessity of standard lobectomy begin to be questioned for small or subsolid types of adenocarcinoma2. Meanwhile, the degree to which MLND is necessary also becomes a challenging issue because the biological behavior of small lung adenocarcinomas showing as subsolid nodules on CT is less aggressive and usually pN0 diseases. These patients are good candidates for the omission of MLND. Therefore, it is very important to establish reliable prediction criteria for pN0 disease preoperatively especially using quantitative radiographic imaging to avoid performing MLND.

In this study, we attempt to analyze the value of radiographic and histopathological features of primary tumors for predicting pN0 disease in cT1N0M0 lung adenocarcinoma to establish an optimal surgical strategy for avoiding unnecessary MLND in these patients. We analyzed both HRCT and PET/CT indexes, and we also proposed three-dimensional computer quantified CT number analysis to calculate the solid component (Sc) proportion and ground glass opacity component (GGOc) proportion.

Results

Patient demographics

The clinicopathological characteristics of the 348 patients with lung adenocarcinoma included in this study are summarized in TableĀ 1. Of the 348 patients, 306 (87.9%) had pN0 disease, and 42 (12.1%) had nodal metastasis. The pN0 disease proportion was 83.6% (122/146) in males and 91.1% (184/202) in females (Pā€‰=ā€‰0.033).

Table 1 Patient characteristics.

Correlation between pN0 disease proportion and pathological features

The pN0 disease proportion in histological Grade 1, Grade 2, Grade 3 and Grade 4 was 100%, 97.7%, 82.5% and 73.3%, respectively (TableĀ 2). The pN0 disease proportion in the well-differentiated group, moderately differentiated group and poorly differentiated group was 97.7%, 80.6% and 65.5%, respectively (TableĀ 3). Decreasing trends were observed for both histological grade and differentiation grade (Pā€‰<ā€‰0.001).

Table 2 Trend of pN0 disease proportion among different tumor grades.
Table 3 Trend of pN0 disease proportion among different tumor differentiation degrees.

Radiographic features analysis

The univariate analysis of radiographic features of node-negative and node-positive disease is presented in TableĀ 4. Tumor size, size stratification, appearance, location, contour, GGOc proportion, Sc proportion and maximum SUV (SUVmax) differed significantly between pN0 disease and node-positive disease (Pā€‰=ā€‰0.002, Pā€‰=ā€‰0.008, Pā€‰<ā€‰0.001, Pā€‰=ā€‰0.032, Pā€‰=ā€‰0.006, Pā€‰<ā€‰0.001, Pā€‰<ā€‰0.001, Pā€‰=ā€‰0.001, respectively). Specific to quantitative radiographic indexes, the pN0 disease had more GGOc proportion, less Sc proportion (Fig.Ā 1) and lower SUVmax (Fig.Ā 2).

Table 4 Correlation between radiographic features and lymph node status.
Figure 1
figure 1

Images of a 59-year-old female patient with invasive adenocarcinoma (aā€“d) and a 45-year-old male patient with invasive adenocarcinoma (eā€“h). (a,b) The lung window and mediastinal window of the axial HRCT showed a partially solid nodule in the left upper lobe. (c) The CT number distribution curve revealed that the solid component proportion of the nodule was 6%. (d) The photomicrograph (hematoxylin and eosin staining; magnificationā€‰Ć—ā€‰20) revealed a lepidic-predominant invasive adenocarcinoma consisting of 80% lepidic pattern and 20% acinar pattern (pT1N0M0). (e,f) The lung window and mediastinal window of the axial HRCT showed a solid nodule in the right upper lobe. (g) The CT number distribution curve revealed that the solid component proportion of the nodule was 80.3%. (h) The photomicrograph (hematoxylin and eosin staining; magnificationā€‰Ć—ā€‰100) revealed an acinar-predominant invasive adenocarcinoma consisting of 80% acinar pattern and 20% solid pattern (pT1N2M0, right superior mediastinal lymph node metastasis, 3/6; subcarinal lymph node metastasis, 1/1).

Figure 2
figure 2

Images of a 59-year-old female patient with invasive adenocarcinoma (aā€“d) and a 46-year-old female patient with invasive adenocarcinoma (eā€“h). (a) The lung window of the axial HRCT showed a solid nodule in the right lower lobe. (b) PET/CT revealed low metabolic activity (SUVmax, 1.5). (c) PET showed no other abnormal FDG uptake in the whole body. (d) The photomicrograph (hematoxylin and eosin staining; magnificationā€‰Ć—ā€‰100) indicated an acinar predominant invasive adenocarcinoma consisting of 60% acinar pattern, 35% lepidic pattern and 5% papillary pattern (pT1N0M0). (e) The lung window of the axial HRCT showed a solid nodule in the right lower lobe; the nodule had a lobulated and speculated margin. (f) PET/CT revealed high metabolic activity (SUVmax, 7.8). (g) PET showed no other abnormal FDG uptake in the whole body. (h) The photomicrograph (hematoxylin and eosin staining; magnificationā€‰Ć—ā€‰100) indicated a solid predominant invasive adenocarcinoma consisting of 70% solid pattern and 30% acinar pattern (pT1N1M0, parabronchial lymph node metastasis of the right middle lobe, 1/3).

Logistic regression analysis demonstrated that the Sc proportion was only significant predictor of pN0 disease (Ī²ā€‰=ā€‰6.153, Pā€‰=ā€‰0.001; ORā€‰=ā€‰470.085, 95% CI: 12.291~17978.987). With the increase of Sc proportion, the chance of pN0 disease is reduced (Fig.Ā 3). The positive predictive value (PPV) of pN0 disease was 100% for patients with a tumor Sc proportion of no more than 21.3%. Of the 348 patients, only 137 underwent PET/CT, and thus SUVmax was removed from the multivariate analysis.

Figure 3
figure 3

The Positive predictive value (PPV) of pN0 disease in cT1N0M0 lung carcinoma based on the tumor solid component (Sc) proportion.

Discussion

A randomized, prospective clinical trial has suggested that MLND can be avoided if there is sufficient evidence of pN0 disease13. Although multiple studies have attempted to identify clinical, pathological and radiographic factors of lymph node status in early-stage adenocarcinoma6,7,14,15,16,17,18, our study differs from these previous studies in several ways. First, this study made cT1N0M0 lung adenocarcinoma a separate object of study. Second, all resected specimens were reviewed pathologically according to the 2011 IASLC/ATS/ERS classification, and a tumor-grading system was established based on histopathological subtypes according to the newly proposed classification19. Third and most important, quantitative radiographic analysis method was introduced, the Sc and GGOc proportion calculated using three-dimensional computer-quantified CT number analysis was adopted to qualitatively reflect the tumor internal appearance, resulting in more objective and convenient selection criteria for clinical application.

Our study results included three main findings: (a) the Sc proportion was an independent predictor of pN0 disease in cT1N0M0 lung adenocarcinoma; (b) all cT1N0M0 lung adenocarcinoma patients with an Sc proportion of no more than 21.3% had no lymph node metastasis; and (c) the pN0 disease proportion exhibited a decreasing trend with both increasing histological grade and differentiation grade. Our results may have clinical implications. With the aid of preoperative quantitative radiographic features, we differentiated cT1N0M0 lung adenocarcinoma patients with or without pN0 disease, thus reducing unnecessary MLND.

In our study, tumors were graded according to four classifications based on histopathological subtypes and prognosis. Sica et al.20 first proposed a 3-tier grading system for lung adenocarcinoma in 2010: Grade 1, a pattern with low metastatic potential (bronchioloalveolar carcinoma, BAC); Grade 2, patterns with intermediate metastatic potential (acinar and papillary); and Grade 3, patterns with high metastatic potential (solid and micropapillary). However, in 2011, the newly proposed lung adenocarcinoma classification discontinued the use of ā€œBACā€, and other changes in uniform terminology, such as adenocarcinoma in situ (AIS) and minimally invasive adenocarcinoma (MIA), and diagnostic criteria were implemented19. Studies21,22 have shown that this new classification system is a stage-independent predictor of survival. Five-year disease-free survival is 100% or near 100% in patients with AIS and MIA, more than 85% in lepidic predominant invasive adenocarcinoma, 65% in acinar and papillary predominant invasive adenocarcinoma, and 40% in solid and micropapillary predominant invasive adenocarcinoma. The tumor-grading system in this study is based on a combination of the classification of Sica et al. and the lung adenocarcinoma classification system proposed in 2011. A higher grade represents a worse prognosis.

A decreasing trend in the pN0 disease proportion with increasing tumor grade and decreasing tumor differentiation was observed in this study, consistent with the worse prognosis of lung adenocarcinoma patients with higher tumor grade and poorer tumor differentiation. In this study, all cT1N0M0 lung adenocarcinoma patients with grade 1 tumors (AIS or MIA) had pN0 disease, whereas only 97.7% of patients with well-differentiated tumors had pN0 disease. Therefore, we believe that tumor grade has stronger predictive ability for lymph node status than tumor differentiation and that tumor grade 1 is a significant predictor of pN0 disease. However, because of histological heterogeneity, AIS and MIA can only be diagnosed in resected specimens and not in small biopsies or cytology specimens23. Consequently, histopathological subtypes have limited significance in predicting node metastasis preoperatively and determining the most appropriate surgical procedures.

In univariate analysis of radiographic features, tumor size, size stratification, location, contour, appearance, Sc proportion, GGOc proportion and SUVmax were associated with pN0 disease. In this study, all 21 patients with a tumor size of no more than 1ā€‰cm had pN0 disease, in accordance with previous study14. We also found that compared with patients with node-positive disease, primary tumors in patients with pN0 disease are more likely to present a smooth margin and peripheral location. Smooth margin and peripheral location were frequently proven to be benign signs in lung nodule diagnosis. Tumor appearance has been proven to be a significant predictor of lymph node status. Researchers have found that most patients with a dominant GGO component on CT have no lymph node metastasis15. In this study, all 52 patients with pGGO had pN0 disease, consistent with previous studies6,16. However, tumor appearance is only used as a subjective and visual characteristic in the past study, and thus is not sufficiently objective and accurate for prediction of pN0 disease. The great progresses of computer make quantitative imaging increasingly applying in modern radiology practice, assisting in the clinical assessment and providing a source of biomarkers for a spectrum of diseases24. In the current study, the Sc proportion and GGOc proportion calculated by three-dimensional computer-quantified CT number analysis were adopted to qualitatively reflect tumor appearance innovatively, resulting in more objective and accurate selection criteria for pN0 disease. The Sc proportion and GGOc proportion were all associated with lymph node status. Primary tumors in patients with pN0 disease had lower Sc proportions and higher GGOc proportions than patients with node-positive disease. Studies7,17,18 have demonstrated that the solid component is a powerful predictor of lymph node metastasis, consistent with our study. SUVmax is a semiquantitative index of FDG uptake of PET/CT, which is associated with glucose metabolism in tissue. SUVmax was significantly lower in patients with pN0 disease than patients with node-positive disease in this study, consistent with previous studies7,25,26,27.

We analyzed all factors except SUVmax by logistic multivariate analysis to assess the joint effects and interactions of the variables on pN0 disease. Our results showed that only the Sc proportion retained statistical significance, whereas tumor size, size stratification, location, contour, appearance and GGOc proportion did not differ significantly in multivariate analysis. We believe that the effects of these factors were substituted by the enrolled index. When Sc proportion ā‰¤21.3% was used as the criterion, the PPV for pN0 disease was 100%. On the basis of our result, patients of cT1N0M0 lung adenocarcinoma who meet this reliable prediction criterion may omit MLND when surgical treatment.

Our study has several limitations. First, this is a single study; further studies at multiple centers are needed to confirm our results, especially the predictive value of Sc. Second, this is a retrospective study, and a prospective study is needed to verify this conclusion. Finally, we did not perform survival analysis because the follow-up time was too short in some cases; advanced study is warranted to analyze patient survival.

In conclusion, we respectively reviewed 348 cT1N0M0 lung adenocarcinoma patients who underwent surgical pulmonary resection with systematic MLND at our hospital in an attempt to investigate independent predictors of pN0 disease. The Sc proportion was identified as statistically significant predictive risk factor. The PPV of pN0 disease was 100% for patients with a tumor Sc proportion of no more than 21.3%. Patients who meet this criterion may be managed with lung resection without MLND.

Materials and Methods

Patients

This prospective study was approved by the institutional review board at the Cancer Hospital of the Chinese Academy of Medical Sciences, and the informed consent was waived due to the retrospective nature of the study and the data were going to be analyzed anonymously. This study was conducted in accordance with the Declaration of Helsinki.

A total of 348 consecutive patients with surgically resected cT1N0M0 lung adenocarcinoma at our hospital between January 2005 and December 2012 were included. The inclusion criteria for our study were as follows: (a) single adenocarcinomas (3ā€‰cm or less in diameter at HRCT) with no evidence of malignant satellite nodules (as proven by prior imaging study or lung biopsy); (b) no hilar or mediastinal lymphadenopathy on imaging study or at mediastinoscopy; a short lymph node diameter greater than 1ā€‰cm on CT axial image with window level in 380 HU and window within 40 HU or maximum SUV (SUVmax) on PET/CT image more than 2.5 was determined as lymphadenopathy; (c) first treatment with surgical pulmonary resection plus systematic MLND (more than 6 hilar or mediastinal lymph nodes resected); (d) either integrated 18F fluorodeoxyglucose (18F-FDG) PET/CT or chest HRCT studies acquired within 1 month before resection; (e) both pathological sections and clinical data were available for review.

Imaging and Interpretation

All 348 patients included in this study had undergone chest HRCT for preoperative staging. HRCT images were obtained with an eight- (LightSpeed Ultra, GE Medical Systems), 16- (ProSpeed or Discovery ST, GE Medical Systems) or 64- (LightSpeed VCT, GE Medical Systems or Toshiba Aqulion, TOSHIBA Medical Systems) slice spiral CT scanner. Unenhanced CT images were obtained at 120ā€‰kVp and 250~350ā€‰mA with a reconstruction kernel with a standard algorithm. The reconstruction thickness was 1.0 or 1.25ā€‰mm, and the interval was 0.8 or 1.0ā€‰mm. A total of 182 patients underwent enhanced HRCT examination, in which 60 to 80ā€‰ml of intravenous contrast was administered at 2.0 to 2.5ā€‰mL/s and enhanced images were obtained 25 to 30ā€‰seconds after contrast infusion.

HRCT images were assessed retrospectively for morphological features and three-dimensional computer-quantified CT number analysis. Two radiologists (M.L. and L.Z., 11 and 8 years of experience, respectively), who were informed that the involved patients had surgically treated adenocarcinoma but were blinded to histological subtypes, analyzed the morphological features on Carestream GCRIS 2.1 PACS workstation (Carestream Health) in consensus. The morphological features included tumor size (the longest tumor diameter on the transverse lung window image in which the largest nodule dimension appeared), appearance (solid nodule [SN], part-solid nodule [PSN], or pure ground-glass-opacity [pGGO])28, location (center or periphery, periphery was defined as within 3ā€‰cm of the pleura)29, and contour (smooth, lobular or spiculated), necrosis (necrosis of tumor was defined as low attenuation in the tumor). Three-dimensional computer-quantified CT number analysis was performed on an ADW 4.6 workstation (GE Medical Systems) to measure the proportion of GGOc and Sc. This measurement was performed in two steps. For the first step, the entire tumor mass was separated from surrounding anatomic structures using computer-aided volume measurement software (Auto Contour in Volume Rendering, GE Medical Systems). The computer generated tumor boundaries were then visually inspected by a radiologist (L.Z.) for accuracy and consistency. If any segmentation results were considered suboptimal, the tumor contours that were superimposed on the original images were edited by the same radiologist (L.Z.), and then the mean CT value of tumor was automatically calculated by computer. Threshold values of GGOc and Sc were obtained by receiver operating characteristic (ROC) curves. As the second step, the entire tumor mass was separated using the same method in the first step, and the computer automatically generated the CT number histogram of the tumor. The proportion of GGOc/Sc was calculated based on the threshold values obtained in the first step. The calculation methods were also described in detail in our previous report30.

Two chest radiologists with PET/CT diagnostic experience prospectively evaluated the integrated PET/CT images. For semiquantitative analysis of 18F-FDG uptakes, a region of interest (ROI) was placed over the tumor site on the hottest trans-axial slice. In some patients, nodular FDG uptake could not be identified on the PET component images of the PET/CT study. In those patients, a ROI was drawn in a presumed nodular location, taking into account the CT component images of PET/CT. The SUVmax within the ROI was used as the reference measurement.

Pathologic Evaluation

All resected specimens were reviewed pathologically by an experienced lung pathologist (S.W., 8 years of experience in lung pathology) according to the 2011 International association for the study of lung cancer/American thoracic society/European respiratory society (IASLC/ATS/ERS) classification19. For difficult cases, histological subtypes were assessed by consultation with pathological experts. Based on the reported prognosis of lung adenocarcinoma31,32,33,34, tumors were graded according to four classifications. Grade 1 included histological subtypes of adenocarcinoma in situ (AIS) and minimally invasive adenocarcinoma (MIA). Grade 2 corresponded to the lepidic pattern of invasive adenocarcinoma. Grade 3 corresponded to tumors that showed acinar or papillary patterns. Grade 4 corresponded to tumors that showed micropapillary or solid patterns. When tumors were composed predominantly of a variant pattern, the tumors were removed from this classification system because they were uncommon and heterogeneous in terms of biological behavior and prognosis. Considering the proximity of tumor cells and normal cells under the microscope, tumors were divided into the well-differentiated group, moderately differentiated group and poorly differentiated group.

Statistical analysis

The trends of pN0 disease proportion among the different tumor grades and tumor differentiation degrees were assessed by chi-square trend test. The prevalences of nominal variables (e.g., tumor appearance, location, contour and intratumoral necrosis) were compared using Fisherā€™s exact test. Differences in the mean values of continuous variables (e.g., tumor size, GGOc proportion, Sc proportion and SUVmax) were compared using the Independent-Studentā€™s t test. Multivariate logistic regression analysis was used to identify significant predictors of pathologic node-negative disease. All statistical analyses were performed using a commercial software package (SPSS, Inc., an IBM Company, Chicago, IL, USA). A P value of less than 0.05 was considered to indicate a significant difference.

Data availability

All data generated or analysed during this study are included in this published article (and its Supplementary Information files).