Connectome spectrum electromagnetic tomography: A method to reconstruct electrical brain source networks at high-spatial resolution.

Connectome spectrum electromagnetic tomography (CSET) combines diffusion MRI-derived structural connectivity data with well-established graph signal processing tools to solve the M/EEG inverse problem. Using simulated EEG signals from fMRI responses, and two EEG datasets on visual-evoked potentials, we provide evidence supporting that (i) CSET captures realistic neurophysiological patterns with better accuracy than state-of-the-art methods, (ii) CSET can reconstruct brain responses more accurately and with more robustness to intrinsic noise in the EEG signal. These results demonstrate that CSET offers high spatio-temporal accuracy, enabling neuroscientists to extend their research beyond the current limitations of low sampling frequency in functional MRI and the poor spatial resolution of M/EEG.

Intra-bin correction and inter-bin compensation of respiratory motion in free-running five-dimensional whole-heart magnetic resonance imaging.

BACKGROUND: Free-running cardiac and respiratory motion-resolved whole-heart five-dimensional (5D) cardiovascular magnetic resonance (CMR) can reduce scan planning and provide a means of evaluating respiratory-driven changes in clinical parameters of interest. However, respiratory-resolved imaging can be limited by user-defined parameters which create trade-offs between residual artifact and motion blur. In this work, we develop and validate strategies for both correction of intra-bin and compensation of inter-bin respiratory motion to improve the quality of 5D CMR. METHODS: Each component of the reconstruction framework was systematically validated and compared to the previously established 5D approach using simulated free-running data (N = 50) and a cohort of 32 patients with congenital heart disease. The impact of intra-bin respiratory motion correction was evaluated in terms of image sharpness while inter-bin respiratory motion compensation was evaluated in terms of reconstruction error, compression of respiratory motion, and image sharpness. The full reconstruction framework (intra-acquisition correction and inter-acquisition compensation of respiratory motion [IIMC] 5D) was evaluated in terms of image sharpness and scoring of image quality by expert reviewers. RESULTS: Intra-bin motion correction provides significantly (p < 0.001) sharper images for both simulated and patient data. Inter-bin motion compensation results in significant (p < 0.001) lower reconstruction error, lower motion compression, and higher sharpness in both simulated (10/11) and patient (9/11) data. The combined framework resulted in significantly (p < 0.001) sharper IIMC 5D reconstructions (End-expiration (End-Exp): 0.45 ± 0.09, End-inspiration (End-Ins): 0.46 ± 0.10) relative to the previously established 5D implementation (End-Exp: 0.43 ± 0.08, End-Ins: 0.39 ± 0.09). Similarly, image scoring by three expert reviewers was significantly (p < 0.001) higher using IIMC 5D (End-Exp: 3.39 ± 0.44, End-Ins: 3.32 ± 0.45) relative to 5D images (End-Exp: 3.02 ± 0.54, End-Ins: 2.45 ± 0.52). CONCLUSION: The proposed IIMC reconstruction significantly improves the quality of 5D whole-heart MRI. This may be exploited for higher resolution or abbreviated scanning. Further investigation of the diagnostic impact of this framework and comparison to gold standards is needed to understand its full clinical utility, including exploration of respiratory-driven changes in physiological measurements of interest.

Trajectory correction enables free-running chemical shift encoded imaging for accurate cardiac proton-density fat fraction quantification at 3T.

BACKGROUND: Metabolic diseases can negatively alter epicardial fat accumulation and composition, which can be probed using quantitative cardiac chemical shift encoded(CSE) MRI by mapping proton-density fat fraction (PDFF). To obtain motion-resolved high-resolution PDFF maps, we proposed a free-running cardiac CSE-MRI framework at 3T. To employ faster bipolar readout gradients, a correction for gradients imperfections was added using the gradient impulse response function (GIRF) and evaluated on intermediate images and PDFF quantification. METHODS: Ten minutes free-running cardiac 3D radial CSE-MRI acquisitions were compared in vitro and in vivo at 3T. Monopolar and bipolar readout gradients schemes provided 8 echoes (TE1/ΔTE = 1.16/1.96ms) and 13 echoes (TE1/ΔTE = 1.12/1.07ms), respectively. Bipolar-gradients free-running cardiac fat and water images and PDFF maps were reconstructed with or without GIRF-correction. PDFF values were evaluated in silico, in vitro on a fat/water phantom, and in vivo in 10 healthy volunteers and three diabetic patients. RESULTS: In monopolar mode, fat-water swaps were demonstrated in silico and confirmed in vitro. Using bipolar readout gradients, PDFF quantification was reliable and accurate with GIRF correction with a mean bias of 0.03% in silico and 0.36% in vitro while it suffered from artifacts without correction, leading to a PDFF bias of 4.9% in vitro and swaps in vivo. Using bipolar readout gradients, in vivo PDFF of epicardial adipose tissue was significantly lower than in subcutaneous fat (80.4±7.1% vs 92.5±4.3%, P<0.0001). CONCLUSION: Aiming for an accurate PDFF quantification, high-resolution free-running cardiac CSE-MRI imaging proved to benefit from bipolar echoes with k-space trajectory correction at 3T. This free-breathing acquisition framework enables to investigate epicardial adipose tissue PDFF in metabolic diseases.

Focused navigation for respiratory-motion-corrected free-running radial 4D flow MRI.

PURPOSE: To validate a respiratory motion correction method called focused navigation (fNAV) for free-running radial whole-heart 4D flow MRI. METHODS: Using fNAV, respiratory signals derived from radial readouts are converted into three orthogonal displacements, which are then used to correct respiratory motion in 4D flow datasets. Hundred 4D flow acquisitions were simulated with non-rigid respiratory motion and used for validation. The difference between generated and fNAV displacement coefficients was calculated. Vessel area and flow measurements from 4D flow reconstructions with (fNAV) and without (uncorrected) motion correction were compared to the motion-free ground-truth. In 25 patients, the same measurements were compared between fNAV 4D flow, 2D flow, navigator-gated Cartesian 4D flow, and uncorrected 4D flow datasets. RESULTS: For simulated data, the average difference between generated and fNAV displacement coefficients was 0.04 ± $$ \pm $$ 0.32 mm and 0.31 ± $$ \pm $$ 0.35 mm in the x and y directions, respectively. In the z direction, this difference was region-dependent (0.02 ± $$ \pm $$ 0.51 mm up to 5.85 ± $$ \pm $$ 3.41 mm). For all measurements (vessel area, net volume, and peak flow), the average difference from ground truth was higher for uncorrected 4D flow datasets (0.32 ± $$ \pm $$ 0.11 cm2 , 11.1 ± $$ \pm $$ 3.5 mL, and 22.3 ± $$ \pm $$ 6.0 mL/s) than for fNAV 4D flow datasets (0.10 ± $$ \pm $$ 0.03 cm2 , 2.6 ± $$ \pm $$ 0.7 mL, and 5.1 ± 0 $$ \pm 0 $$ .9 mL/s, p < 0.05). In vivo, average vessel area measurements were 4.92 ± $$ \pm $$ 2.95 cm2 , 5.06 ± $$ \pm $$ 2.64 cm2 , 4.87 ± $$ \pm $$ 2.57 cm2 , 4.87 ± $$ \pm $$ 2.69 cm2 , for 2D flow and fNAV, navigator-gated and uncorrected 4D flow datasets, respectively. In the ascending aorta, all 4D flow datasets except for the fNAV reconstruction had significantly different vessel area measurements from 2D flow. Overall, 2D flow datasets demonstrated the strongest correlation to fNAV 4D flow for both net volume (r2  = 0.92) and peak flow (r2  = 0.94), followed by navigator-gated 4D flow (r2  = 0.83 and r2  = 0.86, respectively), and uncorrected 4D flow (r2  = 0.69 and r2  = 0.86, respectively). CONCLUSION: fNAV corrected respiratory motion in vitro and in vivo, resulting in fNAV 4D flow measurements that are comparable to those derived from 2D flow and navigator-gated Cartesian 4D flow datasets, with improvements over those from uncorrected 4D flow.

Motion-resolved fat-fraction mapping with whole-heart free-running multiecho GRE and pilot tone.

PURPOSE: To develop a free-running 3D radial whole-heart multiecho gradient echo (ME-GRE) framework for cardiac- and respiratory-motion-resolved fat fraction (FF) quantification. METHODS: (NTE = 8) readouts optimized for water-fat separation and quantification were integrated within a continuous non-electrocardiogram-triggered free-breathing 3D radial GRE acquisition. Motion resolution was achieved with pilot tone (PT) navigation, and the extracted cardiac and respiratory signals were compared to those obtained with self-gating (SG). After extra-dimensional golden-angle radial sparse parallel-based image reconstruction, FF, R2 *, and B0 maps, as well as fat and water images were generated with a maximum-likelihood fitting algorithm. The framework was tested in a fat-water phantom and in 10 healthy volunteers at 1.5 T using NTE = 4 and NTE = 8 echoes. The separated images and maps were compared with a standard free-breathing electrocardiogram (ECG)-triggered acquisition. RESULTS: The method was validated in vivo, and physiological motion was resolved over all collected echoes. Across volunteers, PT provided respiratory and cardiac signals in agreement (r = 0.91 and r = 0.72) with SG of the first echo, and a higher correlation to the ECG (0.1% of missed triggers for PT vs. 5.9% for SG). The framework enabled pericardial fat imaging and quantification throughout the cardiac cycle, revealing a decrease in FF at end-systole by 11.4% ± 3.1% across volunteers (p < 0.0001). Motion-resolved end-diastolic 3D FF maps showed good correlation with ECG-triggered measurements (FF bias of -1.06%). A significant difference in free-running FF measured with NTE = 4 and NTE = 8 was found (p < 0.0001 in sub-cutaneous fat and p < 0.01 in pericardial fat). CONCLUSION: Free-running fat fraction mapping was validated at 1.5 T, enabling ME-GRE-based fat quantification with NTE = 8 echoes in 6:15 min.

Pilot tone navigation for respiratory and cardiac motion-resolved free-running 5D flow MRI.

PURPOSE: In this work, we integrated the pilot tone (PT) navigation system into a reconstruction framework for respiratory and cardiac motion-resolved 5D flow. We tested the hypotheses that PT would provide equivalent respiratory curves, cardiac triggers, and corresponding flow measurements to a previously established self-gating (SG) technique while being independent from changes to the acquisition parameters. METHODS: Fifteen volunteers and 9 patients were scanned with a free-running 5D flow sequence, with PT integrated. Respiratory curves and cardiac triggers from PT and SG were compared across all subjects. Flow measurements from 5D flow reconstructions using both PT and SG were compared to each other and to a reference electrocardiogram-gated and respiratory triggered 4D flow acquisition. Radial trajectories with variable readouts per interleave were also tested in 1 subject to compare cardiac trigger quality between PT and SG. RESULTS: The correlation between PT and SG respiratory curves were 0.95 ± 0.06 for volunteers and 0.95 ± 0.04 for patients. Heartbeat duration measurements in volunteers and patients showed a bias to electrocardiogram measurements of, respectively, 0.16 ± 64.94 ms and 0.01 ± 39.29 ms for PT versus electrocardiogram and of 0.24 ± 63.68 ms and 0.09 ± 32.79 ms for SG versus electrocardiogram. No significant differences were reported for the flow measurements between 5D flow PT and from 5D flow SG. A decrease in the cardiac triggering quality of SG was observed for increasing readouts per interleave, whereas PT quality remained constant. CONCLUSION: PT has been successfully integrated in 5D flow MRI and has shown equivalent results to the previously described 5D flow SG technique, while being completely acquisition-independent.

Free-running cardiac and respiratory motion-resolved 5D whole-heart coronary cardiovascular magnetic resonance angiography in pediatric cardiac patients using ferumoxytol.

BACKGROUND: Coronary cardiovascular magnetic resonance angiography (CCMRA) of congenital heart disease (CHD) in pediatric patients requires accurate planning, adequate sequence parameter adjustments, lengthy scanning sessions, and significant involvement from highly trained personnel. Anesthesia and intubation are commonplace to minimize movements and control respiration in younger subjects. To address the above concerns and provide a single-click imaging solution, we applied our free-running framework for fully self-gated (SG) free-breathing 5D whole-heart CCMRA to CHD patients after ferumoxytol injection. We tested the hypothesis that spatial and motion resolution suffice to visualize coronary artery ostia in a cohort of CHD subjects, both for intubated and free-breathing acquisitions. METHODS: In 18 pediatric CHD patients, non-electrocardiogram (ECG) triggered 5D free-running gradient echo CCMRA with whole-heart 1 mm3 isotropic spatial resolution was performed in seven minutes on a 1.5T CMR scanner. Eleven patients were anesthetized and intubated, while seven were breathing freely without anesthesia. All patients were slowly injected with ferumoxytol (4 mg/kg) over 15 minutes. Cardiac and respiratory motion-resolved 5D images were reconstructed with a fully SG approach. To evaluate the performance of motion resolution, visibility of coronary artery origins was assessed. Intubated and free-breathing patient sub-groups were compared for image quality using coronary artery length and conspicuity as well as lung-liver interface sharpness. RESULTS: Data collection using the free-running framework was successful in all patients in less than 8 min; scan planning was very simple without the need for parameter adjustments, while no ECG lead placement and triggering was required. From the resulting SG 5D motion-resolved reconstructed images, coronary artery origins could be retrospectively extracted in 90% of the cases. These general findings applied to both intubated and free-breathing pediatric patients (no difference in terms of lung-liver interface sharpness), while image quality and coronary conspicuity between both cohorts was very similar. CONCLUSIONS: A simple-to-use push-button framework for 5D whole-heart CCMRA was successfully employed in pediatric CHD patients with ferumoxytol injection. This approach, working without any external gating and for a wide range of heart rates and body sizes provided excellent definition of cardiac anatomy for both intubated and free-breathing patients.

Using 5D flow MRI to decode the effects of rhythm on left atrial 3D flow dynamics in patients with atrial fibrillation.

PURPOSE: This study used a 5D flow framework to explore the influence of arrhythmia on thrombogenic hemodynamic parameters in patients with atrial fibrillation (AF). METHODS: A fully self-gated, 3D radial, highly accelerated free-running 5D flow sequence with interleaved four-point velocity-encoding was acquired using an in vitro arrhythmic flow phantom and in 25 patients with a history of AF (68 ± 8 y, 6 female). Self-gating signals were used to calculate AF burden, bin data, and tag each k-space line with its RRLength . Data were binned as an RR-resolved dataset with four RR-interval bins (RR1-RR4, short-to-long) for compressed sensing reconstruction. AF burden was calculated as interquartile range of all intrascan RR-intervals divided by median RR-interval, and left atrial (LA) stasis as the percent of the cardiac cycle where the velocity was <0.1 m/s. RESULTS: In vitro results demonstrated successful recovery of RR-binned flow curves using RR-resolved 5D flow compared to a real-time PC reference standard. In vivo, 5D flow was acquired in 8:48 minutes. AF burden was significantly correlated with 5D flow-derived peak (PV) and mean (MV) velocity and stasis (|ρ| = 0.54-0.75, P < .001). Sensitivity analyses determined a threshold for low versus high AF burden at 9.7%. High burden patients had increased LA mean stasis (up to +42%, P < .01), and lower MV and PV (-30%, -40.6%, respectively, P < .01). RR4 deviated furthest from respiratory-resolved reconstruction (end-expiration) with increased mean stasis (7.6% ± 14.0%, P = .10) and decreased PV (-12.7 ± 14.2%, P = .09). CONCLUSIONS: RR-resolved 5D flow can capture temporal and RR-resolved 3D hemodynamics in <10 minutes and offers a novel approach to investigate arrhythmias.

Similarity-driven multi-dimensional binning algorithm (SIMBA) for free-running motion-suppressed whole-heart MRA.

PURPOSE: Whole-heart MRA techniques typically target predetermined motion states, address cardiac and respiratory dynamics independently, and require either complex planning or computationally demanding reconstructions. In contrast, we developed a fast data-driven reconstruction algorithm with minimal physiological assumptions and compatibility with ungated free-running sequences. THEORY AND METHODS: We propose a similarity-driven multi-dimensional binning algorithm (SIMBA) that clusters continuously acquired k-space data to find a motion-consistent subset for whole-heart MRA reconstruction. Free-running 3D radial data sets from 12 non-contrast-enhanced scans of healthy volunteers and six ferumoxytol-enhanced scans of pediatric cardiac patients were reconstructed with non-motion-suppressed regridding of all the acquired data ("All Data"), with SIMBA, and with a previously published free-running framework (FRF) that uses cardiac and respiratory self-gating and compressed sensing. Images were compared for blood-myocardium sharpness and contrast ratio, visibility of coronary artery ostia, and right coronary artery sharpness. RESULTS: Both the 20-second SIMBA reconstruction and FRF provided significantly higher blood-myocardium sharpness than All Data in both patients and volunteers (P < .05). The SIMBA reconstruction provided significantly sharper blood-myocardium interfaces than FRF in volunteers (P < .001) and higher blood-myocardium contrast ratio than All Data and FRF, both in volunteers and patients (P < .05). Significantly more ostia could be visualized with both SIMBA (31 of 36) and FRF (34 of 36) than with All Data (4 of 36) (P < .001). Inferior right coronary artery sharpness using SIMBA versus FRF was observed (volunteers: SIMBA 36.1 ± 8.1%, FRF 40.4 ± 8.9%; patients: SIMBA 35.9 ± 7.7%, FRF 40.3 ± 6.1%, P = not significant). CONCLUSION: The SIMBA technique enabled a fast, data-driven reconstruction of free-running whole-heart MRA with image quality superior to All Data and similar to the more time-consuming FRF reconstruction.

Volumetric coronary endothelial function assessment: a feasibility study exploiting stack-of-stars 3D cine MRI and image-based respiratory self-gating.

Abnormal coronary endothelial function (CEF), manifesting as depressed vasoreactive responses to endothelial-specific stressors, occurs early in atherosclerosis, independently predicts cardiovascular events, and responds to cardioprotective interventions. CEF is spatially heterogeneous along a coronary artery in patients with atherosclerosis, and thus recently developed and tested non-invasive 2D MRI techniques to measure CEF may not capture the extent of changes in CEF in a given coronary artery. The purpose of this study was to develop and test the first volumetric coronary 3D MRI cine method for assessing CEF along the proximal and mid-coronary arteries with isotropic spatial resolution and in free-breathing. This approach, called 3D-Stars, combines a 6 min continuous, untriggered golden-angle stack-of-stars acquisition with a novel image-based respiratory self-gating method and cardiac and respiratory motion-resolved reconstruction. The proposed respiratory self-gating method agreed well with respiratory bellows and center-of-k-space methods. In healthy subjects, 3D-Stars vessel sharpness was non-significantly different from that by conventional 2D radial in proximal segments, albeit lower in mid-portions. Importantly, 3D-Stars detected normal vasodilatation of the right coronary artery in response to endothelial-dependent isometric handgrip stress in healthy subjects. Coronary artery cross-sectional areas measured using 3D-Stars were similar to those from 2D radial MRI when similar thresholding was used. In conclusion, 3D-Stars offers good image quality and shows feasibility for non-invasively studying vasoreactivity-related lumen area changes along the proximal coronary artery in 3D during free-breathing.

Compressed sensing with signal averaging for improved sensitivity and motion artifact reduction in fluorine-19 MRI.

Fluorine-19 (19 F) MRI of injected perfluorocarbon emulsions (PFCs) allows for the non-invasive quantification of inflammation and cell tracking, but suffers from a low signal-to-noise ratio and extended scan time. To address this limitation, we tested the hypotheses that a 19 F MRI pulse sequence that combines a specific undersampling regime with signal averaging has both increased sensitivity and robustness against motion artifacts compared with a non-averaged fully sampled pulse sequence, when both datasets are reconstructed with compressed sensing. As a proof of principle, numerical simulations and phantom experiments were performed on selected variable ranges to characterize the point spread function of undersampling patterns, as well as the vulnerability to noise of undersampling and reconstruction parameters with paired numbers of x signal averages and acceleration factor x (NAx-AFx). The numerical simulations demonstrated that a probability density function that uses 25% of the samples to fully sample the k-space central area allowed for an optimal balance between limited blurring and artifact incoherence. At all investigated noise levels, the Dice similarity coefficient (DSC) strongly depended on the regularization parameters and acceleration factor. In phantoms, the motion robustness of an NA8-AF8 undersampling pattern versus NA1-AF1 was evaluated with simulated and real motion patterns. Differences were assessed with the DSC, which was consistently higher for the NA8-AF8 compared with the NA1-AF1 strategy, for both simulated and real cyclic motion patterns (P < 0.001). Both strategies were validated in vivo in mice (n = 2) injected with perfluoropolyether. Here, the images displayed a sharper delineation of the liver with the NA8-AF8 strategy than with the NA1-AF1 strategy. In conclusion, we validated the hypotheses that in 19 F MRI the combination of undersampling and averaging improves both the sensitivity and the robustness against motion artifacts.

Reduced-iodine-dose dual-energy coronary CT angiography: qualitative and quantitative comparison between virtual monochromatic and polychromatic CT images.

OBJECTIVES: To quantitatively evaluate the impact of virtual monochromatic images (VMI) on reduced-iodine-dose dual-energy coronary computed tomography angiography (CCTA) in terms of coronary lumen segmentation in vitro, and secondly to assess the image quality in vivo, compared with conventional CT obtained with regular iodine dose. MATERIALS AND METHODS: A phantom simulating regular and reduced iodine injection was used to determine the accuracy and precision of lumen area segmentation for various VMI energy levels. We retrospectively included 203 patients from December 2017 to August 2018 (mean age, 51.7 ± 16.8 years) who underwent CCTA using either standard (group A, n = 103) or reduced (group B, n = 100) iodine doses. Conventional images (group A) were qualitatively and quantitatively compared with 55-keV VMI (group B). We recorded the location of venous catheters. RESULTS: In vitro, VMI outperformed conventional CT, with a segmentation accuracy of 0.998 vs. 1.684 mm2, respectively (p < 0.001), and a precision of 0.982 vs. 1.229 mm2, respectively (p < 0.001), in simulated overweight adult subjects. In vivo, the rate of diagnostic CCTA in groups A and B was 88.4% (n = 91/103) vs. 89% (n = 89/100), respectively, and noninferiority of protocol B was inferred. Contrast-to-noise ratios (CNR) of lumen versus fat and muscle were higher in group B (p < 0.001) and comparable for lumen versus calcium (p = 0.423). Venous catheters were more often placed on the forearm or hand in group B (p < 0.001). CONCLUSION: In vitro, low-keV VMI improve vessel area segmentation. In vivo, low-keV VMI allows for a 40% iodine dose and injection rate reduction while maintaining diagnostic image quality and improves the CNR between lumen versus fat and muscle. KEY POINTS: • Dual-energy coronary CT angiography is becoming increasingly available and might help improve patient management. • Compared with regular-iodine-dose coronary CT angiography, reduced-iodine-dose dual-energy CT with low-keV monochromatic image reconstructions performed better in phantom-based vessel cross-sectional segmentation and proved to be noninferior in vivo. • Patients receiving reduced-iodine-dose dual-energy coronary CT angiography often had the venous catheter placed on the forearm or wrist without compromising image quality.

Endogenous assessment of myocardial injury with single-shot model-based non-rigid motion-corrected T1 rho mapping.

BACKGROUND: Cardiovascular magnetic resonance T1ρ mapping may detect myocardial injuries without exogenous contrast agent. However, multiple co-registered acquisitions are required, and the lack of robust motion correction limits its clinical translation. We introduce a single breath-hold myocardial T1ρ mapping method that includes model-based non-rigid motion correction. METHODS: A single-shot electrocardiogram (ECG)-triggered balanced steady state free precession (bSSFP) 2D adiabatic T1ρ mapping sequence that collects five T1ρ-weighted (T1ρw) images with different spin lock times within a single breath-hold is proposed. To address the problem of residual respiratory motion, a unified optimization framework consisting of a joint T1ρ fitting and model-based non-rigid motion correction algorithm, insensitive to contrast change, was implemented inline for fast (~ 30 s) and direct visualization of T1ρ maps. The proposed reconstruction was optimized on an ex vivo human heart placed on a motion-controlled platform. The technique was then tested in 8 healthy subjects and validated in 30 patients with suspected myocardial injury on a 1.5T CMR scanner. The Dice similarity coefficient (DSC) and maximum perpendicular distance (MPD) were used to quantify motion and evaluate motion correction. The quality of T1ρ maps was scored. In patients, T1ρ mapping was compared to cine imaging, T2 mapping and conventional post-contrast 2D late gadolinium enhancement (LGE). T1ρ values were assessed in remote and injured areas, using LGE as reference. RESULTS: Despite breath holds, respiratory motion throughout T1ρw images was much larger in patients than in healthy subjects (5.1 ± 2.7 mm vs. 0.5 ± 0.4 mm, P < 0.01). In patients, the model-based non-rigid motion correction improved the alignment of T1ρw images, with higher DSC (87.7 ± 5.3% vs. 82.2 ± 7.5%, P < 0.01), and lower MPD (3.5 ± 1.9 mm vs. 5.1 ± 2.7 mm, P < 0.01). This resulted in significantly improved quality of the T1ρ maps (3.6 ± 0.6 vs. 2.1 ± 0.9, P < 0.01). Using this approach, T1ρ mapping could be used to identify LGE in patients with 93% sensitivity and 89% specificity. T1ρ values in injured (LGE positive) areas were significantly higher than in the remote myocardium (68.4 ± 7.9 ms vs. 48.8 ± 6.5 ms, P < 0.01). CONCLUSIONS: The proposed motion-corrected T1ρ mapping framework enables a quantitative characterization of myocardial injuries with relatively low sensitivity to respiratory motion. This technique may be a robust and contrast-free adjunct to LGE for gaining new insight into myocardial structural disorders.

Measurement accuracy of prototype non-contrast, compressed sensing-based, respiratory motion-resolved whole heart cardiovascular magnetic resonance angiography for the assessment of thoracic aortic dilatation: comparison with computed tomography angiography.

BACKGROUND: Patients with thoracic aortic dilatation who undergo annual computed tomography angiography (CTA) are subject to repeated radiation and contrast exposure. The purpose of this study was to evaluate the feasibility of a non-contrast, respiratory motion-resolved whole-heart cardiovascular magnetic resonance angiography (CMRA) technique against reference standard CTA, for the quantitative assessment of cardiovascular anatomy and monitoring of disease progression in patients with thoracic aortic dilatation.  METHODS: Twenty-four patients (68.6 ± 9.8 years) with thoracic aortic dilatation prospectively underwent clinical CTA and research 1.5T CMRA between July 2017 and November 2018. Scans were repeated in 15 patients 1 year later. A prototype free-breathing 3D radial balanced steady-state free-precession whole-heart CMRA sequence was used in combination with compressed sensing-based reconstruction. Area, circumference, and diameter measurements were obtained at seven aortic levels by two experienced and two inexperienced readers. In addition, area and diameter measurements of the cardiac chambers, pulmonary arteries and pulmonary veins were also obtained. Agreement between the two modalities was assessed with intraclass correlation coefficient (ICC) analysis, Bland-Altman plots and scatter plots. RESULTS: Area, circumference and diameter measurements on a per-level analysis showed good or excellent agreement between CTA and CMRA (ICCs > 0.84). Means of differences on Bland-Altman plots were: area 0.0 cm2 [- 1.7; 1.6]; circumference 1.0 mm [- 10.0; 12.0], and diameter 0.6 mm [- 2.6; 3.6]. Area and diameter measurements of the left cardiac chambers showed good agreement (ICCs > 0.80), while moderate to good agreement was observed for the right chambers (all ICCs > 0.56). Similar good to excellent inter-modality agreement was shown for the pulmonary arteries and veins (ICC range 0.79-0.93), with the exception of the left lower pulmonary vein (ICC < 0.51). Inter-reader assessment demonstrated mostly good or excellent agreement for both CTA and CMRA measurements on a per-level analysis (ICCs > 0.64). Difference in maximum aortic diameter measurements at baseline vs follow up showed excellent agreement between CMRA and CTA (ICC = 0.91). CONCLUSIONS: The radial whole-heart CMRA technique combined with respiratory motion-resolved reconstruction provides comparable anatomical measurements of the thoracic aorta and cardiac structures as the reference standard CTA. It could potentially be used to diagnose and monitor patients with thoracic aortic dilatation without exposing them to radiation or contrast media.

Free-running 5D coronary MR angiography at 1.5T using LIBRE water excitation pulses.

PURPOSE: To implement, optimize, and characterize lipid-insensitive binomial off-resonant RF excitation (LIBRE) pulses for fat-suppressed fully self-gated free-running 5D cardiac MRI. METHODS: Bloch equation simulations were used to optimize LIBRE parameter settings in non-interrupted bSSFP prior to in vitro validation. Thus, optimized LIBRE pulses were subsequently applied to free-running coronary MRA in 20 human adult subjects, where resulting images were quantitatively compared to those obtained with non-fat-suppressing excitation (SP), conventional 1-2-1 water excitation (WE), and a previously published interrupted free-running (IFR) sequence. SAR and scan times were recorded. Respiratory-and-cardiac-motion-resolved images were reconstructed with XD-GRASP, and contrast ratios, coronary artery detection rate, vessel length, and vessel sharpness were computed. RESULTS: The numerically optimized LIBRE parameters were successfully validated in vitro. In vivo, LIBRE had the lowest SAR and a scan time that was similar to that of WE yet 18% shorter than that of IFR. LIBRE improved blood-fat contrast when compared to SP, WE, and IFR, vessel detection relative to SP and IFR, and vessel sharpness when compared to WE and IFR (for example, for the left main and anterior descending coronary artery, 51.5% ± 10.2% [LIBRE] versus 42.1% ± 6.8% [IFR]). Vessel length measurements remained unchanged for all investigated methods. CONCLUSION: LIBRE enabled fully self-gated non-interrupted free-running 5D bSSFP imaging of the heart at 1.5T with suppressed fat signal. Measures of image quality, vessel conspicuity, and scan time compared favorably to those obtained with the more conventional non-interrupted WE and the previously published IFR, while SAR reduction offers added flexibility.

Natively fat-suppressed 5D whole-heart MRI with a radial free-running fast-interrupted steady-state (FISS) sequence at 1.5T and 3T.

PURPOSE: To implement, optimize, and test fast interrupted steady-state (FISS) for natively fat-suppressed free-running 5D whole-heart MRI at 1.5 tesla (T) and 3T. METHODS: FISS was implemented for fully self-gated free-running cardiac- and respiratory-motion-resolved radial imaging of the heart at 1.5T and 3T. Numerical simulations and phantom scans were performed to compare fat suppression characteristics and to determine parameter ranges (number of readouts [NR] per FISS module and TR) for effective fat suppression. Subsequently, free-running FISS data were collected in 10 healthy volunteers and images were reconstructed with compressed sensing. All acquisitions were compared with a continuous balanced steady-state free precession version of the same sequence, and both fat suppression and scan times were analyzed. RESULTS: Simulations demonstrate a variable width and location of suppression bands in FISS that were dependent on TR and NR. For a fat suppression bandwidth of 100 Hz and NR ≤ 8, simulations demonstrated that a TR between 2.2 ms and 3.0 ms is required at 1.5T, whereas a range of 3.0 ms to 3.5 ms applies at 3T. Fat signal increases with NR. These findings were corroborated in phantom experiments. In volunteers, fat SNR was significantly decreased using FISS compared with balanced steady-state free precession (P < 0.05) at both field strengths. After protocol optimization, high-resolution (1.1 mm3 ) 5D whole-heart free-running FISS can be performed with effective fat suppression in under 8 min at 1.5T and 3T at a modest scan time increase compared to balanced steady-state free precession. CONCLUSION: An optimal FISS parameter range was determined enabling natively fat-suppressed 5D whole-heart free-running MRI with a single continuous scan at 1.5T and 3T, demonstrating potential for cardiac imaging and noncontrast angiography.

An automated approach to fully self-gated free-running cardiac and respiratory motion-resolved 5D whole-heart MRI.

PURPOSE: To develop a previously reported, electrocardiogram (ECG)-gated, motion-resolved 5D compressed sensing whole-heart sparse MRI methodology into an automated, optimized, and fully self-gated free-running framework in which external gating or triggering devices are no longer needed. METHODS: Cardiac and respiratory self-gating signals were extracted from raw image data acquired in 12 healthy adult volunteers with a non-ECG-triggered 3D radial golden-angle 1.5 T balanced SSFP sequence. To extract cardiac self-gating signals, central k-space coefficient signal analysis (k0 modulation), as well as independent and principal component analyses were performed on selected k-space profiles. The procedure yielding triggers with the smallest deviation from those of the reference ECG was selected for the automated protocol. Thus, optimized cardiac and respiratory self-gating signals were used for binning in a compressed sensing reconstruction pipeline. Coronary vessel length and sharpness of the resultant 5D images were compared with image reconstructions obtained with ECG-gating. RESULTS: Principal component analysis-derived cardiac self-gating triggers yielded a smaller deviation ( 17.4±6.1ms ) from the reference ECG counterparts than k0 modulation ( 26±7.5ms ) or independent component analysis ( 19.8±5.2ms ). Cardiac and respiratory motion-resolved 5D images were successfully reconstructed with the automated and fully self-gated approach. No significant difference was found for coronary vessel length and sharpness between images reconstructed with the fully self-gated and the ECG-gated approach (all P≥.06 ). CONCLUSION: Motion-resolved 5D compressed sensing whole-heart sparse MRI has successfully been developed into an automated, optimized, and fully self-gated free-running framework in which external gating, triggering devices, or navigators are no longer mandatory. The resultant coronary MRA image quality was equivalent to that obtained with conventional ECG-gating.

Precision and accuracy of cross-sectional area measurements used to measure coronary endothelial function with spiral MRI.

PURPOSE: Coronary endothelial function (CEF) reflects vascular health and conventional invasive CEF measures predict cardiovascular events. MRI can now noninvasively measure CEF by quantifying coronary artery cross-sectional area changes in response to isometric handgrip exercise, an endothelial-dependent stressor. Area changes (10 to 20% in healthy; 2 to -12% in impaired vessels) are only a few imaging voxels because of MRI's limited spatial resolution. Here, with numerical simulations and phantom studies, we test whether Fourier interpolation enables sub-pixel area measurement precision and determine the smallest detectable area change using spiral MRI. METHODS: In vivo coronary SNR with the currently used CEF protocol at 3T was measured in 7 subjects for subsequent in vitro work. Area measurements of circular vessels were simulated by varying partial volume, vessel diameter, voxel size, SNR, and Fourier interpolation factor. A phantom with precision-drilled holes (diameters 3-3.42 mm) was imaged 10 times with the current CEF protocol (voxel size, Δx = 0.89 mm) and a high-resolution protocol (Δx = 0.6 mm) to determine precision, accuracy, and the smallest detectable area changes. RESULTS: In vivo coronary SNR ranged from 30-76. Eight-fold Fourier interpolation improved area measurement precision by a factor 6.5 and 4.9 in the simulations and phantom scans, respectively. The current CEF protocol can detect mean area changes of 4-5% for SNR above 30, and 3-3.5% for SNR above 40 with a higher-resolution protocol. CONCLUSION: Current CEF spiral MRI with in vivo SNR allows detection of a 4-5% area change and Fourier interpolation improves precision several-fold to sub-voxel dimensions.

Accelerated and high-resolution cardiac T2 mapping through peripheral k-space sharing.

PURPOSE: To develop high-spatial-resolution cardiac T2 mapping that allows for a reduced acquisition time while maintaining its precision. We implemented and optimized a new golden-angle radial T2 mapping technique named SKRATCH (Shared k-space Radial T2 Characterization of the Heart) that shares k-space peripheries of T2 -weighted images while preserving their contrasts. METHODS: Six SKRATCH variants (gradient-recalled echo and balanced SSFP, free-breathing and breath-held, with and without a saturation preparation) were implemented, and their precision was compared with a navigator-gated reference technique in phantoms and 22 healthy volunteers at 3 T. The optimal breath-held SKRATCH technique was applied in a small cohort of patients with subacute myocardial infarction. RESULTS: The faster free-breathing SKRATCH technique reduced the acquisition time by 52.4%, while maintaining the precision and spatial resolution of the reference technique. Similarly, the most precise and robust breath-held SKRATCH technique demonstrated homogenous T2 values that did not significantly differ from the navigator-gated reference (T2 = 39.9 ± 3.4 ms versus 39.5 ± 3.4 ms, P > .20, respectively). All infarct patients demonstrated a large T2 elevation in the ischemic regions of the myocardium. CONCLUSION: The optimized SKRATCH technique enabled the accelerated acquisition of high-spatial-resolution T2 maps, was validated in healthy adult volunteers, and was successfully applied to a small initial group of patients.

A black-blood ultra-short echo time (UTE) sequence for 3D isotropic resolution imaging of the lungs.

PURPOSE: Ultra-short echo time MRI is a promising alternative to chest CT for cystic fibrosis patients. Black-blood imaging in particular could help discern small-sized anomalies, such as mucoid plugging, which may otherwise be confused with neighboring blood vessels, particularly when contrast agent is not used. We, therefore, implemented and tested an ultra-short echo time sequence with black-blood preparation. Additionally, this sequence may also be used to generate bright-blood angiograms. METHODS: Using this sequence, data was acquired during free breathing in 10 healthy volunteers to obtain respiratory-motion-resolved 3D volumes covering the entire thorax with an isotropic resolution of (1 mm)3 . The magnitude of signal suppression relative to a bright-blood reference acquisition was quantified and compared with that obtained with a turbo-spin echo (TSE) acquisition. Bright-blood angiograms were also generated by subtraction. Finally, an initial feasibility assessment was performed in 2 cystic fibrosis patients, and images were visually compared with contrast-enhanced images and with CT data. RESULTS: Black-blood preparation significantly decreased the average normalized signal intensity in the vessel lumen (-66%; P < 0.001). Similarly, blood signal was significantly lowered (-60%; P = 0.001) compared with the TSE acquisition. In patients, mucoid plugging could be emphasized in the black-blood datasets. An intercostal artery could also be visualized in the subtraction angiograms. CONCLUSION: Black-blood free-breathing ultra-short echo time imaging was successfully implemented and motion-resolved full volumetric coverage of the lungs with high spatial resolution was achieved, while obtaining an angiogram without contrast agent injection. Encouraging initial results in patients prompt further investigations in a larger cohort.