Skip to main content

ORIGINAL RESEARCH article

Front. Oncol., 19 November 2019
Sec. Cancer Imaging and Image-directed Interventions
This article is part of the Research Topic Novel Methods for Oncologic Imaging Analysis: Radiomics, Machine Learning, and Artificial Intelligence View all 38 articles

CT Morphological Features Integrated With Whole-Lesion Histogram Parameters to Predict Lung Metastasis for Colorectal Cancer Patients With Pulmonary Nodules

\nTingDan Hu,&#x;TingDan Hu1,2ShengPing Wang,&#x;ShengPing Wang1,2Xiangyu EXiangyu E3Ye YuanYe Yuan3Lv HuangLv Huang3JiaZhou WangJiaZhou Wang3DeBing ShiDeBing Shi4Yuan LiYuan Li5WeiJun Peng,WeiJun Peng1,2Tong Tong,
Tong Tong1,2*
  • 1Department of Radiology, Fudan University Shanghai Cancer Center, Shanghai, China
  • 2Department of Oncology, Shanghai Medical College, Fudan University, Shanghai, China
  • 3Department of Radiotherapy, Fudan University Shanghai Cancer Center, Fudan University, Shanghai, China
  • 4Department of Colorectal Surgery, Fudan University Shanghai Cancer Center, Fudan University, Shanghai, China
  • 5Department of Pathology, Fudan University Shanghai Cancer Center, Fudan University, Shanghai, China

Purpose: To retrospectively identify the relationships between both CT morphological features and histogram parameters with pulmonary metastasis in patients with colorectal cancer (CRC) and compare the efficacy of single-slice and whole-lesion histogram analysis.

Methods: Our study enrolled 196 CRC patients with pulmonary nodules (136 in the training dataset and 60 in the validation dataset). Twenty morphological features of contrast-enhanced chest CT were evaluated. The regions of interests were delineated in single-slice and whole-tumor lesions, and 22 histogram parameters were extracted. Stepwise logistic regression analyses were applied to choose the independent factors of lung metastasis in the morphological features model, the single-slice histogram model and whole-lesion histogram model. The areas under the curve (AUC) was applied to quantify the predictive accuracy of each model. Finally, we built a morphological-histogram nomogram for pulmonary metastasis prediction.

Results: The whole-lesion histogram analysis (AUC of 0.888 and 0.865 in the training and validation datasets, respectively) outperformed the single-slice histogram analysis (AUC of 0.872 and 0.819 in the training and validation datasets, respectively) and the CT morphological features model (AUC of 0.869 and 0.845 in the training and validation datasets, respectively). The morphological-histogram model, developed with significant morphological features and whole-lesion histogram parameters, achieved favorable discrimination in both the training dataset (AUC = 0.919) and validation dataset (AUC = 0.895), and good calibration.

Conclusions: CT morphological features in combination with whole-lesion histogram parameters can be used to prognosticate pulmonary metastasis for patients with colorectal cancer.

Introduction

Colorectal cancer (CRC) is the third common cause of morbidity and mortality worldwide (1, 2). Pulmonary is the most common extra-abdominal site of metastasis for those with CRC, with 5–10% of CRC patients developing pulmonary metastasis (PM) (3, 4). The 5-year survival rates after initial colorectal surgery in patients with and without resection for pulmonary metastasis are 68 and 13%, respectively (3). The strong survival benefits of pulmonary metastasectomy make this treatment the generally accepted treatment for patients to achieve long-term survival when there is a definite and clear diagnosis (5, 6). Furthermore, if pulmonary metastasis is diagnosed early and resected aggressively, the survival rate is further improved (7).

However, with chest CT applied as part of preoperative routine examination, an increasing number of CRC patients are being diagnosed with indeterminate pulmonary nodules (IPNs) of unknown nature (8). The reported incidence of IPNs in CRC patients is 25–45.5% (810). Further diagnostic tests can also be problematic as nodules <10 mm in diameter may fall below the threshold of detection for positron emission tomography (PET) (11), and fine-needle aspiration cytology may not be feasible for thoracoscopic localization (12). Therefore, in CRC patients with IPNs, the accurate diagnosis of metastatic disease at an early and surgically treatable stage remains a challenge.

Though early-stage metastatic nodules and benign lesions have similar appearance in images, the importance of morphology should not be underestimated (13). CT imaging allows detailed observation of the morphological features of nodules and lesions, such as their internal density, shape, margin, and other typical characteristics. In recent years, texture analysis has emerged as a valuable methodology for facilitating diagnosis through the deep mining of information from medical images (14, 15). It has achieved great utility in evaluating many kinds of pulmonary diseases, including pulmonary embolisms (16), interstitial lung disease (17), and pulmonary nodules (18, 19). By extracting features of subtle pixel distributions and spatial variations of the gray levels of lesions that are imperceptible to the naked eye, texture analysis provides a complementary method for evaluating subjective and megascopic morphological features.

To date, studies concentrating on the morphological and textural features of IPNs 5–20 mm in diameter on contrast-enhanced CT in CRC patients remain limited. This study sought to determine the morphological characteristics and histogram parameters derived from texture analysis for CRC patients with IPNs and to construct a risk model with a combination of independent predictors to facilitate the accurate diagnosis of pulmonary metastasis.

Materials and Methods

Patients

This retrospective analysis had obtained the ethical approval, and the informed consent requirement was waived. Our study enrolled 196 consecutive colorectal cancer patients (88F/108M; age range, 32–80 years; mean age, 58.49 ± 10.80 years) with lung nodules admitted in our institution between January 2010 and December 2017. The inclusion criteria were as follows: (i) colorectal cancer was histopathologically confirmed; (ii) at least one lung nodule measuring 5–20 mm detected by contrast-enhanced chest CT examination; (iii) available pathology reports with diagnosis of pulmonary metastasis or primary lung cancer for the malignant nodules and at least 2 years follow-up for the benign nodules; and (iv) complete medical history. The exclusion protocol were as follows: (i) with pretreatment 6 months before initial CT examination (including chemotherapy or pneumonectomy); (ii) obsolete nodules detected 6 months before colorectal cancer was detected; (iii) obvious benign nodules with typical imaging characteristics (such as cysts, tuberculosis, or inflammatory nodules); and (iv) adjuvant therapy (including radiation therapy or chemotherapy) applied for no-progress lesions in the process of follow-up. When there are multiple nodules, we choose the largest nodule for morphological and radiomics analysis. Of the 196 people included in the study, 194 of them have been published in our previous research (20).

Nodules were divided into two groups: (i) a pathologically confirmed lung metastasis group (95 PMs; 42F/53M; mean age, 57.46 ± 10.58 years), and (ii) a non-metastasis (NM) group (101 NMs; 46F/55M; mean age, 59.47 ± 10.91 years), including benign nodules (90 cases) with at least 2 years follow-up (88 cases) and pathology confirmation (2 cases) or primary lung cancer confirmed by pathology (11 cases). We used a computer algorithm to randomly divide the patients into a training dataset and a validation dataset at the ratio of 7:3. Figure 1 shows the process of patients' recruitment.

FIGURE 1
www.frontiersin.org

Figure 1. Flow chat of patients' recruitment pathway.

CT Scanning Protocol

Chest CT examinations were performed at our institution with the Sensation 64 scanner (Siemens Healthcare) or the Somatom Definition AS scanner (Siemens Healthcare). The Contrast-enhanced CT scan parameters were as follows: contrast medium, inhexol; tube voltage, 120 kVp; tube current, 250–350 mA; slice thickness, 1.5 mm; slice interval, 1.5 mm; matrix, 512 × 512; field of view (FOV, 35–50 cm; pitch, 1.078; reconstruction algorithm, standard. The arterial phase of the target nodule which was pathologically confirmed or under follow-up was selected for reconstruction.

CT Image Interpretation

The interpretations of CT features are listed in Supplementary Table 1. The CT morphological features were independently evaluated by two operators (SW and TH, with 20 and 3 years of experience in chest CT, respectively). In cases of disagreement, a third radiologist (TT, with 20 years of experience in CT imaging) was consulted, and the majority value was used. Mean values were calculated for continuous variables. The CT images were read with both mediastinal and lung window settings. All of the operators were blinded to the clinical and histologic findings.

Histogram Analysis

Reconstructed images were transferred to the MIM software (v6.6.3; MIM Software Inc.) for histogram analysis. For each patient, regions of interest (ROIs) were first semi-automatically contoured in the largest-cross sectional area of the tumor outline and then manually delineated by an operator and verified by an expert radiologist. Each ROI was propagated to include the entire tumor volume in each consecutive slice using the same contouring method. In the process of delineation, we excluded the border of the lesion and any other irrelevant tissues or regions, such as pleura, normal tissue, air, peripheral vessels, and surrounding organs. Supplementary Figure 1 shows an example of ROI delineation.

The histogram parameters were automatically measured by the software using a volumetric approach on the ROI of the nodule. Single-slice and whole-lesion histogram parameters were extracted and analyzed. From each segmented tumor, we extracted 11 single-slice histogram parameters and 11 whole-lesion histogram features. More information about the methodology used to extract histogram features can be found in Supplementary Material.

Statistical Analysis

R software (version 3.3) was applied for statistical analysis. To measure the agreement of CT morphological features between two readers, intraclass correlation coefficients (ICCs) were calculated (poor: 0.00–0.20; fair: 0.21–0.40; moderate: 0.41–0.60; good: 0.61–0.80; excellent: 0.81–1.00). To compare the proportional differences between the training dataset and the validation dataset, chi-square tests were applied for the categorical variables, and two-sample t-tests were used for the continuous variables. To compare the differences between the PM and NM group, chi-square and two-sample t-tests were applied as appropriate for both the training and validation datasets. Two-sided p < 0.05 was considered significant.

Model Selection

The significant factors were introduced into the stepwise logistic regression to select the independent features for the CT morphological model, the single-slice histogram model and the whole-lesion histogram model. The Akaike information criterion (AIC) was employed as the stopping rule. The validation dataset was used to test the diagnostic performance of the models by applying the multivariable regression formula derived from the training dataset to the patients of the validation dataset, and the probability of metastasis was calculated for each. The area under the receiver characteristic curve (AUC) was calculated to quantify the predictive accuracy of the three models in the training and validation datasets. We also calculated the accuracy, sensitivity, specificity, positive predictive value, and negative predictive value for each model.

We compared the relative strengths of the single-slice and whole-lesion histogram models and then used the more efficient model in combination with the morphological features to construct the morphological-histogram model. A morphological-histogram nomogram was then constructed for clinical application. A receiver operating characteristic (ROC) curve was used to describe the discrimination abilities of the nomogram. An AUC above 0.75 is considered as good (21). Nomogram performance was graphically demonstrated by calibration plots in both the training and validation datasets. Finally, decision curve analysis (DCA) was applied to assess the clinical usefulness of the nomogram.

Results

Patient Characteristics

The patients characteristics and statistically significant CT morphological features are shown in Table 1 (Supplementary Table 3 contains complete morphological features comparison), and the histogram parameters are presented in Table 2. There were no significant differences between the training and validation datasets except in pleural attachment (Supplementary Table 2). The agreement between the two operators was excellent for most characteristics and good for several features (Supplementary Table 4).

TABLE 1
www.frontiersin.org

Table 1. Comparison of morphological features of lung metastasis (PM) and non- metastasis (NM) in the training and validation datasets.

TABLE 2
www.frontiersin.org

Table 2. Comparison of single-slice and whole-lesion histogram parameters of PM and NM in the training and validation datasets.

Significant Morphological Features and Histogram Parameters

Regarding the CT morphological features, the chi-square tests and t-tests revealed that nine CT features were associated with lung metastasis, including lobe location (p = 0.025), size category (p < 0.001), long-axis diameter (p < 0.001), short-axis diameter (p = 0.001), density (p < 0.001), contour (p = 0.043), border (p < 0.001), air bronchogram (p = 0.032), and lymphadenopathy (p = 0.032). After stepwise logistic analysis, long-axis diameter (OR = 1.360, 95%CI: 1.198–1.544, P < 0.001), density (OR = 11.166, 95%CI: 2.721–45.815, P < 0.001) and contour (OR = 0.317, 95%CI: 0.177–0.569, P = 0.001) remained independent predictors in the CT morphological model, as shown in Table 3.

TABLE 3
www.frontiersin.org

Table 3. Comparison of the models by multivariate logistic regression analysis.

Regarding the single-slice histogram parameters (S- means the parameters from the single-slice histogram analysis and W- from the whole-slice histogram), t-tests revealed that the S-average ratio (p < 0.001), S-skewness (p < 0.001), S-area (p < 0.001), S-volume (p < 0.001), S-median (p < 0.001), and S-minimum (p = 0.035) were significant variables related to PM. After stepwise logistic analysis, the S-average ratio (OR = 0.268, 95%CI: 0.111–0.642, P = 0.003), S-area (OR = 559.372, 95%CI: 42.344–7389.333, P < 0.001), and S-median (OR = 1.004, 95%CI: 1.002–1.005, P < 0.001) were selected as independent predictors for the single-slice histogram model.

Regarding the whole-lesion histogram parameters, t-tests revealed that the W-average ratio (p < 0.001), W-mean (p = 0.002), W-skewness (p < 0.001), W-area (p < 0.001), S-volume (p = 0.001), W-median (p < 0.001), and S-minimum (p = 0.035) were significant parameters associated with PM. After stepwise logistic analysis, the W-average ratio (OR = 12.764, 95%CI: 4.653–35.018, P = 0.003), W-mean (OR = 0.977, 95%CI: 0.961–0.994, P = 0.004), and S-median (OR = 1.024, 95%CI: 1.008–1.041, P = 0.009) were selected as independent predictors for the whole-lesion histogram model. Figure 2 shows the distributions of the significant histogram parameters in the training and validation datasets.

FIGURE 2
www.frontiersin.org

Figure 2. Boxplot of the selected histogram parameters in the LM and NM group. (A-C) Boxplot of the S-Average ratio, S-Median, S-Area from the single-slice histogram in the LM and NM group, respectively. (D-F) Boxplot of the W-Average ratio, W-Median, W-Mean from the whole-lesion histogram in the LM and NM group, respectively.

Comparison of Single-Slice and Whole-Lesion Histogram Analyses

The whole-lesion histogram model (AIC = 130.25) had lower AIC value than the single-slice model (AIC = 130.9) and achieved better discrimination. It yielded an AUC of 0.888 for the training dataset and of 0.865 for the validation dataset, exceeding the AUC values of the single-slice model (AUC = 0.872 for the training dataset and AUC = 0.819 for the validation dataset). The ROC curves of the two models are presented in Figure 3.

FIGURE 3
www.frontiersin.org

Figure 3. (A,B) The ROC curves of the single-slice histogram model and whole-lesion histogram model in the training (A) and validation sets (B), respectively. (C,D) The ROC curves of the CT morphological model, the whole-lesion histogram model and the integrated morphological-histogram model in the training (C) and validation cohorts (D), respectively.

Development and Validation of the Morphological-Histogram Nomogram

We subjected the CT morphological features and the whole-lesion histogram parameters to stepwise logistic regression analysis. Density (OR = 5.434, 95%CI: 1.161–25.440, P = 0.032), contour (OR = 0.495, 95%CI: 0.286–0.858, P = 0.012), the W-average ratio (OR = 9.727, 95%CI: 3.538–26.740, P < 0.001), W-mean (OR = 0.977, 95%CI: 0.959–0.995, P = 0.009), and W-median (OR = 1.023, 95%CI: 1.006–1.042, P = 0.013) were identified as independent risk factors in the model. The integrated model also achieved the best performance among the models, with an AUC of 0.919 (95%CI: 0.871–0.968, accuracy: 88.2%, sensitivity: 84.9%, specificity: 92.1%, PPV: 92.5%, NPV: 84.1%) for the training dataset and of 0.895 (95%CI: 0.813–0.977, accuracy: 81.7%, sensitivity: 78.5%, specificity: 84.4%, PPV: 81.5%, NPV: 81.8%) for the validation dataset (Table 4). The ROC curves of the models for both the training and validation datasets are presented in Figure 3.

TABLE 4
www.frontiersin.org

Table 4. Accuracy and predictive value between those models.

The morphological-histogram nomogram was successfully constructed, with good discrimination, based on the morphological-histogram model (Figure 4A). The calibration plots also presented good accordance between the nomogram prediction and actual outcome for PM and NM in both the training and validation datasets (Figures 4B,C). The decision curve analysis demonstrated that given a threshold probability ranging from 0 to 100%, the morphological-histogram model was superior to the treat-all and treat-none schemes in predicting lung metastasis (Figure 4D).

FIGURE 4
www.frontiersin.org

Figure 4. (A) The developed morphological-histogram nomogram for predicting the probability of pulmonary metastases. By summing the scores of each point and locating on the total score scale, the estimated probability of pulmonary metastases could be determined. (B,C) The Calibration curves for predicting pulmonary metastases in the training and validation cohort. The y axis represents the actual rate of LM. The x axis represents the predicted probability of LM. The ideal line represents a perfect prediction by an ideal model. The apparent line represents the performance of the nomogram model, of which a closer fit to the ideal line represents a better prediction. (D) The decision curves analysis for the morphological-histogram nomogram. The red line represents the net benefit of morphological-histogram model. Across the various threshold probabilities, the morphological-histogram curve showed great net benefit.

Discussion

In the present study, we investigated the imaging characteristics of IPNs 5–20 mm in diameter on initial CT in CRC patients and compared the predictive accuracy of whole-lesion and single-slice histogram parameters. We then constructed a morphological-histogram nomogram using a combination of morphological features and whole-lesion histogram parameters for IPNs. This nomogram may be clinically useful for discriminating CRC patients who might benefit from early and curable metastasectomy for metastatic lesions or an appropriate surveillance program.

CT offers direct visualization of lesions and potentially allows a detailed characterization of the morphologic extent of lesions. The careful evaluation of morphologic features is an essential step in pulmonary nodules assessment (13). Although several studies (2224) have sought to identify significant image features for metastatic nodules, there is no consensus regarding the definition of IPNs, which led to slight differences between our results and previously published ones. In our study, we found that significant morphological features associated with pulmonary metastasis were long-axis diameter, density, and contour.

As reported by many other studies, nodule diameter is a reliable indicator of malignant potential (22, 23, 25). We found that solid nodules are more likely to be metastatic lesions. As more than 95% of nodules that originate from colorectal cancer are adenocarcinomas (4), metastatic lesions tend to appear as solid pulmonary nodules (SPN) in CT scans, whereas benign lesions, such as inflammation lesions, or organizing pneumonia/fibrosis consistently present as patchy consolidations or mixed-density regions surrounded by ground-glass opacity (GGO) owing to inflammatory cell infiltration (26). Primary lung cancer consistently evolves from pre-invasive lesions (AIS/AAH) that manifested as pure GGO (27) at the early stage. A post-hoc analysis (24) found that a solid consistency and increasing size were statistically associated with malignancy.

In addition, our study found that metastatic nodules tended to be round or oval, consistent with previous research (28). We speculate that as metastatic nodules often exhibit a largely uniform growth rate and homogenous invasion in all directions, these features contribute to a round or quasi-circular contour, whereas non-metastatic lesions, including benign lesions and primary lung cancer, have irregular shapes due to uneven growth rates at various sites (26). Thus, short-interval CT follow-up is highly recommended for IPNs larger than 5 mm in diameter with solid components and approximately regular margins detected on preoperative chest CT.

In addition to the identification of morphological features, the use of texture analysis is a strength of our study. Previous studies have demonstrated that texture analysis can not only distinguish malignant nodules from benign ones (18) but also differentiate in situ and minimally invasive lung adenocarcinoma subtypes (19). These studies have shown that texture parameters can reveal the underlying histological changes in tissue below the resolution of the given modality or protocol. In this study, we found that the W-average ratio, W-mean, and W-median, which represent the zone of CT attenuation within the ROI, were substantially higher in the metastasis group than in the non-metastasis group. Thus, short-interval CT follow-up is highly recommended for IPNs larger than 5 mm in diameter with solid components and approximately regular margins detected on preoperative chest CT. This speculation is also in line with another finding of our study that vascular convergence was more common and the enhancement degree was higher in the metastasis group than in the non-metastasis group. However, as texture analysis is a mathematical method, the biological mechanisms underlying the textural features are complex and not completely understood (29). In cases where vascular convergence or the enhancement degree is insufficient to differentiate metastatic lesions, the values from the CT attenuation zone might exhibit local variation and more sensitive preservation of spatial information (30).

Another finding of our study was that the whole-lesion texture analysis outperformed the single-slice analysis in evaluating pulmonary nodules, consistent with a previous study (31). Whole-lesion analysis may provide a more comprehensive understanding of the stereo structure of the whole lesion and thereby reflect the integral heterogeneity better than can single-slice analysis. Despite the time-consuming process of the contouring around the whole lesion, it seems more cost-efficient to use this method as it provides improved prediction relative to single-slice analysis and a more definite diagnosis, allowing timely treatment and maximizing the benefits to the patient.

For clinical use, we constructed a risk stratification nomogram for the clinician to predict the risk of PM for an individual CRC patient. As the early and accurate diagnosis of pulmonary metastasis has been recognized as one of the most important steps in treating potential curable lesions with surgery, we propose that patients with a high risk of PM be considered candidates for thoracotomy for resectable lesions to enhance local control and improve the survival rate. We also hope this model can help low-risk patients avoid aggressive follow-up and reduce the burden of radiation exposure. We believe that the clinical use of the nomogram can contribute to reliable diagnoses and help clinicians optimize therapeutic plans for IPNs at an early stage after detection.

Our study has several limitations. First, as a retrospective study, thin-slice contrast-enhanced CT images from our database were used, which limited the number of cases for analysis. And the inclusion and exclusion criteria also limits the implementation of the study in clinical practice. Second, only histogram parameters were extracted in this study. In our previous research (20), 203 radiomic features, including first- and second-order parameters, attained a prognostic value in the differentiation of pulmonary metastasis with an AUC of 0.888, which is slightly higher than that obtained using the histogram parameters (AUC = 0.887). However, the process of extracting radiomic features through MATLAB is intricate and demanding for radiologists and clinicians, which constrains its clinical utilization. The volume histogram analysis performed here allowed the simple, efficient, and automatic acquisition of a density histogram and achieved an accuracy comparable to that of the radiomics analysis. Thus, volume histogram analysis may be more appropriate for imperative clinical decisions, and radiomics analysis can be used as a supplementary method when needed. Another limitation is that the development and validation were performed in a single institution. External validation and multi-center clinical trials are therefore needed for further generalization.

In conclusion, the results of our study demonstrated that histogram parameters may serve as non-invasive imaging biomarkers for differentiating pulmonary metastasis from non-metastatic lesions. When complemented with morphological features, the morphological-histogram nomogram can greatly benefit the diagnosis of pulmonary metastasis in CRC patients.

Data Availability Statement

The datasets generated for this study are available on request to the corresponding author.

Ethics Statement

The studies involving human participants were reviewed and approved by Medical ethics committee of Fudan University Shanghai Cancer Center. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.

Author Contributions

TT and SW carried out the concepts and design of the study. DS provided the patients information. YL confirmed the pathology results. XE, YY, HL, and JW provided assistance for data acquisition and statistical analysis. WP provided the permission of imaging acquisition. SW and TH carried our literature research and manuscript editing. TH and SW contributed equally to this work. All authors have reviewed the final version of the manuscript and approved it for publication.

Funding

This study has received funding by the National Natural Science Foundation of China (Grant No. 81971687) and the Shanghai Committee of Science and Technology of China (Grant No. 19YF1409900).

Conflict of Interest

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

The reviewer DY declared a shared affiliation, though no other collaboration, with the authors to the handling Editor.

Supplementary Material

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

References

1. Siegel RL, Miller KD, Fedewa SA, Ahnen DJ, Meester RGS, Barzi A, et al. Colorectal cancer statistics, 2017. CA Cancer J Clin. (2017) 67:177–93. doi: 10.3322/caac.21395

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Siegel RL, Miller KD, Jemal A. Cancer Statistics, 2017. CA Cancer J Clin. (2017) 67:7–30. doi: 10.3322/caac.21387

CrossRef Full Text | Google Scholar

3. Kobayashi H, Mochizuki H, Sugihara K, Morita T, Kotake K, Teramoto T, et al. Characteristics of recurrence and surveillance tools after curative resection for colorectal cancer: a multicenter study. Surgery. (2007) 141:67–75. doi: 10.1016/j.surg.2006.07.020

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Watanabe K, Saito N, Sugito M, Ito M, Kobayashi A, Nishizawa Y. Incidence and predictive factors for pulmonary metastases after curative resection of colon cancer. Ann Surg Oncol. (2013) 20:1374–80. doi: 10.1245/s10434-012-2747-y

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Hachimaru A, Maeda R, Suda T, Takagi Y. Repeat pulmonary resection for recurrent lung metastases from colorectal cancer: an analysis of prognostic factors. Interact Cardiovasc Thorac Surg. (2016) 22:826–30. doi: 10.1093/icvts/ivv382

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Cao G, Cheng D, Ye L, Pan Y, Yang F, Lyu S. Surgical resection of pulmonary metastases from colorectal cancer: 11 years of experiences. PLoS ONE. (2017) 12:e0175284. doi: 10.1371/journal.pone.0175284

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Lee JH, Ahn BK, Nam YS, Lee KH. Clinical usefulness of preoperative and postoperative chest computed tomography for colorectal cancer. J Korean Soc Coloproctol. (2010) 26:359–64. doi: 10.3393/jksc.2010.26.5.359

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Nordholm-Carstensen A, Wille-Jorgensen PA, Jorgensen LN, Harling H. Indeterminate pulmonary nodules at colorectal cancer staging: a systematic review of predictive parameters for malignancy. Ann Surg Oncol. (2013) 20:4022–30. doi: 10.1245/s10434-013-3062-y

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Grossmann I, Avenarius JK, Mastboom WJ, Klaase JM. Preoperative staging with chest CT in patients with colorectal carcinoma: not as a routine procedure. Ann Surg Oncol. (2010) 17:2045–50. doi: 10.1245/s10434-010-0962-y

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Maithel SK, Ginsberg MS, D'Amico F, DeMatteo RP, Allen PJ, Fong Y, et al. Natural history of patients with subcentimeter pulmonary nodules undergoing hepatic resection for metastatic colorectal cancer. J Am Coll Surg. (2010) 210:31–8. doi: 10.1016/j.jamcollsurg.2009.09.032

PubMed Abstract | CrossRef Full Text | Google Scholar

11. De Wever W, Meylaerts L, De Ceuninck L, Stroobants S, Verschakelen JA. Additional value of integrated PET-CT in the detection and characterization of lung metastases: correlation with CT alone and PET alone. Eur Radiol. (2007) 17:467–73. doi: 10.1007/s00330-006-0362-7

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Yamauchi Y, Izumi Y, Hashimoto K, Inoue M, Nakatsuka S, Kawamura M, et al. Needle-tract seeding after percutaneous cryoablation for lung metastasis of colorectal cancer. Ann Thorac Surg. (2011) 92:e69–71. doi: 10.1016/j.athoracsur.2011.04.072

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Snoeckx A, Reyntiens P, Desbuquoit D, Spinhoven MJ, Van Schil PE, van Meerbeeck JP, et al. Evaluation of the solitary pulmonary nodule: size matters, but do not ignore the power of morphology. Insights Imaging. (2018) 9:73–86. doi: 10.1007/s13244-017-0581-2

CrossRef Full Text | Google Scholar

14. Lambin P, Rios-Velazquez E, Leijenaar R, Carvalho S, van Stiphout RG, Granton P, et al. Radiomics: extracting more information from medical images using advanced feature analysis. Eur J Cancer. (2012) 48:441–6. doi: 10.1016/j.ejca.2011.11.036

PubMed Abstract | CrossRef Full Text | Google Scholar

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

16. Sumikawa H, Johkoh T, Yamamoto S, Yanagawa M, Inoue A, Honda O, et al. Computed tomography values calculation and volume histogram analysis for various computed tomographic patterns of diffuse lung diseases. J Comput Assist Tomogr. (2009) 33:731–8. doi: 10.1097/RCT.0b013e31818da65c

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Koyama H, Ohno Y, Yamazaki Y, Nogami M, Kusaka A, Murase K, et al. Quantitatively assessed CT imaging measures of pulmonary interstitial pneumonia: effects of reconstruction algorithms on histogram parameters. Eur J Radiol. (2010) 74:142–6. doi: 10.1016/j.ejrad.2009.02.006

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Dennie C, Thornhill R, Sethi-Virmani V, Souza CA, Bayanati H, Gupta A, et al. Role of quantitative computed tomography texture analysis in the differentiation of primary lung cancer and granulomatous nodules. Quant Imaging Med Surg. (2016) 6:6–15. doi: 10.3978/j.issn.2223-4292.2016.02.01

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Cohen JG, Reymond E, Medici M, Lederlin M, Lantuejoul S, Laurent F, et al. CT-texture analysis of subsolid nodules for differentiating invasive from in-situ and minimally invasive lung adenocarcinoma subtypes. Diagn Interv Imaging. (2018) 99:291–9. doi: 10.1016/j.diii.2017.12.013

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Hu T, Wang S, Huang L, Wang J, Shi D, Li Y, et al. A clinical-radiomics nomogram for the preoperative prediction of lung metastasis in colorectal cancer patients with indeterminate pulmonary nodules. Eur Radiol. (2018) 29:439–49. doi: 10.1007/s00330-018-5539-3

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Ying GS, Maguire M, Quinn G, Kulp MT, Cyert L, Vision in preschoolers study G. ROC analysis of the accuracy of noncycloplegic retinoscopy, retinomax autorefractor, and suresight vision screener for preschool vision screening. Invest Ophthalmol Vis Sci. (2011) 52:9658–64. doi: 10.1167/iovs.11-8559

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Pomerri F, Pucciarelli S, Maretto I, Perrone E, Pintacuda G, Lonardi S, et al. Significance of pulmonary nodules in patients with colorectal cancer. Eur Radiol. (2012) 22:1680–6. doi: 10.1007/s00330-012-2431-4

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Quyn AJ, Matthews A, Daniel T, Amin AI, Yalamarthi S. The clinical significance of radiologically detected indeterminate pulmonary nodules in colorectal cancer. Colorectal Dis. (2012) 14:828–31. doi: 10.1111/j.1463-1318.2011.02722.x

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Nordholm-Carstensen A, Jorgensen LN, Wille-Jorgensen PA, Hansen H, Harling H. Indeterminate pulmonary nodules in colorectal-cancer: do radiologists agree? Ann Surg Oncol. (2015) 22:543–9. doi: 10.1245/s10434-014-4063-1

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Kim CH, Huh JW, Kim HR, Kim YJ. Indeterminate pulmonary nodules in colorectal cancer: follow-up guidelines based on a risk predictive model. Ann Surg. (2015) 261:1145–52. doi: 10.1097/SLA.0000000000000853

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Xiang W, Xing Y, Jiang S, Chen G, Mao H, Labh K, et al. Morphological factors differentiating between early lung adenocarcinomas appearing as pure ground-glass nodules measuring ≤10 mm on thin-section computed tomography. Cancer Imaging. (2014) 14:33. doi: 10.1186/s40644-014-0033-x

CrossRef Full Text | Google Scholar

27. Lee SM, Park CM, Goo JM, Lee HJ, Wi JY, Kang CH. Invasive pulmonary adenocarcinomas versus preinvasive lesions appearing as ground-glass nodules: differentiation by using CT features. Radiology. (2013) 268:265–73. doi: 10.1148/radiol.13120949

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Yang Q, Wang Y, Ban X, Wu J, Rong D, Zhao Q, et al. Prediction of pulmonary metastasis in pulmonary nodules (ions appearing as ground-glass nodules: differentiation by using CT features. screening. formRadiol Med. (2017) 122:837–49. doi: 10.1007/s11547-017-0790-2

CrossRef Full Text | Google Scholar

29. Suo S, Cheng J, Cao M, Lu Q, Yin Y, Xu J, et al. Assessment of heterogeneity difference between edge and core by using texture analysis: differentiation of malignant from inflammatory pulmonary nodules and masses. Acad Radiol. (2016) 23:1115–22. doi: 10.1016/j.acra.2016.04.009

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Ferreira Junior JR, Koenigkam-Santos M, Cipriano FEG, Fabro AT, Azevedo-Marques PM. RadiomicS-based features for pattern recognition of lung cancer histopathology and metastases. Comput Methods Programs Biomed. (2018) 159:23–30. doi: 10.1016/j.cmpb.2018.02.015

CrossRef Full Text | Google Scholar

31. Ng F, Kozarski R, Ganeshan B, Goh V. Assessment of tumor heterogeneity by CT texture analysis: can the largest crosS-sectional area be used as an alternative to whole tumor analysis? Eur J Radiol. (2013) 82:342–8. doi: 10.1016/j.ejrad.2012.10.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: colorectal cancer, pulmonary metastases, histogram, morphological, morphological features, nomogram

Citation: Hu T, Wang S, E X, Yuan Y, Huang L, Wang J, Shi D, Li Y, Peng W and Tong T (2019) CT Morphological Features Integrated With Whole-Lesion Histogram Parameters to Predict Lung Metastasis for Colorectal Cancer Patients With Pulmonary Nodules. Front. Oncol. 9:1241. doi: 10.3389/fonc.2019.01241

Received: 07 July 2019; Accepted: 29 October 2019;
Published: 19 November 2019.

Edited by:

Chunxiao Guo, University of Texas MD Anderson Cancer Center, United States

Reviewed by:

DingWei Ye, Fudan University Shanghai Cancer Center, China
Emily Ann Thompson, University of Texas MD Anderson Cancer Center, United States

Copyright © 2019 Hu, Wang, E, Yuan, Huang, Wang, Shi, Li, Peng and Tong. 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: Tong Tong, t983352@126.com

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.