3D sodium (23 Na) magnetic resonance fingerprinting for time-efficient relaxometric mapping.

PURPOSE: To develop a framework for 3D sodium (23 Na) MR fingerprinting (MRF), based on irreducible spherical tensor operators with tailored flip angle (FA) pattern and time-efficient data acquisition for simultaneous quantification of T1 , T2l∗ , T2s∗ , and T2∗ in addition to ΔB0 . METHODS: 23 Na-MRF was implemented in a 3D sequence and irreducible spherical tensor operators were exploited in the simulations. Furthermore, the Cramér Rao lower bound was used to optimize the flip angle pattern. A combination of single and double echo readouts was implemented to increase the readout efficiency. A study was conducted to compare results in a multicompartment phantom acquired with MRF and reference methods. Finally, the relaxation times in the human brain were measured in four healthy volunteers. RESULTS: Phantom experiments revealed a mean difference of 1.0% between relaxation times acquired with MRF and results determined with the reference methods. Simultaneous quantification of the longitudinal and transverse relaxation times in the human brain was possible within 32 min using 3D 23 Na-MRF with a nominal resolution of (5 mm)3 . In vivo measurements in four volunteers yielded average relaxation times of: T1,brain = (35.0 ± 3.2) ms, T2l,brain∗ = (29.3 ± 3.8) ms and T2s,brain∗ = (5.5 ± 1.3) ms in brain tissue, whereas T1,CSF = (61.9 ± 2.8) ms and T2,CSF∗ = (46.3 ± 4.5) ms was found in cerebrospinal fluid. CONCLUSION: The feasibility of in vivo 3D relaxometric sodium mapping within roughly ½ h is demonstrated using MRF in the human brain, moving sodium relaxometric mapping toward clinically relevant measurement times.

Sodium relaxometry using 23 Na MR fingerprinting: A proof of concept.

PURPOSE: To evaluate the feasibility of 23 Na MR fingerprinting (MRF) for simultaneous quantification of T1 , T2l∗ , T2s∗ , T2∗ in addition to ΔB0 . METHODS: A framework for sodium relaxometry using MRF at 7T was developed, allowing simultaneous measurement of relaxation times and inhomogeneities in the static field. The technique distinguishes between bi- and monoexponential transverse relaxation and was validated in simulations with respect to the ground truth. In phantom measurements, a resolution of 2 × 2 × 12 mm3 was achieved within 1 h acquisition time, and the resulting parameter maps were compared to results from reference methods. Relaxation times in five healthy volunteers were measured with a resolution of 4 × 4 × 12 mm3 . RESULTS: Phantom experiments revealed an agreement between the relaxation times obtained via 23 Na-MRF and the reference methods. In white matter, a longitudinal relaxation constant of T1 = 38.9 ± 4.8 ms was found, while values of T2l∗ = 29.2 ± 4.9 ms and T2s∗ = 4.7 ± 1.2 ms were found for the long and short component of the transverse relaxation. In cerebrospinal fluid, T1 was 67.7 ± 6.3 ms and T2∗ = 41.5 ± 3.4 ms. CONCLUSION: This work demonstrates the feasibility of 23 Na-MRF for relaxometry in sodium MRI in both phantom and in vivo studies. Simultaneous quantification of T1 , T2l∗ , T2s∗ , T2∗ and ΔB0 was possible within a 1 h measurement time.

Mapping the human brainstem: Brain nuclei and fiber tracts at 3 T and 7 T.

Structural high-resolution imaging of the brainstem can be of high importance in clinical practice. However, ultra-high field magnetic resonance imaging (MRI) is still restricted in use due to limited availability. Therefore, quantitative MRI techniques (quantitative susceptibility mapping [QSM], relaxation measurements [ R2* , R1 ], diffusion tensor imaging [DTI]) and T2 - and proton density (PD)-weighted imaging in the human brainstem at 3 T and 7 T are compared. Five healthy volunteers (mean age: 21.5 ± 1.9 years) were measured at 3 T and 7 T using multi-echo gradient echo sequences for susceptibility mapping and R2* relaxometry, magnetization-prepared 2 rapid acquisition gradient echo sequences for R1 relaxometry, turbo-spin echo sequences for PD- and T2 -weighted imaging and readout-segmented echo planar sequences for DTI. Susceptibility maps were computed using Laplacian-based phase unwrapping, V-SHARP for background field removal and the streaking artifact reduction for QSM algorithm for dipole inversion. Contrast-to-noise ratios (CNRs) were determined at 3 T and 7 T in ten volumes of interest (VOIs). Data acquired at 7 T showed higher CNR. However, in four VOIs, lower CNR was observed for R2* at 7 T. QSM was shown to be the contrast with which the highest number of structures could be identified. The depiction of very fine tracts such as the medial longitudinal fasciculus throughout the brainstem was only possible in susceptibility maps acquired at 7 T. DTI effectively showed the main tracts (crus cerebri, transverse pontine fibers, corticospinal tract, middle and superior cerebellar peduncle, pontocerebellar tract, and pyramid) at both field strengths. Assessing the brainstem with quantitative MRI methods such as QSM, R2* , as well as PD- and T2 -weighted imaging with great detail, is also possible at 3 T, especially when using susceptibility mapping calculated from a gradient echo sequence with a wide range of echo times from 10.5 to 52.5 ms. However, tracing smallest structures strongly benefits from imaging at ultra-high field.

Quantification of patellofemoral cartilage deformation and contact area changes in response to static loading via high-resolution MRI with prospective motion correction.

BACKGROUND: Higher-resolution MRI of the patellofemoral cartilage under loading is hampered by subject motion since knee flexion is required during the scan. PURPOSE: To demonstrate robust quantification of cartilage compression and contact area changes in response to in situ loading by means of MRI with prospective motion correction and regularized image postprocessing. STUDY TYPE: Cohort study. SUBJECTS: Fifteen healthy male subjects. FIELD STRENGTH: 3 T. SEQUENCE: Spoiled 3D gradient-echo sequence augmented with prospective motion correction based on optical tracking. Measurements were performed with three different loads (0/200/400 N). ASSESSMENT: Bone and cartilage segmentation was performed manually and regularized with a deep-learning approach. Average patellar and femoral cartilage thickness and contact area were calculated for the three loading situations. Reproducibility was assessed via repeated measurements in one subject. STATISTICAL TESTS: Comparison of the three loading situations was performed by Wilcoxon signed-rank tests. RESULTS: Regularization using a deep convolutional neural network reduced the variance of the quantified relative load-induced changes of cartilage thickness and contact area compared to purely manual segmentation (average reduction of standard deviation by ∼50%) and repeated measurements performed on the same subject demonstrated high reproducibility of the method. For the three loading situations (0/200/400 N), the patellofemoral cartilage contact area as well as the mean patellar and femoral cartilage thickness were significantly different from each other (P < 0.05). While the patellofemoral cartilage contact area increased under loading (by 14.5/19.0% for loads of 200/400 N), patellar and femoral cartilage thickness exhibited a load-dependent thickness decrease (patella: -4.4/-7.4%, femur: -3.4/-7.1% for loads of 200/400 N). DATA CONCLUSION: MRI with prospective motion correction enables quantitative evaluation of patellofemoral cartilage deformation and contact area changes in response to in situ loading. Regularizing the manual segmentations using a neural network enables robust quantification of the load-induced changes. LEVEL OF EVIDENCE: 2 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2019;50:1561-1570.

Prospective motion correction in functional MRI.

Due to the intrinsic low sensitivity of BOLD-fMRI long scanning is required. Subject motion during fMRI scans reduces statistical significance of the activation maps and increases the prevalence of false activations. Motion correction is therefore an essential tool for a successful fMRI data analysis. Retrospective motion correction techniques are now commonplace and are incorporated into a wide range of fMRI analysis toolboxes. These techniques are advantageous due to robustness, sequence independence and have minimal impact on the fMRI study setup. Retrospective techniques however, do not provide an accurate intra-volume correction, nor can these techniques correct for the spin-history effects. The application of prospective motion correction in fMRI appears to be effective in reducing false positives and increasing sensitivity when compared to retrospective techniques, particularly in the cases of substantial motion. Especially advantageous in this regard is the combination of prospective motion correction with dynamic distortion correction. Nevertheless, none of the recent methods are able to recover activations in presence of motion that are comparable to no-motion conditions, which motivates further research in the area of adaptive dynamic imaging.

Comparative T2 and T1ρ mapping of patellofemoral cartilage under in situ mechanical loading with prospective motion correction.

PURPOSE: To demonstrate improved T2 and T1ρ mapping of patellofemoral cartilage with in situ loading by means of prospective motion correction and to assess load-induced changes in healthy subjects. MATERIALS AND METHODS: Established T2 and T1ρ mapping sequences were augmented with prospective motion correction based on optical tracking. Protocols were optimized for robust imaging of the patellofemoral cartilage at a field strength of 3T. Subjects were positioned in the scanner with knee flexion and in situ loading of the patellofemoral joint was performed with a pneumatic loading device. In a pilot study on a cohort of 10 healthy subjects, load-induced T2 and T1ρ changes were evaluated through measurements with axial loads of 0/20/40 kg. RESULTS: With prospective motion correction and additional lipid saturation, motion artifacts in patellofemoral cartilage magnetic resonance imaging (MRI) with in situ loading could be notably decreased, as demonstrated for T2 mapping. The acquired relaxation maps suggested a T2 /T1ρ decrease in superficial cartilage and a T2 /T1ρ increase in deep cartilage under loading. However, in the quantitative group evaluation of the lateral patellar facet, only T1ρ in superficial cartilage was significantly changed by loading (P ≤ 0.05), while no significant T2 differences for the three loading conditions were observed (P ≥ 0.3). CONCLUSION: Prospective motion correction enables T2 and T1ρ mapping of patellofemoral cartilage with in situ loading and a comparison of the two contrasts in terms of their response to mechanical loading. T1ρ is a more sensitive marker for load-induced patellar cartilage changes than T2 . LEVEL OF EVIDENCE: 3 Technical Efficacy: Stage 1 J. MAGN. RESON. IMAGING 2017;46:452-460.

Left atrial late gadolinium enhancement with water-fat separation: the importance of phase-encoding order.

PURPOSE: To compare two late gadolinium enhancement (LGE) cardiovascular magnetic resonance (CMR) methods: a Dixon LGE sequence with sequential phase-encoding order, reconstructed using water-fat separation, and standard fat-saturated LGE. MATERIALS AND METHODS: We implemented a dual-echo Dixon LGE method for reconstructing water-only images and compared it to fat-saturated LGE in 12 patients prior to their first pulmonary vein isolation (PVI) procedure. Images were analyzed for quality and fat-suppression. Regions of the left atrium were evaluated by a blinded observer (1 = prominent enhancement, 0 = mild or absent enhancement) on two sets of images (fat-saturated and water-only LGE) and agreement was assessed. RESULTS: Water-only LGE showed a trend toward better fat-suppression (P = 0.06), with a significantly more homogeneous blood pool signal and reduced inflow artifacts (both P < 0.01). Agreement between fat-saturated LGE and water-only methods was found in 84% of regions, significantly correlated by chi-squared test (P < 0.001). The kappa value was 0.52 (moderate). The average number of enhancing segments was higher for fat-saturated LGE than water-only LGE (4.2 ± 2.7 vs. 3.2 ± 2.9, P = 0.03). CONCLUSION: The two-point Dixon LGE technique reduces artifacts due to a centric k-space order. A similar enhancement pattern was observed irrespective of the LGE technique, with more enhancement detected by fat-saturated LGE.

Three-dimensional late gadolinium-enhanced MR imaging of the left atrium: a comparison of spiral versus Cartesian k-space trajectories.

PURPOSE: To investigate the feasibility of high-resolution late gadolinium enhancement (LGE) imaging using a three-dimensional (3D) stack of spirals k-space trajectory for the detection of left atrial (LA) ablation lesions. LGE imaging inherently suffers from low SNR, so that improvements in spatial resolution and imaging time are challenging. The spiral trajectory offers greater acquisition efficiency, and this is used for increased spatial resolution. MATERIALS AND METHODS: Nine healthy subjects and 10 pre/post pulmonary vein isolation patients underwent an MRI examination. A Cartesian 3D inversion-recovery gradient-echo sequence was performed followed by a 3D inversion-recovery gradient-echo with a spiral trajectory. Image quality, fat suppression and sharpness were graded by expert cardiologists. RESULTS: No statistical significance was determined for SNR or image quality between the Cartesian and spiral images; however, fat suppression was significantly improved with the spiral approach (P = 0.002). The enhancement in the Cartesian scan was found to have significantly higher CNR (P = 0.02). CONCLUSION: The feasibility of spiral LGE in the left atrium has been demonstrated. Similar image quality and sharpness were observed with both acquisitions. This spiral sequence has sub-millimeter spatial resolution improved fat suppression and maintains a comparable level of SNR compared with the Cartesian scan.

Respiratory bellows-gated late gadolinium enhancement of the left atrium.

PURPOSE: To compare bellows-gated late gadolinium enhancement (LGE) with standard navigator-gated (NAV-gated) LGE for left atrial (LA) imaging, to eliminate the inflow artifacts associated with NAV-gating. MATERIALS AND METHODS: Eleven subjects, including six patients with atrial fibrillation (AF), were imaged with a 3D free-breathing NAV-gated and bellows-gated LGE. Motion compensation was compared by blinded grading of image sharpness and motion ghosting (0 = worst, 2 = best). Inflow artifacts in the right inferior pulmonary vein (RIPV) and right superior PV (RSPV) were characterized on the same scale (0 = none, 2 = prominent). In patients, each PV was divided into four quadrants circumferentially in order to assess agreement about scar presence on both image sets. RESULTS: Respiratory compensation was not different (1.7 ± 0.5 vs. 1.6 ± 0.5, sharpness, 1.6 ± 0.5 vs. 1.6 ± 0.5, ghosting, P = NS) for bellows- and NAV-gated images. For NAV-gated LGE, inflow artifacts were more prominent in the RSPV than the RIPV (1.2 ± 0.8 vs. 0.7 ± 0.5, P = 0.046). Visually, inflow artifacts both obscured and mimicked the true scar. Disagreement on the presence of scar was found in 18% of the assessed quadrants, with 25% disagreement for RSPV quadrants (P = 0.01). CONCLUSION: Bellows-gated LGE provides similar respiratory compensation as NAV-gating, without inflow artifacts, leading to improved assessment of scar presence.

Cardiac MRI to investigate myocardial scar and coronary venous anatomy using a slow infusion of dimeglumine gadobenate in patients undergoing assessment for cardiac resynchronization therapy.

PURPOSE: To evaluate a cardiac MR (CMR) examination with slow infusion of a high-relaxivity contrast agent to visualize coronary venous anatomy (CVA) and myocardial scar in heart failure patients awaiting cardiac resynchronization therapy (CRT). MATERIALS AND METHODS: Fourteen patients awaiting CRT (seven ischemic cardiomyopathy (ICM) and seven non-ICM) and two with normal LV function underwent CMR on a 1.5 Tesla (T) MR scanner. Dimeglumine-gadobenate was slowly infused. Bolus arrival in the LV was measured by a dynamic electrocardiogram (ECG) -triggered inversion recovery (IR) scan subsequent to starting an ECG-triggered respiratory-navigated three-dimensional (3D) SSFP MR scan with IR preparation to acquire systolic whole-heart anatomy for vein visualization. Delayed contrast-enhanced MR scan was performed to assess myocardial scar. CVA obtained by CMR was compared with X-ray venography in 11 patients. CVA and scar were segmented and registered for visual inspection. RESULTS: For all subjects, there was excellent visualization of the CVA. All ICM and one non-ICM patient showed scar. There was excellent correlation between veins seen by CMR and venography. CONCLUSION: We have demonstrated that slow infusion protocol of dimeglumine-gadobenate can be used to assess both CVA and myocardial scar in a single MR examination. Furthermore, an image overlay technique has been used to show the relationship of scar to the CVA.

Advanced image fusion to overlay coronary sinus anatomy with real-time fluoroscopy to facilitate left ventricular lead implantation in CRT.

BACKGROUND: Failure rate for left ventricular (LV) lead implantation in cardiac resynchronization therapy (CRT) is up to 12%. The use of segmentation tools, advanced image registration software, and high-fidelity images from computerized tomography (CT) and cardiac magnetic resonance (CMR) of the coronary sinus (CS) can guide LV lead implantation. We evaluated the feasibility of advanced image registration onto live fluoroscopic images to allow successful LV lead placement. METHODS: Twelve patients (11 male, 59 ± 16.8 years) undergoing CRT had three-dimensional (3D) whole-heart imaging (six CT, six CMR). Eight patients had at least one previously failed LV lead implant. Using segmentation software, anatomical models of the cardiac chambers, CS, and its branches were overlaid onto the live fluoroscopy using a prototype version of the Philips EP Navigator software to guide lead implantation. RESULTS: We achieved high-fidelity segmentations of cardiac chambers, coronary vein anatomy, and accurate registration between the 3D anatomical models and the live fluoroscopy in all 12 patients confirmed by balloon occlusion angiography. The CS was cannulated successfully in every patient and in 11, an LV lead was implanted successfully. (One patient had no acceptable lead values due to extensive myocardial scar). CONCLUSION: Using overlaid 3D segmentations of the CS and cardiac chambers, it is feasible to guide CRT implantation in real time by fusing advanced imaging and fluoroscopy. This enabled successful CRT in a group of patients with previously failed implants. This technology has the potential to facilitate CRT and improve implant success.

Stability of conventional and machine learning-based tumor auto-segmentation techniques using undersampled dynamic radial bSSFP acquisitions on a 0.35 T hybrid MR-linac system.

PURPOSE: Hybrid MRI-linear accelerator systems (MR-linacs) allow for the incorporation of MR images with high soft-tissue contrast into the radiation therapy procedure prior to, during, or post irradiation. This allows not only for the optimization of the treatment planning, but also for real-time monitoring of the tumor position using cine MRI, from which intrafractional motion can be compensated. Fast imaging and accurate tumor tracking are crucial for effective compensation. This study investigates the application of cine MRI with a radial acquisition scheme on a low-field MR-linac to accelerate the acquisition rate and evaluates the effect on tracking accuracy. METHODS: An MR sequence using tiny golden-angle radial k-space sampling was developed and applied to cine imaging on patients with liver tumors on a 0.35 T MR-linac. Tumor tracking was assessed for accuracy and stability from the cine images with increasing k-space undersampling factors. Tracking was achieved using two different auto-segmentation algorithms: a deformable image registration B-spline similar to that implemented on the MR-linac and a convolutional neural network approach known as U-Net. RESULTS: Radial imaging allows for increased temporal resolution with reliable tumor tracking, although tracking robustness decreases as temporal resolution increases. Additional acquisition-based artifacts can be avoided by reducing the angle increment using tiny golden-angles. The U-net algorithm was found to have superior auto-segmentation metrics compared to B-spline. U-net was able to track two well-defined tumors, imaged with just 30 spokes per image (10.6 frames per second), with an average Dice coefficient ≥ 83%, Hausdorff distance ≤ 1.4 pixel, and mean contour distance ≤ 0.5 pixel. CONCLUSIONS: Radial acquisitions are commonplace in dynamic imaging; however, in MR-guided radiotherapy, robust tumor tracking is also required. This study demonstrates the in vivo feasibility of tumor tracking from radially acquired images on a low-field MR-linac. Radial imaging allows for decreased image acquisition times while maintaining robust tracking. The U-net algorithm can track a tumor with higher accuracy in images with undersampling artifacts than a conventional deformable B-spline algorithm and is a promising tool for tracking in MR-guided radiation therapy.

Beyond T2 and 3T: New MRI techniques for clinicians.

Technological advances in Magnetic Resonance Imaging (MRI) in terms of field strength and hybrid MR systems have led to improvements in tumor imaging in terms of anatomy and functionality. This review paper discusses the applications of such advances in the field of radiation oncology with regards to treatment planning, therapy guidance and monitoring tumor response and predicting outcome.

High-resolution FLAIR MRI at 7 Tesla for treatment planning in glioblastoma patients.

Ultra-high field MRI is an emerging technique promising high-resolution images for radiotherapy planning. We compared a 7 Tesla FLAIR sequence with clinical FLAIR imaging at 3 Tesla in glioblastoma patients before radiotherapy. High-resolution 7 Tesla FLAIR imaging may enhance the depiction of organs at risk and possibly modify target volumes.

Pharmacokinetic modeling of delayed gadolinium enhancement in the myocardium.

Delayed contrast-enhanced magnetic resonance imaging (DCE-MRI) provides prognostic information by delineating regions of myocardial scar. The mechanism of this delayed enhancement in myocardial infarctions (MIs) is hypothesized to result from altered kinetics and changes in the volumes of distribution in the myocardium. Pharmacokinetic models with two and three compartments were fitted to the concentration-time curves of dynamic contrast-enhanced MRI data obtained from five patients with known MI. Furthermore, the parameter stability was investigated in simulations for the two different models. The transfer constants and volumes of distribution showed a good correlation with imaging findings on early and delayed contrast-enhanced MRI. The two compartment model showed higher parameter stability. The three compartment model allows a more in-depth quantification of myocardial scarring. These models have the potential to improve the diagnosis of myocardial pathologies involving scar, with differing kinetics and volumes of distribution such as infarction or cardiomyopathy.

Cardiovascular magnetic resonance imaging of scar development following pulmonary vein isolation: a prospective study.

AIMS: Cardiovascular magnetic resonance (MR) provides non-invasive assessment of early (24-hour) edema and injury following pulmonary vein isolation (by ablation) and subsequent scar formation. We hypothesize that 24-hours after ablation, cardiovascular MR would demonstrate a pattern of edema and injury due to ablation and the severity would correlate with subsequent scar. METHODS: Fifteen atrial fibrillation patients underwent cardiovascular MR prior to pulmonary vein isolation, 24-hours post (N = 11) and 30-days post (N = 7) ablation, with T2-weighted (T2W) and late gadolinium enhancement (LGE) imaging. Left atrial wall thickness, edema enhancement ratio and LGE enhancement were assessed at each time point. Volumes of LGE and edema enhancement were measured, and the circumferential presence of injury was assessed at 24-hours, including comparison with LGE enhancement at 30 days. RESULTS: Left atrial wall thickness was increased 24-hours post-ablation (10.7 ± 4.1 mm vs. 7.0 ± 1.8 mm pre-PVI, p<0.05). T2W enhancement at 24-hours showed increased edema enhancement ratio (1.5 ± 0.4 for post-ablation, vs. 0.9 ± 0.2 pre-ablation, p < 0.001). Edema and LGE volumes at 24-hours were correlated with 30-day LGE volume (R = 0.76, p = 0.04, and R = 0.74, p = 0.09, respectively). Using a 16 segment model for assessment, 24-hour T2W had sensitivity, specificity, and accuracy of 82%, 63%, and 79% respectively, for predicting 30-day LGE. 24-hour LGE had sensitivity, specificity, and accuracy of 91%, 47%, and 84%. CONCLUSIONS: Increased left atrial wall thickening and edema were characterized on cardiovascular MR early post-ablation, and found to correlate with 30-day LGE scar.

3-D visualization of acute RF ablation lesions using MRI for the simultaneous determination of the patterns of necrosis and edema.

Catheter ablation using RF energy is a common treatment for atrial arrhythmias. Although this treatment provides a potential cure, currently, there remains a high proportion of patients returning for repeat ablations. Electrophysiologists have little information to verify that a lesion has been created in the myocardium. Temporary electrical block can be created from edema, which will subside. MRI can visualize acute and chronic ablation lesions using delayed-enhancement techniques. However, the ablation patterns cannot be determined from 2-D images alone. Using the combination of T(2)-weighted and delayed-enhancement MRI, ablation lesions can be characterized in terms of necrosis and edema. A novel 3-D visualization technique is presented that projects the image intensity due the lesions onto a 3-D cardiac surface, allowing the complete, simultaneous visualization of the delayed-enhancement and T(2)-weighted ablation patterns. Results show successful visualization of ablation patterns in 18 patients, and an application of this technique is presented in which electroanatomical mapping systems can be validated by overlaying the acquired ablation points onto the cardiac surfaces and assessing the correlation with the lesion maps.