Quantitative MRI Preprocessing: Effects of Tissue-Specific Smoothing Approaches on Statistical inference

Quantitative MRI (qMRI) provides detailed measurements of brain tissue properties that can be used to study aging and microstructural changes. However, traditional spatial smoothing techniques can blur tissue boundaries and introduce partial-volume effects, which may impact both the statistical power and precise localization of results.

In this study, three tissue-specific smoothing methods were evaluated to understand their effects on statistical analyses. These methods included a linear tissue-weighted compensated approach (TWS), a generalized nonlinear tissue-masked approach (gTSPOON), and an intensity-weighted edge-preserving method based on the Smallest Univalue Segment Assimilating Nucleus (SUSAN) algorithm.

The research compared these strategies using a lifespan qMRI dataset of 138 healthy individuals aged 19 to 75 years, analyzing maps of myelin water fraction (MTsat), proton density (PD), R1 and R2* relaxation times. The generalized TSPOON (gTSPOON) method used tissue-specific masks derived from probabilistic segmentation, while all three approaches were parameterized to achieve similar spatial smoothing.

The study found that TWS and gTSPOON produced similar spatial patterns of age-related effects across all qMRI parameters and tissue types. However, TWS tended to identify more significant voxels and clusters, indicating slightly higher statistical sensitivity due to its slightly broader effective smoothing and lower RESEL counts. In contrast, SUSANs generated fewer significant voxels and clusters, reflecting a much lower effective smoothness and substantially higher RESEL counts.

Despite these differences in statistical sensitivity, voxel-wise likelihood (LL) analyses revealed that each smoothing approach had distinct anatomical preferences. TWS provided the best model fit within gray matter (GM) regions, while gTSPOON performed better in homogeneous white matter (WM) areas. SUSANs excelled at the interfaces between GM and WM, particularly at sulcal and gyral transitions, suggesting improved preservation of sharp anatomical gradients.

The results showed that the choice of smoothing method significantly influenced both statistical sensitivity and the quality of voxel-wise model fitting in qMRI analyses. While TWS and gTSPOON yielded consistent results, SUSANs tended to enhance model fit at tissue boundaries. Importantly, voxel-wise log-likelihood mapping indicated that no single smoothing strategy is universally optimal across the entire brain.

Instead, each method performs best in specific anatomical regions, highlighting the need for region-dependent optimization of smoothing strategies tailored to particular neuroanatomical structures and biological processes, including age-related changes in the brain.