Skip to main content

Nomogram using intratumoral and peritumoral radiomics for the preoperative prediction of visceral pleural invasion in clinical stage IA lung adenocarcinoma

Abstract

Background

Accurate prediction of visceral pleural invasion (VPI) in lung adenocarcinoma before operation can provide guidance and help for surgical operation and postoperative treatment. We investigate the value of intratumoral and peritumoral radiomics nomograms for preoperatively predicting the status of VPI in patients diagnosed with clinical stage IA lung adenocarcinoma.

Methods

A total of 404 patients from our hospital were randomly assigned to a training set (n = 283) and an internal validation set (n = 121) using a 7:3 ratio, while 81 patients from two other hospitals constituted the external validation set. We extracted 1218 CT-based radiomics features from the gross tumor volume (GTV) as well as the gross peritumoral tumor volume (GPTV5, 10, 15), respectively, and constructed radiomic models. Additionally, we developed a nomogram based on relevant CT features and the radscore derived from the optimal radiomics model.

Results

The GPTV10 radiomics model exhibited superior predictive performance compared to GTV, GPTV5, and GPTV15, with area under the curve (AUC) values of 0.855, 0.842, and 0.842 in the three respective sets. In the clinical model, the solid component size, pleural indentation, solid attachment, and vascular convergence sign were identified as independent risk factors among the CT features. The predictive performance of the nomogram, which incorporated relevant CT features and the GPTV10-radscore, outperformed both the radiomics model and clinical model alone, with AUC values of 0.894, 0.828, and 0.876 in the three respective sets.

Conclusions

The nomogram, integrating radiomics features and CT morphological features, exhibits good performance in predicting VPI status in lung adenocarcinoma.

Peer Review reports

Background

Lung cancer ranks among the most prevalent and lethal malignant tumors [1], with lung adenocarcinoma (LUAD) constituting the predominant pathological subtype [2]. Visceral pleural invasion (VPI) is strongly related to adverse postoperative prognoses [3]. The presence of VPI elevates the T stage of lung cancer from T1 to T2 and advances the TNM stage from IA to IB [3]. Notably, recent research underscores the significant influence of VPI status on the decision-making process of thoracic surgeons regarding the selection of an appropriate surgical approach for lung cancer [4]. In cases where VPI is absent, sublobectomy serves as a viable option to preserve lung function.3 Conversely, if the tumor does invade the visceral pleura, a lobectomy coupled with a more extensive lymph node dissection may yield a more favorable prognosis [5]. Unfortunately, the conventional intraoperative diagnosis of VPI is time-consuming and prone to inaccuracies, with a reported accuracy of only 56.5% [6]. Consequently, the current gold standard, relying on postoperative pathological diagnosis [2], introduces a delay in determining an appropriate surgical plan, highlighting the need for a reliable preoperative predictive method for VPI status to ensure precise treatment of lung cancer.

Previous studies have used CT features to predict VPI status [7,8,9,10,11,12,13,14,15,16]. Notably, lung cancer presenting as a pure ground-glass nodule (pGGN) on CT or with no contact with the pleura was associated with a VPI-negative status [12,13,14,15]. Conversely, direct tumor-pleural contact, pleural indentation sign, and pleural tags were identified as high-risk CT features for predicting VPI status in LUAD [16]. However, the accuracy of VPI prediction based on different combinations of these high-risk CT signs ranged from 62.7 to 72.3%, with positive predictive values ranging from 44.1 to 56.4% [16]. This suggests a substantial proportion of false positive predictions, highlighting the ongoing challenge in accurately determining whether a part-solid or solid nodule in contact with the pleura has invaded the visceral pleura based on CT morphological features alone.

Radiomics, a method for extracting high-dimensional features from segmented images, including gray level changes and voxel spatial relationships, holds huge promise in achieving accurate diagnoses and prognosis assessments of diseases through feature selection and model establishment [17,18,19]. By leveraging radiomics, the subjective interpretation of CT morphological features by observers can be circumvented, and a wealth of digital information within the image, imperceptible to the human eye, can be comprehensively mined and integrated [20]. Previous pathological studies have indicated the presence of tumor-infiltrating lymphocytes and tumor-associated macrophages on the edge of invasive lesions [21], demonstrating associations with metastasis [22]. The peritumoral region, as an indicator of the tumor microenvironment, holds critical biological significance in reflecting the aggressive behavior of the tumor [23, 24]. Recent studies have demonstrated that incorporating radiomics features from the peritumoral region into modeling analysis can enhance predictive accuracy in the preoperative assessment of pathological invasiveness [17], lymphovascular invasion [18], lymph node metastasis [25], and spread through air space in lung cancer patients, compared to models relying solely on intratumoral features [26]. While previous investigations have successfully applied radiomics for assessing VPI status in early lung cancer [27,28,29,30,31,32], their primary focus has been on intratumoral features, with limited exploration of the potential contribution of peritumoral radiomics features. Additionally, the reliability and reproducibility of these models have not been validated in external sets.

In this study, we hypothesized that high-throughput feature extraction within the volume of interest (VOI), encompassing both intratumoral and peritumoral regions, can not only capture the tumor’s intrinsic heterogeneity but also comprehensively depict the spatial structural characteristics of the tumor and its adjacent tissues in a higher-dimensional space. Consequently, we developed a radiomics model based on the GTV and the incorporation of peritumoral regions at 5 mm, 10 mm, and 15 mm distances. Our objective was to assess the radiomics model’s effectiveness in predicting VPI status in early subpleural LUAD and to investigate whether its diagnostic performance could be further enhanced when combined with traditional CT morphological features.

Methods

Patients

A total of 1146 patients with clinical stage IA LUAD, who underwent surgery between July 2014 and July 2022 at three hospitals, were included retrospectively.

Inclusion criteria were as follows: (1) Clinical stage IA LUAD (cT1N0M0, with the maximum tumor length ≤ 3 cm); (2) Tumors located under the pleura, in direct contact with the pleura or connected to the pleura via lines or strips (pleural tags) [9]; (3) Minimum distance from the lesion to the pleura (DLP) ≤ 10 mm; (4) Thin-slice chest CT (slice thickness ≤ 2 mm) performed within 2 weeks before surgery; (5) Complete pathological report on VPI status.

Exclusion criteria were as follows: (1) Pathologically confirmed atypical adenomatous hyperplasia, or adenocarcinoma in situ (n = 71); (2) Tumors neither in direct contact with the pleural surface nor with pleural tags (n = 278); (3) DLP > 10 mm (n = 65); (4) Prior treatment or biopsy before CT (n = 15); (5) Time interval of more than 2 weeks between CT examination and surgery (n = 10); (6) Poor CT image quality (n = 78); (7) Presence of pGGN (n = 123); (8) Unavailability of VPI status pathological report (n = 21).

Ultimately, 404 patients from our institution (Hospital 1) were included as the internal cohort and were randomly divided into a training set (n = 283) and an internal validation set (n = 121) at a ratio of 7:3. An external validation set was formed, comprising patients from Hospital 2 (n = 25) and Hospital 3 (n = 56) (Fig. 1). Patients were categorized into two groups based on the VPI status determined during pathology: VPI-positive and VPI-negative. Specific details about the collected clinical and pathological data can be found in the Additional file 1. This study was approved by the Ethics Committee of our Hospital (decision number: CZ-20210528-01), and subjects’ informed consent was exempted.

Fig. 1
figure 1

The flowchart of inclusion and exclusion criteria of patients

CT morphological features evaluation and clinical model construction

The study’s workflow is depicted in Fig. 2. CT image acquisition details are provided in the Additional file 2. The DICOM images of non-enhanced chest CT scans were imported into the Radiant DICOM Viewer software (Version 4.2.1, Medixant, Poland) for visualization of tumor CT morphological features through multi-planar reconstruction (MPR) and maximum density projection (MIP). Two thoracic radiologists with 7 years and 9 years of working experience blinded to patient information, conducted qualitative assessments and quantitative measurements of the tumors. Consensus was reached through consultation in cases of differing opinions regarding qualitative indicators. The average of the measurements by the two radiologists was used for quantitative parameters analysis.

Fig. 2
figure 2

The workflow of this study

The assessment included the following CT qualitative indicators: (i) Tumor-pleura signs, such as tumors in direct contact with the pleura (pleural attachment) or connected to the pleura via linear strands (pleural tags) [9] (ii) Pleura signs involving tumor-induced deviation of the pleura from its original position (pleural indentation) [15]; and the presence of a solid component of the tumor in direct contact with the pleura (solid attachment); (iii) Tumor signs, encompassing density type (solid or part-solid), shape (irregular or oval/round), interface (ill-defined or well-defined), two marginal characteristics (lobulation and spiculation), two internal features (vacuole sign, cavity or cyst sign), and three adjacent structural features (vascular convergence sign, air bronchogram signs, and emphysema in the lobe of the tumor) [33, 34].

For CT quantitative parameter measurement: (i) The maximum diameter of the tumor (tumor size, T) and the maximum diameter of the solid component (consolidation size, C) were measured on MPR lung window images, and the consolidation-tumor ratio (CTR) was calculated [35]. (ii) For tumors with pleural tags, the minimum vertical DLP at the lung window on MPR images was measured [7]. (iii) For tumors with pleural attachment and a DLP of 0 mm, the longest interface length of the entire tumor and solid component was measured using a straight line on the lung window of the MPR images [11]. The illustration of the measurements of DLP, pleural contact length and solid pleural contact length are shown in Fig. 3.

Fig. 3
figure 3

The illustration of the measurements of DLP, pleural contact length and solid pleural contact length. (A) A 53-year-old man presented with a mGGN in the upper lobe of the right lung (white arrow), with a DLP of 4.64 mm. The pathological results were invasive LUAD with VPI-positive. (B) A 64-year-old man presented with a mGGN in the lower lobe of the right lung (white arrow), which was in direct contact with the interlobar fissure and accompanied by pleural indentation sign, with solid components in contact with the pleura. The green line was a straight line to measure the pleura contact length (14.2 mm), and the red line was the solid pleural contact length (8.51 mm). The pathological results were invasive LUAD with VPI-positive.

CT images preprocessing and tumor segmentation

Prior to segmentation, standard resampling and grayscale discretization of CT images were performed. The ITK-SNAP software (version 3.80, www.itksnap.org) was used for slice-by-slice delineation of the GTV along the tumor boundary, with the 3D region defined as the VOI. GTV was defined as the entire visible tumor area identified subjectively by the radiologists. The segmentation aimed to exclude vessels, bronchi, and pleura as much as possible during the process.

The intraclass correlation coefficients (ICCs) were calculated to assess agreements between and within observers for the selection of radiomics features with high reproducibility and reliability. Inter-observer agreements were evaluated by two thoracic radiologists with 7 years and 9 years of working experience, who independently segmented 30 randomly selected tumors. Intra-observer agreements were assessed by the 9-year experienced radiologist, who re-segmented 30 tumors one month later. The remaining tumor segmentations were performed solely by the 9-year experienced radiologist, with both radiologists blinded to the patients’ pathological information during the segmentation process.

The VOI segmentation of the GTV was expanded outward by 5 mm, 10 mm, and 15 mm, respectively, and pixel filtering was conducted to automatically exclude peritumoral non-pulmonary tissues (such as peritumoral vessels, soft tissues, and chest wall, mediastinum, and abdominal structures) based on pixel thresholding using a Python platform (version 3.11). The resulting VOI segmentation of the GPTV in Nifty format was generated, named GPTV5, GPTV10, and GPTV15, accordingly.

Radiomics feature extraction and data preprocessing

The extraction of features was carried out using Pyradiomics (version 3.0.1, https://pyradiomics.Read the docs.io/en/latest/changes.html). A total of 1218 radiomics features were extracted, with detailed feature types provided in Additional file 3 and 4. Given the inclusion of CT images from multiple hospitals and diverse CT protocols, the intensities of all radiomics features were normalized using the ComBat compensation method (https://forlhac.shinyapps.io/Shiny_ComBat) and z-score transformation (z = x-µ/σ) [36,37,38].

Radiomics feature selection and model construction

ICCs were employed to evaluate the consistency between the segmented GTV and the corresponding GPTV5, GPTV10, and GPTV15 radiomics features for intra-observer and inter-observer assessments. The “psych” package of R software was used to assess the consistency of the radiomics features. To avoid overfitting, the maximal redundancy minimal relevance (mRMR) algorithm and the least absolute shrinkage and selection operator (LASSO) logistic regression method were applied to features with good consistency (ICC > 0.75) in the training set to limit the dimension of the features. A 10-fold cross-validation process was employed to select the optimal hyperparameter λ. Features with coefficients not equal to zero under the optimal λ were chosen for constructing the radiomics model. The radscore was calculated by summing the selected features, weighted by their coefficients. Wilcoxon tests were used to compare differences between the VPI-positive and VPI-negative groups. Radiomics models were developed based on GTV, GPTV5, GPTV10, and GPTV15, and their diagnostic efficacy was evaluated. The best radiomics model was determined by the highest AUC value in the validation set. The analysis process is depicted in Fig. 2.

Clinical model and nomogram construction

In the training set, a multivariate logistic regression analysis was performed using clinical and CT morphological features with P < 0.1 in univariate logistic analysis. The best variable combinations were selected through a backward stepwise selection process to establish the clinical prediction model. The variables in the clinical model and the radscore from the optimal radiomics model were used as predictive variables to create the combined prediction model and the corresponding nomogram. The predictive performance of each model was assessed in both internal and external validation sets. The model development process is illustrated in Fig. 2.

Statistical analysis

IBM SPSS Statistics (version 20.0, USA) and R software (version 4.2.2, http://www.Rproject.org) were used for statistical analysis. Details of the statistical analysis process are presented in Additional file 5.

Results

Clinical prediction model

Additional file 6 provides the clinical and pathological data for the three respective sets. Regarding CT features, there was strong agreement between the measurements of the two observers (ICC 0.954–0.996), and qualitative evaluation indicators exhibited high consistency (Kappa 0.832-1.000). Additional inter-observer agreement results for each CT feature are detailed in Additional file 7. Univariate analysis revealed significant differences (P < 0.05) between the VPI-positive and VPI-negative groups for several CT features, including tumor size, solid component size, CTR, whole tumor pleural contact length, solid component pleural contact length, density type, pleural indentation sign, solid attachment sign, spiculation sign, vascular convergence sign, and the presence of emphysema. Notably, the VPI-positive group exhibited a higher prevalence of solid nodules, solid attachments, pleural indentation, spiculation, and vascular convergence signs, along with significantly larger tumors and solid components, longer CTR, and longer whole tumor and solid pleural contact lengths (all P < 0.05; Table 1).

Table 1 Univariate analysis of CT features of clinical stage IA lung adenocarcinoma

Multivariate logistic regression analysis was employed to select the optimal combination of predictive variables for constructing the clinical model. The independent risk factors for VPI were the solid component size (OR = 1.23, 95% CI 1.16 ~ 1.30, P < 0.001), pleural indentation sign (OR = 3.36, 95% CI 1.61 ~ 7.02, P = 0.001), solid attachment sign (OR = 2.98, 95% CI 1.56 ~ 5.70, P < 0.001), and vascular convergence sign (OR = 4.51, 95% CI 1.49 ~ 13.69, P = 0.008), as shown in Table 2.

Table 2 Univariate and multivariate logistic regression analysis of the CT features in the training set

Radiomics model

A total of 1218 radiomic features were extracted from the VOI of GTV, GPTV5, GPTV10, and GPTV15, respectively. Among these features, 72.2% (880/1218) of GTV features, 92.3% (1124/1218) of GPTV5 features, 97.1% (1183/1218) of GPTV10 features, and 97.8% (1191/1218) of GPTV15 features demonstrated good repeatability, with inter-class and intra-class ICCs exceeding 0.75. Among the features with ICC > 0.75, the mRMR algorithm was initially used to eliminate redundant and irrelevant features, retaining 30 features in each group. Subsequently, the LASSO regression algorithm was applied to select the optimized feature subset for constructing the final model. A 10-fold cross-validation process was used to determine the optimal hyperparameter λ. The optimal λ values for GTV, GPTV5, GPTV10, and GPTV15 were 0.153, 0.131, 0.068, and 0.100, respectively (see Additional file 8). With these optimal λ values, 2, 5, 7, and 3 features were selected to construct the radiomics models for GTV, GPTV5, GPTV10, and GPTV15, respectively (Fig. 4). The features used for model construction and their ICC details are provided in Additional file 9. The radscore formulas for the four radiomics models can be found in Additional file 10. The radscore for VPI-positive groups in all models was significantly higher than for VPI-negative groups (all P < 0.05), as shown in Additional file 11.

Fig. 4
figure 4

Retained radiomics features and corresponding coefficients of different models after dimensionality reduction by LASSO regression analysis (A) GTV model, (B) GPTV5 model, (C) GPTV10 model, (D) GPTV15 model

Efficacy comparison of radiomics models

The AUC values for the GTV, GPTV5, GPTV10, and GPTV15 models in the training set for predicting VPI status were 0.838, 0.849, 0.855, and 0.841, respectively. In the internal validation set, the corresponding AUC values were 0.808, 0.855, 0.842, and 0.824. Similarly, in the external validation set, the AUC values were 0.809, 0.826, 0.842, and 0.823, respectively. The prediction performance of each radiomics model is summarized in Table 3, and the ROC curves for each radiomics model in the three sets can be found in Additional file 12. The DeLong test indicated that in the training set, the GPTV10 model outperformed the GPTV15 model, and the difference was statistically significant (Z = 2.076, P < 0.05). In the internal validation set, the performance of GPTV5 and GPTV10 models was superior to that of GTV, and the differences were statistically significant (Z = 3.030 and 2.163, both P < 0.05). The radiomics model with the highest AUC value in the external validation set was selected as the best radiomics model. Consequently, a combined model was constructed based on the GPTV10 model’s radscore and CT morphological features.

Table 3 Prediction performance of GTV, GPTV5, GPTV10, and GPTV15 radiomics models in three sets

Efficacy evaluation of combined models

The radscore from the GPTV10 model and CT morphological features included in the clinical model were used as predictive variables to construct a combined model and the corresponding nomogram. The formula for the combined model was as follows:

Nomoscore = (Intercept) * -2.738 + solid component size * 0.076 + pleural indentation sign * 1.169 + the presence of solid component contact pleura * 1.178 + vascular convergence sign * 1.329 + the presence of combined emphysema * 0.650 + GPTV10-radscore * 1.110.

The AUC values of the clinical model were 0.885, 0.814, and 0.838 in the three respective sets. The AUC values of the combined model were 0.894, 0.828, and 0.876 in the three respective sets, as presented in Table 4. The nomogram and examples of its clinical application can be found in Figs. 5 and 6. The ROC curves of the GPTV10-based radiomics model, clinical model, and combined model for predicting VPI in the three sets are shown in Fig. 7.

Table 4 Prediction performance of clinical model, GPTV10-radiomics model and combined model in three sets
Fig. 5
figure 5

The developed nomogram based on GPTV10-Radscore and CT morphological features for predicting VPI status in clinical stage IA LUAD patients

Fig. 6
figure 6

An example of the nomogram in clinical application. This was an axial non-enhanced chest CT image of a 64-year-old male patient with a solid nodule in the upper lobe of the right lung (white arrow) with invasive LUAD and VPI-positive. The factors in the nomogram were analyzed as follows: combined with emphysema = “Yes”, vascular convergence = “Yes”, pleural indentation = “No”, solid component size = 21.1 mm, solid attachment = ”Yes”, GPTV10-Radscore = 0.267, total score was 356 and the probability of VPI-positive was 0.911

Fig. 7
figure 7

ROC curves of GPTV10 radiomics model, clinical model and combined model in three sets (A) training set, (B) internal validation set, (C) external validation set

The DeLong test showed that the combined model outperformed the GPTV10 radiomics model in the training set (Z = 2.987, P < 0.05). In the external validation set, the combined model performed better than the clinical model (Z = 2.348, P < 0.05). The Hosmer-Lemeshow test indicated that the combined model was a good fit in all three sets (all P > 0.05), as shown in Fig. 8. The decision curve analysis (DCA) curves revealed that the combined model achieved a better net benefit in predicting VPI status than the clinical model and the GPTV10 radiomics model, as shown in Fig. 9.

Fig. 8
figure 8

The calibration curves of the combined model in three sets (A) training set, (B) internal validation set, (C) external validation set

Fig. 9
figure 9

DCA for the nomogram, radiomics model, and clinical model in three sets. The y-axis represents the net benefit, and the x-axis represents the threshold probability. (A) training set, (B) internal validation set, (C) external validation set

Discussion

For part-solid or solid nodules in contact with the pleura, selecting the optimal surgical approach remains challenging for thoracic surgeons. In this respect, given that VPI is associated with a high risk of recurrence after sublobectomy [39, 40], it is essential to avoid misdiagnosing lung cancer with VPI. Moreover, in cases where patients have low pulmonary function or high surgical risk, the decision to opt for lobectomy should be made cautiously, especially when the tumor’s location allows for sublobectomy [41]. Therefore, accurately determining the presence of VPI is of great importance. Intraoperative frozen pathology for accurately diagnosing VPI status requires comprehensive histological sampling of the tumor and its adjacent pleura. However, this approach may extend the operation time and lead to unnecessary complications [40]. In this study, we developed a combined prediction model based on the GPTV10 radiomics model and traditional CT signs. In the external validation set, the combined model achieved an AUC value of 0.876, an accuracy of 86.42%, a sensitivity of 72.73%, and a high specificity of 95.83%. Furthermore, we constructed a nomogram that provides a visual representation of the intricate regression equation inherent in the combined model. The nomogram can streamline patient assessment and furnish thoracic surgeons with valuable insights for selecting the optimal surgical strategy.

The clinical model constructed in this study included five CT features: solid component size, solid attachment sign, pleural indentation sign, vascular convergence sign, and the presence of emphysema. The presence of a solid component is often indicative of the more aggressive nature of lung adenocarcinoma, serving as a criterion for T staging in clinical stage IA lung cancer [42]. Larger solid components correlate positively with tumor aggressiveness. Additionally, when a solid component is in contact with the pleura, the risk of visceral pleural invasion increases. As tumor malignancy increases, fibrous hyperplasia occurs within the tumor, leading to more significant traction on adjacent tissues. This traction results in morphological changes in nearby bronchi, vessels, and pleura, which manifest as air bronchogram signs, vascular convergence signs, and pleural indentation signs on CT images [43, 44]. Previous research has confirmed that LUAD accompanied by emphysema exhibits more aggressive characteristics and lower disease-free survival rates compared to cases without emphysema [45, 46] Consistently, our study indicated that LUAD with emphysema is more prone to invade the visceral pleura [47]. The solid component size reflects the invasiveness of the tumor itself, while the other features reflect the relationship between the tumor and the adjacent pleura, vessels, and pulmonary background. The clinical model, which incorporates both tumor and peritumoral CT features, provides a more comprehensive assessment than previous studies [16]. It exhibited promising predictive performance with AUC values of 0.885, 0.814, and 0.838 in the three respective sets.

In this study, precise GTV segmentation formed the basis for automatically expanding three different gradient ranges of the peritumoral region, thereby obtaining the integrated 3D segmentation of GPTV. This exploration marked the first attempt to identify the most efficient radiomics model for predicting VPI status in early LUAD, with a multi-center research providing an evaluation of the model’s generalizability. Emphasizing both the relationship between the tumor and adjacent visceral pleura and the tumor’s inherent aggressiveness, the peritumoral region was not solely utilized as VOI for model construction. The results demonstrated that the predictive efficacy of all GPTV radiomics models was superior to that of GTV, underscoring the peritumoral region’s ability to represent the tumor microenvironment and reflect tumor aggressiveness to a certain extent, as noted in previous studies [17, 23,24,25,26, 48]. In addition, the GPTV10-radiomics model exhibited superior performance compared to other models, likely due to the inclusion of tumors with a DLP ranging from 0 mm to 10 mm, where the 10 mm peritumoral extension range accurately covered various high-order features representing the three-dimensional spatial structure relationship between the tumor and pleura.

Building upon the GPTV10 VOI, seven optimal quantitative radiomics features were selected to indirectly reflect the biological differences between VPI-positive and VPI-negative patients. These features included Large Dependence Low Gray Level Emphasis (LDLGLE), Zone Entropy, Gray Level Variance, Cluster Tendency, and Correlation. LDLGLE, Zone Entropy, and Gray Level Variance are parameters of the Gray Level Size Zone Matrix (GLSZM), while Cluster Tendency and Correlation are parameters of the Gray Level Co-occurrence Matrix (GLCM), all of which are texture feature parameters [49]. It is well-established that GLCM reflects the spatial relationship between pixels, GLDM reflects the gray level relationship between the center pixel and its neighborhood, and GLSZM provides information about the uniformity of each gray level on a three-dimensional image [17]. The higher values of LDLGLE, regional entropy, gray variance, cluster trend, and correlation indicated a more uneven distribution of image texture and more irregular gray changes, signifying higher spatial heterogeneity of the tumor [18], and reflecting its strong invasiveness and increased risk of adjacent VPI.

The clinical model and radiomics model were found to be complementary in this study. Although their diagnostic performance was comparable, their integration exhibited superior diagnostic efficiency. The Delong test results indicated that the combined model outperformed the GPTV10 radiomics model in the training set (P < 0.05) and the clinical model in the external validation set (P < 0.05). The calibration curve demonstrated the model’s predicted probabilities to be in acceptable agreement with the actual probability, while the DCA illustrated that the combined model provided more net benefit than both the clinical model and the radiomics model.

There were several previous studies on predicting VPI status in lung adenocarcinoma by CT-based radiomics features [27,28,29,30,31,32]. Compared with previous studies, our study has some innovative points in patient enrollment and research methods. Firstly, lung tumors with specific VPI-negative on CT presentation such as pGGN or unrelated to the pleura were excluded in our study, therefore, our results were more objective in including cases with potential VPI-positive status to study the predictive performance of radiomics model combined with CT features. Secondly, our study explored the value of radiomics models based on GTV and GPTV with different peritumoral areas for predicting VPI status, no similar research has been reported before. Thirdly, based on the TRIPOD statement, different prediction models were established, and a multi-center dataset was included for internal and external validation of the model to verify the generalization of the models, which was lacking in previous studies.

While this study has notable strengths, it is not devoid of limitations. Firstly, the manual delineation of tumor segmentation was both time-consuming and labor-intensive; however, this method yielded a high level of segmentation accuracy. Moreover, the utilization of Python platform editing code for automatic peritumoral area acquisition, based on semi-automatic and manually precise tumor segmentation, significantly improved segmentation efficiency. The approach included the implementation of a threshold value to eliminate non-lung parenchyma areas, such as soft tissue of the chest wall, ribs, neck, mediastinum, and abdominal regions. This not only preserved the accuracy of artificial tumor delineation but also mitigated subjective errors in the manual delineation of the peritumoral area. Secondly, the sample size in this study was relatively small. Thirdly, the retrospective nature of the study inevitably resulted in selection bias. Additionally, some cases lacked clear PL1 and PL2 pathological grades, precluding the possibility of conducting subgroup analyses for these categories. Further studies are warranted in the future, with a focus on collecting a larger number of cases and more detailed information on pathological grading to address these limitations.

Conclusions

In conclusion, radiomics features-based quantitative analysis provides a noninvasive and accurate diagnostic tool for reflecting the biological behavior of tumors. Combining radiomics features based on segmented GTPV10 with traditional CT morphological signs recognized by radiologists maximizes the diagnostic efficacy for preoperative prediction of VPI in clinical stage IA LUAD, thereby contributing to personalized and accurate treatment strategies.

Data availability

No datasets were generated or analysed during the current study.

Abbreviations

VPI:

Visceral pleural invasion

GTV:

Gross tumor volume

GPTV:

Gross peritumoral tumor volume

AUC:

Area under the curve

LUAD:

Lung adenocarcinoma

pGGN:

Pure ground-glass nodule

VOI:

Volume of interest

DLP:

Distance from the lesion to the pleura

MPR:

Multi-planar reconstruction

CTR:

Consolidation-tumor ratio

ICCs:

Intraclass correlation coefficients

mRMR:

maximal redundancy minimal relevance

LASSO:

The least absolute shrinkage and selection operator

ROC:

Receiver operating characteristic

DCA:

Decision curve analysis

References

  1. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global Cancer statistics 2020: GLOBOCAN estimates of incidence and Mortality Worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021;71(3):209–49. https://doi.org/10.3322/caac.21660.

    Article  CAS  PubMed  Google Scholar 

  2. WHO Classification of Tumours Editorial Board. WHO classification of tumours. Thoracic tumours. 5th ed. Lyon: IARC; 2021.

    Google Scholar 

  3. Goldstraw P, Chansky K, Crowley J, Rami-Porta R, Asamura H, Eberhardt WE, et al. The IASLC Lung Cancer Staging Project: proposals for revision of the TNM Stage groupings in the Forthcoming (Eighth) Edition of the TNM classification for Lung Cancer. J Thorac Oncol. 2016;11(1):39–51. https://doi.org/10.1016/j.jtho.2015.09.009.

    Article  PubMed  Google Scholar 

  4. Zhang T, Zhang JT, Li WF, Lin JT, Liu SY, Yan HH, et al. Visceral pleural invasion in T1 tumors (≤ 3 cm), particularly T1a, in the eighth tumor-node-metastasis classification system for non-small cell lung cancer: a population-based study. J Thorac Dis. 2019;11(7):2754–62. https://doi.org/10.21037/jtd.2019.06.32.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Wo Y, Zhao Y, Qiu T, Li S, Wang Y, Lu T, et al. Impact of visceral pleural invasion on the association of extent of lymphadenectomy and survival in stage I non-small cell lung cancer. Cancer Med. 2019;8(2):669–78. https://doi.org/10.1002/cam4.1990.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Takizawa H, Kondo K, Kawakita N, Tsuboi M, Toba H, Kajiura K, et al. Autofluorescence for the diagnosis of visceral pleural invasion in non-small-cell lung cancer. Eur J Cardiothorac Surg. 2018;53(5):987–92. https://doi.org/10.1093/ejcts/ezx419.

    Article  PubMed  Google Scholar 

  7. Qi LP, Li XT, Yang Y, Chen JF, Wang J, Chen ML et al. Multivariate analysis of Pleural Invasion of Peripheral Non-small Cell Lung Cancer-based computed Tomography features. J Comput Assist Tomogr 2016 Sep-Oct;40(5):757–62. https://doi.org/10.1097/RCT.0000000000000439.

  8. Yang S, Yang L, Teng L, Zhang S, Cui Y, Cao Y, et al. Visceral pleural invasion by pulmonary adenocarcinoma ≤ 3 cm: the pathological correlation with pleural signs on computed tomography. J Thorac Dis. 2018l;10(7):3992–9. https://doi.org/10.21037/jtd.2018.06.125.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Hsu JS, Han IT, Tsai TH, Lin SF, Jaw TS, Liu GC, et al. Pleural tags on CT scans to Predict Visceral Pleural Invasion of Non-small Cell Lung Cancer that does not abut the Pleura. Radiology. 2016;279(2):590–6. https://doi.org/10.1148/radiol.2015151120.

    Article  PubMed  Google Scholar 

  10. Onoda H, Higashi M, Murakami T, Tao H, Yokoyama S, Kunihiro Y, et al. Correlation between pleural tags on CT and visceral pleural invasion of peripheral lung cancer that does not appear touching the pleural surface. Eur Radiol. 2021;31(12):9022–9. https://doi.org/10.1007/s00330-021-07869-y.

    Article  PubMed  Google Scholar 

  11. Heidinger BH, Schwarz-Nemec U, Anderson KR, de Margerie-Mellon C, Monteiro Filho AC, Chen Y, et al. Visceral Pleural Invasion in Pulmonary Adenocarcinoma: differences in CT patterns between solid and subsolid cancers. Radiol Cardiothorac Imaging. 2019;1(3):e190071. https://doi.org/10.1148/ryct.2019190071.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Ahn SY, Park CM, Jeon YK, Kim H, Lee JH, Hwang EJ, et al. Predictive CT features of visceral Pleural Invasion by T1-Sized Peripheral Pulmonary Adenocarcinomas manifesting as Subsolid Nodules. AJR Am J Roentgenol. 2017;209(3):561–6. https://doi.org/10.2214/AJR.16.17280.

    Article  PubMed  Google Scholar 

  13. Zhao Q, Wang JW, Yang L, Xue LY, Lu WW. CT diagnosis of pleural and stromal invasion in malignant subpleural pure ground-glass nodules: an exploratory study. Eur Radiol. 2019;29(1):279–86. https://doi.org/10.1007/s00330-018-5558-0.

    Article  PubMed  Google Scholar 

  14. Shi J, Li F, Yang F, Dong Z, Jiang Y, Nachira D, et al. The combination of computed tomography features and circulating tumor cells increases the surgical prediction of visceral pleural invasion in clinical T1N0M0 lung adenocarcinoma. Transl Lung Cancer Res. 2021;10(11):4266–80. https://doi.org/10.21037/tlcr-21-896.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Kim HJ, Cho JY, Lee YJ, Park JS, Cho YJ, Yoon HI, et al. Clinical significance of Pleural attachment and indentation of Subsolid Nodule Lung Cancer. Cancer Res Treat. 2019;51(4):1540–8. https://doi.org/10.4143/crt.2019.057.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Kim H, Goo JM, Kim YT, Park CM. CT-defined visceral Pleural Invasion in T1 lung adenocarcinoma: lack of relationship to Disease-Free Survival. Radiology. 2019;292(3):741–9. https://doi.org/10.1148/radiol.2019190297.

    Article  PubMed  Google Scholar 

  17. Huang L, Lin W, Xie D, Yu Y, Cao H, Liao G, et al. Development and validation of a preoperative CT-based radiomic nomogram to predict pathology invasiveness in patients with a solitary pulmonary nodule: a machine learning approach, multicenter, diagnostic study. Eur Radiol. 2022;32(3):1983–96. https://doi.org/10.1007/s00330-021-08268-z.

    Article  PubMed  Google Scholar 

  18. Chen Q, Shao J, Xue T, Peng H, Li M, Duan S, et al. Intratumoral and peritumoral radiomics nomograms for the preoperative prediction of lymphovascular invasion and overall survival in non-small cell lung cancer. Eur Radiol. 2023;33(2):947–58. https://doi.org/10.1007/s00330-022-09109-3.

    Article  PubMed  Google Scholar 

  19. Hu Y, Xie C, Yang H, Ho JWK, Wen J, Han L, et al. Assessment of Intratumoral and Peritumoral Computed Tomography Radiomics for Predicting Pathological Complete response to Neoadjuvant Chemoradiation in patients with esophageal squamous cell carcinoma. JAMA Netw Open. 2020;3(9):e2015927. https://doi.org/10.1001/jamanetworkopen.2020.15927.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Lubner MG, Smith AD, Sandrasegaran K, Sahani DV, Pickhardt PJ. CT texture analysis: definitions, applications, Biologic correlates, and challenges. Radiographics. 2017;37(5):1483–503. https://doi.org/10.1148/rg.2017170056.

    Article  PubMed  Google Scholar 

  21. Beig N, Khorrami M, Alilou M, Prasanna P, Braman N, Orooji M, et al. Perinodular and Intranodular Radiomic Features on lung CT images distinguish adenocarcinomas from Granulomas. Radiology. 2019;290(3):783–92. https://doi.org/10.1148/radiol.2018180910.

    Article  PubMed  Google Scholar 

  22. Banat GA, Tretyn A, Pullamsetti SS, Wilhelm J, Weigert A, Olesch C, et al. Immune and Inflammatory Cell Composition of Human Lung Cancer Stroma. PLoS ONE. 2015;28(9):e0139073. https://doi.org/10.1371/journal.pone.0139073.

    Article  CAS  Google Scholar 

  23. Nishino M. Perinodular Radiomic features to assess Nodule Microenvironment: does it help to distinguish malignant versus benign lung. Nodules? Radiol. 2019;290(3):793–5. https://doi.org/10.1148/radiol.2018182619.

    Article  Google Scholar 

  24. Christiansen A, Detmar M. Lymphangiogenesis and cancer. Genes Cancer. 2011;2(12):1146–58. https://doi.org/10.1177/1947601911423028.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Wang X, Zhao X, Li Q, Xia W, Peng Z, Zhang R, et al. Can peritumoral radiomics increase the efficiency of the prediction for lymph node metastasis in clinical stage T1 lung adenocarcinoma on CT? Eur Radiol. 2019;29(11):6049–58. https://doi.org/10.1007/s00330-019-06084-0.

    Article  PubMed  Google Scholar 

  26. Liao G, Huang L, Wu S, Zhang P, Xie D, Yao L, et al. Preoperative CT-based peritumoral and tumoral radiomic features prediction for tumor spread through air spaces in clinical stage I lung adenocarcinoma. Lung Cancer. 2022;163:87–95. https://doi.org/10.1016/j.lungcan.2021.11.017.

    Article  PubMed  Google Scholar 

  27. Wei SH, Zhang JM, Shi B, Gao F, Zhang ZX, Qian LT. The value of CT radiomics features to predict visceral pleural invasion in ≤ 3 cm peripheral type early non-small cell lung cancer. J Xray Sci Technol. 2022;30(6):1115–26. https://doi.org/10.3233/XST-221220.

    Article  PubMed  Google Scholar 

  28. Yuan M, Liu JY, Zhang T, Zhang YD, Li H, Yu TF. Prognostic impact of the findings on thin-section computed Tomography in stage I lung adenocarcinoma with visceral pleural invasion. Sci Rep. 2018;8(1):4743. https://doi.org/10.1038/s41598-018-22853-1.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Zha X, Liu Y, Ping X, Bao J, Wu Q, Hu S, et al. A Nomogram Combined Radiomics and clinical features as imaging biomarkers for prediction of visceral Pleural Invasion in Lung Adenocarcinoma. Front Oncol. 2022;12:876264. https://doi.org/10.3389/fonc.2022.876264.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Cai X, Wang P, Zhou H, Guo H, Yang X, Dai Z, et al. CT-based radiomics nomogram for predicting visceral pleural invasion in peripheral T1-sized solid lung adenocarcinoma. Am J Cancer Res. 2023;13(12):5901–13.

    PubMed  PubMed Central  Google Scholar 

  31. Wang F, Pan X, Zhang T, Zhong Y, Wang C, Li H, et al. Predicting visceral pleural invasion in lung adenocarcinoma presenting as part-solid density utilizing a nomogram model combined with radiomics and clinical features. Thorac Cancer. 2024;15(1):23–34. https://doi.org/10.1111/1759-7714.15151.

    Article  CAS  PubMed  Google Scholar 

  32. Huang S, Xu F, Zhu W, Xie D, Lou K, Huang D, et al. Multi-dimensional radiomics analysis to predict visceral pleural invasion in lung adenocarcinoma of ≤ 3 cm maximum diameter. Clin Radiol. 2023;78(11):e847–55. https://doi.org/10.1016/j.crad.2023.07.014.

    Article  CAS  PubMed  Google Scholar 

  33. Song C, Guo Z, Yu D, Wang Y, Wang Q, Dong Z, et al. A prognostic Nomogram Combining Immune-related gene signature and clinical factors predicts survival in patients with lung adenocarcinoma. Front Oncol. 2020;10:1300. https://doi.org/10.3389/fonc.2020.01300.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Tu W, Sun G, Fan L, Wang Y, Xia Y, Guan Y, et al. Radiomics signature: a potential and incremental predictor for EGFR mutation status in NSCLC patients, comparison with CT morphology. Lung Cancer. 2019;132:28–35. https://doi.org/10.1016/j.lungcan.2019.03.025.

    Article  PubMed  Google Scholar 

  35. Suzuki K, Koike T, Asakawa T, Kusumoto M, Asamura H, Nagai K, et al. A prospective radiological study of thin-section computed tomography to predict pathological noninvasiveness in peripheral clinical IA lung cancer (Japan Clinical Oncology Group 0201). J Thorac Oncol. 2011;6(4):751–6. https://doi.org/10.1097/JTO.0b013e31821038ab.

    Article  PubMed  Google Scholar 

  36. Orlhac F, Lecler A, Savatovski J, Goya-Outi J, Nioche C, Charbonneau F, et al. How can we combat multicenter variability in MR radiomics? Validation of a correction procedure. Eur Radiol. 2021;31(4):2272–80. https://doi.org/10.1007/s00330-020-07284-9.

    Article  PubMed  Google Scholar 

  37. Fortin JP, Cullen N, Sheline YI, Taylor WD, Aselcioglu I, Cook PA, et al. Harmonization of cortical thickness measurements across scanners and sites. NeuroImage. 2018;167:104–20. https://doi.org/10.1016/j.neuroimage.2017.11.024.

    Article  PubMed  Google Scholar 

  38. Beaumont H, Iannessi A, Bertrand AS, Cucchi JM, Lucidarme O. Harmonization of radiomic feature distributions: impact on classification of hepatic tissue in CT imaging. Eur Radiol. 2021;31(8):6059–68. https://doi.org/10.1007/s00330-020-07641-8.

    Article  PubMed  Google Scholar 

  39. Lopez Guerra JL, Gomez DR, Lin SH, Levy LB, Zhuang Y, Komaki R, et al. Risk factors for local and regional recurrence in patients with resected N0-N1 non-small-cell lung cancer, with implications for patient selection for adjuvant radiation therapy. Ann Oncol. 2013;24(1):67–74. https://doi.org/10.1093/annonc/mds274.

    Article  CAS  PubMed  Google Scholar 

  40. Xie H, Su H, Chen D, Xie D, Dai C, Ren Y, et al. Use of Autofluorescence to Intraoperatively diagnose visceral Pleural Invasion from Frozen sections in patients with lung adenocarcinoma 2 cm or less. Am J Clin Pathol. 2019;152(5):608–15. https://doi.org/10.1093/ajcp/aqz081.

    Article  PubMed  Google Scholar 

  41. Ettinger DS, Wood DE, Akerley W, Bazhenova LA, Borghaei H, Camidge DR, et al. Non-small Cell Lung Cancer, Version 6.2015. J Natl Compr Canc Netw. 2015;13(5):515–24. https://doi.org/10.6004/jnccn.2015.0071.

    Article  PubMed  Google Scholar 

  42. Travis WD, Asamura H, Bankier AA, Beasley MB, Detterbeck F, Flieder DB, et al. The IASLC Lung Cancer Staging Project: proposals for Coding T Categories for Subsolid Nodules and Assessment of Tumor Size in Part-Solid tumors in the Forthcoming Eighth Edition of the TNM classification of Lung Cancer. J Thorac Oncol. 2016;11(8):1204–23. https://doi.org/10.1016/j.jtho.2016.03.025.

    Article  PubMed  Google Scholar 

  43. Gruden JF. What is the significance of pleural tags? AJR Am J Roentgenol. 1995;164(2):503–4. https://doi.org/10.2214/ajr.164.2.7840000.

    Article  CAS  PubMed  Google Scholar 

  44. Zhang Y, Qiang JW, Ye JD, Ye XD, Zhang J. High resolution CT in differentiating minimally invasive component in early lung adenocarcinoma. Lung Cancer. 2014;84(3):236–41. https://doi.org/10.1016/j.lungcan.2014.02.008.

    Article  PubMed  Google Scholar 

  45. Lim CG, Shin KM, Lim JK, Kim HJ, Kim WH, Cho SH, et al. Emphysema is associated with the aggressiveness of COPD-related adenocarcinomas. Clin Respir J. 2020;14(4):405–12. https://doi.org/10.1111/crj.13146.

    Article  PubMed  Google Scholar 

  46. Parris BA, O’Farrell HE, Fong KM, Yang IA. Chronic obstructive pulmonary disease (COPD) and lung cancer: common pathways for pathogenesis. J Thorac Dis. 2019;11(Suppl 17):S2155–72. https://doi.org/10.21037/jtd.2019.10.54.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Zuo Z, Li Y, Peng K, Li X, Tan Q, Mo Y, et al. CT texture analysis-based nomogram for the preoperative prediction of visceral pleural invasion in cT1N0M0 lung adenocarcinoma: an external validation cohort study. Clin Radiol. 2022;77(3):e215–21. https://doi.org/10.1016/j.crad.2021.11.008.

    Article  CAS  PubMed  Google Scholar 

  48. Joyce JA, Pollard JW. Microenvironmental regulation of metastasis. Nat Rev Cancer. 2009;9(4):239–52. https://doi.org/10.1038/nrc2618.

    Article  CAS  PubMed  Google Scholar 

  49. van Griethuysen JJM, Fedorov A, Parmar C, Hosny A, Aucoin N, Narayan V, et al. Computational Radiomics System to Decode the Radiographic phenotype. Cancer Res. 2017;77(21):e104–7. https://doi.org/10.1158/0008-5472.CAN-17-0339.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors thank Dr. Zhengwei Zhang of the Department of pathology, Second Affiliated Hospital of Navy Medical University, Shanghai 200003, China for his selfless and valuable assistance in the interpretation of pathological results for this study.

Funding

This research was funded by National Key R&D Program of China (2022YFC2010000, 2022YFC2010002), Key Program of National Natural Science Foundation of China (81930049), National Natural Science Foundation of China (82171926, 82202140, 82271994), Military Logistics Health Care Project (22BJZ07), Shanghai Science and Technology Innovation Action Plan Program (19411951300), Clinical Innovative Project of Shanghai Changzheng Hospital (2020YLCYJ-Y24), Program of Science and Technology Commission of Shanghai Municipality (21DZ2202600).

Author information

Authors and Affiliations

Authors

Contributions

Yun Wang, Deng Lyu and Su Hu conceived and designed the study. Yi Xiao, Li Fan and Shiyuan Liu gave administrative support.Yun Wang, Deng Lyu, Su Hu and Yanqing Ma provided study patients.Shaofeng Duan, Yayuan Geng, Taohu Zhou and Wenting Tu collected data and analyzed them. Yun Wang and Deng Lyu wrote the main manuscript text.All authors reviewed the manuscript.

Corresponding authors

Correspondence to Li Fan or Shiyuan Liu.

Ethics declarations

Ethical approval and consent to participate

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 approved by the Ethics Committee of Changzheng Hospital, Naval Medical University(decision number: CZ-20210528-01), and subjects’ informed consent was exempted.

Consent for publication

N/A as no identifable information (image, name etc.) of the participant is revealed in the submission.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wang, Y., Lyu, D., Hu, S. et al. Nomogram using intratumoral and peritumoral radiomics for the preoperative prediction of visceral pleural invasion in clinical stage IA lung adenocarcinoma. J Cardiothorac Surg 19, 307 (2024). https://doi.org/10.1186/s13019-024-02807-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13019-024-02807-7

Keywords