Search

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 (LME INT), (2) LME that models both site-specific random intercepts and age-related random slopes (LME INT+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.

Background: Posttraumatic stress disorder (PTSD) is accompanied by disrupted cortical neuroanatomy. We investigated alteration in covariance of structural networks associated with PTSD in regions that demonstrate the case-control differences in cortical thickness (CT) and surface area (SA). Methods: Neuroimaging and clinical data were aggregated from 29 research sites in >1300 PTSD cases and >2000 trauma-exposed control subjects (ages 6.2–85.2 years) by the ENIGMA-PGC (Enhancing Neuro Imaging Genetics through Meta Analysis–Psychiatric Genomics Consortium) PTSD working group. Cortical regions in the network were rank ordered by the effect size of PTSD-related cortical differences in CT and SA. The top-n (n = 2–148) regions with the largest effect size for PTSD > non-PTSD formed hypertrophic networks, the largest effect size for PTSD < non-PTSD formed atrophic networks, and the smallest effect size of between-group differences formed stable networks. The mean structural covariance (SC) of a given n-region network was the average of all positive pairwise correlations and was compared with the mean SC of 5000 randomly generated n-region networks. Results: Patients with PTSD, relative to non-PTSD control subjects, exhibited lower mean SC in CT-based and SA-based atrophic networks. Comorbid depression, sex, and age modulated covariance differences of PTSD-related structural networks. Conclusions: Covariance of structural networks based on CT and cortical SA are affected by PTSD and further modulated by comorbid depression, sex, and age. The SC networks that are perturbed in PTSD comport with converging evidence from resting-state functional connectivity networks and networks affected by inflammatory processes and stress hormones in PTSD.

BackgroundProlonged grief disorder (PGD) has recently been recognized as a separate psychiatric diagnosis, despite controversy over the extent to which it is distinctive from posttraumatic stress disorder (PTSD) and major depressive disorder (MDD).MethodsThis study investigated distinctive neural processes underpinning emotion processing in participants with PGD, PTSD, and MDD with functional magnetic resonance study of 117 participants that included PGD (n = 21), PTSD (n = 45), MDD (n = 26), and bereaved controls (BC) (n = 25). Neural responses were measured across the brain while sad, happy, or neutral faces were presented at both supraliminal and subliminal levels.ResultsPGD had greater activation in the pregenual anterior cingulate cortex (pgACC), bilateral insula, bilateral dorsolateral prefrontal cortices and right caudate and also greater pgACC–right pallidum connectivity relative to BC during subliminal processing of happy faces. PGD was distinct relative to both PTSD and MDD groups with greater recruitment of the medial orbitofrontal cortex during supraliminal processing of sad faces. PGD were also distinct relative to MDD (but not PTSD) with greater activation in the left amygdala, caudate, and putamen during subliminal presentation of sad faces. There was no distinction between PGD, PTSD, and MDD during processing of happy faces.ConclusionsThese results provide initial evidence of distinct neural profiles of PGD relative to related psychopathological conditions, and highlight activation of neural regions implicated in reward networks. This pattern of findings validates current models of PGD that emphasize the roles of yearning and appetitive processes in PGD.