Ristretto MRE: A generalized multi‐shot GRE‐MRE sequence

In order to acquire consistent k‐space data in MR elastography, a fixed temporal relationship between the MRI sequence and the underlying period of the wave needs to be ensured. To this end, conventional GRE‐MRE enforces synchronization through repeated triggering of the transducer and forcing the sequence repetition time to be equal to an integer multiple of the wave period. For wave frequencies below 100 Hz, however, this leads to prolonged acquisition times, as the repetition time scales inversely with frequency. A previously developed multi‐shot approach (eXpresso MRE) to multi‐slice GRE‐MRE tackles this issue by acquiring an integer number of slices per wave period, which allows acquisition to be accelerated in typical scenarios by a factor of two or three. In this work, it is demonstrated that the constraints imposed by the eXpresso scheme are overly restrictive. We propose a generalization of the sequence in three steps by incorporating sequence delays into imaging shots and allowing for interleaved wave‐phase acquisition. The Ristretto scheme is compared in terms of imaging shot and total scan duration relative to eXpresso and conventional GRE‐MRE and is validated in three different phantom studies. First, the agreement of measured displacement fields in different stages of the sequence generalization is shown. Second, performance is compared for 25, 36, 40, and 60 Hz actuation frequencies. Third, the performance is assessed for the acquisition of different numbers of slices (13 to 17). In vivo feasibility is demonstrated in the liver and the breast. Here, Ristretto is compared with an optimized eXpresso sequence, leading to scan accelerations of 15% and 5%, respectively, without compromising displacement field and stiffness estimates in general. The Ristretto concept allows us to choose imaging shot durations on a fine grid independent of the number of slices and the wave frequency, permitting 2‐ to 4.5‐fold acceleration of conventional GRE‐MRE acquisitions.     

Measuring changes in muscle stiffness after eccentric exercise using elastography

Muscle stiffness has been reported to increase following eccentric muscle exercise, but to date only indirect methods have been used to measure it. This study aimed to use Magnetic Resonance Elastography (MRE), a noninvasive imaging technique, to assess the time‐course of passive elasticity changes in the medial gastrocnemius and soleus muscles before and after a bout of eccentric exercise. Shear storage modulus (G′) and loss modulus (G′′) measurements were made in eight healthy subjects for both muscles in vivo before, one hour after, 48 hours after and 1 week after eccentric exercise. The results show a 21% increase in medial gastrocnemius storage modulus following eccentric exercise with a peak occurring ~48 hours after exercise (before exercise 1.15 ± 0.23 kPa, 48 hours after 1.38 ± 0.27 kPa). No significant changes in soleus muscle storage modulus were measured for the exercise protocol used in this study, and no significant changes in loss modulus were observed. This study provides the first direct measurements in skeletal muscle before and after eccentric exercise damage and suggests that MRE can be used to detect the time course of changes to muscle properties. Copyright © 2012 John Wiley & Sons, Ltd.     

Vibration safety limits for magnetic resonance elastography

Magnetic resonance elastography (MRE) has been demonstrated to have potential as a clinical tool for assessing the stiffness of tissue in vivo. An essential step in MRE is the generation of acoustic mechanical waves within a tissue via a coupled mechanical driver. Motivated by an increasing volume of human imaging trials using MRE, the objectives of this study were to audit the vibration amplitude of exposure for our IRB-approved human MRE studies, to compare these values to a conservative regulatory standard for vibrational exposure and to evaluate the applicability and implications of this standard for MRE. MRE displacement data were examined from 29 MRE exams, including the liver, brain, kidney, breast and skeletal muscle. Vibrational acceleration limits from a European Union directive limiting occupational exposure to whole-body and extremity vibrations (EU 2002/44/EC) were adjusted for time and frequency of exposure, converted to maximum displacement values and compared to the measured in vivo displacements. The results indicate that the vibrational amplitudes used in MRE studies are below the EU whole-body vibration limit, and the EU guidelines represent a useful standard that could be readily accepted by Institutional Review Boards to define standards for vibrational exposures for MRE studies in humans.     

Rapid acquisition of multifrequency,multislice and multidirectional MR elastography data with a fractionally encoded gradient echo sequence

In MR elastography (MRE), periodic tissue motion is phase encoded using motion‐encoding gradients synchronized to an externally applied periodic mechanical excitation. Conventional methods result in extended scan time for quality phase images, thus limiting the broad application of MRE in the clinic. For practical scan times, researchers have been relying on one‐dimensional or two‐dimensional motion‐encoding, low‐phase sampling and a limited number of slices, and artifact‐prone, single‐shot, echo planar imaging (EPI) readout. Here, we introduce a rapid multislice pulse sequence capable of three‐dimensional motion encoding that is also suitable for simultaneously encoding motion with multiple frequency components. This sequence is based on a gradient‐recalled echo (GRE) sequence and exploits the principles of fractional encoding. This GRE MRE pulse sequence was validated as capable of acquiring full three‐dimensional motion encoding of isotropic voxels in a large volume within less than a minute. This sequence is suitable for monofrequency and multifrequency MRE experiments. In homogeneous paraffin phantoms, the eXpresso sequence yielded similar storage modulus values as those obtained with conventional methods, although with markedly reduced variances (7.11 ± 0.26 kPa for GRE MRE versus 7.16 ± 1.33 kPa for the conventional spin‐echo EPI sequence). The GRE MRE sequence obtained better phase‐to‐noise ratios than the equivalent spin‐echo EPI sequence (matched for identical acquisition time) in both paraffin phantoms and in vivo data in the liver (59.62 ± 11.89 versus 27.86 ± 3.81, 61.49 ± 14.16 versus 24.78 ± 2.48 and 58.23 ± 10.39 versus 23.48 ± 2.91 in the X, Y and Z components, respectively, in the case of liver experiments). Phase‐to‐noise ratios were similar between GRE MRE used in monofrequency or multifrequency experiments (75.39 ± 14.93 versus 86.13 ± 18.25 at 28 Hz, 71.52 ± 24.74 versus 86.96 ± 30.53 at 56 Hz and 95.60 ± 36.96 versus 61.35 ± 26.25 at 84Hz, respectively). Copyright © 2013 John Wiley & Sons, Ltd.     

Retrospective motion gating in cardiac MRI using a simultaneously acquired navigator

A simultaneous acquisition technique of image and navigator signals (simultaneously acquired navigator, SIMNAV) is proposed for cardiac magnetic resonance imaging (CMRI) in Cartesian coordinates. To simultaneously acquire both image and navigator signals, a conventional balanced steady‐state free precession (bSSFP) pulse sequence is modified by adding a radiofrequency (RF) pulse, which excites a supplementary slice for the navigator signal. Alternating phases of the RF pulses make it easy to separate the simultaneously acquired magnetic resonance data into image and navigator signals. The navigator signals of the proposed SIMNAV were compared with those of current gating devices and self‐gating techniques for seven healthy subjects. In vivo experiments demonstrated that SIMNAV could provide cardiac cine images with sufficient image quality, similar to those from electrocardiogram (ECG) gating with breath‐hold. SIMNAV can be used to acquire a cardiac cine image without requiring an ECG device and breath‐hold, whilst maintaining feasible imaging time efficiency.     

Viscoelastic properties of soft gels: comparison of magnetic resonance elastography and dynamic shear testing in the shear wave regime

Magnetic resonance elastography (MRE) is used to quantify the viscoelastic shear modulus, G*, of human and animal tissues. Previously, values of G* determined by MRE have been compared to values from mechanical tests performed at lower frequencies. In this study, a novel dynamic shear test (DST) was used to measure G* of a tissue-mimicking material at higher frequencies for direct comparison to MRE. A closed-form solution, including inertial effects, was used to extract G* values from DST data obtained between 20 and 200 Hz. MRE was performed using cylindrical 'phantoms' of the same material in an overlapping frequency range of 100-400 Hz. Axial vibrations of a central rod caused radially propagating shear waves in the phantom. Displacement fields were fit to a viscoelastic form of Navier's equation using a total least-squares approach to obtain local estimates of G*. DST estimates of the storage G' (Re[G*]) and loss modulus G″ (Im[G*]) for the tissue-mimicking material increased with frequency from 0.86 to 0.97 kPa (20-200 Hz, n = 16), while MRE estimates of G' increased from 1.06 to 1.15 kPa (100-400 Hz, n = 6). The loss factor (Im[G*]/Re[G*]) also increased with frequency for both test methods: 0.06-0.14 (20-200 Hz, DST) and 0.11-0.23 (100-400 Hz, MRE). Close agreement between MRE and DST results at overlapping frequencies indicates that G* can be locally estimated with MRE over a wide frequency range. Low signal-to-noise ratio, long shear wavelengths and boundary effects were found to increase residual fitting error, reinforcing the use of an error metric to assess confidence in local parameter estimates obtained by MRE.     

Magnetic resonance elastography of the lungs: A repeatability and reproducibility study

Lung diseases are one of the leading causes of death worldwide, from which four million people die annually. Lung diseases are associated with changes in the mechanical properties of the lungs. Several studies have shown the feasibility of using magnetic resonance elastography (MRE) to quantify the lungs' shear stiffness. The aim of this study is to investigate the reproducibility and repeatability of lung MRE, and its shear stiffness measurements, obtained using a modified spin echo‐echo planar imaging (SE‐EPI) MRE sequence. In this study, 21 healthy volunteers were scanned twice by repositioning the volunteers to image right lung both at residual volume (RV) and total lung capacity (TLC) to assess the reproducibility of lung shear stiffness measurements. Additionally, 19 out of the 21 volunteers were scanned immediately without moving the volunteers to test the repeatability of the modified SE‐EPI MRE sequence. A paired t‐test was performed to determine the significant difference between stiffness measurements obtained at RV and TLC. Concordance correlation and Bland–Altman's analysis were performed to determine the reproducibility and repeatability of the SE‐EPI MRE‐derived shear stiffness measurements. The SE‐EPI MRE sequence is highly repeatable with a concordance correlation coefficient (CCC) of 0.95 at RV and 0.96 at TLC. Similarly, the stiffness measurements obtained across all volunteers were highly reproducible with a CCC of 0.95 at RV and 0.92 at TLC. The mean shear stiffness of the lung at RV was 0.93 ± 0.22 kPa and at TLC was 1.41 ± 0.41 kPa. TLC showed a significantly higher mean shear stiffness (P = 0.0004) compared with RV. Lung MRE stiffness measurements obtained using the SE‐EPI sequence were reproducible and repeatable, both at RV and TLC. Lung shear stiffness changes across respiratory cycle with significantly higher stiffness at TLC than RV.     

Simultaneous multislice cardiac magnetic resonance fingerprinting using low rank reconstruction

This study introduces a technique for simultaneous multislice (SMS) cardiac magnetic resonance fingerprinting (cMRF), which improves the slice coverage when quantifying myocardial T_1, T₂, and M₀. The single‐slice cMRF pulse sequence was modified to use multiband (MB) RF pulses for SMS imaging. Different RF phase schedules were used to excite each slice, similar to POMP or CAIPIRINHA, which imparts tissues with a distinguishable and slice‐specific magnetization evolution over time. Because of the high net acceleration factor (R = 48 in plane combined with the slice acceleration), images were first reconstructed with a low rank technique before matching data to a dictionary of signal timecourses generated by a Bloch equation simulation. The proposed method was tested in simulations with a numerical relaxation phantom. Phantom and in vivo cardiac scans of 10 healthy volunteers were also performed at 3 T. With single‐slice acquisitions, the mean relaxation times obtained using the low rank cMRF reconstruction agree with reference values. The low rank method improves the precision in T₁ and T₂ for both single‐slice and SMS cMRF, and it enables the acquisition of maps with fewer artifacts when using SMS cMRF at higher MB factors. With this technique, in vivo cardiac maps were acquired from three slices simultaneously during a breathhold lasting 16 heartbeats. SMS cMRF improves the efficiency and slice coverage of myocardial T₁ and T₂ mapping compared with both single‐slice cMRF and conventional cardiac mapping sequences. Thus, this technique is a first step toward whole‐heart simultaneous T₁ and T₂ quantification with cMRF.     

Overhauser‐enhanced magnetic resonance elastography

Magnetic resonance elastography (MRE) is a powerful technique to assess the mechanical properties of living tissue. However, it suffers from reduced sensitivity in regions with short T₂ and T₂* such as in tissue with high concentrations of paramagnetic iron, or in regions surrounding implanted devices. In this work, we exploit the longer T₂* attainable at ultra‐low magnetic fields in combination with Overhauser dynamic nuclear polarization (DNP) to enable rapid MRE at 0.0065 T. A 3D balanced steady‐state free precession based MRE sequence with undersampling and fractional encoding was implemented on a 0.0065 T MRI scanner. A custom‐built RF coil for DNP and a programmable vibration system for elastography were developed. Displacement fields and stiffness maps were reconstructed from data recorded in a polyvinyl alcohol gel phantom loaded with stable nitroxide radicals. A DNP enhancement of 25 was achieved during the MRE sequence, allowing the acquisition of 3D Overhauser‐enhanced MRE (OMRE) images with (1.5 × 2.7 × 9) mm³ resolution over eight temporal steps and 11 slices in 6 minutes. In conclusion, OMRE at ultra‐low magnetic field can be used to detect mechanical waves over short acquisition times. This new modality shows promise to broaden the scope of conventional MRE applications, and may extend the utility of low‐cost, portable MRI systems to detect elasticity changes in patients with implanted devices or iron overload. Copyright © 2016 John Wiley & Sons, Ltd.     

Non-invasive measurement of brain viscoelasticity using magnetic resonance elastography

The purpose of this work was to develop magnetic resonance elastography (MRE) for the fast and reproducible measurement of spatially averaged viscoelastic constants of living human brain. The technique was based on a phase-sensitive echo planar imaging acquisition. Motion encoding was orthogonal to the image plane and synchronized to intracranial shear vibrations at driving frequencies of 25 and 50 Hz induced by a head-rocker actuator. Ten time-resolved phase-difference wave images were recorded within 60 s and analyzed for shear stiffness and shear viscosity. Six healthy volunteers (six men; mean age 34.5 years; age range 25-44 years) underwent 23-39 follow-up MRE studies over a period of 6 months. Interindividual mean +/- SD shear moduli and shear viscosities were found to be 1.17 +/- 0.03 kPa and 3.1 +/- 0.4 Pas for 25 Hz and 1.56 +/- 0.07 kPa and 3.4 +/- 0.2 Pas for 50 Hz, respectively (P < or = 0.01). The intraindividual range of shear modulus data was 1.01-1.31 kPa (25 Hz) and 1.33-1.77 kPa (50 Hz). The observed modulus dispersion indicates a limited applicability of Voigt's model to explain viscoelastic behavior of brain parenchyma within the applied frequency range. The narrow distribution of data within small confidence intervals demonstrates excellent reproducibility of the experimental protocol. The results are necessary as reference data for future comparisons between healthy and pathological human brain viscoelastic data.     

Analysis of errors in diffusion kurtosis imaging caused by slice crosstalk in simultaneous multi‐slice imaging

Simultaneous multi‐slice (SMS) imaging techniques accelerate diffusion MRI data acquisition. However, slice separation is imperfect and results in residual signal leakage between the simultaneously excited slices. The resulting consistent bias may adversely affect diffusion model parameter estimation. Although this bias is usually small and might not affect the simplified diffusion tensor model significantly, higher order diffusion models such as kurtosis are likely to be more susceptible to such effects. In this work, two SMS reconstruction techniques and an alternative acquisition approach were tested to quantify the effects of slice crosstalk on diffusion kurtosis parameters. In reconstruction, two popular slice separation algorithms, slice GRAPPA and split‐slice GRAPPA, are evaluated to determine the effect of slice leakage on diffusion kurtosis metrics. For the alternative acquisition, the slice pairings were varied across diffusion weighted images such that the signal leakage does not come from the same overlapped slice for all diffusion encodings. Simulation results demonstrated the potential benefits of randomizing the slice pairings. However, various experimental factors confounded the advantages of slice pair randomization. In volunteer experiments, region‐of‐interest analyses found high metric errors with each of the SMS acquisitions and reconstructions in the brain white matter.     

Simultaneous multi‐slice (SMS) acquisition enhances the sensitivity of hemodynamic mapping using gas challenges

Hemodynamic mapping using gas inhalation has received increasing interest in recent years. Cerebrovascular reactivity (CVR), which reflects the ability of the brain vasculature to dilate in response to a vasoactive stimulus, can be measured by CO₂ inhalation with continuous acquisition of blood oxygen level‐dependent (BOLD) magnetic resonance images. Cerebral blood volume (CBV) can be measured by O₂ inhalation. These hemodynamic mapping methods are appealing because of their absence of gadolinium contrast agent, their ability to assess both baseline perfusion and vascular reserve, and their utility in calibrating the functional magnetic resonance imaging (fMRI) signal. However, like other functional and physiological indices, a major drawback of these measurements is their poor sensitivity and reliability. Simultaneous multi‐slice echo planar imaging (SMS EPI) is a fast imaging technology that allows the excitation and acquisition of multiple two‐dimensional slices simultaneously, and has been shown to enhance the sensitivity of several MRI applications. To our knowledge, the benefit of SMS in gas inhalation imaging has not been investigated. In this work, we compared the sensitivity of CO₂ and O₂ inhalation data collected using SMS factor 2 (SMS2) and SMS factor 3 (SMS3) with those collected using conventional EPI (SMS1). We showed that the sensitivity of SMS scans was significantly (p = 0.01) higher than that of conventional EPI, although no difference was found between SMS2 and SMS3 (p = 0.3). On a voxel‐wise level, approximately 20–30% of voxels in the brain showed a significant enhancement in sensitivity when using SMS compared with conventional EPI, with other voxels showing an increase, but not reaching statistical significance. When using SMS, the scan duration can be reduced by half, whilst maintaining the sensitivity of conventional EPI. The availability of a sensitive acquisition technique can further enhance the potential of gas inhalation MRI in clinical applications.     

An optimization framework to maximize signal‐to‐noise ratio in simultaneous multi‐slice body imaging

Parallel imaging is essential for the acceleration of abdominal and pelvic 2D multi‐slice imaging, in order to reduce scan time and mitigate motion artifacts. Controlled Aliasing In Parallel Imaging Results IN Higher Acceleration (CAIPIRINHA) accelerated imaging has been shown to increase the signal‐to‐noise ratio (SNR) significantly compared with in‐plane parallel imaging with similar acceleration. We hypothesize that for CAIPIRINHA‐accelerated abdominal imaging the consistency of image quality and SNR is more difficult to achieve due to the subject‐specific coil sensitivity profiles, caused by (1) flexible coil placement; (2) variations in anatomy; and (3) variations in scan coverage along the superior–inferior direction. To test this, a mathematical framework is introduced that calculates the (retained) SNR for in‐plane and simultaneous multi‐slice (SMS)‐accelerated acquisitions. Moreover, this framework was used to optimize the sampling pattern by maximizing the local SNR within a region of interest (ROI) through non‐linear, RF‐induced CAIPIRINHA slice shifts. The framework was evaluated on 14 healthy subjects and the optimized sampling pattern was compared with in‐plane acceleration and CAIPIRINHA acceleration with linear slice shifts, which are primarily used in brain imaging. We demonstrate that the field of view (FOV) in the superior–inferior direction, the coil positioning and the individual anatomy have a large impact on the image SNR (changes up to 50% for varying coil positions and 40% differences between subjects) and image artifacts for simultaneous multi‐slice acceleration. Consequently, sampling patterns have to be optimized for acquisitions employing different FOVs and ideally on an individual basis. Optimization of the sampling pattern, which exploits non‐linear shifts between slices, showed a considerable SNR increase (10–30%) for higher acceleration factors. The framework outlined in this article can be used to optimize sampling patterns for a broad range of accelerated body acquisitions on an individual basis. Copyright © 2015 John Wiley & Sons, Ltd.     

Magnetic Resonance Elastography of kidneys: SE‐EPI MRE reproducibility and its comparison to GRE MRE

The purpose of this study is 1) to demonstrate reproducibility of spin echo‐echo planar imaging (SE‐EPI) magnetic resonance elastography (MRE) to estimate kidney stiffness; and 2) to compare SE‐EPI MRE and gradient recalled echo (GRE) MRE‐derived stiffness estimations in various anatomical regions of the kidney. Kidney MRE was performed on 33 healthy subjects (8 for SE‐EPI MRE reproducibility and 25 for comparison with GRE MRE; age range: 22–66 years) in a 3 T MRI scanner. To demonstrate SE‐EPI MRE reproducibility, subjects were scanned for the first scan and then asked to leave the scan room and repositioned again for the second (repeat) scan. Similar set‐up was used for GRE MRE as well. The displacement data was then processed to obtain overall stiffness estimates of the kidney. Concordance correlation analyses were performed to determine SE‐EPI MRE reproducibility and agreement between GRE MRE and SE‐EPI MRE derived stiffness. A high concordance correlation (ρ_c = 0.95; p‐value<0.0001) was obtained for SE‐EPI MRE reproducibility. Good concordance correlation was observed (ρ_c = 0.84; p < 0.0001 for both kidneys, ρ_c = 0.91; p < 0.0001 for right kidney and ρ_c = 0.78; p < 0.0001 for left kidney) between GRE MRE and SE‐EPI MRE derived stiffness measurements. Paired t‐test results showed that stiffness value of medulla was significantly (p < 0.0001) greater than cortex using SE‐EPI MRE as well as GRE MRE. SE‐EPI MRE was reproducible and good agreement was observed in MRE‐derived stiffness measurements obtained using SE‐EPI and GRE sequences. Therefore, SE‐EPI can be used for kidney MRE applications.     

Analysis and improvement of motion encoding in magnetic resonance elastography

Magnetic resonance elastography (MRE) utilizes phase contrast magnetic resonance imaging (MRI), which is phase locked to externally generated mechanical vibrations, to measure the three‐dimensional wave displacement field. At least four measurements with linear‐independent encoding directions are necessary to correct for spurious phase contributions if effects from imaging gradients are non‐negligible. In MRE, three encoding schemes have been used: unbalanced four‐ and six‐point and balanced four‐point (‘tetrahedral’) encoding. The first two sensitize to motion with orthogonal gradients, with the four‐point method acquiring a single reference scan without motion sensitization, whereas three additional scans with inverted gradients are used with six‐point encoding, leading to two‐fold higher displacement‐to‐noise ratio (DNR) and 50% longer scan duration. Balanced four‐point (tetrahedral) encoding encodes along the four diagonals of a cube, with one direction serving as a reference for the other three encoding directions, similar to four‐point encoding. The objective of this work is to introduce a theoretical framework to compare different motion sensitization strategies with respect to their motion encoding efficiency in two fundamental encoding limits, the gradient strength limit and the dynamic range limit, which are both placed in relation to conventional gradient recalled echo (GRE)‐ and spin echo (SE)‐based MRE sequences. We apply the framework to the three aforementioned schemes and show that the motion encoding efficiency of unbalanced four‐ and six‐point encoding schemes in the gradient‐limited regime can be increased by a factor of 1.5 when using all physical gradient channels concurrently. Furthermore, it is demonstrated that reversing the direction of the reference in balanced four‐point (tetrahedral) encoding results in the Hadamard encoding scheme, which leads to increased DNR by compared with balanced four‐point encoding and 2.8 compared with unbalanced four‐point encoding. As an example, we show that optimal encoding can be utilized to reduce the acquisition time of standard liver MRE in vivo from four to two breath holds.     

Magnetic resonance elastography using 3D gradient echo measurements of steady-state motion

Magnetic resonance elastography (MRE) is an important new method used to measure the elasticity or stiffness of tissues in vivo. While there are many possible applications of MRE, breast cancer detection and classification is currently the most common. Several groups have been developing methods based on MR and ultrasound (US). MR or US is used to estimate the displacements produced by either quasi-static compression or dynamic vibration of the tissue. An important advantage of MRE is the possibility of measuring displacements accurately in all three directions. The central problem in most versions of MRE is recovering elasticity information from the measured displacements. In previous work, we have presented simulation results in two and three dimensions that were promising. In this article, accurate reconstructions of elasticity images from 3D, steady-state experimental data are reported. These results are significant because they demonstrate that the process is truly three-dimensional even for relatively simple geometries and phantoms. Further, they show that the integration of displacement data acquisition and elastic property reconstruction has been successfully achieved in the experimental setting. This process involves acquiring volumetric MR phase images with prescribed phase offsets between the induced mechanical motion and the motion-encoding gradients, converting this information into a corresponding 3D displacement field and estimating the concomitant 3D elastic property distribution through model-based image reconstruction. Fully 3D displacement fields and resulting elasticity images are presented for single and multiple inclusion gel phantoms.     

Single breath‐hold 3D cardiac T1 mapping using through‐time spiral GRAPPA

The quantification of cardiac T₁ relaxation time holds great potential for the detection of various cardiac diseases. However, as a result of both cardiac and respiratory motion, only one two‐dimensional T₁ map can be acquired in one breath‐hold with most current techniques, which limits its application for whole heart evaluation in routine clinical practice. In this study, an electrocardiogram (ECG)‐triggered three‐dimensional Look–Locker method was developed for cardiac T₁ measurement. Fast three‐dimensional data acquisition was achieved with a spoiled gradient‐echo sequence in combination with a stack‐of‐spirals trajectory and through‐time non‐Cartesian generalized autocalibrating partially parallel acquisition (GRAPPA) acceleration. The effects of different magnetic resonance parameters on T₁ quantification with the proposed technique were first examined by simulating data acquisition and T₁ map reconstruction using Bloch equation simulations. Accuracy was evaluated in studies with both phantoms and healthy subjects. These results showed that there was close agreement between the proposed technique and the reference method for a large range of T₁ values in phantom experiments. In vivo studies further demonstrated that rapid cardiac T₁ mapping for 12 three‐dimensional partitions (spatial resolution, 2 × 2 × 8 mm³) could be achieved in a single breath‐hold of ~12 s. The mean T₁ values of myocardial tissue and blood obtained from normal volunteers at 3 T were 1311 ± 66 and 1890 ± 159 ms, respectively. In conclusion, a three‐dimensional T₁ mapping technique was developed using a non‐Cartesian parallel imaging method, which enables fast and accurate T₁ mapping of cardiac tissues in a single short breath‐hold.     

Evaluation of highly accelerated real‐time cardiac cine MRI in tachycardia

Electrocardiogram (ECG)‐gated breath‐hold cine MRI is considered to be the gold standard test for the assessment of cardiac function. However, it may fail in patients with arrhythmia, impaired breath‐hold capacity and poor ECG gating. Although ungated real‐time cine MRI may mitigate these problems, commercially available real‐time cine MRI pulse sequences using parallel imaging typically yield relatively poor spatiotemporal resolution because of their low image acquisition efficiency. As an extension of our previous work, the purpose of this study was to evaluate the diagnostic quality and accuracy of eight‐fold‐accelerated real‐time cine MRI with compressed sensing (CS) for the quantification of cardiac function in tachycardia, where it is challenging for real‐time cine MRI to provide sufficient spatiotemporal resolution. We evaluated the performances of eight‐fold‐accelerated cine MRI with CS, three‐fold‐accelerated real‐time cine MRI with temporal generalized autocalibrating partially parallel acquisitions (TGRAPPA) and ECG‐gated breath‐hold cine MRI in 21 large animals with tachycardia (mean heart rate, 104 beats per minute) at 3T. For each cine MRI method, two expert readers evaluated the diagnostic quality in four categories (image quality, temporal fidelity of wall motion, artifacts and apparent noise) using a Likert scale (1–5, worst to best). One reader evaluated the left ventricular functional parameters. The diagnostic quality scores were significantly different between the three cine pulse sequences, except for the artifact level between CS and TGRAPPA real‐time cine MRI. Both ECG‐gated breath‐hold cine MRI and eight‐fold accelerated real‐time cine MRI yielded all four scores of ≥ 3.0 (acceptable), whereas three‐fold‐accelerated real‐time cine MRI yielded all scores below 3.0, except for artifact (3.0). The left ventricular ejection fraction (LVEF) measurements agreed better between ECG‐gated cine MRI and eight‐fold‐accelerated real‐time cine MRI (mean difference, –1.6%) than between ECG‐gated cine MRI and three‐fold‐accelerated real‐time cine MRI (mean difference, –5.7%). Eight‐fold‐accelerated real‐time cine MRI with CS yields acceptable diagnostic quality and relatively accurate LVEF measurements in the challenging setting of tachycardia. Copyright © 2013 John Wiley & Sons, Ltd.     

T1 measurement using a short acquisition period for quantitative cardiac applications

Myocardial signal intensity curves for myocardial perfusion studies may be made quantitative by the use of T1 measurements made after the first-pass of contrast agent. A short data acquisition method for T1 mapping is presented in which all data for each T1 map are acquired in a short breath hold, and the slice geometry and timing in the cardiac cycle exactly match that of the dynamic first-pass perfusion sequence. This allows accurate image registration of the T1 map with the first-pass series of images. The T1 method is based on varying the preparation-pulse delay time of a saturation recovery sequence, and in this implementation employs an ECG-triggered, single-shot, spoiled gradient echo technique with SENSE reconstruction. The method allows T1 estimates of three slices to be made in fifteen heartbeats. For a range of samples with T1 values equivalent to those found in the myocardium during the first-pass of contrast agent, T1 estimates were accurate to within 6%, and the variation between slices was 2% or less.     

Accelerated high‐frame‐rate mouse heart cine‐MRI using compressed sensing reconstruction

We introduce a new protocol to obtain very high‐frame‐rate cinematographic (Cine) MRI movies of the beating mouse heart within a reasonable measurement time. The method is based on a self‐gated accelerated fast low‐angle shot (FLASH) acquisition and compressed sensing reconstruction. Key to our approach is that we exploit the stochastic nature of the retrospective triggering acquisition scheme to produce an undersampled and random k–t space filling that allows for compressed sensing reconstruction and acceleration. As a standard, a self‐gated FLASH sequence with a total acquisition time of 10 min was used to produce single‐slice Cine movies of seven mouse hearts with 90 frames per cardiac cycle. Two times (2×) and three times (3×) k–t space undersampled Cine movies were produced from 2.5‐ and 1.5‐min data acquisitions, respectively. The accelerated 90‐frame Cine movies of mouse hearts were successfully reconstructed with a compressed sensing algorithm. The movies had high image quality and the undersampling artifacts were effectively removed. Left ventricular functional parameters, i.e. end‐systolic and end‐diastolic lumen surface areas and early‐to‐late filling rate ratio as a parameter to evaluate diastolic function, derived from the standard and accelerated Cine movies, were nearly identical. Copyright © 2012 John Wiley & Sons, Ltd.