Inversion-recovery ultrashort-echo-time (IR-UTE) MRI-based detection of radiation dose heterogeneity in gynecologic cancer patients treated with HDR brachytherapy.

PURPOSE: To evaluate the relationship between delivered radiation (RT) and post-RT inversion-recovery ultrashort-echo-time (IR-UTE) MRI signal-intensity (SI) in gynecologic cancer patients treated with high-dose-rate (HDR) brachytherapy (BT). METHODS: Seven patients underwent whole-pelvis RT (WPRT) followed by BT to the high-risk clinical target volume (HR-CTV). MR images were acquired at three time-points; pre-RT, post-WPRT/pre-BT, and 3-6 months post-BT. Diffuse-fibrosis (FDiffuse) was imaged with a non-contrast dual-echo IR (inversion time [TI] = 60 ms) UTE research application, with image-subtraction of the later echo, only retaining the ultrashort-echo SI. Dense-fibrosis (FDense) imaging utilized single-echo Late-Gadolinium-Enhanced IR-UTE, acquired ∼ 15 min post-Gadavist injection. Resulting FDiffuse and FDense SI were normalized to the corresponding gluteal-muscle SI. Images were deformably registered between time-points based on normal tissue anatomy. The remnant tumor at both time-points was segmented using multi-parametric MRI. Contours corresponding to the 50%, 100%, 150%, and 200% isodose lines (IDLs) of the prescription BT-dose were created. Mean FDiffuse and FDense SI within (i) each IDL contour and (ii) the remnant tumor were calculated. Post-BT FDiffuse and FDense SI were correlated with prescribed BT-dose. To determine the relationship between BT-dose and IR-UTE SI, the differences in the post-BT FDense across IDLs was determined using paired t-tests with Bonferroni correction. RESULTS: FDense was higher in regions of higher dose for 6/7 patients, with mean ± SD values of 357 ± 103% and 331 ± 97% (p = .03) in the 100% and 50% IDL, respectively. FDense was higher in regions of higher dose in the responsive regions with mean ± SD values of 380 ± 122% and 356 ± 135% (p = .03) in the 150% and 50% IDL, respectively. Within the segmented remnant tumor, an increase in prescribed dose correlated with an increase in FDense post-BT (n = 5, r = .89, p = .04). Post-BT FDiffuse inversely correlated (n = 7, r = -.83, p = .02) with prescribed BT-dose within the 100% IDL. CONCLUSIONS: Results suggest that FDense SI 3-6 months post-BT is a sensitive measure of tissue response to heterogeneous BT radiation-dose. Future studies will validate whether FDiffuse and FDense are accurate biomarkers of fibrotic radiation response.

The Peaks and Valleys of Photon Versus Proton Spatially Fractionated Radiotherapy.

Spatially-fractionated radiotherapy (SFRT) delivers high doses to small areas of tumor while sparing adjacent tissue, including intervening disease. In this review, we explore the evolution of SFRT technological advances, contrasting approaches with photon and proton beam radiotherapy. We discuss unique dosimetric considerations and physical properties of SFRT, as well as review the preclinical literature that provides an emerging understanding of biological mechanisms. We emphasize crucial areas of future study and highlight clinical trials that are underway to assess SFRT's safety and efficacy, with a focus on immunotherapeutic synergies. The review concludes with practical considerations for SFRT's clinical application, advocating for strategies that leverage its unique dosimetric and biological properties for improved patient outcomes.

Overview and Recommendations for Prospective Multi-institutional Spatially Fractionated Radiation Therapy Clinical Trials.

PURPOSE: The highly heterogeneous dose delivery of spatially fractionated radiation therapy (SFRT) is a profound departure from standard radiation planning and reporting approaches. Early SFRT studies have shown excellent clinical outcomes. However, prospective multi-institutional clinical trials of SFRT are still lacking. This NRG Oncology/American Association of Physicists in Medicine working group consensus aimed to develop recommendations on dosimetric planning, delivery, and SFRT dose reporting to address this current obstacle toward the design of SFRT clinical trials. METHODS AND MATERIALS: Working groups consisting of radiation oncologists, radiobiologists, and medical physicists with expertise in SFRT were formed in NRG Oncology and the American Association of Physicists in Medicine to investigate the needs and barriers in SFRT clinical trials. RESULTS: Upon reviewing the SFRT technologies and methods, this group identified challenges in several areas, including the availability of SFRT, the lack of treatment planning system support for SFRT, the lack of guidance in the physics and dosimetry of SFRT, the approximated radiobiological modeling of SFRT, and the prescription and combination of SFRT with conventional radiation therapy. CONCLUSIONS: Recognizing these challenges, the group further recommended several areas of improvement for the application of SFRT in cancer treatment, including the creation of clinical practice guidance documents, the improvement of treatment planning system support, the generation of treatment planning and dosimetric index reporting templates, and the development of better radiobiological models through preclinical studies and through conducting multi-institution clinical trials.

Clinical Characterization of a Table Mounted Range Shifter Board for Synchrotron-Based Intensity Modulated Proton Therapy for Pediatric Craniospinal Irradiation.

Purpose: To report our design, manufacturing, commissioning and initial clinical experience with a table-mounted range shifter board (RSB) intended to replace the machine-mounted range shifter (MRS) in a synchrotron-based pencil beam scanning (PBS) system to reduce penumbra and normal tissue dose for image-guided pediatric craniospinal irradiation (CSI). Methods: A custom RSB was designed and manufactured from a 3.5 cm thick slab of polymethyl methacrylate (PMMA) to be placed directly under patients, on top of our existing couch top. The relative linear stopping power (RLSP) of the RSB was measured using a multi-layer ionization chamber, and output constancy was measured using an ion chamber. End-to-end tests were performed using the MRS and RSB approaches using an anthropomorphic phantom and radiochromic film measurements. Cone beam CT (CBCT) and 2D planar kV X-ray image quality were compared with and without the RSB present using image quality phantoms. CSI plans were produced using MRS and RSB approaches for two retrospective pediatric patients, and the resultant normal tissue doses were compared. Results: The RLSP of the RSB was found to be 1.163 and provided computed penumbra of 6.9 mm in the phantom compared to 11.8 mm using the MRS. Phantom measurements using the RSB demonstrated errors in output constancy, range, and penumbra of 0.3%, -0.8%, and 0.6 mm, respectively. The RSB reduced mean kidney and lung dose compared to the MRS by 57.7% and 46.3%, respectively. The RSB decreased mean CBCT image intensities by 86.8 HU but did not significantly impact CBCT or kV spatial resolution providing acceptable image quality for patient setup. Conclusions: A custom RSB for pediatric proton CSI was designed, manufactured, modeled in our TPS, and found to significantly reduce lateral proton beam penumbra compared to a standard MRS while maintaining CBCT and kV image-quality and is in routine use at our center.

Treatment Planning of Bulky Tumors Using Pencil Beam Scanning Proton GRID Therapy.

Purpose: To compare spatially fractionated radiation therapy (GRID) treatment planning techniques using proton pencil-beam-scanning (PBS) and photon therapy. Materials and Methods: PBS and volumetric modulated arc therapy (VMAT) GRID plans were retrospectively generated for 5 patients with bulky tumors. GRID targets were arranged along the long axis of the gross tumor, spaced 2 and 3 cm apart, and treated with a prescription of 18 Gy. PBS plans used 2- to 3-beam multiple-field optimization with robustness evaluation. Dosimetric parameters including peak-to-edge ratio (PEDR), ratio of dose to 90% of the valley to dose to 10% of the peak VPDR(D90/D10), and volume of normal tissue receiving at least 5 Gy (V5) and 10 Gy (V10) were calculated. The peak-to-valley dose ratio (PVDR), VPDR(D90/D10), and organ-at-risk doses were prospectively assessed in 2 patients undergoing PBS-GRID with pretreatment quality assurance computed tomography (QACT) scans. Results: PBS and VMAT GRID plans were generated for 5 patients with bulky tumors. Gross tumor volume values ranged from 826 to 1468 cm3. Peak-to-edge ratio for PBS was higher than for VMAT for both spacing scenarios (2-cm spacing, P = .02; 3-cm spacing, P = .01). VPDR(D90/D10) for PBS was higher than for VMAT (2-cm spacing, P = .004; 3-cm spacing, P = .002). Normal tissue V5 was lower for PBS than for VMAT (2-cm spacing, P = .03; 3-cm spacing, P = .02). Normal tissue mean dose was lower with PBS than with VMAT (2-cm spacing, P = .03; 3-cm spacing, P = .02). Two patients treated using PBS GRID and assessed with pretreatment QACT scans demonstrated robust PVDR, VPDR(D90/D10), and organs-at-risk doses. Conclusions: The PEDR was significantly higher for PBS than VMAT plans, indicating lower target edge dose. Normal tissue mean dose was significantly lower with PBS than VMAT. PBS GRID may result in lower normal tissue dose compared with VMAT plans, allowing for further dose escalation in patients with bulky disease.

Planning and Treatment Recommendations for Breast Proton Therapy From a Single Center's Experience.

Purpose: Proton therapy use for breast cancer has grown due to advantages in coverage and potentially reduced late toxicities compared with conventional radiation therapy. We aimed to provide recommendations for robustness criteria, daily imaging, and quality assurance computed tomography (QA CT) frequency for these patients. Methods and Materials: All patients treated for localized breast cancer at the Johns Hopkins Proton Center between November 2019 and February 2022 were eligible for inclusion. Daily shift information was extracted and examined through control charts. If an adaptive plan was used, the time to replan was recorded. Three and 5 mm setup uncertainty was used to calculate robustness. Robust evaluation of QA CTs was compared with initial robustness range for breast/chest wall and lymph node target coverage. Results: Sixty-six patients were included: 19 with intact breast, 25 with non-reconstructed chest wall, and 22 with chest wall plus expanders or implants. Sixteen percent, 13%, and 41% of breast, chest wall, and expander/implant patients had a replan. Only patients with expanders or implants required 2 adaptive plans. Daily shift data showed large variation and did not correlate with plan adaptation. Patients without adaptive plans had QA CTs with dose-volume histogram metrics within robustness more frequently than those with adaptive plans. Using 3 mm robustness for patients who did not require an adaptive plan, 91% to 100% of patients had QA CTs within robustness, while 55% to 60% of patients with an adaptive plan had QA CTs within robustness for the axilla, internal mammary nodes, and supraclavicular nodes. Five millimeter setup uncertainty did not significantly improve this. Conclusions: We recommend using daily cone beam CT because of the large variation in daily setup with 3 mm setup uncertainty in robustness analysis. If daily cone beam CT imaging is not available, then larger setup uncertainty should be used. Two QA CTs should be conducted during treatment if the patient has expanders or implants; otherwise, one QA CT is sufficient.

Functional Lung Avoidance for Individualized Radiation Therapy: Results of a Double-Masked, Randomized Controlled Trial.

PURPOSE: To determine whether functional lung avoidance based on 3He magnetic resonance imaging (MRI) improves quality of life (QOL) for patients undergoing concurrent chemoradiotherapy (CCRT) for advanced non-small cell lung cancer. METHODS AND MATERIALS: Patients with stage III non-small cell lung cancer (or oligometastatic disease treated with curative intent) undergoing CCRT with at least a 10 pack-year smoking history were eligible. Patients underwent pretreatment 3He MRI to measure lung ventilation and had 2 radiation therapy (RT) plans created before randomization: a standard plan, which did not make use of the 3He MRI, and an avoidance plan, preferentially sparing well-ventilated lung. All participants were masked to assignment except the physicist responsible for exporting the selected plan. The primary end point was patient-reported QOL measured at 3-months post-RT by the FACT-L lung cancer subscale (LCS); secondary end points included other QOL metrics, toxicity, and survival outcomes. Target accrual was 64. RESULTS: Twenty-seven patients were randomized before the trial was closed due to slower-than-expected accrual, with 11 randomized to the standard arm and 16 to the avoidance arm. Baseline patient characteristics were well-balanced. At 3 months post-RT, the mean ± SD LCS scores were 17.4 ± 2.8 versus 17.3 ± 6.1 for the standard and avoidance arms, respectively (P = .485). A clinically meaningful, prespecified decline of ≥3 points in the LCS score was observed in 50% (4/8) in the standard arm and 33% (4/12) in the avoidance arm (P = .648). Two patients in each arm developed grade ≥2 radiation pneumonitis, with no grade ≥4 toxicities. CONCLUSIONS: Although this trial did not reach full accrual, QOL scores were very similar between arms. Due to the scarcity of 3He MRI, other, more commonly available methods to measure functional lung, such as 4-dimensional computed tomography ventilation mapping, may be considered in the assessment of functional lung avoidance RT, and a larger, multicenter approach would be needed to accrue sufficient patients to test such approaches.

Evaluating Proton Dose and Associated Range Uncertainty Using Daily Cone-Beam CT.

Purpose: This study aimed to quantitatively evaluate the range uncertainties that arise from daily cone-beam CT (CBCT) images for proton dose calculation compared to CT using a measurement-based technique. Methods: For head and thorax phantoms, wedge-shaped intensity-modulated proton therapy (IMPT) treatment plans were created such that the gradient of the wedge intersected and was measured with a 2D ion chamber array. The measured 2D dose distributions were compared with 2D dose planes extracted from the dose distributions using the IMPT plan calculated on CT and CBCT. Treatment plans of a thymoma cancer patient treated with breath-hold (BH) IMPT were recalculated on 28 CBCTs and 9 CTs, and the resulting dose distributions were compared. Results: The range uncertainties for the head phantom were determined to be 1.2% with CBCT, compared to 0.5% for CT, whereas the range uncertainties for the thorax phantom were 2.1% with CBCT, compared to 0.8% for CT. The doses calculated on CBCT and CT were similar with similar anatomy changes. For the thymoma patient, the primary source of anatomy change was the BH uncertainty, which could be up to 8 mm in the superior-inferior (SI) direction. Conclusion: We developed a measurement-based range uncertainty evaluation method with high sensitivity and used it to validate the accuracy of CBCT-based range and dose calculation. Our study demonstrated that the CBCT-based dose calculation could be used for daily dose validation in selected proton patients.

Dosimetric evaluation of cone-beam CT-based synthetic CTs in pediatric patients undergoing intensity-modulated proton therapy.

PURPOSE: To evaluate dosimetric changes detected using synthetic computed tomography (sCT) derived from online cone-beam CTs (CBCT) in pediatric patients treated using intensity-modulated proton therapy (IMPT). METHODS: Ten pediatric patients undergoing IMPT and aligned daily using proton gantry-mounted CBCT were identified for retrospective analysis with treated anatomical sites fully encompassed in the CBCT field of view. Dates were identified when the patient received both a CBCT and a quality assurance CT (qCT) for routine dosimetric evaluation. sCTs were generated based on a deformable registration between the initial plan CT (pCT) and CBCT. The clinical IMPT plans were re-computed on the same day qCT and sCT, and dosimetric changes due to tissue change or response from the initial plan were computed using each image. Linear regression analysis was performed to determine the correlation between dosimetric changes detected using the qCT and the sCT. Gamma analysis was also used to compare the dose distributions computed on the qCT and sCT. RESULTS: The correlation coefficients (p-values) between qCTs and sCTs for changes detected in target coverage, overall maximum dose, and organ at risk dose were 0.97 (< .001), 0.84 (.002) and 0.91 (< .001), respectively. Mean ± SD gamma pass rates of the sCT-based dose compared to the qCT-based dose at 3%/3 mm, 3%/2 mm, and 2%/2 mm criteria were 96.5%±4.5%, 93.2%±6.3%, and 91.3%±7.8%, respectively. Pass rates tended to be lower for targets near lung. CONCLUSION: While insufficient for re-planning, sCTs provide approximate dosimetry without administering additional imaging dose in pediatric patients undergoing IMPT. Dosimetric changes detected using sCTs are correlated with changes detected using clinically-standard qCTs; however, residual differences in dosimetry remain a limitation. Further improvements in sCT image quality may both improve online dosimetric evaluation and reduce imaging dose for pediatric patients by reducing the need for routine qCTs.

Radiomic modeling to predict risk of vertebral compression fracture after stereotactic body radiation therapy for spinal metastases.

OBJECTIVE: In the treatment of spinal metastases with stereotactic body radiation therapy (SBRT), vertebral compression fracture (VCF) is a common and potentially morbid complication. Better methods to identify patients at high risk of radiation-induced VCF are needed to evaluate prophylactic measures. Radiomic features from pretreatment imaging may be employed to more accurately predict VCF. The objective of this study was to develop and evaluate a machine learning model based on clinical characteristics and radiomic features from pretreatment imaging to predict the risk of VCF after SBRT for spinal metastases. METHODS: Vertebral levels C2 through L5 containing metastases treated with SBRT were included if they were naive to prior surgery or radiation therapy, target delineation was based on consensus guidelines, and 1-year follow-up data were available. Clinical features, including characteristics of the patient, disease, and treatment, were obtained from chart review. Radiomic features were extracted from the planning target volume (PTV) on pretreatment CT and T1-weighted MRI. Clinical and radiomic features selected by least absolute shrinkage and selection operator (LASSO) regression were included in random forest classification models, which were trained to predict VCF within 1 year after SBRT. Model performance was assessed with leave-one-out cross-validation. RESULTS: Within 1 year after SBRT, 15 of 95 vertebral levels included in the analysis demonstrated new or progressive VCF. Selected clinical features included BMI, performance status, total prescription dose, dose to 99% of the PTV, lumbar location, and 2 components of the Spine Instability Neoplastic Score (SINS): lytic tumor character and spinal misalignment. Selected radiomic features included 5 features from CT and 3 features from MRI. The best-performing classification model, derived from a combination of selected clinical and radiomic features, demonstrated a sensitivity of 0.844, specificity of 0.800, and area under the receiver operating characteristic (ROC) curve (AUC) of 0.878. This model was significantly more accurate than alternative models derived from only selected clinical features (AUC = 0.795, p = 0.048) or only components of the SINS (AUC = 0.579, p < 0.0001). CONCLUSIONS: In the treatment of spinal metastases with SBRT, a machine learning model incorporating both clinical features and radiomic features from pretreatment imaging predicted VCF at 1 year after SBRT with excellent sensitivity and specificity, outperforming models developed from clinical features or components of the SINS alone. If validated, these findings may allow more judicious selection of patients for prophylactic interventions.

Radiation Versus Immune Checkpoint Inhibitor Associated Pneumonitis: Distinct Radiologic Morphologies.

BACKGROUND: Patients with non-small cell lung cancer may develop pneumonitis after thoracic radiotherapy (RT) and immune checkpoint inhibitors (ICIs). We hypothesized that distinct morphologic features are associated with different pneumonitis etiologies. MATERIALS AND METHODS: We systematically compared computed tomography (CT) features of RT- versus ICI-pneumonitis. Clinical and imaging features were tested for association with pneumonitis severity. Lastly, we constructed an exploratory radiomics-based machine learning (ML) model to discern pneumonitis etiology. RESULTS: Between 2009 and 2019, 82 patients developed pneumonitis: 29 after thoracic RT, 23 after ICI, and 30 after RT + ICI. Fifty patients had grade 2 pneumonitis, 22 grade 3, and 7 grade 4. ICI-pneumonitis was more likely bilateral (65% vs. 28%; p = .01) and involved more lobes (66% vs. 45% involving at least three lobes) and was less likely to have sharp border (17% vs. 59%; p = .004) compared with RT-pneumonitis. Pneumonitis morphology after RT + ICI was heterogeneous, with 47% bilateral, 37% involving at least three lobes, and 40% sharp borders. Among all patients, risk factors for severe pneumonitis included poor performance status, smoking history, worse lung function, and bilateral and multifocal involvement on CT. An ML model based on seven radiomic features alone could distinguish ICI- from RT-pneumonitis with an area under the receiver-operating curve of 0.76 and identified the predominant etiology after RT + ICI concordant with multidisciplinary consensus. CONCLUSION: RT- and ICI-pneumonitis exhibit distinct spatial features on CT. Bilateral and multifocal lung involvement is associated with severe pneumonitis. Integrating these morphologic features in the clinical management of patients who develop pneumonitis after RT and ICIs may improve treatment decision-making. IMPLICATIONS FOR PRACTICE: Patients with non-small cell lung cancer often receive thoracic radiation and immune checkpoint inhibitors (ICIs), both of which can cause pneumonitis. This study identified similarities and differences in pneumonitis morphology on computed tomography (CT) scans among pneumonitis due to radiotherapy (RT) alone, ICI alone, and the combination of both. Patients who have bilateral CT changes involving at least three lobes are more likely to have ICI-pneumonitis, whereas those with unilateral CT changes with sharp borders are more likely to have radiation pneumonitis. After RT and/or ICI, severe pneumonitis is associated with bilateral and multifocal CT changes. These results can help guide clinicians in triaging patients who develop pneumonitis after radiation and during ICI treatment.

Predicting acute radiation induced xerostomia in head and neck Cancer using MR and CT Radiomics of parotid and submandibular glands.

PURPOSE: To analyze baseline CT/MR-based image features of salivary glands to predict radiation-induced xerostomia 3-months after head-and-neck cancer (HNC) radiotherapy. METHODS: A retrospective analysis was performed on 266 HNC patients who were treated using radiotherapy at our institution between 2009 and 2018. CT and T1 post-contrast MR images along with NCI-CTCAE xerostomia grade (3-month follow-up) were prospectively collected at our institution. CT and MR images were registered on which parotid/submandibular glands were contoured. Image features were extracted for ipsilateral/contralateral parotid and submandibular glands relative to the location of the primary tumor. Dose-volume-histogram (DVH) parameters were also acquired. Features were pre-selected based on Spearman correlation before modelling by examining the correlation with xerostomia (p < 0.05). A shrinkage regression analysis of the pre-selected features was performed using LASSO. The internal validity of the variable selection was estimated by repeating the entire variable selection procedure using a leave-one-out-cross-validation. The most frequently selected variables were considered in the final model. A generalized linear regression with repeated ten-fold cross-validation was developed to predict radiation-induced xerostomia at 3-months after radiotherapy. This model was tested in an independent dataset (n = 50) of patients who were treated at the same institution in 2017-2018. We compared the prediction performances under eight conditions (DVH-only, CT-only, MR-only, CT + MR, DVH + CT, DVH + CT + MR, Clinical+CT + MR, and Clinical+DVH + CT + MR) using the area under the receiver operating characteristic curve (ROC-AUC). RESULTS: Among extracted features, 7 CT, 5 MR, and 2 DVH features were selected. The internal cohort (n = 216) ROC-AUC values for DVH, CT, MR, and Clinical+DVH + CT + MR features were 0.73 ± 0.01, 0.69 ± 0.01, 0.70 ± 0.01, and 0.79 ± 0.01, respectively. The validation cohort (n = 50) ROC-AUC values for DVH, CT, MR, and Clinical+DVH + CT + MR features were 0.63, 0.57, 0.66, and 0.68, respectively. The DVH-ROC was not significantly different than the CT-ROC (p = 0.8) or MR-ROC (p = 0.4). However, the CT + MR-ROC was significantly different than the CT-ROC (p = 0.03), but not the Clinical+DVH + CT + MR model (p = 0.5). CONCLUSION: Our results suggest that baseline CT and MR image features may reflect baseline salivary gland function and potential risk for radiation injury. The integration of baseline image features into prediction models has the potential to improve xerostomia risk stratification with the ultimate goal of truly personalized HNC radiotherapy.

Comparison of treatment planning approaches for spatially fractionated irradiation of deep tumors.

PURPOSE: The purpose of this work was to compare the dosimetry and delivery times of 3D-conformal (3DCRT)-, volumetric modulated arc therapy (VMAT)-, and tomotherapy-based approaches for spatially fractionated radiation therapy for deep tumor targets. METHODS: Two virtual GRID phantoms were created consisting of 7 "target" cylinders (1-cm diameter) aligned longitudinally along the tumor in a honey-comb pattern, mimicking a conventional GRID block, with 2-cm center-to-center spacing (GRID2 cm ) and 3-cm center-to-center spacing (GRID3 cm ), all contained within a larger cylinder (8 and 10 cm in diameter for the GRID2 cm and GRID3 cm , respectively). In a single patient, a GRID3 cm structure was created within the gross tumor volume (GTV). Tomotherapy, VMAT (6 MV + 6 MV-flattening-filter-free) and multi-leaf collimator segment 3DCRT (6 MV) plans were created using commercially available software. Two tomotherapy plans were created with field widths (TOMO2.5 cm ) 2.5 cm and (TOMO5 cm ) 5 cm. Prescriptions for all plans were set to deliver a mean dose of 15 Gy to the GRID targets in one fraction. The mean dose to the GRID target and the heterogeneity of the dose distribution (peak-to-valley and peak-to-edge dose ratios) inside the GRID target were obtained. The volume of normal tissue receiving 7.5 Gy was determined. RESULTS: The peak-to-valley ratios for GRID2 cm /GRID3 cm /Patient were 2.1/2.3/2.8, 1.7/1.5/2.8, 1.7/1.9/2.4, and 1.8/2.0/2.8 for the 3DCRT, VMAT, TOMO5 cm , and TOMO2.5 cm plans, respectively. The peak-to-edge ratios for GRID2 cm /GRID3 cm /Patient were 2.8/3.2/5.4, 2.1/1.8/5.4, 2.0/2.2/3.9, 2.1/2.7/5.2 and for the 3DCRT, VMAT, TOMO5 cm , and TOMO2.5 cm plans, respectively. The volume of normal tissue receiving 7.5 Gy was lowest in the TOMO2.5 cm plan (GRID2 cm /GRID3 cm /Patient = 54 cm3 /19 cm3 /10 cm3 ). The VMAT plans had the lowest delivery times (GRID2 cm /GRID3 cm /Patient = 17 min/8 min/9 min). CONCLUSION: Our results present, for the first time, preliminary evidence comparing IMRT-GRID approaches which result in high-dose "islands" within a target, mimicking what is achieved with a conventional GRID block but without high-dose "tail" regions outside of the target. These approaches differ modestly in their ability to achieve high peak-to-edge ratios and also differ in delivery times.

Dose/Volume histogram patterns in Salivary Gland subvolumes influence xerostomia injury and recovery.

Xerostomia is a common consequence of radiotherapy in head and neck cancer. The objective was to compare the regional radiation dose distribution in patients that developed xerostomia within 6 months of radiotherapy and those recovered from xerostomia within 18 months post-radiotherapy. We developed a feature generation pipeline to extract dose volume histogram features from geometrically defined ipsilateral/contralateral parotid glands, submandibular glands, and oral cavity surrogates for each patient. Permutation tests with multiple comparisons were performed to assess the dose difference between injury vs. non-injury and recovery vs. non-recovery. Ridge logistic regression models were applied to predict injury and recovery using clinical features along with dose features (D10-D90) of the subvolumes extracted from oral cavity and salivary gland contours + 3 mm peripheral shell. Model performances were assessed by the area under the receiver operating characteristic curve (AUC) using nested cross-validation. We found that different regional dose/volume metrics patterns exist for injury vs. recovery. Compared to injury, recovery has increased importance to the subvolumes receiving lower dose. Within the subvolumes, injury tends to have increased importance towards D10 from D90. This suggests that different threshold for xerostomia injury and recovery. Injury is induced by the subvolumes receiving higher dose, and the ability to recover can be preserved by further reducing the dose to subvolumes receiving lower dose.

Distinguishing True Progression From Radionecrosis After Stereotactic Radiation Therapy for Brain Metastases With Machine Learning and Radiomics.

PURPOSE: Treatment effect or radiation necrosis after stereotactic radiosurgery (SRS) for brain metastases is a common phenomenon often indistinguishable from true progression. Radiomics is an emerging field that promises to improve on conventional imaging. In this study, we sought to apply a radiomics-based prediction model to the problem of diagnosing treatment effect after SRS. METHODS AND MATERIALS: We included patients in the Johns Hopkins Health System who were treated with SRS for brain metastases who subsequently underwent resection for symptomatic growth. We also included cases of likely treatment effect in which lesions grew but subsequently regressed spontaneously. Lesions were segmented semiautomatically on preoperative T1 postcontrast and T2 fluid-attenuated inversion recovery magnetic resonance imaging, and radiomic features were extracted with software developed in-house. Top-performing features on univariate logistic regression were entered into a hybrid feature selection/classification model, IsoSVM, with parameter optimization and further feature selection performed using leave-one-out cross-validation. Final model performance was assessed by 10-fold cross-validation with 100 repeats. All cases were independently reviewed by a board-certified neuroradiologist for comparison. RESULTS: We identified 82 treated lesions across 66 patients, with 77 lesions having pathologic confirmation. There were 51 radiomic features extracted per segmented lesion on each magnetic resonance imaging sequence. An optimized IsoSVM classifier based on top-ranked radiomic features had sensitivity and specificity of 65.38% and 86.67%, respectively, with an area under the curve of 0.81 on leave-one-out cross-validation. Only 73% of cases were classifiable by the neuroradiologist, with a sensitivity of 97% and specificity of 19%. CONCLUSIONS: Radiomics holds promise for differentiating between treatment effect and true progression in brain metastases treated with SRS. A predictive model built on radiomic features from an institutional cohort performed well on cross-validation testing. These results warrant further validation in independent datasets. Such work could prove invaluable for guiding management of individual patients and assessing outcomes of novel interventions.

Development of a pulmonary imaging biomarker pipeline for phenotyping of chronic lung disease.

We designed and generated pulmonary imaging biomarker pipelines to facilitate high-throughput research and point-of-care use in patients with chronic lung disease. Image processing modules and algorithm pipelines were embedded within a graphical user interface (based on the .NET framework) for pulmonary magnetic resonance imaging (MRI) and x-ray computed-tomography (CT) datasets. The software pipelines were generated using C++ and included: (1) inhaled He3/Xe129 MRI ventilation and apparent diffusion coefficients, (2) CT-MRI coregistration for lobar and segmental ventilation and perfusion measurements, (3) ultrashort echo-time H1 MRI proton density measurements, (4) free-breathing Fourier-decomposition H1 MRI ventilation/perfusion and free-breathing H1 MRI specific ventilation, (5) multivolume CT and MRI parametric response maps, and (6) MRI and CT texture analysis and radiomics. The image analysis framework was implemented on a desktop workstation/tablet to generate biomarkers of regional lung structure and function related to ventilation, perfusion, lung tissue texture, and integrity as well as multiparametric measures of gas trapping and airspace enlargement. All biomarkers were generated within 10 min with measurement reproducibility consistent with clinical and research requirements. The resultant pulmonary imaging biomarker pipeline provides real-time and automated lung imaging measurements for point-of-care and high-throughput research.

Pulmonary MRI morphometry modeling of airspace enlargement in chronic obstructive pulmonary disease and alpha-1 antitrypsin deficiency.

PURPOSE: We generated lung morphometry measurements using single-breath diffusion-weighted MRI and three different acinar duct models in healthy participants and patients with emphysema stemming from chronic obstructive lung disease (COPD) and alpha-1 antitrypsin deficiency (AATD). METHODS: Single-breath-inhaled 3 He MRI with five diffusion sensitizations (b-value = 0, 1.6, 3.2, 4.8, and 6.4 s/cm2 ) was used, and signal intensities were fit using a cylindrical and single-compartment acinar-duct model to estimate MRI-derived mean linear intercept (Lm ) and surface-to-volume ratio (S/V). A stretched exponential model was also developed to estimate the mean airway length and Lm . RESULTS: We evaluated 42 participants, including 15 elderly never-smokers (69 ± 5 years), 12 ex-smokers without COPD (67 ± 11 years), 9 COPD ex-smokers (80 ± 6 years), and 6 AATD patients (59 ± 6 years). In the never- and ex-smokers, the diffusing capacity of the lung for carbon monoxide (DLCO ) and computed tomography relative area of less than -950 Hounsfield units (RA950 ) were normal, but these were abnormal in the COPD and AATD patients, which is reflective of emphysema. Although cylindrical and stretched-exponential-model estimates of Lm and S/V were not significantly different, the single-compartment-model estimates were significantly different (P < 0.05) for the never- and ex-smoker subgroups. All models estimated significantly worse Lm and S/V in the AATD and COPD subgroups compared with the never- and ex-smokers without emphysema. CONCLUSIONS: Differences in airspace enlargement may be estimated using Lm and S/V, generated using MRI and a stretched-exponential or cylindrical model of the acinar ducts. Magn Reson Med 79:439-448, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Pulmonary 3He Magnetic Resonance Imaging Biomarkers of Regional Airspace Enlargement in Alpha-1 Antitrypsin Deficiency.

RATIONALE AND OBJECTIVES: Thoracic x-ray computed tomography (CT) and hyperpolarized 3He magnetic resonance imaging (MRI) provide quantitative measurements of airspace enlargement in patients with emphysema. For patients with panlobular emphysema due to alpha-1 antitrypsin deficiency (AATD), sensitive biomarkers of disease progression and response to therapy have been difficult to develop and exploit, especially those biomarkers that correlate with outcomes like quality of life. Here, our objective was to generate and compare CT and diffusion-weighted inhaled-gas MRI measurements of emphysema including apparent diffusion coefficient (ADC) and MRI-derived mean linear intercept (Lm) in patients with AATD, chronic obstructive pulmonary disease (COPD) ex-smokers, and elderly never-smokers. MATERIALS AND METHODS: We enrolled patients with AATD (n = 8; 57 ± 7 years), ex-smokers with COPD (n = 8; 77 ± 6 years), and a control group of never-smokers (n = 5; 64 ± 2 years) who underwent thoracic CT, MRI, spirometry, plethysmography, the St. George's Respiratory Questionnaire, and the 6-minute walk test during a single 2-hour visit. MRI-derived ADC, Lm, surface-to-volume ratio, and ventilation defect percent were generated for the apical, basal, and whole lung as was CT lung area ≤-950 Hounsfield units (RA950), low attenuating clusters, and airway count. RESULTS: In patients with AATD, there was a significantly different MRI-derived ADC (P = .03), Lm (P < .0001), and surface-to-volume ratio (P < .0001), but not diffusing capacity of carbon monoxide, residual volume or total lung capacity, or CT RA950 (P > .05) compared to COPD ex-smokers with a significantly different St. George's Respiratory Questionnaire. CONCLUSIONS: In this proof-of-concept demonstration, we evaluated CT and MRI lung emphysema measurements and observed significantly worse MRI biomarkers of emphysema in patients with AATD compared to patients with COPD, although CT RA950 and diffusing capacity of carbon monoxide were not significantly different, underscoring the sensitivity of MRI measurements of AATD emphysema.

Free-breathing Functional Pulmonary MRI: Response to Bronchodilator and Bronchoprovocation in Severe Asthma.

RATIONALE AND OBJECTIVES: Ventilation heterogeneity is a hallmark feature of asthma. Our objective was to evaluate ventilation heterogeneity in patients with severe asthma, both pre- and post-salbutamol, as well as post-methacholine (MCh) challenge using the lung clearance index, free-breathing pulmonary 1H magnetic resonance imaging (FDMRI), and inhaled-gas MRI ventilation defect percent (VDP). MATERIALS AND METHODS: Sixteen severe asthmatics (49 ± 10 years) provided written informed consent to an ethics board-approved protocol. Spirometry, plethysmography, and multiple breath nitrogen washout to measure the lung clearance index were performed during a single visit within 15 minutes of MRI. Inhaled-gas MRI and FDMRI were performed pre- and post-bronchodilator to generate VDP. For asthmatics with forced expiratory volume in 1 second (FEV1) >70%predicted, MRI was also performed before and after MCh challenge. Wilcoxon signed-rank tests, Spearman correlations, and a repeated-measures analysis of variance were performed. RESULTS: Hyperpolarized 3He (P = .02) and FDMRI (P = .02) VDP significantly improved post-salbutamol and for four asthmatics who could perform MCh (n = 4). 3He and FDMRI VDP significantly increased at the provocative concentration of MCh, resulting in a 20% decrease in FEV1 (PC20) and decreased post-bronchodilator (P = .02), with a significant difference between methods (P = .01). FDMRI VDP was moderately correlated with 3He VDP (ρ = .61, P = .01), but underestimated VDP relative to 3He VDP (-6 ± 9%). Whereas 3He MRI VDP was significantly correlated with the lung clearance index, FDMRI was not (ρ = .49, P = .06). CONCLUSIONS: FDMRI VDP generated in free-breathing asthmatic patients was correlated with static inspiratory breath-hold 3He MRI VDP but underestimated VDP relative to 3He MRI VDP. Although less sensitive to salbutamol and MCh, FDMRI VDP may be considered for asthma patient evaluations at centers without inhaled-gas MRI.

Thoracic CT-MRI coregistration for regional pulmonary structure-function measurements of obstructive lung disease.

PURPOSE: Recent pulmonary imaging research has revealed that in patients with chronic obstructive pulmonary disease (COPD) and asthma, structural and functional abnormalities are spatially heterogeneous. This novel information may help optimize treatment in individual patients, monitor interventional efficacy, and develop new treatments. Moreover, by automating the measurement of regional biomarkers for the 19 different anatomical lung segments, there is an opportunity to embed imaging biomarkers into clinically acceptable clinical workflows and improve lung disease clinical care. Therefore, to exploit the regional structure-function information provided by thoracic imaging, and as a first step toward this goal, our objective was to develop a fully automated registration pipeline for thoracic x-ray computed tomography (CT) and inhaled gas functional magnetic resonance imaging (MRI) whole lung and segmental structure-function biomarkers. METHODS: Thirty-five patients including 15 severe, poorly controlled asthmatics and 20 COPD patients [classified according to the global initiative for chronic obstructive lung disease (GOLD) criteria)] provided written informed consent to a study protocol approved by Health Canada and underwent pulmonary function tests, MRI, and CT during a single 2-hour visit. Using this diverse patient dataset, we developed and evaluated a joint deformable registration approach to simultaneously coregister CT with both 1 H and 3 He MRI by enforcing the similarity of the deformation fields from the two individual registrations. We derived a simpler model that was equivalent to the original challenging optimization problem through variational analysis and the simpler model gave rise to an efficient numerical solver that was parallelized on a graphics processing unit. The coregistered CT-3 He MRI and whole lung/segmental lung masks were used to generate whole lung and segmental 3 He MRI ventilation defect percent (VDP). To estimate fiducial localization reproducibility, a single observer manually identified 109 pairs of CT and 3 He MRI fiducials for 35 patient images on five separate occasions and determined the fiducial localization error (FLE). CT-3 He MRI registration accuracy was evaluated using the target registration error (TRE). Whole lung VDP generated using the algorithm was compared with VDP generated using a previously validated semiautomated approach and computational efficiency was evaluated using run time. RESULTS: In 35 patients including 15 with severe asthma and 20 with COPD, mean forced expiratory volume in 1 s (FEV1 ) was 63±24%pred and FEV1 /forced vital capacity (FVC) was 54 ± 17%. FLE was 0.16 mm and 0.34 mm for 3 He MRI and CT, respectively. TRE was 4.5 ± 2.0 mm, 4.0 ± 1.7 mm, 4.8 ± 2.3 mm for asthma, COPD GOLD II, and GOLD III groups, respectively, with a mean of 4.4 ± 2.0 mm for the entire dataset. TRE was significantly improved for joint CT-1 H/3 He MRI registration compared with CT-1 H MRI rigid registration (P < 0.0001). Whole lung VDP generated using the pipeline was not significantly different (P = 0.37) compared to a semiautomated method with which it was strongly correlated (r = 0.93, P < 0.0001). The fully automated pipeline required 11 ± 0.4 min to generate whole lung and segmental VDP. CONCLUSIONS: For a diverse group of patients with COPD and asthma, whole lung and segmental VDP was measured using an automated lung image analysis pipeline which provides a way to incorporate lung functional biomarkers into clinical research and patient care.