High Quality Preclinical 4D-Flow Phase Contrast Imaging
Moritz Braig1, Jochen Leupold1, Ko Cheng-Wen2, Marius Menza1, Juergen Hennig1, Jan Korvink3, and Dominik von Elverfeldt1

1University Medical Center Freiburg, Freiburg, Germany, 2Dept. Computer Science and Engineering, National Sun Yat-sen University, Kaohsiung, Taiwan, 3Institute of Microstructure Technology, Karlsruhe Institute of Technology, Karlsruhe, Germany

Synopsis

So far preclinical 4D-Flow MRI has not been able to deliver an analysis of complex flow due to low resolution. The presented framework and improvements allow high quality data acquisitions with a reduced measurement time and the possibility to visualize regional flow abnormalities. An automatic magnitude segmentation in every timeframe combines anatomic information with the underlying blood flow showing even small vessels. It will draw new conclusions in mouse models of cardiovascular diseases as a valuable tool for preclinical researchers.

Purpose

4D-Flow Phase Contrast Magnetic Resonance Imaging has been shown useful for fundamental research and clinical applications. This technique could derive useful insights about the cardiovascular function and function changes due to diseases in murine models. So far long scan times, low SNR, low resolution and image artefacts compromise image quality and therefore the ability to use this method. We present an improved steady state sequence design in combination with a cryogenic coil and a sophisticated reconstruction process. It allows to characterize aortic blood flow in comparable quality to humans and to derive common quantitative parameters.

Methods

For the preclinical MR imaging a Bruker Biospec 70/20 USR equipped with a two channel transmit/receive cryogenic (26 K helium cooled, Bruker) mouse head surface coil was used. This coil shows a factor 3-5 improvement in SNR compared to conventional room temperature coils.1 Physiological observation and sequence gating was done with a SA Instruments 1030 (Stony Brook, NY 11790) monitoring device. All animals (n = 3) were measured under the approval of the ethics committee (Ref: G-14-91). A balanced four point phase contrast sequence from the manufacturer was redesigned to keep the magnetization in a steady state, by an online trigger decision using a double gating strategy. The gating device was modified to generate two trigger outputs, one to determine the respiration state and the other for the R-peak of the ECG. The steady state is interrupted a few millisecond when detecting the next R-peak from the ECG. The data was reconstructed using Matlab. In order to be able to segment and display smaller vessels a POCS algorithm was implemented in the magnitude reconstruction. A multilevel wavelet based segmentation algorithm with a hard threshold was developed. It automatically segments the aorta in 3D in every cine frame from the magnitude data and detects static tissue that is used for eddy current correction. Results with Daubechies wavelets showed a good agreement with a manual segmentation. The phase offset in the velocity data was corrected by a 3D fitting algorithm based on Legendre Polynoms.2 The data was then fed in our default flow processing pipeline used for human flow post-processing as described in literature.3,4 The data was visualized in Ensight (CEI Inc). Parameters: FOV = 20x16x14 mm³, TE/TR = 2/5; Matrix = 60x50x50 (interpolated to 125x100x88), VENC = 150 cm/s; FA = 15°, Acquisition time in vivo ~ 45 min, POCS iterations = 20.

Results

It was possible to acquire the blood flow in the aorta with a very high resolution showing complex flow patterns even in the aortic branches. The images show an uniform signal intensity over the cardiac cycle. The wavelet segmentation allowed a magnitude segmentation mask in every time frame and could differentiate static tissue from moving tissue. It was possible to remove the phase offsets by means of the 3D phase offset correction based on Legendre Polynomials. Common parameters to quantify the blood flow were derived. For example peak speed 127 ± 29 cm/s, flow through a plane in the ascending aorta was equal to 62 ± 15 ml/min.

Discussion

With the new steady state sequence design it was possible to remove commonly neglected steady state artefacts and to achieve an uniform signal intensity over the cardiac cycle. This facilitates the segmentation process as it is based on hard thresholding. Higher SNR due to the cryogenic coil, redundantizes signal averaging. Thus measurement time is reduced by a factor of three compared to published data5 (at equal resolutions). The wavelet segmentation allowed a magnitude segmentation mask in every time frame and a superior static tissue detection compared to conventional 2D methods (based on standard deviation over time). Even the carotid arteries are well segmented which are not covered with thresholded speed-sum-squares masking methods from the velocity information. Several improvements in sequence design, image reconstruction, post-processing in combination with a highly sensitive coil enables to achieve a sufficient image quality for murine high resolution 4D-Flow. Acceleration methods like GRAPPA, SENSE could be used to speed up the acquisition. However the resulting loss of SNR could impede flow quantification in small vessels.

Conclusion

The presented framework - to our knowledge- allows for the first time to generate data quality high enough to investigate preclinical questions. Our framework -to our knowlege- allows for the first time to study complex flow patterns at high spatial resolutions. It will permit to study flow changes non-invasively in mouse models. Thus allowing for example to investigate the development of atherosclerotic plaques or valvular heart disease, which could provide a better understanding of these lethal developments.

Acknowledgements

German Research Foundation (DFG), grant number EL 534/6-1

References

1. Wagenhaus, B. et al. Functional and Morphological Cardiac Magnetic Resonance Imaging of Mice Using a Cryogenic Quadrature Radiofrequency Coil. PLoS ONE 7, (2012).
2. Testud, F & Zaitsev, M. B0 Field Monitoring by Air-Matched Phantoms. in Proceedings of the 17rd Annual Meeting of Intl. Soc. Magn. Reson. Med., Honolulu, Hawaii 2791 (2009).
3. Stalder, A. F. et al. Quantitative 2D and 3D phase contrast MRI: optimized analysis of blood flow and vessel wall parameters. Magn. Reson. Med. 60, 1218–1231 (2008).
4. Bock, J., Kreher, W., Hennig, J. & Markl, M. Optimized pre-processing of time-resolved 2d and 3d phase contrast MRI data. in In Proceedings of 15th Scientific Meeting of International Society for Magnetic Resonance in Medicine, (2007).
5. Bovenkamp, P. R. et al. Velocity mapping of the aortic flow at 9.4 T in healthy mice and mice with induced heart failure using time-resolved three-dimensional phase-contrast MRI (4D PC MRI). Magma N. Y. N (2014). doi:10.1007/s10334-014-0466-z
6. Feintuch, A. et al. Hemodynamics in the mouse aortic arch as assessed by MRI, ultrasound, and numerical modeling. Am. J. Physiol. Heart Circ. Physiol. 292, H884–892 (2007).

Figures

Figure 1:

Streamline visualization in systole overlain by an automated magnitude segmentation. There is a visible signal decrease in the descending aorto due to the use of a surface coil.


Figure 2:

Automatic magnitude segmentation. Besides the aorta the main branches are well defined. A simple rectangle (shown in left image) is used to remove remaining unwanted tissue. This allows a clear visualization of the region of interest (right image).


Figure 3:

Exemplary result: flow through a plane perpendicular to the ascending aorta over the cardiac cycle. Previous results6 could be confirmed.




Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
0324