An MRI multi‐scanner neuroimaging data harmonization study using RAVEL and ComBat

Abstract

Background Large‐scale multi‐site neuroimaging studies provide higher power for statistical analyses. However, these aggregated datasets are susceptible to unwanted variability resulting from scanner and acquisition protocol differences. To address this problem, a group of intensity normalization and harmonization methods have been developed recently. This study used a paired 1.5T‐3T MRI dataset to evaluate an intensity normalization procedure, RAVEL (Fortin et al., 2016), a data harmonization method, ComBat (Johnson et al., 2007), and a pipeline combining the two. Method Sixteen healthy elderly participants (median 77.5 years (range=70‐79), 25% (n=4) males) underwent 1.5T and 3T T1‐weighted MRI, three months apart. These images (RAW) were used as input for RAVEL (an intensity normalization technique) and/or ComBat (a harmonization procedure) to characterize the cross‐scanner harmonization of FreeSurfer‐extracted summary measures (cortical thickness and/or volume) relevant to Alzheimer’s disease (Zhang et al., 2001). In setting‐1 (RAVEL), measures were extracted from RAVEL‐normalized images. In setting‐2 (ComBat), measures were extracted from RAW images and then harmonized using ComBat. In setting‐3 (RAVEL‐ComBat), ComBat was applied to the measures resulted from setting‐1. Paired t‐tests were used to evaluate cross‐scanner differences of measures in each of these sets. The measures with statistically significant differences were considered with potential scanner differences and decreasing such differences to lose significance, was considered to achieve harmonization. Result Table 1 shows the means of the differences and their 95% CI for statistical significance, indicating 14 measures with potential scanner differences in RAW (blue/yellow under RAW). Out of these 14 measures, RAVEL‐Combat harmonized 13 measures (green under RAVEL‐ComBat) compared to 8 and 9 harmonized measures by RAVEL and ComBat, respectively. For example, for the cortical thickness of the Inferior Temporal Left, the mean of differences for RAW was 0.25mm, 95%CI:(0.16mm, 0.34mm), indicating a statistically significant difference (blue). This difference has been reduced to 0.08mm, 95%CI (‐0.05mm, 0.21mm) using RAVEL‐ComBat (green) and no longer statistically significant. We observed that, both, RAVEL‐ComBat and RAVEL slightly increased the variability of the differences. Conclusion A combination of intensity normalization (RAVEL) and data harmonization (ComBat) performs best in terms of harmonizing data across 1.5T and 3T MRI within the same cohort.

Cerebrovascular lesion loads and accelerated brain aging: insights into the cognitive spectrum

“…In MRI pre-processing stage, we employed a validated software, specifically FreeSurfer , renowned for its suitability in multi-center and multi-scanner studies (Dadar et al, 2020 ). Moreover, in order to reduce any potential impact of various scanner manufacturers on FreeSurfer measurements, we incorporated a validated harmonization technique known as ComBat (Fortin et al, 2018 ; Torbati et al, 2021 ). To ensure the separation of brain features extracted by FreeSurfer from WMH (Dadar et al, 2021b ), only cortical brain features were used in our brain age prediction model.…”

Machine Learning for Alzheimer’s Disease and Related Dementias

“…Using empirical Bayes to improve the estimation of the site, this model can be used to correct several imaging modalities while associating relevant clinical and demographic information. It was originally developed to correct gene expression microarray data [159], being later extended to correct DTI maps [158], cortical thickness measurements [160], or structural MRI [161]. Additional extensions include longitudinal data [162], site effects due to covariance [163], and a generalized additive model in order to handle nonlinear trajectories over the life span [164].…”