Variability of Image Features Computed from Conventional and Respiratory-Gated PET/CT Images of Lung Cancer

AbstractRadiomics is being explored for potential applications in radiation therapy. How various imaging protocols affect quantitative image features is currently a highly active area of research. To assess the variability of image features derived from conventional [three-dimensional (3D)] and respiratory-gated (RG) positron emission tomography (PET)/computed tomography (CT) images of lung cancer patients, image features were computed from 23 lung cancer patients. Both protocols for each patient were acquired during the same imaging session. PET tumor volumes were segmented using an adaptive technique which accounted for background. CT tumor volumes were delineated with a commercial segmentation tool. Using RG PET images, the tumor center of mass motion, length, and rotation were calculated. Fifty-six image features were extracted from all images consisting of shape descriptors, first-order features, and second-order texture features. Overall, 26.6% and 26.2% of total features demonstrated less than 5% difference between 3D and RG protocols for CT and PET, respectively. Between 10 RG phases in PET, 53.4% of features demonstrated percent differences less than 5%. The features with least variability for PET were sphericity, spherical disproportion, entropy (first and second order), sum entropy, information measure of correlation 2, Short Run Emphasis (SRE), Long Run Emphasis (LRE), and Run Percentage (RPC); and those for CT were minimum intensity, mean intensity, Root Mean Square (RMS), Short Run Emphasis (SRE), and RPC. Quantitative analysis using a 3D acquisition versus RG acquisition (to reduce the effects of motion) provided notably different image feature values. This study suggests that the variability between 3D and RG features is mainly due to the impact of respiratory motion

Adaptive hypofractionted and stereotactic body radiotherapy for lung tumors with real-time MRI guidance

The treatment of central and ultracentral lung tumors with radiotherapy remains an ongoing clinical challenge. The risk of Grade 5 toxicity with ablative radiotherapy doses to these high-risk regions is significant as shown in recent prospective studies. Magnetic resonance (MR) image-guided adaptive radiotherapy (MRgART) is a new technology and may allow the delivery of ablative radiotherapy to these high-risk regions safely. MRgART is able to achieve this by utilizing small treatment margins, real-time gating/tracking and on-table plan adaptation to maintain dose to the tumor but limit dose to critical structures. The process of MRgART is complex and has nuances and challenges for the treatment of lung tumors. We outline the critical steps needed for appropriate delivery of MRgART for lung tumors safely and effectively

A proliferation saturation index to predict radiation response and personalize radiotherapy fractionation

BACKGROUND: Although altered protocols that challenge conventional radiation fractionation have been tested in prospective clinical trials, we still have limited understanding of how to select the most appropriate fractionation schedule for individual patients. Currently, the prescription of definitive radiotherapy is based on the primary site and stage, without regard to patient-specific tumor or host factors that may influence outcome. We hypothesize that the proportion of radiosensitive proliferating cells is dependent on the saturation of the tumor carrying capacity. This may serve as a prognostic factor for personalized radiotherapy (RT) fractionation. METHODS: We introduce a proliferation saturation index (PSI), which is defined as the ratio of tumor volume to the host-influenced tumor carrying capacity. Carrying capacity is as a conceptual measure of the maximum volume that can be supported by the current tumor environment including oxygen and nutrient availability, immune surveillance and acidity. PSI is estimated from two temporally separated routine pre-radiotherapy computed tomography scans and a deterministic logistic tumor growth model. We introduce the patient-specific pre-treatment PSI into a model of tumor growth and radiotherapy response, and fit the model to retrospective data of four non-small cell lung cancer patients treated exclusively with standard fractionation. We then simulate both a clinical trial hyperfractionation protocol and daily fractionations, with equal biologically effective dose, to compare tumor volume reduction as a function of pretreatment PSI. RESULTS: With tumor doubling time and radiosensitivity assumed constant across patients, a patient-specific pretreatment PSI is sufficient to fit individual patient response data (R(2) = 0.98). PSI varies greatly between patients (coefficient of variation >128 %) and correlates inversely with radiotherapy response. For this study, our simulations suggest that only patients with intermediate PSI (0.45–0.9) are likely to truly benefit from hyperfractionation. For up to 20 % uncertainties in tumor growth rate, radiosensitivity, and noise in radiological data, the absolute estimation error of pretreatment PSI is <10 % for more than 75 % of patients. CONCLUSIONS: Routine radiological images can be used to calculate individual PSI, which may serve as a prognostic factor for radiation response. This provides a new paradigm and rationale to select personalized RT dose-fractionation

Gonadal-sparing total body irradiation with the use of helical tomotherapy for nonmalignant indications

Background: The aim was to demonstrate the feasibility and technique of gonadal sparing total body irradiation (TBI) with helical tomotherapy. Total body irradiation is a common part of the conditioning regimen prior to allogeneic stem cell transplantation. Shielding or dose-reduction to the gonads is often desired to preserve fertility, particularly in young patients undergoing transplant for non-malignant indications. Helical tomotherapy (HT) has been shown to be superior to traditional TBI delivery for organ at risk (OAR) doses and dose homogeneity. Materials and methods: We present two representative cases (one male and one female) to illustrate the feasibility of this technique, each of whom received 3Gy in a single fraction prior to allogeneic stem cell transplant for benign indications. The planning target volume (PTV) included the whole body with a subtraction of OARs including the lungs, heart, and brain (each contracted by 1cm) as well as the gonads (testicles expanded by 5 cm and ovaries expanded by 0.5 cm). Results: For the male patient we achieved a homogeneity index of 1.35 with a maximum and median planned dose to the testes of 0.53 Gy and 0.35 Gy, respectively. In-vivo dosimetry demonstrated an actual received dose of 0.48 Gy. For the female patient we achieved a homogeneity index of 1.13 with a maximum and median planned dose to the ovaries of 1.66 Gy and 0.86 Gy, respectively. Conclusion: Gonadal sparing TBI is feasible and deliverable using HT in patients with non-malignant diseases requiring TBI as part of a pre-stem cell transplant conditioning regimen

Assessment of the Dependence of Ventilation Image Calculation from 4D-CT on Deformation and Ventilation Algorithms

Ventilation imaging using 4D-CT is a convenient and cost effective functional imaging methodology which might be of value in radiotherapy treatment planning to spare functional lung volumes. To calculate ventilation imaging from 4D-CT we must use deformable image registration (DIR). This study validates the DIR methods and investigates the dependence of calculated ventilation on DIR methods and ventilation algorithms. The first hypothesis is if ventilation algorithms are robust then they will be insensitive to the precise DIR used provided the DIR is accurate. The second hypothesis is that the change in Houndsfield Unit (HU) method is less dependent on the DIR used and depends more on the CT image quality due to the inherent noise of HUs in normal CT imaging. DIR of the normal end expiration and inspiration phases of the 4D-CT images was used to correlate the voxels between the two respiratory phases. All DIR algorithms were validated using a 4D pixel-based and point-validated breathing thorax model, consisting of a 4D-CT image data set along with associated landmarks. Three different DIR algorithms, Optical Flow (OF), Diffeomorphic Demons (DD) and Diffeomorphic Morphons (DM), were retrospectively applied to the same group of 10 esophagus and 10 lung cancer cases all of which had associated 4D-CT image sets that encompassed the entire lung volume. Three different ventilation calculation algorithms were compared (Jacobian, ΔV, and HU) using the Dice similarity coefficient comparison. In the validation of the DIR algorithms, the average target registration errors with one standard deviation for the DIR algorithms were 1.6 ± 0.7 mm, maximum 3.1 mm for OF, 1.3 ± 0.6 mm, maximum 3.3 mm for DM, 1.3 ± 0.6 mm, maximum 2.8 mm for DD, indicating registration errors were within 2 voxels. Dependence of ventilation images on the DIR was greater for the ΔV and the Jacobian methods than for the HU method. The Dice similarity coefficient for 20% of low ventilation volume for ΔV was 0.33 ± 0.03 between OF and DM, 0.44 ± 0.05 between OF and DD and 0.51 ± 0.04 between DM and DD. The similarity comparisons for Jacobian was 0.32 ± 0.03, 0.44 ± 0.05 and 0.51 ± 0.04 respectively, and for HU 0.53 ± 0.03, 0.56 ± 0.03 and 0.76 ± 0.04 respectively. Dependence of ventilation images on the ventilation method used showed good agreement between the ΔV and Jacobian methods but differences between these two and the HU method were significantly greater. Dice similarity coefficient for using OF as DIR was 0.86 ± 0.01 between ΔV and Jacobian, 0.28 ± 0.04 between ΔV and HU and 0.28 ± 0.04 between Jacobian and HU respectively. When using DM or DD as DIR, similar values were obtained when comparing the different ventilation calculation methods. The similarity values for 20% of the high ventilation volume were close to those found for the 20% low ventilation volume. Mean target registration error for all three DIR methods was within one voxel suggesting that the registration done by either of the methods is quite accurate. Ventilation calculation from 4D-CT demonstrates some degree of dependency on the DIR algorithm employed. Similarities between ΔV and Jacobian are higher than between ΔV and HU and Jacobian and HU. This shows that ΔV and Jacobian are very similar, but HU is a very different ventilation calculation method

Validation of Three Deformable Image Registration Algorithms for the Thorax

Deformable image registration (DIR) has been proposed for lung ventilation calculation using 4D CT. Spatial accuracy of DIR can be evaluated using expert landmark correspondences. Additionally, image differences between the deformed and the target images give a degree of accuracy of DIR algorithms for the same image modality registration. DIR of the normal end-expiration (50%), end- inspiration (0%), midexpiration (30%), and midinspiration image (70%) phases of the 4D CT images was used to correlate the voxels between the respiratory phases. Three DIR algorithms, optical flow (OF), diffeomorphic morphons (DM), and diffeomorphic demons (DD) were validated using a 4D thorax model, consisting of a 4D CT image dataset, along with associated landmarks delineated by a radiologist. Image differences between the deformed and the target images were used to evaluate the degree of registration accuracy of the three DIR algorithms. In the validation of the DIR algorithms, the average target registration error (TRE) for normal end-expiration-to-end-inspiration registration with one standard deviation (SD) for the DIR algorithms was (maximum 3.1 mm) for OF, (maximum 3.3 mm) for DM, and (maximum 3.3 mm) for DD, indicating registration errors were within two voxels. As a reference, the median value of TRE between 0 and 50% phases with rigid registration only was 5.0 mm with one SD of 2.5 mm and the maximum value of 12.0 mm. For the OF algorithm, 81% of voxels were within a difference of 50 HU, and 93% of the voxels were within 100 HU. For the DM algorithm, 69% of voxels were within 50 HU, and 87% within 100 HU. For the DD algorithm, 71% of the voxels were within 50 HU, and 87% within a difference of 100 HU. These data suggest that the three DIR methods perform accurate registrations in the thorax region. The mean TRE for all three DIR methods was less than two voxels suggesting that the registration performed by all methods are equally accurate in the thorax

Validation of Three Deformable Image Registration Algorithms for the Thorax

Deformable image registration (DIR) has been proposed for lung ventilation calculation using 4D CT. Spatial accuracy of DIR can be evaluated using expert landmark correspondences. Additionally, image differences between the deformed and the target images give a degree of accuracy of DIR algorithms for the same image modality registration. DIR of the normal end-expiration (50%), end- inspiration (0%), midexpiration (30%), and midinspiration image (70%) phases of the 4D CT images was used to correlate the voxels between the respiratory phases. Three DIR algorithms, optical flow (OF), diffeomorphic morphons (DM), and diffeomorphic demons (DD) were validated using a 4D thorax model, consisting of a 4D CT image dataset, along with associated landmarks delineated by a radiologist. Image differences between the deformed and the target images were used to evaluate the degree of registration accuracy of the three DIR algorithms. In the validation of the DIR algorithms, the average target registration error (TRE) for normal end-expiration-to-end-inspiration registration with one standard deviation (SD) for the DIR algorithms was (maximum 3.1 mm) for OF, (maximum 3.3 mm) for DM, and (maximum 3.3 mm) for DD, indicating registration errors were within two voxels. As a reference, the median value of TRE between 0 and 50% phases with rigid registration only was 5.0 mm with one SD of 2.5 mm and the maximum value of 12.0 mm. For the OF algorithm, 81% of voxels were within a difference of 50 HU, and 93% of the voxels were within 100 HU. For the DM algorithm, 69% of voxels were within 50 HU, and 87% within 100 HU. For the DD algorithm, 71% of the voxels were within 50 HU, and 87% within a difference of 100 HU. These data suggest that the three DIR methods perform accurate registrations in the thorax region. The mean TRE for all three DIR methods was less than two voxels suggesting that the registration performed by all methods are equally accurate in the thorax

Monte Carlo study of radiation dose enhancement by gadolinium in megavoltage and high dose rate radiotherapy.

MRI is often used in tumor localization for radiotherapy treatment planning, with gadolinium (Gd)-containing materials often introduced as a contrast agent. Motexafin gadolinium is a novel radiosensitizer currently being studied in clinical trials. The nanoparticle technologies can target tumors with high concentration of high-Z materials. This Monte Carlo study is the first detailed quantitative investigation of high-Z material Gd-induced dose enhancement in megavoltage external beam photon therapy. BEAMnrc, a radiotherapy Monte Carlo simulation package, was used to calculate dose enhancement as a function of Gd concentration. Published phase space files for the TrueBeam flattening filter free (FFF) and conventional flattened 6MV photon beams were used. High dose rate (HDR) brachytherapy with Ir-192 source was also investigated as a reference. The energy spectra difference caused a dose enhancement difference between the two beams. Since the Ir-192 photons have lower energy yet, the photoelectric effect in the presence of Gd leads to even higher dose enhancement in HDR. At depth of 1.8 cm, the percent mean dose enhancement for the FFF beam was 0.38±0.12, 1.39±0.21, 2.51±0.34, 3.59±0.26, and 4.59±0.34 for Gd concentrations of 1, 5, 10, 15, and 20 mg/mL, respectively. The corresponding values for the flattened beam were 0.09±0.14, 0.50±0.28, 1.19±0.29, 1.68±0.39, and 2.34±0.24. For Ir-192 with direct contact, the enhanced were 0.50±0.14, 2.79±0.17, 5.49±0.12, 8.19±0.14, and 10.80±0.13. Gd-containing materials used in MRI as contrast agents can also potentially serve as radiosensitizers in radiotherapy. This study demonstrates that Gd can be used to enhance radiation dose in target volumes not only in HDR brachytherapy, but also in 6 MV FFF external beam radiotherapy, but higher than the currently used clinical concentration (>5 mg/mL) would be needed

Voxel Size and Gray Level Normalization of CT Radiomic Features in Lung Cancer

Radiomic features are potential imaging biomarkers for therapy response assessment in oncology. However, the robustness of features with respect to imaging parameters is not well established. Previously identified potential imaging biomarkers were found to be intrinsically dependent on voxel size and number of gray levels (GLs) in a recent texture phantom investigation. Here, we validate the voxel size and GL in-phantom normalizations in lung tumors. Eighteen patients with non-small cell lung cancer of varying tumor volumes were analyzed. To compare with patient data, phantom scans were acquired on eight different scanners. Twenty four previously identified features were extracted from lung tumors. The Spearman rank (rs) and interclass correlation coefficient (ICC) were used as metrics. Eight out of 10 features showed high (rs \u3e 0.9) and low (rs \u3c 0.5) correlations with number of voxels before and after normalizations, respectively. Likewise, texture features were unstable (ICC \u3c 0.6) and highly stable (ICC \u3e 0.8) before and after GL normalizations, respectively. We conclude that voxel size and GL normalizations derived from a texture phantom study also apply to lung tumors. This study highlights the importance and utility of investigating the robustness of radiomic features with respect to CT imaging parameters in radiomic phantoms

Normalization of ventilation data from 4D-CT to facilitate comparison between datasets acquired at different times.

PURPOSE: The 4D-CT data used for comparing a patient's ventilation distributions before and after lung radiotherapy are acquired at different times. As a result, an additional variable--the tidal volume (TV)--can alter the results. Therefore, in this paper we propose to normalize the ventilation to the same TV to eliminate that uncertainty. METHODS: Absolute ventilation (AV) data were generated for 6 stereotactic body radiation therapy (SBRT) cases before and after treatment, using the direct geometric algorithm and diffeomorphic morphons deformable image registration (DIR). Each pair of AV distributions was converted to TV-normalized, percentile ventilation (PV) and low-dose well-ventilated-normalized ventilation (LDWV) distributions. The ventilation change was calculated in various dose regions based on the treatment plans using the DIR-registered before and after treatment data sets. The ventilation change based on TV-normalized ventilation was compared with the AV as well as the data normalized by PV and LDWV. RESULTS: AV change may be misleading when the TV differs before and after treatment, which was found to be up to 6.7%. All three normalization methods produced a similar trend in ventilation change: the higher the dose to a region of lung, the greater the degradation in ventilation. In low dose regions (<5 Gy), ventilation appears relatively improved after treatment due to the relative nature of the normalized ventilation. However, the LDWV may not be reliable when the ventilation in the low-dose regions varies. PV exhibited a similar ventilation change trend compared to the TV-normalized in all cases. However, by definition, the ventilation distribution in the PV is significantly different from the original distribution. CONCLUSION: Normalizing ventilation distributions by the TV is a simple and reliable method for evaluation of ventilation changes