Search

Background and Purpose: Radiation-induced pneumonitis (RP), diagnosed 6-12 weeks after treatment, is a complication of lung tumor radiotherapy. So far, clinical and dosimetric parameters have not been reliable in predicting RP. We propose using non-contrast enhanced magnetic resonance imaging (MRI) based functional parameters acquired over the treatment course for patient stratification for improved follow-up., Materials and Methods: 23 lung tumor patients received MR-guided hypofractionated stereotactic body radiation therapy at a 0.35T MR-Linac. Ventilation- and perfusion-maps were generated from 2D-cine MRI-scans acquired after the first and last treatment fraction (Fx) using non-uniform Fourier decomposition. The relative differences in ventilation and perfusion between last and first Fx in three regions (planning target volume (PTV), lung volume receiving more than 20Gy (V20) excluding PTV, whole tumor-bearing lung excluding PTV) and three dosimetric parameters (mean lung dose, V20, mean dose to the gross tumor volume) were investigated. Univariate receiver operating characteristic curve - area under the curve (ROC-AUC) analysis was performed (endpoint RP grade≥1) using 5000 bootstrapping samples. Differences between RP and non-RP patients were tested for statistical significance with the non-parametric Mann-Whitney U test (α=0.05)., Results: 14/23 patients developed RP of grade≥1 within 3 months. The dosimetric parameters showed no significant differences between RP and non-RP patients. In contrast, the functional parameters, especially the relative ventilation difference in the PTV, achieved a p-value<0.05 and an AUC value of 0.84., Conclusion: MRI-based functional parameters extracted from 2D-cine MRI-scans were found to be predictive of RP development in lung tumor patients., Competing Interests: Declaration of competing interest The Department of Radiation Oncology of the LMU University Hospital, Munich has research agreements with Elekta, Brainlab and C-RAD., (Copyright © 2024 The Authors. Published by Elsevier B.V. All rights reserved.)

Introduction: This study evaluated the performance of magnetic resonance thermometry (MRT) during deep-regional hyperthermia (HT) in pelvic and lower-extremity soft-tissue sarcomas., Materials and Methods: 17 pelvic (45 treatments) and 16 lower-extremity (42 treatments) patients underwent standard regional HT and chemotherapy. Pairs of double-echo gradient-echo scans were acquired during the MR protocol 1.4 s apart. For each pair, precision was quantified using phase data from both echoes ('dual-echo') or only one ('single-echo') in- or excluding body fat pixels in the field drift correction region of interest. The precision of each method was compared to that of the MRT approach using a built-in clinical software tool (SigmaVision). Accuracy was assessed in three lower-extremity patients (six treatments) using interstitial temperature probes. The Jaccard coefficient quantified pretreatment motion; receiver operating characteristic analysis assessed its predictability for acceptable precision (<1 °C) during HT., Results: Compared to the built-in dual-echo approach, single-echo thermometry improved the mean temporal precision from 1.32 ± 0.40 °C to 1.07 ± 0.34 °C (pelvis) and from 0.99 ± 0.28 °C to 0.76 ± 0.23 °C (lower extremities). With body fat-based field drift correction, single-echo mean accuracy improved from 1.4 °C to 1.0 °C. Pretreatment bulk motion provided excellent precision prediction with an area under the curve of 0.80-0.86 (pelvis) and 0.81-0.83 (lower extremities), compared to gastrointestinal air motion (0.52-0.58)., Conclusion: Single-echo MRT exhibited better precision than dual-echo MRT. Body fat-based field-drift correction significantly improved MRT accuracy. Pretreatment bulk motion showed improved prediction of acceptable MRT temporal precision over gastrointestinal air motion.

Background: Hybrid devices that combine radiation therapy and MR-imaging have been introduced in the clinical routine for the treatment of lung cancer. This opened up not only possibilities in terms of accurate tumor tracking, dose delivery and adapted treatment planning, but also functional lung imaging. The aim of this study was to show the feasibility of Non-uniform Fourier Decomposition (NuFD) MRI at a 0.35 T MR-Linac as a potential treatment response assessment tool, and propose two signal normalization strategies for enhancing the reproducibility of the results., Methods: Ten healthy volunteers (median age 28 ± 8 years, five female, five male) were repeatedly scanned at a 0.35 T MR-Linac using an optimized 2D+t balanced steady-state free precession (bSSFP) sequence for two coronal slice positions. Image series were acquired in normal free breathing with breaks inside and outside the scanner as well as deep and shallow breathing. Ventilation- and perfusion-weighted maps were generated for each image series using NuFD. For intra-volunteer ventilation map reproducibility, a normalization factor was defined based on the linear correlation of the ventilation signal and diaphragm position of each scan as well as the diaphragm motion amplitude of a reference scan. This allowed for the correction of signal dependency on the diaphragm motion amplitude, which varies with breathing patterns. The second strategy, which can be used for ventilation and perfusion, eliminates the dependency on the signal amplitude by normalizing the ventilation/perfusion maps with the average ventilation/perfusion signal within a selected region-of-interest (ROI). The position and size dependency of this ROI was analyzed. To evaluate the performance of both approaches, the normalized ventilation/perfusion-weighted maps were compared and the deviation of the mean ventilation/perfusion signal from the reference was calculated for each scan. Wilcoxon signed-rank tests were performed to test whether the normalization methods can significantly improve the reproducibility of the ventilation/perfusion maps., Results: The ventilation- and perfusion-weighted maps generated with the NuFD algorithm demonstrated a mostly homogenous distribution of signal intensity as expected for healthy volunteers regardless of the breathing maneuver and slice position. Evaluation of the ROI's size and position dependency showed small differences in the performance. Applying both normalization strategies improved the reproducibility of the ventilation by reducing the median deviation of all scans to 9.1%, 5.7% and 8.6% for the diaphragm-based, the best and worst performing ROI-based normalization, respectively, compared to 29.5% for the non-normalized scans. The significance of this improvement was confirmed by the Wilcoxon signed rank test with [Formula: see text] at [Formula: see text]. A comparison of the techniques against each other revealed a significant difference in the performance between best ROI-based normalization and worst ROI ([Formula: see text]) and between best ROI-based normalization and scaling factor ([Formula: see text]), but not between scaling factor and worst ROI ([Formula: see text]). Using the ROI-based approach for the perfusion-maps, the uncorrected deviation of 10.2% was reduced to 5.3%, which was shown to be significant ([Formula: see text])., Conclusions: Using NuFD for non-contrast enhanced functional lung MRI at a 0.35 T MR-Linac is feasible and produces plausible ventilation- and perfusion-weighted maps for volunteers without history of chronic pulmonary diseases utilizing different breathing patterns. The reproducibility of the results in repeated scans significantly benefits from the introduction of the two normalization strategies, making NuFD a potential candidate for fast and robust early treatment response assessment of lung cancer patients during MR-guided radiotherapy., (© 2023. The Author(s).)

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., (© 2022. The Author(s).)

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.

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.

Purpose: In photon radiotherapy, respiratory-induced target motion can be accounted for by internal target volumes (ITV) or mid-ventilation target volumes (midV) defined on the basis of four-dimensional computed tomography (4D-CT). Intrinsic limitations of these approaches can result in target volumes that are not representative for the gross tumor volume (GTV) motion over the course of treatment. To address these limitations, we propose a novel patient-specific ITV definition method based on real-time 4D magnetic resonance imaging (rt-4DMRI)., Methods: Three lung cancer patients underwent weekly rt-4DMRI scans. A total of 24 datasets were included in this retrospective study. The GTV was contoured on breath-hold MR images and propagated to all rt-4DMRI images by deformable image registration. Different targets were created for the first (reference) imaging sessions: ITVs encompassing all GTV positions over the complete (ITV 80 s ) or partial acquisition time ( ITV 10 s ), ITVs including only voxels with a GTV probability-of-presence (POP) of at least 5% ( ITV 5 % ) or 10% ( ITV 10 % ), and the mid-ventilation GTV position. Reference planning target volumes ( PTV r ) were created by adding margins around the ITVs and midV target volumes. The geometrical overlap of the PTV r with ITV n 5 % from the six to eight subsequent imaging sessions on days n was quantified in terms of the Dice similarity coefficient (DSC), sensitivity [SE: ( PTV r ∩ ITV n 5 % )/ ITV n 5 % ] and precision [PRE: ( PTV r ∩ ITV n 5 % )/ PTV r ] as surrogates for target coverage and normal tissue sparing., Results: Patient-specific analysis yielded a high variance of the overlap values of PTV r 10 s , when different periods within the reference imaging session were sampled. The mid-ventilation-based PTVs were smaller than the ITV-based PTVs. While the SE was high for patients with small breathing pattern variations, changes of the median breathing amplitudes in different imaging sessions led to inferior SE values for the mid-ventilation PTV for one patient. In contrast, PTV r 5 % and PTV r 10 % showed higher SE values with a higher robustness against interfractional changes, at the cost of larger target volumes., Conclusions: The results indicate that rt-4DMRI could be valuable for the definition of target volumes based on the GTV POP to achieve a higher robustness against interfractional changes than feasible with today's 4D-CT-based target definition concepts., (© 2020 The Authors. Medical Physics published by Wiley Periodicals, Inc. on behalf of American Association of Physicists in Medicine.)

Introduction: The recent developments of magnetic resonance (MR) based adaptive strategies for photon and, potentially for proton therapy, require a fast and reliable conversion of MR images to X-ray computed tomography (CT) values. CT values are needed for photon and proton dose calculation. The improvement of conversion results employing a 3D deep learning approach is evaluated. Material and methods: A database of 89 T1-weighted MR head scans with about 100 slices each, including rigidly registered CTs, was created. Twenty-eight validation patients were randomly sampled, and four patients were selected for application. The remaining patients were used to train a 2D and a 3D U-shaped convolutional neural network (Unet). A stack size of 32 slices was used for 3D training. For all application cases, volumetric modulated arc therapy photon and single-field uniform dose pencil-beam scanning proton plans at four different gantry angles were optimized for a generic target on the CT and recalculated on 2D and 3D Unet-based pseudoCTs. Mean (absolute) error (MAE/ME) and a gradient sharpness estimate were used to quantify the image quality. Three-dimensional gamma and dose difference analyses were performed for photon (gamma criteria: 1%, 1 mm) and proton dose distributions (gamma criteria: 2%, 2 mm). Range (80% fall off) differences for beam's eye view profiles were evaluated for protons. Results: Training 36 h for 1000 epochs in 3D (6 h for 200 epochs in 2D) yielded a maximum MAE of 147 HU (135 HU) for the application patients. Except for one patient gamma pass rates for photon and proton dose distributions were above 96% for both Unets. Slice discontinuities were reduced for 3D training at the cost of sharpness. Conclusions: Image analysis revealed a slight advantage of 2D Unets compared to 3D Unets. Similar dose calculation performance was reached for the 2D and 3D network.

A magnetic resonance (MR)-only radiotherapy workflow can reduce cost, radiation exposure and uncertainties introduced by CT-MRI registration. A crucial prerequisite is generating the so called pseudo-CT (pCT) images for accurate dose calculation and planning. Many pCT generation methods have been proposed in the scope of photon radiotherapy. This work aims at verifying for the first time whether a commercially available photon-oriented pCT generation method can be employed for accurate intensity-modulated proton therapy (IMPT) dose calculation. A retrospective study was conducted on ten prostate cancer patients. For pCT generation from MR images, a commercial solution for creating bulk-assigned pCTs, called MR for Attenuation Correction (MRCAT), was employed. The assigned pseudo-Hounsfield Unit (HU) values were adapted to yield an increased agreement to the reference CT in terms of proton range. Internal air cavities were copied from the CT to minimise inter-scan differences. CT- and MRCAT-based dose calculations for opposing beam IMPT plans were compared by gamma analysis and evaluation of clinically relevant target and organ at risk dose volume histogram (DVH) parameters. The proton range in beam's eye view (BEV) was compared using single field uniform dose (SFUD) plans. On average, a [Formula: see text] mm) gamma pass rate of 98.4% was obtained using a [Formula: see text] dose threshold after adaptation of the pseudo-HU values. Mean differences between CT- and MRCAT-based dose in the DVH parameters were below 1 Gy ([Formula: see text]). The median proton range difference was [Formula: see text] mm, with on average 96% of all BEV dose profiles showing a range agreement better than 3 mm. Results suggest that accurate MR-based proton dose calculation using an automatic commercial bulk-assignment pCT generation method, originally designed for photon radiotherapy, is feasible following adaptation of the assigned pseudo-HU values.

Purpose: Magnetic resonance imaging (MRI)-guided radiotherapy enables adaptive treatment plans based on daily anatomical changes and accurate organ visualization. However, the bias field artifact can compromise image quality, affecting diagnostic accuracy and quantitative analyses. This study aims to assess the impact of bias field correction on 0.35 T pelvis MRIs by evaluating clinical anatomy visualization and generative adversarial network (GAN) autosegmentation performance. Materials and methods: 3D simulation MRIs from 60 prostate cancer patients treated on MR-Linac (0.35 T) were collected and preprocessed with the N4ITK algorithm for bias field correction. A 3D GAN architecture was trained, validated, and tested on 40, 10, and 10 patients, respectively, to auto-segment the organs at risk (OARs) rectum and bladder. The GAN was trained and evaluated either with the original or the bias-corrected MRIs. The Dice similarity coefficient (DSC) and 95th percentile Hausdorff distance (HD95th) were computed for the segmented volumes of each patient. The Wilcoxon signed-rank test assessed the statistical difference of the metrics within OARs, both with and without bias field correction. Five radiation oncologists blindly scored 22 randomly chosen patients in terms of overall image quality and visibility of boundaries (prostate, rectum, bladder, seminal vesicles) of the original and bias-corrected MRIs. Bennett's S score and Fleiss' kappa were used to assess the pairwise interrater agreement and the interrater agreement among all the observers, respectively. Results: In the test set, the GAN trained and evaluated on original and biascorrected MRIs showed DSC/HD95th of 0.92/5.63 mm and 0.92/5.91 mm for the bladder and 0.84/10.61 mm and 0.83/9.71 mm for the rectum. No statistical differences in the distribution of the evaluation metrics were found neither for the bladder (DSC: p = 0.07; HD95th: p = 0.35) nor for the rectum (DSC: p = 0.32; HD95th: p = 0.63). From the clinical visual grading assessment, the biascorrected MRI resulted mostly in either no change or an improvement of the image quality and visualization of the organs' boundaries compared with the original MRI. Conclusion: The bias field correction did not improve the anatomy visualization from a clinical point of view and the OARs' auto-segmentation outputs generated by the GAN. [ABSTRACT FROM AUTHOR]

Background: Current commercially available hybrid magnetic resonance linear accelerators (MR‐Linac) use 2D+t cine MR imaging to provide intra‐fractional motion monitoring. However, given the limited temporal resolution of cine MR imaging, target intra‐frame motion deterioration effects, resulting in effective time latency and motion artifacts in the image domain, can be appreciable, especially in the case of fast breathing. Purpose: The aim of this work is to investigate intra‐frame motion deterioration effects in MR‐guided radiotherapy (MRgRT) by simulating the motion‐corrupted image acquisition, and to explore the feasibility of deep‐learning‐based compensation approaches, relying on the intra‐frame motion information which is spatially and temporally encoded in the raw data (k‐space). Methods: An intra‐frame motion model was defined to simulate motion‐corrupted MR images, with 4D anthropomorphic digital phantoms being exploited to provide ground truth 2D+t cine MR sequences. A total number of 10 digital phantoms were generated for lung cancer patients, with randomly selected eight patients for training or validation and the remaining two for testing. The simulation code served as the data generator, and a dedicated motion pattern perturbation scheme was proposed to build the intra‐frame motion database, where three degrees of freedom were designed to guarantee the diversity of intra‐frame motion trajectories, enabling a thorough exploration in the domain of the potential anatomical structure positions. U‐Nets with three types of loss functions: L1 or L2 loss defined in image or Fourier domain, referred to as NNImgLoss‐L1, NNFloss‐L1 and NNL2‐Loss were trained to extract information from the motion‐corrupted image and used to estimate the ground truth final‐position image, corresponding to the end of the acquisition. Images before and after compensation were evaluated in terms of (i) image mean‐squared error (MSE) and mean absolute error (MAE), and (ii) accuracy of gross tumor volume (GTV) contouring, based on optical‐flow image registration. Results: Image degradation caused by intra‐frame motion was observed: for a linearly and fully acquired Cartesian readout k‐space trajectory, intra‐frame motion resulted in an imaging latency of approximately 50% of the acquisition time; in comparison, the motion artifacts exhibited only a negligible contribution to the overall geometric errors. All three compensation models led to a decrease in image MSE/MAE and GTV position offset compared to the motion‐corrupted image. In the investigated testing dataset for GTV contouring, the average dice similarity coefficients (DSC) improved from 88% to 96%, and the 95th percentile Hausdorff distance (HD95) dropped from 4.8 mm to 2.1 mm. Different models showed slight performance variations across different intra‐frame motion amplitude categories: NNImgLoss‐L1 excelled for small/medium amplitudes, whereas NNFloss‐L1 demonstrated higher DSC median values at larger amplitudes. The saliency maps of the motion‐corrupted image highlighted the major contribution of the later acquired k‐space data, as well as the edges of the moving anatomical structures at their final positions, during the model inference stage. Conclusions: Our results demonstrate the deep‐learning‐based approaches have the potential to compensate for intra‐frame motion by utilizing the later acquired data to drive the convergence of the earlier acquired k‐space components. [ABSTRACT FROM AUTHOR]

Background: Magnetic resonance imaging (MRI)‐guided radiotherapy with multileaf collimator (MLC)‐tracking is a promising technique for intra‐fractional motion management, achieving high dose conformality without prolonging treatment times. To improve beam‐target alignment, the geometric error due to system latency should be reduced by using temporal prediction. Purpose: To experimentally compare linear regression (LR) and long‐short‐term memory (LSTM) motion prediction models for MLC‐tracking on an MRI‐linac using multiple patient‐derived traces with different complexities. Methods: Experiments were performed on a prototype 1.0 T MRI‐linac capable of MLC‐tracking. A motion phantom was programmed to move a target in superior‐inferior (SI) direction according to eight lung cancer patient respiratory motion traces. Target centroid positions were localized from sagittal 2D cine MRIs acquired at 4 Hz using a template matching algorithm. The centroid positions were input to one of four motion prediction models. We used (1) a LSTM network which had been optimized in a previous study on patient data from another cohort (offline LSTM). We also used (2) the same LSTM model as a starting point for continuous re‐optimization of its weights during the experiment based on recent motion (offline+online LSTM). Furthermore, we implemented (3) a continuously updated LR model, which was solely based on recent motion (online LR). Finally, we used (4) the last available target centroid without any changes as a baseline (no‐predictor). The predictions of the models were used to shift the MLC aperture in real‐time. An electronic portal imaging device (EPID) was used to visualize the target and MLC aperture during the experiments. Based on the EPID frames, the root‐mean‐square error (RMSE) between the target and the MLC aperture positions was used to assess the performance of the different motion predictors. Each combination of motion trace and prediction model was repeated twice to test stability, for a total of 64 experiments. Results: The end‐to‐end latency of the system was measured to be (389 ± 15) ms and was successfully mitigated by both LR and LSTM models. The offline+online LSTM was found to outperform the other models for all investigated motion traces. It obtained a median RMSE over all traces of (2.8 ± 1.3) mm, compared to the (3.2 ± 1.9) mm of the offline LSTM, the (3.3 ± 1.4) mm of the online LR and the (4.4 ± 2.4) mm when using the no‐predictor. According to statistical tests, differences were significant (p‐value <0.05) among all models in a pair‐wise comparison, but for the offline LSTM and online LR pair. The offline+online LSTM was found to be more reproducible than the offline LSTM and the online LR with a maximum deviation in RMSE between two measurements of 10%. Conclusions: This study represents the first experimental comparison of different prediction models for MRI‐guided MLC‐tracking using several patient‐derived respiratory motion traces. We have shown that among the investigated models, continuously re‐optimized LSTM networks are the most promising to account for the end‐to‐end system latency in MRI‐guided radiotherapy with MLC‐tracking. [ABSTRACT FROM AUTHOR]

Background and purpose: Magnetic resonance imaging guided radiotherapy (MRgRT) offers treatment plan adaptation to the anatomy of the day. In the current MRgRT workflow, this requires the time consuming and repetitive task of manual delineation of organs-at-risk (OARs), which is also prone to inter- and intra-observer variability. Therefore, deep learning autosegmentation (DLAS) is becoming increasingly attractive. No investigation of its application to OARs in thoracic magnetic resonance images (MRIs) from MRgRT has been done so far. This study aimed to fill this gap. Materials and methods: 122 planning MRIs from patients treated at a 0.35 T MR-Linac were retrospectively collected. Using an 80/19/23 (training/validation/test) split, individual 3D U-Nets for segmentation of the left lung, right lung, heart, aorta, spinal canal and esophagus were trained. These were compared to the clinically used contours based on Dice similarity coefficient (DSC) and Hausdorff distance (HD). They were also graded on their clinical usability by a radiation oncologist. Results: Median DSC was 0.96, 0.96, 0.94, 0.90, 0.88 and 0.78 for left lung, right lung, heart, aorta, spinal canal and esophagus, respectively. Median 95th percentile values of the HD were 3.9, 5.3, 5.8, 3.0, 2.6 and 3.5 mm, respectively. The physician preferred the network generated contours over the clinical contours, deeming 85 out of 129 to not require any correction, 25 immediately usable for treatment planning, 15 requiring minor and 4 requiring major corrections. Conclusions: We trained 3D U-Nets on clinical MRI planning data which produced accurate delineations in the thoracic region. DLAS contours were preferred over the clinical contours. [ABSTRACT FROM AUTHOR]

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 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 mm2). Median (95 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. [ABSTRACT FROM AUTHOR]

Simple Summary: Stereotactic body radiotherapy is used in the treatment of liver tumors. However, adjacent organs at risk (OAR) frequently limit the applicable dose to the target volume. The introduction of hybrid magnetic resonance imaging (MRI)-guided radiotherapy systems may allow dose escalation strategies with better OAR sparing due to improved soft tissue visualization, adaptive treatment planning and real-time motion management. Here we report the feasibility and early results of online adaptive MR-guided radiotherapy of primary and secondary liver tumors in eleven patients. The treatment was feasible and successfully completed in all patients. After a median follow-up of five months, no local failure occurred and no ≥ grade 2 toxicity was observed. The technique should be compared to conventional SBRT in further studies to assess the advantages of the technique. Purpose: To assess the feasibility and early results of online adaptive MR-guided radiotherapy (oMRgRT) of liver tumors. Methods: We retrospectively examined consecutive patients with primary or secondary liver lesions treated at our institution using a 0.35T hybrid MR-Linac (Viewray Inc., Mountain View, CA, USA). Online-adaptive treatment planning was used to account for interfractional anatomical changes, and real-time intrafractional motion management using online 2D cine MRI was performed using a respiratory gating approach. Treatment response and toxicity were assessed during follow-up. Results: Eleven patients and a total of 15 lesions were evaluated. Histologies included cholangiocarcinomas and metastases of neuroendocrine tumors, colorectal carcinomas, sarcomas and a gastrointestinal stroma tumor. The median BED10 of the PTV prescription doses was 84.4 Gy (range 59.5–112.5 Gy) applied in 3–5 fractions and the mean GTV BED10 was in median 147.9 Gy (range 71.7–200.5 Gy). Online plan adaptation was performed in 98% of fractions. The median overall treatment duration was 53 min. The treatment was feasible and successfully completed in all patients. After a median follow-up of five months, no local failure occurred and no ≥ grade two toxicity was observed. OMRgRT resulted in better PTV coverage and fewer OAR constraint violations. Conclusion: Early results of MR-linac based oMRgRT for the primary and secondary liver tumors are promising. The treatment was feasible in all cases and well tolerated with minimal toxicity. The technique should be compared to conventional SBRT in further studies to assess the advantages of the technique. [ABSTRACT FROM AUTHOR]

• Comparison of three AI algorithms for real-time prediction of future tumor contours. • Usage of clinical cine MRI data from low-field MR-linacs from two institutions. • Prediction times compatible with reported MLC-tracking latencies on MR-linacs. Magnetic resonance imaging guided radiotherapy (MRgRT) with deformable multileaf collimator (MLC) tracking would allow to tackle both rigid displacement and tumor deformation without prolonging treatment. However, the system latency must be accounted for by predicting future tumor contours in real-time. We compared the performance of three artificial intelligence (AI) algorithms based on long short-term memory (LSTM) modules for the prediction of 2D-contours 500 ms into the future. Models were trained (52 patients, 3.1 h of motion), validated (18 patients, 0.6 h) and tested (18 patients, 1.1 h) with cine MRs from patients treated at one institution. Additionally, we used three patients (2.9 h) treated at another institution as second testing set. We implemented 1) a classical LSTM network (LSTM-shift) predicting tumor centroid positions in superior-inferior and anterior-posterior direction which are used to shift the last observed tumor contour. The LSTM-shift model was optimized both in an offline and online fashion. We also implemented 2) a convolutional LSTM model (ConvLSTM) to directly predict future tumor contours and 3) a convolutional LSTM combined with spatial transformer layers (ConvLSTM-STL) to predict displacement fields used to warp the last tumor contour. The online LSTM-shift model was found to perform slightly better than the offline LSTM-shift and significantly better than the ConvLSTM and ConvLSTM-STL. It achieved a 50% Hausdorff distance of 1.2 mm and 1.0 mm for the two testing sets, respectively. Larger motion ranges were found to lead to more substantial performance differences across the models. LSTM networks predicting future centroids and shifting the last tumor contour are the most suitable for tumor contour prediction. The obtained accuracy would allow to reduce residual tracking errors during MRgRT with deformable MLC-tracking. [ABSTRACT FROM AUTHOR]