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.

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,343 cases with posttraumatic stress disorder (PTSD) (6.2-81.8 years old) and 2,067 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 were more sensitive to the detection of significant case-control differences in regional cortical thickness (X2(3) = 34.339, p < 0.001), and case-control differences in age-related cortical thinning (X2(3) = 15.128, p = 0.002). Specifically, ComBat-GAM led to larger effect size estimates of cortical thickness reductions (corrected p-values < 0.001), smaller age-appropriate declines (corrected p-values < 0.001), and lower female to male contrast (corrected p-values < 0.001) in cases compared to controls relative to other harmonization methods. Harmonization with ComBat-GAM also led to greater estimates of age-related declines in cortical thickness (corrected p-values < 0.001) in both cases and controls compared to other harmonization methods. Our results support the use of ComBat-GAM for harmonizing cortical thickness data aggregated from multiple sites and scanners to minimize confounds and increase statistical power.

Posttraumatic stress disorder (PTSD) is a complex psychiatric condition that has generated much attention in the neuroimaging literature. A neurocircuitry model supporting fronto-limbic dysfunction as a major player in facilitating clinical symptoms of PTSD is well-characterized; however, recent literature suggests that network-based approaches may provide additional insight into neural dysfunction in PTSD. Our analysis uses resting-state neuroimaging scans of 1063 adults from the PGC-ENIGMA PTSD Consortium to investigate a network-based model of functional connectivity in PTSD. With a novel, resolution limit-free community detection approach, 16 communities corresponding to functionally meaningful networks were detected with high quality. After group-level community detection, participants were classified into three groups (PTSD, n=418, trauma-exposed controls without PTSD, n=434, and non-trauma exposed healthy controls, n=211). Individual network connectivity metrics were calculated, including whole-brain, default mode network, and central executive network participation coefficient and connectivity strength. Linear mixed effects models revealed group differences in the whole-brain, default mode, and central executive network participation coefficient and connectivity strength such that individuals with PTSD demonstrated overall greater values. We also described sex differences such that males demonstrate greater whole-brain participation coefficient vs. females and females demonstrate greater default mode network connectivity strength vs. males. Our results suggest that PTSD in adults is associated with reduced specialization and enhanced inter-module communication throughout the brain network, which may contribute to inefficient information processing and poor emotional regulation. This study presents a novel use of resolution limit-free community detection in a large PTSD sample, revealing robust differences in resting-state network topology.

BACKGROUND: Depersonalization/derealization disorder (DPD) is a chronic and distressing condition characterized by detachment from oneself and/or the external world. Neuroimaging studies have associated DPD with structural and functional alterations in a variety of distinct brain regions. Such local neuronal changes might be mediated by altered interregional white matter connections. However, to our knowledge, no research on network characteristics in this patient population exists to date.METHODS: We explored the structural connectome in 23 individuals with DPD and 23 matched, healthy controls by applying graph theory to diffusion tensor imaging data. Mean interregional fractional anisotropy (FA) was used to define the network weights. Group differences were assessed using network-based statistics and a link-based controlling procedure.RESULTS: Our main finding refers to lower FA values within left temporal and right temporoparietal regions in individuals with DPD than in healthy controls when using a link-based controlling procedure. These links were also associated with dissociative symptom severity and could not be explained by anxiety or depression scores. Using network-based statistics, no significant results emerged. However, we found a trend for 1 subnetwork that may support the model of frontolimbic dysbalance suggested to underlie DPD symptomatology.LIMITATIONS: To ensure ecological validity, patients with certain comorbidities or psychotropic medication were included in the study. Confirmatory replications are necessary to corroborate the results of this explorative investigation.CONCLUSION: In patients with DPD, the structural connectivity between brain regions crucial for multimodal integration and emotion regulation may be altered. Aberrations in fibre tract communication seem to be not solely a secondary effect of local grey matter volume loss, but may present a primary pathophysiology in patients with DPD.

BackgroundPTSD and depression commonly co-occur and have been associated with smaller hippocampal volumes compared to healthy and trauma-exposed controls. However, the hippocampus is heterogeneous, with subregions that may be uniquely affected in individuals with PTSD and depression.MethodsWe used random effects regressions and a harmonized neuroimaging protocol based on FreeSurfer (v6.0) to identify sub-structural hippocampal markers of current PTSD (C-PTSD), depression, and the interaction of these conditions across 31 cohorts worldwide (N=3,115;Mage=38.9±13.9 years). Secondary analyses tested these associations by sex and after modeling the simultaneous effects of remitted PTSD, childhood trauma, mild traumatic brain injury, and alcohol use disorder.ResultsA significant negative main effect of depression (n=800, vs. no depression, n=1456) was observed in the hippocampal tail (ß=−0.13) and CA1 (ß=−0.09) after adjusting for covariates and multiple testing (adjusted p’s (q)=0.028). A main effect of C-PTSD (n=1042 vs. control, n=1359) was not significant, but an interaction between C-PTSD and depression was significant in the CA1 (ß=−0.24, q=0.044). Pairwise comparisons revealed significantly smaller CA1 volumes in individuals with C-PTSD+Depression than controls (ß=−0.12, q=0.012), C-PTSD-only (ß=−0.17, q=0.001), and Depression-only (ß=−0.18, q=0.023). Follow-up analyses revealed sex effects in the hippocampal tail of depressed females, and an interaction effect of C-PTSD and depression in the fimbria of males.ConclusionsCollectively our results suggest that depression is a stronger predictor of hippocampal volumetry than PTSD, particularly in the CA1, and provide compelling evidence of more pronounced hippocampal phenotypes in comorbid PTSD and depression compared to either condition alone.