Mapping registration sensitivity in MR mouse brain images

Nonlinear registration algorithms provide a way to estimate structural (brain) differences based on magnetic resonance images. Their ability to align images of different individuals and across modalities has been well-researched, but the bounds of their sensitivity with respect to the recovery of salient morphological differences between groups are unclear. Here we develop a novel approach to simulate deformations on MR brain images to evaluate the ability of two registration algorithms to extract structural differences corresponding to biologically plausible atrophy and expansion. We show that at a neuroanatomical level registration accuracy is influenced by the size and compactness of structures, but do so differently depending on how much change is simulated. The size of structures has a small influence on the recovered accuracy. There is a trend for larger structures to be recovered more accurately, which becomes only significant as the amount of simulated change is large. More compact structures can be recovered more accurately regardless of the amount of simulated change. Both tested algorithms underestimate the full extent of the simulated atrophy and expansion. Finally we show that when multiple comparisons are corrected for at a voxelwise level, a very low rate of false positives is obtained. More interesting is that true positive rates average around 40%, indicating that the simulated changes are not fully recovered. Simulation experiments were run using two fundamentally different registration algorithms and we identified the same results, suggesting that our findings are generalizable across different classes of nonlinear registration algorithms.