Deep learning based time-to-event analysis with PET, CT and joint PET/CT for head and neck cancer prognosis

OBJECTIVES

Recent studies have shown that deep learning based on pre-treatment positron emission tomography (PET) or computed tomography (CT) is promising for distant metastasis (DM) and overall survival (OS) prognosis in head and neck cancer (HNC). However, lesion segmentation is typically required, resulting in a predictive power susceptible to variations in primary and lymph node gross tumor volume (GTV) segmentation. This study aimed at achieving prognosis without GTV segmentation, and extending single modality prognosis to joint PET/CT to allow investigating the predictive performance of combined- compared to single-modality inputs.

METHODS

We employed a 3D-Resnet combined with a time-to-event outcome model to incorporate censoring information. We focused on the prognosis of DM and OS for HNC patients. For each clinical endpoint, five models with PET and/or CT images as input were compared: PET-GTV, PET-only, CT-GTV, CT-only, and PET/CT-GTV models, where -GTV indicates that the corresponding images were masked using the GTV contour. Publicly available delineated CT and PET scans from 4 different Canadian hospitals (293) and the MAASTRO clinic (74) were used for training by 3-fold cross-validation (CV). For independent testing, we used 110 patients from a collaborating institution. The predictive performance was evaluated via Harrell's Concordance Index (HCI) and Kaplan-Meier curves.

RESULTS

In a 5-year time-to-event analysis, all models could produce CV HCIs with median values around 0.8 for DM and 0.7 for OS. The best performance was obtained with the PET-only model, achieving a median testing HCI of 0.82 for DM and 0.69 for OS. Compared with the PET/CT-GTV model, the PET-only still had advantages of up to 0.07 in terms of testing HCI. The Kaplan-Meier curves and corresponding log-rank test results also demonstrated significant stratification capability of our models for the testing cohort.

CONCLUSION

Deep learning-based DM and OS time-to-event models showed predictive capability and could provide indications for personalized RT. The best predictive performance achieved by the PET-only model suggested GTV segmentation might be less relevant for PET-based prognosis.

Business exit during the COVID-19 pandemic: Non-traditional measures in historical context

Lags in official data releases have forced economists and policymakers to leverage "alternative" or "non-traditional" data to measure business exit resulting from the COVID-19 pandemic. We first review official data on business exit in recent decades to place the alternative measures of exit within historical context. For the U.S., business exit is fairly common, with about 7.5 percent of firms exiting annually in recent years. The high level of exit is driven by very small firms and establishments. We then explore a range of alternative measures of business exit, including novel measures based on paycheck issuance and phone-tracking data, which indicate exit was elevated in certain sectors during the first year of the pandemic. That said, we find many industries have likely seen lower-than-usual exit rates, and exiting businesses do not appear to represent a large share of U.S. employment. As a result, exit appears lower than widespread expectations from early in the pandemic.

Offline and online LSTM networks for respiratory motion prediction in MR-guided radiotherapy

Objective. Gated beam delivery is the current clinical practice for respiratory motion compensation in MR-guided radiotherapy, and further research is ongoing to implement tracking. To manage intra-fractional motion using multileaf collimator tracking the total system latency needs to be accounted for in real-time. In this study, long short-term memory (LSTM) networks were optimized for the prediction of superior–inferior tumor centroid positions extracted from clinically acquired 2D cine MRIs. Approach. We used 88 patients treated at the University Hospital of the LMU Munich for training and validation (70 patients, 13.1 h), and for testing (18 patients, 3.0 h). Three patients treated at Fondazione Policlinico Universitario Agostino Gemelli were used as a second testing set (1.5 h). The performance of the LSTMs in terms of root mean square error (RMSE) was compared to baseline linear regression (LR) models for forecasted time spans of 250 ms, 500 ms and 750 ms. Both the LSTM and the LR were trained with offline (offline LSTM and offline LR) and online schemes (offline+online LSTM and online LR), the latter to allow for continuous adaptation to recent respiratory patterns. Main results. We found the offline+online LSTM to perform best for all investigated forecasts. Specifically, when predicting 500 ms ahead it achieved a mean RMSE of 1.20 mm and 1.00 mm, while the best performing LR model achieved a mean RMSE of 1.42 mm and 1.22 mm for the LMU and Gemelli testing set, respectively. Significance. This indicates that LSTM networks have potential as respiratory motion predictors and that continuous online re-optimization can enhance their performance.

Evaluation of an anthropomorphic ion chamber and 3D gel dosimetry head phantom at a 0.35 T MR-linac using separate 1.5 T MR-scanners for gel readout

PURPOSE

To date, no universally accepted technique for the evaluation of the overall dosimetric performance of hybrid integrated magnetic resonance imaging (MR) - linear accelerators (linacs) is available. We report on the suitability and reliability of a novel phantom with modular inserts for combined polymer gel (PG) and ionisation chamber (IC) measurements at a 0.35 T MR-linac.

METHODS

Three 3D-printed, modular head phantoms, based on real patient anatomy, were used for repeated (2 times) PG irradiations of cranial treatment plans on a 0.35 T MR-linac. The PG readout was performed on two 1.5 T diagnostic MR-scanners to reduce scanning time. The PG dose volumes were normalised to the IC dose (normalised dose N1) and to the median planning target volume dose (normalised dose N2). Linearity of the PG dose response was validated and dose profiles, centres of mass (COM) of the 95% isodoses and dose volume histograms (DVH) were compared between planned and measured dose distributions and a 3D gamma analysis was performed.

RESULTS

Dose linearity of the PG was good (R²> 0.99 for all linear fit functions). High agreement was found between planned and measured dose volumes in the dose profiles and DVHs. The largest dose deviation was found in the intermediate dose region (mean dose deviation 0.2Gy; 5.6%). A mean COM offset of 1.2mm indicated high spatial accuracy. Mean 3D gamma passing rates (2%, 2mm) of 83.3% for N1 and 91.6% for N2 dose distributions were determined. When comparing repeated PG measurements to each other, a mean gamma passing rate of 95.7% was found.

CONCLUSION

The new modular phantom was found practical for use at a 0.35 T MR-linac. In contrast to the high dose region, larger mean deviations were found in the mid dose range. The PG measurements showed high reproducibility. The MR-linac performed well in a non-adaptive setting in terms of spatial and dosimetric accuracy.

Dosimetric benefit of MR-guided online adaptive radiotherapy in different tumor entities: liver, lung, abdominal lymph nodes, pancreas and prostate

Background

Hybrid magnetic resonance (MR)-Linac systems have recently been introduced into clinical practice. The systems allow online adaption of the treatment plan with the aim of compensating for interfractional anatomical changes. The aim of this study was to evaluate the dose volume histogram (DVH)-based dosimetric benefits of online adaptive MR-guided radiotherapy (oMRgRT) across different tumor entities and to investigate which subgroup of plans improved the most from adaption.

Methods

Fifty patients treated with oMRgRT for five different tumor entities (liver, lung, multiple abdominal lymph nodes, pancreas, and prostate) were included in this retrospective analysis. Various target volume (gross tumor volume GTV, clinical target volume CTV, and planning target volume PTV) and organs at risk (OAR) related DVH parameters were compared between the dose distributions before and after plan adaption.

Results

All subgroups clearly benefited from online plan adaption in terms of improved PTV coverage. For the liver, lung and abdominal lymph nodes cases, a consistent improvement in GTV coverage was found, while many fractions of the prostate subgroup showed acceptable CTV coverage even before plan adaption. The largest median improvements in GTV near-minimum dose (D_98%) were found for the liver (6.3%, p < 0.001), lung (3.9%, p < 0.001), and abdominal lymph nodes (6.8%, p < 0.001) subgroups. Regarding OAR sparing, the largest median OAR dose reduction during plan adaption was found for the pancreas subgroup (-87.0%). However, in the pancreas subgroup an optimal GTV coverage was not always achieved because sparing of OARs was prioritized.

Conclusion

With online plan adaptation, it was possible to achieve significant improvements in target volume coverage and OAR sparing for various tumor entities and account for interfractional anatomical changes.

Dosimetric impact of deep learning-based CT auto-segmentation on radiation therapy treatment planning for prostate cancer

Background

The evaluation of automatic segmentation algorithms is commonly performed using geometric metrics. An analysis based on dosimetric parameters might be more relevant in clinical practice but is often lacking in the literature. The aim of this study was to investigate the impact of state-of-the-art 3D U-Net-generated organ delineations on dose optimization in radiation therapy (RT) for prostate cancer patients.

Methods

A database of 69 computed tomography images with prostate, bladder, and rectum delineations was used for single-label 3D U-Net training with dice similarity coefficient (DSC)-based loss. Volumetric modulated arc therapy (VMAT) plans have been generated for both manual and automatic segmentations with the same optimization settings. These were chosen to give consistent plans when applying perturbations to the manual segmentations. Contours were evaluated in terms of DSC, average and 95% Hausdorff distance (HD). Dose distributions were evaluated with the manual segmentation as reference using dose volume histogram (DVH) parameters and a 3%/3 mm gamma-criterion with 10% dose cut-off. A Pearson correlation coefficient between DSC and dosimetric metrics, i.e. gamma index and DVH parameters, has been calculated.

Results

3D U-Net-based segmentation achieved a DSC of 0.87 (0.03) for prostate, 0.97 (0.01) for bladder and 0.89 (0.04) for rectum. The mean and 95% HD were below 1.6 (0.4) and below 5 (4) mm, respectively. The DVH parameters, V\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_{60/65/70\,{\mathrm{Gy}}}$$\end{document}60/65/70Gy for the bladder and V\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_{50/65/70\,{\mathrm{Gy}}}$$\end{document}50/65/70Gy for the rectum, showed agreement between dose distributions within \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pm \, 5\%$$\end{document}±5% and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pm \,2\%$$\end{document}±2%, respectively. The D\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_{98/2\%}$$\end{document}98/2% and V\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_{95\%}$$\end{document}95%, for prostate and its 3 mm expansion (surrogate clinical target volume) showed agreement with the reference dose distribution within 2% and 3 Gy with the exception of one case. The average gamma pass-rate was 85%. The comparison between geometric and dosimetric metrics showed no strong statistically significant correlation.

Conclusions

The 3D U-Net developed for this work achieved state-of-the-art geometrical performance. Analysis based on clinically relevant DVH parameters of VMAT plans demonstrated neither excessive dose increase to OARs nor substantial under/over-dosage of the target in all but one case. Yet the gamma analysis indicated several cases with low pass rates. The study highlighted the importance of adding dosimetric analysis to the standard geometric evaluation.

Risk Stratification Using 18F-FDG PET/CT and Artificial Neural Networks in Head and Neck Cancer Patients Undergoing Radiotherapy

This study retrospectively analyzed the performance of artificial neural networks (ANN) to predict overall survival (OS) or locoregional failure (LRF) in HNSCC patients undergoing radiotherapy, based on 2-[¹⁸F]FDG PET/CT and clinical covariates. We compared predictions relying on three different sets of features, extracted from 230 patients. Specifically, (i) an automated feature selection method independent of expert rating was compared with (ii) clinical variables with proven influence on OS or LRF and (iii) clinical data plus expert-selected SUV metrics. The three sets were given as input to an artificial neural network for outcome prediction, evaluated by Harrell’s concordance index (HCI) and by testing stratification capability. For OS and LRF, the best performance was achieved with expert-based PET-features (0.71 HCI) and clinical variables (0.70 HCI), respectively. For OS stratification, all three feature sets were significant, whereas for LRF only expert-based PET-features successfully classified low vs. high-risk patients. Based on 2-[¹⁸F]FDG PET/CT features, stratification into risk groups using ANN for OS and LRF is possible. Differences in the results for different feature sets confirm the relevance of feature selection, and the key importance of expert knowledge vs. automated selection.

Validation of proton dose calculation on scatter corrected 4D cone beam computed tomography using a porcine lung phantom

Proton therapy treatment for lungs remains challenging as images enabling the detection of inter- and intra-fractional motion, which could be used for proton dose adaptation, are not readily available. 4D computed tomography (4DCT) provides high image quality but is rarely available in-room, while in-room 4D cone beam computed tomography (4DCBCT) suffers from image quality limitations stemming mostly from scatter detection. This study investigated the feasibility of using virtual 4D computed tomography (4DvCT) as a prior for a phase-per-phase scatter correction algorithm yielding a 4D scatter corrected cone beam computed tomography image (4DCBCT_cor), which can be used for proton dose calculation. 4DCT and 4DCBCT scans of a porcine lung phantom, which generated reproducible ventilation, were acquired with matching breathing patterns. Diffeomorphic Morphons, a deformable image registration algorithm, was used to register the mid-position 4DCT to the mid-position 4DCBCT and yield a 4DvCT. The 4DCBCT was reconstructed using motion-aware reconstruction based on spatial and temporal regularization (MA-ROOSTER). Successively for each phase, digitally reconstructed radiographs of the 4DvCT, simulated without scatter, were exploited to correct scatter in the corresponding CBCT projections. The 4DCBCT_corwas then reconstructed with MA-ROOSTER using the corrected CBCT projections and the same settings and deformation vector fields as those already used for reconstructing the 4DCBCT. The 4DCBCT_corand the 4DvCT were evaluated phase-by-phase, performing proton dose calculations and comparison to those of a ground truth 4DCT by means of dose-volume-histograms (DVH) and gamma pass-rates (PR). For accumulated doses, DVH parameters deviated by at most 1.7% in the 4DvCT and 2.0% in the 4DCBCT_corcase. The gamma PR for a (2%, 2 mm) criterion with 10% threshold were at least 93.2% (4DvCT) and 94.2% (4DCBCT_cor), respectively. The 4DCBCT_cortechnique enabled accurate proton dose calculation, which indicates the potential for applicability to clinical 4DCBCT scans.

Dosimetric comparison of MR-linac-based IMRT and conventional VMAT treatment plans for prostate cancer

Background

The aim of this study was to evaluate and compare the performance of intensity modulated radiation therapy (IMRT) plans, planned for low-field strength magnetic resonance (MR) guided linear accelerator (linac) delivery (labelled IMRT MRL plans), and clinical conventional volumetric modulated arc therapy (VMAT) plans, for the treatment of prostate cancer (PCa). Both plans used the original planning target volume (PTV) margins. Additionally, the potential dosimetric benefits of MR-guidance were estimated, by creating IMRT MRL plans using smaller PTV margins.

Materials and methods

20 PCa patients previously treated with conventional VMAT were considered. For each patient, two different IMRT MRL plans using the low-field MR-linac treatment planning system were created: one with original (orig.) PTV margins and the other with reduced (red.) PTV margins. Dose indices related to target coverage, as well as dose-volume histogram (DVH) parameters for the target and organs at risk (OAR) were compared. Additionally, the estimated treatment delivery times and the number of monitor units (MU) of each plan were evaluated.

Results

The dose distribution in the high dose region and the target volume DVH parameters (D_98%, D_50%, D_2% and V_95%) were similar for all three types of treatment plans, with deviations below 1% in most cases. Both IMRT MRL plans (orig. and red. PTV margins) showed similar homogeneity indices (HI), however worse values for the conformity index (CI) were also found when compared to VMAT. The IMRT MRL plans showed similar OAR sparing when the orig. PTV margins were used but a significantly better sparing was feasible when red. PTV margins were applied. Higher number of MU and longer predicted treatment delivery times were seen for both IMRT MRL plans.

Conclusions

A comparable plan quality between VMAT and IMRT MRL plans was achieved, when applying the same PTV margin. However, online MR-guided adaptive radiotherapy allows for a reduction of PTV margins. With a red. PTV margin, better sparing of the surrounding tissues can be achieved, while maintaining adequate target coverage. Nonetheless, longer treatment delivery times, characteristic for the IMRT technique, have to be expected.

Measurement-based range evaluation for quality assurance of CBCT-based dose calculations in adaptive proton therapy

PURPOSE

The implementation of volumetric in-room imaging for online adaptive radiotherapy makes extensive testing of this image data for treatment planning necessary. Especially for proton beams the higher sensitivity to stopping power properties of the tissue results in more stringent requirements. Current approaches mainly focus on recalculation of the plans on the new image data, lacking experimental verification, and ignoring the impact on the plan re-optimization process. The aim of this study was to use gel and film dosimetry coupled with a three-dimensional (3D) printed head phantom (based on the planning CT of the patient) for 3D range verification of intensity-corrected cone beam computed tomography (CBCT) image data for adaptive proton therapy.

METHODS

Single field uniform dose pencil beam scanning proton plans were optimized for three different patients on the patients' planning CT (planCT) and the patients' intensity-corrected CBCT (scCBCT) for the same target volume using the same optimization constraints. The CBCTs were corrected on projection level using the planCT as a prior. The dose optimized on planCT and recalculated on scCBCT was compared in terms of proton range differences (80% distal fall-off, recalculation). Moreover, the dose distribution resulting from recalculation of the scCBCT-optimized plan on the planCT and the original planCT dose distribution were compared (simulation). Finally, the two plans of each patient were irradiated on the corresponding patient-specific 3D printed head phantom using gel dosimetry inserts for one patient and film dosimetry for all three patients. Range differences were extracted from the measured dose distributions. The measured and the simulated range differences were corrected for range differences originating from the initial plans and evaluated.

RESULTS

The simulation approach showed high agreement with the standard recalculation approach. The median values of the range differences of these two methods agreed within 0.1 mm and the interquartile ranges (IQRs) within 0.3 mm for all three patients. The range differences of the film measurement were accurately matching with the simulation approach in the film plane. The median values of these range differences deviated less than 0.1 mm and the IQRs less than 0.4 mm. For the full 3D evaluation of the gel range differences, the median value and IQR matched those of the simulation approach within 0.7 and 0.5 mm, respectively. scCBCT- and planCT-based dose distributions were found to have a range agreement better than 3 mm (median and IQR) for all considered scenarios (recalculation, simulation, and measurement).

CONCLUSIONS

The results of this initial study indicate that an online adaptive proton workflow based on scatter-corrected CBCT image data for head irradiations is feasible. The novel presented measurement- and simulation-based method was shown to be equivalent to the standard literature recalculation approach. Additionally, it has the capability to catch effects of image differences on the treatment plan optimization. This makes the measurement-based approach particularly interesting for quality assurance of CBCT-based online adaptive proton therapy. The observed uncertainties could be kept within those of the registration and positioning. The proposed validation could also be applied for other alternative in-room images, e.g. for MR-based pseudoCTs.

Porcine lung phantom-based validation of estimated 4D-MRI using orthogonal cine imaging for low-field MR-Linacs

Real-time motion monitoring of lung tumors with low-field magnetic resonance imaging-guided linear accelerators (MR-Linacs) is currently limited to sagittal 2D cine magnetic resonance imaging (MRI). To provide input data for improved intrafractional and interfractional adaptive radiotherapy, the 4D anatomy has to be inferred from data with lower dimensionality. The purpose of this study was to experimentally validate a previously proposed propagation method that provides continuous time-resolved estimated 4D-MRI based on orthogonal cine MRI for a low-field MR-Linac. Ex vivo porcine lungs were injected with artificial nodules and mounted in a dedicated phantom that allows for the simulation of periodic and reproducible breathing motion. The phantom was scanned with a research version of a commercial 0.35 T MR-Linac. Respiratory-correlated 4D-MRI were reconstructed and served as ground truth images. Series of interleaved orthogonal slices in sagittal and coronal orientation, intersecting the injected targets, were acquired at 7.3 Hz. Estimated 4D-MRI at 3.65 Hz were created in post-processing using the propagation method and compared to the ground truth 4D-MRI. Eight datasets at different breathing frequencies and motion amplitudes were acquired for three porcine lungs. The overall median (95[Formula: see text] percentile) deviation between ground truth and estimated deformation vector fields was 2.3 mm (5.7 mm), corresponding to 0.7 (1.6) times the in-plane imaging resolution (3.5 × 3.5 mm²). Median (95[Formula: see text] percentile) estimated nodule position errors were 1.5 mm (3.8 mm) for nodules intersected by orthogonal slices and 2.1 mm (7.1 mm) for nodules located more than 2 cm away from either of the orthogonal slices. The estimation error depended on the breathing phase, the motion amplitude and the location of the estimated position with respect to the orthogonal slices. By using the propagation method, the 4D motion within the porcine lung phantom could be accurately and robustly estimated. The method could provide valuable information for treatment planning, real-time motion monitoring, treatment adaptation, and post-treatment evaluation of MR-guided radiotherapy treatments.

Dosimetric impact of geometric distortions in an MRI-only proton therapy workflow for lung, liver and pancreas

In a radiation therapy workflow based on Magnetic Resonance Imaging (MRI), dosimetric errors may arise due to geometric distortions introduced by MRI. The aim of this study was to quantify the dosimetric effect of system-dependent geometric distortions in an MRI-only workflow for proton therapy applied at extra-cranial sites. An approach was developed, in which computed tomography (CT) images were distorted using an MRI displacement map, which represented the MR distortions in a spoiled gradient-echo sequence due to gradient nonlinearities and static magnetic field inhomogeneities. A retrospective study was conducted on 4DCT/MRI digital phantoms and 18 4DCT clinical datasets of the thoraco-abdominal site. The treatment plans were designed and separately optimized for each beam in a beam specific Planning Target Volume on the distorted CT, and the final dose distribution was obtained as the average. The dose was then recalculated in undistorted CT using the same beam geometry and beam weights. The analysis was performed in terms of Dose Volume Histogram (DVH) parameters. No clinically relevant dosimetric impact was observed on organs at risk, whereas in the target structure, geometric distortions caused statistically significant variations in the planned dose DVH parameters and dose homogeneity index (DHI). The dosimetric variations in the target structure were smaller in abdominal cases (ΔD_2%, ΔD_98%, and ΔD_mean all below 0.1% and ΔDHI below 0.003) compared to the lung cases. Indeed, lung patients with tumors isolated inside lung parenchyma exhibited higher dosimetric variations (ΔD_2%≥0.3%, ΔD_98%≥15.9%, ΔD_mean≥3.3% and ΔDHI≥0.102) than lung patients with tumor close to soft tissue (ΔD_2%≤0.4%, ΔD_98%≤5.6%, ΔD_mean≤0.9% and ΔDHI≤0.027) potentially due to higher density variations along the beam path. Results suggest the potential applicability of MRI-only proton therapy, provided that specific analysis is applied for isolated lung tumors.

MR-guided proton therapy: a review and a preview

Background

The targeting accuracy of proton therapy (PT) for moving soft-tissue tumours is expected to greatly improve by real-time magnetic resonance imaging (MRI) guidance. The integration of MRI and PT at the treatment isocenter would offer the opportunity of combining the unparalleled soft-tissue contrast and real-time imaging capabilities of MRI with the most conformal dose distribution and best dose steering capability provided by modern PT. However, hybrid systems for MR-integrated PT (MRiPT) have not been realized so far due to a number of hitherto open technological challenges. In recent years, various research groups have started addressing these challenges and exploring the technical feasibility and clinical potential of MRiPT. The aim of this contribution is to review the different aspects of MRiPT, to report on the status quo and to identify important future research topics.

Methods

Four aspects currently under study and their future directions are discussed: modelling and experimental investigations of electromagnetic interactions between the MRI and PT systems, integration of MRiPT workflows in clinical facilities, proton dose calculation algorithms in magnetic fields, and MRI-only based proton treatment planning approaches.

Conclusions

Although MRiPT is still in its infancy, significant progress on all four aspects has been made, showing promising results that justify further efforts for research and development to be undertaken. First non-clinical research solutions have recently been realized and are being thoroughly characterized. The prospect that first prototype MRiPT systems for clinical use will likely exist within the next 5 to 10 years seems realistic, but requires significant work to be performed by collaborative efforts of research groups and industrial partners.

Medical physics challenges in clinical MR-guided radiotherapy

The integration of magnetic resonance imaging (MRI) for guidance in external beam radiotherapy has faced significant research and development efforts in recent years. The current availability of linear accelerators with an embedded MRI unit, providing volumetric imaging at excellent soft tissue contrast, is expected to provide novel possibilities in the implementation of image-guided adaptive radiotherapy (IGART) protocols. This study reviews open medical physics issues in MR-guided radiotherapy (MRgRT) implementation, with a focus on current approaches and on the potential for innovation in IGART.Daily imaging in MRgRT provides the ability to visualize the static anatomy, to capture internal tumor motion and to extract quantitative image features for treatment verification and monitoring. Those capabilities enable the use of treatment adaptation, with potential benefits in terms of personalized medicine. The use of online MRI requires dedicated efforts to perform accurate dose measurements and calculations, due to the presence of magnetic fields. Likewise, MRgRT requires dedicated quality assurance (QA) protocols for safe clinical implementation.Reaction to anatomical changes in MRgRT, as visualized on daily images, demands for treatment adaptation concepts, with stringent requirements in terms of fast and accurate validation before the treatment fraction can be delivered. This entails specific challenges in terms of treatment workflow optimization, QA, and verification of the expected delivered dose while the patient is in treatment position. Those challenges require specialized medical physics developments towards the aim of fully exploiting MRI capabilities. Conversely, the use of MRgRT allows for higher confidence in tumor targeting and organs-at-risk (OAR) sparing.The systematic use of MRgRT brings the possibility of leveraging IGART methods for the optimization of tumor targeting and quantitative treatment verification. Although several challenges exist, the intrinsic benefits of MRgRT will provide a deeper understanding of dose delivery effects on an individual basis, with the potential for further treatment personalization.