Dataset Information

Harmonization of cortical thickness measurements across scanners and sites.

ABSTRACT: With the proliferation of multi-site neuroimaging studies, there is a greater need for handling non-biological variance introduced by differences in MRI scanners and acquisition protocols. Such unwanted sources of variation, which we refer to as "scanner effects", can hinder the detection of imaging features associated with clinical covariates of interest and cause spurious findings. In this paper, we investigate scanner effects in two large multi-site studies on cortical thickness measurements across a total of 11 scanners. We propose a set of tools for visualizing and identifying scanner effects that are generalizable to other modalities. We then propose to use ComBat, a technique adopted from the genomics literature and recently applied to diffusion tensor imaging data, to combine and harmonize cortical thickness values across scanners. We show that ComBat removes unwanted sources of scan variability while simultaneously increasing the power and reproducibility of subsequent statistical analyses. We also show that ComBat is useful for combining imaging data with the goal of studying life-span trajectories in the brain.

SUBMITTER: Fortin JP

PROVIDER: S-EPMC5845848 | biostudies-literature | 2018 Feb

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Harmonization of cortical thickness measurements across scanners and sites.

Fortin Jean-Philippe JP Cullen Nicholas N Sheline Yvette I YI Taylor Warren D WD Aselcioglu Irem I Cook Philip A PA Adams Phil P Cooper Crystal C Fava Maurizio M McGrath Patrick J PJ McInnis Melvin M Phillips Mary L ML Trivedi Madhukar H MH Weissman Myrna M MM Shinohara Russell T RT

NeuroImage 20171117

With the proliferation of multi-site neuroimaging studies, there is a greater need for handling non-biological variance introduced by differences in MRI scanners and acquisition protocols. Such unwanted sources of variation, which we refer to as "scanner effects", can hinder the detection of imaging features associated with clinical covariates of interest and cause spurious findings. In this paper, we investigate scanner effects in two large multi-site studies on cortical thickness measurements ...[more]

PMID: 29155184

Similar Datasets

Project description:Results of neuroimaging datasets aggregated from multiple sites may be biased by site-specific profiles in participants' demographic and clinical characteristics, as well as MRI acquisition protocols and scanning platforms. We compared the impact of four different harmonization methods on results obtained from analyses of cortical thickness data: (1) linear mixed-effects model (LME) that models site-specific random intercepts (LMEINT), (2) LME that models both site-specific random intercepts and age-related random slopes (LMEINT+SLP), (3) ComBat, and (4) ComBat with a generalized additive model (ComBat-GAM). Our test case for comparing harmonization methods was cortical thickness data aggregated from 29 sites, which included 1,340 cases with posttraumatic stress disorder (PTSD) (6.2-81.8 years old) and 2,057 trauma-exposed controls without PTSD (6.3-85.2 years old). We found that, compared to the other data harmonization methods, data processed with ComBat-GAM was more sensitive to the detection of significant case-control differences (Χ2(3) = 63.704, p < 0.001) as well as case-control differences in age-related cortical thinning (Χ2(3) = 12.082, p = 0.007). Both ComBat and ComBat-GAM outperformed LME methods in detecting sex differences (Χ2(3) = 9.114, p = 0.028) in regional cortical thickness. ComBat-GAM also led to stronger estimates of age-related declines in cortical thickness (corrected p-values < 0.001), stronger estimates of case-related cortical thickness reduction (corrected p-values < 0.001), weaker estimates of age-related declines in cortical thickness in cases than controls (corrected p-values < 0.001), stronger estimates of cortical thickness reduction in females than males (corrected p-values < 0.001), and stronger estimates of cortical thickness reduction in females relative to males in cases than controls (corrected p-values < 0.001). Our results support the use of ComBat-GAM to minimize confounds and increase statistical power when harmonizing data with non-linear effects, and the use of either ComBat or ComBat-GAM for harmonizing data with linear effects.

Project description:The cerebral cortex is a distinctive part of the mammalian nervous system, displaying a spatial variety in cyto-, chemico-, and myelinoarchitecture. As part of a rich history of histological findings, pioneering anatomists von Economo and Koskinas provided detailed mappings on the cellular structure of the human cortex, reporting on quantitative aspects of cytoarchitecture of cortical areas. Current day investigations into the structure of human cortex have embraced technological advances in Magnetic Resonance Imaging (MRI) to assess macroscale thickness and organization of the cortical mantle in vivo. However, direct comparisons between current day MRI estimates and the quantitative measurements of early anatomists have been limited. Here, we report on a simple, but nevertheless important cross-analysis between the histological reports of von Economo and Koskinas on variation in thickness of the cortical mantle and MRI derived measurements of cortical thickness. We translated the von Economo cortical atlas to a subdivision of the commonly used Desikan-Killiany atlas (as part of the FreeSurfer Software package and a commonly used parcellation atlas in studies examining MRI cortical thickness). Next, values of "width of the cortical mantle" as provided by the measurements of von Economo and Koskinas were correlated to cortical thickness measurements derived from high-resolution anatomical MRI T1 data of 200+ subjects of the Human Connectome Project (HCP). Cross-correlation revealed a significant association between group-averaged MRI measurements of cortical thickness and histological recordings (r?=?0.54, P?<?0.001). Further validating such a correlation, we manually segmented the von Economo parcellation atlas on the standardized Colin27 brain dataset and applied the obtained three-dimensional von Economo segmentation atlas to the T1 data of each of the HCP subjects. Highly consistent with our findings for the mapping to the Desikan-Killiany regions, cross-correlation between in vivo MRI cortical thickness and von Economo histology-derived values of cortical mantle width revealed a strong positive association (r?=?0.62, P?<?0.001). Linking today's state-of-the-art T1-weighted imaging to early histological examinations our findings indicate that MRI technology is a valid method for in vivo assessment of thickness of human cortex.

Dataset Information

Harmonization of cortical thickness measurements across scanners and sites.

Publications

Harmonization of cortical thickness measurements across scanners and sites.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets