Optimization of the clinical management of screen-detected lung nodules is needed to avoid unnecessary diagnostic interventions. Herein we demonstrate the potential value of a novel radiomics-based ...approach for the classification of screen-detected indeterminate nodules.
Independent quantitative variables assessing various radiologic nodule features such as sphericity, flatness, elongation, spiculation, lobulation and curvature were developed from the NLST dataset using 726 indeterminate nodules (all ≥ 7 mm, benign, n = 318 and malignant, n = 408). Multivariate analysis was performed using least absolute shrinkage and selection operator (LASSO) method for variable selection and regularization in order to enhance the prediction accuracy and interpretability of the multivariate model. The bootstrapping method was then applied for the internal validation and the optimism-corrected AUC was reported for the final model.
Eight of the originally considered 57 quantitative radiologic features were selected by LASSO multivariate modeling. These 8 features include variables capturing Location: vertical location (Offset carina centroid z), Size: volume estimate (Minimum enclosing brick), Shape: flatness, Density: texture analysis (Score Indicative of Lesion/Lung Aggression/Abnormality (SILA) texture), and surface characteristics: surface complexity (Maximum shape index and Average shape index), and estimates of surface curvature (Average positive mean curvature and Minimum mean curvature), all with P<0.01. The optimism-corrected AUC for these 8 features is 0.939.
Our novel radiomic LDCT-based approach for indeterminate screen-detected nodule characterization appears extremely promising however independent external validation is needed.
Accurate assessment of prognosis in idiopathic pulmonary fibrosis remains elusive due to significant individual radiological and physiological variability. We hypothesised that short-term ...radiological changes may be predictive of survival. We explored the use of CALIPER (Computer-Aided Lung Informatics for Pathology Evaluation and Rating), a novel software tool developed by the Biomedical Imaging Resource Laboratory at the Mayo Clinic Rochester (Rochester, MN, USA) for the analysis and quantification of parenchymal lung abnormalities on high-resolution computed tomography. We assessed baseline and follow-up (time-points 1 and 2, respectively) high-resolution computed tomography scans in 55 selected idiopathic pulmonary fibrosis patients and correlated CALIPER-quantified measurements with expert radiologists' assessments and clinical outcomes. Findings of interval change (mean 289 days) in volume of reticular densities (hazard ratio 1.91, p=0.006), total volume of interstitial abnormalities (hazard ratio 1.70, p=0.003) and per cent total interstitial abnormalities (hazard ratio 1.52, p=0.017) as quantified by CALIPER were predictive of survival after a median follow-up of 2.4 years. Radiologist interpretation of short-term global interstitial lung disease progression, but not specific radiological features, was also predictive of mortality. These data demonstrate the feasibility of quantifying interval short-term changes on high-resolution computed tomography and their possible use as independent predictors of survival in idiopathic pulmonary fibrosis.
The mechanisms underlying airflow obstruction in COPD cannot be distinguished by standard spirometry. We ascertain whether mathematical modeling of airway biomechanical properties, as assessed from ...spirometry, could provide estimates of emphysema presence and severity, as quantified by computed tomography (CT) metrics and CT-based radiomics.
We quantified presence and severity of emphysema by standard CT metrics (VIDA) and co-registration analysis (ImbioLDA) of inspiratory-expiratory CT in 194 COPD patients who underwent pulmonary function testing. According to percentages of low attenuation area below - 950 Hounsfield Units (%LAA
) patients were classified as having no emphysema (NE) with %LAA
< 6, moderate emphysema (ME) with %LAA
≥ 6 and < 14, and severe emphysema (SE) with %LAA
≥ 14. We also obtained stratified clusters of emphysema CT features by an automated unsupervised radiomics approach (CALIPER). An emphysema severity index (ESI), derived from mathematical modeling of the maximum expiratory flow-volume curve descending limb, was compared with pulmonary function data and the three CT classifications of emphysema presence and severity as derived from CT metrics and radiomics.
ESI mean values and pulmonary function data differed significantly in the subgroups with different emphysema degree classified by VIDA, ImbioLDA and CALIPER (p < 0.001 by ANOVA). ESI differentiated NE from ME/SE CT-classified patients (sensitivity 0.80, specificity 0.85, AUC 0.86) and SE from ME CT-classified patients (sensitivity 0.82, specificity 0.87, AUC 0.88).
Presence and severity of emphysema in patients with COPD, as quantified by CT metrics and radiomics can be estimated by mathematical modeling of airway function as derived from standard spirometry.
Diffuse parenchymal lung diseases (DPLDs) are characterized by widespread pathological changes within the pulmonary tissue that impair the elasticity and gas exchange properties of the lungs. ...Clinical-radiological diagnosis of these diseases remains challenging and their clinical course is characterized by variable disease progression. These challenges have hindered the introduction of robust objective biomarkers for patient-specific prediction based on specific phenotypes in clinical practice for patients with DPLD. Therefore, strategies facilitating individualized clinical management, staging and identification of specific phenotypes linked to clinical disease outcomes or therapeutic responses are urgently needed. A classification schema consistently reflecting the radiological, clinical (lung function and clinical outcomes) and pathological features of a disease represents a critical need in modern pulmonary medicine. Herein, we report a quantitative stratification paradigm to identify subsets of DPLD patients with characteristic radiologic patterns in an unsupervised manner and demonstrate significant correlation of these self-organized disease groups with clinically accepted surrogate endpoints. The proposed consistent and reproducible technique could potentially transform diagnostic staging, clinical management and prognostication of DPLD patients as well as facilitate patient selection for clinical trials beyond the ability of current radiological tools. In addition, the sequential quantitative stratification of the type and extent of parenchymal process may allow standardized and objective monitoring of disease, early assessment of treatment response and mortality prediction for DPLD patients.
To investigate the correlations between densitometric and CALIPER-derived indices of pulmonary emphysema and their change in the short-term period for groups of patients with different smoking ...habits.
This retrospective study included 284 subjects from the ITALUNG trial (198 men and 86 women; mean age±sd, 60±4) who underwent low-dose chest CT at baseline and 2-year follow-up. Subjects were divided in four groups (persistent smokers, restarters, quitters, and former smokers) according to their smoking habit at baseline and follow-up. Densitometric and texture analyses were performed, using CALIPER software. A correlation analysis was conducted between CALIPER-derived low-attenuation-areas (LAAs) and densitometric indices, including the 15th percentile of the whole-lung attenuation histogram (Perc15) and the relative areas below -950HU (RA950). Densitometric indices and LAAs were evaluated at baseline and variation assessed longitudinally with comparisons between groups with different smoking habit. Further analysis of parenchymal changes per pulmonary zone was performed.
LAAs were strongly correlated with Perc15 (rs=0.81; p<0.001) and RA950 (rs=0.905; p<0.001). At baseline, the group of smokers showed higher Perc15, lower RA950, lower LAAs (particularly mild sub-class of LAAs) than the group of ex-smokers (p<0.001). At 2-year follow-up, densitometric indices and LAAs increased in persistent smokers, former smokers and quitters (p<0.05). The progression was larger and statistically more significant in quitters (p<0.001).
CALIPER texture analysis provides an objective measure comparable to traditional density/histogram features to assess the lung parenchymal changes in relation to different smoking habit.
PURPOSE:High-resolution chest computed tomography (HRCT) is essential in the characterization of interstitial lung disease. The HRCT features of some diseases can be diagnostic. Longitudinal ...monitoring with HRCT can assess progression of interstitial lung disease; however, subtle changes in the volume and character of abnormalities can be difficult to assess. Accuracy of diagnosis can be dependent on expertise and experience of the radiologist, pathologist, or clinician. Quantitative analysis of thoracic HRCT has the potential to determine the extent of disease reproducibly, classify the types of abnormalities, and automate the diagnostic process.
MATERIALS AND METHODS:Novel software that utilizes histogram signatures to characterize pulmonary parenchyma was used to analyze chest HRCT data, including retrospective processing of clinical CT scans and research data from the Lung Tissue Research Consortium. Additional information including physiological, pathologic, and semiquantitative radiologist assessment was available to allow comparison of quantitative results, with visual estimates of the disease, physiological parameters, and measures of disease outcome.
RESULTS:Quantitative analysis results were provided in regional volumetric quantities for statistical analysis and a graphical representation. These results suggest that quantitative HRCT analysis can serve as a biomarker with physiological, pathologic, and prognostic significance.
CONCLUSIONS:It is likely that quantitative analysis of HRCT can be used in clinical practice as a means to aid in identifying a probable diagnosis, stratifying prognosis in early disease, and consistently determining progression of the disease or response to therapy. Further optimization of quantitative techniques and longitudinal analysis of well-characterized subjects would be helpful in validating these methods.
Lung adenocarcinoma (ADC), the most common lung cancer type, is recognized increasingly as a disease spectrum. To guide individualized patient care, a non-invasive means of distinguishing indolent ...from aggressive ADC subtypes is needed urgently. Computer-Aided Nodule Assessment and Risk Yield (CANARY) is a novel computed tomography (CT) tool that characterizes early ADCs by detecting nine distinct CT voxel classes, representing a spectrum of lepidic to invasive growth, within an ADC. CANARY characterization has been shown to correlate with ADC histology and patient outcomes. This study evaluated the inter-observer variability of CANARY analysis. Three novice observers segmented and analyzed independently 95 biopsy-confirmed lung ADCs from Vanderbilt University Medical Center/Nashville Veterans Administration Tennessee Valley Healthcare system (VUMC/TVHS) and the Mayo Clinic (Mayo). Inter-observer variability was measured using intra-class correlation coefficient (ICC). The average ICC for all CANARY classes was 0.828 (95% CI 0.76, 0.895) for the VUMC/TVHS cohort, and 0.852 (95% CI 0.804, 0.901) for the Mayo cohort. The most invasive voxel classes had the highest ICC values. To determine whether nodule size influenced inter-observer variability, an additional cohort of 49 sub-centimeter nodules from Mayo were also segmented by three observers, with similar ICC results. Our study demonstrates that CANARY ADC classification between novice CANARY users has an acceptably low degree of variability, and supports the further development of CANARY for clinical application.
Most computed tomography (CT)-detected lung cancers are adenocarcinomas (ACs), representing lesions with variable tissue invasion, aggressiveness, and clinical outcome. Visual radiologic ...characterization of AC pulmonary nodules is both inconsistent and inadequate to confidently predict histopathologic classification or prognosis. Comprehensive pathologic interpretation requires full nodule resection. We have described a computerized scoring system for AC detected on CT scans that can noninvasively estimate the degree of histologic invasion and simultaneously predict patient survival.
The Computer-Aided Nodule Assessment and Risk Yield has been validated to characterize CT-detected nodules across the spectrum of AC. With the use of unsupervised clustering, nine natural exemplars were identified as basic radiographic features of AC nodules. We now introduce the Score Indicative of Lung Cancer Aggression (SILA), which is a cumulative aggregate of normalized distributions of ordered Computer-Aided Nodule Assessment and Risk Yield exemplars. The SILA values for each of 237 unique nodules in AC were compared with the histopathologically defined maximum linear extent of tumor invasion. With use of the SILA, Kaplan-Meier survival and Cox proportionality analysis were performed on patients with stage I AC, who comprised a subset of our cohort.
The SILA discriminated between indolent and invasive AC (p < 0.0001). In addition, prediction of linear extent of histopathologic tumor invasion was possible. In stage I AC, three separate SILA prognosis groups were identified: indolent, intermediate, and poor, with 5-year survival rates of 100%, 79%, 58%, respectively. Cox proportionality hazard modeling predicted a 50% increase in mortality, for a 0.1 unit increase in the SILA over a median follow-up time of 3.6 years (p < 0.0002).
The SILA is a computer-based analytic measure allowing noninvasive approximation of histologic invasion and prediction of patient survival in CT-detected AC nodules.
Implementation of low-dose chest computed tomography (CT) lung cancer screening and the ever-increasing use of cross-sectional imaging are resulting in the identification of many screen- and ...incidentally detected indeterminate pulmonary nodules. While the management of nodules with low or high pre-test probability of malignancy is relatively straightforward, those with intermediate pre-test probability commonly require advanced imaging or biopsy. Noninvasive risk stratification tools are highly desirable.
We previously developed the BRODERS classifier (Benign
aggRessive nODule Evaluation using Radiomic Stratification), a conventional predictive radiomic model based on eight imaging features capturing nodule location, shape, size, texture and surface characteristics. Herein we report its external validation using a dataset of incidentally identified lung nodules (Vanderbilt University Lung Nodule Registry) in comparison to the Brock model. Area under the curve (AUC), as well as sensitivity, specificity, negative and positive predictive values were calculated.
For the entire Vanderbilt validation set (n=170, 54% malignant), the AUC was 0.87 (95% CI 0.81-0.92) for the Brock model and 0.90 (95% CI 0.85-0.94) for the BRODERS model. Using the optimal cut-off determined by Youden's index, the sensitivity was 92.3%, the specificity was 62.0%, the positive (PPV) and negative predictive values (NPV) were 73.7% and 87.5%, respectively. For nodules with intermediate pre-test probability of malignancy, Brock score of 5-65% (n=97), the sensitivity and specificity were 94% and 46%, respectively, the PPV was 78.4% and the NPV was 79.2%.
The BRODERS radiomic predictive model performs well on an independent dataset and may facilitate the management of indeterminate pulmonary nodules.