Non-linear registration improves statistical power to detect hippocampal atrophy in aging and dementia

• Published in: NeuroImage clinical 2019-01, Vol.23, p.101902, Article 101902
• Main Authors: Bartel, F, Visser, M, de Ruiter, M, Belderbos, J, Barkhof, F, Vrenken, H, de Munck, J C, van Herk, M
• Format: Article
• Language: English
• Subjects: Aged; Aged, 80 and over; Aging - pathology; Alzheimer Disease - diagnostic imaging; Alzheimer Disease - pathology; Atrophy - pathology; Cognitive Dysfunction - diagnostic imaging; Cognitive Dysfunction - pathology; Databases, Factual; Female; Hippocampus - diagnostic imaging; Hippocampus - pathology; Humans; Image Processing, Computer-Assisted - methods; Magnetic Resonance Imaging - methods; Male; Neuroimaging - methods; Regular
• ISSN: 2213-1582; 2213-1582
• DOI: 10.1016/j.nicl.2019.101902

To compare the performance of different methods for determining hippocampal atrophy rates using longitudinal MRI scans in aging and Alzheimer's disease (AD). Quantifying hippocampal atrophy caused by neurodegenerative diseases is important to follow the course of the disease. In dementia, the efficacy of new therapies can be partially assessed by measuring their effect on hippocampal atrophy. In radiotherapy, the quantification of radiation-induced hippocampal volume loss is of interest to quantify radiation damage. We evaluated plausibility, reproducibility and sensitivity of eight commonly used methods to determine hippocampal atrophy rates using test-retest scans. Manual, FSL-FIRST, FreeSurfer, multi-atlas segmentation (MALF) and non-linear registration methods (Elastix, NiftyReg, ANTs and MIRTK) were used to determine hippocampal atrophy rates on longitudinal T1-weighted MRI from the ADNI database. Appropriate parameters for the non-linear registration methods were determined using a small training dataset (N = 16) in which two-year hippocampal atrophy was measured using test-retest scans of 8 subjects with low and 8 subjects with high atrophy rates. On a larger dataset of 20 controls, 40 mild cognitive impairment (MCI) and 20  AD patients, one-year hippocampal atrophy rates were measured. A repeated measures ANOVA analysis was performed to determine differences between controls, MCI and AD patients. For each method we calculated effect sizes and the required sample sizes to detect one-year volume change between controls and MCI (N ) and between controls and AD (N ). Finally, reproducibility of hippocampal atrophy rates was assessed using within-session rescans and expressed as an average distance measure D , which expresses the difference in atrophy rate, averaged over all subjects. The same D was used to determine the agreement between different methods. Except for MALF, all methods detected a significant group difference between CTRL and AD, but none could find a significant difference between the CTRL and MCI. FreeSurfer and MIRTK required the lowest sample sizes (FreeSurfer: N  = 115, N  = 17 with D  = 3.26%; MIRTK: N  = 97, N  = 11 with D  = 3.76%), while ANTs was most reproducible (N  = 162, N  = 37 with D  = 1.06%), followed by Elastix (N  = 226, N  = 15 with D  = 1.78%) and NiftyReg (N  = 193, N  = 14 with D  = 2.11%). Manually measured hippocampal atrophy rates required largest sample sizes to detect volume change and were poorly reproduced (