Spatial profiling of in vivo diffusion-weighted MRI parameters in the healthy human kidney.

OBJECTIVE: Diffusion-weighted MRI is a technique that can infer microstructural and microcirculatory features from biological tissue, with particular application to renal tissue. There is extensive literature on diffusion tensor imaging (DTI) of anisotropy in the renal medulla, intravoxel incoherent motion (IVIM) measurements separating microstructural from microcirculation effects, and combinations of the two. However, interpretation of these features and adaptation of more specific models remains an ongoing challenge. One input to this process is a whole organ distillation of corticomedullary contrast of diffusion metrics, as has been explored for other renal biomarkers. MATERIALS AND METHODS: In this work, we probe the spatial dependence of diffusion MRI metrics with concentrically layered segmentation in 11 healthy kidneys at 3 T. The metrics include those from DTI, IVIM, a combined approach titled "REnal Flow and Microstructure AnisotroPy (REFMAP)", and a multiply encoded model titled "FC-IVIM" providing estimates of fluid velocity and branching length. RESULTS: Fractional anisotropy decreased from the inner kidney to the outer kidney with the strongest layer correlation in both parenchyma (including cortex and medulla) and medulla with Spearman correlation coefficients and p-values (r, p) of (0.42, <0.001) and (0.37, <0.001), respectively. Also, dynamic parameters derived from the three models significantly decreased with a high correlation from the inner to the outer parenchyma or medulla with (r, p) ranges of (0.46-0.55, <0.001). CONCLUSIONS: These spatial trends might find implications for indirect assessments of kidney physiology and microstructure using diffusion MRI.

Incidental coronary artery calcification on non-gated CT thorax correlates with risk of cardiovascular events and death.

OBJECTIVES: To assess coronary artery calcification (CAC) on non-contrast non-ECG-gated CT thorax (NC-NECG-CTT) and to evaluate its correlation with short-term risk of cardiovascular disease (CVD) events and death. METHODS: Single-institution retrospective study including all patients 40-70 years old who underwent NC-NECG-CTT over a period of 6 months. Individuals with known CVD were excluded. The presence of CAC was assessed and quantified by the Agatston score (CACS). CAC severity was defined as mild (< 100), moderate (100-400), or severe (> 400). CVD events (including CVD death, myocardial infarction, revascularisation procedures, ischaemic stroke, acute peripheral atherosclerotic ischaemia), and all-cause mortality over a median of 3.5 years were recorded. Cox proportional-hazards regression modelling was performed including CACS, age, gender and CVD risk factors (smoking, hypertension, diabetes mellitus, dyslipidaemia, and family history of CVD). RESULTS: Of the total 717 eligible cases, 325 (45%) had CAC. In patients without CAC, there was only one CVD event, compared to 26 CVD events including 5 deaths in patients with CAC. The presence and severity of CAC correlated with CVD events (p < 0.001). A CACS > 100 was significantly associated with both CVD events, hazard ratio (HR) 5.74, 95% confidence interval: 2.19-15.02; p < 0.001, and all-cause mortality, HR 1.7, 95% CI: 1.08-2.66; p = 0.02. Ever-smokers with CAC had a significantly higher risk for all-cause mortality compared to never-smokers (p = 0.03), but smoking status was not an independent predictor for CVD events in any subgroup category of CAC severity. CONCLUSIONS: The presence and severity of CAC assessed on NC-NECG-CTT correlates with short-term cardiovascular events and death. KEY POINTS: • Patients aged 40-70 years old without known CVD but with CAC on NC-NECG-CTT have a higher risk of CVD events compared to those without CAC. • CAC (Agatston) score above 100 confers a 5.7-fold increase in the risk of short-term CVD events in these patients. • The presence and severity of CAC on NC-NECG-CTT may have prognostic and therapeutic implications.

Echo time dependence of biexponential and triexponential intravoxel incoherent motion parameters in the liver.

PURPOSE: Intravoxel incoherent motion (IVIM) studies are performed with different acquisition protocols. Comparing them requires knowledge of echo time (TE) dependencies. The TE-dependence of the biexponential perfusion fraction f is well-documented, unlike that of its triexponential counterparts f1 and f2 and the biexponential and triexponential pseudodiffusion coefficients D* , D1∗ , and D2∗ . The purpose was to investigate the TE-dependence of these parameters and to check whether the triexponential pseudodiffusion compartments are associated with arterial and venous blood. METHODS: Fifteen healthy volunteers (19-58 y; mean: 24.7 y) underwent diffusion-weighted imaging of the abdomen with 24 b-values (0.2-800 s/mm2 ) at TEs of 45, 60, 75, and 90 ms. Regions of interest (ROIs) were manually drawn in the liver. One set of bi- and triexponential IVIM parameters per volunteer and TE was determined. The TE-dependence was assessed with the Kruskal-Wallis test. RESULTS: TE-dependence was observed for f (P < .001), f1 (P = .001), and f2 (P < .001). Their median values at the four measured TEs were: f: 0.198/0.240/0.274/0.359, f1 : 0.113/0.139/0.146/0.205, f2 : 0.115/0.155/0.182/0.194. D, D* , D1∗ , and D2∗ showed no significant TE-dependence. Their values were: diffusion coefficient D (10-4 mm2 /s): 9.45/9.63/9.75/9.41, biexponential D* (10-2 mm2 /s): 5.26/5.52/6.13/5.82, triexponential D1∗ (10-2 mm2 /s): 1.73/2.91/2.25/2.51, triexponential D2∗ (mm2 /s): 0.478/1.385/0.616/0.846. CONCLUSION: f1 and f2 show similar TE-dependence as f, ie, increase with rising TE; an effect that must be accounted for when comparing different studies. The diffusion and pseudodiffusion coefficients might be compared without TE correction. Because of the similar TE-dependence of f1 and f2 , the triexponential pseudodiffusion compartments are most probably not associated to venous and arterial blood.

The future of MRI in radiation therapy: Challenges and opportunities for the MR community.

Radiation therapy is a major component of cancer treatment pathways worldwide. The main aim of this treatment is to achieve tumor control through the delivery of ionizing radiation while preserving healthy tissues for minimal radiation toxicity. Because radiation therapy relies on accurate localization of the target and surrounding tissues, imaging plays a crucial role throughout the treatment chain. In the treatment planning phase, radiological images are essential for defining target volumes and organs-at-risk, as well as providing elemental composition (e.g., electron density) information for radiation dose calculations. At treatment, onboard imaging informs patient setup and could be used to guide radiation dose placement for sites affected by motion. Imaging is also an important tool for treatment response assessment and treatment plan adaptation. MRI, with its excellent soft tissue contrast and capacity to probe functional tissue properties, holds great untapped potential for transforming treatment paradigms in radiation therapy. The MR in Radiation Therapy ISMRM Study Group was established to provide a forum within the MR community to discuss the unmet needs and fuel opportunities for further advancement of MRI for radiation therapy applications. During the summer of 2021, the study group organized its first virtual workshop, attended by a diverse international group of clinicians, scientists, and clinical physicists, to explore our predictions for the future of MRI in radiation therapy for the next 25 years. This article reviews the main findings from the event and considers the opportunities and challenges of reaching our vision for the future in this expanding field.

An optimized b-value distribution for triexponential intravoxel incoherent motion (IVIM) in the liver.

PURPOSE: To find an optimized b-value distribution for reproducible triexponential intravoxel incoherent motion (IVIM) exams in the liver. METHODS: A numeric optimization of b-value distributions was performed using the triexponential IVIM equation and 27 different IVIM parameter sets. Starting with an initially optimized distribution of 6 b-values, the number of b-values was increased stepwise. Each new b-value was chosen from a set of 64 predefined b-values based on the computed summed relative mean error of the fitted triexponential IVIM parameters. This process was repeated for up to 100 b-values. In simulations and in vivo measurements, optimized b-value distributions were compared to 4 representative distributions found in literature. RESULTS: The first 16 optimized b-values were 0, 0.3, 0.3, 70, 200, 800, 70, 1, 3.5, 5, 70, 1.2, 6, 45, 1.5, and 60 in units of s/mm2 . Low b-values were much more frequent than high b-values. The optimized b-value distribution resulted in a higher fit stability compared to distributions used in literature in both, simulation and in vivo measurements. Using more than 6 b-values, ideally 16 or more, increased the fit stability considerably. CONCLUSION: Using optimized b-values, the fit uncertainty in triexponential IVIM can be largely reduced. Ideally, 16 or more b-values should be acquired.

Optimal acquisition scheme for flow-compensated intravoxel incoherent motion diffusion-weighted imaging in the abdomen: An accurate and precise clinically feasible protocol.

PURPOSE: Flow-compensated (FC) diffusion-weighted MRI (DWI) for intravoxel-incoherent motion (IVIM) modeling allows for a more detailed description of tissue microvasculature than conventional IVIM. The long acquisition time of current FC-IVIM protocols, however, has prohibited clinical application. Therefore, we developed an optimized abdominal FC-IVIM acquisition with a clinically feasible scan time. METHODS: Precision and accuracy of the FC-IVIM parameters were assessed by fitting the FC-IVIM model to signal decay curves, simulated for different acquisition schemes. Diffusion-weighted acquisitions were added subsequently to the protocol, where we chose the combination of b-value, diffusion time and gradient profile (FC or bipolar) that resulted in the largest improvement to its accuracy and precision. The resulting two optimized FC-IVIM protocols with 25 and 50 acquisitions (FC-IVIMopt25 and FC-IVIMopt50 ), together with a complementary acquisition consisting of 50 diffusion-weighting (FC-IVIMcomp ), were acquired in repeated abdominal free-breathing FC-IVIM imaging of seven healthy volunteers. Intersession and intrasession within-subject coefficient of variation of the FC-IVIM parameters were compared for the liver, spleen, and kidneys. RESULTS: Simulations showed that the performance of FC-IVIM improved in tissue with larger perfusion fraction and signal-to-noise ratio. The scan time of the FC-IVIMopt25 and FC-IVIMopt50 protocols were 8 and 16 min. The best in vivo performance was seen in FC-IVIMopt50 . The intersession within-subject coefficients of variation of FC-IVIMopt50 were 11.6%, 16.3%, 65.5%, and 36.0% for FC-IVIM model parameters diffusivity, perfusion fraction, characteristic time and blood flow velocity, respectively. CONCLUSIONS: We have optimized the FC-IVIM protocol, allowing for clinically feasible scan times (8-16 min).

Twice-refocused stimulated echo diffusion imaging: Measuring diffusion time dependence at constant T1 weighting.

PURPOSE: Diffusion times longer than 50 ms are typically probed with stimulated-echo sequences. Varying the diffusion time in stimulated-echo sequences affects the T1 weighting of subcompartments, complicating the analysis of diffusion time dependence. Although inversion recovery preparation could be used to change the T1 weighting, it cannot ensure equal T1 weighting at arbitrary mixing times. In this article, a sequence that ensures constant T1 weighting over a wide range of diffusion times is presented. METHODS: The proposed sequence features 2 independent longitudinal storage periods: TM1 and TM2 . Diffusion encoding is performed during TM1 , effectively coupling the diffusion time and TM1 . Equal T1 weighting at arbitrary diffusion times is realized by keeping the total mixing time TM1 + TM2 constant. The sequence was compared with conventional stimulated-echo measurements of diffusion in a 2-compartment phantom consisting of distilled water and paraffinum perliquidum. Additionally, in vivo DTI of the brain was carried out for 8 healthy volunteers with diffusion times ranging from 50 to 500 ms. RESULTS: Diffusion time dependence of the axial and radial diffusivity was detected in the brain. Both sequences resulted in almost identical diffusivities in white matter. In regions containing partial volumes of gray and white matter, a dependency on T1 weighting was observed. CONCLUSION: In accordance with previous studies, little variance of T1 values appeared to be present in healthy white matter. However, this is likely different in diseased tissue. Here, the proposed sequence can be effective in differentiating between diffusion time dependence and T1 weighting effects.

On the Field Strength Dependence of Bi- and Triexponential Intravoxel Incoherent Motion (IVIM) Parameters in the Liver.

BACKGROUND: Studies on intravoxel incoherent motion (IVIM) imaging are carried out with different acquisition protocols. PURPOSE: To investigate the dependence of IVIM parameters on the B0 field strength when using a bi- or triexponential model. STUDY TYPE: Prospective. STUDY POPULATION: 20 healthy volunteers (age: 19-28 years). FIELD STRENGTH/SEQUENCE: Volunteers were examined at two field strengths (1.5 and 3T). Diffusion-weighted images of the abdomen were acquired at 24 b-values ranging from 0.2 to 500 s/mm2 . ASSESSMENT: ROIs were manually drawn in the liver. Data were fitted with a bi- and a triexponential IVIM model. The resulting parameters were compared between both field strengths. STATISTICAL TESTS: One-way analysis of variance (ANOVA) and Kruskal-Wallis test were used to test the obtained IVIM parameters for a significant field strength dependency. RESULTS: At b-values below 6 s/mm2 , the triexponential model provided better agreement with the data than the biexponential model. The average tissue diffusivity was D = 1.22/1.00 µm2 /msec at 1.5/3T. The average pseudodiffusion coefficients for the biexponential model were D* = 308/260 µm2 /msec at 1.5/3T; and for the triexponential model D1* = 81.3/65.9 µm2 /msec, D2* = 2453/2333 µm2 /msec at 1.5/3T. The average perfusion fractions for the biexponential model were f = 0.286/0.303 at 1.5/3T; and for the triexponential model f1 = 0.161/0.174 and f2 = 0.152/0.159 at 1.5/3T. A significant B0 dependence was only found for the biexponential pseudodiffusion coefficient (ANOVA/KW P = 0.037/0.0453) and tissue diffusivity (ANOVA/KW: P < 0.001). DATA CONCLUSION: Our experimental results suggest that triexponential pseudodiffusion coefficients and perfusion fractions obtained at different field strengths could be compared across different studies using different B0 . However, it is recommended to take the field strength into account when comparing tissue diffusivities or using the biexponential IVIM model. Considering published values for oxygenation-dependent transversal relaxation times of blood, it is unlikely that the two blood compartments of the triexponential model represent venous and arterial blood. LEVEL OF EVIDENCE: 1 Technical Efficacy Stage: 2 J. Magn. Reson. Imaging 2019;50:1883-1892.

Characterization of the diffusion coefficient of blood.

PURPOSE: To characterize the diffusion coefficient of human blood for accurate results in intravoxel incoherent motion imaging. METHODS: Diffusion-weighted MRI of blood samples from 10 healthy volunteers was acquired with a single-shot echo-planar-imaging sequence at body temperature. Effects of gradient profile (monopolar or flow-compensated), diffusion time (40-100 ms), and echo time (60-200 ms) were investigated. RESULTS: Although measured apparent diffusion coefficients of blood were larger for flow-compensated than for monopolar gradients, no dependence of the apparent diffusion coefficient on the diffusion time was found. Large differences between individual samples were observed, with results ranging from 1.26 to 1.66 µm2 /ms for flow-compensated and 0.94 to 1.52 µm2 /ms for monopolar gradients. Statistical analysis indicates correlations of the flow-compensated apparent diffusion coefficient with hematocrit (P = 0.007) and hemoglobin (P = 0.017), but not with mean corpuscular volume (P = 0.64). Results of Monte-Carlo simulations support the experimental observations. CONCLUSIONS: Measured blood apparent diffusion coefficient values depend on hematocrit/hemoglobin concentration and applied gradient profile due to non-Gaussian diffusion. Because in vivo measurement is delicate, an estimation based on blood count results could be an alternative. For intravoxel incoherent motion modeling, the use of a blood self-diffusion constant Db = 1.54 ± 0.12 µm2 /ms for flow-compensated and Db = 1.30 ± 0.18 µm2 /ms for monopolar encoding is suggested. Magn Reson Med 79:2752-2758, 2018. © 2017 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Quantitative susceptibility mapping and 23 Na imaging-based in vitro characterization of blood clotting kinetics.

Blood clotting is a fundamental biochemical process in post-hemorrhagic hemostasis. Although the varying appearance of coagulating blood in T1 - and T2 -weighted images is widely used to qualitatively determine bleeding age, the technique permits only a rough discrimination of coagulation stages, and it remains difficult to distinguish acute and chronic hemorrhagic stages because of low T1 - and T2 -weighted signal intensities in both instances. To investigate new biomedical parameters for magnetic resonance imaging-based characterization of blood clotting kinetics, sodium imaging and quantitative susceptibility mapping (QSM) were compared with conventional T1 - and T2 -weighted imaging, as well as with biochemical hemolysis parameters. For this purpose, a blood-filled spherical agar phantom was investigated daily for 14 days, as well as after 24 days at 7 T after initial preparation with fresh blood. T1 - and T2 -weighted sequences, a three-dimensional (3D) gradient echo sequence and a density-adapted 3D radial projection reconstruction pulse sequence for 23 Na imaging were applied. For hemolysis estimations, free hemoglobin and free potassium concentrations were measured photometrically and with the direct ion-selective electrode method, respectively, in separate heparinized whole-blood samples along the same timeline. Initial mean susceptibility was low (0.154 ± 0.020 ppm) and increased steadily during the course of coagulation to reach up to 0.570 ± 0.165 ppm. The highest total sodium (NaT) values (1.02 ± 0.06 arbitrary units) in the clot were observed initially, dropped to 0.69 ± 0.13 arbitrary units after one day and increased again to initial values. Compartmentalized sodium (NaS) showed a similar signal evolution, and the NaS/NaT ratio steadily increased over clot evolution. QSM depicts clot evolution in vitro as a process associated with hemoglobin accumulation and transformation, and enables the differentiation of the acute and chronic coagulation stages. Sodium imaging visualizes clotting independent of susceptibility and seems to correspond to clot integrity. A combination of QSM and sodium imaging may enhance the characterization of hemorrhage.

Eddy current compensated double diffusion encoded (DDE) MRI.

PURPOSE: Eddy currents might lead to image distortions in diffusion-weighted echo planar imaging. A method is proposed to reduce their effects on double diffusion encoding (DDE) MRI experiments and the thereby derived microscopic fractional anisotropy (µFA). METHODS: The twice-refocused spin echo scheme was adapted for DDE measurements. To assess the effect of individual diffusion encodings on the image distortions, measurements of a grid of plastic rods in water were performed. The effect of eddy current compensation on µFA measurements was evaluated in the brains of six healthy volunteers. RESULTS: The use of an eddy current compensation reduced the signal variation. As expected, the distortions caused by the second encoding were larger than those of the first encoding, entailing a stronger need to compensate for them. For an optimal result, however, both encodings had to be compensated. The artifact reduction strongly improved the measurement of the µFA in ventricles and gray matter by reducing the overestimation. An effect of the compensation on absolute µFA values in white matter was not observed. CONCLUSION: It is advisable to compensate both encodings in DDE measurements for eddy currents. Magn Reson Med 77:328-335, 2017. © 2015 Wiley Periodicals, Inc.

4D respiratory motion-compensated image reconstruction of free-breathing radial MR data with very high undersampling.

PURPOSE: To develop four-dimensional (4D) respiratory time-resolved MRI based on free-breathing acquisition of radial MR data with very high undersampling. METHODS: We propose the 4D joint motion-compensated high-dimensional total variation (4D joint MoCo-HDTV) algorithm, which alternates between motion-compensated image reconstruction and artifact-robust motion estimation at multiple resolution levels. The algorithm is applied to radial MR data of the thorax and upper abdomen of 12 free-breathing subjects with acquisition times between 37 and 41 s and undersampling factors of 16.8. Resulting images are compared with compressed sensing-based 4D motion-adaptive spatio-temporal regularization (MASTeR) and 4D high-dimensional total variation (HDTV) reconstructions. RESULTS: For all subjects, 4D joint MoCo-HDTV achieves higher similarity in terms of normalized mutual information and cross-correlation than 4D MASTeR and 4D HDTV when compared with reference 4D gated gridding reconstructions with 8.4 ± 1.1 times longer acquisition times. In a qualitative assessment of artifact level and image sharpness by two radiologists, 4D joint MoCo-HDTV reveals higher scores (P < 0.05) than 4D HDTV and 4D MASTeR at the same undersampling factor and the reference 4D gated gridding reconstructions, respectively. CONCLUSIONS: 4D joint MoCo-HDTV enables time-resolved image reconstruction of free-breathing radial MR data with undersampling factors of 16.8 while achieving low-streak artifact levels and high image sharpness. Magn Reson Med 77:1170-1183, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

On contrast mechanisms in p-space imaging.

PURPOSE: The effects of microstructure on susceptibility mapping have recently received increasing attention. In this work, the capability of p-space imaging to resolve subvoxel structure and susceptibility is assessed. METHODS: In a simulation study, the p-space contrast of axon bundles comprised of hollow cylinders is investigated. Various axon and susceptibility distributions are considered and compared with the corresponding case of a voxel with homogeneous substructure of the order of the voxel size. RESULTS: MR signal behavior for p-space imaging of a voxel containing axon bundles and a voxel containing homogeneous substructure are nearly identical. CONCLUSION: p-Space imaging resolves subvoxel structure of the order of the voxel size. Due to dephasing effects, strong susceptibility variations alter p-space contrast. However, p-space contrast is not directly linked to the susceptibility of the axon compartment. Magn Reson Med 75:2526-2533, 2016. © 2015 Wiley Periodicals, Inc.

Regional Lung Ventilation Analysis Using Temporally Resolved Magnetic Resonance Imaging.

OBJECTIVE: We propose a computer-aided method for regional ventilation analysis and observation of lung diseases in temporally resolved magnetic resonance imaging (4D MRI). METHODS: A shape model-based segmentation and registration workflow was used to create an atlas-derived reference system in which regional tissue motion can be quantified and multimodal image data can be compared regionally. Model-based temporal registration of the lung surfaces in 4D MRI data was compared with the registration of 4D computed tomography (CT) images. A ventilation analysis was performed on 4D MR images of patients with lung fibrosis; 4D MR ventilation maps were compared with corresponding diagnostic 3D CT images of the patients and 4D CT maps of subjects without impaired lung function (serving as reference). RESULTS: Comparison between the computed patient-specific 4D MR regional ventilation maps and diagnostic CT images shows good correlation in conspicuous regions. Comparison to 4D CT-derived ventilation maps supports the plausibility of the 4D MR maps. Dynamic MRI-based flow-volume loops and spirograms further visualize the free-breathing behavior. CONCLUSIONS: The proposed methods allow for 4D MR-based regional analysis of tissue dynamics and ventilation in spontaneous breathing and comparison of patient data. The proposed atlas-based reference coordinate system provides an automated manner of annotating and comparing multimodal lung image data.

Flow-compensated intravoxel incoherent motion diffusion imaging.

PURPOSE: The pseudo-diffusion coefficient D* in intravoxel incoherent motion (IVIM) imaging was found difficult to seize. Flow-compensated diffusion gradients were used to test the validity of the commonly assumed biexponential limit and to determine not only D*, but also characteristic timescale τ and velocity v of the incoherent motion. THEORY AND METHODS: Bipolar and flow-compensated diffusion gradients were inserted into a flow-compensated single-shot EPI sequence. Images were obtained from a pipe-shaped flow phantom and from healthy volunteers. To calculate the IVIM signal outside the biexponential limit, a formalism based on normalized phase distributions was developed. RESULTS: The flow-compensated diffusion gradients caused less signal attenuation than the bipolar ones. A signal dependence on the duration of the flow-compensated gradients was found at low b-values in the volunteer datasets. The characteristic IVIM parameters were estimated to be v = 4.60 ± 0.34 mm/s and τ = 144 ± 10 ms for liver and v = 3.91 ± 0.54 mm/s and τ = 224 ± 47 ms for pancreas. CONCLUSION: Our results strongly indicate that the biexponential limit does not adequately model the diffusion signal in liver and pancreas. By using both bipolar and flow-compensated diffusion gradients of different duration, the characteristic timescale and velocity of the incoherent motion can be determined.

Real-time motion-including dose estimation of simulated multi-leaf collimator-tracked magnetic resonance-guided radiotherapy.

BACKGROUND: Real-time dose estimation is a key-prerequisite to enable online intra-fraction treatment adaptation in magnetic resonance (MR)-guided radiotherapy (MRgRT). It is an essential component for the assessment of the dosimetric benefits and risks of online adaptive treatments, such as multi-leaf collimator (MLC)-tracking. PURPOSE: We present a proof-of-concept for a software workflow for real-time dose estimation of MR-guided adaptive radiotherapy based on real-time data-streams of the linac delivery parameters and target positions. METHODS: A software workflow, combining our in-house motion management software DynaTrack, a real-time dose calculation engine that connects to a research version of the treatment planning software (TPS) Monaco (v.6.09.00, Elekta AB, Stockholm, Sweden) was developed and evaluated. MR-guided treatment delivery on the Elekta Unity MR-linac was simulated with and without MLC-tracking for three prostate patients, previously treated on the Elekta Unity MR-linac (36.25 Gy/five fractions). Three motion scenarios were used: no motion, regular motion, and erratic prostate motion. Accumulated monitor units (MUs), centre of mass target position and MLC-leaf positions, were forwarded from DynaTrack at a rate of 25 Hz to a Monte Carlo (MC) based dose calculation engine which utilises the research GPUMCD-library (Elekta AB, Stockholm, Sweden). A rigid isocentre shift derived from the selected motion scenarios was applied to a bulk density-assigned session MR-image. The respective electron density used for treatment planning was accessed through the research Monaco TPS. The software workflow including the online dose reconstruction was validated against offline dose reconstructions. Our investigation showed that MC-based real-time dose calculations that account for all linac states (including MUs, MLC positions and target position) were infeasible, hence states were randomly sampled and used for calculation as follows; Once a new linac state was received, a dose calculation with 106 photons was started. Linac states that arrived during the time of the ongoing calculation were put into a queue. After completion of the ongoing calculation, one new linac state was randomly picked from the queue and assigned the MU accumulated from the previous state until the last sample in the queue. The queue was emptied, and the process repeated throughout treatment simulation. RESULTS: On average 27% (23%-30%) of received samples were used in the real-time calculation, corresponding to a calculation time for one linac state of 148 ms. Median gamma pass rate (2%/3 mm local) was 100.0% (99.9%-100%) within the PTV volume and 99.1% (90.1%-99.4.0%) with a 15% dose cut off. Differences in PTVDmean , CTVDmean , RectumD2% , and BladderD2% (offline-online, % of prescribed dose) were below 0.64%. Beam-by-beam comparisons showed deviations below 0.07 Gy. Repeated simulations resulted in standard deviations below 0.31% and 0.12 Gy for the investigated volume and dose criteria respectively. CONCLUSIONS: Real-time dose estimation was successfully performed using the developed software workflow for different prostate motion traces with and without MLC-tracking. Negligible dosimetric differences were seen when comparing online and offline reconstructed dose, enabling online intra-fraction treatment decisions based on estimates of the delivered dose.

Respiratory motion modelling for MR-guided lung cancer radiotherapy: model development and geometric accuracy evaluation.

Objective.Respiratory motion of lung tumours and adjacent structures is challenging for radiotherapy. Online MR-imaging cannot currently provide real-time volumetric information of the moving patient anatomy, therefore limiting precise dose delivery, delivered dose reconstruction, and downstream adaptation methods.Approach.We tailor a respiratory motion modelling framework towards an MR-Linac workflow to estimate the time-resolved 4D motion from real-time data. We develop a multi-slice acquisition scheme which acquires thick, overlapping 2D motion-slices in different locations and orientations, interleaved with 2D surrogate-slices from a fixed location. The framework fits a motion model directly to the input data without the need for sorting or binning to account for inter- and intra-cycle variation of the breathing motion. The framework alternates between model fitting and motion-compensated super-resolution image reconstruction to recover a high-quality motion-free image and a motion model. The fitted model can then estimate the 4D motion from 2D surrogate-slices. The framework is applied to four simulated anthropomorphic datasets and evaluated against known ground truth anatomy and motion. Clinical applicability is demonstrated by applying our framework to eight datasets acquired on an MR-Linac from four lung cancer patients.Main results.The framework accurately reconstructs high-quality motion-compensated 3D images with 2 mm3isotropic voxels. For the simulated case with the largest target motion, the motion model achieved a mean deformation field error of 1.13 mm. For the patient cases residual error registrations estimate the model error to be 1.07 mm (1.64 mm), 0.91 mm (1.32 mm), and 0.88 mm (1.33 mm) in superior-inferior, anterior-posterior, and left-right directions respectively for the building (application) data.Significance.The motion modelling framework estimates the patient motion with high accuracy and accurately reconstructs the anatomy. The image acquisition scheme can be flexibly integrated into an MR-Linac workflow whilst maintaining the capability of online motion-management strategies based on cine imaging such as target tracking and/or gating.

Accelerating 4D image reconstruction for magnetic resonance-guided radiotherapy.

Background and purpose: Physiological motion impacts the dose delivered to tumours and vital organs in external beam radiotherapy and particularly in particle therapy. The excellent soft-tissue demarcation of 4D magnetic resonance imaging (4D-MRI) could inform on intra-fractional motion, but long image reconstruction times hinder its use in online treatment adaptation. Here we employ techniques from high-performance computing to reduce 4D-MRI reconstruction times below two minutes to facilitate their use in MR-guided radiotherapy. Material and methods: Four patients with pancreatic adenocarcinoma were scanned with a radial stack-of-stars gradient echo sequence on a 1.5T MR-Linac. Fast parallelised open-source implementations of the extra-dimensional golden-angle radial sparse parallel algorithm were developed for central processing unit (CPU) and graphics processing unit (GPU) architectures. We assessed the impact of architecture, oversampling and respiratory binning strategy on 4D-MRI reconstruction time and compared images using the structural similarity (SSIM) index against a MATLAB reference implementation. Scaling and bottlenecks for the different architectures were studied using multi-GPU systems. Results: All reconstructed 4D-MRI were identical to the reference implementation (SSIM > 0.99). Images reconstructed with overlapping respiratory bins were sharper at the cost of longer reconstruction times. The CPU  + GPU implementation was over 17 times faster than the reference implementation, reconstructing images in 60 ± 1 s and hyper-scaled using multiple GPUs. Conclusion: Respiratory-resolved 4D-MRI reconstruction times can be reduced using high-performance computing methods for online workflows in MR-guided radiotherapy with potential applications in particle therapy.

Joint radial trajectory correction for accelerated T2 * mapping on an MR-Linac.

BACKGROUND: T2 * mapping can characterize tumor hypoxia, which may be associated with resistance to therapy. Acquiring T2 * maps during MR-guided radiotherapy could inform treatment adaptation by, for example, escalating the dose to resistant sub-volumes. PURPOSE: The purpose of this work is to demonstrate the feasibility of the accelerated T2 * mapping technique using model-based image reconstruction with integrated trajectory auto-correction (TrACR) for MR-guided radiotherapy on an MR-Linear accelerator (MR-Linac). MATERIALS AND METHODS: The proposed method was validated in a numerical phantom, where two T2 * mapping approaches (sequential and joint) were compared for different noise levels (0,0.1,0.5,1) and gradient delays ([1, -1] and [1, -2] in units of dwell time for x- and y-axis, respectively). Fully sampled k-space was retrospectively undersampled using two different undersampling patterns. Root mean square errors (RMSEs) were calculated between reconstructed T2 * maps and ground truth. In vivo data was acquired twice weekly in one prostate and one head and neck cancer patient undergoing treatment on a 1.5 T MR-Linac. Data were retrospectively undersampled and T2 * maps reconstructed, with and without trajectory corrections were compared. RESULTS: Numerical simulations demonstrated that, for all noise levels, T2 * maps reconstructed with a joint approach demonstrated less error compared to an uncorrected and sequential approach. For a noise level of 0.1, uniform undersampling and gradient delay [1, -1] (in units of dwell time for x- and y-axis, respectively), RMSEs for sequential and joint approaches were 13.01 and 9.32 ms, respectively, which reduced to 10.92 and 5.89 ms for a gradient delay of [1, 2]. Similarly, for alternate undersampling and gradient delay [1, -1], RMSEs for sequential and joint approaches were 9.80 and 8.90 ms, respectively, which reduced to 9.10 and 5.40 ms for gradient delay [1, 2]. For in vivo data, T2 * maps reconstructed with our proposed approach resulted in less artifacts and improved visual appearance compared to the uncorrected approach. For both prostate and head and neck cancer patients, T2 * maps reconstructed from different treatment fractions showed changes within the planning target volume (PTV). CONCLUSION: Using the proposed approach, a retrospective data-driven gradient delay correction can be performed, which is particularly relevant for hybrid devices, where full information on the machine configuration is not available for image reconstruction. T2 * maps were acquired in under 5 min and can be integrated into MR-guided radiotherapy treatment workflows, which minimizes patient burden and leaves time for additional imaging for online adaptive radiotherapy on an MR-Linac.