Optimizing an Automated Processing Pipeline for Regional Macromolecular Proton Fraction and Volume Measurements in the Brain

There is increasing demand for quantitative measurements in neuroimaging research, including absolute quantification of regional myelin density. Fast whole-brain macromolecular proton fraction (MPF) imaging enables acquisition and reconstruction of quantitative maps of MPF, a biomarker of myelin, and has been used in many research applications. To date, no standard protocols have been published for parcellating MPF maps. FreeSurfer is a widely used software package for automated brain parcellation. It is designed for use with high-resolution T1-weighted anatomical images, typically acquired with an MPRAGE sequence. MPF maps inherently provide high contrast between white and gray matter and can potentially be used as input images for parcellation, enabling simultaneous quantitative tissue characterization and volumetric assessment. However, it is unclear how the use of MPF maps as input images with FreeSurfer software affect reliability and repeatability of regional MPF and volume estimates. We evaluated three FreeSurfer parcellation workflows for this purpose. The workflows differed in the amount of neuroimaging data required and computational intensity. MPF and MPRAGE data were acquired in two separate sessions for 11 adults. Reliability was evaluated using mean relative differences in estimates, Bland-Altman analysis, and intra-class correlation coefficients (ICCs). Repeatability was assessed using within-subject coefficients of variation (CVws). All three workflows produced similarly high within-subject repeatability across scans. However, reliability in regional MPF estimates was lower in gray matter than in white matter for workflows that used MPF maps as input. Workflows that use MPF maps as FreeSurfer input may be adequate for estimating mean MPF in GM and WM of cortical parcels and in subcortical regions. However, the MPRAGE-based workflow is recommended when reliable individual-level estimates or volume estimates are required.