Skip to main content

ORIGINAL RESEARCH article

Front. Oncol., 10 August 2023
Sec. Thoracic Oncology

Dual-energy CT-based radiomics for predicting invasiveness of lung adenocarcinoma appearing as ground-glass nodules

Yuting Zheng,&#x;Yuting Zheng1,2†Xiaoyu Han,&#x;Xiaoyu Han1,2†Xi Jia,Xi Jia1,2Chengyu DingChengyu Ding3Kailu Zhang,Kailu Zhang1,2Hanting Li,Hanting Li1,2Xuexiang CaoXuexiang Cao4Xiaohui ZhangXiaohui Zhang4Xin Zhang,*Xin Zhang1,2*Heshui Shi,*Heshui Shi1,2*
  • 1Department of Radiology, Union Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, Hubei, China
  • 2Hubei Province Key Laboratory of Molecular Imaging, Wuhan, China
  • 3ShuKun (BeiJing) Technology Co., Ltd., Beijing, China
  • 4Clinical Solution, Philips Healthcare, Shanghai, China

Objectives: To explore the value of radiomics based on Dual-energy CT (DECT) for discriminating preinvasive or MIA from IA appearing as GGNs before surgery.

Methods: The retrospective study included 92 patients with lung adenocarcinoma comprising 30 IA and 62 preinvasive-MIA, which were further divided into a training (n=64) and a test set (n=28). Clinical and radiographic features along with quantitative parameters were recorded. Radiomics features were derived from virtual monoenergetic images (VMI), including 50kev and 150kev images. Intraclass correlation coefficients (ICCs), Pearson’s correlation analysis and least absolute shrinkage and selection operator (LASSO) penalized logistic regression were conducted to eliminate unstable and redundant features. The performance of the models was evaluated by area under the curve (AUC) and the clinical utility was assessed using decision curve analysis (DCA).

Results: The DECT-based radiomics model performed well with an AUC of 0.957 and 0.865 in the training and test set. The clinical-DECT model, comprising sex, age, tumor size, density, smoking, alcohol, effective atomic number, and normalized iodine concentration, had an AUC of 0.929 in the training and 0.719 in the test set. In addition, the radiomics model revealed a higher AUC value and a greater net benefit to patients than the clinical-DECT model.

Conclusion: DECT-based radiomics features were valuable in predicting the invasiveness of GGNs, yielding a better predictive performance than the clinical-DECT model.

Introduction

Lung cancer is the main cause of mortality among malignant tumors. Non-small cell lung cancer represents approximately 85% of all lung cancer cases, and adenocarcinoma is the primary histological type (1, 2). Ground-glass nodules (GGNs) refer to nodules with a slight increase in density on computed tomography (CT) that do not cover the bronchial structures or vascular edges (3). According to solid components, GGNs can be divided into pure GGNs (pGGNs) and part-solid nodules. When a GGN is persistent on HRCT, it usually indicates the existence of lung adenocarcinoma (LUAD) or its precursors, including invasive adenocarcinomas (IA), minimally invasive adenocarcinoma (MIA), adenocarcinoma in situ (AIS) and atypical adenomatous hyperplasia (AAH) (46), which can be difficult to distinguish without any additional investigation or more invasive measures.

Preinvasive-MIA and IA required different surgical approaches, postoperative treatments and had distinct prognoses. Limited wedge resection or segmental resection is typically performed for preinvasive-MIA to preserve maximum functional lung parenchyma, while lobectomy is always performed for IA to reduce tumor recurrence (7, 8). Thus, accurately predicting the invasiveness of lung adenocarcinoma before surgical decision-making is crucial for selecting appropriate surgical approach and improving prognosis.

Although some CT findings, such as maximal tumor diameter and density, have previously been shown to aid in identifying tumor invasion (9, 10), evaluating these characteristics relies on radiologists’ experience. Dual-energy CT (DECT) has emerged as a potential clinical diagnostic tool (1113), enabling low-dose scanning, obtaining diagnostic images and various quantitative parameters unavailable with conventional CT. Recent studies (14, 15) have explored using DECT parameters, including CT60 keV values, virtual HU, and normalized iodine concentration (NIC), to assess the invasiveness of LUAD.

Radiomics can extract numerous quantitative features from medical images and convert them into minable data for subsequent analysis to support decision-making (16). Several studies showed the potential of radiomics to predict LUAD invasiveness based on conventional CT (17, 18). The developed radiomics models performed moderately well in predicting tumor invasion with areas under the curves (AUCs) of 0.77 and 0.79, respectively. However, to our knowledge, no previous radiomics analysis have examined the invasion of LUAD manifesting as GGNs based on DECT. It is of great interest whether DECT-based radiomics can improve the evaluation of tumor invasiveness.

Consequently, the aim of the study was to identify the value of DECT-based radiomics in differentiating preinvasive-MIA from IA characterized as GGNs and compared its predictive value with the clinical-DECT model.

Materials and methods

Patients

From January 2021 to February 2022, we included 220 patients with pulmonary nodules who performed DECT-enhanced scanning before surgery at our institution. 128 participants were removed due to (1): solid nodules (n=41) (2); large lesions over 3cm (n=6) (3); inflammatory lesions (n=36) (4); non-adenocarcinoma (n=29) (5); preoperative adjuvant therapy (n=16). Thus, 92 patients were included with 62 cases of preinvasive-MIA and 30 cases of IA and then were randomly assigned to a training (64 patients) and test cohort (28 patients), with a ratio of 7:3 (Figure 1).

FIGURE 1
www.frontiersin.org

Figure 1 Study flowchart. MIA, minimally invasive adenocarcinoma; IA, invasive adenocarcinomas.

Pathological analysis

All surgical specimens were fixed in formalin and stained with hematoxylin-eosin as routinely prescribed by our hospital. Two pathologists examined the specimens and recorded the pathological subtype of each tumor following the International Association for the Study of Lung Cancer (IASLC)/American Thoracic Society (ATS)/European Respiratory Society (ERS) (19). All GGNs were separated into preinvasive-MIA (AAH, AIS and MIA) and IA groups.

DECT acquisition

Following the standard procedure, the CT scans were taken using a dual-layer spectral CT (IQon; Philips Healthcare). The acquisition parameters were as follows:120kVp, with the automatic regulation of the tube current; 512 × 512 matrix; collimation 64 × 0.625 mm; reconstructed slice thickness and interval 1.5 mm/1.5 mm. A contrast medium (Iohexol, 320 mg/ml) was injected intravenously at a rate of 2–2.5 ml/s. Conventional images were reconstructed using the iDose4 (Philips Healthcare) algorithm, and spectral-based images (SBIs) were reconstructed using a spectral reconstruction algorithm.

Image analysis

Two radiologists (HSS and YTZ, with 31 and 2 years of experience in thoracic radiology, respectively) blind to the clinical and histologic findings analyzed CT images and measured the quantitative parameters in the artery phase (AP) and venous phase (VP). Tumor size and density (part-solid nodules and pGGNs) were chosen to help identify tumor invasion in patients with LUAD.

Spectral CT images were quantitatively analyzed by commercially available tools (IntelliSpace Portal v. 10.1, Philips Healthcare). Reconstructed images were as follows: iodine-based material decomposition images, effective atomic number (Zeff) images, and 101 sets of virtual monochromatic images (VMIs). The rules for measuring GGNs were shown below (20, 21): (I) Select a region-of-interest (ROI) to cover over 70% of the tumor area in the largest slice, avoiding major vessels and bubbles. (II) Keep the ROIs consistent in size and position across different image types and phases. (III) Take the average values of 2 independent measurement results for analysis and comparison. The following imaging parameters were acquired: CT attenuation values of nodules from VMIs (including 40 kev and 100 kev), iodine concentration (IC), and Zeff at AP and VP. Since the body weight, circulation status and body composition could affect IC values, they were standardized to that of the aorta in the same slice to obtain the normalized iodine concentration (NIC): NIC = IC Ca/IC aorta, as previously described (22). Mean IC values were recorded for lesions and aorta ROIs. The slope of the spectral curve [slope(k)] was computed: slope(k) = |CT number (40 keV) − CT number (100 keV) |/60.

DECT images segmentation

Studies (23, 24) have shown that VMIs generated by spectral CT can improve reproducibility in measuring part-solid nodules. The monochromatic reconstructions at 50-55 keV yield the best image quality for lung parenchyma. Thus, we chose to delineate tumors on 50kev monoenergetic images. Two junior radiologists (XH and YZ with 6 and 2 years of experience in thoracic imaging, respectively) delineated the ROI on DECT images semi-automatically layer by layer using 3D-slicer software. The 3D-slicer software was then used to automatically reconstruct the three-dimensional volumes of interest (3D-VOIs) (Figure 2). To evaluate intra-and inter-observer repeatability of feature extraction, images from 30 patients were randomly chosen and independently segmented on the target slice by the same radiologist (YZ) after one week and another radiologist (SHS) with 31 years of experience, respectively. Intra-and inter-observer agreement in the imaging evaluation was assessed by the intraclass correlation coefficient (ICC) (25).

FIGURE 2
www.frontiersin.org

Figure 2 Radiomics workflow.

Radiomics feature extraction

Studies (26, 27) have shown that venous phase images can better reflect tumor microcirculation and result in high clinical usefulness. Low-energy images (40–70 keV) in VMI can substantially increase vascular contrast, while high-energy images(80-190kev) can decrease the metallic artifact to improve the detection of lung nodules (28). As a result, the data in venous phase were transmitted to workstations to generate VMIs with 50 keV and 150 keV energy levels for radiomics investigation. A Radiomics system (ultra scholar, ShuKun (BeiJing) Technology) was used for the extraction of radiomics features, and a third-party Python library called Pyradiomics was used for extracting all the radiomics features in the feature extraction module (29). Before radiomics examination, each 3D CT image was resampled to a spacing of (0.7, 0.7, 1.5) mm using a B-spline curve interpolation algorithm. 1037 3D-radiomics features were extracted for each VOI (types and numbers are displayed in Figure 2, and details are available at pyradiomics.readthedocs.io/en/latest/features.html).

Delta- and mean-radiomics features

For DECT scans with 50 keV and 150 keV energy levels, radiomics features (RFs) were extracted for these two energy levels, respectively. The delta-RFs referred to the relative net change of RFs between two energy levels (30):

Relative Net Change =(Feature150kev  Feature50kev)/Feature50kev

And mean-RFs referred to the mean values of RFs between two energy levels:

Mean RF =(Feature150kev + Feature50kev)/2

Statistical analysis

All statistical analyses were carried out using R (version 4.0.2; http://www.Rproject.org) and the SPSS software (SPSS, version 21, IBM, Chicago, IL, USA). The ‘glmnet’ package was used for LASSO binary logistic regression. The ‘rms’ package was used for multivariate binary logistic regression. The ‘pROC’ program was employed to conduct receiver operating characteristic (ROC) analysis. Categorical variables were represented as frequency (percentage) and Fisher’s exact X2 test was applied to compare categorical variables. We used the Shapiro-Wilk test to check the normality assumption of the continuous data. When the data were normally distributed, they were expressed as mean ± standard deviation and the independent sample t test was applied to compare them. If the data distribution was not normal, they were represented as median (interquartile range, IQR) and the Mann-Whitney U test was employed for the nonparametric analysis. P<0.05 (two-tailed) was regarded as statistically significant.

All radiomics features were normalized to the z-score. To eliminate redundant and unstable features, ICC calculation and Pearson’s correlation analysis were conducted (r>0.8, ICC>0.75). Least absolute shrinkage and selection operator (LASSO) analysis was carried out to identify characteristics for further evaluation (31). The maximum area under the curve (AUC) and five-fold cross-validation served as features filtering criteria. A multivariate logistic regression (LR) algorithm established a classification model based on the selected features. Two models were then constructed: a clinical-DECT model incorporating clinical and DECT quantitative parameters and a radiomics model based on the radiomics features. The accuracy, sensitivity and specificity were calculated for each model. The confounder matrix illustrated the predictive ability. The DeLong test was performed to verify the efficiency of different models for diagnosis. Moreover, we conducted a decision curve analysis (DCA) to evaluate the clinical value of the models by estimating the net benefits under different threshold probabilities.

Results

Clinical characteristics and quantitative DECT results of adenocarcinoma

The study included 62 preinvasive-MIA and 30 IA patients. Table 1 summarizes the clinical characteristics, main CT features (tumor size and tumor density type) and quantitative results (IC, Zeff and slope[k]) of lesions. Significant differences could be found in age, sex, smoking and alcohol history (all p <0.05). Considering CT features, IA was typically larger than preinvasive-MIA (15.5 mm vs. 10 mm, p < 0.001). Tumor density differed between the IA and preinvasive-MIA group (p < 0.001). Most IA manifested as part-solid nodules (25/30, 83.3%; Figure 3). Regarding quantitative parameters, there was significant difference in Zeff at AP and VP (p = 0.008, p = 0.003, respectively). NIC was larger in IA than in preinvasive-MIA at VP (p = 0.04). However, no significant difference was found in slope (k) values between the groups.

TABLE 1
www.frontiersin.org

Table 1 Comparison of clinical characteristics and DECT parameters of GGNs between preinvasive-MIA and IA.

FIGURE 3
www.frontiersin.org

Figure 3 (A, B): Invasive adenocarcinoma (IA) in a 44-year-old man with a part-solid nodule (A). Photomicrograph (H&E staining, ×100) confirmed IA (B). (C, D): Minimal invasive adenocarcinoma (MIA) in a 66-year-old woman with a part-solid nodule (C). Photomicrograph (H&E staining, ×100) confirmed MIA (D).

Afterwards, the clinical-DECT model was established by LR algorithm (Table 2), including sex, age, smoking, alcohol history, tumor size, density, Zeff, and NIC, with an AUC of 0.929 (sensitivity of 90.5% and specificity of 69.8%) in the training cohort and an AUC of 0.719 (sensitivity of 55.6% and specificity of 73.7%) in the test cohort (Figure 4A). Table 3 illustrated the distribution of selected clinical and DECT characteristics in the training and test set.

TABLE 2
www.frontiersin.org

Table 2 Features included in clinical-DECT model and their coefficients.

FIGURE 4
www.frontiersin.org

Figure 4 Receiver operating characteristic curves of the clinical-DECT model (A) and the radiomics model (B) in the training and test sets.

TABLE 3
www.frontiersin.org

Table 3 Comparison of selected clinical characteristics and DECT parameters of GGNs between preinvasive-MIA and IA in the training and test cohorts.

Radiomics model building and validation

According to Pearson’s correlation analysis results, 332 delta-radiomics features were selected to build the delta-radiomics model. Then, 138 features were preserved from 50kev images by ICC and PCC. Based on the mean-RFs mentioned above, 276 features were derived combined with the features of 150kev images. According to LASSO penalized logistic regression analysis, seven features (three first order and four second order parameters, including GLCM and GLRLM features) revealed a significant correlation between radiomics and adenocarcinoma invasiveness (Figures 5A, B). Then, these radiomics features were taken into the LR model to establish a radiomics model. Table 4 illustrated the features included in the radiomics model and their coefficients. The radiomics model performed well in the training (AUC: 0.957, sensitivity: 95.2%, and specificity: 79.1%) and test cohorts (AUC:0.865; sensitivity:88.9%, and specificity: 73.7%) (Figure 4B).

FIGURE 5
www.frontiersin.org

Figure 5 Radiomics feature selection using least absolute shrinkage and selection operator (LASSO) logistic regression. LASSO logistic regression of radiomics features (A) and the AUC versus the regularization parameter lambda (B).

TABLE 4
www.frontiersin.org

Table 4 Features included in radiomics model and their coefficients.

The rad-score of each lesion was computed according to the following formula:

Rad-Score = -1.54+1.99*delta_wavelet.LLH_firstorder_Median   +0.64*mean_wavelet.LHL_glrlm_LongRunEmphasis   +0.94*mean_wavelet.HLL_glrlm_LongRunHighGrayLevelEmphasis   +1.42*mean_wavelet.LLL_firstorder_RobustMeanAbsoluteDeviation    -1.65*mean_log.sigma.6.0.mm.3D_firstorder_10Percentile   -0.87*mean_log.sigma.6.0.mm.3D_firstorder_InterquartileRange   +0.10*mean_log.sigma.6.0.mm.3D_glcm_ClusterShade

Figure S1 displayed the rad scores significantly higher in IA than in preinvasive-MIA (p <0.001) in the training and test cohorts.

Performance comparison between different models

In comparison with the clinical-DECT model, the radiomics model exhibited improved predictive ability in both the training and test cohorts (AUC, 0.957 vs 0.929, p = 0.341; 0.865 vs 0.719, p = 0.09, respectively). Table 5 summarized the findings of confounder matrix analysis. The radiomics model exhibited better accuracy, sensitivity, and specificity than the clinical-DECT model in the training and test cohorts.

TABLE 5
www.frontiersin.org

Table 5 Confounder matrix for the training and test sets in the two models.

Clinical use

DCA assessed the clinical utility of the two predictive models (Figure S2). Compared to the treat-all and treat-none models, both the clinical-DECT and radiomics models provided a net benefit to patients, with the radiomics model showing superior benefits. Figure S3 showed waterfall plots of the rad-score for the test and training sets, indicating that the model can effectively distinguish IA from preinvasive-MIA.

Discussion

The aim of the study was to investigate the predictive value of DECT-based radiomics for identifying the invasiveness of GGNs. We evaluated the main radiographic features, measured quantitative parameters, and extracted radiomics features from VMI (50kev and 150kev images). Our findings confirmed that the DECT-based radiomics model was superior to the clinical-DECT model in distinguishing IA from preinvasive-MIA, providing a simple and reliable prediction method for accurately determining tumor invasiveness before surgery.

As expected, our study found that IA tended to be larger and part-solid nodules on CT images, similar to the findings of previous studies (9, 32). For quantitative parameters from DECT in our study, Zeff and NIC were considered significant differentiators of IA from preinvasive-MIA. Zeff represents the composite atom, which can be used to identify material composition (12). Li et al. (33) found that tumors with a predominant pattern of solid/micropapillary had lower Zeff than those with a predominant pattern of lepidic/acinar/papillary. Consistent with the previous study, the Zeff of IA was smaller than that of preinvasive-MIA in our research. This finding may be because preinvasive-MIA tends to have lepidic growth patterns, whereas IA can have the more aggressive adenocarcinoma subtype (solid and micropapillary patterns) (34). The result confirmed the value of Zeff in differentiating IA from preinvasive-MIA. In contrast to the finding of Zhang et al. (35), the NIC in IA was higher than in preinvasive-MIA in our study, suggesting that IA had more underlying microvascular and tumor angiogenesis. The observation can be explained as we included the part-solid nodules, and carefully avoided vessels when outlining the ROIs. In addition, we observed no statistical difference in slope(k) at both VP and AP, consistent with Zhang et al. (35). It indicated that different histological subtypes of adenocarcinomas manifested as GGNs had similar changes in iodine concentration.

Radiomics analysis could use the information from VMI datasets to expand new horizons for noninvasively assessing tumor signatures. One study (36) found that the radiomics model based on 70keV images showed good performance in predicting head and neck squamous carcinoma differentiation. Another study (37) demonstrated that radiomics derived from VMIs have high diagnostic efficiencies for differentiating low- and high-grade renal cell carcinoma. In our study, the radiomics model based on 50 keV and 150 kev images could help predict tumor invasiveness with an accuracy of 78.6% in the test set. Among the features included in the radiomics model, the firstorder_Median (media gray level intensity within the tumor) was closely associated with LUAD invasiveness with the highest estimate coefficient (1.99). We assumed that the difference in heterogeneity between the preinvasive-MIA and IA may be associated with their different density. Our results suggested that 50kev images combined with 150kev images can provide helpful information for predicting tumor invasiveness.

Several studies (38, 39) have analyzed radiomics for predicting the invasiveness of GGNs in LUAD. Zhang et al. (38) enrolled 65 patients with pGGNs and found that the radiomics model performed well in identifying invasive lung lesions with an AUC of 0.82. Weng et al. (39) enrolled 119 patients with part-solid nodules and found that the radiomics model had superior predictive performance compared to the radiographic model (AUC, 0.81 vs. 0.76). Previous studies mainly explored radiographic features and radiomics based on traditional CT. The current study, for the first time, used DECT-based radiomics to establish a predictive model (AUC, 0.96 and 0.87 in the training and test set, respectively) of LUAD invasiveness manifesting as GGNs. In addition to the different imaging techniques, the difference among the three studies might also be attributable to inclusion and grouping standards of the patients. Weng et al. included part-solid nodules, while Zhang et al. analyzed pure GGNs. In addition, Zhang et al. divided the nodules into preinvasive lesions (AAH and AIS) and invasive lesions (MIA and IA), whereas Weng et al. classified nodules into MIA and IA. Thus, larger samples and better designs are required in future studies to validate the present results. Despite the differences among various studies, these results still indicate that models from DECT can well estimate the invasiveness of LUAD. Furthermore, our radiomics model did not significantly outperform the clinical-DECT model in predictive ability. This may be because some of the imaging features we evaluated were probably contained in the radiomics features.

Study limitations

There were several limitations to this study. First, this was a single-center retrospective study and the sample size was relatively small, which may constrain the generalizability of our findings. Second, we only included the CT features that have been reported and validated. The results may not reflect the complete CT morphological features of tumors. Third, the main purpose of the study was to identify tumor invasiveness of GGNs, and the density of nodules was not further subdivided. Finally, as the patients were examined from 2021 to 2022 with a limited follow-up time, the influence of IA and preinvasive-MIA on patient outcomes was not assessed.

Conclusions

In summary, the DECT-based radiomics model showed satisfactory predictive performance in preoperatively differentiating IA and preinvasive-MIA. Further prospective multicenter studies are necessary to assess the utility of the model for clinical application.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.

Ethics statement

The studies involving human participants were reviewed and approved by Ethics Committee of Wuhan Union Hospital. The ethics committee waived the requirement of written informed consent for participation.

Author contributions

Conception and design: YZ, XH, XZ, HS. Administrative support: HS, XZ. Provision of study materials or patients: XJ, XC, XHZ. Collection and assembly of data: YZ, XH, KZ, HL. Data analysis and interpretation: CD, XC, XHZ. Manuscript writing: all authors. Final approval of manuscript: all authors.

Conflict of interest

Author CD was employed by the company ShuKun BeiJing Technology Co., Ltd. XHZ and XC were employed by Philips Healthcare.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2023.1208758/full#supplementary-material

Supplementary Figure 1 | The Rad-score between groups in the training set (A) and test set (B).

Supplementary Figure 2 | Decision curves for the two models.

Supplementary Figure 3 | Waterfall plot of the rad-score showing the predicted probabilities of the training (A) and test(B) sets.

References

1. Siegel RL, Miller KD, Jemal A. Cancer statistics, 2020. CA Cancer J Clin (2020) 70:7–30. doi: 10.3322/caac.21590

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Ma J, Ward EM, Smith R, Jemal A. Annual number of lung cancer deaths potentially avertable by screening in the United States. Cancer (2013) 119:1381–5. doi: 10.1002/cncr.27813

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Austin JH, Müller NL, Friedman PJ, Hansell DM, Naidich DP, Remy-Jardin M, et al. Glossary of terms for CT of the lungs: recommendations of the Nomenclature Committee of the Fleischner Society. Radiology (1996) 200:327–31. doi: 10.1148/radiology.200.2.8685321

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Lim HJ, Ahn S, Lee KS, Han J, Shim YM, Woo S, et al. Persistent pure ground-glass opacity lung nodules ≥ 10 mm in diameter at CT scan: histopathologic comparisons and prognostic implications. Chest (2013) 144:1291–9. doi: 10.1378/chest.12-2987

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Son JY, Lee HY, Kim JH, Han J, Jeong JY, Lee KS, et al. Quantitative CT analysis of pulmonary ground-glass opacity nodules for distinguishing invasive adenocarcinoma from non-invasive or minimally invasive adenocarcinoma: the added value of using iodine mapping. Eur Radiol (2016) 26:43–54. doi: 10.1007/s00330-015-3816-y

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Son JY, Lee HY, Lee KS, Kim JH, Han J, Jeong JY, et al. Quantitative CT analysis of pulmonary ground-glass opacity nodules for the distinction of invasive adenocarcinoma from pre-invasive or minimally invasive adenocarcinoma. PloS One (2014) 9:e104066. doi: 10.1371/journal.pone.0104066

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Zhang Y, Ma X, Shen X, Wang S, Li Y, Hu H, et al. Surgery for pre- and minimally invasive lung adenocarcinoma. J Thorac Cardiovasc Surg (2022) 163:456–64. doi: 10.1016/j.jtcvs.2020.11.151

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Liu S, Wang R, Zhang Y, Li Y, Cheng C, Pan Y, et al. Precise diagnosis of intraoperative frozen section is an effective method to guide resection strategy for peripheral small-sized lung adenocarcinoma. J Clin Oncol (2016) 34:307–13. doi: 10.1200/jco.2015.63.4907

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Ohde Y, Nagai K, Yoshida J, Nishimura M, Takahashi K, Suzuki K, et al. The proportion of consolidation to ground-glass opacity on high resolution CT is a good predictor for distinguishing the population of non-invasive peripheral adenocarcinoma. Lung Cancer (2003) 42:303–10. doi: 10.1016/j.lungcan.2003.07.001

PubMed Abstract | CrossRef Full Text | Google Scholar

10. 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:236–41. doi: 10.1016/j.lungcan.2014.02.008

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Li M, Zhang L, Tang W, Jin YJ, Qi LL, Wu N. Identification of epidermal growth factor receptor mutations in pulmonary adenocarcinoma using dual-energy spectral computed tomography. Eur Radiol (2019) 29:2989–97. doi: 10.1007/s00330-018-5756-9

PubMed Abstract | CrossRef Full Text | Google Scholar

12. González-Pérez V, Arana E, Barrios M, Bartrés A, Cruz J, Montero R, et al. Differentiation of benign and Malignant lung lesions: Dual-Energy Computed Tomography findings. Eur J Radiol (2016) 85:1765–72. doi: 10.1016/j.ejrad.2016.07.019

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Deniffel D, Sauter A, Fingerle A, Rummeny EJ, Makowski MR, Pfeiffer D. Improved differentiation between primary lung cancer and pulmonary metastasis by combining dual-energy CT-derived biomarkers with conventional CT attenuation. Eur Radiol (2021) 31:1002–10. doi: 10.1007/s00330-020-07195-9

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Yang Y, Li K, Sun D, Yu J, Cai Z, Cao Y, et al. Invasive pulmonary adenocarcinomas versus preinvasive lesions appearing as pure ground-glass nodules: differentiation using enhanced dual-source dual-energy CT. AJR Am J Roentgenol (2019) 213:W114–22. doi: 10.2214/AJR.19.21245

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Wang S, Liu G, Fu Z, Jiang Z, Qiu J. Predicting pathological invasiveness of lung adenocarcinoma manifesting as GGO-predominant nodules: a combined prediction model generated from DECT. Acad Radiol (2021) 28:509–16. doi: 10.1016/j.acra.2020.03.007

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Gillies RJ, Kinahan PE, Hricak H. Radiomics: images are more than pictures, they are data. Radiology (2016) 278:563–77. doi: 10.1148/radiol.2015151169

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Sun Y, Li C, Jin L, Gao P, Zhao W, Ma W, et al. Radiomics for lung adenocarcinoma manifesting as pure ground-glass nodules: invasive prediction. Eur Radiol (2020) 30:3650–9. doi: 10.1007/s00330-020-06776-y

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Xue X, Yang Y, Huang Q, Cui F, Lian Y, Zhang S, et al. Use of a radiomics model to predict tumor invasiveness of pulmonary adenocarcinomas appearing as pulmonary ground-glass nodules. BioMed Res Int (2018) 2018:6803971. doi: 10.1155/2018/6803971

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Zugazagoitia J, Enguita AB, Nuñez JA, Iglesias L, Ponce S. The new IASLC/ATS/ERS lung adenocarcinoma classification from a clinical perspective: current concepts and future prospects. J Thorac Dis (2014) 6:S526–36. doi: 10.3978/j.issn.2072-1439.2014.01.27

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Henschke CI, Yankelevitz DF, Mirtcheva R, McGuinness G, McCauley D, Miettinen OS, et al. CT screening for lung cancer: frequency and significance of part-solid and nonsolid nodules. AJR Am J Roentgenol (2002) 178:1053–7. doi: 10.2214/ajr.178.5.1781053

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Kim HY, Shim YM, Lee KS, Han J, Yi CA, Kim YK. Persistent pulmonary nodular ground-glass opacity at thin-section CT: histopathologic comparisons. Radiology (2007) 245:267–75. doi: 10.1148/radiol.2451061682

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Deniffel D, Sauter A, Dangelmaier J, Fingerle A, Rummeny EJ, Pfeiffer D. Differentiating intrapulmonary metastases from different primary tumors via quantitative dual-energy CT based iodine concentration and conventional CT attenuation. Eur J Radiol (2019) 111:6–13. doi: 10.1016/j.ejrad.2018.12.015

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Ohana M, Labani A, Severac F, Jeung MY, Gaertner S, Caspar T, et al. Single source dual energy CT: What is the optimal monochromatic energy level for the analysis of the lung parenchyma? Eur J Radiol (2017) 88:163–70. doi: 10.1016/j.ejrad.2017.01.016

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Kim J, Lee KH, Kim J, Shin YJ, Lee KW. Improved repeatability of subsolid nodule measurement in low-dose lung screening with monoenergetic images: a phantom study. Quant Imaging Med Surg (2019) 9:171–9. doi: 10.21037/qims.2018.10.06

CrossRef Full Text | Google Scholar

25. Zou GY. Sample size formulas for estimating intraclass correlation coefficients with precision and assurance. Stat Med (2012) 31:3972–81. doi: 10.1002/sim.5466

CrossRef Full Text | Google Scholar

26. Zegadło A, Żabicka M, Kania-Pudło M, Maliborski A, Różyk A, Sośnicki W. Assessment of solitary pulmonary nodules based on virtual monochrome images and iodine-dependent images using a single-source dual-energy CT with fast kVp switching. J Clin Med (2020) 9:2514. doi: 10.3390/jcm9082514

CrossRef Full Text | Google Scholar

27. Zhang Z, Zou H, Yuan A, Jiang F, Zhao B, Liu Y, et al. A single enhanced dual-energy CT scan may distinguish lung squamous cell carcinoma from adenocarcinoma during the venous phase. Acad Radiol (2020) 27:624–9. doi: 10.1016/j.acra.2019.07.018

CrossRef Full Text | Google Scholar

28. D’Angelo T, Cicero G, Mazziotti S, Ascenti G, Albrecht MH, Martin SS, et al. Dual energy computed tomography virtual monoenergetic imaging: technique and clinical applications. Br J Radiol (2019) 92:20180546. doi: 10.1259/bjr.20180546

CrossRef Full Text | Google Scholar

29. 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:e104–7. doi: 10.1158/0008-5472.CAN-17-0339

CrossRef Full Text | Google Scholar

30. Fave X, Zhang L, Yang J, Mackin D, Balter P, Gomez D, et al. Delta-radiomics features for the prediction of patient outcomes in non-small cell lung cancer. Sci Rep (2017) 7:588. doi: 10.1038/s41598-017-00665-z

CrossRef Full Text | Google Scholar

31. Huang Y, Liu Z, He L, Chen X, Pan D, Ma Z, et al. Radiomics signature: A potential biomarker for the prediction of disease-free survival in early-stage (I or II) non-small cell lung cancer. Radiology (2016) 281:947–57. doi: 10.1148/radiol.2016152234

CrossRef Full Text | Google Scholar

32. Yu Y, Cheng JJ, Li JY, Zhang Y, Lin LY, Zhang F, et al. Determining the invasiveness of pure ground-glass nodules using dual-energy spectral computed tomography. Transl Lung Cancer Res (2020) 9:484–95. doi: 10.21037/tlcr.2020.03.33

CrossRef Full Text | Google Scholar

33. Li Q, Li X, Li XY, He XQ, Chu ZG, Luo TY. Histological subtypes of solid-dominant invasive lung adenocarcinoma: differentiation using dual-energy spectral CT. Clin Radiol (2021) 76:77. doi: 10.1016/j.crad.2020.08.034

CrossRef Full Text | Google Scholar

34. Weichert W, Warth A. Early lung cancer with lepidic pattern: adenocarcinoma in situ, minimally invasive adenocarcinoma, and lepidic predominant adenocarcinoma. Curr Opin Pulm Med (2014) 20:309–16. doi: 10.1097/MCP.0000000000000065

CrossRef Full Text | Google Scholar

35. Zhang Y, Tang J, Xu J, Cheng J, Wu H. Analysis of pulmonary pure ground-glass nodule in enhanced dual energy CT imaging for predicting invasive adenocarcinoma: comparing with conventional thin-section CT imaging. J Thorac Dis (2017) 9:4967–78. doi: 10.21037/jtd.2017.11.04

CrossRef Full Text | Google Scholar

36. Li Z, Liu Z, Guo Y, Wang S, Qu X, Li Y, et al. Dual-energy CT-based radiomics nomogram in predicting histological differentiation of head and neck squamous carcinoma: a multicenter study. Neuroradiology (2022) 64:361–9. doi: 10.1007/s00234-021-02860-2

CrossRef Full Text | Google Scholar

37. Han D, Yu Y, He T, Yu N, Dang S, Wu H, et al. Effect of radiomics from different virtual monochromatic images in dual-energy spectral CT on the WHO/ISUP classification of clear cell renal cell carcinoma. Clin Radiol (2021) 76:627.e23–627.e29. doi: 10.1016/j.crad.2021.02.033

CrossRef Full Text | Google Scholar

38. Zhang T, Li X, Liu J. Prediction of the invasiveness of ground-glass nodules in lung adenocarcinoma by radiomics analysis using high-resolution computed tomography imaging. Cancer Control (2022) 29:10732748221089408. doi: 10.1177/10732748221089408

CrossRef Full Text | Google Scholar

39. Weng Q, Zhou L, Wang H, Hui J, Chen M, Pang P, et al. A radiomics model for determining the invasiveness of solitary pulmonary nodules that manifest as part-solid nodules. Clin Radiol (2019) 74:933–43. doi: 10.1016/j.crad.2019.07.026

CrossRef Full Text | Google Scholar

Keywords: dual-energy spectral computed tomography, radiomics, prediction, lung adenocarcinoma, ground-glass nodules, invasion

Citation: Zheng Y, Han X, Jia X, Ding C, Zhang K, Li H, Cao X, Zhang X, Zhang X and Shi H (2023) Dual-energy CT-based radiomics for predicting invasiveness of lung adenocarcinoma appearing as ground-glass nodules. Front. Oncol. 13:1208758. doi: 10.3389/fonc.2023.1208758

Received: 19 April 2023; Accepted: 24 July 2023;
Published: 10 August 2023.

Edited by:

Zhongxiang Ding, Zhejiang University, China

Reviewed by:

Xiaohu Li, First Affiliated Hospital of Anhui Medical University, China
Asli Suner, Ege University, Türkiye

Copyright © 2023 Zheng, Han, Jia, Ding, Zhang, Li, Cao, Zhang, Zhang and Shi. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Xin Zhang, 286938051@qq.com; Heshui Shi, heshuishi@hust.edu.cn

These authors have contributed equally to this work

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.