Ronald A. Karwoski

Automated quantification of radiological patterns predicts survival in idiopathic pulmonary fibrosis

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. Short-term quantified CT changes are predictive of survival in IPF http://ow.ly/qmbjd

Quantitative computed tomography imaging of interstitial lung diseases.

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.

System for the analysis of whole-bone strength from helical CT images

Osteoporosis affects an estimated 44 million Americans. This condition results from bone loss, but the measured change in bone mass does not fully account for the marked decrease in whole-bone structural integrity seen in osteoporosis. In order to study structural changes in bone mineral distribution due to normal ageing and osteoporosis, we have developed a method for progressive analysis of whole-bone mechanical integrity from helical CT images. The system provides rapid semi-automated alignment of femur and vertebrae volume images into standard anatomic reference planes, and calculates bone mineral density in any selected 3D sections of bone. Mineral density measures are obtained using both full-width-half-max contours and threshold-derived masks, and are obtained for cortical bone and trabecular bone separately. Biomechanical properties of the bone cross-section are also assessed, including the 2-D bending moment of the cortical bone region and the integrated flexural rigidity of the cortical region or whole-bone region in arbitrary planes. This method facilitates progressive refinement of the analysis protocol by separating the labor-intensive alignment and landmark selection process from the analysis process. As the analysis protocol evolves to include new measures, previously analyzed images can be automatically reanalyzed, using the image regions originally specified. Initial results show inverse correlation of indices of biomechanical bone strength with age, greater loss of bone strength in the lumbar spine than in the femoral neck, and more trabecular than cortical bone loss at both sites.

Mortality prediction in idiopathic pulmonary fibrosis: evaluation of computer-based CT analysis with conventional severity measures

Computer-based computed tomography (CT) analysis can provide objective quantitation of disease in idiopathic pulmonary fibrosis (IPF). A computer algorithm, CALIPER, was compared with conventional CT and pulmonary function measures of disease severity for mortality prediction. CT and pulmonary function variables (forced expiratory volume in 1 s, forced vital capacity, diffusion capacity of the lung for carbon monoxide, transfer coefficient of the lung for carbon monoxide and composite physiologic index (CPI)) of 283 consecutive patients with a multidisciplinary diagnosis of IPF were evaluated against mortality. Visual and CALIPER CT features included total extent of interstitial lung disease, honeycombing, reticular pattern, ground glass opacities and emphysema. In addition, CALIPER scored pulmonary vessel volume (PVV) while traction bronchiectasis and consolidation were only scored visually. A combination of mortality predictors was compared with the Gender, Age, Physiology model. On univariate analyses, all visual and CALIPER-derived interstitial features and functional indices were predictive of mortality to a 0.01 level of significance. On multivariate analysis, visual CT parameters were discarded. Independent predictors of mortality were CPI (hazard ratio (95% CI) 1.05 (1.02–1.07), p<0.001) and two CALIPER parameters: PVV (1.23 (1.08–1.40), p=0.001) and honeycombing (1.18 (1.06–1.32), p=0.002). A three-group staging system derived from this model was powerfully predictive of mortality (2.23 (1.85–2.69), p<0.0001). CALIPER-derived parameters, in particular PVV, are more accurate prognostically than traditional visual CT scores. Quantitative tools such as CALIPER have the potential to improve staging systems in IPF. CALIPER-derived parameters such as pulmonary vessel volume are more accurate prognostically than visual CT scores http://ow.ly/2b5G304exlA

Automated Quantitative Computed Tomography Versus Visual Computed Tomography Scoring in Idiopathic Pulmonary Fibrosis: Validation Against Pulmonary Function.

Purpose: The aim of the study was to determine whether a novel computed tomography (CT) postprocessing software technique (CALIPER) is superior to visual CT scoring as judged by functional correlations in idiopathic pulmonary fibrosis (IPF). Materials and Methods: A total of 283 consecutive patients with IPF had CT parenchymal patterns evaluated quantitatively with CALIPER and by visual scoring. These 2 techniques were evaluated against: forced expiratory volume in 1 second (FEV1), forced vital capacity (FVC), diffusing capacity for carbon monoxide (DLco), carbon monoxide transfer coefficient (Kco), and a composite physiological index (CPI), with regard to extent of interstitial lung disease (ILD), extent of emphysema, and pulmonary vascular abnormalities. Results: CALIPER-derived estimates of ILD extent demonstrated stronger univariate correlations than visual scores for most pulmonary function tests (PFTs): (FEV1: CALIPER R2=0.29, visual R2=0.18; FVC: CALIPER R2=0.41, visual R2=0.27; DLco: CALIPER R2=0.31, visual R2=0.35; CPI: CALIPER R2=0.48, visual R2=0.44). Correlations between CT measures of emphysema extent and PFTs were weak and did not differ significantly between CALIPER and visual scoring. Intriguingly, the pulmonary vessel volume provided similar correlations to total ILD extent scored by CALIPER for FVC, DLco, and CPI (FVC: R2=0.45; DLco: R2=0.34; CPI: R2=0.53). Conclusions: CALIPER was superior to visual scoring as validated by functional correlations with PFTs. The pulmonary vessel volume, a novel CALIPER CT parameter with no visual scoring equivalent, has the potential to be a CT feature in the assessment of patients with IPF and requires further exploration.

Noninvasive Characterization of the Histopathologic Features of Pulmonary Nodules of the Lung Adenocarcinoma Spectrum using Computer-Aided Nodule Assessment and Risk Yield (CANARY)—A Pilot Study

Introduction: Pulmonary nodules of the adenocarcinoma spectrum are characterized by distinctive morphological and radiologic features and variable prognosis. Noninvasive high-resolution computed tomography–based risk stratification tools are needed to individualize their management. Methods: Radiologic measurements of histopathologic tissue invasion were developed in a training set of 54 pulmonary nodules of the adenocarcinoma spectrum and validated in 86 consecutively resected nodules. Nodules were isolated and characterized by computer-aided analysis, and data were analyzed by Spearman correlation, sensitivity, and specificity and the positive and negative predictive values. Results: Computer-aided nodule assessment and risk yield (CANARY) can noninvasively characterize pulmonary nodules of the adenocarcinoma spectrum. Unsupervised clustering analysis of high-resolution computed tomography data identified nine unique exemplars representing the basic radiologic building blocks of these lesions. The exemplar distribution within each nodule correlated well with the proportion of histologic tissue invasion, Spearman R = 0.87, p < 0.0001 and 0.89 and p < 0.0001 for the training and the validation set, respectively. Clustering of the exemplars in three-dimensional space corresponding to tissue invasion and lepidic growth was used to develop a CANARY decision algorithm that successfully categorized these pulmonary nodules as “aggressive” (invasive adenocarcinoma) or “indolent” (adenocarcinoma in situ and minimally invasive adenocarcinoma). Sensitivity, specificity, positive predictive value, and negative predictive value of this approach for the detection of aggressive lesions were 95.4, 96.8, 95.4, and 96.8%, respectively, in the training set and 98.7, 63.6, 94.9, and 87.5%, respectively, in the validation set. Conclusion: CANARY represents a promising tool to noninvasively risk stratify pulmonary nodules of the adenocarcinoma spectrum.

Noninvasive Computed Tomography–based Risk Stratification of Lung Adenocarcinomas in the National Lung Screening Trial

RATIONALE Screening for lung cancer using low-dose computed tomography (CT) reduces lung cancer mortality. However, in addition to a high rate of benign nodules, lung cancer screening detects a large number of indolent cancers that generally belong to the adenocarcinoma spectrum. Individualized management of screen-detected adenocarcinomas would be facilitated by noninvasive risk stratification. OBJECTIVES To validate that Computer-Aided Nodule Assessment and Risk Yield (CANARY), a novel image analysis software, successfully risk stratifies screen-detected lung adenocarcinomas based on clinical disease outcomes. METHODS We identified retrospective 294 eligible patients diagnosed with lung adenocarcinoma spectrum lesions in the low-dose CT arm of the National Lung Screening Trial. The last low-dose CT scan before the diagnosis of lung adenocarcinoma was analyzed using CANARY blinded to clinical data. Based on their parametric CANARY signatures, all the lung adenocarcinoma nodules were risk stratified into three groups. CANARY risk groups were compared using survival analysis for progression-free survival. MEASUREMENTS AND MAIN RESULTS A total of 294 patients were included in the analysis. Kaplan-Meier analysis of all the 294 adenocarcinoma nodules stratified into the Good, Intermediate, and Poor CANARY risk groups yielded distinct progression-free survival curves (P < 0.0001). This observation was confirmed in the unadjusted and adjusted (age, sex, race, and smoking status) progression-free survival analysis of all stage I cases. CONCLUSIONS CANARY allows the noninvasive risk stratification of lung adenocarcinomas into three groups with distinct post-treatment progression-free survival. Our results suggest that CANARY could ultimately facilitate individualized management of incidentally or screen-detected lung adenocarcinomas.

Noninvasive risk stratification of lung adenocarcinoma using quantitative computed tomography

Introduction: Lung cancer remains the leading cause of cancer-related deaths in the United States and worldwide. Adenocarcinoma is the most common type of lung cancer and encompasses lesions with widely variable clinical outcomes. In the absence of noninvasive risk stratification, individualized patient management remains challenging. Consequently a subgroup of pulmonary nodules of the lung adenocarcinoma spectrum is likely treated more aggressively than necessary. Methods: Consecutive patients with surgically resected pulmonary nodules of the lung adenocarcinoma spectrum (lesion size ⩽3 cm, 2006–2009) and available presurgical high-resolution computed tomography (HRCT) imaging were identified at Mayo Clinic Rochester. All cases were classified using an unbiased Computer-Aided Nodule Assessment and Risk Yield (CANARY) approach based on the quantification of presurgical HRCT characteristics. CANARY-based classification was independently correlated to postsurgical progression-free survival. Results: CANARY analysis of 264 consecutive patients identified three distinct subgroups. Independent comparisons of 5-year disease-free survival (DFS) between these subgroups demonstrated statistically significant differences in 5-year DFS, 100%, 72.7%, and 51.4%, respectively (p = 0.0005). Conclusions: Noninvasive CANARY-based risk stratification identifies subgroups of patients with pulmonary nodules of the adenocarcinoma spectrum characterized by distinct clinical outcomes. This technique may ultimately improve the current expert opinion-based approach to the management of these lesions by facilitating individualized patient management.

Evaluation of computer-based computer tomography stratification against outcome models in connective tissue disease-related interstitial lung disease: a patient outcome study.

BackgroundTo evaluate computer-based computer tomography (CT) analysis (CALIPER) against visual CT scoring and pulmonary function tests (PFTs) when predicting mortality in patients with connective tissue disease-related interstitial lung disease (CTD-ILD). To identify outcome differences between distinct CTD-ILD groups derived following automated stratification of CALIPER variables.MethodsA total of 203 consecutive patients with assorted CTD-ILDs had CT parenchymal patterns evaluated by CALIPER and visual CT scoring: honeycombing, reticular pattern, ground glass opacities, pulmonary vessel volume, emphysema, and traction bronchiectasis. CT scores were evaluated against pulmonary function tests: forced vital capacity, diffusing capacity for carbon monoxide, carbon monoxide transfer coefficient, and composite physiologic index for mortality analysis. Automated stratification of CALIPER-CT variables was evaluated in place of and alongside forced vital capacity and diffusing capacity for carbon monoxide in the ILD gender, age physiology (ILD-GAP) model using receiver operating characteristic curve analysis.ResultsCox regression analyses identified four independent predictors of mortality: patient age (P < 0.0001), smoking history (P = 0.0003), carbon monoxide transfer coefficient (P = 0.003), and pulmonary vessel volume (P < 0.0001). Automated stratification of CALIPER variables identified three morphologically distinct groups which were stronger predictors of mortality than all CT and functional indices. The Stratified-CT model substituted automated stratified groups for functional indices in the ILD-GAP model and maintained model strength (area under curve (AUC) = 0.74, P < 0.0001), ILD-GAP (AUC = 0.72, P < 0.0001). Combining automated stratified groups with the ILD-GAP model (stratified CT-GAP model) strengthened predictions of 1- and 2-year mortality: ILD-GAP (AUC = 0.87 and 0.86, respectively); stratified CT-GAP (AUC = 0.89 and 0.88, respectively).ConclusionsCALIPER-derived pulmonary vessel volume is an independent predictor of mortality across all CTD-ILD patients. Furthermore, automated stratification of CALIPER CT variables represents a novel method of prognostication at least as robust as PFTs in CTD-ILD patients.

Comparison of manual and semiautomated techniques for analyzing gastric volumes with MRI in humans

Gastric emptying, accommodation, and motility can be quantified with magnetic resonance imaging (MRI). The first step in image analysis entails segmenting the stomach from surrounding structures, usually by a time-consuming manual process. We have developed a semiautomated process to segment and measure gastric volumes with MRI. Gastric images were acquired with a three-dimensional gradient echo MRI sequence at 5, 10, 20, and 30 min after ingestion of a liquid nutrient (Ensure, 296 ml) labeled with gadolinium in 20 healthy volunteers and 29 patients with dyspeptic symptoms. The agreement between gastric volumes measured by manual segmentation and our new semiautomated algorithm was assessed with Lins concordance correlation coefficient (CCC) and the Bland Altman test. At 5 min after a meal, food volumes measured by manual (352 ± 4 ml) and semiautomated (346 ± 4 ml) techniques were correlated {CCC[95% confidence interval (CI)] 0.70 (0.52, 0.81)}; air volumes measured by manual (88 ± 6 ml) and semiautomated (84 ± 6 ml) techniques were also correlated [CCC (95% CI) 0.89 (0.82, 0.94)]. Findings were similar at subsequent time points. The Bland Altman test was not significant. The time required for semiautomated segmentation ranged from an average of 204 s for the 5-min images to 233 s for the 20-min images. These times were appreciably smaller than the typical times of many tens of minutes, even hours, required for manual segmentation. To conclude, a semiautomated process can measure gastric food and air volume using MRI with comparable accuracy and far better efficiency than a manual process.