Real-time fMRI neurofeedback of the anterior insula using arterial spin labelling.

Arterial spin labelling (ASL) is the only non-invasive technique that allows absolute quantification of perfusion and is increasingly used in brain activation studies. Contrary to the blood oxygen level-dependent (BOLD) effect ASL measures the cerebral blood flow (CBF) directly. However, the ASL signal has a lower signal-to-noise ratio (SNR), than the BOLD signal, which constrains its utilization in neurofeedback studies. If successful, ASL neurofeedback can be used to aid in the rehabilitation of health conditions with impaired blood flow, for example, stroke. We provide the first ASL-based neurofeedback study incorporating a double-blind, sham-controlled design. A pseudo-continuous ASL (pCASL) approach with background suppression and 3D GRASE readout was combined with a real-time post-processing pipeline. The real-time pipeline allows to monitor the ASL signal and provides real-time feedback on the neural activity to the subject. In total 41 healthy adults (19-56 years) divided into three groups underwent a neurofeedback-based emotion imagery training of the left anterior insula. Two groups differing only in the explicitness level of instruction received real training and a third group received sham feedback. Only those participants receiving real feedback with explicit instruction showed significantly higher absolute CBF values in the trained region during neurofeedback than participants receiving sham feedback. However, responder analyses of percent signal change values show no differences in activation between the three groups. Persisting limitations, such as the lower SNR, confounding effects of arterial transit time and partial volume effects still impact negatively the implementation of ASL neurofeedback.

Fully Automated Hippocampus Segmentation using T2-informed Deep Convolutional Neural Networks.

Hippocampal atrophy (tissue loss) has become a fundamental outcome parameter in clinical trials on Alzheimer's disease. To accurately estimate hippocampus volume and track its volume loss, a robust and reliable segmentation is essential. Manual hippocampus segmentation is considered the gold standard but is extensive, time-consuming, and prone to rater bias. Therefore, it is often replaced by automated programs like FreeSurfer, one of the most commonly used tools in clinical research. Recently, deep learning-based methods have also been successfully applied to hippocampus segmentation. The basis of all approaches are clinically used T1-weighted whole-brain MR images with approximately 1 mm isotropic resolution. However, such T1 images show low contrast-to-noise ratios (CNRs), particularly for many hippocampal substructures, limiting delineation reliability. To overcome these limitations, high-resolution T2-weighted scans are suggested for better visualization and delineation, as they show higher CNRs and usually allow for higher resolutions. Unfortunately, such time-consuming T2-weighted sequences are not feasible in a clinical routine. We propose an automated hippocampus segmentation pipeline leveraging deep learning with T2-weighted MR images for enhanced hippocampus segmentation of clinical T1-weighted images based on a series of 3D convolutional neural networks and a specifically acquired multi-contrast dataset. This dataset consists of corresponding pairs of T1- and high-resolution T2-weighted images, with the T2 images only used to create more accurate manual ground truth annotations and to train the segmentation network. The T2-based ground truth labels were also used to evaluate all experiments by comparing the masks visually and by various quantitative measures. We compared our approach with four established state-of-the-art hippocampus segmentation algorithms (FreeSurfer, ASHS, HippoDeep, HippMapp3r) and demonstrated a superior segmentation performance. Moreover, we found that the automated segmentation of T1-weighted images benefits from the T2-based ground truth data. In conclusion, this work showed the beneficial use of high-resolution, T2-based ground truth data for training an automated, deep learning-based hippocampus segmentation and provides the basis for a reliable estimation of hippocampal atrophy in clinical studies.

Improved quantification in CEST-MRI by joint spatial total generalized variation.

PURPOSE: In this work, the use of joint Total Generalized Variation (TGV) regularization to improve Multipool-Lorentzian fitting of chemical exchange saturation transfer (CEST) Spectra in terms of stability and parameter signal-to-noise ratio (SNR) was investigated. THEORY AND METHODS: The joint TGV term was integrated into the nonlinear parameter fitting problem. To increase convergence and weight the gradients, preconditioning using a voxel-wise singular value decomposition was applied to the problem, which was then solved using the iteratively regularized Gauss-Newton method combined with a Primal-Dual splitting algorithm. The TGV method was evaluated on simulated numerical phantoms, 3T phantom data and 7T in vivo data with respect to systematic errors and robustness. Three reference methods were also implemented: The standard nonlinear fitting, a method using a nonlocal-means filter for denoising and the pyramid scheme, which uses downsampled images to acquire accurate start values. RESULTS: The proposed regularized fitting method showed significantly improved robustness (compared to the reference methods). In testing, over a range of SNR values the TGV fit outperformed the other methods and showed accurate results even for large amounts of added noise. Parameter values found were closer or comparable to the ground truth. For in vivo datasets, the added regularization increased the parameter map SNR and prevented instabilities. CONCLUSION: The proposed fitting method using TGV regularization leads to improved results over a range of different data-sets and noise levels. Furthermore, it can be applied to all Z-spectrum data, with different amounts of pools, where the improved SNR and stability can increase diagnostic confidence.

Enhanced and robust contrast in CEST MRI: Saturation pulse shape design via optimal control.

PURPOSE: To employ optimal control for the numerical design of Chemical Exchange Saturation Transfer (CEST) saturation pulses to maximize contrast and stability against B 0 $$ {\mathrm{B}}_0 $$ inhomogeneities. THEORY AND METHODS: We applied an optimal control framework for the design pulse shapes for CEST saturation pulse trains. The cost functional minimized both the pulse energy and the discrepancy between the corresponding CEST spectrum and the target spectrum based on a continuous radiofrequency (RF) pulse. The optimization is subject to hardware limitations. In measurements on a 7 T preclinical scanner, the optimal control pulses were compared to continuous-wave and Gaussian saturation methods. We conducted a comparison of the optimal control pulses with Gaussian, block pulse trains, and adiabatic spin-lock pulses. RESULTS: The optimal control pulse train demonstrated saturation levels comparable to continuous-wave saturation and surpassed Gaussian saturation by up to 50 % in phantom measurements. In phantom measurements at 3 T the optimized pulses not only showcased the highest CEST contrast, but also the highest stability against field inhomogeneities. In contrast, block pulse saturation resulted in severe artifacts. Dynamic Bloch-McConnell simulations were employed to identify the source of these artifacts, and underscore the B 0 $$ {\mathrm{B}}_0 $$ robustness of the optimized pulses. CONCLUSION: In this work, it was shown that a substantial improvement in pulsed saturation CEST imaging can be achieved by using Optimal Control design principles. It is possible to overcome the sensitivity of saturation to B0 inhomogeneities while achieving CEST contrast close to continuous wave saturation.

Robust dual-angle T 1 measurement in magnetization transfer spectroscopy by time-optimal control.

Magnetization transfer spectroscopy relies heavily on the robust determination of T 1 relaxation times of nuclei participating in metabolic exchange. Challenges arise due to the use of surface RF coils for transmission (high B 1 + variation) and the broad resonance band of most X nuclei. These challenges are particularly pronounced when fast T 1 mapping methods, such as the dual-angle method, are employed. Consequently, in this work, we develop resonance offset and B 1 + robust excitation RF pulses for 31P magnetization transfer spectroscopy at 7T through ensemble-based time-optimal control. In our approach, we introduce a cost functional for designing robust pulses, incorporating the full Bloch equations as constraints, which are solved using symmetric operator splitting techniques. The optimal control design of the RF pulses developed demonstrates improved accuracy, desired phase properties, and reduced RF power when applied to dual-angle T 1 mapping, thereby improving the precision of exchange-rate measurements, as demonstrated in a preclinical in vivo study quantifying brain creatine kinase activity.

Whole-Cerebrum distortion-free three-dimensional pseudo-continuous arterial spin labeling at 7T.

Fulfilling potentials of ultrahigh field for pseudo-Continuous Arterial Spin Labeling (pCASL) has been hampered by B1/B0 inhomogeneities that affect pCASL labeling, background suppression (BS), and the readout sequence. This study aimed to present a whole-cerebrum distortion-free three-dimensional (3D) pCASL sequence at 7T by optimizing pCASL labeling parameters, BS pulses, and an accelerated Turbo-FLASH (TFL) readout. A new set of pCASL labeling parameters (Gave = 0.4 mT/m, Gratio = 14.67) was proposed to avoid interferences in bottom slices while achieving robust labeling efficiency (LE). An OPTIM BS pulse was designed based on the range of B1/B0 inhomogeneities at 7T. A 3D TFL readout with 2D-CAIPIRINHA undersampling (R = 2 × 2) and centric ordering was developed, and the number of segments (Nseg) and flip angle (FA) were varied in simulation to achieve the optimal trade-off between SNR and spatial blurring. In-vivo experiments were performed on 19 subjects. The results showed that the new set of labeling parameters effectively achieved whole-cerebrum coverage by eliminating interferences in bottom slices while maintaining a high LE. The OPTIM BS pulse achieved 33.3% higher perfusion signal in gray matter (GM) than the original BS pulse with a cost of 4.8-fold SAR. Incorporating a moderate FA (8°) and Nseg (2), whole-cerebrum 3D TFL-pCASL imaging was achieved with a 2 × 2 × 4 mm3 resolution without distortion and susceptibility artifacts compared to 3D GRASE-pCASL. In addition, 3D TFL-pCASL showed a good to excellent test-retest repeatability and potential of higher resolution (2 mm isotropic). The proposed technique also significantly improved SNR when compared to the same sequence at 3T and simultaneous multislice TFL-pCASL at 7T. By combining a new set of labeling parameters, OPTIM BS pulse, and accelerated 3D TFL readout, we achieved high resolution pCASL at 7T with whole-cerebrum coverage, detailed perfusion and anatomical information without distortion, and sufficient SNR.

Whole-Cerebrum distortion-free three-dimensional pseudo-Continuous Arterial Spin Labeling at 7T.

Fulfilling potentials of ultrahigh field for pseudo-Continuous Arterial Spin Labeling (pCASL) has been hampered by B1/B0 inhomogeneities that affect pCASL labeling, background suppression (BS), and the readout sequence. This study aimed to present a whole-cerebrum distortion-free three-dimensional (3D) pCASL sequence at 7T by optimizing pCASL labeling parameters, BS pulses, and an accelerated Turbo-FLASH (TFL) readout. A new set of pCASL labeling parameters (Gave=0.4mT/m, Gratio=14.67) was proposed to avoid interferences in bottom slices while achieving robust labeling efficiency (LE). An OPTIM BS pulse was designed based on the range of B1/B0 inhomogeneities at 7T. A 3D TFL readout with 2D-CAIPIRINHA undersampling (R=2×2) and centric ordering was developed, and the number of segments (Nseg) and flip angle (FA) were varied in simulation to achieve the optimal trade-off between SNR and spatial blurring. In-vivo experiments were performed on 19 subjects. The results showed that the new set of labeling parameters effectively achieved whole-cerebrum coverage by eliminating interferences in bottom slices while maintaining a high LE. The OPTIM BS pulse achieved 33.3% higher perfusion signal in gray matter (GM) than the original BS pulse with a cost of 4.8-fold SAR. Incorporating a moderate FA (8 ° ) and Nseg (2), whole-cerebrum 3D TFL-pCASL imaging was achieved with a 2×2×4 mm 3 resolution without distortion and susceptibility artifacts compared to 3D GRASE-pCASL. In addition, 3D TFL-pCASL showed a good to excellent test-retest repeatability and potential of higher resolution (2 mm isotropic). The proposed technique also significantly improved SNR when compared to the same sequence at 3T and simultaneous multislice TFL-pCASL at 7T. By combining a new set of labeling parameters, OPTIM BS pulse, and accelerated 3D TFL readout, we achieved high resolution pCASL at 7T with whole-cerebrum coverage, detailed perfusion and anatomical information without distortion, and sufficient SNR.

Model-based reconstructions for intravoxel incoherent motion and diffusion tensor imaging parameter map estimations.

Intravoxel incoherent motion (IVIM) imaging and diffusion tensor imaging (DTI) facilitate noninvasive quantification of tissue perfusion and diffusion. Both are promising biomarkers in various diseases and a combined acquisition is therefore desirable. This comes with challenges, including noisy parameter maps and long scan times, especially for the perfusion fraction f and pseudo-diffusion coefficient D*. A model-based reconstruction has the potential to overcome these challenges. As a first step, our goal was to develop a model-based reconstruction framework for IVIM and combined IVIM-DTI parameter estimation. The IVIM and IVIM-DTI models were implemented in the PyQMRI model-based reconstruction framework and validated with simulations and in vivo data. Commonly used voxel-wise nonlinear least-squares fitting was used as the reference. Simulations with the IVIM and IVIM-DTI models were performed with 100 noise realizations to assess accuracy and precision. Diffusion-weighted data were acquired for IVIM reconstruction in the liver (n = 5), as well as for IVIM-DTI in the kidneys (n = 5) and lower-leg muscles (n = 6) of healthy volunteers. The median and interquartile range (IQR) values of the IVIM and IVIM-DTI parameters were compared to assess bias and precision. With model-based reconstruction, the parameter maps exhibited less noise, which was most pronounced in the f and D* maps, both in the simulations and in vivo. The bias values in the simulations were comparable between model-based reconstruction and the reference method. The IQR was lower with model-based reconstruction compared with the reference for all parameters. In conclusion, model-based reconstruction is feasible for IVIM and IVIM-DTI and improves the precision of the parameter estimates, particularly for f and D* maps.

Interpretable brain disease classification and relevance-guided deep learning.

Deep neural networks are increasingly used for neurological disease classification by MRI, but the networks' decisions are not easily interpretable by humans. Heat mapping by deep Taylor decomposition revealed that (potentially misleading) image features even outside of the brain tissue are crucial for the classifier's decision. We propose a regularization technique to train convolutional neural network (CNN) classifiers utilizing relevance-guided heat maps calculated online during training. The method was applied using T1-weighted MR images from 128 subjects with Alzheimer's disease (mean age = 71.9 ± 8.5 years) and 290 control subjects (mean age = 71.3 ± 6.4 years). The developed relevance-guided framework achieves higher classification accuracies than conventional CNNs but more importantly, it relies on less but more relevant and physiological plausible voxels within brain tissue. Additionally, preprocessing effects from skull stripping and registration are mitigated. With the interpretability of the decision mechanisms underlying CNNs, these results challenge the notion that unprocessed T1-weighted brain MR images in standard CNNs yield higher classification accuracy in Alzheimer's disease than solely atrophy.

Advanced design of MRI inversion pulses for inhomogeneous field conditions by optimal control.

Non-selective inversion pulses find widespread use in MRI applications, where requirements on them are increasingly demanding. With the use of high and ultra-high field strength systems, robustness to Δ B 0 and B 1 + inhomogeneities, while tackling SAR and hardware limitations, has rapidly become important. In this work, we propose a time-optimal control framework for the optimization of Δ B 0 - and B 1 + -robust inversion pulses. Robustness is addressed by means of ensemble formulations, while allowing inclusion of hardware and energy limitations. The framework is flexible and performs excellently for various optimization goals. The optimization results are analyzed extensively in numerical experiments. Furthermore, they are validated, and compared with adiabatic RF pulses, in various phantom and in vivo measurements on a 3 T MRI system.

Periventricular magnetisation transfer abnormalities in early multiple sclerosis.

OBJECTIVE: Recent studies suggested that CSF-mediated factors contribute to periventricular (PV) T2-hyperintense lesion formation in multiple sclerosis (MS) and this in turn correlates with cortical damage. We thus investigated if such PV-changes are observable microstructurally in early-MS and if they correlate with cortical damage. METHODS: We assessed the magnetisation transfer ratio (MTR) in PV normal-appearing white matter (NAWM) and in MS lesions in 44 patients with a clinically isolated syndrome (CIS) suggestive of MS and 73 relapsing-remitting MS (RRMS) patients. Band-wise MTR values were related to cortical mean thickness (CMT) and compared with 49 healthy controls (HCs). For each band, MTR changes were assessed relative to the average MTR values of all HCs. RESULTS: Relative to HCs, PV-MTR was significantly reduced up to 2.63% in CIS and 5.37% in RRMS (p < 0.0001). The MTR decreased towards the lateral ventricles with 0.18%/mm in CIS and 0.31%/mm in RRMS patients, relative to HCs. In RRMS, MTR-values adjacent to the ventricle and in PV-lesions correlated positively with CMT and negatively with EDSS. CONCLUSION: PV-MTR gradients are present from the earliest stage of MS, consistent with more pronounced microstructural WM-damage closer to the ventricles. The positive association between reduced CMT and lower MTR in PV-NAWM suggests a common pathophysiologic mechanism. Together, these findings indicate the potential use of multimodal MRI as refined marker for MS-related tissue changes.

Automated vortical blood flow-based estimation of mean pulmonary arterial pressure from 4D flow MRI.

PURPOSE: Elevated mean pulmonary arterial pressure (mPAP), or pulmonary hypertension (PH), is associated with vortical blood flow along the main pulmonary artery. We present and validate a method for automated detection and tracking of the PH-related vortex from magnetic resonance 4D flow data that allows estimation of mPAP. METHODS: The proposed method detects the presence of a PH-related vortex in the main pulmonary artery based on geometrical properties of swirling streamlines and estimates mPAP from the PH-related vortex duration (tvortex) using a previously established model. 4D flow data of 32 subjects (19/13 with/without PH) who underwent right heart catheterization (RHC) for mPAP measurement and diagnosis of PH (mPAP >20 mmHg) were used to compare visual and automated PH-related vortex detection and to validate estimated mPAP against RHC-derived results. RESULTS: Visually and automatically determined tvortex values correlated strongly (r = 0.98); they yielded no bias, and the standard deviation of differences between them was small (5.9% of the cardiac interval). mPAP estimates from visual and automated analyses both allowed diagnosis of PH with an area under the curve of 1.00 [0.89,1.00]. For subjects with PH, neither visually nor automatically estimated mPAP differed from mPAP measured by RHC, while the standard deviation between estimated and invasively measured mPAP was lower with visual estimation (3.1 mmHg vs. 5.3 mmHg). CONCLUSION: An automated method for PH-related vortex detection and tracking from magnetic resonance 4D flow data was introduced, which demonstrated very good agreement with visual analysis and accurate estimation of elevated mPAP.

Accuracy and performance analysis for Bloch and Bloch-McConnell simulation methods.

PURPOSE: To introduce new solution methods for the Bloch and Bloch-McConnell equations and compare them quantitatively to different known approaches. THEORY AND METHODS: A new exact solution per time step is derived by means of eigenvalues and generalized eigenvectors. Fast numerical solution methods based on asymmetric and symmetric operator splitting, which are already known for the Bloch equations, are extended to the Bloch-McConnell equations. Those methods are compared to other numerical methods including spin domain, one-step and multi-step methods, and matrix exponential. Error metrics are introduced based on the exact solution method, which allows to assess the accuracy of each solution method quantitatively for arbitrary example data. RESULTS: Accuracy and performance properties for nine different solution methods are analyzed and compared in extensive numerical experiments including various examples for non-selective and slice-selective MR imaging applications. The accuracy of the methods heavily varies, in particular for short relaxation times and long pulse durations. CONCLUSION: In absence of relaxation effects, the numerical results confirm the rotation matrices approach as accurate and computationally efficient Bloch solution method. Otherwise, as well as for the Bloch-McConnell equations, symmetric operator splitting methods are recommended due to their excellent numerical accuracy paired with efficient run time.

Joint multi-field T1 quantification for fast field-cycling MRI.

PURPOSE: Recent developments in hardware design enable the use of fast field-cycling (FFC) techniques in MRI to exploit the different relaxation rates at very low field strength, achieving novel contrast. The method opens new avenues for in vivo characterizations of pathologies but at the expense of longer acquisition times. To mitigate this, we propose a model-based reconstruction method that fully exploits the high information redundancy offered by FFC methods. METHODS: The proposed model-based approach uses joint spatial information from all fields by means of a Frobenius - total generalized variation regularization. The algorithm was tested on brain stroke images, both simulated and acquired from FFC patients scans using an FFC spin echo sequences. The results are compared to three non-linear least squares fits with progressively increasing complexity. RESULTS: The proposed method shows excellent abilities to remove noise while maintaining sharp image features with large signal-to-noise ratio gains at low-field images, clearly outperforming the reference approach. Especially patient data show huge improvements in visual appearance over all fields. CONCLUSION: The proposed reconstruction technique largely improves FFC image quality, further pushing this new technology toward clinical standards.

Non-linear fitting with joint spatial regularization in arterial spin labeling.

Multi-Delay single-shot arterial spin labeling (ASL) imaging provides accurate cerebral blood flow (CBF) and, in addition, arterial transit time (ATT) maps but the inherent low SNR can be challenging. Especially standard fitting using non-linear least squares often fails in regions with poor SNR, resulting in noisy estimates of the quantitative maps. State-of-the-art fitting techniques improve the SNR by incorporating prior knowledge in the estimation process which typically leads to spatial blurring. To this end, we propose a new estimation method with a joint spatial total generalized variation regularization on CBF and ATT. This joint regularization approach utilizes shared spatial features across maps to enhance sharpness and simultaneously improves noise suppression in the final estimates. The proposed method is evaluated at three levels, first on synthetic phantom data including pathologies, followed by in vivo acquisitions of healthy volunteers, and finally on patient data following an ischemic stroke. The quantitative estimates are compared to two reference methods, non-linear least squares fitting and a state-of-the-art ASL quantification algorithm based on Bayesian inference. The proposed joint regularization approach outperforms the reference implementations, substantially increasing the SNR in CBF and ATT while maintaining sharpness and quantitative accuracy in the estimates.

Impact of the Choice of Native T1 in Pixelwise Myocardial Blood Flow Quantification.

BACKGROUND: Quantification of myocardial blood flow (MBF) from dynamic contrast-enhanced (DCE) MRI can be performed using a signal intensity model that incorporates T1 values of blood and myocardium. PURPOSE: To assess the impact of T1 values on pixelwise MBF quantification, specifically to evaluate the influence of 1) study population-averaged vs. subject-specific, 2) diastolic vs. systolic, and 3) regional vs. global myocardial T1 values. STUDY TYPE: Prospective. SUBJECTS: Fifteen patients with chronic coronary heart disease. FIELD STRENGTH/SEQUENCE: 3T; modified Look-Locker inversion recovery for T1 mapping and saturation recovery gradient echo for DCE imaging, both acquired in a mid-ventricular short-axis slice in systole and diastole. ASSESSMENT: MBF was estimated using Fermi modeling and signal intensity nonlinearity correction with different T1 values: study population-averaged blood and myocardial, subject-specific systolic and diastolic, and segmental T1 values. Myocardial segments with perfusion deficits were identified visually from DCE series. STATISTICAL TESTS: The relationships between MBF parameters derived by different methods were analyzed by Bland-Altman analysis; corresponding mean values were compared by t-test. RESULTS: Using subject-specific diastolic T1 values, global diastolic MBF was 0.61 ± 0.13 mL/(min·g). It did not differ from global MBF derived from the study population-averaged T1 (P = 0.88), but the standard deviation of differences was large (0.07 mL/(min·g), 11% of mean MBF). Global diastolic and systolic MBF did not differ (P = 0.12), whereas global diastolic MBF using systolic (0.62 ± 0.13 mL/(min·g)) and diastolic T1 values differed (P < 0.05). If regional instead of global T1 values were used, segmental MBF was lower in segments with perfusion deficits (bias = -0.03 mL/(min·g), -7% of mean MBF, P < 0.05) but higher in segments without perfusion deficits (bias = 0.01 mL/(min·g), 1% of mean MBF, P < 0.05). DATA CONCLUSION: Whereas cardiac phase-specific T1 values have a minor impact on MBF estimates, subject-specific and myocardial segment-specific T1 values substantially affect MBF quantification. LEVEL OF EVIDENCE: 3 TECHNICAL EFFICACY STAGE: 3.

Adaptive slice-specific z-shimming for 2D spoiled gradient-echo sequences.

PURPOSE: To reduce the misbalance between compensation gradients and macroscopic field gradients, we introduce an adaptive slice-specific z-shimming approach for 2D spoiled multi-echo gradient-echoe sequences in combination with modeling of the signal decay. METHODS: Macroscopic field gradients were estimated for each slice from a fast prescan (15 seconds) and then used to calculate slice-specific compensation moments along the echo train. The coverage of the compensated field gradients was increased by applying three positive and three negative moments. With a forward model, which considered the effect of the slice profile, the z-shim moment, and the field gradient, R2∗ maps were estimated. The method was evaluated in phantom and in vivo measurements at 3 T and compared with a spoiled multi-echo gradient-echo and a global z-shimming approach without slice-specific compensation. RESULTS: The proposed method yielded higher SNR in R2∗ maps due to a broader range of compensated macroscopic field gradients compared with global z-shimming. In global white matter, the mean interquartile range, proxy for SNR, could be decreased to 3.06 s-1 with the proposed approach, compared with 3.37 s-1 for global z-shimming and 3.52 s-1 for uncompensated multi-echo gradient-echo. CONCLUSION: Adaptive slice-specific compensation gradients between echoes substantially improved the SNR of R2∗ maps, and the signal could also be rephased in anatomical areas, where it has already been completely dephased.

Automated mitral valve vortex ring extraction from 4D-flow MRI.

PURPOSE: To present and validate a method for automated extraction and analysis of the temporal evolution of the mitral valve (MV) vortex ring from MR 4D-flow data. METHODS: The proposed algorithm uses the divergence-free part of the velocity vector field for Q criterion-based identification and tracking of MV vortex ring core and region within the left ventricle (LV). The 4D-flow data of 20 subjects (10 healthy controls, 10 patients with ischemic heart disease) were used to validate the algorithm against visual analysis as well as to assess the method's sensitivity to manual LV segmentation. Quantitative MV vortex ring parameters were analyzed with respect to both their differences between healthy subjects and patients and their correlation with transmitral peak velocities. RESULTS: The algorithm successfully extracted MV vortex rings throughout the entire cardiac cycle, which agreed substantially with visual analysis (Cohen's kappa = 0.77). Furthermore, vortex cores and regions were robustly detected even if a static end-diastolic LV segmentation mask was applied to all frames (Dice coefficients 0.82 ± 0.08 and 0.94 ± 0.02 for core and region, respectively). Early diastolic MV vortex ring vorticity, kinetic energy and circularity index differed significantly between healthy controls and patients. In contrast to vortex shape parameters, vorticity and kinetic energy correlated strongly with transmitral peak velocities. CONCLUSION: An automated method for temporal MV vortex ring extraction demonstrating robustness with respect to LV segmentation strategies is introduced. Quantitative vortex parameter analysis indicates importance of the MV vortex ring for LV diastolic (dys)function.

T2 and T2∗ mapping in ex situ porcine myocardium: myocardial intravariability, temporal stability and the effects of complete coronary occlusion.

Diagnosis of ischaemia-related sudden cardiac death in the absence of microscopic and macroscopic ischaemic lesions remains a challenge for medical examiners. Medical imaging techniques increasingly provide support in post-mortem examinations by detecting and documenting internal findings prior to autopsy. Previous studies have characterised MR relaxation times to investigate post-mortem signs of myocardial infarction in forensic cohorts. In this prospective study based on an ex situ porcine heart model, we report fundamental findings related to intramyocardial variability and temporal stability of T2 as well as the effects of permanent coronary occlusion on T2 and T2∗ relaxation in post-mortem myocardium. The ex situ porcine hearts included in this study (n= 19) were examined in two groups (Ss, n= 11 and Si, n= 8). All magnetic resonance imaging (MRI) examinations were performed ex situ, at room temperature and at 3 T. In the Ss group, T2 mapping was performed on slaughterhouse porcine hearts at different post-mortem intervals (PMI) between 7 and 26 h. Regarding the intramyocardial variability, no statistically significant differences in T2 were observed between myocardial segments (p= 0.167). Assessment of temporal stability indicated a weak negative correlation (r=- 0.21) between myocardial T2 and PMI. In the Si group, animals underwent ethanol-induced complete occlusion of the left anterior descending artery. T2 and T2∗ mapping were performed within 3 h of death. Differences between the expected ischaemic and remote regions were statistically significant for T2 (p= 0.007), however not for T2∗ (p= 0.062). Our results provide important information for future assessment of the diagnostic potential of quantitative MRI in the post-mortem detection of early acute myocardial infarction.

Assessment and correction of macroscopic field variations in 2D spoiled gradient-echo sequences.

PURPOSE: To model and correct the dephasing effects in the gradient-echo signal for arbitrary RF excitation pulses with large flip angles in the presence of macroscopic field variations. METHODS: The dephasing of the spoiled 2D gradient-echo signal was modeled using a numerical solution of the Bloch equations to calculate the magnitude and phase of the transverse magnetization across the slice profile. Additionally, regional variations of the transmit RF field and slice profile scaling due to macroscopic field gradients were included. Simulations, phantom, and in vivo measurements at 3 T were conducted for R2∗ and myelin water fraction (MWF) mapping. RESULTS: The influence of macroscopic field gradients on R2∗ and myelin water fraction estimation can be substantially reduced by applying the proposed model. Moreover, it was shown that the dephasing over time for flip angles of 60° or greater also depends on the polarity of the slice-selection gradient because of phase variation along the slice profile. CONCLUSION: Substantial improvements in R2∗ accuracy and myelin water fraction mapping coverage can be achieved using the proposed model if higher flip angles are required. In this context, we demonstrated that the phase along the slice profile and the polarity of the slice-selection gradient are essential for proper modeling of the gradient-echo signal in the presence of macroscopic field variations.