Dataset Information

Cross-scanner and cross-protocol multi-shell diffusion MRI data harmonization: Algorithms and results.

ABSTRACT: Cross-scanner and cross-protocol variability of diffusion magnetic resonance imaging (dMRI) data are known to be major obstacles in multi-site clinical studies since they limit the ability to aggregate dMRI data and derived measures. Computational algorithms that harmonize the data and minimize such variability are critical to reliably combine datasets acquired from different scanners and/or protocols, thus improving the statistical power and sensitivity of multi-site studies. Different computational approaches have been proposed to harmonize diffusion MRI data or remove scanner-specific differences. To date, these methods have mostly been developed for or evaluated on single b-value diffusion MRI data. In this work, we present the evaluation results of 19 algorithms that are developed to harmonize the cross-scanner and cross-protocol variability of multi-shell diffusion MRI using a benchmark database. The proposed algorithms rely on various signal representation approaches and computational tools, such as rotational invariant spherical harmonics, deep neural networks and hybrid biophysical and statistical approaches. The benchmark database consists of data acquired from the same subjects on two scanners with different maximum gradient strength (80 and 300 mT/m) and with two protocols. We evaluated the performance of these algorithms for mapping multi-shell diffusion MRI data across scanners and across protocols using several state-of-the-art imaging measures. The results show that data harmonization algorithms can reduce the cross-scanner and cross-protocol variabilities to a similar level as scan-rescan variability using the same scanner and protocol. In particular, the LinearRISH algorithm based on adaptive linear mapping of rotational invariant spherical harmonics features yields the lowest variability for our data in predicting the fractional anisotropy (FA), mean diffusivity (MD), mean kurtosis (MK) and the rotationally invariant spherical harmonic (RISH) features. But other algorithms, such as DIAMOND, SHResNet, DIQT, CMResNet show further improvement in harmonizing the return-to-origin probability (RTOP). The performance of different approaches provides useful guidelines on data harmonization in future multi-site studies.

SUBMITTER: Ning L

PROVIDER: S-EPMC10243465 | biostudies-literature | 2020 Nov

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Cross-scanner and cross-protocol multi-shell diffusion MRI data harmonization: Algorithms and results.

Ning Lipeng L Bonet-Carne Elisenda E Grussu Francesco F Sepehrband Farshid F Kaden Enrico E Veraart Jelle J Blumberg Stefano B SB Khoo Can Son CS Palombo Marco M Kokkinos Iasonas I Alexander Daniel C DC Coll-Font Jaume J Scherrer Benoit B Warfield Simon K SK Karayumak Suheyla Cetin SC Rathi Yogesh Y Koppers Simon S Weninger Leon L Ebert Julia J Merhof Dorit D Moyer Daniel D Pietsch Maximilian M Christiaens Daan D Gomes Teixeira Rui Azeredo RA Tournier Jacques-Donald JD Schilling Kurt G KG Huo Yuankai Y Nath Vishwesh V Hansen Colin C Blaber Justin J Landman Bennett A BA Zhylka Andrey A Pluim Josien P W JPW Parker Greg G Rudrapatna Umesh U Evans John J Charron Cyril C Jones Derek K DK Tax Chantal M W CMW

NeuroImage 20200713

Cross-scanner and cross-protocol variability of diffusion magnetic resonance imaging (dMRI) data are known to be major obstacles in multi-site clinical studies since they limit the ability to aggregate dMRI data and derived measures. Computational algorithms that harmonize the data and minimize such variability are critical to reliably combine datasets acquired from different scanners and/or protocols, thus improving the statistical power and sensitivity of multi-site studies. Different computat ...[more]

PMID: 32673745

Similar Datasets

Project description:The clinical usefulness MRI biomarkers for aging and dementia studies relies on precise brain morphological measurements; however, scanner and/or protocol variations may introduce noise or bias. One approach to address this is post-acquisition scan harmonization. In this work, we evaluate deep learning (neural style transfer, CycleGAN and CGAN), histogram matching, and statistical (ComBat and LongComBat) methods. Participants who had been scanned on both GE and Siemens scanners (cross-sectional participants, known as Crossover (n = 113), and longitudinally scanned participants on both scanners (n = 454)) were used. The goal was to match GE MPRAGE (T1-weighted) scans to Siemens improved resolution MPRAGE scans. Harmonization was performed on raw native and preprocessed (resampled, affine transformed to template space) scans. Cortical thicknesses were measured using FreeSurfer (v.7.1.1). Distributions were checked using Kolmogorov-Smirnov tests. Intra-class correlation (ICC) was used to assess the degree of agreement in the Crossover datasets and annualized percent change in cortical thickness was calculated to evaluate the Longitudinal datasets. Prior to harmonization, the least agreement was found at the frontal pole (ICC = 0.72) for the raw native scans, and at caudal anterior cingulate (0.76) and frontal pole (0.54) for the preprocessed scans. Harmonization with NST, CycleGAN, and HM improved the ICCs of the preprocessed scans at the caudal anterior cingulate (>0.81) and frontal poles (>0.67). In the Longitudinal raw native scans, over- and under-estimations of cortical thickness were observed due to the changing of the scanners. ComBat matched the cortical thickness distributions throughout but was not able to increase the ICCs or remove the effects of scanner changeover in the Longitudinal datasets. CycleGAN and NST performed slightly better to address the cortical thickness variations between scanner change. However, none of the methods succeeded in harmonizing the Longitudinal dataset. CGAN was the worst performer for both datasets. In conclusion, the performance of the methods was overall similar and region dependent. Future research is needed to improve the existing approaches since none of them outperformed each other in terms of harmonizing the datasets at all ROIs. The findings of this study establish framework for future research into the scan harmonization problem.

Project description:The clinical translation of diffusion magnetic resonance imaging (dMRI)-derived quantitative contrasts hinges on robust reproducibility, minimizing both same-scanner and cross-scanner variability. As multi-site data sets, including multi-shell dMRI, expand in scope, enhancing reproducibility across variable MRI systems and MRI protocols becomes crucial. This study evaluates the reproducibility of diffusion kurtosis imaging (DKI) metrics (beyond conventional diffusion tensor imaging (DTI)), at the voxel and region-of-interest (ROI) levels on magnitude and complex-valued dMRI data, using denoising with and without harmonization. We compared same-scanner, cross-scanner, and cross-protocol variability for a multi-shell dMRI protocol (2-mm isotropic resolution, b = 0, 1000, 2000 s/mm2) in 20 subjects. We first evaluated the effectiveness of Marchenko-Pastur Principal Component Analysis (MPPCA) based denoising strategies for both magnitude and complex data to mitigate noise-induced bias and variance, to improve dMRI parametric maps and reproducibility. Next, we examined the impact of denoising under different population analysis approaches, specifically comparing voxel-wise versus region of interest (ROI)-based methods. We also evaluated the role of denoising when harmonizing dMRI across scanners and protocols. The results indicate that DTI and DKI maps visually improve after MPPCA denoising, with noticeably fewer outliers in kurtosis maps. Denoising, either using magnitude or complex dMRI, enhances voxel-wise reproducibility, with test-retest variability of kurtosis indices reduced from 15%-20% without denoising to 5%-10% after denoising. Complex dMRI denoising reduces the noise floor by up to 60%. Denoising not only reduced variability across scans and protocols, but also increased statistical power for low SNR voxel-wise comparisons when comparing cross sectional groups. In conclusion, MPPCA denoising, either over magnitude or complex dMRI data, enhances the reproducibility and precision of higher-order diffusion metrics across same-scanner, cross-scanner, and cross-protocol assessments. The enhancement in data quality and precision facilitates the broader application and acceptance of these advanced imaging techniques in both clinical practice and large-scale neuroimaging studies.

Project description:Parkinson's disease (PD) is a progressive neurodegenerative disorder that affects extensive regions of the central nervous system. In this work, we evaluated the structural connectome of patients with PD, as mapped by diffusion-weighted MRI tractography and a multi-shell, multi-tissue (MSMT) constrained spherical deconvolution (CSD) method to increase the precision of tractography at tissue interfaces. The connectome was mapped with probabilistic MSMT-CSD in 21 patients with PD and in 21 age- and gender-matched controls. Mapping was also performed by deterministic single-shell, single tissue (SSST)-CSD tracking and probabilistic SSST-CSD tracking for comparison. A support vector machine was trained to predict diagnosis based on a linear combination of graph metrics. We showed that probabilistic MSMT-CSD could detect significantly reduced global strength, efficiency, clustering, and small-worldness, and increased global path length in patients with PD relative to healthy controls; by contrast, probabilistic SSST-CSD only detected the difference in global strength and small-worldness. In patients with PD, probabilistic MSMT-CSD also detected a significant reduction in local efficiency and detected clustering in the motor, frontal temporoparietal associative, limbic, basal ganglia, and thalamic areas. The network-based statistic identified a subnetwork of reduced connectivity by MSMT-CSD and probabilistic SSST-CSD in patients with PD, involving key components of the cortico-basal ganglia-thalamocortical network. Finally, probabilistic MSMT-CSD had superior diagnostic accuracy compared with conventional probabilistic SSST-CSD and deterministic SSST-CSD tracking. In conclusion, probabilistic MSMT-CSD detected a greater extent of connectome pathology in patients with PD, including those with cortico-basal ganglia-thalamocortical network disruptions. Connectome analysis based on probabilistic MSMT-CSD may be useful when evaluating the extent of white matter connectivity disruptions in PD.

Dataset Information

Cross-scanner and cross-protocol multi-shell diffusion MRI data harmonization: Algorithms and results.

Publications

Cross-scanner and cross-protocol multi-shell diffusion MRI data harmonization: Algorithms and results.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets