A dataset of multi-contrast unbiased average MRI templates of a Parkinson’s disease population

Subjects

All Parkinson’s patients and healthy controls were recruited through the registry of the Quebec Parkinson Network from 2018 onward and are included in the Montreal Neurological Institute’s (MNI) Open Science Clinical Biospecimen Imaging and Genetic (C-BIG) Repository [9]. After rigorous and detailed informed consent, 126 PD patients and 17 controls were scanned, and their age, sex, Unified Parkinson’s Disease Rating Scale Part III (UPDRS-III) scores, and Hoen and Yahr (H&Y) scores were recorded. Table 1 shows the demographics of each of the population cohorts.

Image acquisition

Each subject was scanned on a Siemens 3T Prisma magnetic resonance imaging (MRI) scanner with the following sequences and accompanying parameters listed in Table 2. From the multi-echo GRE data, several imaging contrasts were created: T2*w image, R2* (i.e., 1/T2*), and CLEAR-SWI [3]. All image contrasts were acquired simultaneously at the same imaging session, and are thus inherently co-registered, per subject.

Image processing and template creation

For each of the imaging modalities, brain masks were generated using BEaST: brain extraction based on nonlocal segmentation technique [10] using the T1w image. After non-local means denoising [7], non-uniformity correction [8], and intensity normalization, each image was registered to stereotaxic space using an in-house PPMI template (which is in ICBM152 stereotaxic space) as the stereotaxic registration target [4].

T1-T2* fusion MR images were created for each subject following methods from [11]. Briefly, T1w and T2*w images from the same subject are added together, using spatially varying weights calculated from the R2* volume (i.e., 1/T2*), such that the T1w structure near the cortex and T2*w deep grey nuclei are preserved. Improved susceptibility-weighted images (CLEAR-SWI) were created from methods described in [3]. All other imaging modalities (T1, T2*w, R2*, T2w, PDw, FLAIR, and NM) were ready for template creation.

The T1-T2* fusion images were used to create the transformations that would drive the unbiased average template creation for all imaging modalities. Template creation follows the methods from [4], where hierarchical registration using Automatic Nonlinear Image Matching and Anatomical labelling (ANIMAL) [12] is employed to register an initial template model (here, the MNI PD25 template [1] [2] was used) to each of the subjects. This process is performed iteratively, correcting for residual mis-registration at each step, until convergence is reached. Once convergence is reached and the unbiased average T1-T2* fusion template is created, the individual subject-to-template transformations are saved and used to create the templates for the remaining contrasts (T1, T2*w, R2*, T2w, PDw, FLAIR, NM, and CLEAR-SWI). The same transformations were used to create 0.5 mm and 0.3 mm templates, resampling the pre-processed MRI volumes to the appropriate resolution in stereotaxic space.