Impact Factor 4.848 | CiteScore 3.5
More on impact ›

Original Research ARTICLE

Front. Oncol., 20 October 2020 | https://doi.org/10.3389/fonc.2020.580809

The Growth Trend Predictions in Pulmonary Ground Glass Nodules Based on Radiomic CT Features

Chen Gao1,2, Jing Yan1,2, Yifan Luo1,2, Linyu Wu1,2, Peipei Pang3, Ping Xiang1,2 and Maosheng Xu1,2*
  • 1Department of Radiology, The First Affiliated Hospital of Zhejiang Chinese Medical University (Zhejiang Provincial Hospital of Traditional Chinese Medicine), Hangzhou, China
  • 2The First Clinical Medical College of Zhejiang Chinese Medical University, Hangzhou, China
  • 3GE Healthcare Life Sciences, Hangzhou, China

Background: The management of ground glass nodules (GGNs) remains a distinctive challenge. This study is aimed at comparing the predictive growth trends of radiomic features against current clinical features for the evaluation of GGNs.

Methods: A total of 110 GGNs in 85 patients were included in this retrospective study, in which follow up occurred over a span ≥2 years. A total of 396 radiomic features were manually segmented by radiologists and quantitatively analyzed using an Analysis Kit software. After feature selection, three models were developed to predict the growth of GGNs. The performance of all three models was evaluated by a receiver operating characteristic (ROC) curve. The best performing model was also assessed by calibration and clinical utility.

Results: After using a stepwise multivariate logistic regression analysis and dimensionality reduction, the diameter and five specific radiomic features were included in the clinical model and the radiomic model. The rad-score [odds ratio (OR) = 5.130; P < 0.01] and diameter (OR = 1.087; P < 0.05) were both considered as predictive indicators for the growth of GGNs. Meanwhile, the area under the ROC curve of the combined model reached 0.801. The high degree of fitting and favorable clinical utility was detected using the calibration curve with the Hosmer-Lemeshow test and the decision curve analysis was utilized for the nomogram.

Conclusions: A combined model using the current clinical features alongside the radiomic features can serve as a powerful tool to assist clinicians in guiding the management of GGNs.

Introduction

The detection rate of pulmonary nodules has been significantly increased since the introduction of low dose CT screening, especially for the ground glass nodule (GGN) (1, 2). The GGN, which includes pure and part-solid GGN, is defined as a hazy region of increased opacity on lung windows without obscurity to bronchial and vascular structures (3). The pathophysiology of GGN is based on the accumulation of fluid, cells or amorphous material in the alveoli itself, or thickening of the alveolar walls and septal interstitium (4). The GGN is observed in many lesions, such as malignant tumors and benign lesions which include inflammatory lesions, interstitial lung disease, and so on (57). In 2017, the Fleischner society released new guidelines for GGN management, with more aggressive guidelines toward follow up (8). Although radiologists were able to observe changes of GGN in follow-up CT examinations, most GGNs progress at a slow rate, particularly persistent GGN (9). Therefore, long windows of follow-up are often required. This is a source of great anxiety for patients and their families. Furthermore, the increased duration of follow up often increases the rate of no-shows. Therefore, several studies have sought to provide a greater diagnostic indicator for the growth of GGNs through the analysis of traditional imaging features, such as diameter and CT attenuation (1013). For example, Matsuguma et al. showed there were significant differences in diameter between rapidly growing and non-growing pure GGNs (12). However, distinguishing the growing GGNs from static GGNs using traditional quantitative CT imaging remains a distinctive challenge.

With the advancements in imaging technology, many radiologists attach importance to these parameters, such as the Gray-Level Co-occurrence Matrix (GLCM) and the Gray Level Run-length Matrix (RLM) (14). While there have been many radiomic pulmonary studies in recent years, there have been no studies comparing and contrasting radiomic features with clinical features to predict the growth of GGNs (15, 16).

Therefore, the purpose of the current study is to compare the performance of clinical signatures and radiomic features in predicting the growth of GGNs and to build a clinical-radiomic nomogram to accurately predict the growth of GGNs.

Materials and Methods

Patients

This retrospective study was approved by our institutional review board and informed consent requirement was waived.

Between September 2012 and December 2018, a total of 85 patients with 110 GGNs were involved in this study which included 68 patients with a single GGN, 11 patients with two GGNs, 5 patients with three GGNs, and 1 patient with five GGNs. The inclusion criteria were as follows: (1) detected pulmonary nodules showed GGNs on non-enhanced CT thin-sectioned images; (2) the GGN diameter between 5 and 30 mm in the initial CT image; (3) there were more than two follow-up, thin-section CT examinations and the follow-up interval was longer than 2 years.

The exclusion criteria were as follows: (1) biopsy, radiotherapy, chemotherapy or surgical resection during any follow-up; (2) severe respiratory artifacts on CT images; (3) a history of lung surgery; (4) the first or final CT examinations were low dose CT examinations.

The patients were randomly divided into a training cohort and a validation cohort. A flow chart of patients who were selected is presented in Figure 1.

FIGURE 1
www.frontiersin.org

Figure 1. Flow diagram of GGN selection.

CT Examination Acquisition

All images were obtained with a Siemens Somatom sensation, 64 slice, CT scanner (Siemens Healthcare). The imaging parameters were as follows: tube voltage, 120 kVp; automatic tube current modulation; collimation, 0.6 mm*64; matrix, 512*512; rotation time, 0.37 s, reconstruction slice thickness, 0.75 mm with a 0.5 mm interval, reconstruction kernel, B31f.

Nodule Selection and Growth Definition

The GGNs were assessed by two radiologists (one with 11 years of experience and another with 4 years of experience in pulmonary radiology) based on thin-section unenhanced CT images. All GGNs were confirmed by two radiologists with consensus agreement. If the two radiologists cannot reach a consensus, the GGNs were assessed by a third-party professor with 28 years of experience in pulmonary radiology. The diameter of GGN was defined as the maximum length of on the transverse lung window in thin-section CT images. The solid diameter was defined as the maximum size of the solid portion of GGNs on the transverse lung window in thin-section CT images. All measurements in the initial and final CT images were constructed from transverse sections by two radiologists to reach a consensus. To eliminate measuring error, growth was defined as an increase in diameter or the size of the solid portion ≥2 mm, or an emerging solid portion (17).

Region of Interest Segmentation and Feature Extraction

All regions of interest (ROI) were manually segmentation by a radiologist with 4 years of experience in pulmonary radiology on initial thin-section CT images by using ITK-SNAP 3.6.0 (www.itksnap.org), and further verified by another radiologist with 11 years of experience in pulmonary radiology. For situations of a discrepancy between the two radiologists, the segmentation patterns were evaluated by a professor of radiology with 28 years of experience in pulmonary radiology. Normal structures within or around the GGNs, such as vessels and pleura, were not included in ROIs. A total of 396 radiomic features were quantitatively extracted using Analysis Kit software (AK, GE Healthcare). These features included single-order (histograms and morphologic features) alongside higher-order parameters (Supplementary Descriptions). The higher-order parameters were described as “texture” features, such as GLCM and RLM. Texture features described statistical interrelationships between voxels with similar (or dissimilar) contrast values. The values of each feature for all GGNs were normalized with Z-scores ((x – μ)/σ) for the purpose of removing the unit limits of each feature before being applied to a machine learning model for classification. For the model parameters, x represents the value of the feature, μ indicates the average feature values in all GGNs within the cohort, and σ represents the corresponding standard deviation.

Radiomic Feature Selection and Construction of Rad-Score

The minimum Redundancy Maximum Relevance (mRMR) and the least absolute shrinkage and selection operator (LASSO) algorithm were used to select radiomic features (18). Initially, mRMR was performed to eliminate redundant and irrelevant features. The top 20 features were selected, and the optimized subset of features was chosen by LASSO to construct the final model. After the number of features was determined, the most predictive feature subset was chosen and the corresponding coefficients were calculated. Rad-score was composed by summing the chosen features, weighted by their coefficients, and comparing it between the training group and test group. The performance of the model was then evaluated by ROC analysis.

Establishment of the Clinical and Combined Model

In the establishment of the clinical model, a univariate logistic regression was used to evaluate the clinicopathological factors. Factors with a P-value < 0.05 were considered in the stepwise multivariate logistic regression analysis. Meanwhile, the Akaike information criterion (AIC) was used to determine a stopping rule. A combination of the clinical signatures from the clinical model and rad-score were used to develop the combined model with multivariate logistic regression. Afterward a test process was implemented.

Model Comparison and Nomogram Establishment

The predictive accuracy of the three models was assessed by the area under the curve (AUC) in both the training and validation group. The AIC of the clinical model was applied to identify the most appropriate clinical model. The probability of growth for each GGN was analyzed by logistic regression, and GGNs were grouped into growth and non-growth cohorts based on the highest Youden index. According to the actual growth results, the accuracy, specificity, sensitivity, negative-predictive value (NPV) and positive-predictive value (PPV) for the three models were calculated in the training and test group. Then, the nomogram of the most appropriate model was established. According to the reference of Iasonos et al. and Stephenson et al., the usefulness of a nomogram is that it maps the predicted probabilities into points on a scale from 0 to 100 in a picture and the total points accumulated by the various factors correspond to the predicted probability for a patient (19, 20). Meanwhile, the calibration curves measured the consistency between the actual growth probability and the predicted growth probability to evaluate the running characteristic of the nomogram. The degree of fit of the prediction models was also evaluated by the Hosmer-Lemeshow (H-L) test.

Development of Decision Curve Analysis (DCA)

To assess the added value of radiomic features to clinical in the prediction of the growth of GGNs, three DCA was performed based on clinical diagnosis, radiomics, and the combined model. The clinical application of said model could be verified by quantifying the net benefits for a range of threshold probabilities.

Statistical Analysis

R statistical software (version 3.5.1) was used for all statistical tests. A chi-square test or Fisher's exact test was used for the categorical variable. A student's t-test, Mann-Whitney U-test or Kruskal-Wallis H-test was used for the continuous variable. The “mRMRe” package was used to perform the mRMR model, and the “glmnet” package was used to conduct the LASSO model. The “pROC” package was used to plot the ROC curves and the “rms” package was used to build nomogram and perform calibration curves. The ROC curve analysis was performed using the “ROC.TEST” packages. Meanwhiles, the “generalhoslem” package and the “dca.R.” package were used to conduct the H-L test and DCA, respectively. A two-sided P < 0.05 was considered statistically significant.

Results

Clinical Characteristics and Development of the Clinical Model

The baseline of GGNs is shown in Tables 1, 2. There were no significant differences between the training and validation cohorts (Table 1). The difference in clinical characteristics between non-growth group and growth group was shown in Table 2. The univariate logistic regression analysis demonstrated the type, diameter, solid-diameter, and volume as risk signatures for the growth of GGNs. However, only the diameter was still considered as a viable predictive indicator in the clinical model after using the stepwise multivariate logistic regression analysis as shown in Table 3.

TABLE 1
www.frontiersin.org

Table 1. Characteristics of the GGNs in the training and validation group.

TABLE 2
www.frontiersin.org

Table 2. Characteristics of the non-growth and growth cohorts.

TABLE 3
www.frontiersin.org

Table 3. Stepwise multivariate logistic regression analysis.

Features Selection and Construction of Model

By using the mRMR method, 20 features were retained. Then, five features after regression by the LASSO model. These features were presented in the rad-score and were calculated by using the following formula: rad-score = 0.088*RunLengthNonuniformity_AllDirection_offset7_SD-0.367*SurfaceVolumeRatio-0.214*LongRunLowGreyLevelEmphasis_angle0_offset1+0.03*ShortRunEmphasis_AllDirection_offset1_SD+0.227*VolumeCC-0.753. As seen in Figure 2, rad-score was significantly different between the growth and non-growth groups in both the primary and test cohort (P < 0.01) when using the Wilcoxon rank-sum test. The rad-score [odds ratio (OR) = 5.130; 95%CI: 0.948–37.835; P < 0.01] and diameter (OR = 1.087; 95%CI: 0.785–1.564; P < 0.05) were both considered as predictive indicators for the growth of GGNs by using the multivariate logistic regression analysis as seen in Table 4. In the combined model, the rad-score was the key predictive indicator of the growth of GGNs.

FIGURE 2
www.frontiersin.org

Figure 2. Feature selection and the performance of rad-score and three models. (A) The least absolute shrinkage and selection operator (LASSO) regression was used to choose the features to construct the final model. (B,C) The rad-score from non-growth (class 0) and growth (class 1) on the training group (B) and test group (C) were compared, respectively. (D,E) Receiver operating characteristic (ROC) curve for predicting the growth of GGNs in the training group (D) and test group (E).

TABLE 4
www.frontiersin.org

Table 4. Risk factors for the growth of GGNs.

Model Comparison and Construction of the Nomogram

For the training group, the AUC of the combined model, the radiomic model, and the clinical model were 0.801, 0.803, and 0.741, respectively. For the validation cohort, the AUC of the combined model, the radiomic model, and the clinical model were 0.782, 0.791, and 0.686, respectively. There were significant differences between the ROC of the combined model and clinical model (Z = 1.987, P = 0.047). No significant differences were found between the ROC of the combined model and radiomics model (Z = −0.490, P = 0.624). Meanwhile, the combined model also showed the greatest accuracy (accuracy: 80.8%; sensitivity: 86.7%; specificity: 79.4%; PPV: 50.0%; NPV: 96.2%) in the prediction of the growth of GGNs (Table 5). Therefore, the nomogram was generated based on the combined model (Figure 3). Compared to the diameter, the rad-score made up a high proportion of the nomogram.

TABLE 5
www.frontiersin.org

Table 5. Accuracy and predictive value between three models.

FIGURE 3
www.frontiersin.org

Figure 3. The evaluation of the degree of fitting for the combined model and comparison of clinical utility of three models. (A) The combined nomogram based on clinical factors and rad-score for predicting the growth of GGNs. (B,C) Calibration curves for prediction of the growth of GGNs based on the combined model in two cohorts. The X-axis represents the predicted probability of GGN growth based on the combined model and the Y-axis is the actual probability for the growth of GGNs. (D) The X-axis represents high-risk threshold and the Y-axis represents net benefit. The green line represents the clinical model. The blue and red line, respectively represent the radiomic model and the combined model. The black line represents a hypothetical GGN growth. The yellow line a non-growing GGN.

As the calibration plots illustrates, there is a high degree of consistency between actual observation and the combined model in both the training and the validation cohort (Figure 3). Meanwhile, the results of the H-L test were non-significant statistics in the training group (P = 0.6305) and test group (P = 0.6698), which represented a good fitting model.

Clinical Use of DCA

The DCA based on three models was shown in Figure 3. However, the DCA based on the radiomic model showed a greater benefit in the prediction of GGN growth in the 10–60% threshold probabilities as opposed to the clinical model. In essence, the diagnostic utility of the rad-score surpassed that of the clinical model within this threshold range. Combining the clinical features and rad-score allowed for a DCA similar to that of the radiomic model.

Discussion

In this study, the predictive value of radiomic features was analyzed to determine their utility in predicting the growth of GGNs. For the combined model, an improvement in diagnostic utility was also evident through the rad-score for predicting the growth of GGNs. The combined model showed the performance with an AUC of 0.801. Then, the results of DCA showed that the combined model and the radiomic model demonstrated greater performance over the clinical model in predicting the growth of GGNs (Figure 3).

Before this study, the majority of investigations used conventional CT features to predict the growth of GGNs. Parameters, such as the diameter, CT attenuation, volume, and shape were commonly used (21, 22). For example, Masaya Tamura et al. retrospectively analyzed the potential value of conventional CT features including diameter, and mean CT attenuation in the prediction of the growth of GGNs (22). Their result showed the mean CT attenuation (OR = 7.572; P = 0.0023) was the best predictor of the growth of GGNs by using multivariate analysis. However, in our study, the diameter of GGNs was found to be significant indicator of GGN growth in the clinical model after using a stepwise multivariate logistic regression analysis. A potential explanation is that conventional CT features require visual inspection and measurement on a macro scale. In essence, features less than the resolution of the equipment are missed. CT radiomic features offer the advantage of extracting high-throughput data from CT images (23). Therefore, many studies using radiomic methods to analyze GGNs have been reported (2426). A retrospective study by Li Fan et al. analyzed the value of radiomic features as a potential biomarker for the distinction between pre-invasive lesions and invasive lung adenocarcinoma appearing as GGNs (25). A separate study by Q. Sun et al. sought to investigate the relationship between CT texture features and the growth trends of GGNs (27). Features quantitatively extracted from 89 GGNs including the mean value, uniformity, entropy, and energy, were used in their analysis. Their results showed that there was a significant relationship between uniformity and volume doubling time for pure GGNs (P = 0.022). In our study, the radiomic model were compared alongside the clinical and combined model for the evaluation of GGNs. The SurfaceVolumeRatio and VolumeCC were considered key parameters for the accuracy of the proposed model according to its corresponding coefficients. Specifically, the SurfaceVolumeRatio was defined as the surface area of lesions in square millimeters divided by the volume of lesions in cubic centimeters using AK. The SurfaceVolumeRatio was inversely related to the rad-score using corresponding coefficients, which suggested that the value of SurfaceVolumeRatio may be large in non-growth GGNs. According to the general pathophysiology, the SurfaceVolumeRatio is relatively small in irregular GGN which may indicate heterogeneity and malignancy (28). Therefore, the SurfaceVolumeRatio is inversely related to the degree of malignancy, further confirming the findings made herein. According to definitions of these radiomic features (Supplementary Descriptions), the RunLengthNonuniformity, LongRunLowGreyLevelEmphasis, and ShortRunEmphasis also reflect the heterogeneity of lesions. The definition of VolumeCC was set as the volume of GGNs in cubic centimeters using AK. Based on the corresponding coefficients, the value of VolumeCC was positively correlated with the rad-score in our study, suggesting that the value of VolumeCC may be larger in the growth group. Results by Jacob Scharcanski et al. indicated that the growth pattern of pulmonary nodules is exponential (29). Therefore, GGNs with large volumes will have a faster growth rate over smaller nodules. These findings are in line with those obtained herein. Additionally, the mRMR method was used to select radiomic features, which can then evaluate the correlation between the radiomic features and results and the correlation between different radiomic features (30). The aim was to select the features most relevant to the results and remove redundant features. In the current study, the top 20 features selected from 396 total features were used in the mRMR model. By doing so, the accuracy of feature selection was greatly increased. Once the nomogram was developed, a calibration curve and H-L test were performed to evaluate the predictive model (31). The results of both studies suggested that the combined model showed a good correlation with the actual data.

Despite the findings presented herein, several limitations are of note. Firstly, the design of the retrospective study meant that sample sizes were small as a result of the strict inclusion criteria. There was also a lack of external validation as data was gathered from a single institution. Secondly, the growth of GGNs was only measured in two dimensions. Yet, a volume change may better reflect the growth of GGNs as their growth can often be in asymmetric axis (32). Several studies have also investigated the use of volume double time and mass double time to reflect the growth rate of GGNs (33, 34). Yet, current controversy exists as to the best method to measure the natural history of GGNs. The use of manual segmentation in the current study also predisposes results to observer bias. Furthermore, clinical confounding variables, such as patients' social and family history were not included in this analysis, as information regarding such variables is insufficient in the current long-term study. Specifically, the smoking history of patients was not elicited, despite smoking history is a potential confounding factor in the growth of GGNs (3537). Finally, pathological data from GGNs were not obtained, meaning that one cannot infer as to the relationship between the results presented and the types of GGNs with their respective growth trends.

In conclusion, this study has developed and compared three models for predicting the growth of GGNs; the current clinical model, the radiomic model and a combination of the two. The results suggested that the radiomic model and the combined model showed increased utility in predicting the growth of GGNs as opposed to the clinical model. The nomogram studies conducted herein suggested the combined model as offering the greatest diagnostic value. Therefore, this study indicates the utility and versatility of the combined model in guiding the management of GGNs.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics Statement

The studies involving human participants were reviewed and approved by the Ethics Committee of the First Affiliated Hospital of Zhejiang Chinese Medical University. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.

Author Contributions

CG: investigation, formal analysis, data curation, and writing–original draft. JY, YL, and LW: investigation and visualization. PP: software, validation, formal analysis, and visualization. PX: investigation. MX: conceptualization, methodology, resources, writing–reviewing and editing, supervision, and project administration. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by the Zhejiang Provincial Natural Science Foundation of China (LSY19H180003) and the Medical Health Science and Technology Project of Zhejiang Provincial Health Commission (2019KY117, WKJ-ZJ-2039).

Conflict of Interest

PP was employed by GE Healthcare Life Sciences.

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.

Acknowledgments

We thank Michael C. Lin for language editing. In addition, CG especially wants to thank Starbucks for their Latte and Cold Brew that have given him strong support.

Supplementary Material

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

Abbreviations

AIC, akaike information criterion; AK, Analysis Kit; AUC, area under the curve; CI, confidence interval; CT, computed tomography; DCA, decision curve analysis; GGNs, ground glass nodules; GLCM, gray-level co-occurrence matrix; H-L, Hosmer-Lemeshow; LASSO, least absolute shrinkage and selection operator; LLL, left lower lobe; LUL, left upper lobe; MRMR, minimum redundancy maximum relevance; NA, not available; NPV, negative-predictive value; OR, odds ratio; PPV, positive-predictive value; RLL, right lower lobe; RLM, gray level run-length matrix; RML, right middle lobe; ROC, receiver operating characteristic; ROI, regions of interest; RUL, right upper lobe.

References

1. Kakinuma R, Ohmatsu H, Kaneko M, Kusumoto M, Yoshida J, Nagai K, et al. Progression of focal pure ground-glass opacity detected by low-dose helical computed tomography screening for lung cancer. J Comput Assist Tomogr. (2004) 28:17–23. doi: 10.1097/00004728-200401000-00003

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Kovalchik SA, Tammemagi M, Berg CD, Caporaso NE, Riley TL, Korch M, et al. Targeting of low-dose CT screening according to the risk of lung-cancer death. N Engl J Med. (2013) 369:245–54. doi: 10.1056/NEJMoa1301851

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Collins J, Stern EJ. Ground-glass opacity at CT: the ABCs. Am J Roentgenol. (1997) 169:355–67. doi: 10.2214/ajr.169.2.9242736

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Engeler CE, Tashjian JH, Trenkner SW, Walsh JW. Ground-glass opacity of the lung parenchyma: a guide to analysis with high-resolution CT. Am J Roentgenol. (1993) 160:249–51. doi: 10.2214/ajr.160.2.8424326

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Park CM, Goo JM, Lee HJ, Lee CH, Chun EJ, Im J-G. Nodular ground-glass opacity at thin-section CT: histologic correlation and evaluation of change at follow-up. Radiographics. (2007) 27:391–408. doi: 10.1148/rg.272065061

PubMed Abstract | CrossRef Full Text | Google Scholar

6. 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

7. Diederich S. Pulmonary nodules: do we need a separate algorithm for non-solid lesions? Cancer Imaging. (2009) 9:S126–8. doi: 10.1102/1470-7330.2009.9046

PubMed Abstract | CrossRef Full Text | Google Scholar

8. MacMahon H, Naidich DP, Goo JM, Lee KS, Leung ANC, Mayo JR, et al. Guidelines for management of incidental pulmonary nodules detected on CT images: from the Fleischner Society 2017. Radiology. (2017) 284:228–43. doi: 10.1148/radiol.2017161659

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Lee HJ, Goo JM, Lee CH, Yoo C-G, Kim YT, Im J-G. Nodular ground-glass opacities on thin-section CT: size change during follow-up and pathological results. Korean J Radiol. (2007) 8:22–31. doi: 10.3348/kjr.2007.8.1.22

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Eguchi T, Kondo R, Kawakami S, Matsushita M, Yoshizawa A, Hara D, et al. Computed tomography attenuation predicts the growth of pure ground-glass nodules. Lung Cancer. (2014) 84:242–7. doi: 10.1016/j.lungcan.2014.03.009

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Lee SW, Leem CS, Kim TJ, Lee KW, Chung JH, Jheon S, et al. The long-term course of ground glass opacities detected on thin-section computed tomography. Respir Med. (2013) 107:904–10. doi: 10.1016/j.rmed.2013.02.014

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Matsuguma H, Mori K, Nakahara R, Suzuki H, Kasai T, Kamiyama Y, et al. Characteristics of subsolid pulmonary nodules showing growth during follow-up with CT scanning. Chest. (2013) 143:436–43. doi: 10.1378/chest.11-3306

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Gao C, Li J, Wu L, Kong D, Xu M, Zhou C. The natural growth of subsolid nodules predicted by quantitative initial CT features: a systematic review. Front Oncol. (2020) 10:318. doi: 10.3389/fonc.2020.00318

PubMed Abstract | CrossRef Full Text | Google Scholar

14. 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

15. Tunali I, Gray JE, Qi J, Abdalah M, Jeong DK, Guvenis A, et al. Novel clinical and radiomic predictors of rapid disease progression phenotypes among lung cancer patients treated with immunotherapy: an early report. Lung Cancer. (2019) 129:75–9. doi: 10.1016/j.lungcan.2019.01.010

PubMed Abstract | CrossRef Full Text | Google Scholar

16. 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:783–92. doi: 10.1148/radiol.2018180910

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Hiramatsu M, Inagaki T, Inagaki T, Matsui Y, Satoh Y, Okumura S, et al. Pulmonary ground glass opacity (GGO) lesions-large size and a history of lung cancer are risk factors for growth. J Thorac Oncol. (2008) 3:1245–50. doi: 10.1097/JTO.0b013e318189f526

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Unler A, Murat A, Chinnam RB. mr2PSO: a maximum relevance minimum redundancy feature selection method based on swarm intelligence for support vector machine classification. Inf Sci. (2011) 181:4625–41. doi: 10.1016/j.ins.2010.05.037

CrossRef Full Text | Google Scholar

19. Iasonos A, Schrag D, Raj GV, Panageas KS. How to build and interpret a nomogram for cancer prognosis. J Clin Oncol. (2008) 26:1364–70. doi: 10.1200/JCO.2007.12.9791

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Stephenson AJ, Scardino PT, Eastham JA, Bianco FJ Jr, Dotan ZA, DiBlasio CJ, et al. Postoperative nomogram predicting the 10-year probability of prostate cancer recurrence after radical prostatectomy. J Clin Oncol. (2005) 23:7005–12. doi: 10.1200/JCO.2005.01.867

CrossRef Full Text | Google Scholar

21. Bak SH, Lee HY, Kim JH, Um SW, Kwon OJ, Han J, et al. Quantitative CT scanning analysis of pure ground-glass opacity nodules predicts further CT scanning change. Chest. (2016) 149:180–91. doi: 10.1378/chest.15-0034

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Tamura M, Shimizu Y, Yamamoto T, Yoshikawa J, Hashizume Y. Predictive value of one-dimensional mean computed tomography value of ground-glass opacity on high-resolution images for the possibility of future change. J Thorac Oncol. (2014) 9:469–72. doi: 10.1097/JTO.0000000000000117

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Thawani R, Mclane M, Beig N, Ghose S, Prasanna P, Velcheti V, et al. Radiomics and radiogenomics in lung cancer : a review for the clinician. Lung Cancer. (2018) 115:34–41. doi: 10.1016/j.lungcan.2017.10.015

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Chae H-D, Park CM, Park SJ, Lee SM, Kim KG, Goo JM. Computerized texture analysis of persistent part-solid ground-glass nodules: differentiation of preinvasive lesions from invasive pulmonary adenocarcinomas. Radiology. (2014) 273:285–93. doi: 10.1148/radiol.14132187

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Fan L, Fang M, Li Z, Tu W, Wang S, Chen W, et al. Radiomics signature: a biomarker for the preoperative discrimination of lung invasive adenocarcinoma manifesting as a ground-glass nodule. Eur Radiol. (2019) 29:889–97. doi: 10.1007/s00330-018-5530-z

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Gao C, Xiang P, Ye J, Pang P, Wang S, Xu M. Can texture features improve the differentiation of infiltrative lung adenocarcinoma appearing as ground glass nodules in contrast-enhanced CT? Eur J Radiol. (2019) 117:126–31. doi: 10.1016/j.ejrad.2019.06.010

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Sun Q, Huang Y, Wang J, Zhao S, Zhang L, Tang W, et al. Applying CT texture analysis to determine the prognostic value of subsolid nodules detected during low-dose CT screening. Clin Radiol. (2019) 74:59–66. doi: 10.1016/j.crad.2018.07.103

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Niu R, Shao X, Shao X, Wang J, Jiang Z, Wang Y. Lung adenocarcinoma manifesting as ground-glass opacity nodules 3 cm or smaller: evaluation with combined high-resolution CT and PET/CT modality. Am J Roentgenol. (2019) 213:W236–45. doi: 10.2214/AJR.19.21382

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Scharcanski J, da Silva LS, Koff D, Wong A. Interactive modeling and evaluation of tumor growth. J Digit Imaging. (2010) 23:755–68. doi: 10.1007/s10278-009-9234-4

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Parmar C, Grossmann P, Rietveld D, Rietbergen MM, Lambin P, Aerts HJWL. Radiomic machine-learning classifiers for prognostic biomarkers of head and neck cancer. Front Oncol. (2015) 5:272. doi: 10.3389/fonc.2015.00272

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Xu L, Yang P, Liang W, Liu W, Wang W, Luo C, et al. A radiomics approach based on support vector machine using MR images for preoperative lymph node status evaluation in intrahepatic cholangiocarcinoma. Theranostics. (2019) 9:5374–85. doi: 10.7150/thno.34149

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Yanagawa M, Tanaka Y, Kusumoto M, Watanabe S, Tsuchiya R, Honda O, et al. Automated assessment of malignant degree of small peripheral adenocarcinomas using volumetric CT data: correlation with pathologic prognostic factors. Lung Cancer. (2010) 70:286–94. doi: 10.1016/j.lungcan.2010.03.009

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Borghesi A, Farina D, Michelini S, Ferrari M, Benetti D, Fisogni S, et al. Pulmonary adenocarcinomas presenting as ground-glass opacities on multidetector CT: three-dimensional computer-assisted analysis of growth pattern and doubling time. Diagn Interv Radiol. (2016) 22:525–33. doi: 10.5152/dir.2016.16110

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Oda S, Awai K, Murao K, Ozawa A, Utsunomiya D, Yanaga Y, et al. Volume-doubling time of pulmonary nodules with ground glass opacity at multidetector CT: Assessment with computer-aided three-dimensional volumetry. Acad Radiol. (2011) 18:63–9. doi: 10.1016/j.acra.2010.08.022

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Kobayashi Y, Sakao Y, Deshpande GA, Fukui T, Mizuno T, Kuroda H, et al. The association between baseline clinical-radiological characteristics and growth of pulmonary nodules with ground-glass opacity. Lung Cancer. (2014) 83:61–6. doi: 10.1016/j.lungcan.2013.10.017

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Chang B, Hwang JH, Choi YH, Chung MP, Kim H, Kwon OJ, et al. Natural history of pure ground-glass opacity lung nodules detected by low-dose CT scan. Chest. (2013) 143:172–8. doi: 10.1378/chest.11-2501

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Seki N, Sawada S, Nakata M, Inoue T, Nishimura R, Segawa Y, et al. Lung cancer with localized ground-glass attenuation represents early-stage adenocarcinoma in nonsmokers. J Thorac Oncol. (2008) 3:483–90. doi: 10.1097/JTO.0b013e31816a4994

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: machine learning, growth, solitary pulmonary nodule, tomography, X-ray computed, nomograms

Citation: Gao C, Yan J, Luo Y, Wu L, Pang P, Xiang P and Xu M (2020) The Growth Trend Predictions in Pulmonary Ground Glass Nodules Based on Radiomic CT Features. Front. Oncol. 10:580809. doi: 10.3389/fonc.2020.580809

Received: 07 July 2020; Accepted: 10 September 2020;
Published: 20 October 2020.

Edited by:

Jiuquan Zhang, Chongqing University, China

Reviewed by:

Chen Chen, Central South University, China
Zhongxiang Ding, Zhejiang University, China

Copyright © 2020 Gao, Yan, Luo, Wu, Pang, Xiang and Xu. 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: Maosheng Xu, xums166@zcmu.edu.cn