An efficient rectangular optimization method for sparse orthogonal collimator based small animal irradiation

Objective.Intensity-modulated radiotherapy (IMRT) is widely used in clinical radiotherapy, treating varying malignancies with conformal doses. As the test field for clinical translation, preclinical small animal experiments need to mimic the human radiotherapy condition, including IMRT. However, small animal IMRT is a systematic challenge due to the lack of corresponding hardware and software for miniaturized targets.Approach.The sparse orthogonal collimators (SOC) based on the direct rectangular aperture optimization (RAO) substantially simplified the hardware for miniaturization. This study investigates and evaluates a significantly improved RAO algorithm for complex mouse irradiation using SOC. Because the Kronecker product representation of the rectangular aperture is the main limitation of the computational performance, we reformulated matrix multiplication in the data fidelity term using multiplication with small matrices instead of the Kronecker product of the dose loading matrices. Solving the optimization problem was further accelerated using the Fast Iterative Shrinkage-Thresholding Algorithm (FISTA).Main results.Four mouse cases, including a liver, a brain tumor, a concave U-target, and a complex total marrow irradiation (TMI) case, were included in this study with manually delineated targets and OARs. Seven coplanar-field SOC IMRT (sIMRT) plans were compared with idealistic fluence map based IMRT (iIMRT) plans. For the first three cases with simpler and smaller targets, the differences between sIMRT plans and iIMRT plans in the planning target volumes (PTV) statistics are within 1%. For the TMI case, the sIMRT plans are superior in reducing hot spots (also termedDmax) of PTV, kidneys, lungs, heart, and bowel by 20.5%, 31.5%, 24.67%, 20.13%, and 17.78%, respectively. On average, in four cases in this study, the sIMRT plan conformity is comparable to that of the iIMRT's with lightly increased R50 and Integral Dose by 2.23% and 2.78%.Significance.The significantly improved sIMRT optimization method allows fast plan creation in under 1 min for smaller targets and makes complex TMI planning feasible while achieving comparable dosimetry to idealistic IMRT with fluence map optimization.

Feasibility of a Novel Sparse Orthogonal Collimator-Based Preclinical Total Marrow Irradiation for Enhanced Dosimetric Conformality

Total marrow irradiation (TMI) has significantly improved radiation conditioning for hematopoietic cell transplantation in hematologic diseases by reducing conditioning-induced toxicities and improving survival outcomes in relapsed/refractory patients. Recently, preclinical three-dimensional image-guided TMI has been developed to enhance mechanistic understanding of the role of TMI and to support the development of experimental therapeutics. However, a dosimetric comparison between preclinical and clinical TMI reveals that the preclinical TMI treatment lacks the ability to reduce the dose to some of the vital organs that are very close to the skeletal system and thus limits the ability to evaluate radiobiological relevance. To overcome this limit, we introduce a novel Sparse Orthogonal Collimator (SOC)-based TMI and evaluate its ability to enhance dosimetric conformality. The SOC-TMI-based dose modulation technique significantly improves TMI treatment planning by reducing radiation exposures to critical organs that are close to the skeletal system that leads to reducing the gap between clinical and preclinical TMI.

Urethral Interfractional Geometric and Dosimetric Variations of Prostate Cancer Patients: A Study Using an Onboard MRI

Purpose

For a cohort of prostate cancer patients treated on an MR-guided radiotherapy (MRgRT) system, we retrospectively analyzed urethral interfractional geometric and dosimetric variations based on onboard MRIs acquired at different timepoints and evaluated onboard prostatic urethra visualization for urethra-focused online adaptive RT.

Methods

Twenty-six prostate cancer patients were prospectively scanned on a 0.35-T MRgRT system using an optimized T2-weighted HASTE sequence at simulation and final fraction. Two radiation oncologists (RO1 and RO2) contoured the urethras on all HASTE images. The simulation and final fraction HASTE images were rigidly registered, and urethral interobserver and interfractional geometric variation was evaluated using the 95th percentile Hausdorff distance (HD95), mean distance to agreement (MDA), center-of-mass shift (COMS), and DICE coefficient. For dosimetric analysis, simulation and final fraction HASTE images were registered to the 3D bSSFP planning MRI and 3D bSSFP final setup MRI, respectively. Both ROs’ urethra contours were transferred from HASTE images for initial treatment plan optimization and final fraction dose estimation separately. Stereotactic body radiotherapy (SBRT) plans, 40 Gy in 5 fractions, were optimized to meet clinical constraints, including urethral V42Gy ≤0.03 cc, on the planning MRI. The initial plan was then forward calculated on the final setup MRI to estimate urethral dose on the final fraction and evaluate urethral dosimetric impact due to anatomy change.

Results

The average interobserver HD95, MDA, COMS, and DICE were 2.85 ± 1.34 mm, 1.02 ± 0.36 mm, 3.16 ± 1.61 mm, and 0.58 ± 0.15, respectively. The average interfractional HD95, MDA, COMS, and DICE were 3.26 ± 1.54 mm, 1.29 ± 0.54 mm, 3.34 ± 2.01 mm, and 0.49 ± 0.18, respectively. All patient simulation MRgRT plans met all clinical constraints. For RO1 and RO2, 23/26 (88%) and 21/26 (81%) patients’ final fraction estimated urethral dose did not meet the planned constraint. The average urethral V42Gy change was 0.48 ± 0.58 cc.

Conclusion

Urethral interfractional motion and anatomic change can result in daily treatment violating urethral constraints. Onboard MRI with good visualization of the prostatic urethra can be a valuable tool to help better protect the urethra through patient setup or online adaptive RT.

Voxelwise Prediction of Recurrent High-Grade Glioma via Proximity Estimation-Coupled Multidimensional Support Vector Machine

PURPOSE

To provide early and localized glioblastoma (GBM) recurrence prediction, we introduce a novel postsurgery multiparametric magnetic resonance-based support vector machine (SVM) method coupling with stem cell niche (SCN) proximity estimation.

METHODS AND MATERIALS

This study used postsurgery magnetic resonance imaging (MRI) scans from 50 patients with recurrent GBM, obtained approximately 2 months before clinically diagnosed recurrence. The main prediction pipeline consisted of a proximity-based estimator to identify regions with high risk of recurrence (HRRs) and an SVM classifier to provide voxelwise prediction in HRRs. The HRRs were estimated using the weighted sum of inverse distances to 2 possible origins of recurrence-the SCN and the tumor cavity. Subsequently, multiparametric voxels (from T1, T1 contrast-enhanced, fluid-attenuated inversion recovery, T2, and apparent diffusion coefficient) within the HRR were grouped into recurrent (warped from the clinical diagnosis) and nonrecurrent subregions and fed into the proximity estimation-coupled SVM classifier (SVM_PE). The cohort was randomly divided into 40% and 60% for training and testing, respectively. The trained SVM_PE was then extrapolated to an earlier time point for earlier recurrence prediction. As an exploratory analysis, the SVM_PE predictive cluster sizes and the image intensities from the 5 magnetic resonance sequences were compared across time to assess the progressive subclinical traces.

RESULTS

On 2-month prerecurrence MRI scans from 30 test cohort patients, the SVM_PE classifier achieved a recall of 0.80, a precision of 0.69, an F1-score of 0.73, and a mean boundary distance of 7.49 mm. Exploratory analysis at early time points showed spatially consistent but significantly smaller subclinical clusters and significantly increased T1 contrast-enhanced and apparent diffusion coefficient values over time.

CONCLUSIONS

We demonstrated a novel voxelwise early prediction method, SVM_PE, for GBM recurrence based on clinical follow-up MR scans. The SVM_PE is promising in localizing subclinical traces of recurrence 2 months ahead of clinical diagnosis and may be used to guide more effective personalized early salvage therapy.

Reformulated McNamara RBE-weighted beam orientation optimization for intensity modulated proton therapy

PURPOSE

Empirical relative biological effectiveness (RBE) models have been used to estimate the biological dose in proton therapy but do not adequately capture the factors influencing RBE values for treatment planning. We reformulate the McNamara RBE model such that it can be added as a linear biological dose fidelity term within our previously developed sensitivity-regularized and heterogeneity-weighted beam orientation optimization (SHBOO) framework.

METHODS

Based on our SHBOO framework, we formulated the biological optimization problem to minimize total McNamara RBE dose to OARs. We solve this problem using two optimization algorithms: FISTA (McNam-FISTA) and Chambolle-Pock (McNam-CP). We compare their performances with a physical dose optimizer assuming RBE = 1.1 in all structures (PHYS-FISTA) and an LET-weighted dose model (LET-FISTA). Three head and neck patients were planned with the four techniques and compared on dosimetry and robustness.

RESULTS

Compared to Phys-FISTA, McNam-CP was able to match CTV [HI, Dmax, D95%, D98%] by [0.00, 0.05%, 1.4%, 0.8%]. McNam-FISTA and McNam-CP were able to significantly improve overall OAR [Dmean, Dmax] by an average of [36.1%,26.4%] and [29.6%, 20.3%], respectively. Regarding CTV robustness, worst [Dmax, V95%, D95%, D98%] improvement of [-6.6%, 6.2%, 6.0%, 4.8%] was reported for McNam-FISTA and [2.7%, 2.7%, 5.3%, -4.3%] for McNam-CP under combinations of range and setup uncertainties. For OARs, worst [Dmax, Dmean] were improved by McNam-FISTA and McNam-CP by an average of [25.0%, 19.2%] and [29.5%, 36.5%], respectively. McNam-FISTA considerably improved dosimetry and CTV robustness compared to LET-FISTA, which achieved better worst-case OAR doses.

CONCLUSION

The four optimization techniques deliver comparable biological doses for the head and neck cases. Besides modest CTV coverage and robustness improvement, OAR biological dose and robustness were substantially improved with both McNam-FISTA and McNam-CP, showing potential benefit for directly incorporating McNamara RBE in proton treatment planning.

Fixed Beamline Optimization for Intensity Modulated Carbon-Ion Therapy

A major obstacle for the adoption of heavy ion therapy is the cost and technical difficulties to construct and maintain a rotational gantry. Many heavy ion treatment facilities instead choose to construct fixed beamlines as a compromise, which we propose to mitigate with optimized treatment couch angle. We formulate the integrated beam orientation and scanning spot optimization problem as a quadratic cost function with a group sparsity regularization term. The optimization problem is efficiently solved using fast iterative shrinkage-thresholding algorithm (FISTA). To test the method, we created the fixed beamline plans with couch rotation (FBCR) and without couch rotation (FB) for intensity modulated carbon-ion therapy (IMCT) and compared with the ideal scenario where both the couch and gantry have 360 degrees of freedom (GCR). FB, FBCR, and GCR IMCT plans were compared for ten pancreas cases. The FBCR plans show comparable PTV coverage and OAR doses for each pancreas case. In conclusion, the dosimetric limitation of fixed beams in heavy ion radiotherapy may be largely mitigated with integrated beam orientation optimization of the couch rotation.

Magnetic resonance imaging-guided versus computed tomography-guided stereotactic body radiotherapy for prostate cancer (MIRAGE): Interim analysis of a phase III randomized trial

255

Background: Magnetic resonance imaging (MRI) guidance offers several theoretical advantages over computed tomography (CT) guidance in the context of stereotactic body radiotherapy (SBRT) for intact prostate cancer. Here we report the results of an interim analysis of the phase III MIRAGE trial, which directly compared MRI- and CT-guidance with a pragmatic primary endpoint of acute grade ≥2 genitourinary (GU) toxicity. Methods: MIRAGE is a single center, randomized phase 3 trial. Men undergoing SBRT for localized prostate cancer were randomly assigned to either CT-guidance or MRI-guidance. Planning margins of 4 mm (CT-arm) and 2 mm (MRI-arm) were placed around the prostate and proximal seminal vesicles, and this volume received 40 Gy in five fractions. Elective nodal radiotherapy and rectal spacers were allowed per physician discretion. The primary outcome was the incidence of acute (i.e., within 90 days of SBRT) grade ≥2 GU physician-reported toxicity (by CTCAE version 4.03). Secondary outcomes of interest included the incidence of acute grade ≥2 GU physician-reported toxicity, changes in IPSS scores at 1 and 3 months, and changes in EPIC-26 bowel domain summary scores at 1 and 3 months. A pre-specified efficacy analysis was planned once the 100th patient was eligible for evaluation of the primary endpoint. Results: On 9/1/2021, 100 patients became eligible for evaluation for the interim analysis (51 CT arm, 49 MRI arm). Acute grade ≥2 GU toxicity was significantly reduced in men receiving MRI-guided SBRT (incidence of 24 (47.1%) vs. 11 (22.4%), p = 0.01). Acute grade ≥2 GI toxicity was also significantly reduced in men receiving MRI-guided SBRT (incidence of 7 (13.7%) vs. 0 (0%), p = 0.01.). The increase in IPSS scores from baseline was significantly higher in men receiving CT-guided SBRT at 1 month post-SBRT (median change of 10 vs. 6, p = 0.03), but not at 3 months (median change of 3 vs. 2, p = 0.3). The decrement in EPIC-26 bowel domain scores was significantly greater at 1 month in men receiving CT-guided SBRT (median change of -8.3 vs. 0, p = 0.03), but not at 3 months (median change of -2.3 vs. 0, p = 0.4). Given the large primary endpoint signal seen, our protocol was amended to reduce the projected sample size to 154 while still maintaining 89% power to detect a difference. Conclusions: This interim analysis demonstrates a statistically significant reduction in acute grade ≥2 GU toxicity with MRI-guidance versus CT-guidance in the context of prostate SBRT. Patient-reported urinary and bowel function metrics are also better preserved at the 1 month time point with MRI-guidance, though this difference dissipates (potentially due to side-effect management) at the 3 month time point. Accrual has been completed as of October 2021 and a final analysis for the primary endpoint is anticipated in early 2022. Clinical trial information: NCT04384770.

A Prospective Phase II Study of Automated Non-Coplanar VMAT for Recurrent Head and Neck Cancer: Initial Report of Feasibility, Safety, and Patient-Reported Outcomes

Simple Summary

The delivery of higher radiation doses has been shown to increase local control, and ultimately survival, for head and neck cancer patients, but highly conformal dose distributions are necessary to minimize normal tissue toxicity. Varian’s HyperArc non-coplanar automated treatment planning and delivery technique has been shown to improve dose conformity for intracranial treatment, but its safety and efficacy for head and neck cancer treatment has yet to be verified. This study evaluates the initial results of a prospective clinical trial using HyperArc for recurrent head and neck cancer patients. We demonstrated that HyperArc can enable significant tumor dose escalation compared to conventional volumetric modulated arc therapy (VMAT) planning while minimizing the dose to organs at risk. Treatment delivery was feasible and safe, with minimal treatment-related toxicities and positive patient-reported quality of life measures.

Abstract

This study reports the initial results for the first 15 patients on a prospective phase II clinical trial exploring the safety, feasibility, and efficacy of the HyperArc technique for recurrent head and neck cancer treatment. Eligible patients were simulated and planned with both conventional VMAT and HyperArc techniques and the plan with superior dosimetry was selected for treatment. Dosimetry, delivery feasibility and safety, treatment-related toxicity, and patient-reported quality of life (QOL) were all evaluated. HyperArc was chosen over conventional VMAT for all 15 patients and enabled statistically significant increases in dose conformity (R50% reduced by 1.2 ± 2.1, p < 0.05) and mean PTV and GTV doses (by 15.7 ± 4.9 Gy, p < 0.01 and 17.1 ± 6.0 Gy, p < 0.01, respectively). The average HyperArc delivery was 2.8 min longer than conventional VMAT (p < 0.01), and the mean intrafraction motion was ≤ 0.5 ± 0.4 mm and ≤0.3 ± 0.1°. With a median follow-up of 12 months, treatment-related toxicity was minimal (only one grade 3 acute toxicity above baseline) and patient-reported QOL metrics were favorable. HyperArc enabled superior dosimetry and significant target dose escalation compared to conventional VMAT planning, and treatment delivery was feasible, safe, and well-tolerated by patients.

Bladder surface dose modeling in prostate cancer radiotherapy: An analysis of motion-induced variations and the cumulative dose across the treatment

PURPOSE

To introduce a novel surface-based dose mapping method to improve quantitative bladder dosimetric assessment in prostate cancer (PC) radiotherapy.

METHODS

Based on the planning and daily pre and postfraction MRIs of 12 PC patients, bladder surface models (SMs) were generated on manually delineated contours and regionally aligned via surface-based registration. Subsequently, bladder surface dose models (SDMs) were created using face-wise dose sampling. To determine the bladder intrafractional and interfractional motion and dose variation, we performed a pose analysis between pre and postfraction bladder SMs, as well as surface mapping for fractional SMs. Discrepancies between the received dose, accumulated from daily SDMs, and the planned dose were then assessed on the corresponding SDMs. Complementary to the surface dose mapping, dose surface histogram (DSH)-based comparisons were also performed.

RESULTS

The intrafraction pose analysis revealed a significant (p < 0.05) bladder expansion, as well as an anterior/superior drift during the treatment. The intrafraction motion substantially altered dose to mid-bladder body, but not the bladder surface areas distal to or contiguous with the target. A similar pattern of dose variations was also detected by interfraction comparisons. With surface registration to the common SM, the cumulative bladder dose significantly differs from the planned dose. The discrepancy is evident in the mid-posterior range that corresponds to a mid- to high-dose region. The received DSH significantly differs from the planned DSH after permutation correction (p = 0.0122), while the overall surface-based comparison after multiple comparison correction is nonsignificant (p = 0.0800).

CONCLUSIONS

We developed a novel surface-based intra and interdose mapping framework applied to a unique daily MR dataset for image-guided radiotherapy. The framework identified significant intrafraction bladder positional changes, localized the intra and interfraction variations, and quantified planned versus received dose differences on the bladder surface. The result indicates the importance of adopting the motion-integrated bladder SDM for bladder dose management.

Quantitative Characterization of Tumor Proximity to Stem Cell Niches: Implications on Recurrence and Survival in GBM Patients

PURPOSE

Emerging evidence has linked glioblastoma multiforme (GBM) recurrence and survival to stem cell niches (SCNs). However, the traditional tumor-ventricle distance is insufficiently powered for an accurate prediction. We aimed to use a novel inverse distance map for improved prediction.

METHODS AND MATERIALS

Two T1-magnetic resonance imaging data sets were included for a total of 237 preoperative scans for prognostic stratification and 55 follow-up scans for recurrent pattern identification. SCN, including the subventricular zone (SVZ) and subgranular zone (SGZ), were manually defined on a standard template. A proximity map was generated using the summed inverse distances to all SCN voxels. The mean and maximum proximity scores (PS_m-SCN and PS_max-SCN) were calculated for each primary/recurrent tumor, deformably transformed into the template. The prognostic capacity of proximity score (PS)-derived metrics was assessed using Cox regression and log-rank tests. To evaluate the impact of SCNs on recurrence patterns, we performed group comparisons of PS-derived metrics between the primary and recurrent tumors. For comparison, the same analyses were conducted on PS derived from SVZ alone and traditional edge/center-to-ventricle metrics.

RESULTS

Among all SCN-derived features, PS_m-SCN was the strongest survival predictor (P < .0001). PS_max-SCN was the best in risk stratification, using either evenly sorted (P = .0001) or k-means clustering methods (P = .0045). PS metrics based on SVZ only also correlated with overall survival and risk stratification, but to a lesser degree of significance. In contrast, edge/center-to-ventricle metrics showed weak to no prediction capacities in either task. Moreover, PS_m-SCN,PS_m-SVZ, and center-to-ventricle metrics revealed a significantly closer SCN distribution of recurrence than primary tumors.

CONCLUSIONS

We introduced a novel inverse distance-based metric to comprehensively capture the anatomic relationship between GBM tumors and SCN zones. The derived metrics outperformed traditional edge or center distance-based measurements in overall survival prediction, risk stratification, and recurrent pattern differentiation. Our results reveal the potential role of SGZ in recurrence aside from SVZ.

Using neural networks to extend cropped medical images for deformable registration among images with differing scan extents

PURPOSE

Missing or discrepant imaging volume is a common challenge in deformable image registration (DIR). To minimize the adverse impact, we train a neural network to synthesize cropped portions of head and neck CT's and then test its use in DIR.

METHODS

Using a training dataset of 409 head and neck CT's, we trained a generative adversarial network to take in a cropped 3D image and output an image with synthesized anatomy in the cropped region. The network used a 3D U-Net generator along with Visual Geometry Group (VGG) deep feature losses. To test our technique, for each of the 53 test volumes, we used Elastix to deformably register combinations of a randomly cropped, full, and synthetically full volume to a single cropped, full, and synthetically full target volume. We additionally tested our method's robustness to crop extent by progressively increasing the amount of cropping, synthesizing the missing anatomy using our network, and then performing the same registration combinations. Registration performance was measured using 95% Hausdorff distance across 16 contours.

RESULTS

We successfully trained a network to synthesize missing anatomy in superiorly and inferiorly cropped images. The network can estimate large regions in an incomplete image, far from the cropping boundary. Registration using our estimated full images was not significantly different from registration using the original full images. The average contour matching error for full image registration was 9.9 mm, whereas our method was 11.6, 12.1, and 13.6 mm for synthesized-to-full, full-to-synthesized, and synthesized-to-synthesized registrations, respectively. In comparison, registration using the cropped images had errors of 31.7 mm and higher. Plotting the registered image contour error as a function of initial preregistered error shows that our method is robust to registration difficulty. Synthesized-to-full registration was statistically independent of cropping extent up to 18.7 cm superiorly cropped. Synthesized-to-synthesized registration was nearly independent, with a -0.04 mm of change in average contour error for every additional millimeter of cropping.

CONCLUSIONS

Different or inadequate in scan extent is a major cause of DIR inaccuracies. We address this challenge by training a neural network to complete cropped 3D images. We show that with image completion, the source of DIR inaccuracy is eliminated, and the method is robust to varying crop extent.

Clinical assessment of geometric distortion for a 0.35T MR-guided radiotherapy system

Purpose

To estimate the overall spatial distortion on clinical patient images for a 0.35 T MR‐guided radiotherapy system.

Methods

Ten patients with head‐and‐neck cancer underwent CT and MR simulations with identical immobilization. The MR images underwent the standard systematic distortion correction post‐processing. The images were rigidly registered and landmark‐based analysis was performed by an anatomical expert. Distortion was quantified using Euclidean distance between each landmark pair and tagged by tissue interface: bone‐tissue, soft tissue, or air‐tissue. For baseline comparisons, an anthropomorphic phantom was imaged and analyzed.

Results

The average spatial discrepancy between CT and MR landmarks was 1.15 ± 1.14 mm for the phantom and 1.46 ± 1.78 mm for patients. The error histogram peaked at 0–1 mm. 66% of the discrepancies were <2 mm and 51% <1 mm. In the patient data, statistically significant differences (p‐values < 0.0001) were found between the different tissue interfaces with averages of 0.88 ± 1.24 mm, 2.01 ± 2.20 mm, and 1.41 ± 1.56 mm for the air/tissue, bone/tissue, and soft tissue, respectively. The distortion generally correlated with the in‐plane radial distance from the image center along the longitudinal axis of the MR.

Conclusion

Spatial distortion remains in the MR images after systematic distortion corrections. Although the average errors were relatively small, large distortions observed at bone/tissue interfaces emphasize the need for quantitative methods for assessing and correcting patient‐specific spatial distortions.

Interfractional Geometric Variations and Dosimetric Benefits of Stereotactic MRI Guided Online Adaptive Radiotherapy (SMART) of Prostate Bed after Radical Prostatectomy: Post-Hoc Analysis of a Phase II Trial

PURPOSE

To evaluate geometric variations of patients receiving stereotactic body radiotherapy (SBRT) after radical prostatectomy and the dosimetric benefits of stereotactic MRI guided adaptive radiotherapy (SMART) to compensate for these variations.

MATERIALS/METHODS

The CTV and OAR were contoured on 55 MRI setup scans of 11 patients treated with an MR-LINAC and enrolled in a phase II trial of post-prostatectomy SBRT. All patients followed institutional bladder and rectum preparation protocols and received five fractions of 6-6.8 Gy to the prostate bed. Interfractional changes in volume were calculated and shape deformation was quantified by the Dice similar coefficient (DSC). Changes in CTV-V95%, bladder and rectum maximum dose, V32.5Gy and V27.5Gy were predicted by recalculating the initial plan on daily MRI. SMART was retrospectively simulated if the predicted dose exceeded pre-set criteria.

RESULTS

The CTV volume and shape remained stable with a median volumetric change of 3.0% (IQR -3.0% to 11.5%) and DSC of 0.83 (IQR 0.79 to 0.88). Relatively large volumetric changes in bladder (median -24.5%, IQR -34.6% to 14.5%) and rectum (median 5.4%, IQR - 9.7% to 20.7%) were observed while shape changes were moderate (median DSC of 0.79 and 0.73, respectively). The median CTV-V95% was 98.4% (IQR 94.9% to 99.6%) for the predicted doses. However, SMART would have been deemed beneficial for 78.2% of the 55 fractions based on target undercoverage (16.4%), exceeding OAR constraints (50.9%), or both (10.9%). Simulated SMART improved the dosimetry and met dosimetric criteria in all fractions. Moderate correlations were observed between the CTV-V95% and target DSC (R2 = 0.73) and bladder mean dose versus volumetric changes (R2 = 0.61).

CONCLUSIONS

Interfractional dosimetric variations resulting from anatomic deformation are commonly encountered with post-prostatectomy RT and can be mitigated with SMART.

Magnetic resonance imaging-guided stereotactic body radiotherapy for prostate cancer (mirage): a phase iii randomized trial

Background

Stereotactic body radiotherapy (SBRT) is becoming increasingly used in treating localized prostate cancer (PCa), with evidence showing similar toxicity and efficacy profiles when compared with longer courses of definitive radiation. Magnetic resonance imaging (MRI)-guided radiotherapy has multiple potential advantages over standard computed tomography (CT)-guided radiotherapy, including enhanced prostate visualization (abrogating the need for fiducials and MRI fusion), enhanced identification of the urethra, the ability to track the prostate in real-time, and the capacity to perform online adaptive planning. However, it is unknown whether these potential advantages translate into improved outcomes. This phase III randomized superiority trial is designed to prospectively evaluate whether toxicity is lower after MRI-guided versus CT-guided SBRT.

Methods

Three hundred men with localized PCa will be randomized in a 1:1 ratio to SBRT using CT or MRI guidance. Randomization will be stratified by baseline International Prostate Symptom Score (IPSS) (≤15 or > 15) and prostate gland volume (≤50 cc or > 50 cc). Five fractions of 8 Gy will be delivered to the prostate over the course of fourteen days, with or without hormonal therapy and elective nodal radiotherapy (to a dose of 5 Gy per fraction) as per the investigator’s discretion. The primary endpoint is the incidence of physician-reported acute grade ≥ 2 genitourinary (GU) toxicity (during the first 90 days after SBRT), as assessed by the CTCAE version 4.03 scale. Secondary clinical endpoints include incidence of acute grade ≥ 2 gastrointestinal (GI) toxicity, 5-year cumulative incidences of physician-reported late grade ≥ 2 GU and GI toxicity, temporal changes in patient-reported quality of life (QOL) outcomes, 5-year biochemical recurrence-free survival and the proportion of fractions of MRI-guided SBRT in which online adaptive radiotherapy is used.

Discussion

The MIRAGE trial is the first randomized trial comparing MRI-guided with standard CT-guided SBRT for localized PCa. The primary hypothesis is that MRI-guided SBRT will lead to an improvement in the cumulative incidence of acute grade ≥ 2 GU toxicity when compared to CT-guided SBRT. The pragmatic superiority design focused on an acute toxicity endpoint will allow an early comparison of the two technologies.

Trial registration

Clinicaltrials.gov identifier: NCT04384770. Date of registration: May 12, 2020.

https://clinicaltrials.gov/ct2/show/NCT04384770

Protocol version

Version 2.1, Aug 28, 2020.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-021-08281-x.

Automated Non-Coplanar VMAT for Dose Escalation in Recurrent Head and Neck Cancer Patients

Simple Summary

The ability to escalate the radiation dose to head and neck tumors has been shown to offer improved local control, and consequently, survival for recurrent head and neck cancer (rHNC) patients. This study evaluates the HyperArc automated non-coplanar planning technique (originally developed for intracranial treatment) for 20 rHNC patients, and compares this technique to conventional planning methods. HyperArc enables significant tumor dose escalation, with average increases in mean target dose of over 11.5 Gy (26%), while maintaining clinically-equivalent doses to nearby organs. Our results show that the average probability of tumor control is 23% higher for HyperArc than conventional techniques.

Abstract

This study evaluates the potential for tumor dose escalation in recurrent head and neck cancer (rHNC) patients with automated non-coplanar volumetric modulated arc therapy (VMAT) stereotactic body radiation therapy (SBRT) planning (HyperArc). Twenty rHNC patients are planned with conventional VMAT SBRT to 40 Gy while minimizing organ-at-risk (OAR) doses. They are then re-planned with the HyperArc technique to match these minimal OAR doses while escalating the target dose as high as possible. Then, we compare the dosimetry, tumor control probability (TCP), and normal tissue complication probability (NTCP) for the two plan types. Our results show that the HyperArc technique significantly increases the mean planning target volume (PTV) and gross tumor volume (GTV) doses by 10.8 ± 4.4 Gy (25%) and 11.5 ± 5.1 Gy (26%) on average, respectively. There are no clinically significant differences in OAR doses, with maximum dose differences of <2 Gy on average. The average TCP is 23% (± 21%) higher for HyperArc than conventional plans, with no significant differences in NTCP for the brainstem, cord, mandible, or larynx. HyperArc can achieve significant tumor dose escalation while maintaining minimal OAR doses in the head and neck—potentially enabling improved local control for rHNC SBRT patients without increased risk of treatment-related toxicities.

Weak Magnetic Fields Enhance the Efficacy of Radiation Therapy

Purpose

The clinical efficacy of radiation therapy is mechanistically linked to ionization-induced free radicals that cause cell and tissue injury through direct and indirect mechanisms. Free radical reaction dynamics are influenced by many factors and can be manipulated by static weak magnetic fields (WMF) that perturb singlet-triplet state interconversion. Our study exploits this phenomenon to directly increase ionizing radiation (IR) dose absorption in tumors by combining WMF with radiation therapy as a new and effective method to improve treatment.

Methods and Materials

Coils were custom made to produce both homogeneous and gradient magnetic fields. The gradient coil enabled simultaneous in vitro assessment of free radical/reactive oxygen species reactivity across multiple field strengths from 6 to 66 G. First, increases in IR-induced free radical concentrations using oxidant-sensitive fluorescent dyes in a cell-free system were measured and verified. Next, human and murine cancer cell lines were evaluated in in vitro and in vivo models after exposure to clinically relevant doses of IR in combination with WMF.

Results

Cellular responses to IR and WMF were field strength and cell line dependent. WMF was able to enhance IR effects on reactive oxygen species formation, DNA double-strand break formation, cell death, and tumor growth.

Conclusions

We demonstrate that the external presence of a magnetic field enhances radiation-induced cancer cell injury and death in vitro and in vivo. The effect extends beyond the timeframe when free radicals are induced in the presence of radiation into the window when endogenous free radicals are produced and therefore extends the applicability of this novel adjunct to cancer therapy in the context of radiation treatment.

Treating Glioblastoma Multiforme (GBM) with super hyperfractionated radiation therapy: Implication of temporal dose fractionation optimization including cancer stem cell dynamics

PURPOSE

A previously developed ordinary differential equation (ODE) that models the dynamic interaction and distinct radiosensitivity between cancer stem cells (CSC) and differentiated cancer cells (DCC) was used to explain the definitive treatment failure in Glioblastoma Multiforme (GBM) for conventionally and hypo-fractionated treatments. In this study, optimization of temporal dose modulation based on the ODE equation is performed to explore the feasibility of improving GBM treatment outcome.

METHODS

A non-convex optimization problem with the objective of minimizing the total cancer cell number while maintaining the normal tissue biological effective dose (BEDnormal) at 100 Gy, equivalent to the conventional 2 Gy × 30 dosing scheme was formulated. With specified total number of dose fractions and treatment duration, the optimization was performed using a paired simulated annealing algorithm with fractional doses delivered to the CSC and DCC compartments and time intervals between fractions as variables. The recurrence time, defined as the time point at which the total tumor cell number regrows to 2.8×109 cells, was used to evaluate optimization outcome. Optimization was performed for conventional treatment time frames equivalent to currently and historically utilized fractionation schemes, in which limited improvement in recurrence time delay was observed. The efficacy of a super hyperfractionated approach with a prolonged treatment duration of one year was therefore tested, with both fixed regular and optimized variable time intervals between dose fractions corresponding to total number of fractions equivalent to weekly, bi-weekly, and monthly deliveries (n = 53, 27, 13). Optimization corresponding to BEDnormal of 150 Gy was also obtained to evaluate the possibility in further recurrence delay with dose escalation.

RESULTS

For the super hyperfractionated schedules with dose fraction number equivalent to weekly, bi-weekly, and monthly deliveries, the recurrence time points were found to be 430.5, 423.9, and 413.3 days, respectively, significantly delayed compared with the recurrence time of 250.3 days from conventional fractionation. Results show that optimal outcome was achieved by first delivering infrequent fractions followed by dense once per day fractions in the middle and end of the treatment course, with sparse and low dose treatments in the between. The dose to the CSC compartment was held relatively constant throughout while larger dose fractions to the DCC compartment were observed in the beginning and final fractions that preceded large time intervals. Dose escalation to BEDnormal of 150 Gy was shown capable of further delaying recurrence time to 452 days.

CONCLUSION

The development and utilization of a temporal dose fractionation optimization framework in the context of CSC dynamics have demonstrated that substantial delay in GBM local tumor recurrence could be achieved with a super hyperfractionated treatment approach. Preclinical and clinical studies are needed to validate the efficacy of this novel treatment delivery method.

ROAD: ROtational direct Aperture optimization with a Decoupled ring-collimator for FLASH radiotherapy

Ultra-high dose rate in radiotherapy (FLASH) has been shown to increase the therapeutic index with markedly reduced normal tissue toxicity and the same or better tumor cell killing. The challenge to achieve FLASH using x-rays, besides developing a high output linac, is to intensity-modulate the high-dose-rate x-rays so that the biological gain is not offset by the lack of physical dose conformity. In this study, we develop the ROtational direct Aperture optimization with a Decoupled ring-collimator (ROAD) to achieve simultaneous ultrafast delivery and complex dose modulation. The ROAD design includes a fast-rotating slip-ring linac and a decoupled collimator-ring with 75 pre-shaped multi-leaf-collimator (MLC) modules. The ring-source rotates at 1 rotation per second (rps) clockwise while the ring-collimator is either static or rotating at 1 rps counterclockwise, achieving 75 (ROAD-75) or 150 (ROAD-150) equal-angular beams for one full arc. The Direct Aperture Optimization (DAO) for ROAD was formulated to include a least-square dose fidelity, an anisotropic total variation term, and a single segment term. The FLASH dose (FD) and FLASH biological equivalent dose (FBED) were computed voxelwise, with the latter using a spatiotemporal model accounting for radiolytic oxygen depletion. ROAD was compared with clinical volumetric modulated arc therapy (VMAT) on a brain, a lung, a prostate, and a head and neck cancer patient. The mean dose rate of ROAD-75 and ROAD-150 are 76.2 Gy s^-1 and 112 Gy s^-1 respectively to deliver 25 Gy single-fraction dose in 1 s. With improved PTV homogeneity, ROAD-150 reduced (max, mean) OAR physical dose by (4.8 Gy, 6.3 Gy). The average R50 and integral dose of (VMAT, ROAD-75, ROAD-150) are (4.8, 3.2, 3.2) and (89, 57, 56) Gy×Liter, respectively. The FD and FBED showed model dependent FLASH effects. The novel ROAD design achieves ultrafast dose delivery and improves physical dosimetry compared with clinical VMAT, providing a potentially viable engineering solution for x-ray FLASH radiotherapy.

DeepMC: a deep learning method for efficient Monte Carlo beamlet dose calculation by predictive denoising in magnetic resonance-guided radiotherapy

Emerging magnetic resonance (MR) guided radiotherapy affords significantly improved anatomy visualization and, subsequently, more effective personalized treatment. The new therapy paradigm imposes significant demands on radiation dose calculation quality and speed, creating an unmet need for the acceleration of Monte Carlo (MC) dose calculation. Existing deep learning approaches to denoise the final plan MC dose fail to achieve the accuracy and speed requirements of large-scale beamlet dose calculation in the presence of a strong magnetic field for online adaptive radiotherapy planning. Our deep learning dose calculation method, DeepMC, addresses these needs by predicting low-noise dose from extremely noisy (but fast) MC-simulated dose and anatomical inputs, thus enabling significant acceleration. DeepMC simultaneously reduces MC sampling noise and predicts corrupted dose buildup at tissue-air material interfaces resulting from MR-field induced electron return effects. Here we demonstrate our model's ability to accelerate dose calculation for daily treatment planning by a factor of 38 over traditional low-noise MC simulation with clinically meaningful accuracy in deliverable dose and treatment delivery parameters. As a post-processing approach, DeepMC provides compounded acceleration of large-scale dose calculation when used alongside established MC acceleration techniques in variance reduction and graphics processing unit-based MC simulation.

Multi-task edge-recalibrated network for male pelvic multi-organ segmentation on CT images

Automated male pelvic multi-organ segmentation on CT images is highly desired for applications, including radiotherapy planning. To further improve the performance and efficiency of existing automated segmentation methods, in this study, we propose a multi-task edge-recalibrated network (MTER-Net), which aims to overcome the challenges, including blurry boundaries, large inter-patient appearance variations, and low soft-tissue contrast. The proposed MTER-Net is equipped with the following novel components. (a) To exploit the saliency and stability of femoral heads, we employed a light-weight localization module to locate the target region and efficiently remove the complex background. (b) We add an edge stream to the regular segmentation stream to focus on processing the edge-related information, distinguish the organs with blurry boundaries, and then boost the overall segmentation performance. Between the regular segmentation stream and edge stream, we introduce an edge recalibration module at each resolution level to connect the intermediate layers and deliver the higher-level activations from the regular stream to the edge stream to denoise the irrelevant activations. (c) Finally, using a 3D Atrous Spatial Pyramid Pooling (ASPP) feature fusion module, we fuse the features at different scales in the regular stream and the predictions from the edge stream to form the final segmentation result. The proposed segmentation network was evaluated on 200 prostate cancer patient CT images with manually delineated contours of bladder, rectum, seminal vesicle, and prostate. The segmentation performance of the proposed method was quantitatively evaluated using three metrics including Dice similarity coefficient (DSC), average surface distance (ASD), and 95% surface distance (95SD). The proposed MTER-Net achieves average DSC of 86.35%, ASD of 1.09 mm, and 95SD of 3.53 mm on the four organs, which outperforms the state-of-the-art segmentation networks by a large margin. Specifically, the quantitative DSC evaluation results of the four organs are 96.49% (bladder), 86.39% (rectum), 76.38% (seminal vesicle), and 86.14% (prostate), respectively. In conclusion, we demonstrate that the proposed MTER-Net efficiently attains superior performance to state-of-the-art pelvic organ segmentation methods.

Automatic detection and segmentation of multiple brain metastases on magnetic resonance image using asymmetric UNet architecture

Detection of brain metastases is a paramount task in cancer management due both to the number of high-risk patients and the difficulty of achieving consistent detection. In this study, we aim to improve the accuracy of automated brain metastasis (BM) detection methods using a novel asymmetric UNet (asym-UNet) architecture. An end-to-end asymmetric 3D-UNet architecture, with two down-sampling arms and one up-sampling arm, was constructed to capture the imaging features. The two down-sampling arms were trained using two different kernels (3 × 3 × 3 and 1 × 1 × 3, respectively) with the kernel (1 × 1 × 3) dominating the learning. As a comparison, vanilla single 3D UNets were trained with different kernels and evaluated using the same datasets. Voxel-based Dice similarity coefficient (DSC_v), sensitivity (S _v), precision (P _v), BM-based sensitivity (S _BM), and false detection rate (F _BM) were used to evaluate model performance. Contrast-enhanced T1 MR images from 195 patients with a total of 1034 BMs were solicited from our institutional stereotactic radiosurgery database. The patient cohort was split into training (160 patients, 809 lesions), validation (20 patients, 136 lesions), and testing (15 patients, 89 lesions) datasets. The lesions in the testing dataset were further divided into two subgroups based on the diameters (small S = 1-10 mm, large L = 11-26 mm). In the testing dataset, there were 72 and 17 BMs in the S and L sub-groups, respectively. Among all trained networks, asym-UNet achieved the highest DSC_v of 0.84 and lowest F _BM of 0.24. Although vanilla 3D-UNet with a single 1 × 1 × 3 kernel achieved the highest sensitivities for the S group, it resulted in the lowest precision and highest false detection rate. Asym-UNet was shown to balance sensitivity and false detection rate as well as keep the segmentation accuracy high. The novel asym-UNet segmentation network showed overall competitive segmentation performance and more pronounced improvement in hard-to-detect small BMs comparing to the vanilla single 3D UNet.