Altered White Matter Microstructural Organization in Post-Traumatic Stress Disorder across 3,047 Adults: Results from the PGC-ENIGMA PTSD Consortium

A growing number of studies have examined alterations in white matter organization in people with posttraumatic stress disorder (PTSD) using diffusion MRI (dMRI), but the results have been mixed, which may be partially due to relatively small sample sizes among studies. Altered structural connectivity may be both a neurobiological vulnerability for, and a result of, PTSD. In an effort to find reliable effects, we present a multi-cohort analysis of dMRI metrics across 3,047 individuals from 28 cohorts currently participating in the PGC-ENIGMA PTSD working group (a joint partnership between the Psychiatric Genomics Consortium and the Enhancing NeuroImaging Genetics through Meta-Analysis consortium). Comparing regional white matter metrics across the full brain in 1,426 individuals with PTSD and 1,621 controls (2174 males/873 females) between ages 18–83, 92% of whom were trauma-exposed, we report associations between PTSD and disrupted white matter organization measured by lower fractional anisotropy (FA) in the tapetum region of the corpus callosum (Cohen’s d=−0.11, p=0.0055). The tapetum connects the left and right hippocampus, structures for which structure and function have been consistently implicated in PTSD. Results remained significant/similar after accounting for the effects of multiple potentially confounding variables: childhood trauma exposure, comorbid depression, history of traumatic brain injury, current alcohol abuse or dependence, and current use of psychotropic medications. Our results show that PTSD may be associated with alterations in the broader hippocampal network.


Introduction
Posttraumatic stress disorder (PTSD) is a debilitating mental health condition with a lifetime prevalence varying globally between 1-9% 1 , with higher rates in women. Rates of PTSD are higher in populations exposed to greater levels of trauma, such as combat veterans 2 and civilians in conflict zones 3 . In addition to trauma type, genetics, and other sociological, psychological, and biological factors, individual differences in brain structure and function may explain vulnerability to developing PTSD following exposure to trauma, may result from trauma, or may be exacerbated by PTSD 4 . There is a lack of mechanistic evidence on the effects of stress and trauma on white matter structure. Exposure to trauma could lead to white matter damage, as excessive glucocorticoid levels can be neurotoxic and can impact myelination 5,6 . Diffusion MRI (dMRI) is able to model white matter tracts and assess microstructural organization 7 . Fractional anisotropy (FA) is the most commonly used metric of microstructural organization, reflecting the degree to which water is diffusing along the axon (axially) as compared with across it (radially). Greater FA can reflect higher outlier detection. Details on ENIGMA-DTI methods 39 , inclusion/exclusion criteria, and clinical information may be found in Supplementary Note 1, Supplementary Table 1, and Supplementary Note 2, respectively.

Image Acquisition and Processing
The acquisition parameters for each cohort are provided in Supplementary Table 2.
Preprocessing, including eddy current correction, echo-planar imaging-induced distortion correction and tensor fitting, was carried out at each site. Recommended protocols and quality control procedures are available on the ENIGMA-DTI and NITRC (Neuroimaging Informatics Tools and Resources Clearinghouse) webpages. Harmonization of preprocessing schemes was not enforced across sites to accommodate site-and acquisition-specific pipelines. Once tensors were estimated, they were mapped to the ENIGMA DTI template and projected onto the ENIGMA-DTI template and were averaged within ROIs (http:// enigma.ini.usc.edu/protocols/dti-protocols/). Further details and ROI abbreviations can be seen in Supplementary Note 1. In a subset of sites we extracted motion parameters from the Eddy current correction step to determine if motion played a role in our case control findings. We examined rotation and translation averaged across the X, Y, and Z axes and found a marginally significant difference between PTSD and control groups in average rotation (d=0.12 p=0.030). Follow-up analyses including motion parameters as covariates yielded results that were consistent with our main analyses (for details see Supplementary Note 1 and Supplementary Figure 2).

Statistical Analysis
For each cohort/study, a linear model was fit using the ppcor and matrixStats packages in R 3.1.3, with the ROI FA as the response variable and PTSD and covariates as predictors. For cohorts/studies including more than one data collection site, site was included as a fixed dummy variable in the site-level analysis. As in prior ENIGMA disease working group meta-analyses 39 , a random-effects inverse-variance weighted meta-analysis was conducted at a central coordinating site (the University of Southern California Imaging Genetics Center) in R (metafor package, version 1.99-118 http://www.metafor-project.org/) to combine individual cohort estimated effect sizes. Cohen's d for the main effect of group and unstandardized β coefficients (regression parameters) for continuous predictors were computed with 95% confidence intervals. We used the Cohen's d calculation that accounts for covariates in the fixed effects model, using the following equation: Heterogeneity scores (I 2 ) for each test were computed, indicating the percent of total variance in effect size explained by heterogeneity across cohorts. Bilaterally averaged FA was the primary imaging measure, with corresponding MD, RD, and AD examined post hoc when FA was significant for an effect of diagnosis. Lateralized ROIs were examined post hoc when a significant association was found with the bilateral average. The corticospinal tract was not analyzed as it has poor reliability 38 . The average correlation in FA between all ROIs was r=0.61. A Bonferroni correction is overly conservative when there are correlations between the multiple dependent measures being tested. Therefore, we calculated the effective number of independent tests based on the observed correlation structure between the alternate responses using two different methods. The equation proposed by Nyholt 45 yielded an estimated 13.3529 effective tests (Meff), giving a significance threshold of p<0.05/13.3529=0.0037. The equation of Li and Ji 46 yielded Veff=9, giving a significance threshold of p<0.05/9=0.0057. Recent empirical work based on simulations concludes that when the number of effective tests were known, the Nyholt calculation is less accurate than the method of Li and Ji because it overestimates the number of effective variables 47 . Given the lack of consensus in the field along with initial evidence that the Li and Ji threshold is more accurate, we adopted this method but also note whether results survived the Nyholt threshold.

Code Availability
All analyses were conducted using generalizable scripts available on the ENIGMA-GitHub https://github.com/ENIGMA-git/ENIGMA/tree/master/WorkingGroups/ EffectSize_and_GLM. Individual ROI level data was shared with the central site and processed using a set of R scripts with regressions customized for the current PGC ENIGMA-PTSD DTI analysis, publicly available on a set of Google Spreadsheet configuration files.
Non-linear age term-We first conducted analyses to examine whether a non-linear age term should be included in statistical models along with age and sex, as age has a non-linear effect on FA 48 . As this analysis did reveal a significant effect of non-linear age above and beyond linear age, age 2 was included in all subsequent analyses. Demographic details are shown in Supplementary Figure 3.
Primary -group comparison-We compared PTSD cases to all controls (both trauma exposed and unexposed), PTSD cases to trauma-exposed controls only, and trauma-exposed to trauma un-exposed controls. Given the significant impact of trauma-exposure on brain structure 49 , we compared trauma-exposed to un-exposed controls to examine whether potential differences may be attributed to trauma exposure rather than PTSD. We repeated analyses separately in sites acquiring 30 or fewer diffusion directions, sites acquiring more than 30 directions, sites using isotropic voxel dimensions, and sites using anisotropic voxel dimensions. These are discussed in further detail in Supplementary Note 3, Supplementary Table 3, and Supplementary Figure 4. Secondary -subgroups: We examined PTSD associations in males and females separately, and in military and civilian samples separately. These results may be found in Supplementary Note 4. Secondary -interactions: We examined potential interactions between PTSD and age or sex. These results can be seen in Supplementary Note 5.
Secondary -additional covariates-We tested a model including ancestry, but as this was a meta-analysis and most cohorts were primarily composed of participants of white non Hispanic descent, this had a very minimal impact, and we did not include this variable as a covariate in our analysis. We examined the impact of five potentially confounding covariates on the associations of PTSD with FA -childhood trauma, depression, alcohol dependence/ abuse, traumatic brain injury (TBI, of any severity), and use of psychotropic medications. We compared the white matter microstructure of individuals with PTSD to that of controls with each covariate included individually in the model, and in the subset of sites that collected data on childhood trauma, depression, alcohol use disorders, TBI, or medication without that covariate in the model to determine whether differences in results were due to the inclusion of the covariate or the reduction in sample size. There were not enough participants with all five variables to simultaneously model these potential confounds in a single model. Details of these methods and results are provided in Supplementary Note 6. Briefly, binary variables were created for depression, TBI, and medication use. As depression was assessed using a variety of measures, we used published cut-offs to recode the data as categorical depression (see Supplementary Note 2 for more details). Alcohol use disorders and childhood trauma were coded as three-level ordinal variables based on evidence of dose-dependent effects on brain structure and clinical severity, respectively 50,51 : Alcohol use disorders: 0=no alcohol use disorder, 1=alcohol abuse, 2=alcohol dependence, as measured by the SCID or AUDIT 52 ; childhood trauma (as measured by the Childhood Trauma Questionnaire): 0=no reported childhood trauma, 1=one type of childhood trauma exposure, 2=two or more types of childhood trauma exposure; PTSD severity: To examine PTSD severity, we conducted linear regressions on CAPS-4 score in the PTSD group for sites that collected CAPS-4. We examined linear associations with CAPS-4 score covarying for childhood trauma, depression, alcohol use disorders, TBI, and medication use, and we tested associations with CAPS-4 separately in military veterans and civilians as well as males and females (see Supplementary Note 6 for more details).

Group differences
We found significantly lower FA in the PTSD group in the tapetum of the corpus callosum (d=−0.11, p=0.0055) when comparing PTSD (n=1,1377) and all controls (n=1,620). This result did not survive the more conservative but potentially less accurate Nyholt threshold. Post hoc analysis revealed a larger effect in the left than in the right tapetum (left d=−0.13, p=0.00090; right d=−0.067, p=0.042). Post hoc analysis also revealed higher RD in the tapetum in the PTSD group (bilateral d=0.09, p=0.027; left d=0.11, p=0.0038).
In the analysis comparing participants with PTSD (n=1,319) to trauma-exposed controls (n=1,498), we found lower FA in the left tapetum in PTSD (d=−0.13, p=0.0014), although the effect in bilateral tapetum did not survive correction for multiple comparisons but had comparable effect size (d=−0.10, p=0.011) (see Table 2 and Figure 1, see Figure 2 for site-specific effects, and Figure 3 for a neuroanatomical figure). PTSD participants from cohorts that only included trauma-unexposed controls were not included.
Comparing trauma-exposed (n=200) to trauma un-exposed controls (n=93) from 6 sites, we found an effect that did not withstand correction for multiple comparisons but had a comparable effect size of lower FA in exposed controls in the tapetum, splenium of corpus callosum, and fornix/stria-terminalis (d=−0.41, p=0.014; d=−0.48, p=0.0042; d=−0.36, p=0.019, respectively), along with significantly higher MD and RD in the splenium (d=0.48, p=0.0017; d=0.56, p=0.00023, respectively) and an effect that did not survive correction for multiple comparisons but had a comparable effect size of higher RD in the tapetum (d=0.32, p=0.036). Lastly, for a subset of cohorts we had data on lifetime PTSD (444 lifetime PTSD vs. 230 controls across 7 cohorts). We did not find any significant results for this comparison.

Subgroups
We examined military vs. civilian cohorts, and male vs. female participants separately. We found significantly lower FA in males in the SFO, along with associations that did not survive correction for multiple comparisons with tapetum FA in the military-only and male-only subgroups separately (see Supplementary Note 4 and Supplementary Figure 5 for more details). Results of group-by-sex and group-by-age interactions were not significant and are shown in Supplementary Note 5 and Supplementary Figure 6.

Additional Covariates
The role of potentially confounding variables on the association between PTSD and the tapetum was tested in several post hoc analyses focused on left, right, and bilateral tapetum FA (see Supplementary Figure 7). As these analyses were considered post hoc and limited to the tapetum, we used a test-wise significance threshold of p<0. 05 PTSD vs. 598 controls) did not yield any significant results, but neither did the analysis in the reduced sample without CT in the model, suggesting that the sample reduction impacted these results. To control for covariate-and cohort-dependent changes in sample size, each analysis was repeated in a smaller sample that corresponded to omitting the relevant covariate. The tapetum results remained consistent in nearly all reduced sample analyses -significant effects survived covariate adjustment and effects that disappeared (such as with childhood trauma) were also absent in the reduced sample. Thus, covariates had minimal impact beyond the reduction in sample size (see Supplementary Note 6 and Supplementary Figures 8-12). A table showing how many participants at each site had information on these potentially confounding variables may be found in Supplementary  Table 4.

PTSD Severity
Examining PTSD symptom severity in subjects assessed with CAPS-4 (N=1764 from 17 sites), we found an association, which did not survive correction for multiple comparisons, between higher CAPS-4 score and lower FA in the tapetum and SFO (tapetum: β=−9.0×10 −5 , p=0.024; SFO: β=−6.6×10 −5 , p=0.028; Figure 4). PTSD symptom severity in the PTSD group (measured by the CAPS-4, N=956 from 17 sites) was not associated with FA ( Figure 4). Subgroup analyses yielded associations that did not survive correction for multiple comparisons between CAPS-4 score and tapetum FA in the military and female cohorts and between CAPS-4 score and superior fronto-occipital fasciculus (SFO) FA in the male subgroup. Results were similar when potentially confounding variables were included, with no significant associations, although there was an effect in tapetum that did not withstand correction for multiple comparisons when psychotropic medication use was included. Detailed analyses of PTSD severity and covariates within subgroups are in Supplementary Note 6 and Supplementary Figures 13 and 14.

Discussion
We present DTI results from a multi-cohort study conducted by the PGC-ENIGMA PTSD consortium. In a meta-analysis of 3,047 participants from 28 sites, we found lower FA and higher RD in the tapetum among adults with PTSD (neuroanatomical figure -Figure 3), which remained after accounting for several potentially confounding factors. The tapetum is a major tract within the corpus callosum that serves as a conduit between right and left hippocampus. Prior studies of white matter disruption in PTSD have found alterations in other hippocampal tracts, but were generally hindered by small sample sizes leading to inconsistent findings across studies. Our results add to the existing literature in identifying structural disruptions that compromise putative hippocampal functions, which are known to play a central role in PTSD symptomatology 53,54 .
The tapetum is a small segment of the corpus callosum that connects the temporal lobes, in particular the left and right hippocampus 55 . It is one of the last corpus callosum segments to develop and experiences rapid growth around age 14, which may make it vulnerable to the effects of trauma for a longer period of time 56 . Structural and functional alterations in the hippocampus are frequently reported in PTSD, with smaller volumes 37 , decreased activation, and disrupted functional connectivity with the medial and lateral prefrontal cortices 57,58 . Disrupted functional connectivity between the left and right hippocampus has also been previously reported in PTSD 59 . Here we report microstructural evidence that structural connectivity between the left and right hippocampus may also be disrupted in PTSD.
We did not find PTSD-related effects in the cingulum, uncinate, fornix, and corpus callosum, as some prior studies of smaller, more homogeneous populations have shown [9][10][11][12][13][14]16,[18][19][20][21][23][24][25][26][27][28][29]31,32,36,60 . This may reflect a limitation in the meta-analysis of heterogeneous study samples, but this is an unlikely explanation given our low estimates of heterogeneity (I 2 ) in these tracts (< 10%). However, these effects were not entirely absent in our subgroup analyses, suggesting that they may be more specific to population, trauma type, etc. While many studies have reported that PTSD is associated with alterations in the cingulum bundle 9-13,16,18,21,23-29,31,32,36 , which has a hippocampal component, the tapetum has not yet emerged for several possible reasons. Many prior studies took an ROI approach, which limited analyses to pre-determined regions that frequently omitted the tapetum, a small region often grouped with other tracts such as the splenium or posterior thalamic radiation. Critically, in 2013, an error was uncovered in the JHU atlas used as part of the TBSS pipeline, with the uncinate incorrectly identified as the inferior fronto-occipital fasciculus, and the tapetum incorrectly identified as the uncinate 61 . Thus, the tapetum was simply not examined in prior studies, with one very recent exception showing that tapetum abnormalities are associated with lower major depressive disorder remission 62 . There is further data suggesting that the absence of tapetum findings may partly relate to methodological limitations of prior studies; in a recent DTI meta-analysis of attention deficit hyperactivity disorder (ADHD), lower tapetum FA was found to characterize ADHD, despite the fact that this specific region of the CC has not frequently emerged in ADHD studies 63 . Finally, the precise role of the tapetum in connecting the left and right hippocampus was only recently elucidated by mapping the subcortical connectome with exquisitely high-resolution mapping capable of discerning the intermingling of tapetum and other corpus callosum fibers 55 .
Childhood trauma is the greatest single risk factor for future vulnerability to PTSD 64 ; numerous studies show significant alterations in brain structure and function in individuals who experience significant early life stress 51,65 . Some of these alterations likely contribute to a higher risk for psychopathology, but childhood trauma exposure alone did not explain the association between tapetum white matter disruption and PTSD that we report here. While it is likely that childhood trauma exposure plays some role in tapetum organization, these analyses indicate that our PTSD results are not solely due to confounding effects between PTSD and exposure to childhood trauma. Depression is frequently comorbid with PTSD 66 and is associated with disrupted white matter organization, although the affected tracts are broadly distributed 67 . Accounting for depression in group comparisons did not significantly alter our results, suggesting that tapetum white matter disruption is not confounded by depression. The UK Biobank sample revealed FA reductions associated with depression in commissural and association fibers, thalamic radiations, left superior longitudinal fasciculus (SLF), superior and forceps major 68 . The ENIGMA MDD Working Group reports alterations in the corpus callosum and hippocampal cingulum but not in the tapetum (labeled in that paper as the uncinate, p=0.172) 69 . Particularly in military populations, which formed the majority of our sample, PTSD is often comorbid with traumatic brain injury (TBI) 70 . White matter is particularly vulnerable to TBI, which produces stretching and shearing of axons and altered neurometabolism 71 . Accounting for TBI also had some effect on our results but the left tapetum remained significant in post hoc analyses, indicating that TBI is not associated with white matter damage in the tapetum despite numerous reports mild TBI is associated with white matter damage generally. Preliminary results from the ENIGMA Military Brain Injury group have not shown significant differences in the tapetum 72 . Psychotropic medications are another potential confound, given their neurotrophic and neuroprotective effects 73 . The result in the tapetum persisted after covarying for psychotropic medication, indicating that our findings are unlikely to be explained by medication. Lastly, PTSD can be comorbid with alcohol use disorders, which have a poorer clinical prognosis 74,75 . Alcohol use disorders have been associated with significant changes in white matter organization 76,77 but did not influence the present results. Although none of these potentially confounding variables fully explained our reported results, it is possible that altered tapetum microstructure could represent a more general response to trauma or stress. As shown in Figure 1, altered tapetum FA is present in trauma-exposed controls compared to unexposed controls, although the sample size for this comparison was around 10% of the total sample size. The post hoc analysis of the tapetum was still significant comparing PTSD to trauma-exposed controls, suggesting a dose response, with some alterations related to trauma that does not lead to psychopathology but more in PTSD. Trauma is one form of severe stress, but we do not know how more minor stress may affect tapetum organization. Smaller hippocampal volume has been reported across disorders, suggesting that it may be a non-specific marker of disease 78 .
We found an association that did not withstand correction for multiple comparisons of PTSD with FA in the tapetum in male and military subgroups separately. Although results were non-significant in female or civilian subgroups, the effect size was slightly larger and in the same direction. The female and civilian subgroups were smaller and therefore the analyses had lower power than in male and military subgroups. Most prior dMRI studies in civilians report lower FA in PTSD 9,10,[12][13][14]16,17,[21][22][23][24][25][30][31][32]60 . Studies of military cohorts have been mixed, reporting higher FA [26][27][28][29]36 , lower FA 8,15,[18][19][20] , and null results [33][34][35] . This discrepancy may be due to differences in age, chronicity, and type of trauma exposure, although military personnel often also experience civilian trauma. Combat-related PTSD is often comorbid with TBI, which is also associated with white matter disruption, constituting a potentially confounding factor for studies 72 .
In the absence of longitudinal data, our analysis cannot make causal inferences nor make conclusions about the direction of the relationship between PTSD and tapetum white matter organization. Disrupted white matter of the tapetum may represent a vulnerability that predates the onset of PTSD, or a pathological response to trauma. In twins discordant for exposure to combat stress, the unexposed twins of combat veterans with PTSD have smaller hippocampal volume than the unexposed twins of combat veterans without PTSD 79 . Individuals with two risk alleles of the FKBP5 gene have demonstrated lower cingulum FA above and beyond the association of cingulum FA with PTSD 11,60 . These studies suggest that heritable differences in brain structure may influence risk of developing PTSD. Evidence that alterations are caused by PTSD was observed in Israeli Defense Force recruits with reduced structural connectivity between the hippocampus and ventromedial prefrontal cortex, but only after exposure to military stress 80 . With the varying developmental trajectories of brain structure, function, and connectivity, along with the varying distribution of stress hormone receptors in the brain, the complex question of vulnerability vs. consequence will require prospective longitudinal neuroimaging studies. Some evidence indicates that high FA is a marker of resilience to the effects of stress 81,82 . A putative marker of resilience is the ability to attenuate stress-induced increases in corticotropin-releasing hormone and glucocorticoids through an elaborate negative feedback system, and to modulate the expression of brain-derived neurotrophic factor (BDNF) 83,84 . BDNF has myriad functions including supporting neuronal differentiation, maturation, and survival 85,86 . In particular, hippocampal BDNF is implicated in the development of neural circuits that promote stress adaptations 83 . These stress adaptation circuits involve white matter in the fornix and other fronto-limbic connections 87 .

Limitations
Our study has several limitations. One limitation of TBSS studies is the inability to fully attribute results to particular fiber bundles, especially given that the tapetum is a relatively small structure. Future studies may benefit by using tractography to more reliably identify the affected bundles, but this is difficult across sites. Second, not all participants classified as PTSD received clinician-administered interview (such as the CAPS) to confirm diagnoses. Third, we could not reliably measure chronicity across different cohorts. Other variables that we could not examine given the heterogeneity across sites include treatment effects, symptom clusters, trauma types, and lifetime as opposed to current PTSD diagnosis. Given that effects have not been previously reported in the tapetum, it would be helpful to know if altered tapetum organization is linked to particular symptoms or trauma type and we encourage researchers to examine this possibility in their datasets. Although we analyzed data from over 3,000 participants, we may have been underpowered to examine group-by sex interactions, as 55% of our sample came from cohorts including only males or only females or samples that were >90% male. Diffusion metrics are not scanner invariant, and can vary even in scanners of the same model. For this reason, we were limited to a meta-analytic approach which may have lower power than mega-analysis. However, studies including both meta-and mega-analyses of brain volume in other ENIGMA groups have found minimal differences 88,89 . Lastly, our cross-sectional design did not allow us to investigate to what degree WM differences predate the onset of PTSD.
Future studies should further investigate the tapetum using high spatial and angular resolution tractography to replicate our findings. Future and existing studies with more in-depth phenotyping than was possible here could examine how alterations in the tapetum vary with trauma type, chronicity, treatment, and whether they are associated with specific symptom clusters. The current study excluded pediatric cases, so additional research on white matter disruption in pediatric trauma and PTSD is warranted. Lastly, while we considered comorbidities as potential confounding variables, we did not examine their association with dMRI metrics. Future collaborations with the ENIGMA Brain Injury, MDD, and Addiction working groups will provide opportunities to separate general neuroimaging biomarkers of psychopathology and disorder-specific effects.

Conclusions
Here we presented results from the PGC-ENIGMA PTSD working group, reporting poorer white matter organization in the tapetum in individuals currently suffering from PTSD. We present the largest DTI study in PTSD to date and the first to use harmonized image processing across sites, increasing our power to detect subtle effects. While future studies need to confirm the involvement of the tapetum specifically, our results add to the existing literature implicating the hippocampus and associated white matter connections as neural markers of structural disruptions in PTSD.  Results from the group comparisons. (A) Results comparing PTSD to all controls; (B) results comparing PTSD to trauma-exposed controls only; (C) results comparing trauma exposed to unexposed participants; (D) lifetime PTSD compared to controls. Cohen's d statistics are shown across all bilateral and midline ROIs and average FA, with bars indicating the 95% confidence interval. The ROI abbreviations are explained in Supplementary Note 1. As PTSD was coded "1" and control "0", negative statistics indicate lower FA in PTSD. Total N is listed for each comparison. Dark orange bars indicate significance (p<0.0057) and light orange bars indicate results that did not withstand correction for multiple comparisons (0.05>p>0.0057). Error bars are 95% CI. Site effects for tapetum result. Forest plot shows the effect sizes (Cohen's d) for each of the 25 cohorts, scaled by sample size, with bars for 95% CI. The effect size and 95% CI of the meta-analysis is shown at the bottom of the figure, along with effect sizes and 95% CI for the subgroup analyses of the military cohorts, civilian cohorts, male cohorts, and female cohorts. Tapetum displayed on the ENIGMA template FA. The skeleton is shown in red, the left tapetum (green) and right tapetum (blue) ROIs are displayed. Left in image is right in brain. Linear association with CAPS-4 across the whole sample (left) and within the PTSD cases only (right). Meta-regression unstandardized β statistics are shown across all bilateral and midline ROIs and average FA, with bars indicating the 95% confidence interval. The ROI abbreviations are explained in Supplementary Note 1. Light orange bars indicate results that did not withstand correction for multiple comparisons (0.05>p>0.0057). Error bars are 95% CI. Dennis    Results from the group comparisons. (A) Results comparing PTSD to all controls, (B) results comparing PTSD to trauma-exposed controls only. Cohen's d values, uncorrected p-values, the 95% confidence interval for the d statistic, and the I 2 (heterogeneity) are shown for the group comparisons. Bolded results are significant when corrected for multiple comparisons, italicized results are marginally significant (based on the Li and Ji adjusted Bonferroni correction).

Site
(A) PTSD vs all controls (B) PTSD vs trauma-exposed controls