The effect of subject motion on fractional anisotropy estimates: A simulation study of angular bias
Szabolcs David1, Anneriet Heemskerk1, Max Viergever1, and Alexander Leemans1

1Image Sciences Institute, UMC Utrecht, Utrecht, Netherlands

Synopsis

Subject motion during diffusion weighted MRI acquisition is a well-known confounding factor which can affect diffusion tensor metrics. In this simulation study, we investigated whether this confound is related to the change in distribution of diffusion gradient orientations due to the required “B-matrix rotation” during subject motion. According to our findings, subject motion can indeed generate a significant angular bias in addition to the noise bias, suggesting that subject motion itself may cause group differences in diffusion metrics.

Introduction

In many studies, diffusion tensor imaging (DTI) measures, such as fractional anisotropy (FA), of white matter structures are compared between groups, such as different age groups or healthy controls vs patients. The observed differences in these measures, however, may not always be due to biologically related processes, but may originate from group differences in physiological characteristics during acquisition, something already well-studied in the fMRI society [1]. While for diffusion MRI several processing steps have been optimized [2] [3], a recent study with real data showed that FA group differences could still be observed simply based on differences in subject motion between the groups [4]. In this simulation study, we explored whether this confound can be explained by the concomitant change in distribution of diffusion gradient orientations due to the required “B-matrix rotation” during subject motion [5]. We also investigated how the introduced angular bias due to subject motion compares to another well-known confound, the noise bias.

Methods

Subject motion trajectories were simulated along all three coordinate axes to mimic the slow drift of head motion. These trajectories were defined as described previously in [5], where the average absolute angular difference (represented by D), a measure that quantifies the amount of subject motion (rotation only), was varied from 0° (i.e., no subject motion) to 5° (typically considered as a lot of subject motion) and subsequently quantified in steps of 0.5° (500 trajectories per step and per axis). Other simulation settings were defined in line with current acquisition settings for typical DTI studies: 60 diffusion gradient directions distributed isotropically on the unit half-sphere [6]; the mean diffusivity (MD) of the diffusion tensor is 0.7×10-3 mm2/s; the FA was set between 0.2 to 0.8; and the b-value was set at 1000 s/mm2. For each motion trajectory, Rician noise was added (100 instances) with SNR levels equal to 10, 20, 30, 40, and 50. To investigate the relative contribution of the noise bias with respect to the angular bias induced by subject motion, we also simulated diffusion signals as described above, but then without adding this angular bias (i.e., as if there was no subject motion). Mean, standard deviation, coefficient of variation (COV), and percentage difference of FA estimates are computed and evaluated as function of SNR level and amount of subject motion (D).

Results

Figure 1 shows that for a small amount of subject motion (average D = 1°), the introduced angular bias does not add significantly to the noise bias, even at high SNR levels. By contrast, for higher degrees of subject motion (average D = 4°) the angular bias becomes more apparent (Figure 2), especially at higher SNR levels (e.g., see Figure 2a and 2d). Figure 3 shows the estimated FA percentage difference as a function of D for a simulated diffusion tensor with FA = 0.8. One can clearly appreciate that for higher amounts of subject motion, the angular bias becomes more dominant, with a net result of nearly 2% difference for the largest values of D compared to the case without subject motion. Figure 4 shows the COV of the FA percentage difference demonstrating the difference in relative variability between the estimated FA values with and without the induced angular bias.

Discussion and Conclusion

Our main finding is that the underestimation and variance of FA due to motion depends on the underlying FA value and the SNR level. At low SNR levels, the noise will swamp the effects of diffusion gradient imperfections, whereas at higher SNR levels, these disturbances become more apparent and will act as a confounding factor by introducing an additional increase in bias and variability, becoming more dominant with larger amounts of subject motion. Note that this finding is in line with a similarly observed underestimation of FA as reported in work on non-perfect gradient distributions [7]. Our study demonstrates that confounds in FA can be explained by associated changes in the distribution of diffusion gradient orientations due to subject motion.

Acknowledgements

This research is supported by VIDI Grant 639.072.411 from the Netherlands Organisation for Scientific Research (NWO).

References

[1] Power, J.D., Barnes, K.A., Snyder, A.Z., Schlaggar, B.L., Petersen, S.E., Spurious but systematic correlations in functional connectivity MRI networks arise from subject motion Neuroimage. 2012 Feb 1; 59(3):2142-54

[2] Andersson, J.L.R., Sotiropoulos, S.N., An integrated approach to correction for off-resonance effects and subject movement in diffusion MR imaging, Neuroimage. 2015 Oct 20. pii: S1053-8119(15)00920-9.

[3] Andersson, J.L.R., Skare, S., Ashburner, J., How to correct susceptibility distortions in spin-echo echo-planar images: application to diffusion tensor imaging. Neuroimage. 2003 Oct; 20(2):870-88.

[4] Yendiki, A., Koldewyn, K., Kakunoori, S., Kanwisher, N., Fischl, B., Spurious group differences due to head motion in a diffusion MRI study. Neuroimage. 2013 Nov 21;88C:79-90

[5] Leemans, A. and Jones, D. K. (2009), The B-matrix must be rotated when correcting for subject motion in DTI data. Magn Reson Med, 61: 1336–1349.

[6] Jones, D.K., Horsfield, M.A., Simmons, A., Optimal strategies for measuring diffusion in anisotropic systems by magnetic resonance imaging Magn Reson Med. 1999 Sep;42(3):515-25.

[7] Skare, S., Hedehus, M., Moseley, M.E., and Li T.Q., Condition Number as a Measure of Noise Performance of Diffusion Tensor Data Acquisition Schemes with MRI. J Magn Reson. 2000 Dec; 147(2):340-52.

Figures

Estimated FA percentage difference as a function of simulated FA for D = 1° for four different SNR levels (error bars represent standard deviations).

Estimated FA percentage difference as a function of simulated FA for D = 4° for four different SNR levels (error bars represent standard deviations).

Estimated FA percentage difference (means and standard deviations) with and without adding angular bias for three different SNR levels with a simulated FA of 0.8.

Coefficient of variation of estimated FA percentage difference with and without adding angular bias. The subject motion introduces an offset which is more pronounced for high FA values and higher SNR levels as indicated by the shift on the color scale (equal range for both).



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