Noninvasive pulmonary nodule elastometry by CT and deformable image registration

Assessing four-dimensional CT stress maps derived from patient-specific biomechanical models of the lung with pulmonary function test data in lung cancer patients

OBJECTIVE

This study aims to retrospectively compare the stress map of the lung with pulmonary function test (PFT) results in lung cancer patients and to evaluate the potential of the stress map as an imaging biomarker for chronic obstructive pulmonary disease (COPD).

METHODS

25 lung cancer patients with pre-treatment four-dimensional CT (4DCT) and PFT data were retrospectively analysed. PFT metrics were used to diagnose obstructive lung disease. For each patient, forced expiratory volume in 1 s (FEV1 % predicted) and the ratio of FEV1 and forced vital capacity (FEV1/FVC) were recorded. 4DCT and biomechanical model-deformable image registration (BM-DIR) were used to obtain the lung stress map. The relationship between the mean of the total lung stress and PFT data was evaluated, and the COPD classification grade was also evaluated.

RESULTS

The mean values of the total lung stress and FEV1 % predicted showed a significant strong correlation [R = 0.833, (p < 0.001)]. The mean values and FEV1/FVC showed a significant strong correlation [R = 0.805, (p < 0.001)]. For the total lung stress, the area under the curve and the optimal cut-off value were 0.94 and 510.8 Pa for the classification of normal or abnormal lung function, respectively.

CONCLUSION

This study has demonstrated the potential of lung stress maps based on BM-DIR to accurately assess lung function by comparing them with PFT data.

ADVANCES IN KNOWLEDGE

The derivation of stress map directly from 4DCT is novel method. The BM-DIR-based lung stress map can provide an accurate assessment of lung function.

A Feasibility Study of Single-inhalation, Single-energy Xenon-enhanced CT for High-resolution Imaging of Regional Lung Ventilation in Humans

RATIONALE AND OBJECTIVES

The objective of this study was to assess the feasibility of single-inhalation xenon-enhanced computed tomography (XeCT) to provide clinically practical, high-resolution pulmonary ventilation imaging to clinics with access to only a single-energy computed tomography scanner, and to reduce the subject's overall exposure to xenon by utilizing a higher (70%) concentration for a much shorter time than has been employed in prior studies.

MATERIALS AND METHODS

We conducted an institutional review board-approved prospective feasibility study of XeCT for 15 patients undergoing thoracic radiotherapy. For XeCT, we acquired two breath-hold single-energy computed tomography images of the entire lung with a single inhalation each of 100% oxygen and a mixture of 70% xenon and 30% oxygen, respectively. A video biofeedback system for coached patient breathing was used to achieve reproducible breath holds. We assessed the technical success of XeCT acquisition and side effects. We then used deformable image registration to align the breath-hold images with each other to accurately subtract them, producing a map of lung xenon distribution. Additionally, we acquired ventilation single-photon emission computed tomography-computed tomography (V-SPECT-CT) images for 11 of the 15 patients. For a comparative analysis, we partitioned each lung into 12 sectors, calculated the xenon concentration from the Hounsfield unit enhancement in each sector, and then correlated this with the corresponding V-SPECT-CT counts.

RESULTS

XeCT scans were tolerated well overall, with a mild (grade 1) dizziness as the only side effect in 5 of the 15 patients. Technical failures in five patients occurred because of inaccurate breathing synchronization with xenon gas delivery, leaving seven patients analyzable for XeCT and single-photon emission computed tomography correlation. Sector-wise correlations were strong (Spearman coefficient >0.75, Pearson coefficient >0.65, P value <.002) for two patients for whom ventilation deficits were visibly pronounced in both scans. Correlations were nonsignificant for the remaining five who had more homogeneous XeCT ventilation maps, as well as strong V-SPECT-CT imaging artifacts attributable to airway deposition of the aerosolized imaging agent. Qualitatively, XeCT demonstrated higher resolution and no central airway deposition artifacts compared to V-SPECT-CT.

CONCLUSIONS

In this pilot study, single-breath XeCT ventilation imaging was generally feasible for patients undergoing thoracic radiotherapy, using an imaging protocol that is clinically practical and potentially widely available. In the future, the xenon delivery failures can be addressed by straightforward technical improvements to the patient biofeedback coaching system.

Feasibility of deriving a novel imaging biomarker based on patient-specific lung elasticity for characterizing the degree of COPD in lung SBRT patients

OBJECTIVE:

Lung tissue elasticity is an effective spatial representation for Chronic Obstructive Pulmonary Disease phenotypes and pathophysiology. We investigated a novel imaging biomarker based on the voxel-by-voxel distribution of lung tissue elasticity. Our approach combines imaging and biomechanical modeling to characterize tissue elasticity.

METHODS:

We acquired 4DCT images for 13 lung cancer patients with known COPD diagnoses based on GOLD 2017 criteria. Deformation vector fields (DVFs) from the deformable registration of end-inhalation and end-exhalation breathing phases were taken to be the ground-truth. A linear elastic biomechanical model was assembled from end-exhalation datasets with a density-guided initial elasticity distribution. The elasticity estimation was formulated as an iterative process, where the elasticity was optimized based on its ability to reconstruct the ground-truth. An imaging biomarker (denoted YM_1-3) derived from the optimized elasticity distribution, was compared with the current gold standard, RA₉₅₀ using confusion matrix and area under the receiver operating characteristic (AUROC) curve analysis.

RESULTS:

The estimated elasticity had 90 % accuracy when representing the ground-truth DVFs. The YM_1-3 biomarker had higher diagnostic accuracy (86% vs 71 %), higher sensitivity (0.875 vs 0.5), and a higher AUROC curve (0.917 vs 0.875) as compared to RA₉₅₀. Along with acting as an effective spatial indicator of lung pathophysiology, the YM_1-3 biomarker also proved to be a better indicator for diagnostic purposes than RA₉₅₀.

CONCLUSIONS:

Overall, the results suggest that, as a biomarker, lung tissue elasticity will lead to new end points for clinical trials and new targeted treatment for COPD subgroups.

ADVANCES IN KNOWLEDGE:

The derivation of elasticity information directly from 4DCT imaging data is a novel method for performing lung elastography. The work demonstrates the need for a mechanics-based biomarker for representing lung pathophysiology.

Prognostic factors associated with the accuracy of deformable image registration in lung cancer patients treated with stereotactic body radiotherapy

We evaluated the accuracy of an in-house program in lung stereotactic body radiation therapy (SBRT) cancer patients, and explored the prognostic factors associated with the accuracy of deformable image registrations (DIRs). The accuracy of the 3 programs which implement the free-form deformation and the B-spline algorithm was compared regarding the structures on 4-dimensional computed tomography (4DCT) image datasets between the peak-inhale and peak-exhale phases. The dice similarity coefficient (DSC) and normalized DSC (NDSC) were measured for the gross tumor volumes from 19 lung SBRT patients. We evaluated the accuracy of DIR using gross tumor volume, magnitude of displacement from 0% phase to 50% phase, whole lung volume in the 50% phase image, and status of tumor pleural attachment. The median NDSC values using the NiftyReg, MIM Maestro and Velocity AI programs were 1.027, 1.005, and 0.946, respectively, indicating that NiftyReg and MIM Maestro programs had similar accuracy with an uncertainty of < 1 mm. Larger uncertainty of 1 to 2 mm was observed using the Velocity AI program. The NiftyReg and the MIM programs provided higher NDSC values than the median values when the gross tumor volume was attached to the pleura (p <0.05). However, it showed a different trend in using the Velocity AI program. All software programs provided unexpected results, and there is a possibility that such results would reduce the accuracy of 4D treatment planning and adaptive radiotherapy. The unexpected results may be because the tumors are surrounded by other tissues, and there are differences regarding the region of interest for rigid and nonrigid registration. Furthermore, our results indicated that the pleural attachment status might be an important predictor of DIR accuracy for thoracic images, indicating that there is a potentially large dose distribution discrepancy concerning 4D treatment planning and adaptive radiotherapy.

Dynamic CT imaging of volumetric changes in pulmonary nodules correlates with physical measurements of stiffness

BACKGROUND AND PURPOSE

A major challenge in CT screening for lung cancer is limited specificity when distinguishing between malignant and non-malignant pulmonary nodules (PN). Malignant nodules have different mechanical properties and tissue characteristics ('stiffness') from non-malignant nodules. This study seeks to improve CT specificity by demonstrating in rats that measurements of volumetric ratios in PNs with varying composition can be determined by respiratory-gated dynamic CT imaging and that these ratios correlate with direct physical measurements of PN stiffness.

METHODS AND MATERIALS

Respiratory-gated MicroCT images acquired at extreme tidal volumes of 9 rats with PNs from talc, matrigel and A549 human lung carcinoma were analyzed and their volumetric ratios (δ) derived. PN stiffness was determined by measuring the Young's modulus using atomic force microscopy (AFM) for each nodule excised immediately after MicroCT imaging.

RESULTS

There was significant correlation (p=0.0002) between PN volumetric ratios determined by respiratory-gated CT imaging and the physical stiffness of the PNs determined from AFM measurements.

CONCLUSION

We demonstrated proof of concept that PN volume changes measured non-invasively correlate with direct physical measurements of stiffness. These results may translate clinically into a means of improving the specificity of CT screening for lung cancer and/or improving individual prognostic assessments based on lung tumor stiffness.