Detecting liver fibrosis using a machine learning-based approach to the quantification of the heart-induced deformation in tagged MR images.

Liver disease causes millions of deaths per year worldwide, and approximately half of these cases are due to cirrhosis, which is an advanced stage of liver fibrosis that can be accompanied by liver failure and portal hypertension. Early detection of liver fibrosis helps in improving its treatment and prevents its progression to cirrhosis. In this work, we present a novel noninvasive method to detect liver fibrosis from tagged MRI images using a machine learning-based approach. Specifically, coronal and sagittal tagged MRI imaging are analyzed separately to capture cardiac-induced deformation of the liver. The liver is manually delineated and a novel image feature, namely, the histogram of the peak strain (HPS) value, is computed from the segmented liver region and is used to classify the liver as being either normal or fibrotic. Classification is achieved using a support vector machine algorithm. The in vivo study included 15 healthy volunteers (10 males; age range 30-45 years) and 22 patients (15 males; age range 25-50 years) with liver fibrosis verified and graded by transient elastography, and 10 patients only had a liver biopsy and were diagnosed with a score of F3-F4. The proposed method demonstrates the usefulness and efficiency of extracting the HPS features from the sagittal slices for patients with moderate fibrosis. Cross-validation of the method showed an accuracy of 83.7% (specificity = 86.6%, sensitivity = 81.8%).

Imaging sequence for joint myocardial T1 mapping and fat/water separation.

PURPOSE: To develop and evaluate an imaging sequence to simultaneously quantify the epicardial fat volume and myocardial T1 relaxation time. METHODS: We introduced a novel simultaneous myocardial T1 mapping and fat/water separation sequence (joint T1 -fat/water separation). Dixon reconstruction is performed on a dual-echo data set to generate water/fat images. T1 maps are computed using the water images, whereas the epicardial fat volume is calculated from the fat images. A phantom experiment using vials with different T1 /T2 values and a bottle of oil was performed. Additional phantom experiment using vials of mixed fat/water was performed to show the potential of this sequence to mitigate the effect of intravoxel fat on estimated T1 maps. In vivo evaluation was performed in 17 subjects. Epicardial fat volume, native myocardial T1 measurements and precision were compared among slice-interleaved T1 mapping, Dixon, and the proposed sequence. RESULTS: In the first phantom, the proposed sequence separated oil from water vials and there were no differences in T1 of the fat-free vials (P = .1). In the second phantom, the T1 error decreased from 22%, 36%, 57%, and 73% to 8%, 9%, 16%, and 26%, respectively. In vivo there was no difference between myocardial T1 values (1067 ± 17 ms versus 1077 ± 24 ms, P = .6). The epicardial fat volume was similar for both sequences (54.3 ± 33 cm3 versus 52.4 ± 32 cm3 , P = .8). CONCLUSION: The proposed sequence provides simultaneous quantification of native myocardial T1 and epicardial fat volume. This will eliminate the need for an additional sequence in the cardiac imaging protocol if both measurements are clinically indicated.

Improved dark blood late gadolinium enhancement (DB-LGE) imaging using an optimized joint inversion preparation and T2 magnetization preparation.

PURPOSE: To develop a dark blood-late gadolinium enhancement (DB-LGE) sequence that improves scar-blood contrast and delineation of scar region. METHODS: The DB-LGE sequence uses an inversion pulse followed by T2 magnetization preparation to suppress blood and normal myocardium. Time delays inserted after preparation pulses and T2 -magnetization-prep duration are used to adjust tissue contrast. Selection of these parameters was optimized using numerical simulations and phantom experiments. We evaluated DB-LGE in 9 swine and 42 patients (56 ± 14 years, 33 male). Improvement in scar-blood contrast and overall image quality was subjectively evaluated by two independent readers (1 = poor, 4 = excellent). The signal ratios among scar, blood, and myocardium were compared. RESULTS: Simulations and phantom studies demonstrated that simultaneous nulling of myocardium and blood can be achieved by selecting appropriate timing parameters. The scar-blood contrast score was significantly higher for DB-LGE (P < 0.001) with no significant difference in overall image quality (P > 0.05). Scar-blood signal ratios for DB-LGE versus LGE were 5.0 ± 2.8 versus 1.5 ± 0.5 (P < 0.001) for patients, and 2.2 ± 0.7 versus 1.0 ± 0.4 (P = 0.0023) for animals. Scar-myocardium signal ratios were 5.7 ± 2.9 versus 6.3 ± 2.6 (P = 0.35) for patients, and 3.7 ± 1.1 versus 4.1 ± 2.0 (P = 0.60) for swine. CONCLUSIONS: The DB-LGE sequence simultaneously reduces normal myocardium and blood signal intensity, thereby enhancing scar-blood contrast while preserving scar-myocardium contrast. Magn Reson Med 79:351-360, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Gray blood late gadolinium enhancement cardiovascular magnetic resonance for improved detection of myocardial scar.

BACKGROUND: Low scar-to-blood contrast in late gadolinium enhanced (LGE) MRI limits the visualization of scars adjacent to the blood pool. Nulling the blood signal improves scar detection but results in lack of contrast between myocardium and blood, which makes clinical evaluation of LGE images more difficult. METHODS: GB-LGE contrast is achieved through partial suppression of the blood signal using T2 magnetization preparation between the inversion pulse and acquisition. The timing parameters of GB-LGE sequence are determined by optimizing a cost-function representing the desired tissue contrast. The proposed 3D GB-LGE sequence was evaluated using phantoms, human subjects (n = 45) and a swine model of myocardial infarction (n = 5). Two independent readers subjectively evaluated the image quality and ability to identify and localize scarring in GB-LGE compared to black-blood LGE (BB-LGE) (i.e., with complete blood nulling) and conventional (bright-blood) LGE. RESULTS: GB-LGE contrast was successfully generated in phantoms and all in-vivo scans. The scar-to-blood contrast was improved in GB-LGE compared to conventional LGE in humans (1.1 ± 0.5 vs. 0.6 ± 0.4, P < 0.001) and in animals (1.5 ± 0.2 vs. -0.03 ± 0.2). In patients, GB-LGE detected more tissue scarring compared to BB-LGE and conventional LGE. The subjective scores of the GB-LGE ability for localizing LV scar and detecting papillary scar were improved as compared with both BB-LGE (P < 0.024) and conventional LGE (P < 0.001). In the swine infarction model, GB-LGE scores for the ability to localize LV scar scores were consistently higher than those of both BB-LGE and conventional-LGE. CONCLUSION: GB-LGE imaging improves the ability to identify and localize myocardial scarring compared to both BB-LGE and conventional LGE. Further studies are warranted to histologically validate GB-LGE.

Clinical performance of high-resolution late gadolinium enhancement imaging with compressed sensing.

PURPOSE: To evaluate diagnostic image quality of 3D late gadolinium enhancement (LGE) with high isotropic spatial resolution (∼1.4 mm3 ) images reconstructed from randomly undersampled k-space using LOw-dimensional-structure Self-learning and Thresholding (LOST). MATERIALS AND METHODS: We prospectively enrolled 270 patients (181 men; 55 ± 14 years) referred for myocardial viability assessment. 3D LGE with isotropic spatial resolution of 1.4 ± 0.1 mm3 was acquired at 1.5T using a LOST acceleration rate of 3 to 5. In a subset of 121 patients, 3D LGE or phase-sensitive LGE were acquired with parallel imaging with an acceleration rate of 2 for comparison. Two readers evaluated image quality using a scale of 1 (poor) to 4 (excellent) and assessed for scar presence. The McNemar test statistic was used to compare the proportion of detected scar between the two sequences. We assessed the association between image quality and characteristics (age, gender, torso dimension, weight, heart rate), using generalized linear models. RESULTS: Overall, LGE detection proportions for 3D LGE with LOST were similar between readers 1 and 2 (16.30% vs. 18.15%). For image quality, readers gave 85.9% and 80.0%, respectively, for images categorized as good or excellent. Overall proportion of scar presence was not statistically different from conventional 3D LGE (28% vs. 33% [P = 0.17] for reader 1 and 26% vs. 31% [P = 0.37] for reader 2). Increasing subject heart rate was associated with lower image quality (estimated slope = -0.009 (P = 0.001)). CONCLUSION: High-resolution 3D LGE with LOST yields good to excellent image quality in >80% of patients and identifies patients with LV scar at the same rate as conventional 3D LGE. LEVEL OF EVIDENCE: 2 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2017;46:1829-1838.

Joint myocardial T1 and T2 mapping using a combination of saturation recovery and T2 -preparation.

PURPOSE: To develop a heart-rate independent breath-held joint T1 -T2 mapping sequence for accurate simultaneous estimation of coregistered myocardial T1 and T2 maps. METHODS: A novel preparation scheme combining both a saturation pulse and T2 -preparation in a single R-R interval is introduced. The time between these two pulses, as well as the duration of the T2 -preparation is varied in each heartbeat, acquiring images with different T1 and T2 weightings, and no magnetization dependence on previous images. Inherently coregistered T1 and T2 maps are calculated from these images. Phantom imaging is performed to compare the proposed maps with spin echo references. In vivo imaging is performed in ten subjects, comparing the accuracy and precision of the proposed technique to existing myocardial T1 and T2 mapping sequences of the same duration. RESULTS: Phantom experiments show that the proposed technique provides accurate quantification of T1 and T2 values over a wide-range (T1 : 260 ms to 1460 ms, T2 : 40 ms to 200 ms). In vivo imaging shows that the proposed sequence quantifies T1 and T2 values similar to a saturation-based T1 mapping and a conventional breath-hold T2 mapping sequence, respectively. CONCLUSION: The proposed sequence allows joint estimation of accurate and coregistered quantitative myocardial T1 and T2 maps in a single breath-hold. Magn Reson Med 76:888-896, 2016. © 2015 Wiley Periodicals, Inc.

Free-breathing slice-interleaved myocardial T2 mapping with slice-selective T2 magnetization preparation.

PURPOSE: To develop and evaluate a free-breathing slice-interleaved T2 mapping sequence by proposing a new slice-selective T2 magnetization preparation (T2 prep) sequence that allows interleaved data acquisition for different slices in subsequent heartbeats. METHODS: We developed a slice-selective T2 prep for myocardial T2 mapping by adding slice-selective gradients to a conventional single-slice T2 prep sequence. In this sequence, five slices are acquired during five consecutive heartbeats, each using a slice-selective T2 prep. The scheme was repeated four times using different T2 prep echo times. We compared the performance of the proposed slice-interleaved T2 mapping sequence and the conventional single-slice T2 mapping sequence in term of accuracy, precision, and reproducibility using phantom experiments and in vivo imaging in 10 healthy subjects. We also evaluated the feasibility of the proposed sequence in 28 patients with cardiovascular disease, and the quality of the maps was scored subjectively. Furthermore, we investigated the impact of through-plane motion by comparing T2 measurements acquired during end-systole versus mid-diastole. RESULTS: T2 measurements using a slice-interleaved T2 mapping sequence were correlated with a spin echo (r(2) = 0.88) and single-slice T2 mapping sequence (r(2) = 0.98). The mean myocardial T2 values were correlated between slice-interleaved (48 ms) and single-slice (51 ms) T2 mapping sequences. Subjective scores of T2 map quality were good to excellent in 81% of the maps in patients. There was no difference in T2 measurements between end-systole versus mid-diastole. CONCLUSIONS: The proposed free-breathing slice-interleaved T2 mapping sequence allows T2 measurements of five left ventricular slices in 20 heartbeats with similar reproducibility and precision as the single-slice T2 mapping sequence but with a four-fold reduction in acquisition time. Magn Reson Med 76:555-565, 2016. © 2015 Wiley Periodicals, Inc.

Comparison of spoiled gradient echo and steady-state free-precession imaging for native myocardial T1 mapping using the slice-interleaved T1 mapping (STONE) sequence.

Cardiac T1 mapping allows non-invasive imaging of interstitial diffuse fibrosis. Myocardial T1 is commonly calculated by voxel-wise fitting of the images acquired using balanced steady-state free precession (SSFP) after an inversion pulse. However, SSFP imaging is sensitive to B1 and B0 imperfection, which may result in additional artifacts. A gradient echo (GRE) imaging sequence has been used for myocardial T1 mapping; however, its use has been limited to higher magnetic field to compensate for the lower signal-to-noise ratio (SNR) of GRE versus SSFP imaging. A slice-interleaved T1 mapping (STONE) sequence with SSFP readout (STONE-SSFP) has been recently proposed for native myocardial T1 mapping, which allows longer recovery of magnetization (>8 R-R) after each inversion pulse. In this study, we hypothesize that a longer recovery allows higher SNR and enables native myocardial T1 mapping using STONE with GRE imaging readout (STONE-GRE) at 1.5T. Numerical simulations and phantom and in vivo imaging were performed to compare the performance of STONE-GRE and STONE-SSFP for native myocardial T1 mapping at 1.5T. In numerical simulations, STONE-SSFP shows sensitivity to both T2 and off resonance. Despite the insensitivity of GRE imaging to T2 , STONE-GRE remains sensitive to T2 due to the dependence of the inversion pulse performance on T2 . In the phantom study, STONE-GRE had inferior accuracy and precision and similar repeatability as compared with STONE-SSFP. In in vivo studies, STONE-GRE and STONE-SSFP had similar myocardial native T1 times, precisions, repeatabilities and subjective T1 map qualities. Despite the lower SNR of the GRE imaging readout compared with SSFP, STONE-GRE provides similar native myocardial T1 measurements, precision, repeatability, and subjective image quality when compared with STONE-SSFP at 1.5T.

Reproducibility of myocardial T1 and T2 relaxation time measurement using slice-interleaved T1 and T2 mapping sequences.

PURPOSE: To assess measurement reproducibility and image quality of myocardial T1 and T2 maps using free-breathing slice-interleaved T1 and T2 mapping sequences at 1.5 Tesla (T). MATERIALS AND METHODS: Eleven healthy subjects (33 ± 16 years; 6 males) underwent a slice-interleaved T1 and T2 mapping test/retest cardiac MR study at 1.5T on 2 days. For each day, subjects were imaged in two sessions with removal out of the magnet and repositioning before the subsequent session. We studied measurement reproducibility as well as the required sample size for sufficient statistical power to detect a predefined change in T1 and T2 . In a separate prospective study, we assessed T1 and T2 map image quality in 241 patients (54 ± 15 years; 73 women) with known/suspected cardiovascular disease referred for clinical cardiac MR. A subjective quality score was used to assess a segment-based image quality. RESULTS: In the healthy cohort, the slice-interleaved T1 measurements were highly reproducible, with global coefficients of variation (CVs) of 2.4% between subjects, 2.1% between days, and 1.7% between sessions. Slice-interleaved T2 mapping sequences provided similar reproducibility with global CVs of 7.2% between subjects, 6.3% between days, and 5.0 between sessions. A lower variability resulted in a reduction of the required number of subjects to achieve a certain statistical power when compared with other T1 mapping sequences. In the subjective image quality assessment, >80% of myocardial segments had interpretable data. CONCLUSION: Slice-interleaved T1 and T2 mapping sequences yield highly reproducible T1 and T2 measurements with >80% of interpretable myocardial segments. J. Magn. Reson. Imaging 2016;44:1159-1167.

Relationship between native papillary muscle T1 time and severity of functional mitral regurgitation in patients with non-ischemic dilated cardiomyopathy.

BACKGROUND: Functional mitral regurgitation is one of the severe complications of non-ischemic dilated cardiomyopathy (DCM). Non-contrast native T1 mapping has emerged as a non-invasive method to evaluate myocardial fibrosis. We sought to evaluate the potential relationship between papillary muscle T1 time and mitral regurgitation in DCM patients. METHODS: Forty DCM patients (55 ± 13 years) and 20 healthy adult control subjects (54 ± 13 years) were studied. Native T1 mapping was performed using a slice interleaved T1 mapping sequence (STONE) which enables acquisition of 5 slices in the short-axis plane within a 90 s free-breathing scan. We measured papillary muscle diameter, length and shortening. DCM patients were allocated into 2 groups based on the presence or absence of functional mitral regurgitation. RESULTS: Papillary muscle T1 time was significantly elevated in DCM patients with mitral regurgitation (n = 22) in comparison to those without mitral regurgitation (n = 18) (anterior papillary muscle: 1127 ± 36 msec vs 1063 ± 16 msec, p < 0.05; posterior papillary muscle: 1124 ± 30 msec vs 1062 ± 19 msec, p < 0.05), but LV T1 time was similar (1129 ± 38 msec vs 1134 ± 58 msec, p = 0.93). Multivariate linear regression analysis showed that papillary muscle native T1 time (ß = 0.10, 95 % CI: 0.05-0.17, p < 0.05) is significantly correlated with mitral regurgitant fraction. Elevated papillary muscle T1 time was associated with larger diameter, longer length and decreased papillary muscle shortening (all p values <0.05). CONCLUSIONS: In DCM, papillary muscle native T1 time is significantly elevated and related to mitral regurgitant fraction.

Improved quantitative myocardial T2 mapping: Impact of the fitting model.

PURPOSE: To develop an improved T2 prepared (T2 prep) balanced steady-state free-precession (bSSFP) sequence and signal relaxation curve fitting method for myocardial T2 mapping. METHODS: Myocardial T2 mapping is commonly performed by acquisition of multiple T2 prep bSSFP images and estimating the voxel-wise T2 values using a two-parameter fit for relaxation. However, a two-parameter fit model does not take into account the effect of imaging pulses in a bSSFP sequence or other imperfections in T2 prep RF pulses, which may decrease the robustness of T2 mapping. Therefore, we propose a novel T2 mapping sequence that incorporates an additional image acquired with saturation preparation, simulating a very long T2 prep echo time. This enables the robust estimation of T2 maps using a 3-parameter fit model, which captures the effect of imaging pulses and other imperfections. Phantom imaging is performed to compare the T2 maps generated using the proposed 3-parameter model with the conventional two-parameter model, as well as a spin echo reference. In vivo imaging is performed on eight healthy subjects to compare the different fitting models. RESULTS: Phantom and in vivo data show that the T2 values generated by the proposed 3-parameter model fitting do not change with different choices of the T2 prep echo times, and are not statistically different than the reference values for the phantom (P = 0.10 with three T2 prep echoes). The two-parameter model exhibits dependence on the choice of T2 prep echo times and are significantly different than the reference values (P = 0.01 with three T2 prep echoes). CONCLUSION: The proposed imaging sequence in combination with a three-parameter model allows accurate measurement of myocardial T2 values, which is independent of number and duration of T2 prep echo times. Magn Reson Med 74:93-105, 2015. © 2014 Wiley Periodicals, Inc.

Free-breathing multislice native myocardial T1 mapping using the slice-interleaved T1 (STONE) sequence.

PURPOSE: To develop a novel pulse sequence for free-breathing, multislice, native myocardial T1 mapping. METHODS: The slice-interleaved T1 (STONE) sequence consists of multiple sets of single-shot images of different slices, acquired after a single nonselective inversion pulse. Each slice is only selectively excited once after each inversion pulse to allow sampling of the unperturbed longitudinal magnetization in the adjacent slices. For respiratory motion, a prospective slice-tracking respiratory navigator is used to decrease through-plane motion followed by a retrospective image registration to reduce in-plane motion. STONE T1 maps were calculated using both a two-parameter and three-parameter fit model. The accuracy and precision of the STONE sequence for different T1 , T2 , and inversion pulse efficiency were studied using numerical simulations and phantom experiments. T1 maps from 14 subjects were acquired with the STONE sequence and T1 s were compared to the MOdified Look-Locker Inversion recovery sequence (MOLLI). RESULTS: In numerical simulations and phantom experiments, the STONE sequence using a two-parameter fit model yields more accurate T1 times compared to MOLLI, with similar high precision. The three-parameter fit model further improves the accuracy, but with a reduced precision. The native myocardial T1 times were higher in the STONE sequence using two- or three-parameter fit compared to MOLLI. The standard deviation of the T1 times was lower in the STONE T1 maps with a two-parameter fit compared with MOLLI or a three-parameter fit. CONCLUSION: The STONE sequence allows accurate and precise quantification of native myocardial T1 times with the additional benefit of covering the entire ventricle. Magn Reson Med 74:115-124, 2015. © 2014 Wiley Periodicals, Inc.

Free-breathing combined three-dimensional phase sensitive late gadolinium enhancement and T1 mapping for myocardial tissue characterization.

PURPOSE: To develop a novel MR sequence for combined three-dimensional (3D) phase-sensitive (PS) late gadolinium enhancement (LGE) and T1 mapping to allow for simultaneous assessment of focal and diffuse myocardial fibrosis. METHODS: In the proposed sequence, four 3D imaging volumes are acquired with different T1 weightings using a combined saturation and inversion preparation, after administration of a gadolinium contrast agent. One image is acquired fully sampled with the inversion time selected to null the healthy myocardial signal (the LGE image). The other three images are three-fold under-sampled and reconstructed using compressed sensing. An acquisition scheme with two interleaved imaging cycles and joint navigator-gating of those cycles ensures spatial registration of the imaging volumes. T1 maps are generated using all four imaging volumes. The signal-polarity in the LGE image is restored using supplementary information from the T1 fit to generate PS-LGE images. The accuracy of the proposed method was assessed with respect to a inversion-recovery spin-echo sequence. In vivo T1 maps and LGE images were acquired with the proposed sequence and quantitatively compared with 2D multislice Modified Look-Locker inversion recovery (MOLLI) T1 maps. Exemplary images in a patient with focal scar were compared with conventional LGE imaging. RESULTS: The deviation of the proposed method and the spin-echo reference was < 11 ms in phantom for T1 times between 250 and 600 ms, regardless of the inversion time selected in the LGE image. There was no significant difference in the in vivo T1 times of the proposed sequence and the 2D MOLLI technique (myocardium: 292 ± 75 ms versus 310 ± 49 ms, blood-pools: 191 ± 75 ms versus 182.0 ± 33). The LGE images showed proper nulling of the healthy myocardium in all subjects and clear depiction of scar in the patient. CONCLUSION: The proposed sequence enables simultaneous acquisition of 3D PS-LGE images and spatially registered 3D T1 maps in a single scan.

Free-breathing post-contrast three-dimensional T1 mapping: Volumetric assessment of myocardial T1 values.

PURPOSE: To develop a three-dimensional (3D) free-breathing myocardial T1 mapping sequence for assessment of left ventricle diffuse fibrosis after contrast administration. METHODS: In the proposed sequence, multiple 3D inversion recovery images are acquired in an interleaved manner. A mixed prospective/retrospective navigator scheme is used to obtain the 3D Cartesian k-space data with fully sampled center and randomly undersampled outer k-space. The resulting undersampled 3D k-space data are then reconstructed using compressed sensing. Subsequently, T1 maps are generated by voxel-wise curve fitting of the individual interleaved images. In a phantom study, the accuracy of the 3D sequence was evaluated against two-dimensional (2D) modified Look-Locker inversion recovery (MOLLI) and spin-echo sequences. In vivo T1 times of the proposed method were compared with 2D multislice MOLLI T1 mapping. Subsequently, the feasibility of high-resolution 3D T1 mapping with spatial resolution of 1.7 × 1.7 × 4 mm(3) was demonstrated. RESULTS: The proposed method shows good agreement with 2D MOLLI and the spin-echo reference in phantom. No significant difference was found in the in vivo T1 times estimated using the proposed sequence and the 2D MOLLI technique (myocardium, 330 ± 66 ms versus 319 ± 93 ms; blood pools, 211 ± 68 ms versus 210 ± 98 ms). However, improved homogeneity, as measured using standard deviation of the T1 signal, was observed in the 3D T1 maps. CONCLUSION: The proposed sequence enables high-resolution 3D T1 mapping after contrast injection during free-breathing with volumetric left ventricle coverage.

Accelerated three-dimensional cine phase contrast imaging using randomly undersampled echo planar imaging with compressed sensing reconstruction.

The aim of this study was to implement and evaluate an accelerated three-dimensional (3D) cine phase contrast MRI sequence by combining a randomly sampled 3D k-space acquisition sequence with an echo planar imaging (EPI) readout. An accelerated 3D cine phase contrast MRI sequence was implemented by combining EPI readout with randomly undersampled 3D k-space data suitable for compressed sensing (CS) reconstruction. The undersampled data were then reconstructed using low-dimensional structural self-learning and thresholding (LOST). 3D phase contrast MRI was acquired in 11 healthy adults using an overall acceleration of 7 (EPI factor of 3 and CS rate of 3). For comparison, a single two-dimensional (2D) cine phase contrast scan was also performed with sensitivity encoding (SENSE) rate 2 and approximately at the level of the pulmonary artery bifurcation. The stroke volume and mean velocity in both the ascending and descending aorta were measured and compared between two sequences using Bland-Altman plots. An average scan time of 3 min and 30 s, corresponding to an acceleration rate of 7, was achieved for 3D cine phase contrast scan with one direction flow encoding, voxel size of 2 × 2 × 3 mm(3) , foot-head coverage of 6 cm and temporal resolution of 30 ms. The mean velocity and stroke volume in both the ascending and descending aorta were statistically equivalent between the proposed 3D sequence and the standard 2D cine phase contrast sequence. The combination of EPI with a randomly undersampled 3D k-space sampling sequence using LOST reconstruction allows a seven-fold reduction in scan time of 3D cine phase contrast MRI without compromising blood flow quantification.

Impact of motion correction on reproducibility and spatial variability of quantitative myocardial T2 mapping.

BACKGROUND: To evaluate and quantify the impact of a novel image-based motion correction technique in myocardial T2 mapping in terms of measurement reproducibility and spatial variability. METHODS: Twelve healthy adult subjects were imaged using breath-hold (BH), free breathing (FB), and free breathing with respiratory navigator gating (FB + NAV) myocardial T2 mapping sequences. Fifty patients referred for clinical CMR were imaged using the FB + NAV sequence. All sequences used a T2 prepared (T2prep) steady-state free precession acquisition. In-plane myocardial motion was corrected using an adaptive registration of varying contrast-weighted images for improved tissue characterization (ARCTIC). DICE similarity coefficient (DSC) and myocardial boundary errors (MBE) were measured to quantify the motion estimation accuracy in healthy subjects. T2 mapping reproducibility and spatial variability were evaluated in healthy subjects using 5 repetitions of the FB + NAV sequence with either 4 or 20 T2prep echo times (TE). Subjective T2 map quality was assessed in patients by an experienced reader using a 4-point scale (1-non diagnostic, 4-excellent). RESULTS: ARCTIC led to increased DSC in BH data (0.85 ± 0.08 vs. 0.90 ± 0.02, p = 0.007), FB data (0.78 ± 0.13 vs. 0.90 ± 0.21, p < 0.001), and FB + NAV data (0.86 ± 0.05 vs. 0.90 ± 0.02, p = 0.002), and reduced MBE in BH data (0.90 ± 0.40 vs. 0.64 ± 0.19 mm, p = 0.005), FB data (1.21 ± 0.65 vs. 0.63 ± 0.10 mm, p < 0.001), and FB + NAV data (0.81 ± 0.21 vs. 0.63 ± 0.08 mm, p < 0.001). Improved reproducibility (4TE: 5.3 ± 2.5 ms vs. 4.0 ± 1.5 ms, p = 0.016; 20TE: 3.9 ± 2.3 ms vs. 2.2 ± 0.5 ms, p = 0.002), reduced spatial variability (4TE: 12.8 ± 3.5 ms vs. 10.3 ± 2.5 ms, p < 0.001; 20TE: 9.7 ± 3.5 ms vs. 7.5 ± 1.4 ms) and improved subjective score of T2 map quality (3.43 ± 0.79 vs. 3.69 ± 0.55, p < 0.001) were obtained using ARCTIC. CONCLUSIONS: The ARCTIC technique substantially reduces spatial mis-alignment among T2-weighted images and improves the reproducibility and spatial variability of in-vivo T2 mapping.

Accelerated isotropic sub-millimeter whole-heart coronary MRI: compressed sensing versus parallel imaging.

PURPOSE: To enable accelerated isotropic sub-millimeter whole-heart coronary MRI within a 6-min acquisition and to compare this with a current state-of-the-art accelerated imaging technique at acceleration rates beyond what is used clinically. METHODS: Coronary MRI still faces major challenges, including lengthy acquisition time, low signal-to-noise-ratio (SNR), and suboptimal spatial resolution. Higher spatial resolution in the sub-millimeter range is desirable, but this results in increased acquisition time and lower SNR, hindering its clinical implementation. In this study, we sought to use an advanced B1-weighted compressed sensing technique for highly accelerated sub-millimeter whole-heart coronary MRI, and to compare the results to parallel imaging, the current-state-of-the-art, where both techniques were used at acceleration rates beyond what is used clinically. Two whole-heart coronary MRI datasets were acquired in seven healthy adult subjects (30.3 ± 12.1 years; 3 men), using prospective 6-fold acceleration, with random undersampling for the proposed compressed sensing technique and with uniform undersampling for sensitivity encoding reconstruction. Reconstructed images were qualitatively compared in terms of image scores and perceived SNR on a four-point scale (1 = poor, 4 = excellent) by an experienced blinded reader. RESULTS: The proposed technique resulted in images with clear visualization of all coronary branches. Overall image quality and perceived SNR of the compressed sensing images were significantly higher than those of parallel imaging (P = 0.03 for both), which suffered from noise amplification artifacts due to the reduced SNR. CONCLUSION: The proposed compressed sensing-based reconstruction and acquisition technique for sub-millimeter whole-heart coronary MRI provides 6-fold acceleration, where it outperforms parallel imaging with uniform undersampling.

Localized spatio-temporal constraints for accelerated CMR perfusion.

PURPOSE: To develop and evaluate an image reconstruction technique for cardiac MRI (CMR) perfusion that uses localized spatio-temporal constraints. METHODS: CMR perfusion plays an important role in detecting myocardial ischemia in patients with coronary artery disease. Breath-hold k-t-based image acceleration techniques are typically used in CMR perfusion for superior spatial/temporal resolution and improved coverage. In this study, we propose a novel compressed sensing-based image reconstruction technique for CMR perfusion, with applicability to free-breathing examinations. This technique uses local spatio-temporal constraints by regularizing image patches across a small number of dynamics. The technique was compared with conventional dynamic-by-dynamic reconstruction, and sparsity regularization using a temporal principal-component (pc) basis, as well as zero-filled data in multislice two-dimensional (2D) and three-dimensional (3D) CMR perfusion. Qualitative image scores were used (1 = poor, 4 = excellent) to evaluate the technique in 3D perfusion in 10 patients and five healthy subjects. On four healthy subjects, the proposed technique was also compared with a breath-hold multislice 2D acquisition with parallel imaging in terms of signal intensity curves. RESULTS: The proposed technique produced images that were superior in terms of spatial and temporal blurring compared with the other techniques, even in free-breathing datasets. The image scores indicated a significant improvement compared with other techniques in 3D perfusion (x-pc regularization, 2.8 ± 0.5 versus 2.3 ± 0.5; dynamic-by-dynamic, 1.7 ± 0.5; zero-filled, 1.1 ± 0.2). Signal intensity curves indicate similar dynamics of uptake between the proposed method with 3D acquisition and the breath-hold multislice 2D acquisition with parallel imaging. CONCLUSION: The proposed reconstruction uses sparsity regularization based on localized information in both spatial and temporal domains for highly accelerated CMR perfusion with potential use in free-breathing 3D acquisitions.

Free-breathing cardiac MR stress perfusion with real-time slice tracking.

PURPOSE: To develop a free-breathing cardiac MR perfusion sequence with slice tracking for use after physical exercise. METHODS: We propose to use a leading navigator, placed immediately before each 2D slice acquisition, for tracking the respiratory motion and updating the slice location in real-time. The proposed sequence was used to acquire CMR perfusion datasets in 12 healthy adult subjects and 8 patients. Images were compared with the conventional perfusion (i.e., without slice tracking) results from the same subjects. The location and geometry of the myocardium were quantitatively analyzed, and the perfusion signal curves were calculated from both sequences to show the efficacy of the proposed sequence. RESULTS: The proposed sequence was significantly better compared with the conventional perfusion sequence in terms of qualitative image scores. Changes in the myocardial location and geometry decreased by 50% in the slice tracking sequence. Furthermore, the proposed sequence had signal curves that are smoother and less noisy. CONCLUSION: The proposed sequence significantly reduces the effect of the respiratory motion on the image acquisition in both rest and stress perfusion scans.

Free-breathing phase contrast MRI with near 100% respiratory navigator efficiency using k-space-dependent respiratory gating.

PURPOSE: To investigate the efficacy of a novel respiratory motion scheme, where only the center of k-space is gated using respiratory navigators, versus a fully respiratory-gated acquisition for three-dimensional flow imaging. METHODS: Three-dimensional flow images were acquired axially using a gradient echo sequence in a volume, covering the ascending and descending aorta, and the pulmonary artery bifurcation in 12 healthy subjects (33.2 ± 15.8 years; five men). For respiratory motion compensation, two gating and tracking strategies were used with a 7-mm gating window: (1) All of k-space acquired within the gating window (fully gated) and (2) central k-space acquired within the gating window, and the remainder of k-space acquired without any gating (center gated). Each scan was repeated twice. Stroke volume, mean flow, peak velocity, and signal-to-noise-ratio measurements were performed both on the ascending and on the descending aorta for all acquisitions, which were compared using a linear mixed-effects model and Bland-Altman analysis. RESULTS: There were no statistical differences between the fully gated and the center-gated strategies for the quantification of stroke volume, peak velocity, and mean flow, as well as the signal-to-noise-ratio measurements. Furthermore, the proposed center-gated strategy had significantly shorter acquisition time compared to the fully gated strategy (13:19 ± 3:02 vs. 19:35 ± 5:02, P < 0.001). CONCLUSIONS: The proposed novel center-gated strategy for three-dimensional flow MRI allows for markedly shorter acquisition time without any systematic variation in quantitative flow measurements in this small group of healthy volunteers.