Dataset Information

Harmonization of diffusion MRI data sets with adaptive dictionary learning.

ABSTRACT: Diffusion magnetic resonance imaging can indirectly infer the microstructure of tissues and provide metrics subject to normal variability in a population. Potentially abnormal values may yield essential information to support analysis of controls and patients cohorts, but subtle confounds could be mistaken for purely biologically driven variations amongst subjects. In this work, we propose a new harmonization algorithm based on adaptive dictionary learning to mitigate the unwanted variability caused by different scanner hardware while preserving the natural biological variability of the data. Our harmonization algorithm does not require paired training data sets, nor spatial registration or matching spatial resolution. Overcomplete dictionaries are learned iteratively from all data sets at the same time with an adaptive regularization criterion, removing variability attributable to the scanners in the process. The obtained mapping is applied directly in the native space of each subject toward a scanner-space. The method is evaluated with a public database which consists of two different protocols acquired on three different scanners. Results show that the effect size of the four studied diffusion metrics is preserved while removing variability attributable to the scanner. Experiments with alterations using a free water compartment, which is not simulated in the training data, shows that the modifications applied to the diffusion weighted images are preserved in the diffusion metrics after harmonization, while still reducing global variability at the same time. The algorithm could help multicenter studies pooling their data by removing scanner specific confounds, and increase statistical power in the process.

SUBMITTER: St-Jean S

PROVIDER: S-EPMC7555079 | biostudies-literature | 2020 Aug

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Harmonization of diffusion MRI data sets with adaptive dictionary learning.

St-Jean Samuel S Viergever Max A MA Leemans Alexander A

Human brain mapping 20200826 16

Diffusion magnetic resonance imaging can indirectly infer the microstructure of tissues and provide metrics subject to normal variability in a population. Potentially abnormal values may yield essential information to support analysis of controls and patients cohorts, but subtle confounds could be mistaken for purely biologically driven variations amongst subjects. In this work, we propose a new harmonization algorithm based on adaptive dictionary learning to mitigate the unwanted variability ca ...[more]

PMID: 32851729

Similar Datasets

Project description: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.

Project description:ObjectivesAcquisition-related differences in diffusion magnetic resonance imaging (dMRI) hamper pooling of multicentre data to achieve large sample sizes. A promising solution is to harmonize the raw diffusion signal using rotation invariant spherical harmonic (RISH) features, but this has not been tested in elderly subjects. Here we aimed to establish if RISH harmonization effectively removes acquisition-related differences in multicentre dMRI of elderly subjects with cerebral small vessel disease (SVD), while preserving sensitivity to disease effects.MethodsFive cohorts of patients with SVD (N = 397) and elderly controls (N = 175) with 3 Tesla MRI on different systems were included. First, to establish effectiveness of harmonization, the RISH method was trained with data of 13 to 15 age and sex-matched controls from each site. Fractional anisotropy (FA) and mean diffusivity (MD) were compared in matched controls between sites using tract-based spatial statistics (TBSS) and voxel-wise analysis, before and after harmonization. Second, to assess sensitivity to disease effects, we examined whether the contrast (effect sizes of FA, MD and peak width of skeletonized MD - PSMD) between patients and controls within each site remained unaffected by harmonization. Finally, we evaluated the association between white matter hyperintensity (WMH) burden, FA, MD and PSMD using linear regression analyses both within individual cohorts as well as with pooled scans from multiple sites, before and after harmonization.ResultsBefore harmonization, significant differences in FA and MD were observed between matched controls of different sites (p < 0.05). After harmonization these site-differences were removed. Within each site, RISH harmonization did not alter the effect sizes of FA, MD and PSMD between patients and controls (relative change in Cohen's d = 4 %) nor the strength of association with WMH volume (relative change in R2 = 2.8 %). After harmonization, patient data of all sites could be aggregated in a single analysis to infer the association between WMH volume and FA (R2 = 0.62), MD (R2 = 0.64), and PSMD (R2 = 0.60).ConclusionsWe showed that RISH harmonization effectively removes acquisition-related differences in dMRI of elderly subjects while preserving sensitivity to SVD-related effects. This study provides proof of concept for future multicentre SVD studies with pooled datasets.

Project description:BackgroundPseudotime estimation from dynamic single-cell transcriptomic data enables characterisation and understanding of the underlying processes, for example developmental processes. Various pseudotime estimation methods have been proposed during the last years. Typically, these methods start with a dimension reduction step because the low-dimensional representation is usually easier to analyse. Approaches such as PCA, ICA or t-SNE belong to the most widely used methods for dimension reduction in pseudotime estimation methods. However, these methods usually make assumptions on the derived dimensions, which can result in important dataset properties being missed. In this paper, we suggest a new dictionary learning based approach, dynDLT, for dimension reduction and pseudotime estimation of dynamic transcriptomic data. Dictionary learning is a matrix factorisation approach that does not restrict the dependence of the derived dimensions. To evaluate the performance, we conduct a large simulation study and analyse 8 real-world datasets.ResultsThe simulation studies reveal that firstly, dynDLT preserves the simulated patterns in low-dimension and the pseudotimes can be derived from the low-dimensional representation. Secondly, the results show that dynDLT is suitable for the detection of genes exhibiting the simulated dynamic patterns, thereby facilitating the interpretation of the compressed representation and thus the dynamic processes. For the real-world data analysis, we select datasets with samples that are taken at different time points throughout an experiment. The pseudotimes found by dynDLT have high correlations with the experimental times. We compare the results to other approaches used in pseudotime estimation, or those that are method-wise closely connected to dictionary learning: ICA, NMF, PCA, t-SNE, and UMAP. DynDLT has the best overall performance for the simulated and real-world datasets.ConclusionsWe introduce dynDLT, a method that is suitable for pseudotime estimation. Its main advantages are: (1) It presents a model-free approach, meaning that it does not restrict the dependence of the derived dimensions; (2) Genes that are relevant in the detected dynamic processes can be identified from the dictionary matrix; (3) By a restriction of the dictionary entries to positive values, the dictionary atoms are highly interpretable.

Dataset Information

Harmonization of diffusion MRI data sets with adaptive dictionary learning.

Publications

Harmonization of diffusion MRI data sets with adaptive dictionary learning.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets