To read the full PDF, click here
A Comparison of Methods to Harmonize
Cross-Platform Neuroimaging Data
Delin Sun 1,2, Rajendra A. Morey 1,2 *.
1 Brain Imaging and Analysis Center, Duke University, Durham, NC, USA.
2 Department of Veteran Affairs (VA) Mid-Atlantic Mental Illness Research, Education and Clinical Center, Durham, NC, USA.
* Corresponding Author
Rajendra A. Morey, M.D.
40 Duke Medicine Circle, Room 414
Durham, NC 27710 USA
Phone: 919-286-0411 ext. 6425
Facsimile: 919-416-5912
E-mail: rajendra.morey@ekc23
Methods for Harmonization of Neuroimaging Data
The ENIGMA (Enhancing Neuro Imaging Genetics through Meta-Analysis; http://enigma.ini.usc.edu) Consortium is a large international cooperation project that aggregates neuroimaging, demographic, and clinical data in participants from multiple sites, aiming at increasing statistical power and leading to findings that are more representative of the general population (Dennis et al., 2020; Thompson et al., 2014). Neuroimaging studies based on the ENIGMA cohorts may be more robust and more representative of the population with PTSD given that the same preprocessing protocol was applied in all sites participating in an ENIGMA project. However, the results may still be biased by site-related differences in participants’ demographic and clinical information, acquisition protocols, and scanning platforms (Radua et al., 2020). Different approaches have been developed to harmonize neuroimaging data across sites. The linear mixed-effects (LME) model considers study site as a random factor and produces results with higher statistical power than the widely used random-effects meta-analysis approach (Boedhoe et al., 2017; Favre et al., 2019; van Rooij et al., 2018). The ComBat approach was developed to remove the site differences while preserving variations due to biologically relevant covariates such as age, sex, and diagnosis (Fortin et al., 2018). ComBat overcomes the limitations of the LME method by assuming varying normal distributions at different sites. ComBat has been successfully used in recently published studies to harmonize neuroimaging data including cortical thickness (Fortin et al., 2018), surface area and subcortical volume (Radua et al., 2020), diffusion tensor imaging (Fortin et al., 2017), and resting-state functional connectivity (Yu et al., 2018). The newly published method by Pomponio et al. (2020) further improves the ComBat approach by allowing for non-linear trend in the data using Generalized Additive Models (GAMs; we thus call this approach ComBat-GAM) and assuming different distributions for the multiplicative (scale, or variance) and additive (location, or mean) effects, respectively. The ComBat-GAM method outperformed the traditional ComBat approach in harmonizing cortical and subcortical grey matter volumes in a large sample of 10,477 healthy subjects (3-96 years old) aggregated from 18 sites (Pomponio et al., 2020).
Test Case for Comparing Methods
We select cortical thickness data in trauma-exposed people with and without PTSD to compare cross platform harmonization. PTSD is associated with anatomical and functional alternations in widely distributed regions that play a role in fear learning, threat detection, executive function and emotion regulation, and contextual processing (Shalev, Liberzon, & Marmar, 2017). These include inferior temporal regions (Clausen et al., 2020), default model network (DMN), salience network (SN), central executive network (CEN) and the others (Sripada et al., 2012; Yehuda et al., 2015; Zandvakili et al., 2020). Inconsistent results in PTSD concerning cortical thickness point to an increase (Li et al., 2016), decrease (Lindemer, Salat, Leritz, McGlinchey, & Milberg, 2013; Wrocklage et al., 2017), and even no significant change (Rinne-Albers et al., 2020). Furthermore, little is known about the age-related trajectories of cortical thickness in PTSD, whereas age-appropriate cortical thickness decline is typical in human brain development (Frangou et al., 2020). Military service members with traumatic brain injury and PTSD have faster declines in cortical thickness (Santhanam, Wilson, Oakes, & Weaver, 2019; Savjani, Taylor, Acion, Wilde, & Jorge, 2017) and surface area (Santhanam, Wilson, Mulatya, Oakes, & Weaver, 2019). PTSD-related changes in cortical thickness may be associated with faster chronological age (Katrinli et al., 2020) and DNA methylation age (Wolf, Logue, et al., 2016) (however, see Connolly et al. (2018)). It is possible that PTSD may enhance the association between cortical thickness reduction and metabolic syndrome including cardiometabolic pathology, neurodegeneration, premature aging, and therefore leads to significant medical and cognitive decline (Wolf, Sadeh, et al., 2016).
- Research Aims and Hypotheses
Aims. In the present study, we aim to compare results from three different harmonization approaches (1) LME, (2) ComBat, and (3) ComBat-GAM. We propose comparing methods using cortical thickness data from 29 ENIGMA-PTSD sites as a test case to investigate age-related trajectories of cortical thickness in participants with PTSD.
Hypotheses. We hypothesize more statistically significant between-group differences in age-related cortical thinning through utilization of ComBat-GAM relative to LME and ComBat.
- Analyses plan
Participants. All data were aggregated by the PGC-ENIGMA PTSD Working Group from 29 sites located in five countries (PTSD/Non-PTSD N = 1,344/2,073). Most sites used clinician-administered measures such as the Structured Clinical Interview (SCID) (First, 2015) or Clinician-Administered PTSD Scale (CAPS) (F. W. Weathers et al., 2018; F. W. Weathers, Keane, & Davidson, 2001) to ascertain PTSD diagnosis. One site used a psychiatrist diagnosis for PTSD, and a few additional sites used self-report scales such as the PTSD Checklist (PCL) (F.W. Weathers et al., 2013). While the majority used DSM-IV criteria, a small subset of sites used DSM-5 criteria. Severity of PTSD symptoms was derived from the same measures used for diagnosis, except when a clinician-administered measure (e.g., SCID) lacked severity information, the PCL or Davidson Trauma Scale (DTS) (Davidson et al., 1997) was used. The measures reflect PTSD diagnosis and symptoms in the past month before scanning. For depression, sites used a mix of clinician-administered and self-report instruments to diagnose depression. We will harmonize depression data by assigning participants to major depressive disorder (MDD) or control groups based on a standardized depression severity cut-off. The majority of sites reported depression severity using the Beck Depression Inventory-II (BDI-II) (Beck, Steer, & Brown, 1996), while other commonly used scales were the Center for Epidemiologic Studies Depression Scale (CES-D) (Radloff, 1977), Hamilton Depression Rating Scale (HAM-D) (Hamilton, 1960), Depression Anxiety Stress Scales (DASS) (Lovibond & Lovibond, 1995), and Children’s Depression Inventory (CDI) (Kovacs, 1985). All study sites obtained approval from local institutional review boards or ethics committees. All participants provided written informed consent.
Imaging Data Preprocessing. Anatomical brain images were preprocessed at Duke University through a standardized neuroimaging and QC pipeline developed by the ENIGMA Consortium (http://enigma.ini.usc.edu/protocols/imaging-protocols/) (Logue et al., 2018). Cortical thickness measurements were generated by FreeSurfer software (https://surfer.nmr.mgh.harvard.edu) based on the Destrieux atlas (Destrieux, Fischl, Dale, & Halgren, 2010) that contains 74 regions per hemisphere. In brief, white matter surfaces were deformed toward the gray matter boundary at each surface vertex. Cortical thickness was calculated based on the average distance between the parcellated portions of white and pial surfaces within each region. The estimates of cortical thickness for each region and participant were entered into further respective analyses. The image quality was estimated by the Euler number, which is a measure of the topological complexity of the reconstructed cortical surface (Dale, Fischl, & Sereno, 1999).
ComBat Harmonization. The ComBat approach removes the effects of study sites and/or scanners while preserving inherent biological associations in the data (Fortin et al., 2018). It achieves harmonization by first modeling expected imaging features as linear combinations of the biological variables and site effects whose error term is further modulated by site-specific scaling factors. After that, it applies empirical Bayes to improve the estimation of site parameters for small samples. In the present study, PTSD diagnosis, age, and sex will be designated as biological variables. The ComBat approach will be implemented using the R scripts (https://github.com/Jfortin1/ComBatHarmonization) running on RStudio (ver. 1.3.1073) and R (ver. 4.0.2).
ComBat-GAM Harmonization. Compared to the traditional ComBat approach, the ComBat-GAM method specifies covariates with generic nonlinear effects using the generalized additive models (GAM) (Pomponio et al., 2020). In the present study, we specify age as a nonlinear term in the model for our research purpose. The ComBat-GAM approach will be implemented using the Python scripts (https://github.com/rpomponio/neuroHarmonize) running on Anaconda (ver. 4.8.2) and Python (ver. 3.7.6).
Statistical Models. For the tests based on the LME model, the estimates of cortical thickness from all the participants will be entered into a linear mixed-effects model per cortical region with study site as a random factor, and age, sex, PTSD diagnosis, and age by PTSD diagnosis interaction as covariates. For the analyses based on ComBat and ComBat-GAM approaches, the estimates of cortical thickness from all the participants will be first harmonized through the ComBat or ComBat-GAM approach, and then entered a linear model per cortical region with age, sex, PTSD diagnosis and age by PTSD diagnosis interaction as covariates. The covariate of depression diagnosis will also be added to test the robustness of the findings. The PTSD symptom severity will be utilized to replace PTSD diagnosis in the models to investigate the associations between PTSD severity and age. The false discovery rate (FDR) method (Benjamini & Hochberg, 1995) with a q-value threshold of 5% will be employed to correct for multiple comparisons across the 148 cortical regions. The R packages lme4, lmerTest and sjPlot will be utilized for statistical analyses.
- Analyses Personnel
Delin Sun, PhD
- Resources Needed
1. Preprocessed cortical thickness data by the ENIGMA-PGC PTSD Consortium.
2. Demographic information including age, sex and sites/scanners, and clinical variables including PTSD diagnosis, symptom severity, and comorbid depression.
- Timeline
First results report is anticipated within two months.
- Collaboration
Follow the standard PGC policy regarding secondary analyses, ENIGMA-PGC investigators who wish to substantively contribute to this project are welcome to discuss with and join the proposing group.
- Authorship
We will follow the authorship policy of the PGC-PTSD (https://pgc-ptsd.com/wpcontent/uploads/2017/06/Authorship-Guidelines-PGC-PTSD.pdf)
a) Are you following the authorship policies of the groups involved? YES
b) Will there be a writing group and if so, who will be included? Delin Sun, Rajendra A. Morey, Joaquim Radua, and Aris Sotiras will be included in the writing group.
c) What groups or individuals will be listed as authors? Authors will include the writing group and individuals and group contributors of data and analysis from each site (generally 2-5 co-authors from each site).
d) Will PGC members not listed as named authors be listed at the end of the manuscript? All individuals who meet the criteria established in the PGC-PTSD authorship policy will be co-authors. Other PGC members will not be listed at the end of the manuscript.
e) Will PGC members or groups be listed as “collaborators” on the PubMed abstract page? All individuals and groups who meet the authorship criteria of the PGC-PTSD authorship policy will be listed as collaborators on the PubMed abstract page. No other individuals or groups will be listed.
f) How will funding sources be handled or acknowledged? All funding sources that supported data collection and analysis will be listed in the manuscript.
References
Beck, A. T., Steer, R. A., & Brown, G. K. (1996). Manual for the Beck Depression Inventory-II. San Antonio, TX: Psychological Corporation.
Benjamini, Y., & Hochberg, Y. (1995). Controlling the False Discovery Rate – a Practical and Powerful Approach to Multiple Testing. Journal of the Royal Statistical Society Series B-Statistical Methodology, 57(1), 289-300. doi:DOI 10.1111/j.2517-6161.1995.tb02031.x
Boedhoe, P. S., Schmaal, L., Abe, Y., Ameis, S. H., Arnold, P. D., Batistuzzo, M. C., . . . van den Heuvel, O. A. (2017). Distinct Subcortical Volume Alterations in Pediatric and Adult OCD: A Worldwide Meta- and Mega-Analysis. Am J Psychiatry, 174(1), 60-69. doi:10.1176/appi.ajp.2016.16020201
Clausen, A. N., Clarke, E., Phillips, R. D., Haswell, C., Morey, R. A., Brancu, M., . . . Workgrp, V. M.-A. M. (2020). Combat exposure, posttraumatic stress disorder, and head injuries differentially relate to alterations in cortical thickness in military Veterans. Neuropsychopharmacology, 45(3), 491-498. doi:10.1038/s41386-019-0539-9
Connolly, S. L., Stoop, T. B., Logue, M. W., Orr, E. H., De Vivo, I., Miller, M. W., & Wolf, E. J. (2018). Posttraumatic Stress Disorder Symptoms, Temperament, and the Pathway to Cellular Senescence. Journal of Traumatic Stress, 31(5), 676-686. doi:10.1002/jts.22325
Dale, A. M., Fischl, B., & Sereno, M. I. (1999). Cortical surface-based analysis – I. Segmentation and surface reconstruction. Neuroimage, 9(2), 179-194. doi:DOI 10.1006/nimg.1998.0395
Davidson, J. R., Book, S. W., Colket, J. T., Tupler, L. A., Roth, S., David, D., . . . Feldman, M. E. (1997). Assessment of a new self-rating scale for post-traumatic stress disorder. Psychol Med, 27(1), 153-160. doi:10.1017/s0033291796004229
Dennis, E. L., Baron, D., Bartnik-Olson, B., Caeyenberghs, K., Esopenko, C., Hillary, F. G., . . . Wilde, E. A. (2020). ENIGMA brain injury: Framework, challenges, and opportunities. Human Brain Mapping. doi:10.1002/hbm.25046
Destrieux, C., Fischl, B., Dale, A., & Halgren, E. (2010). Automatic parcellation of human cortical gyri and sulci using standard anatomical nomenclature. Neuroimage, 53(1), 1-15. doi:10.1016/j.neuroimage.2010.06.010
Favre, P., Pauling, M., Stout, J., Hozer, F., Sarrazin, S., Abe, C., . . . Group, E. B. D. W. (2019). Widespread white matter microstructural abnormalities in bipolar disorder: evidence from mega- and meta-analyses across 3033 individuals. Neuropsychopharmacology, 44(13), 2285-2293. doi:10.1038/s41386-019-0485-6
First, M. B. (2015). Structured Clinical Interview for the DSM (SCID). In The Encyclopedia of Clinical Psychology (pp. 1-6).
Fortin, J. P., Cullen, N., Sheline, Y. I., Taylor, W. D., Aselcioglu, I., Cook, P. A., . . . Shinohara, R. T. (2018). Harmonization of cortical thickness measurements across scanners and sites. Neuroimage, 167, 104-120. doi:10.1016/j.neuroimage.2017.11.024
Fortin, J. P., Parker, D., Tunc, B., Watanabe, T., Elliott, M. A., Ruparel, K., . . . Shinohara, R. T. (2017). Harmonization of multi-site diffusion tensor imaging data. Neuroimage, 161, 149-170. doi:10.1016/j.neuroimage.2017.08.047
Frangou, S., Modabbernia, A., Doucet, G. E., Papachristou, E., Williams, S. C., Agartz, I., . . . Dima, D. (2020). Cortical Thickness Trajectories across the Lifespan: Data from 17,075 healthy individuals aged 3-90 years. bioRxiv, 2020.2005.2005.077834. doi:10.1101/2020.05.05.077834
Hamilton, M. (1960). A rating scale for depression. J Neurol Neurosurg Psychiatry, 23, 56-62. doi:10.1136/jnnp.23.1.56
Katrinli, S., Stevens, J., Wani, A. H., Lori, A., Kilaru, V., van Rooij, S. J. H., . . . Smith, A. K. (2020). Evaluating the impact of trauma and PTSD on epigenetic prediction of lifespan and neural integrity. Neuropsychopharmacology, 45(10), 1609-1616. doi:10.1038/s41386-020-0700-5
Kovacs, M. (1985). The Children’s Depression, Inventory (CDI). Psychopharmacol Bull, 21(4), 995-998. Retrieved from https://www.ncbi.nlm.nih.gov/pubmed/4089116
Li, S. G., Huang, X. Q., Li, L. J., Du, F., Li, J., Bi, F., . . . Gong, Q. Y. (2016). Posttraumatic Stress Disorder: Structural Characterization with 3-T MR Imaging. Radiology, 280(2), 537-544. doi:10.1148/radiol.2016150477
Lindemer, E. R., Salat, D. H., Leritz, E. C., McGlinchey, R. E., & Milberg, W. P. (2013). Reduced cortical thickness with increased lifetime burden of PTSD in OEF/OIF Veterans and the impact of comorbid TBI. Neuroimage-Clinical, 2, 601-611. doi:10.1016/j.nicl.2013.04.009
Logue, M. W., van Rooij, S. J. H., Dennis, E. L., Davis, S. L., Hayes, J. P., Stevens, J. S., . . . Morey, R. A. (2018). Smaller Hippocampal Volume in Posttraumatic Stress Disorder: A Multisite ENIGMA-PGC Study: Subcortical Volumetry Results From Posttraumatic Stress Disorder Consortia. Biological Psychiatry, 83(3), 244-253. doi:10.1016/j.biopsych.2017.09.006
Lovibond, P. F., & Lovibond, S. H. (1995). The structure of negative emotional states: comparison of the Depression Anxiety Stress Scales (DASS) with the Beck Depression and Anxiety Inventories. Behav Res Ther, 33(3), 335-343. doi:10.1016/0005-7967(94)00075-u
Pomponio, R., Erus, G., Habes, M., Doshi, J., Srinivasan, D., Mamourian, E., . . . Davatzikos, C. (2020). Harmonization of large MRI datasets for the analysis of brain imaging patterns throughout the lifespan. Neuroimage, 208. doi:ARTN 116450
10.1016/j.neuroimage.2019.116450
Radloff, L. S. (1977). The CES-D Scale: A Self-Report Depression Scale for Research in the General Population. Applied Psychological Measurement, 1(385-401). doi:10.1177/014662167700100306
Radua, J., Vieta, E., Shinohara, R., Kochunov, P., Quide, Y., Green, M. J., . . . collaborators, E. C. (2020). Increased power by harmonizing structural MRI site differences with the ComBat batch adjustment method in ENIGMA. Neuroimage, 218, 116956. doi:10.1016/j.neuroimage.2020.116956
Rinne-Albers, M. A., Boateng, C. P., van der Werff, S. J., Lamers-Winkelman, F., Rombouts, S. A., Vermeiren, R. R., & van der Wee, N. J. (2020). Preserved cortical thickness, surface area and volume in adolescents with PTSD after childhood sexual abuse. Scientific Reports, 10(1). doi:ARTN 3266
10.1038/s41598-020-60256-3
Santhanam, P., Wilson, S. H., Mulatya, C., Oakes, T. R., & Weaver, L. K. (2019). Age-Accelerated Reduction in Cortical Surface Area in United States Service Members and Veterans with Mild Traumatic Brain Injury and Post-Traumatic Stress Disorder. Journal of Neurotrauma, 36(20), 2922-2929. doi:10.1089/neu.2018.6242
Santhanam, P., Wilson, S. H., Oakes, T. R., & Weaver, L. K. (2019). Accelerated age-related cortical thinning in mild traumatic brain injury. Brain and Behavior, 9(1). doi:ARTN e01161
10.1002/brb3.1161
Savjani, R. R., Taylor, B. A., Acion, L., Wilde, E. A., & Jorge, R. E. (2017). Accelerated Changes in Cortical Thickness Measurements with Age in Military Service Members with Traumatic Brain Injury. Journal of Neurotrauma, 34(22), 3107-3116. doi:10.1089/neu.2017.5022
Shalev, A., Liberzon, I., & Marmar, C. (2017). Post-Traumatic Stress Disorder. New England Journal of Medicine, 376(25), 2459-2469. doi:10.1056/NEJMra1612499
Sripada, R. K., King, A. P., Welsh, R. C., Garfinkel, S. N., Wang, X., Sripada, C. S., & Liberzon, I. (2012). Neural Dysregulation in Posttraumatic Stress Disorder: Evidence for Disrupted Equilibrium Between Salience and Default Mode Brain Networks. Psychosomatic Medicine, 74(9), 904-911. doi:10.1097/PSY.0b013e318273bf33
Thompson, P. M., Stein, J. L., Medland, S. E., Hibar, D. P., Vasquez, A. A., Renteria, M. E., . . . Grp, S. Y. S. S. (2014). The ENIGMA Consortium: large-scale collaborative analyses of neuroimaging and genetic data. Brain Imaging and Behavior, 8(2), 153-182. doi:10.1007/s11682-013-9269-5
van Rooij, D., Anagnostou, E., Arango, C., Auzias, G., Behrmann, M., Busatto, G. F., . . . Buitelaar, J. K. (2018). Cortical and Subcortical Brain Morphometry Differences Between Patients With Autism Spectrum Disorder and Healthy Individuals Across the Lifespan: Results From the ENIGMA ASD Working Group. Am J Psychiatry, 175(4), 359-369. doi:10.1176/appi.ajp.2017.17010100
Weathers, F. W., Blake, D. D., Schnurr, P. P., Kaloupek, D. G., Marx, B. P., & Keane, T. M. (2013). The PTSD Checklist for DSM-5 (PCL-5).
Weathers, F. W., Bovin, M. J., Lee, D. J., Sloan, D. M., Schnurr, P. P., Kaloupek, D. G., . . . Marx, B. P. (2018). The Clinician-Administered PTSD Scale for DSM-5 (CAPS-5): Development and Initial Psychometric Evaluation in Military Veterans. Psychological Assessment, 30(3), 383-395. doi:10.1037/pas0000486
Weathers, F. W., Keane, T. M., & Davidson, J. R. T. (2001). Clinician-administered PTSD scale: A review of the first ten years of research. Depression and Anxiety, 13(3), 132-156. doi:DOI 10.1002/da.1029
Wolf, E. J., Logue, M. W., Hayes, J. P., Sadeh, N., Schichman, S. A., Stone, A., . . . Miller, M. W. (2016). Accelerated DNA methylation age: Associations with PTSD and neural integrity. Psychoneuroendocrinology, 63, 155-162. doi:10.1016/j.psyneuen.2015.09.020
Wolf, E. J., Sadeh, N., Leritz, E. C., Logue, M. W., Stoop, T. B., McGlinchey, R., . . . Miller, M. W. (2016). Posttraumatic Stress Disorder as a Catalyst for the Association Between Metabolic Syndrome and Reduced Cortical Thickness. Biological Psychiatry, 80(5), 363-371. doi:10.1016/j.biopsych.2015.11.023
Wrocklage, K. M., Averill, L. A., Scottd, J. C., Averill, C. L., Schweinsburg, B., Trejo, M., . . . Abdallah, C. G. (2017). Cortical thickness reduction in combat exposed US veterans with and without PTSD. European Neuropsychopharmacology, 27(5), 515-525. doi:10.1016/j.euroneuro.2017.02.010
Yehuda, R., Hoge, C. W., McFarlane, A. C., Vermetten, E., Lanius, R., Nievergelt, C. M., . . . Hyman, S. E. (2015). Post-traumatic stress disorder. Nature Reviews Disease Primers, 1. doi:ARTN 15057
10.1038/nrdp.2015.57
Yu, M. C., Linn, K. A., Cook, P. A., Phillips, M. L., McInnis, M., Fava, M., . . . Sheline, Y. I. (2018). Statistical harmonization corrects site effects in functional connectivity measurements from multi-site fMRI data. Human Brain Mapping, 39(11), 4213-4227. doi:10.1002/hbm.24241
Zandvakili, A., Barredo, J., Swearingen, H. R., Aiken, E. M., Berlow, Y. A., Greenberg, B. D., . . . Philip, N. S. (2020). Mapping PTSD symptoms to brain networks: a machine learning study. Translational Psychiatry, 10(1). doi:10.1038/s41398-020-00879-2