Search

Background: Longitudinal assessments of apparent diffusion coefficients (ADCs) derived from diffusion-weighted imaging (DWI) during intracranial radiotherapy at magnetic resonance imaging-guided linear accelerators (MR-linacs) could enable early response assessment by tracking tumor diffusivity changes. However, DWI pulse sequences are currently unavailable in clinical practice at low-field MR-linacs. Quantifying the in vivo repeatability of ADC measurements is a crucial step towards clinical implementation of DWI sequences but has not yet been reported on for low-field MR-linacs. This study assessed ADC measurement repeatability in a phantom and in vivo at a 0.35 T MR-linac. Methods: Eleven volunteers and a diffusion phantom were imaged on a 0.35 T MR-linac. Two echo-planar imaging DWI sequence variants, emphasizing high spatial resolution ("highRes") and signal-to-noise ratio ("highSNR"), were investigated. A test–retest study with an intermediate outside-scanner-break was performed to assess repeatability in the phantom and volunteers' brains. Mean ADCs within phantom vials, cerebrospinal fluid (CSF), and four brain tissue regions were compared to literature values. Absolute relative differences of mean ADCs in pre- and post-break scans were calculated for the diffusion phantom, and repeatability coefficients (RC) and relative RC (relRC) with 95% confidence intervals were determined for each region-of-interest (ROI) in volunteers. Results: Both DWI sequence variants demonstrated high repeatability, with absolute relative deviations below 1% for water, dimethyl sulfoxide, and polyethylene glycol in the diffusion phantom. RelRCs were 7% [5%, 12%] (CSF; highRes), 12% [9%, 22%] (CSF; highSNR), 9% [8%, 12%] (brain tissue ROIs; highRes), and 6% [5%, 7%] (brain tissue ROIs; highSNR), respectively. ADCs measured with the highSNR variant were consistent with literature values for volunteers, while smaller mean values were measured for the diffusion phantom. Conversely, the highRes variant underestimated ADCs compared to literature values, indicating systematic deviations. Conclusions: High repeatability of ADC measurements in a diffusion phantom and volunteers' brains were measured at a low-field MR-linac. The highSNR variant outperformed the highRes variant in accuracy and repeatability, at the expense of an approximately doubled voxel volume. The observed high in vivo repeatability confirms the potential utility of DWI at low-field MR-linacs for early treatment response assessment. [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: Stereotactic body radiation therapy (SBRT) of central lung tumors with photon or proton therapy has a risk of increased toxicity. Treatment planning studies comparing accumulated doses for state‐of‐the‐art treatment techniques, such as MR‐guided radiotherapy (MRgRT) and intensity modulated proton therapy (IMPT), are currently lacking. Purpose: We conducted a comparison of accumulated doses for MRgRT, robustly optimized non‐adaptive IMPT, and online adaptive IMPT for central lung tumors. A special focus was set on analyzing the accumulated doses to the bronchial tree, a parameter linked to high‐grade toxicities. Methods: Data of 18 early‐stage central lung tumor patients, treated at a 0.35 T MR‐linac in eight or five fractions, were analyzed. Three gated treatment scenarios were compared: (S1) online adaptive MRgRT, (S2) non‐adaptive IMPT, and (S3) online adaptive IMPT. The treatment plans were recalculated or reoptimized on the daily imaging data acquired during MRgRT, and accumulated over all treatment fractions. Accumulated dose‐volume histogram (DVH) parameters of the gross tumor volume (GTV), lung, heart, and organs‐at‐risk (OARs) within 2 cm of the planning target volume (PTV) were extracted for each scenario and compared in Wilcoxon signed‐rank tests between S1 & S2, and S1 & S3. Results: The accumulated GTV D98% was above the prescribed dose for all patients and scenarios. Significant reductions (p < 0.05) of the mean ipsilateral lung dose (S2: –8%; S3: –23%) and mean heart dose (S2: –79%; S3: –83%) were observed for both proton scenarios compared to S1. The bronchial tree D0.1cc was significantly lower for S3 (S1: 48.1 Gy; S3: 39.2 Gy; p = 0.005), but not significantly different for S2 (S2: 45.0 Gy; p = 0.094), compared to S1. The D0.1cc for S2 and S3 compared to S1 was significantly (p < 0.05) smaller for OARs within 1–2 cm of the PTV (S1: 30.2 Gy; S2: 24.6 Gy; S3: 23.1 Gy), but not significantly different for OARs within 1 cm of the PTV. Conclusions: A significant dose sparing potential of non‐adaptive and online adaptive proton therapy compared to MRgRT for OARs in close, but not direct proximity of central lung tumors was identified. The near‐maximum dose to the bronchial tree was not significantly different for MRgRT and non‐adaptive IMPT. Online adaptive IMPT achieved significantly lower doses to the bronchial tree compared to MRgRT. [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]

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. 23 lung tumor patients received MR-guided hypofractionated stereotactic body radiation therapy at a 0.35 T 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 20 Gy (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). 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. MRI-based functional parameters extracted from 2D-cine MRI-scans were found to be predictive of RP development in lung tumor patients. [Display omitted] • Non-contrast enhanced functional imaging during lung SBRT at low-field MR-Linac. • Non-uniform Fourier decomposition to extract ventilation and perfusion over treatment. • Ventilation changes in high-dose region predictive of radiation-induced pneumonitis. • Higher predictive performance of functional than dosimetric parameters. [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]