ABSTRACT
CONCLUSION
The nomogram integrating the SUVmax and the CT features could non-invasively predict VPI status before surgery in subpleural clinical stage IA lung adenocarcinoma patients.
RESULTS
The size of the solid component, the consolidation-to-tumor ratio, the solid component pleural contact length, the SUVmax, the density type, the pleural indentation, the spiculation, and the vascular convergence sign demonstrated significant differences between VPI-positive (n = 40) and VPI-negative (n = 58) cases on univariate analysis in the training set. A multivariate logistic regression model incorporated the SUVmax [odds ratio (OR): 1.753, P = 0.002], the solid component pleural contact length (OR: 1.101, P = 0.034), the pleural indentation (OR: 5.075, P = 0.041), and the vascular convergence sign (OR: 13.324, P = 0.025) as the best combination of predictors, which were all independent risk factors for VPI in the training group. The nomogram indicated promising discrimination, with an area under the curve value of 0.892 [95% confidence interval (CI), 0.813–0.946] in the training set and 0.885 (95% CI, 0.748–0.962) in the validation set. The calibration curve demonstrated that its predicted probabilities were in acceptable agreement with the actual probability. The decision curve analysis illustrated that the current nomogram would add more net benefit.
METHODS
A total of 140 patients with subpleural clinical stage IA peripheral lung adenocarcinoma were recruited and divided into a training set (n = 98) and a validation set (n = 42), according to the positron emission tomography/CT examination temporal sequence, with a 7:3 ratio. Next, VPI-positive and VPI-negative groups were formed based on the pathological results. In the training set, the clinical information, the SUVmax, the relationship between the tumor and the pleura, and the CT features were analyzed using univariate analysis. The variables with significant differences were included in the multivariate analysis to construct a prediction model. A nomogram based on multivariate analysis was developed, and its predictive performance was verified in the validation set.
PURPOSE
Preoperative prediction of visceral pleural invasion (VPI) is important because it enables thoracic surgeons to choose appropriate surgical plans. This study aimed to develop and validate a multivariate logistic regression model incorporating the maximum standardized uptake value (SUVmax) and valuable computed tomography (CT) signs for the non-invasive prediction of VPI status in subpleural clinical stage IA lung adenocarcinoma patients before surgery.
Main points
• 18F-fluorodeoxyglucose positron emission tomography/computed tomography is valuable in the non-invasive diagnosis of visceral pleural invasion (VPI) in subpleural clinical stage IA lung adenocarcinoma patients.
• The maximum standardized uptake value (SUVmax), the solid component pleural contact length, the pleural indentation, and the vascular convergence sign are independent predictors of VPI.
• A nomogram incorporating SUVmax and independent computed tomography features improves predictive performance.
The Global Cancer Statistics 2020 estimate that lung cancer incidence is 11.4%, with lung cancer being the world’s second most prevalent type of malignant tumor and leading cause of cancer death, with an estimated 1.8 million deaths (18%) worldwide in 2020.1 Lung adenocarcinoma accounts for the largest proportion of non-small cell lung cancer (NSCLC).1
Visceral pleural invasion (VPI) in lung cancer is defined as a tumor that invades beyond the visceral pleural elastic layer (PL1) and exposes the pleural surface (PL2); it is one of the poor prognostic factors.2 According to the eighth edition of the tumor-node-metastasis (TNM) staging criteria, when cT1N0M0 lung adenocarcinoma has pathologically positive VPI, the tumor (T) stage is upgraded from T1 to T2, and the TNM stage is upgraded from IA to IB.2
According to recent studies, lobectomy surgery results in a better prognosis and lower intrathoracic recurrence for patients with clinical T1 stage tumors (including T1a) and pathologically confirmed positive VPI.3 For patients with T1-sized VPI-positive lung cancer, more extensive lymph node dissection rather than lymph node sampling is required.4 However, the conventional intraoperative diagnosing of VPI is time-consuming and inexact, with a reported accuracy of 56.5%.5 The gold standard for diagnosis still relies on postoperative elastic fiber staining to evaluate VPI,2 but postoperative pathology cannot guide preoperative surgical planning. Therefore, if the VPI status can be accurately assessed using preoperative imaging methods, it may significantly impact the surgical plan.
Previous studies have explored using computed tomography (CT) findings to predict VPI status.6,7,8,9,10,11,12,13,14,15 Subpleural lung adenocarcinoma presenting as pure ground-glass nodules (pGGNs) do not enter the visceral pleura due to their low invasiveness.11,12,13 In addition, VPI-positive status is never observed in lung cancer without a relationship with the pleura.13,14 However, it is challenging to accurately determine the VPI status of lung cancers related to the adjacent pleura on CT, such as pleural tags or those in direct contact with the pleura. Positive predictive values range from 44.1% to 56.4%, indicating that approximately half of the CT-based predictions are false positives.15 Therefore, finding novel and quickly available methods to improve efficiency is critical.
18F-fluorodeoxyglucose (18F-FDG) positron emission tomography (PET)/CT can reflect the glucose metabolism, the morphological characteristics of the tumors, and the anatomical relationship between adjacent structures, which is of great value for the TNM staging of tumors.16 Recent studies have indicated that 18F-FDG PET/CT has high accuracy in predicting the invasiveness,17 the lymph node metastasis,18 and the spread through air spaces19 of clinical T1 stage NSCLC before surgery. However, establishing and validating a multivariate logistic regression model incorporating the maximum standardized uptake value (SUVmax) and the relevant CT signs for predicting the VPI status of clinical stage IA lung adenocarcinoma, excluding pGGNs and tumors unrelated to the pleura, is rare. A nomogram is based on multivariate regression analysis, integrating multiple predictors to transform complex regression equations into visual graphics, making the results of the prediction model more readable and convenient for evaluating patients. Because a nomogram is intuitive and easy to understand, it has been widely used in preoperative nodal staging, predicting invasiveness, and giving a lung cancer prognosis.20,21,22 Therefore, this study aims to explore whether tumor metabolism information and CT features could accurately predict VPI in clinical stage IA lung adenocarcinoma based on a nomogram.
Methods
This retrospective study was approved by the Shanghai Changzheng Hospital Ethical Review Board of the hospital (decision number: CZ-20210528-01). Due to the retrospective nature of this study, informed consent was waived.
Patients
Between February 2015 and June 2022, 328 patients with clinical stage IA lung adenocarcinoma who underwent 18F-FDG PET/CT examination before surgery and surgical resection were recruited. The inclusion criteria were as follows: (i) tumors clinically diagnosed as clinical T1 stage with a tumor size smaller than 30 mm23 and (ii) tumors without pathological lymph node or distant metastasis. The exclusion criteria were as follows: (i) postoperative pathological diagnosis of atypical adenomatous hyperplasia (n = 1) or adenocarcinoma in situ (n = 1); (ii) no report of VPI status (n = 20); (iii) pGGNs (n = 14); (iv) neither directly in contact with the pleural surface nor had pleural tags (n = 107); (v) the minimum distance between the lesion and the pleura (DLP) >10 mm (n = 24); (vi) had treatment prior to 18F-FDG PET/CT examination (n = 5); (vii) had a surgery and an examination interval of longer than 14 days (n = 6); and (viii) poor image quality or incomplete clinical data (n = 10). A total of 140 patients were included in our research (Figure 1). The patients were classified into a training and validation set with a 7:3 ratio based on the PET/CT examination temporal sequence. The training set included 98 patients between February 2015 and January 2020; the temporal validation set comprised 42 patients between April 2020 and June 2022.
Equipment and parameters
Whole-body PET/CT tomography was performed using a Siemens Biograph Truepoint 64 PET/CT. The fasting blood glucose level of the patient was less than 10 mmol/L before the examination. After measuring the patient’s body weight, the patient was given an intravenous injection of 18F-FDG (Shanghai Atomic Kexing Pharmaceutical Co., Ltd.) 3.70–5.55 MBq/kg (0.10–0.15 mCi/kg) body weight with a radiochemical purity >95%, followed by a 300 mL water drink. The patient was instructed to lie down and rest for 60 min in a dark room. After emptying the bladder, the body scan ranged from the base of the skull to the middle of the femur, scanning 5–7 beds, 2–3 min/bed, with a reconstruction matrix of 192 × 192. The PET images were attenuated by CT images and reconstructed iteratively. The body CT scanning parameters were as follows: tube voltage 120 kV, tube current 160 mAs, scanning layer thickness 3.75 mm, reconstruction matrix 512 × 512, and pitch 0.8 s. The chest high-resolution computed tomography (HRCT) scan parameters were as follows: tube voltage 120 kV, tube current 150 mAs, scanning layer thickness 5 mm, reconstruction layer thickness and layer interval 1 mm, reconstruction matrix 512 × 512, and pitch 0.8 s. The lung window images were reconstructed by a high-resolution algorithm (B70f), and the mediastinal window images were reconstructed using a standard algorithm (B40f).
Clinical and pathological data collection
Patients’ clinical data, including age, gender, smoking history, preoperative carcinoembryonic antigen (CEA) level, surgical type, histological subtype, and tumor location, were reviewed.
The resected tissues were stained with hematoxylin and eosin, and the pathological diagnoses were performed by two pathologists with at least 10 years of experience. The specific elastic fiber stain was performed if the VPI status could not be determined accurately. Additionally, VPI was classified as no pleural invasion beyond the elastic layer (PL0), tumor invasion beyond the elastic layer (PL1), and tumor invasion to the surface of the visceral pleura (PL2), with PL1 and PL2 indicating the presence of VPI.2
Image evaluation
The 18F-FDG PET/CT images were imported into the software (RadiAnt DICOM Viewer 4.2.1, Medixant, Poland) and analyzed by two independent radiologists with seven years of experience who were blinded to the pathological information. The lung window [width: 1.500 Hounsfield scale (HU), level: −500 HU], mediastinal window (width: 300 HU, level: 50 HU), multiplanar reformation (MPR), and maximal intensity projection were used to analyze the lesion. For quantitative indicators, the average measurements of two independent radiologists were used as the final data. For qualitative analysis, disagreements were discussed until a consensus was reached.
First, the tumor density type was classified as solid or part solid. The tumor size (the longest length of the tumor, T) and the solid component size (the longest length of the consolidation part, C) were measured at the lung window on the MPR images, and the proportion of the consolidation part was calculated (C/T ratio, CTR).24
Second, the tumor–pleura relationship was classified as a pleural attachment (directly contacting the pleura) or pleural tags (without abutting the pleura). Pleural tags were defined as one or multiple high-density linear strands connecting the tumor margin and the pleura (Figures 2-4).8 The minimum vertical DLP was measured on the MPR images at the lung window for tumors with pleural tags (Figures 5, 6).6 The pleural indentation sign was defined as the deviation of the pleura from its original position due to tumor traction at the lung window, which can be observed in tumors with pleural tags or pleural attachment (Figures 2, 7).14 The longest interface length of the whole tumor and solid component was measured for tumors with pleural attachment by drawing a straight line at the lung window on the MPR images (Figure 7).10 The solid interface length was 0 mm for a part-solid nodule without the solid component contacting the tumor.
Third, the presence or absence of lobulation, spiculation, vascular convergence, and air bronchogram signs were analyzed for all tumors. Lobulation is defined as a petaloid or wavy appearance at the tumor’s margins. Spiculation refers to short, thin linear strands radiating around the surface of the tumor. Vascular convergence is the convergence of pulmonary vessels around the tumor towards the lesion (Figure 3).22 Air bronchogram refers to air-filled bronchus manifesting as natural, dilated/distorted, or cut-off within the lesions.25
Fourth, the 18F-FDG metric was measured by setting the region of interest covering the tumor on the PET/CT fusion images slice by slice and automatically generating the SUVmax.
Statistical analysis
R software (version 4.1.0, http://www.Rproject.org) and IBM SPSS Statistics (version 20.0, USA) software were utilized for the statistical analysis. The Shapiro–Wilk test was used to examine the normality of numeric variables. The normally distributed numerical variables were represented as the mean ± standard deviation, and the comparison between the two groups was carried out by using the two independent samples t-test. Non-normally distributed data were described as the median and the 25% and 75% quartiles, and the Mann–Whitney U test was performed. Pearson’s chi-squared test or Fisher’s exact test was used for categorical variables analysis. Variables significantly different (P < 0.05) in the univariate analysis were involved in the multivariate analysis of logistic regression, and a backward stepwise selection was applied by using the likelihood-ratio test with the Akaike information criterion (AIC) as the stopping rule to select the best combination of variables to build the prediction model in the training set and the corresponding nomogram. Spearman’s rank correlation was used to analyze the correlation between the tumor size, the solid component size, the CTR, and the SUVmax. The interobserver agreement of numeric and categorical variables was assessed using the intraclass correlation coefficient (ICC) and κ-statistic, respectively. The receiver operating characteristic (ROC) curve with the corresponding area under the curve (AUC) value was used to evaluate the discrimination ability of the prediction model and each risk factor in predicting VPI in the training and validation set. The calibration curve and Hosmer–Lemeshow test were used to evaluate the goodness-of-fit of the prediction model, and a P value of greater than 0.05 indicated a good goodness-of-fit. The decision curve analysis (DCA) was used to evaluate the clinical utility of the nomogram. Multivariate binary logistic regression, nomogram, validation, and calibration plots were done with the “rms” package of R software. The ROC was performed by the “pROC” package, and the DCA was performed with the function of “ggDCA”.
Results
Baseline characteristics of the study cohorts
Among the 140 clinical stage IA lung adenocarcinomas, 57 cases were VPI-positive, and 83 were VPI-negative. The baseline characteristics of the training and validation sets were similar (P > 0.05) (Table 1).
Clinical and 18F-fluorodeoxyglucose positron emission tomography/computed tomography features by visceral pleural invasion status in the training group
Clinical features showed no significant differences in age (P = 0.460), gender (P = 0.359), tumor location (P = 0.453), smoking status (P = 0.349), and CEA level (P = 1.000) between the VPI-negative and VPI-positive groups in the training cohort (Table 2).
For 18F-FDG PET/CT characteristics, the consistency of measurements between the two observers was good (ICCs ranged from 0.962 to 0.998), and the consistency of qualitative evaluation indicators was strong (a Kappa value of 0.879 to 0.971). The interobserver agreement assessment results of each index are shown in Supplementary File 1. The Spearman correlation analysis showed no linear correlation between the tumor size and the SUVmax. The solid component size was positively correlated with the SUVmax (rs: 0.721, P < 0.001), while the CTR was positively correlated with the SUVmax (rs: 0.742, P < 0.001). The univariate analysis showed that the SUVmax, the CTR, the solid component size, the solid pleural contact length, the density type, the pleural indentation, the spiculation, and the vascular convergence sign significantly differed between the VPI-positive and the VPI-negative groups in the training set (P < 0.05). The solid nodules, the pleural indentation, the spiculation, and the vascular convergence signs were more common in the VPI-positive group (P < 0.001). The VPI-positive group presented a significantly higher SUVmax, larger solid component size, greater CTR, and longer solid pleural contact length than the VPI-negative group (P < 0.05; Table 2).
Nomogram development and evaluation
Variables significantly different in the univariate analysis were involved in the multivariate logistic regression analysis. Based on the principle of AIC value minimization, a prediction model was constructed by a backward stepwise selection of variables, including the SUVmax, the solid component pleural contact length, the pleural indentation, and the vascular convergence sign as the best combination of prediction variables. The SUVmax [odds ratio (OR): 1.753, 95% confidence interval (CI) 1.232–2.496, P = 0.002], the solid component pleural contact length (OR: 1.101, 95% CI 1.007–1.204, P = 0.034), the pleural indentation (OR: 5.075, 95% CI 1.065–24.172, P = 0.041), and the vascular convergence sign (OR: 13.324, 95% CI 1.379–128.691, P = 0.025) were independent risk factors for VPI (Table 3). Figures 2 to 8 show the representative cases.
Based on the regression coefficients of the variables included in this model, a nomogram was constructed to evaluate the VPI risk intuitively (Figure 9). The sensitivity, specificity, accuracy, and the AUC for the prediction model in the training set were 82.5%, 79.31%, 80.61%, and 0.892 (95% CI 0.813–0.946), respectively, using the optimal cut-off value of 0.35. The prediction model achieved good discrimination performance in the validation set with sensitivity, specificity, accuracy, and AUC values of 100%, 76.00%, 85.71%, and 0.885 (95% CI, 0.748–0.962), respectively. Figures 10 and 11 present the ROCs for the training and validation cohorts. Table 4 shows the cut-off values and predictive performance of each independent risk factor for predicting VPI in the training and validation cohort. The calibration curve demonstrated that the predicted probabilities were in acceptable agreement with the actual probabilities for the training and validation cohorts, and the Hosmer–Lemeshow test showed good goodness of fit, with P values of 0.648 and 0.051, respectively (Figures 12, 13). The DCA showed that the prediction model adds more net benefit than the “treat all” or “treat none” approach (Figures 14, 15). Supplementary Figures 1, 2 give an example of the clinical application of this nomogram.
Discussion
The present study developed a prediction model as a non-invasive method to detect VPI status in clinical stage IA lung adenocarcinoma. The nomogram incorporating the SUVmax, the solid component pleural contact length, the pleural indentation, and the vascular convergence sign provided a scoring system to help predict VPI status before surgery and had excellent discrimination and acceptable calibration, with an AUC of 0.892 in the training set and 0.885 in the validation set.
The visceral pleura contains abundant lymphatic channels, and the subpleural lymphatic system (peripheral lymphatic system) communicates with the axial lymphatic system distributed around the bronchi and pulmonary vessels.26 Therefore, lung cancer patients with VPI-positive status have a higher probability of hilar and mediastinal lymph node metastasis via the axial lymphatic system27 and are more prone to skip N2 lymph node metastasis.28 Thus, lobectomy and extensive lymph node dissection are required for VPI-positive patients.3,4 Therefore, accurately predicting VPI before surgery is of great clinical significance.
For qualitative indicators, the present study found that the spiculation sign was more common in the VPI-positive group. Spiculation is associated with tumor cell infiltration into adjacent blood and lymphatic vessels, suggesting lung adenocarcinoma is more aggressive.29 However, the spiculation sign was not an independent risk factor for VPI prediction in the present study, similar to the previous literature report.30 As previously reported,6,31 VPI was significantly related to the pleural indentation in the authors’ study. The mechanism could be that as the tumor grows, the development of reactive fibrous hyperplasia of the tumor, to a certain extent, may cause the adjacent pleura to retract, thereby increasing the risk of VPI.32 In addition to the pleural indentation sign, the vascular convergence sign was an independent risk factor in the present study, and no previous study reported similar results. It was speculated that it might be related to the pathological basis of the vascular convergence sign, which was also caused by the reactive fibrous hyperplasia of the tumor.33 As tumor invasiveness increased, reactive fibrous hyperplasia increased, and traction on adjacent pulmonary vessels or the pleura became more obvious.
For quantitative indicators, there was no linear correlation between the tumor size and the SUVmax, according to the Spearman correlation analysis, while the solid component size and the CTR had a relatively high degree of positive correlation with the SUVmax. This may be related to the larger proportion of part-solid nodules (78.57%) in the training set included in the present study. The tumor cells in the solid component area of the part-solid nodules were denser than in the ground glass area, so the 18F-FDG uptake value was greater in the solid component area. Additionally, the SUVmax, the solid component size, the CTR, and the solid pleural contact length were significantly greater in the VPI-positive group (P < 0.05) than in the VPI-negative group, but the tumor size and the whole tumor pleural contact length did not show a difference between the two groups. The results indicate that the solid component represents the more aggressive part of the tumor, consistent with clinical T staging depending on the solid component’s size rather than the tumor’s overall size, in tumors less than 3 cm.23 In this study, the analysis found that the SUVmax and the solid pleural contact length were independent risk factors for predicting VPI among quantitative indicators, similar to previous studies.34,35 This result suggests that whether a solid component contacts the pleura and its contact length should also be an important consideration when determining VPI status. However, solid components in tumors may also correspond to fibroblast proliferation, alveolar collapse, inflammatory cell infiltration, and tumor-secreted mucus.36 In some cases in the authors’ study, although the solid component accounted for a large proportion, the SUVmax value was low, and the adjacent visceral pleura was not invaded (Supplementary Figures 1 to 2). Therefore, diagnostic efficacy for VPI status is limited when only considering the solid component size and the CTR. The SUVmax reflected the 18F-FDG metabolism of tumors, representing the activity of tumor cell proliferation, and was a better indicator of tumor invasiveness in previous studies.16,17,18,19 Furthermore, previous results found no VPI-positive status in clinical IA lung adenocarcinoma with the SUVmax less than 1.3,34 indicating that an extremely low SUVmax can represent less invasiveness, which has a high negative predictive value.34 Similar to the previous study, the authors’ results also found that VPI-positive status was not observed in patients with an SUVmax of less than 1.4 in the training group. However, contrary to previous reports,6,30 there was no significant difference in the DLP between the two groups. This may be because the proportion of tumors directly in contact with the pleura (pleura attachment) included in this study was comparable between the VPI-positive and VPI-negative groups in the training set, while the type of tumors in direct contact with the pleura (DLP = 0 mm) was more prevalent in the VPI-positive group in previous studies.6,30 The different results are due to selection bias, and further verification is required after collecting large sample cases in the future.
There were two previous studies on predicting VPI status in lung adenocarcinoma based on 18F-FDG PET/CT, but both of them only analyzed the risk factors related to VPI-positive status and did not validate the model. Tanaka et al.34 studied the value of HRCT and 18F-FDG PET/CT in predicting the pleural invasion of lung adenocarcinoma in direct contact with the pleura. In the subgroup analysis, including whole tumors size less than 3 cm, the multivariate analysis showed that the SUVmax and the whole tumor contact length with the pleura were independent risk factors for predicting pleural invasion. The predictive performance of only the SUVmax (AUC value 0.844) was comparable to the multivariate model (AUC value 0.845–0.857). Unlike the present study, the previous study included pGGNs and 11 cases with pathological parietal pleural invasion. The study by Chen et al.35 focused on subsolid nodules, including pGGNs, and the multivariate analysis revealed that the DLP and the SUVmax were independent risk factors for predicting VPI. The AUC value of the combined model, including the two variables, was 0.90. However, Chen et al.35 did not indicate whether the tumors included were related to the pleura, especially in the VPI-negative group. The present study excluded tumors presenting as pGGNs and tumors unrelated to the pleura that do not invade the visceral pleura, unlike in previous reports.11,12,13 Therefore, the authors’ results were more objective in including cases with potential VPI-positive status to study the predictive performance of the SUVmax combined with the CT features. The different results may be due to the different inclusion criteria of the study.
A previous study showed that 18F-FDG liver metrics (the SUVmax and the SUVmean) obtained from PET/CT studies of the same patient scanned on the same machine at different time periods were highly repeatable and might reliably be used in following patients longitudinally.37 The SUVmax values in the present study were all obtained from a single 18F-FDG PET/CT device in the same hospital; standardized image acquisition and post-processing were adopted, avoiding the differences brought by the hardware and post-processing software of the PET/CT equipment. It is worth noting that standardized image acquisition and post-processing of 18F-FDG PET/CT and accurate measurement of the SUVmax are key factors in the preoperative assessment of whether lung cancer has invaded the visceral pleura.
This study has several limitations. First, the study was retrospective, resulting in inevitable selection bias. Second, although internal validation of the model yielded good discrimination, the generalization of this nomogram needs to be verified on external data. Third, VPI status in some cases lacks clear PL1 and PL2 pathological grades. Therefore, the clinical application and generalization of the model still need to be further improved and validated by multicenter studies.
In conclusion, the nomogram incorporating the SUVmax and the independent CT features (including the solid component pleural contact length, the pleural indentation, and the vascular convergence sign) showed good value for predicting VPI status in clinical stage IA lung adenocarcinoma in this preliminary study. As a non-invasive quantitative method, 18F-FDG PET/CT holds the potential to help in the surgical decision-making process and additional treatment of clinical stage IA lung adenocarcinoma.