Impact of acquisition parameters on cortical thickness and volume derived from Multi-Echo MPRAGE scans
Ross W. Mair1,2, Martin Reuter2,3, and Andre J. van der Kouwe2

1Center for Brain Science, Harvard University, Cambridge, MA, United States, 2A.A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Charlestown, MA, United States, 3Department of Mechanical Engineering, Massachusetts Institute of Technology, Cambridge, MA, United States

Synopsis

The multi-echo MPRAGE (MEMPRAGE) sequence was implemented to reduce signal distortion by acquiring at a higher bandwidth and averaging multiple echoes to recover SNR while providing additional T2* information that can enhance cortical segmentation. A rapid 2-minute MEMPRAGE protocol has been implemented for large multi-center studies. Here, we investigate the impact on morphometric results for the cortex by systematically varying all the parameters modified between the rapid 2-minute scan and a conventional 6-minute structural scan. Small but significant differences in cortical thickness and gray matter volume result from a combination of the use of partial fourier acquisition and lowering the spatial resolution to 1.2mm.

Introduction

Automated MRI-derived measurements of human brain volumes from anatomical scans provide novel insights into normal and abnormal neuroanatomy, but few studies have probed the effects of sequence-dependent parameters on these measurements.1 The multi-echo MPRAGE (MEMPRAGE) sequence reduces signal distortion by using a higher bandwidth and averaging multiple echoes to recover SNR while using variable T2* decays to enhance contrast, and hence, cortical segmentation.2 For maximal efficiency, a rapid 2-minute MEMPRAGE protocol has been implemented for large multi-center studies.3 This yields quantitatively repeatable morphometric information on sub-cortical/white-matter structures in comparison to scans employing a higher-resolution MEMPRAGE protocol with lower image acceleration.4 However, small but significant differences in cortical thickness and gray matter volume were observed between the two protocols.5 Here, we investigate the impact on morphometric results for the cortex by systematically varying all the parameters modified between the rapid 2-minute and the conventional 6-minute MEMPRAGE scans.

Methods

All measurements were performed using a 3.0 T MRI scanner (Siemens Tim Trio). 10 subjects (mean: 28.9 years, 5 male) were scanned using the 12-channel coil. Each session included a standard MEMPRAGE (6:03 min, TR/TI = 2530/1200 ms, matrix 256×256×176, resolution = 1 mm iso (no partial fourier), parallel imaging (PI) acceleration = 2 (p2)); and a rapid MEMPRAGE (2:12 min, TR/TI = 2200/1100 ms, matrix 230×230×144, resolution = 1.2 mm iso (6/8 partial fourier), PI = 4 (p4)). TR/TI; matrix/resolution/partial fourier; and PI acceleration were systematically varied between standard and rapid scans, resulting in a total of eight MEMPRAGE scans for each subject (See Fig. 1). Images were analyzed using FreeSurfer,6 after the eight scans from each subject were aligned using the FreeSurfer robust registration tool.7 An automated parcellation of the cortex, subcortical and white matter structures was performed. The 33 cortical regions of the Desikan-Killiany atlas were combined into five principal cortical lobes.8 Correlation and Bland-Altman difference plots were made for the thickness and volume of each principal cortical lobe determined from each scan.

Results

Three difference-analyses were performed: all 1.0mm scans for all subjects vs. all corresponding 1.2mm scans for all subjects; all p2 scans for all subjects vs. all p4 scans; and all TR2530 vs. all TR2200 scans. Example Bland-Altman plots are shown in Fig. 2. The difference analyses are summarized in Fig. 3. The resolution change resulted in significant differences in cortical thickness and volume, similar to those observed previously.5 TR/TI variation had no effect, while PI acceleration had no significant effect. Bland-Altman slopes were ~ -0.2 – 0.2, indicating no bias with measurement value. Correlation coefficients are summarized in Fig. 4.

Discussion

Small but significant differences in cortical thickness and gray matter volume were previously seen between 6-min and 2-min MEPRAGE scans.5 This could result from either slight blurring of the GM/WM boundary in the rapid scan from the use of partial fourier acquisition, or from the resampling of the 1.2mm data to 1.0mm by the standard Freesurfer processing stream. Surprisingly, the reduced TR/TI from the Freesurfer recommended values had no effect on cortical data, while PI acceleration had no significant effect. It has recently been shown that motion by subjects can result in reduced estimates of cortical thickness.9 While longer scans are more likely to be impacted by subject motion, it is the shorter scans here which show a reduced cortical thickness, indicating motion is likely not responsible for the effects seen in this study. Rather, blurring of the GM/WM boundary must happen at the expense of pure GM, resulting in cortical thinning in the case of the faster scans. For a rapid MEMPRAGE with no impact on cortical thickness compared to conventional scans, a 1.0 mm spatial resolution scan employing the other techniques to speed the acquisition is optimal. Such a scan would take 3 mins 13 sec.

Acknowledgements

NIH Grants P41-RR14075, U24-RR021382.

References

1. Wonderlick, J.S. et al., ‘Reliability of MRI-derived cortical and subcortical morphometric measures: Effects of pulse sequence, voxel geometry, and parallel imaging’, NeuroImage, 2009:44;1324-1333.

2. van der Kouwe, A.J.W. et al., ‘Brain morphometry with multiecho MPRAGE’, NeuroImage, 2008:40:559–569.

3. Holmes, A.J. et al., 'Brain Genomics Superstruct Project initial data release with structural, functional, and behavioral measures', Scientific Data, 2015;2:150031.

4. Mair, R.W. et al., ‘Quantitative Comparison of Morphometric Data from Multi-Echo MPRAGE Variable Acceleration and Different Head Coils’, Proceedings ISMRM, 2013;21:947.

5. Mair, R.W. et al., ‘Validation of Cortical Thickness/Volume Data from Multi-Echo MPRAGE Scans with Variable Acceleration in Young and Elderly Populations’, Proceedings ISMRM, 2014;22:1795.

6. Fischl, B, ‘Measuring the thickness of the human cerebral cortex from magnetic resonance images’, Proc. Natl. Acad. Sci. USA, 2000;97:11050-11055.

7. Reuter, M., ‘Highly Accurate Inverse Consistent Registration: A Robust Approach’, NeuroImage, 2010;53;1181-1196.

8. surfer.nmr.mgh.harvard.edu/fswiki/CorticalParcellation

9. Reuter, M. et al. (2015), ‘Head motion during MRI acquisition reduces gray matter volume and thickness estimates’, NeuroImage, 2015;107;105-115.

Figures

Table showing parameter variations in the eight different MEMPRAGE scans acquired from each subject. Acquisition order was randomized among subjects.

Example Bland-Altman plots showing cortical thickness difference (a, c, e) and gray matter volume difference (b, d, f) in the frontal lobe from Freesurfer analysis of pairs of MEMPRAGE scans. All scans from all subjects with 1.0 mm spatial resolution were compared with all scans from all subjects with 1.2 mm resolution (a,b). Similarly, (c,d) show differences between all scans with TR 2530 and TR 2200ms, and (e,f) show differences between all scans with p2 and p4 acceleration.

Summary difference plots showing cortical thickness difference (a, c, e) and gray matter volume difference (b, d, f) for the 5 principal cortical lobes in both hemi-spheres from Freesurfer analysis of pairs of MEMPRAGE scans. Differences in morphometric results plotted for all scans from all subjects with 1.0 and 1.2 mm spatial resolution (a,b), all scans with TR 2530 ms and TR 2200 ms (c, d), and all scans with p2 and p4 acceleration (e,f).

Correlation coefficients (R2) for gray matter volume of the 5 principal cortical lobes in the right hemisphere, measured with MEMPRAGE scans with either varying resolution, TR or PI acceleration. Results are similar for the left hemisphere. High correlations for GM volume, even when comparing the different resolutions, indicates the difference is systematic, and not random. Automated parcellation in the occipital region failed for one subject, resulting in lower R2.



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