Fast cardiac T1ρ,adiab mapping using slice-selective adiabatic spin-lock preparation pulses.

PURPOSE: Myocardial T1ρ mapping techniques commonly acquire multiple images in one breathhold to calculate a single-slice T1ρ map. Recently, non-selective adiabatic pulses have been used for robust spin-lock preparation (T1ρ,adiab). The objective of this study was to develop a fast multi-slice myocardial T1ρ,adiab mapping approach. METHODS: The proposed-sequence reduces the number of breathholds required for whole-heart 2D T1ρ,adiab mapping by acquiring multiple interleaved slices in each breathhold using slice-selective T1ρ,adiab preparation pulses. The proposed-sequence was implemented with two interleaved slices per breathhold scan and was quantitatively evaluated in phantom experiments and 10 healthy-volunteers against a single-slice T1ρ,adiab mapping sequence. The sequence was demonstrated in two patients with myocardial scar. RESULTS: The phantom experiments showed the proposed-sequence had slice-to-slice variation of 1.62% ± 1.05% and precision of 4.51 ± 0.68 ms. The healthy volunteer cohort subject-wise mean relaxation time was lower for the proposed-sequence than the single-slice sequence (137.7 ± 5.3 ms vs. 148.4 ± 8.3 ms, p < 0.001), and spatial-standard-deviation was better (18.7 ± 1.8 ms vs. 21.8 ± 3.4 ms, p < 0.018). The mean within-subject, coefficient of variation was 5.93% ± 1.57% for the proposed-sequence and 6.31% ± 1.92% for the single-slice sequence (p = 0.35) and the effect of slice variation (0.81 ± 4.87 ms) was not significantly different to zero (p = 0.61). In both patient examples increased T1ρ,adiab (maximum American Heart Association-segment mean = 174 and 197 ms) was measured within the myocardial scar. CONCLUSION: The proposed sequence provides a twofold acceleration for myocardial T1ρ,adiab mapping using a multi-slice approach. It has no significant difference in within-subject variability, and significantly better precision, compared to a 2D T1ρ,adiab mapping sequence based on non-selective adiabatic spin-lock preparations.

Simultaneous 3D T 1 $$ {\mathrm{T}}_1 $$ , T 2 $$ {\mathrm{T}}_2 $$ , and fat-signal-fraction mapping with respiratory-motion correction for comprehensive liver tissue characterization at 0.55 T.

PURPOSE: To develop a framework for simultaneous three-dimensional (3D) mapping of T 1 $$ {\mathrm{T}}_1 $$ , T 2 $$ {\mathrm{T}}_2 $$ , and fat signal fraction in the liver at 0.55 T. METHODS: The proposed sequence acquires four interleaved 3D volumes with a two-echo Dixon readout. T 1 $$ {\mathrm{T}}_1 $$ and T 2 $$ {\mathrm{T}}_2 $$ are encoded into each volume via preparation modules, and dictionary matching allows simultaneous estimation of T 1 $$ {\mathrm{T}}_1 $$ , T 2 $$ {\mathrm{T}}_2 $$ , and M 0 $$ {M}_0 $$ for water and fat separately. 2D image navigators permit respiratory binning, and motion fields from nonrigid registration between bins are used in a nonrigid respiratory-motion-corrected reconstruction, enabling 100% scan efficiency from a free-breathing acquisition. The integrated nature of the framework ensures the resulting maps are always co-registered. RESULTS: T 1 $$ {\mathrm{T}}_1 $$ , T 2 $$ {\mathrm{T}}_2 $$ , and fat-signal-fraction measurements in phantoms correlated strongly (adjusted r 2 > 0 . 98 $$ {r}^2>0.98 $$ ) with reference measurements. Mean liver tissue parameter values in 10 healthy volunteers were 427 ± 22 $$ 427\pm 22 $$ , 47 . 7 ± 3 . 3 ms $$ 47.7\pm 3.3\;\mathrm{ms} $$ , and 7 ± 2 % $$ 7\pm 2\% $$ for T 1 $$ {\mathrm{T}}_1 $$ , T 2 $$ {\mathrm{T}}_2 $$ , and fat signal fraction, giving biases of 71 $$ 71 $$ , - 30 . 0 ms $$ -30.0\;\mathrm{ms} $$ , and - 5 $$ -5 $$ percentage points, respectively, when compared to conventional methods. CONCLUSION: A novel sequence for comprehensive characterization of liver tissue at 0.55 T was developed. The sequence provides co-registered 3D T 1 $$ {\mathrm{T}}_1 $$ , T 2 $$ {\mathrm{T}}_2 $$ , and fat-signal-fraction maps with full coverage of the liver, from a single nine-and-a-half-minute free-breathing scan. Further development is needed to achieve accurate proton-density fat fraction (PDFF) estimation in vivo.

Real-time automatic image-based slice tracking of gadolinium-filled balloon wedge catheter during MR-guided cardiac catheterization: A proof-of-concept study.

PURPOSE: MR-guided cardiac catheterization procedures currently use passive tracking approaches to follow a gadolinium-filled catheter balloon during catheter navigation. This requires frequent manual tracking and repositioning of the imaging slice during navigation. In this study, a novel framework for automatic real-time catheter tracking during MR-guided cardiac catheterization is presented. METHODS: The proposed framework includes two imaging modes (Calibration and Runtime). The sequence starts in Calibration mode, in which the 3D catheter coordinates are determined using a stack of 10-20 contiguous saturated slices combined with real-time image processing. The sequence then automatically switches to Runtime mode, where three contiguous slices (acquired with partial saturation), initially centered on the catheter balloon using the Calibration feedback, are acquired continuously. The 3D catheter balloon coordinates are estimated in real time from each Runtime slice stack using image processing. Each Runtime stack is repositioned to maintain the catheter balloon in the central slice based on the prior Runtime feedback. The sequence switches back to Calibration mode if the catheter is not detected. This framework was evaluated in a heart phantom and 3 patients undergoing MR-guided cardiac catheterization. Catheter detection accuracy and rate of catheter visibility were evaluated. RESULTS: The automatic detection accuracy for the catheter balloon during the Calibration/Runtime mode was 100%/95% in phantom and 100%/97 ± 3% in patients. During Runtime, the catheter was visible in 82% and 98 ± 2% of the real-time measurements in the phantom and patients, respectively. CONCLUSION: The proposed framework enabled real-time continuous automatic tracking of a gadolinium-filled catheter balloon during MR-guided cardiac catheterization.

Referenceless Nyquist ghost correction outperforms standard navigator-based method and improves efficiency of in vivo diffusion tensor cardiovascular magnetic resonance.

PURPOSE: The study aims to assess the potential of referenceless methods of EPI ghost correction to accelerate the acquisition of in vivo diffusion tensor cardiovascular magnetic resonance (DT-CMR) data using both computational simulations and data from in vivo experiments. METHODS: Three referenceless EPI ghost correction methods were evaluated on mid-ventricular short axis DT-CMR spin echo and STEAM datasets from 20 healthy subjects at 3T. The reduced field of view excitation technique was used to automatically quantify the Nyquist ghosts, and DT-CMR images were fit to a linear ghost model for correction. RESULTS: Numerical simulation showed the singular value decomposition (SVD) method is the least sensitive to noise, followed by Ghost/Object method and entropy-based method. In vivo experiments showed significant ghost reduction for all correction methods, with referenceless methods outperforming navigator methods for both spin echo and STEAM sequences at b = 32, 150, 450, and 600 smm - 2 $$ {\mathrm{smm}}^{-2} $$ . It is worth noting that as the strength of the diffusion encoding increases, the performance gap between the referenceless method and the navigator-based method diminishes. CONCLUSION: Referenceless ghost correction effectively reduces Nyquist ghost in DT-CMR data, showing promise for enhancing the accuracy and efficiency of measurements in clinical practice without the need for any additional reference scans.

Non-rigid motion-compensated 3D whole-heart T2 mapping in a hybrid 3T PET-MR system.

PURPOSE: Simultaneous PET-MRI improves inflammatory cardiac disease diagnosis. However, challenges persist in respiratory motion and mis-registration between free-breathing 3D PET and 2D breath-held MR images. We propose a free-breathing non-rigid motion-compensated 3D T2 -mapping sequence enabling whole-heart myocardial tissue characterization in a hybrid 3T PET-MR system and provides non-rigid respiratory motion fields to correct also simultaneously acquired PET data. METHODS: Free-breathing 3D whole-heart T2 -mapping was implemented on a hybrid 3T PET-MRI system. Three datasets were acquired with different T2 -preparation modules (0, 28, 55 ms) using 3-fold undersampled variable-density Cartesian trajectory. Respiratory motion was estimated via virtual 3D image navigators, enabling multi-contrast non-rigid motion-corrected MR reconstruction. T2 -maps were computed using dictionary-matching. Approach was tested in phantom, 8 healthy subjects, 14 MR only and 2 PET-MR patients with suspected cardiac disease and compared with spin echo reference (phantom) and clinical 2D T2 -mapping (in-vivo). RESULTS: Phantom results show a high correlation (R2 = 0.996) between proposed approach and gold standard 2D T2 mapping. In-vivo 3D T2 -mapping average values in healthy subjects (39.0 ± 1.4 ms) and patients (healthy tissue) (39.1 ± 1.4 ms) agree with conventional 2D T2 -mapping (healthy = 38.6 ± 1.2 ms, patients = 40.3 ± 1.7 ms). Bland-Altman analysis reveals bias of 1.8 ms and 95% limits of agreement (LOA) of -2.4-6 ms for healthy subjects, and bias of 1.3 ms and 95% LOA of -1.9 to 4.6 ms for patients. CONCLUSION: Validated efficient 3D whole-heart T2 -mapping at hybrid 3T PET-MRI provides myocardial inflammation characterization and non-rigid respiratory motion fields for simultaneous PET data correction. Comparable T2 values were achieved with both 3D and 2D methods. Improved image quality was observed in the PET images after MR-based motion correction.

Characterization of quantitative susceptibility mapping in the left ventricular myocardium.

BACKGROUND: Myocardial quantitative susceptibility mapping (QSM) may offer better specificity to iron than conventional T2* imaging in the assessment of cardiac diseases, including intra-myocardial hemorrhage. However, the precision and repeatability of cardiac QSM have not yet been characterized. The aim of this study is to characterize these key metrics in a healthy volunteer cohort and show the feasibility of the method in patients. METHODS: Free breathing respiratory-navigated multi-echo 3D gradient echo images were acquired, from which QSM maps were reconstructed using the Morphology Enhanced Dipole Inversion toolbox. This technique was first evaluated in a susceptibility phantom containing tubes with known concentrations of gadolinium. In vivo characterization of myocardial QSM was then performed in a cohort of 10 healthy volunteers where each subject was scanned twice. Mean segment susceptibility, precision (standard deviation of voxel magnetic susceptibilities within one segment), and repeatability (absolute difference in segment mean susceptibility between repeats) of QSM were calculated for each American Heart Association (AHA) myocardial segment. Finally, the feasibility of the method was shown in 10 patients, including four with hemorrhagic infarcts. RESULTS: The phantom experiment showed a strong linear relationship between measured and predicted susceptibility shifts (R2 > 0.99). For the healthy volunteer cohort, AHA segment analysis showed the mean segment susceptibility was 0.00 ± 0.02 ppm, the mean precision was 0.05 ± 0.04 ppm, and the mean repeatability was 0.02 ± 0.02 ppm. Cardiac QSM was successfully performed in all patients. Focal iron deposits were successfully visualized in the patients with hemorrhagic myocardial infarctions. CONCLUSION: The precision and repeatability of cardiac QSM were successfully characterized in phantom and in vivo experiments. The feasibility of the technique was also successfully demonstrated in patients. While challenges still remain, further clinical evaluation of the technique is now warranted. TRIAL REGISTRATION: This work does not report on a health care intervention.

Highly efficient free-breathing 3D whole-heart imaging in 3-min: single center study in adults with congenital heart disease.

BACKGROUND: Three dimensional, whole-heart (3DWH) MRI is an established non-invasive imaging modality in patients with congenital heart disease (CHD) for the diagnosis of cardiovascular morphology and for clinical decision making. Current techniques utilise diaphragmatic navigation (dNAV) for respiratory motion correction and gating and are frequently limited by long acquisition times. This study proposes and evaluates the diagnostic performance of a respiratory gating-free framework, which considers respiratory image-based navigation (iNAV), and highly accelerated variable density Cartesian sampling in concert with non-rigid motion correction and low-rank patch-based denoising (iNAV-3DWH-PROST). The method is compared to the clinical dNAV-3DWH sequence in adult patients with CHD. METHODS: In this prospective single center study, adult patients with CHD who underwent the clinical dNAV-3DWH MRI were also scanned with the iNAV-3DWH-PROST. Diagnostic confidence (4-point Likert scale) and diagnostic accuracy for common cardiovascular lesions was assessed by three readers. Scan times and diagnostic confidence were compared using the Wilcoxon-signed rank test. Co-axial vascular dimensions at three anatomic landmarks were measured, and agreement between the research and the corresponding clinical sequence was assessed with Bland-Altman analysis. RESULTS: The study included 60 participants (mean age ± [SD]: 33 ± 14 years; 36 men). The mean acquisition time of iNAV-3DWH-PROST was significantly lower compared with the conventional clinical sequence (3.1 ± 0.9 min vs 13.9 ± 3.9 min, p < 0.0001). Diagnostic confidence was higher for the iNAV-3DWH-PROST sequence compared with the clinical sequence (3.9 ± 0.2 vs 3.4 ± 0.8, p < 0.001), however there was no significant difference in diagnostic accuracy. Narrow limits of agreement and mean bias less than 0.08 cm were found between the research and the clinical vascular measurements. CONCLUSIONS: The iNAV-3DWH-PROST framework provides efficient, high quality and robust 3D whole-heart imaging in significantly shorter scan time compared to the standard clinical sequence.

Combining generative modelling and semi-supervised domain adaptation for whole heart cardiovascular magnetic resonance angiography segmentation.

BACKGROUND: Quantification of three-dimensional (3D) cardiac anatomy is important for the evaluation of cardiovascular diseases. Changes in anatomy are indicative of remodeling processes as the heart tissue adapts to disease. Although robust segmentation methods exist for computed tomography angiography (CTA), few methods exist for whole-heart cardiovascular magnetic resonance angiograms (CMRA) which are more challenging due to variable contrast, lower signal to noise ratio and a limited amount of labeled data. METHODS: Two state-of-the-art unsupervised generative deep learning domain adaptation architectures, generative adversarial networks and variational auto-encoders, were applied to 3D whole heart segmentation of both conventional (n = 20) and high-resolution (n = 45) CMRA (target) images, given segmented CTA (source) images for training. An additional supervised loss function was implemented to improve performance given 10%, 20% and 30% segmented CMRA cases. A fully supervised nn-UNet trained on the given CMRA segmentations was used as the benchmark. RESULTS: The addition of a small number of segmented CMRA training cases substantially improved performance in both generative architectures in both standard and high-resolution datasets. Compared with the nn-UNet benchmark, the generative methods showed substantially better performance in the case of limited labelled cases. On the standard CMRA dataset, an average 12% (adversarial method) and 10% (variational method) improvement in Dice score was obtained. CONCLUSIONS: Unsupervised domain-adaptation methods for CMRA segmentation can be boosted by the addition of a small number of supervised target training cases. When only few labelled cases are available, semi-supervised generative modelling is superior to supervised methods.

Evaluation of Deep Learning Estimation of Whole Heart Anatomy from Automated Cardiovascular Magnetic Resonance Short- and Long-Axis Analyses in UK Biobank.

AIMS: Standard methods of heart chamber volume estimation in cardiovascular magnetic resonance (CMR) typically utilize simple geometric formulae based on a limited number of slices. We aimed to evaluate whether an automated deep learning neural network prediction of 3D anatomy of all four chambers would show stronger associations with cardiovascular risk factors and disease than standard volume estimation methods in the UK Biobank. METHODS: A deep learning network was adapted to predict 3D segmentations of left and right ventricles (LV, RV) and atria (LA, RA) at ∼1mm isotropic resolution from CMR short and long axis 2D segmentations obtained from a fully automated machine learning pipeline in 4723 individuals with cardiovascular disease (CVD) and 5733 without in the UK Biobank. Relationships between volumes at end-diastole (ED) and end-systole (ES) and risk/disease factors were quantified using univariate, multivariate and logistic regression analyses. Strength of association between deep learning volumes and standard volumes was compared using the area under the receiving operator characteristic curve (AUC). RESULTS: Univariate and multivariate associations between deep learning volumes and most risk and disease factors were stronger than for standard volumes (higher R2 and more significant P values), particularly for sex, age, and body mass index. AUC for all logistic regressions were higher for deep learning volumes than standard volumes (p<0.001 for all four chambers at ED and ES). CONCLUSIONS: Neural network reconstructions of whole heart volumes had significantly stronger associations with cardiovascular disease and risk factors than standard volume estimation methods in an automatic processing pipeline.

Magnetic resonance elastography resolving all gross anatomical segments of the kidney during controlled hydration.

Introduction: Magnetic resonance elastography (MRE) is a non-invasive method to quantify biomechanical properties of human tissues. It has potential in diagnosis and monitoring of kidney disease, if established in clinical practice. The interplay of flow and volume changes in renal vessels, tubule, urinary collection system and interstitium is complex, but physiological ranges of in vivo viscoelastic properties during fasting and hydration have never been investigated in all gross anatomical segments simultaneously. Method: Ten healthy volunteers underwent two imaging sessions, one following a 12-hour fasting period and the second after a drinking challenge of >10 mL per kg body weight (60-75 min before the second examination). High-resolution renal MRE was performed using a novel driver with rotating eccentric mass placed at the posterior-lateral wall to couple waves (50 Hz) to the kidney. The biomechanical parameters, shear wave speed (cs in m/s), storage modulus (Gd in kPa), loss modulus (Gl in kPa), phase angle (Υ=2πatanGlGd) and attenuation (α in 1/mm) were derived. Accurate separation of gross anatomical segments was applied in post-processing (whole kidney, cortex, medulla, sinus, vessel). Results: High-quality shear waves coupled into all gross anatomical segments of the kidney (mean shear wave displacement: 163 ± 47 µm, mean contamination of second upper harmonics <23%, curl/divergence: 4.3 ± 0.8). Regardless of the hydration state, median Gd of the cortex and medulla (0.68 ± 0.11 kPa) was significantly higher than that of the sinus and vessels (0.48 ± 0.06 kPa), and consistently, significant differences were found in cs, Υ, and Gl (all p < 0.001). The viscoelastic parameters of cortex and medulla were not significantly different. After hydration sinus exhibited a small but significant reduction in median Gd by -0.02 ± 0.04 kPa (p = 0.01), and, consequently, the cortico-sinusoidal-difference in Gd increased by 0.04 ± 0.07 kPa (p = 0.05). Only upon hydration, the attenuation in vessels became lower (0.084 ± 0.013 1/mm) and differed significantly from the whole kidney (0.095 ± 0.007 1/mm, p = 0.01). Conclusion: High-resolution renal MRE with an innovative driver and well-defined 3D segmentation can resolve all renal segments, especially when including the sinus in the analysis. Even after a prolonged hydration period the approach is sensitive to small hydration-related changes in the sinus and in the cortico-sinusoidal-difference.

Fast reconstruction of SMS bSSFP myocardial perfusion images using noise map estimation network (NoiseMapNet): a head-to-head comparison with parallel imaging and iterative reconstruction.

Background: Simultaneous multi-slice (SMS) bSSFP imaging enables stress myocardial perfusion imaging with high spatial resolution and increased spatial coverage. Standard parallel imaging techniques (e.g., TGRAPPA) can be used for image reconstruction but result in high noise level. Alternatively, iterative reconstruction techniques based on temporal regularization (ITER) improve image quality but are associated with reduced temporal signal fidelity and long computation time limiting their online use. The aim is to develop an image reconstruction technique for SMS-bSSFP myocardial perfusion imaging combining parallel imaging and image-based denoising using a novel noise map estimation network (NoiseMapNet), which preserves both sharpness and temporal signal profiles and that has low computational cost. Methods: The proposed reconstruction of SMS images consists of a standard temporal parallel imaging reconstruction (TGRAPPA) with motion correction (MOCO) followed by image denoising using NoiseMapNet. NoiseMapNet is a deep learning network based on a 2D Unet architecture and aims to predict a noise map from an input noisy image, which is then subtracted from the noisy image to generate the denoised image. This approach was evaluated in 17 patients who underwent stress perfusion imaging using a SMS-bSSFP sequence. Images were reconstructed with (a) TGRAPPA with MOCO (thereafter referred to as TGRAPPA), (b) iterative reconstruction with integrated motion compensation (ITER), and (c) proposed NoiseMapNet-based reconstruction. Normalized mean squared error (NMSE) with respect to TGRAPPA, myocardial sharpness, image quality, perceived SNR (pSNR), and number of diagnostic segments were evaluated. Results: NMSE of NoiseMapNet was lower than using ITER for both myocardium (0.045 ± 0.021 vs. 0.172 ± 0.041, p < 0.001) and left ventricular blood pool (0.025 ± 0.014 vs. 0.069 ± 0.020, p < 0.001). There were no significant differences between all methods for myocardial sharpness (p = 0.77) and number of diagnostic segments (p = 0.36). ITER led to higher image quality than NoiseMapNet/TGRAPPA (2.7 ± 0.4 vs. 1.8 ± 0.4/1.3 ± 0.6, p < 0.001) and higher pSNR than NoiseMapNet/TGRAPPA (3.0 ± 0.0 vs. 2.0 ± 0.0/1.3 ± 0.6, p < 0.001). Importantly, NoiseMapNet yielded higher pSNR (p < 0.001) and image quality (p < 0.008) than TGRAPPA. Computation time of NoiseMapNet was only 20s for one entire dataset. Conclusion: NoiseMapNet-based reconstruction enables fast SMS image reconstruction for stress myocardial perfusion imaging while preserving sharpness and temporal signal profiles.

Inversion recovery and saturation recovery pulmonary vein MR angiography using an image based navigator fluoro trigger and variable-density 3D cartesian sampling with spiral-like order.

Contrast enhanced pulmonary vein magnetic resonance angiography (PV CE-MRA) has value in atrial ablation pre-procedural planning. We aimed to provide high fidelity, ECG gated PV CE-MRA accelerated by variable density Cartesian sampling (VD-CASPR) with image navigator (iNAV) respiratory motion correction acquired in under 4 min. We describe its use in part during the global iodinated contrast shortage. VD-CASPR/iNAV framework was applied to ECG-gated inversion and saturation recovery gradient recalled echo PV CE-MRA in 65 patients (66 exams) using .15 mmol/kg Gadobutrol. Image quality was assessed by three physicians, and anatomical segmentation quality by two technologists. Left atrial SNR and left atrial/myocardial CNR were measured. 12 patients had CTA within 6 months of MRA. Two readers assessed PV ostial measurements versus CTA for intermodality/interobserver agreement. Inter-rater/intermodality reliability, reproducibility of ostial measurements, SNR/CNR, image, and anatomical segmentation quality was compared. The mean acquisition time was 3.58 ± 0.60 min. Of 35 PV pre-ablation datasets (34 patients), mean anatomical segmentation quality score was 3.66 ± 0.54 and 3.63 ± 0.55 as rated by technologists 1 and 2, respectively (p = 0.7113). Good/excellent anatomical segmentation quality (grade 3/4) was seen in 97% of exams. Each rated one exam as moderate quality (grade 2). 95% received a majority image quality score of good/excellent by three physicians. Ostial PV measurements correlated moderate to excellently with CTA (ICCs range 0.52-0.86). No difference in SNR was observed between IR and SR. High quality PV CE-MRA is possible in under 4 min using iNAV bolus timing/motion correction and VD-CASPR.