Methods 1. All T1-weighted images were segmented with FSL/FIRST.

2. Hippocampal surfaces were reconstructed from the segmentations. 3. Conformal parameterizations of hippocampal surfaces were computed with holomorphic 1-form. 4. Point-to-point correspondences between hippocampal surfaces were computed by the inverse consistent fluid registration method, which was extended to work with general surfaces [1]. 5. Surface deformations were measured by new multivariate statistics [2] consisting of radial distance and multivariate tensor-based morphometry (mTBM). 6. Permutation tests were used to perform group comparisons