Structural brain correlates of childhood trauma with replication across two large, independent community-based samples ====================================================================================================================== * Rebecca A. Madden * Kimberley Atkinson * Xueyi Shen * Claire Green * Robert F. Hillary * Emma Hawkins * Anca-Larisa Sandu * Gordon Waiter * Christopher McNeil * Mathew Harris * Archie Campbell * David Porteous * Jennifer A. Macfarlane * Alison Murray * Douglas Steele * Liana Romaniuk * Stephen M. Lawrie * Andrew M. McIntosh * Heather C. Whalley ## Abstract **Introduction** Childhood trauma and adversity are common across societies and have strong associations with physical and psychiatric morbidity throughout the life-course. One mechanism through which childhood trauma may predispose individuals to poor psychiatric outcomes, such as raised risk of lifetime depression, could be via associations with brain structure. This study aimed to elucidate the associations between childhood trauma scores and brain structure across two large, independent community cohorts. **Methods** The two samples comprised (i) a subsample of individuals from Generation Scotland with imaging and in-depth phenotyping, including the CTQ-28 (n=1,024); and (ii) individuals from UK Biobank with imaging and a modified summary CTQ measure (n=27,202). This comprised n=28,226 for mega-analysis. Scans were processed using *FreeSurfer* image processing software, providing cortical and subcortical as well as global brain metrics. Regression models were used to determine associations between these metrics and childhood trauma measures. Associations between childhood trauma measures and psychiatric phenotypes were also explored. **Results** Childhood trauma measures associated with lifetime risk of depression diagnosis with similar ORs across cohorts (OR 1.06, 1.23 GS and UKB respectively), which also related to earlier onset and more recurrent course within both samples. There was also evidence for associations between childhood trauma measures and a range of brain structures. Replicated findings included reduced global brain volumes, reduced cortical surface area but not thickness, with highest effects at mega-analysis seen in the frontal (β=-0.0385, SE=0.0048, p(FDR)=5.43×10−15) and parietal lobes (β=-0.0387, SE=0.005, p(FDR)=1.56×10−14). At a regional level, one subcortical regional volume in particular – the ventral diencephalon (VDc) – displayed significant associations with childhood trauma measures across the two cohorts and at mega-analysis (β=-0.0232, SE=0.0039, p(FDR)=2.91×10−8). There was also evidence for associations with reduced hippocampus, thalamus, and nucleus accumbens volumes, however these were not as consistent across cohorts. **Discussion** There was strong evidence for associations between childhood trauma and reduced global and regional brain volumes across cohorts. In particular, the presence of an association between childhood trauma and the volume of the VDc (which includes the hypothalamic area), with replication, provides further evidence of the importance of neuroendocrine stress response pathways in links between early life stress and clinical outcomes. ## Introduction Adverse childhood experiences (ACEs) are reported to affect around 50% of the UK population (1, 2), with similar proportions worldwide (3). The classification of ACEs varies, but typically include abuse and neglect, witnessing domestic or neighbourhood violence, family substance misuse, and parental divorce. ACEs have been found to act in a dose-dependent manner to raise the risk of a broad range of adverse outcomes in adulthood, including poor psychiatric outcomes (4). The childhood trauma questionnaire (CTQ) describes a narrower range of early life experiences, focussing specifically on abuse and neglect, mainly in the home environment (5). It has been repeatedly demonstrated that experiences that fall under this definition of ‘childhood trauma’ (CT) have stronger associations with poor psychiatric outcomes than other types of childhood adversity (CA). It has been shown, for instance, that interpersonal trauma (adolescent sexual or physical abuse) is more strongly associated with lifetime SCID diagnoses than traumatic bereavement during adolescence (6); and that physical neglect, emotional neglect, and sexual abuse show stronger associations with symptom severity among adult psychiatric inpatients than other ACEs such as peer bullying (7). The biological underpinnings of the relationship between CA and psychiatric disorders are likely to be multifaceted; including chronic dysregulation of the hypothalamic-pituitary-adrenal axis (8), disrupted attachment and emotional development (9), and deviations from typical brain development (10). The latter hypothesis posits that the experience of trauma during sensitive periods of neurodevelopment may lead to structural abnormalities which predispose the adult individual to pathological symptoms (11). A great deal of previous work has focused on describing the neurostructural consequences of CA; the literature has, however, proven to be inconsistent where results fail to replicate study-to-study. The field to date has often taken a hypothesis led approach; while there can be a strong rationale for examining specific regions of interest (ROIs) in relation to CA, this *a priori* approach can also contribute to bias and inconsistencies, whereby meaningful effects in other regions may be missed. The hippocampus, for instance, has received a lot of attention due to evidence for smaller hippocampal volumes in depression (12), alongside findings from preclinical work indicating the hippocampus is particularly vulnerable to early stressful life events (13). These region-of-interest (ROI) studies generally find reduced hippocampal grey matter (GM) volume (14-20) in CA-exposed groups; however other studies report no association (21-24). A meta-analysis including 17 of these ROI studies looking specifically at the hippocampus did find evidence, however, for reduced volumes of the left, right, and total hippocampus (hedges g=-0.642, -0.616, -0.517 respectively), despite significant between-study heterogeneity, and publication bias in the latter (Egger’s t(df=33)=3.44, p=0.001)(25). The amygdala has been another focus of ROI analyses of the neurostructural effects of CA, given its role in emotion and fear processing (19, 24, 26). Two meta-analyses of amygdala-specific ROI analyses are however contradictory, one reporting reductions in amygdala volume from meta-analysis of 19 studies (25); the other finding no effect in a meta-analysis of 15 studies (16). Other commonly used ROIs have included the caudate nucleus (22, 23, 27), nucleus accumbens (21, 22), orbitofrontal cortex (18, 22), and anterior cingulate cortex (18, 22-24), with mixed findings. Despite multiple strands of evidence that traumatic events can affect the structure and function of the thalamus, this region is rarely included as an ROI in *a priori* studies of CA. After traumatic events in adulthood, structural and functional changes have been observed in the thalamus (28-30). Thalamic structural deviations after CA have also been observed (31, 32), and childhood trauma has been associated with functional hyperconnectivity of the thalamus in a transdiagnostic sample (33). In addition to these structural and functional imaging findings, endocrine work has repeatedly shown evidence of dysregulation in the hypothalamic-pituitary-adrenal (HPA) axis in CT-exposed adult populations (34-37); implicating associated thalamic structures as well. The lack of inclusion of the thalamus and hypothalamus in *a priori* structural imaging studies of CA demonstrates the potential to miss interesting regions by taking an ROI-based approach. While brain-wide (*a priori*) analyses have been conducted, they are rarer and often lack robust sample sizes. One such study revealed an association between CTQ scores and right insular surface area, (38). Meta-analyses of voxel-based morphometry (VBM) approaches to brain-wide studies of CA report reduced GM volumes in the dorsolateral prefrontal cortex, right hippocampus, and right postcentral gyrus (across 19 studies with 1095 adult subjects (25)); and GM volume reductions in temporal, frontal, and midbrain regions, and the left postcentral gyrus (across 12 studies with a total of 693 children, adolescents, and adults (39)). The ENIGMA consortium (40) has recently produced larger-sample, *a priori* analyses of the impact of CA on adult brain structure based on meta-analysed data from multiple cohort studies. In an analysis across 7 subcortical grey matter structures defined by the FreeSurfer image processing software, no significant main effects of childhood adversity were found (41). A similar approach taken in cortical thickness and surface area measures found main effects of CA score in the thickness of the banks of the superior temporal sulcus and the supramarginal gyrus, and the surface area of the middle temporal gyrus (42). They also found various interaction effects with age, sex, and MDD diagnosis. It is evident that there is limited consistency in determining the relationships between childhood trauma and adult structural imaging findings from whole-brain analyses. These inconsistencies may result from between-study variability. In particular, differences in the detail and definition of childhood trauma/adversity; in the demographics; and in the use of healthy or psychiatric populations (with varying diagnoses and severities of clinical status). With the growing opportunities provided by neuro-imaged population cohorts, large sample sizes are available within a single study protocol, enabling the analysis of neurostructural sequalae of childhood trauma more homogenous samples. In this study, therefore, we describe psychiatric phenotyping, whole-brain analysis, and regional analysis of childhood trauma (CT) in a subsample of the Generation Scotland: Scottish Family Health Study (GS) cohort with neuroimaging data. We based our initial analyses in this cohort due to its characteristically deep phenotyping, including the 28-item Childhood Trauma Questionnaire. We then sought to replicate in the larger UK Biobank cohort, which utilised a much shorter questionnaire assessing CT. In addition, we report imaging mega-analysis from the combined cohorts, with a total sample size of n=28,226. These large samples make the findings presented here well-powered, and more generalisable to the population than smaller clinically oriented samples. ## Methods The analysis comprised two independent populations: a subsample from The Generation Scotland (n=1,024) and the UK Biobank cohort (n=27,202), totalling a population of n=28,226. In both GS and UKB, freesurfer derived measures of cortical metrics and subcortical volumes were analysed. The analysis focussed initially on global cortical measures, lobar regions, and finally individual regional measures. ### The Generation Scotland cohort Generation Scotland: Scottish Family Health Study (GS hereafter) is a population-based cohort of over 24,000 individuals with in-depth phenotyping recruited between 2006-2011. A subsample of participants was re-contacted in 2015-2019 for neuroimaging and further data collection. This subsample was used as the basis for the current study, and is described in detail elsewhere (43-46). Cognitive assessments, blood sampling, physical measurements, and clinical questionnaires were also collected from GS participants, including the 28-item Childhood Trauma Questionnaire (CTQ; (5)). The subsample included n=1,153 participants with CTQ and depression phenotyping, and n=1,024 also had MRI imaging data. The 28-item CTQ is a questionnaire validated for self-report of abuse and neglect during childhood (5). The questionnaire is made up of five subscales – Emotional Abuse (EA), Emotional Neglect (EN), Physical Abuse (PA), Physical Neglect (PN), and Sexual Abuse (SA) – with five items for each subscale, scored on a five-point Likert scale rating frequency of each experience (1-5). Emotional and Physical Neglect items are reverse-scored. The minimisation and denial scale makes up the remaining three items, a scale devised to help detect under-reporting of childhood trauma (5). The analyses reported here focused mainly on the summed score of the 25 subscale items of this questionnaire, scored on a scale 25-125, which provides an index of the cumulative traumatic experience of an individual, consistent with previous literature (e.g. (47, 48)). The questionnaire can also be used to break scores down into severity categories of none-low, low-moderate, moderate-severe, and severe-extreme, for both total CTQ score and each subscale (49, 50). All analyses were performed using scaled scores for the five childhood trauma subscales captured by the CTQ (and, in UKB, CTM) measures – emotional neglect, emotional abuse, physical neglect, physical abuse, and sexual abuse – as well as scaled composite ‘abuse’ and ‘neglect’ scores created by summing scores on the abuse and neglect items, respectively. Major Depressive Disorder (MDD) diagnoses were derived from the Structured Clinical Interview for DSM-IV (SCID; (43, 51)) taken at GS baseline recruitment; a binary variable describing presence/absence of a lifetime diagnosis encompasses single-episode, chronic, post-partum onset, and depression with manic/hypomanic features. Two individuals with manic/hypomanic episodes alone were excluded from further analyses. The Quick Inventory of Depressive Symptomatology (52) was used as measure of continuous (and current) depressive symptom severity at the time of the imaging clinic visit. ### MRI Imaging in GS Imaging was conducted at two sites, Aberdeen and Dundee. In Aberdeen, brain magnetic resonance imaging (MRI) data were acquired on a 3T Philips Achieva TX-series MRI system (Philips Healthcare, Best, Netherlands) with a 32-channel phased-array head coil with a back facing mirror (software version 5.1.7; gradients with maximum amplitude 80 mT/m and maximum slew rate 100 T/m/s). In Dundee, images were acquired on a Siemens 3T Prisma-FIT (Siemens Healthineers, Erlangen, Germany) with 20 channel head and neck coil and a back facing mirror (software version VE11, gradient with max amplitude 80 mT/m and maximum slew rate 200 T/m/s). The current study uses T1 structural data, though other sequences were also acquired (53, 54). Structural measures were derived in-house from raw images using FreeSurfer version 5.3 (55-57). The 1,070 images were segmented into grey matter, white matter, and cerebrospinal fluid. Grey matter was further segmented into cortical and subcortical grey matter. The cortex was divided into 34 regions per hemisphere according to the Desikan-Killany atlas (58). Visual quality control was undertaken for exclusion of participants with major output errors in segmentation or parcellation. There were 424 subjects for whom a degree of manual editing was required to make small corrections to the images. Participants (n=12) were excluded if errors could not be corrected in this way. For the 68 bilateral cortical regions, measures of mean thickness, surface area, and volume were calculated. For 8 bilateral subcortical grey matter structures, volume alone was calculated. Coordinates of head position within the scanner were included as covariates. ### The UK Biobank cohort UK Biobank (UKB) is a large, UK-based population cohort of over n=500,000 adults recruited between 2006-2010 (59). Data was collected over several instances of postal questionnaires and clinic visits. MRI imaging was performed on n=100,000, which was acquired in-clinic during the second assessment instance (59, 60); at the time of writing, MRI data were available for n=41,985 (61). An online follow-up questionnaire was disseminated to participants after completion of the imaging appointment, comprising the Mental Health Questionnaire (MHQ; which included childhood trauma questions). Five items relating to childhood trauma were included in the MHQ, derived from the 28-item CTQ with each item scored on a five-point likert scale (0-4). Hence only one item related to each of the five subscales (as described in the fuller 28-item CTQ) was included in this cohort, making this data less sensitive than the GS and necessitating a focus throughout the study on total scores (with subscale results presented in Appendix 3). For the purposes of this study, childhood trauma items were used on a 0-20 scale of the summed total of responses to these items. Any individuals with incomplete responses to the childhood trauma items were excluded from analyses, leaving a maximum sample of n=153,650 for this report, of which n=27,202 had MRI imaging data available. MDD diagnoses were derived from ICD10 diagnostic criteria assessed at the initial assessment centre visit (62), and CIDI-SF (63) derived questionnaire in the MHQ at online follow up. A modified form of the Generalised Anxiety Disorder Assessment (GAD-7; (64)) was delivered in the MHQ, as well, to indicate recent anxiety symptoms. The 4-item Patient Health Questionnaire (PHQ-4 (62)) was used to take current mood level at the assessment centre appointment. ### MRI Scanning in UKB Imaging for the 40,000 data release in January 2020 was conducted at three sites: Stockport, Newcastle, and Reading. These sites are equipped with identical MRI scanners - 3T Siemens Skyra (software VD13), using the standard Siemens 32-channel head coil (65). The current study uses T1 structural scans, although other sequences were also collected. Structural measures were derived for the 40,000 data release using FreeSurfer software (61). See supplementary information p.3 ### Statistical analyses All statistical analyses were performed in R v.3.6.4. All analyses were corrected for multiple comparisons using the false-discovery rate (FDR) method. Linear and binomial regressions were used to describe phenotypic relationships between scaled childhood trauma and demographic information, as well as features of depressive symptomatology in both cohorts. Numbers for each model varied in both cohorts, as different measures had different numbers of non-completers. These models adjusted for age and sex as minimum covariates. To determine the relationships between scaled childhood trauma total and subscale scores and structural brain imaging features, linear regression models were conducted using the NLME package in R. For cortical structures, volume, surface area, and thickness measures were available; for subcortical structures volume data were available. Covariates included sex, age, age squared (to account for non-linear effects of age), scaled intracranial volume, hemisphere, scan site. For 5 summed lobar measures (lobar cortical volume, lobar cortical surface area, and mean lobar cortical thickness), the same covariates were used. For global measures (such as whole brain volume (WBV), total grey matter volume, and total white matter volume), hemisphere was not a required covariate and ICV was also excluded at a covariate due to high collinearity with the measures. Imaging analyses were performed in n=1,024 in GS, and n=27,202 (n=26,639 for WBV) in UKB. These numbers represent those with complete imaging data as well as complete childhood trauma data. ### Mega-analysis of imaging samples For the mega-analysis, the total and subscale childhood trauma scores were scaled within each cohort, and then re-scaled across the merged mega-analysis sample. Structural imaging data including ICV was scaled within cohort. Analyses were conducted, as before, using the same covariates. ## Results ### Demographics characteristics #### Generation Scotland The GS study population (n=1,153) was 58.54% female, with a mean age of 59.52yrs (26-84yrs). At SCID interview 30.79% met criteria for lifetime depression (Table 1). View this table: [Table 1:](http://medrxiv.org/content/early/2022/06/07/2022.06.07.22276081/T1) Table 1: Demographic characteristics and psychiatric traits of the GS study population, incorporating participants with complete CTQ responses and SCID interviews. Prevalence of the 5 CTQ subscales are also shown. ** For the n=355 reporting lifetime depression* Using the CTQ score cut-offs (see supplementary materials, page 1), 26.02% of the study population reported experiencing ‘low-moderate’ levels of childhood trauma, or above. Using the subscale-specific cut-offs, 51% of the study population reported a score of ‘low-moderate’ or above on at least one CTQ subscale. The most frequent subtype of trauma reported was emotional neglect, and the least frequent were physical and sexual abuse (Table 1). Total CTQ scores were significantly higher in females (mean=34.75) than males (mean=32.94; t(df=1,138.8)=2.79, p=0.0054), and were associated with higher BMI (β=0.97, SE=0.17, p(FDR)=2.31×10−8), but were not significantly associated with age (β=-4.04×10−4, SE=0.0029, p(FDR) =0.89). Participants with lifetime MDD diagnoses reported higher total CTQ scores (mean=39.59) than healthy controls (mean=31.52; t(df=442.54)=-9.11, p<0.001). A significant positive association was found between total CTQ score and risk of lifetime depression (β=0.68, OR=1.98, p(FDR)=6.45×10−18). Depression risk was positively associated with total scores on all trauma subscales (Appendix 1). Higher total CTQ scores also predicted younger age-of-onset for depression (β=-3.2, 95% CI=-4.22 – - 2.17, p(FDR)=1.44×10−8), and higher odds for a recurrent course of depressive illness (β=0.36, OR=1.43, p(FDR)=0.0022). #### UK BioBank The UKB study population (n=153,650) was 56.31% female, with a mean age of 55.91 years (38-72yrs). Of n=94,379 participants for whom ICD-10 coded diagnostic information was available, 2.68% report lifetime depression. For the n=123,355 participants who completed a CIDI diagnostic questionnaire, however, 29.47% met criteria for lifetime depression (Table 2). View this table: [Table 2:](http://medrxiv.org/content/early/2022/06/07/2022.06.07.22276081/T2) Table 2: Demographic characteristics and psychiatric traits of the UKB study population, incorporating participants with complete childhood trauma responses. ** Data comes from initial assessment centre visit, not tied to a diagnostic questionnaire*. ** * Data comes from online follow-up to imaging appointment, not tied to a diagnostic questionnaire*. Due to the less detailed nature of the childhood trauma metric (CTM) items in UKB’s MHQ, the severity cut-offs inherent to the CTQ could not be replicated. The proportion of participants responding higher than ‘Rarely True’ (or lower than ‘Very Often True’ for reverse-coded items) to any CTM item was 59.28% (Table 2). Total CTM scores in UKB were significantly higher in females (mean=1.85) than in males (mean=1.64; t(df=152,796)=16.85, p<0.001). They showed a positive association with BMI (β=0.37, SE=0.012, p(FDR)=1.86×10−220), and a small but significant negative association with age (β=-0.0059, SE=3.3×10−4, p(FDR)=3.28×10−72). Participants with CIDI-coded lifetime MDD had higher CTM scores (mean=2.57) than healthy controls (mean=1.32; t(df=48,247)=-72.35, p<0.001); similarly, those with ICD-10 coded lifetime MDD had higher CTM scores (mean=3.34) than healthy controls (mean=1.67; t(df=2582)=-22.18, p<0.001). Significant positive associations were found between total CTM scores and risk of lifetime depression for CIDI (β=0.49, OR=1.63, p(FDR)<1×10−314), and ICD-10 data (β=0.4, OR=1.49, p(FDR)=8.52×10−186). ICD-10- and CIDI-diagnosed depression risk both positively associated with responses to each of the CTM subscale items, in addition (Appendix 1). Total CTM scores predicted younger age of depression onset (β=-2.1, SE=0.046, p(FDR)<1×10−314), and higher odds of a recurrent course of depressive illness (β=0.22, OR=1.24, p(FDR)=2.66×10−51). ### Neuroimaging findings We took the approach of analysing relationships between scaled CT scores and brain metrics in increasing levels of detail, starting with global brain and grey/white matter volumes, then lobar-level metrics, then across all regions. #### Global brain measures The relationships between CT scores and global brain metrics were analysed in GS (n=1,024), UKB (n=26,639), and the mega-analysis sample (n=28,226). In GS, significant negative associations were revealed between total CTQ score and global white matter volume (β=-0.0852, SE=0.0258, p(FDR)=0.0015), global grey matter volume (β=-0.0682, SE=0.0219, p(FDR)=0.0019), and whole brain volume (β=-0.0844, SE=0.0239, p(FDR)=0.0013; Appendix 2). In UKB, CTM scores showed significant negative associations with global white matter volume (β=-0.0311, SE=0.0049, p(FDR)=2.67×10−10), global grey matter volume (β=-0.0433, SE=0.0049, p(FDR)=5.27×10−18), and whole brain volume (β=-0.0399, SE=0.0048, p(FDR)=1.75×10−16; Appendix 2). In the mega-analysis population, the whole-brain analyses followed much the same pattern, with significant negative associations between childhood trauma scores and global white matter volume (β=-0.0330, SE=0.0048, p(FDR)=9.15×10−12), global grey matter volume (β=-0.0442, SE=0.0048, p(FDR)=1.74×10−19), and whole brain volume (β=-0.0417, SE=0.0047, p(FDR)=1.51×10−18; Appendix 2). #### Lobar measures The relationships between scaled CT scores and the whole-lobe cortical volume, whole-lobe cortical surface area, and mean cortical thickness of all five lobar regions were examined and analysed in GS (n=1,024), UKB (n=27,202), and the mega-analysis sample (n=28,226). In GS, there were significant negative associations between total CTQ scores and both cortical surface areas and cortical volumes of the frontal, parietal, temporal, occipital, and cingulate lobes (Appendix 2). There were no significant associations for lobar cortical thicknesses. In UKB, we observed significant negative associations between total CTM score and cortical volumes and cortical surface areas of all five lobes, but no associations with cortical thicknesses were observed (Appendix 2). In the mega-analysis sample, relationships between childhood trauma score and cortical volume and surface areas of all five lobes were significant (Appendix 2). Again, no associations between childhood trauma scores and whole-lobe cortical thickness were observed. #### Regional measures Regional analyses examined the relationships between CT scores and volumes of 34 cortical and 8 subcortical regions, as well as thickness and surface area for the same 34 cortical regions in GS (n=1,024). In UKB (n=27,202) and mega-analysis (n=28,226), 33 cortical regions were examined on the same metrics, plus 8 subcortical regional volumes. In GS, across all regional metrics analysed three regions were significantly associated after FDR correction (**Figure 2**, Appendix 2). These were the volumes of the hippocampus (β=-0.0582, SE=0.0247, p(FDR)=0.050), nucleus accumbens (β=-0.0654, SE=0.0231, p(FDR)=0.037), and ventral diencephalon (VDc; **Supplementary Figure 1**; β=-0.0526, SE=0.0218, p(FDR)=0.050). All three regions are subcortical areas, and all were negatively associated with total CTQ score, with higher CTQ score predicting lower volume of these regions. ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/06/07/2022.06.07.22276081/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2022/06/07/2022.06.07.22276081/F1) Figure 1: Analysis structure of the present study, demonstrating the initial discovery analysis in GS due to deeper phenotyping of childhood trauma, with replication in UKB and mega-analysis for maximum sample size. ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/06/07/2022.06.07.22276081/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2022/06/07/2022.06.07.22276081/F2) Figure 2: Lollipop plots showing –log10 of p(FDR) for each region and metric in the regional analysis for (A) GS, (B) UKB, and (C) the mega-analysis; and in the global and lobar analyses for (D) GS, (E) UKB, and (F) the meg-analysis. The –log10 of p(FDR)=0.05 is represented by the dotted grey line, any points exceeding this line achieved statistical significance after FDR correction. Regions of higher significance are labelled, the order in which other regions are presented can be found in Appendix 2. *In this figure only: Th = Cortical thickness, SA = Cortical surface area, CV = Cortical volume, ScV = Subcortical Volume*. In UKB, many significant associations were found (**Figure 2;** Appendix 2), due to the higher power afforded by the larger sample size in this data set. Of note was the replication of a significant negative association between total CTM score and the volume of the VDc (β=-0.0220, SE=0.0040, p(FDR)=2.84×10−7). Significant associations between CTM scores and cortical regions were demonstrated – again showing more consistent associations with SA and volume than cortical thickness. Regions with the strongest effect sizes were the volume of the inferior parietal lobule (β=-0.0223, SE=0.0042, p(FDR)=4.78×10−6), and the surface area of the pericalcarine cortex (β=-0.0248, SE=0.0054, p(FDR)=1.61×10−4). All findings from the UKB analyses were replicated in the mega-analysis, which also revealed some novel significant associations (**Figure 2, Figure 3**, Appendix 2). Two of the three findings in the GS regional analysis were also replicated in the mega-analysis; the significant negative associations in the hippocampus (β=-0.0119, SE=0.0044, p(FDR)=0.017), and the ventral diencephalon (β=-0.0232, SE=0.0039, p(FDR)=2.91×10−8). The volume of the thalamus was additionally significant at mega-analysis (β=-0.0134, SE=0.0039, p(FDR)=0.0026). ![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/06/07/2022.06.07.22276081/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2022/06/07/2022.06.07.22276081/F3) Figure 3: Map of showing Beta values for regions significantly associated with childhood trauma in the mega-analysis, for cortical thickness, cortical surface area, and cortical volume. The red-toned colours represent a negative association and the purple-toned represent a positive association, with the strength of that associated denoted by the shade of the colour. Images created using the Cox & Liewald *Heatmapper* programme (66) #### Childhood trauma subscales Across all three global brain metrics in GS, significant associations were demonstrated with Emotional Abuse, Physical Abuse, Sexual Abuse, and Abuse composite score. This was replicated in the UKB and mega-analysis. In GS, however, no significant association was indicated for Emotional Neglect, Physical Neglect, or Neglect Composite score. This was not replicated in the UKB and mega-analysis; EN, PN, and Neglect composite were significantly associated with all three global brain metrics in these analyses (Appendix 3). Similarly, at the lobar level, imaging samples demonstrated associations with some subscale scores with particularly consistent results between PA, SA, and abuse composite score and lobar surface areas; and between PA and abuse composite score and lobar volumes. Very few significant associations were found with lobar thicknesses, excepting significant associations with both abuse and neglect composite scores in the UKB analysis only (Appendix 3). At the regional level in GS, the only significant associations observed were between the physical abuse score and the surface area of superior parietal cortex (β=-0.0771, SE=0.0235, p(FDR)=0.036), the cortical volume of the entorhinal cortex (β=-0.0733, SE=0.0224, p(FDR)=0.037), and 6 of the 8 subcortical volumes: the hippocampus, VDc, thalamus, amygdala, nucleus accumbens, and pallidum (Appendix 3). No other subscale or composite score showed any significant associations with any regional metric in GS, except for the neglect composite score which associated significantly with the volume of the nucleus accumbens (β=-0.0805, SE=0.0229, p(FDR)=0.0037). In UKB and in the mega-analysis, many associations were found between the different trauma subtypes and regional imaging metrics, including some associations between the emotional abuse items and cortical thickness measures (Appendix 3). The abuse composite score demonstrated significant associations with the volumes of the hippocampus, thalamus, and ventral diencephalon in both the UKB analysis and mega-analysis – the same three subcortical regions with which total childhood trauma scores associated in the mega-analysis (Appendix 3). The only region for which a significant association was found across all three analyses was the VDc (GS: β=-0.0676, SE=0.0217, p(FDR)=0.005; UKB: β=-0.0212, SE=0.004, p(FDR)=1.01×10−6; Mega-analysis: β=-0.023, SE=0.0039, p(FDR)=4.19×10−8). ## Discussion This study reports psychiatric and structural brain correlates of childhood trauma in a mega-analysis of two independent UK cohorts. The mega-analysis sample size of n=28,226 improves significantly upon the sample sizes of similar multicentre studies (41, 42). The data presented here provide further validation for proposed relationships between childhood trauma and depressive pathology in adulthood. These links to poor mental health outcomes have been demonstrated previously (67), as have links to poor physical health outcomes (4). The study also highlights the high prevalence of traumatic childhood experiences in our population, with 51% and 59% of participants reporting having experienced some kind of adverse childhood experience in GS and UKB, respectively. In the GS CTQ data, where more detailed analysis was possible, we observed what might be considered a ‘clinically significant’ level of trauma (a total CTQ score over the ‘low-moderate’ threshold) in 26% of the study population. This fits with previous estimates of childhood trauma prevalence (1, 2, 4). The association between higher CT score and younger age of onset for MDD is strong in both cohorts (GS: β=-3.2, 95% CI=-4.22 – -2.17, p(FDR)=1.44×10−8; UKB: β=-2.1, SE=0.046, p(FDR)<1×10−314). This replicated indicates that CT is having early effects on the development of psychopathology, and given the links between age of onset and poorer long-term outcomes in MDD could also suggest a link with lifetime severity (68). The severity of traumatic childhood experiences associated with global brain volume, global grey and white matter volumes, and the volumes and surface areas of all five major brain lobes at mega-analysis with effect sizes ranging from -0.0442 – -0.0204 (Appendix 2). These results point to a whole-brain effect of adverse childhood events on brain development, perhaps with exacerbated results in particular regions. Interestingly, in the global, lobar, and regional analyses very few associations are seen between CT scores and cortical thickness metrics. The evidence of this work suggests that CT severity may be more strongly associated with cortical surface areas and volumes, especially in front and temporo-parietal regions (**Figure 3**). Specific reductions in cortical surface area rather than thickness have also been reported in the context of adolescent depression (69), this could reflect a common feature associated with early onset of psychiatric morbidity. Cortical surface area has been shown to be established earlier in development within a shorter window than cortical thickness (70); adverse experiences within this earlier window could be interrupting the developmental process resulting in the surface area findings presented here. Cortical thickness is established over a much longer time period, well into adulthood (70), and may therefore be less susceptible to interruption by early life events. The effect sizes found in the global and lobar analyses presented here are larger than those found in the regional analyses, but there were nevertheless multiple significant associations demonstrated. Total CT scores was significantly associated with multiple cortical and subcortical regional volumes at mega-analysis, with effect sizes ranging from -0.0262 – 0.0208. The strongest effect sizes were seen in the negative relationship between total childhood trauma scores and the volume of the ventral diencephalon (β=-0.0232, SE=0.0039, p(FDR)=2.91×10−8), the volume of the inferior parietal cortex (β=-0.0225, SE=0.0042, p(FDR)=2.43×10−6), and the surface area of the pericalcarine cortex (β=-0.0262, SE=0.0053, p(FDR)=2.75×10−5); as well as the positive association with the thickness of the lingual cortex (β=0.0208, SE=0.0053, p(FDR)=0.0026). The subcortical regions the hippocampus (β=-0.0119, SE=0.0044, p(FDR)=0.017) and the thalamus (β=-0.0134, SE=0.0039, p(FDR)=0.0026) were also significantly associated with total CT scores at mega-analysis, and in the GS analysis the nucleus accumbens was significant with a relatively strong effect size (β=-0.0654, SE=0.023, p(FDR)=0.037). We have previously reported an association between elevated hair glucocorticoids and the volume of the nucleus accumbens (as well as with MDD and CT) in the GS imaging sample, suggesting a possible link to the HPA axis and stress reactivity (71). The ventral diencephalon was the only regional metric in which a consistent effect was seen across both GS and UK Biobank, as well as in the mega-analysis. The VDc is a *FreeSurfer* derived region comprising several structures, including the hypothalamus – a region crucial to the regulation of the stress response – as well as subthalamic nuclei, the substantia nigra, the medial and lateral geniculate nuclei and associated white matter structures (**Figure 2**; as described (72, 73)). A significant relationship was also observed between CT and the volume of the thalamus (β=-0.0134, SE=0.0039, p(FDR)=0.0026). This could implicate involvement of thalamic/hypothalamic function after CT, supporting previous work (31-37). Future work could employ finer-grained approaches to probe associations with this region more deeply. It is well established that different types of trauma may affect the brain in different ways (10). The subscale mega-analysis presented here generally found an influence of physical abuse on a greater number of regional surface areas and cortical volumes than other subscales, although physical neglect also associated with a high number of regional surface areas (Appendix 3). Emotional abuse was the only subscale which had any associations with regional cortical thicknesses (Appendix 3). Subcortical regional volumes were affected by all subscales of childhood trauma, with the abuse composite measure showing the same pattern of effects as the total CT score – negatively associating with the hippocampus (β=-0.0124, SE=0.044, p(FDR)=0.013), thalamus (β=-0.0166, SE=0.0040, p(FDR)=1.03××10−4), and VDc (β=-0.0249, SE=0.0040, p(FDR)=2.42×10−9) – indicating that abuse may be more strongly associated with structural brain changes after CT than neglect. It is important, however, to treat the UKB and mega-analysis subscales with caution due to the single-factor nature of the subscale items in UKB. In GS alone, due to the lower power of the cohort in comparison, few significant associations were found in the subscale analyses. Notably, physical abuse associated significantly with the surface area of the superior parietal cortex (β=-0.0771, SE=0.024, p(FDR)=0.036), and the volume of the entorhinal cortex (β=-0.0733, SE=0.022, p(FDR)=0.037), as well as 6 of 8 subcortical regions (Appendix 3). No other subscale associated with any cortical measures in GS, although physical and emotional neglect both associated with the volume of the nucleus accumbens (β=-0.0681, SE=0.023, p(FDR)=0.025; β=-0.0743, SE=0.023, p(FDR)=0.0099 respectively), for which the neglect composite score was accordingly significant (β=-0.0805, SE=0.023, p(FDR)=0.0037). In the global and lobar level analyses in GS, physical and sexual abuse associated with most lobar surface areas, and physical abuse associated with most lobar volumes (Appendix 3). The physical, emotional, and sexual abuse items and the abuse composite score each associated with the global brain volume, and global grey and white matter measures (effect sizes ranging from -0.114 --0480; Appendix 3); meanwhile, neither neglect item not the composite score associated with any global measure in GS alone. This suggests that the effect of childhood trauma on global brain development may be more influenced by experiences of abuse than neglect, with physical abuse seeming to display the strongest influence. ### Limitations Despite the very strong sample size derived from mega-analysis of two independent UK cohorts, enabling well-powered brain-wide analysis and conservative correction for multiple testing, the multi-sample nature of this study introduces unavoidable limitations. Chiefly the use of different childhood trauma measures, upon which this work relies heavily. The GS subsample employed the full CTQ-28, while the UKB study used a five-item CT metric comprising one item taken from each subscale of the CTQ-28. While this does mean that the two measures share features, they are not precisely equivalent. The cruder nature of the UKB measure is evident in the high report of emotional neglect in UKB (47.59% responding ‘rarely true’ or above, compared to 31.92% reporting a ‘low-moderate’ score or above in GS). The single item used to denote emotional neglect in UKB was ‘When I was growing up… I felt loved’, score ‘never true’ to ‘very often true’. This item carries a high emotional valence, which may have led respondents to over-report compared to the scores derived from a multi-factorial scale in the full CTQ-28. The subscale analyses presented here should be treated particularly cautiously as a result, in both UKB and in the mega-analysis sample due to the heavy representation of UKB participants there. Other metrics used in the two cohorts – such as psychiatric diagnostic tools – differed, leading again to a degree of mismatch in the data derived from each cohort. The imaging protocols, too, were conducted at different sites and scanners and according to different protocols, although the ‘scan site’ variable employed in all imaging analyses attempts to account for this. Despite these limitations, however, the two samples were relatively homogenous compared to previous multi-cohort studies, as can be seen by comparing the study population demographics presented for each cohort (**Table 1, Table 2**). Further limitations are evident in the nature of the CTQ in itself – retrospective report of traumatic childhood experiences has been shown to identify a different population than prospective report (74). While it is difficult to conceive of a study on this scale that does not rely on retrospective report, it is important to consider these results in light of this fact. Neither cohort provides any information on the timing of reported trauma, either. This has been shown to affect outcomes, with different types of trauma affecting brain development in different ways according to the timing of these events (75). An additional concern relevant to any population imaging cohort is evidence that ‘healthy bias’ may exist in the group of participants who were able and willing to attend MRI scanning appointments as compared to the un-scanned members of the cohort. Healthy bias has been demonstrated in UKB (76) but may well be present in GS, thus participants with more extreme symptoms following CT may not have been included in the study populations reported here. ### Conclusions This large mega-analysis study reports associations between childhood trauma and altered structural metrics across the adult brain, with analyses performed brain-wide in a non-hypothesis driven manner. These effects can be seen at the global, lobar, and regional levels. A negative association between childhood trauma metric scores and the ventral diencephalon – an area comprising the hypothalamus and other thalamus-associated regions – replicated across all three analyses. The effects on global brain measures are particularly noteworthy, and support evidence of an involvement of traumatic childhood experiences in global brain development. ## Supporting information Suplementary Information [[supplements/276081_file02.docx]](pending:yes) Appendix 3 [[supplements/276081_file03.xlsx]](pending:yes) Appendix 2 [[supplements/276081_file04.xlsx]](pending:yes) Appendix 1 [[supplements/276081_file05.xlsx]](pending:yes) ## Data Availability All data produced in the present study are available upon reasonable request to the authors ## List of Abbreviations CA : Childhood adversity ACEs : Adverse childhood experiences CT : Childhood trauma CTQ : Childhood trauma questionnaire CTM : childhood trauma metric (shorthand for UKB childhood trauma items) SCID : Structured clinical interview for DSM-IV VDc : Ventral diencephalon EN : Emotional Neglect EA : Emotional Abuse PN : Physical Neglect PA : Physical Abuse SA : Sexual Abuse * Received June 7, 2022. * Revision received June 7, 2022. * Accepted June 7, 2022. * © 2022, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/) ## References 1. 1.Bellis MA, Ashton K, Hughes K, Ford K, Bishop J, Paranjothy S. Adverse childhood experiences and their impact on health-harming behaviours in the Welsh adult population: Public Health Wales NHS Trust; 2016. 2. 2.Bellis MA, Lowey H, Leckenby N, Hughes K, Harrison D. Adverse childhood experiences: retrospective study to determine their impact on adult health behaviours and health outcomes in a UK population. Journal of public health. 2013;36(1):81–91. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23587573&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F07%2F2022.06.07.22276081.atom) 3. 3.Felitti VJ, Anda RF, Nordenberg D, Williamson DF, Spitz AM, Edwards V, et al. Relationship of Childhood Abuse and Household Dysfunction to Many of the Leading Causes of Death in Adults. American Journal of Preventive Medicine. 1998;14(4). 4. 4.Felitti VJ, Anda RF, Nordenberg D, Williamson DF, Spitz AM, Edwards V, et al. Relationship of childhood abuse and household dysfunction to many of the leading causes of death in adults: The Adverse Childhood Experiences (ACE) Study. American journal of preventive medicine. 2019;56(6):774–86. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.amepre.2019.04.001&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F07%2F2022.06.07.22276081.atom) 5. 5.Bernstein DP, Stein JA, Newcomb MD, Walker E, Pogge D, Ahluvalia T, et al. Development and validation of a brief screening version of the Childhood Trauma Questionnaire. Child abuse & neglect. 2003;27(2):169–90. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0145-2134(02)00541-0&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12615092&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F07%2F2022.06.07.22276081.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000181130900005&link_type=ISI) 6. 6.Krupnick JL, Green BL, Stockton P, Goodman L, Corcoran C, Petty R. Mental health effects of adolescent trauma exposure in a female college sample: Exploring differential outcomes based on experiences of unique trauma types and dimensions. Psychiatry: Interpersonal and Biological Processes. 2004;67(3):264–79. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1521/psyc.67.3.264.48986&link_type=DOI) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000224686900010&link_type=ISI) 7. 7.Schalinski I, Teicher MH, Nischk D, Hinderer E, Müller O, Rockstroh B. Type and timing of adverse childhood experiences differentially affect severity of PTSD, dissociative and depressive symptoms in adult inpatients. BMC Psychiatry. 2016;16(1):295. 8. 8.Heim C, Newport DJ, Mletzko T, Miller AH, Nemeroff CB. The link between childhood trauma and depression: insights from HPA axis studies in humans. Psychoneuroendocrinology. 2008;33(6):693–710. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.psyneuen.2008.03.008&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18602762&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F07%2F2022.06.07.22276081.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000258356200001&link_type=ISI) 9. 9.Cook A, Spinazzola J, Ford J, Lanktree C, Blaustein M, Cloitre M, et al. Complex trauma in children and adolescents. Psychiatric annals. 2017;35(5):390–8. 10. 10.Cassiers LLM, Sabbe BGC, Schmaal L, Veltman DJ, Penninx BWJH, Van Den Eede F. Structural and Functional Brain Abnormalities Associated With Exposure to Different Childhood Trauma Subtypes: A Systematic Review of Neuroimaging Findings. Front Psychiatry. 2018;9:329-. 11. 11.Teicher MH, Andersen SL, Polcari A, Anderson CM, Navalta CP, Kim DM. The neurobiological consequences of early stress and childhood maltreatment. Neuroscience & Biobehavioral Reviews. 2003;27(1):33–44. 12. 12.Videbech P, Ravnkilde B. Hippocampal volume and depression: a meta-analysis of MRI studies. American Journal of Psychiatry. 2004;161(11):1957–66. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1176/appi.ajp.161.11.1957&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15514393&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F07%2F2022.06.07.22276081.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000224768400003&link_type=ISI) 13. 13.Sapolsky RM, Uno H, Rebert CS, Finch CE. Hippocampal damage associated with prolonged glucocorticoid exposure in primates. Journal of Neuroscience. 1990;10(9):2897–902. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Njoiam5ldXJvIjtzOjU6InJlc2lkIjtzOjk6IjEwLzkvMjg5NyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzA2LzA3LzIwMjIuMDYuMDcuMjIyNzYwODEuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 14. 14.Malhi GS, Das P, Outhred T, Dobson-Stone C, Irwin L, Gessler D, et al. Effect of stress gene-by-environment interactions on hippocampal volumes and cortisol secretion in adolescent girls. Australian & New Zealand Journal of Psychiatry. 2019;53(4):316–25. 15. 15.Dannlowski U, Stuhrmann A, Beutelmann V, Zwanzger P, Lenzen T, Grotegerd D, et al. Limbic Scars: Long-Term Consequences of Childhood Maltreatment Revealed by Functional and Structural Magnetic Resonance Imaging. Biological Psychiatry. 2012;71(4):286–93. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.biopsych.2011.10.021&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22112927&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F07%2F2022.06.07.22276081.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000299314800006&link_type=ISI) 16. 16.Calem M, Bromis K, McGuire P, Morgan C, Kempton MJ. Meta-analysis of associations between childhood adversity and hippocampus and amygdala volume in non-clinical and general population samples. NeuroImage: Clinical. 2017;14:471–9. 17. 17.Colle R, Segawa T, Chupin M, Dong MNTKT, Hardy P, Falissard B, et al. Early life adversity is associated with a smaller hippocampus in male but not female depressed in-patients: a case–control study. BMC psychiatry. 2017;17(1):71. 18. 18.Chaney A, Carballedo A, Amico F, Fagan A, Skokauskas N, Meaney J, et al. Effect of childhood maltreatment on brain structure in adult patients with major depressive disorder and healthy participants. J Psychiatry Neurosci. 2014;39(1):50–9. 19. 19.Lawson GM, Camins JS, Wisse L, Wu J, Duda JT, Cook PA, et al. Childhood socioeconomic status and childhood maltreatment: Distinct associations with brain structure. PLOS ONE. 2017;12(4):e0175690. 20. 20.Frodl T, Reinhold E, Koutsouleris N, Reiser M, Meisenzahl EM. Interaction of childhood stress with hippocampus and prefrontal cortex volume reduction in major depression. Journal of Psychiatric Research. 2010;44(13):799–807. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jpsychires.2010.01.006&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20122698&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F07%2F2022.06.07.22276081.atom) 21. 21.Mielke EL, Neukel C, Bertsch K, Reck C, Möhler E, Herpertz SC. Alterations of brain volumes in women with early life maltreatment and their associations with oxytocin. Hormones and Behavior. 2018;97:128–36. 22. 22.Saleh A, Potter GG, McQuoid DR, Boyd B, Turner R, MacFall JR, et al. Effects of early life stress on depression, cognitive performance and brain morphology. Psychological medicine. 2017;47(1):171–81. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S0033291716002403&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27682320&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F07%2F2022.06.07.22276081.atom) 23. 23.Cohen RA, Grieve S, Hoth KF, Paul RH, Sweet L, Tate D, et al. Early Life Stress and Morphometry of the Adult Anterior Cingulate Cortex and Caudate Nuclei. Biological Psychiatry. 2006;59(10):975–82. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.biopsych.2005.12.016&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16616722&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F07%2F2022.06.07.22276081.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000237676100012&link_type=ISI) 24. 24.Rokita KI, Holleran L, Dauvermann MR, Mothersill D, Holland J, Costello L, et al. Childhood trauma, brain structure and emotion recognition in patients with schizophrenia and healthy participants. Social Cognitive and Affective Neuroscience. 2020;15(12):1325–39. 25. 25.Paquola C, Bennett MR, Lagopoulos J. Understanding heterogeneity in grey matter research of adults with childhood maltreatment—A meta-analysis and review. Neuroscience & Biobehavioral Reviews. 2016;69:299–312. 26. 26.LeDoux J. The emotional brain, fear, and the amygdala. Cellular and molecular neurobiology. 2003;23(4-5):727-38. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1023/A:1025048802629&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=14514027&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F07%2F2022.06.07.22276081.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000184589000023&link_type=ISI) 27. 27.Janiri D, Sani G, Rossi PD, Piras F, Iorio M, Banaj N, et al. Amygdala and hippocampus volumes are differently affected by childhood trauma in patients with bipolar disorders and healthy controls. Bipolar disorders. 2017;19(5):353–62. 28. 28.Bruno F, Splendiani A, Tommasino E, Conson M, Quarantelli M, Saporito G, et al. Multimodal MRI Assessment of Thalamic Structural Changes in Earthquake Survivors. Diagnostics. 2021;11(1):70. 29. 29.Shucard JL, Cox J, Shucard DW, Fetter H, Chung C, Ramasamy D, et al. Symptoms of posttraumatic stress disorder and exposure to traumatic stressors are related to brain structural volumes and behavioral measures of affective stimulus processing in police officers. Psychiatry Research: Neuroimaging. 2012;204(1):25–31. 30. 30.Nilsen AS, Blix I, Leknes S, Ekeberg Ø, Skogstad L, Endestad T, et al. Brain Activity in Response to Trauma-specific, Negative, and Neutral Stimuli. A fMRI Study of Recent Road Traffic Accident Survivors. Frontiers in Psychology. 2016;7(1173). 31. 31.Ancelin M-L, Carrière I, Artero S, Maller JJ, Meslin C, Dupuy A-M, et al. Structural brain alterations in older adults exposed to early-life adversity. Psychoneuroendocrinology. 2021;129:105272. 32. 32.Kumari V, Gudjonsson GH, Raghuvanshi S, Barkataki I, Taylor P, Sumich A, et al. Reduced thalamic volume in men with antisocial personality disorder or schizophrenia and a history of serious violence and childhood abuse. European Psychiatry. 2013;28(4):225–34. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.eurpsy.2012.03.002&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22944337&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F07%2F2022.06.07.22276081.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000319548700005&link_type=ISI) 33. 33.Philip NS, Tyrka AR, Albright SE, Sweet LH, Almeida J, Price LH, et al. Early life stress predicts thalamic hyperconnectivity: A transdiagnostic study of global connectivity. Journal of Psychiatric Research. 2016;79:93–100. 34. 34.Carvalho Fernando S, Beblo T, Schlosser N, Terfehr K, Otte C, Löwe B, et al. Associations of childhood trauma with hypothalamic-pituitary-adrenal function in borderline personality disorder and major depression. Psychoneuroendocrinology. 2012;37(10):1659–68. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.psyneuen.2012.02.012&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22444624&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F07%2F2022.06.07.22276081.atom) 35. 35.Videlock EJ, Adeyemo M, Licudine A, Hirano M, Ohning G, Mayer M, et al. Childhood Trauma Is Associated With Hypothalamic-Pituitary-Adrenal Axis Responsiveness in Irritable Bowel Syndrome. Gastroenterology. 2009;137(6):1954–62. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1053/j.gastro.2009.08.058&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19737564&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F07%2F2022.06.07.22276081.atom) 36. 36.Schäfer I, Teske L, Schulze-Thüsing J, Homann K, Reimer J, Haasen C, et al. Impact of Childhood Trauma on Hypothalamus-Pituitary-Adrenal Axis Activity in Alcohol-Dependent Patients. European Addiction Research. 2010;16(2):108–14. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20224278&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F07%2F2022.06.07.22276081.atom) 37. 37.Kempke S, Luyten P, De Coninck S, Van Houdenhove B, Mayes LC, Claes S. Effects of early childhood trauma on hypothalamic–pituitary–adrenal (HPA) axis function in patients with Chronic Fatigue Syndrome. Psychoneuroendocrinology. 2015;52:14–21. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.psyneuen.2014.10.027&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25459889&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F07%2F2022.06.07.22276081.atom) 38. 38.Opel N, Redlich R, Dohm K, Zaremba D, Goltermann J, Repple J, et al. Mediation of the influence of childhood maltreatment on depression relapse by cortical structure: a 2-year longitudinal observational study. The Lancet Psychiatry. 2019;6(4):318–26. 39. 39.Lim L, Radua J, Rubia K. Gray matter abnormalities in childhood maltreatment: a voxel-wise meta-analysis. American Journal of Psychiatry. 2014;171(8):854–63. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1176/appi.ajp.2014.13101427&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24781447&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F07%2F2022.06.07.22276081.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000339988900011&link_type=ISI) 40. 40.Thompson PM, Stein JL, Medland SE, Hibar DP, Vasquez AA, Renteria ME, et al. The ENIGMA Consortium. large-scale collaborative analyses of neuroimaging and genetic data. 2014;N/A(N/A):N/A. 41. 41.Frodl T, Janowitz D, Schmaal L, Tozzi L, Dobrowolny H, Stein DJ, et al. Childhood adversity impacts on brain subcortical structures relevant to depression. Journal of Psychiatric Research. 2017;86:58–65. 42. 42.Tozzi L, Garczarek L, Janowitz D, Stein DJ, Wittfeld K, Dobrowolny H, et al. Interactive impact of childhood maltreatment, depression, and age on cortical brain structure: mega-analytic findings from a large multi-site cohort. Psychological Medicine. 2020;50(6):1020–31. 43. 43.Habota T, Sandu A, Waiter G, McNeil C, Steele J, Macfarlane J, et al. Cohort profile for the STratifying Resilience and Depression Longitudinally (STRADL) study: A depression-focused investigation of Generation Scotland, using detailed clinical, cognitive, and neuroimaging assessments [version 1; peer review: 1 approved, 1 not approved]. Wellcome Open Research. 2019;4(185). 44. 44.Navrady LB, Wolters MK, MacIntyre DJ, Clarke T-K, Campbell AI, Murray AD, et al. Cohort Profile: Stratifying Resilience and Depression Longitudinally (STRADL): a questionnaire follow-up of Generation Scotland: Scottish Family Health Study (GS:SFHS). International Journal of Epidemiology. 2017;47(1):13–4g. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ije/dyx115&link_type=DOI) 45. 45.Stolicyn A, Harris MA, Shen X, Barbu MC, Adams MJ, Hawkins EL, et al. Automated classification of depression from structural brain measures across two independent community-based cohorts. Human Brain Mapping. 2020;41(14):3922–37. 46. 46.Green C, Shen X, Stevenson AJ, Conole ELS, Harris MA, Barbu MC, et al. Structural brain correlates of serum and epigenetic markers of inflammation in major depressive disorder. Brain Behav Immun. 2021;92:39–48. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/J.BBI.2020.11.024&link_type=DOI) 47. 47.Pereira J, Vickers K, Atkinson L, Gonzalez A, Wekerle C, Levitan R. Parenting stress mediates between maternal maltreatment history and maternal sensitivity in a community sample. Child Abuse & Neglect. 2012;36(5):433–7. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.chiabu.2012.01.006&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22633056&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F07%2F2022.06.07.22276081.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000305367700006&link_type=ISI) 48. 48.Gerdner A, Allgulander C. Psychometric properties of the Swedish version of the Childhood Trauma Questionnaire—Short Form (CTQ-SF). Nordic Journal of Psychiatry. 2009;63(2):160–70. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/08039480802514366&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19021077&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F07%2F2022.06.07.22276081.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000264482200010&link_type=ISI) 49. 49.Schmidt MR, Narayan AJ, Atzl VM, Rivera LM, Lieberman AF. Childhood Maltreatment on the Adverse Childhood Experiences (ACEs) Scale versus the Childhood Trauma Questionnaire (CTQ) in a Perinatal Sample. Journal of Aggression, Maltreatment & Trauma. 2020;29(1):38–56. 50. 50.Paivio SC, Cramer KM. Factor structure and reliability of the Childhood Trauma Questionnaire in a Canadian undergraduate student sample. Child Abuse & Neglect. 2004;28(8):889–904. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.chiabu.2004.01.011&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15350772&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F07%2F2022.06.07.22276081.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000223876200007&link_type=ISI) 51. 51.First MB. Structured clinical interview for the DSM (SCID). The encyclopedia of clinical psychology. 2014:1–6. 52. 52.Rush AJ, Trivedi MH, Ibrahim HM, Carmody TJ, Arnow B, Klein DN, et al. The 16-Item Quick Inventory of Depressive Symptomatology (QIDS), clinician rating (QIDS-C), and self-report (QIDS-SR): a psychometric evaluation in patients with chronic major depression. Biological psychiatry. 2003;54(5):573–83. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0006-3223(02)01866-8&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12946886&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F07%2F2022.06.07.22276081.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000184906100009&link_type=ISI) 53. 53.Rupprechter S, Romaniuk L, Series P, Hirose Y, Hawkins E, Sandu A-L, et al. Blunted medial prefrontal cortico-limbic reward-related effective connectivity and depression. Brain. 2020;143(6):1946–56. 54. 54.Romaniuk L, Sandu A-L, Waiter GD, McNeil CJ, Xueyi S, Harris MA, et al. The Neurobiology of Personal Control During Reward Learning and Its Relationship to Mood. Biological Psychiatry: Cognitive Neuroscience and Neuroimaging. 2019;4(2):190–9. 55. 55.Dale AM, Fischl B, Sereno MI. Cortical surface-based analysis: I. Segmentation and surface reconstruction. Neuroimage. 1999;9(2):179–94. 56. 56.Fischl B, Salat DH, Van Der Kouwe AJ, Makris N, Ségonne F, Quinn BT, et al. Sequence-independent segmentation of magnetic resonance images. Neuroimage. 2004;23:S69–S84. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2004.07.016&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15501102&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F07%2F2022.06.07.22276081.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000225374100007&link_type=ISI) 57. 57.Fischl B, Sereno MI, Dale AM. Cortical surface-based analysis: II: inflation, flattening, and a surface-based coordinate system. Neuroimage. 1999;9(2):195–207. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1006/nimg.1998.0396&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9931269&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F07%2F2022.06.07.22276081.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000078608900002&link_type=ISI) 58. 58.Desikan RS, Ségonne F, Fischl B, Quinn BT, Dickerson BC, Blacker D, et al. An automated labeling system for subdividing the human cerebral cortex on MRI scans into gyral based regions of interest. Neuroimage. 2006;31(3):968–80. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2006.01.021&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16530430&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F07%2F2022.06.07.22276081.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000238358500003&link_type=ISI) 59. 59.Sudlow C, Gallacher J, Allen N, Beral V, Burton P, Danesh J, et al. UK Biobank: An Open Access Resource for Identifying the Causes of a Wide Range of Complex Diseases of Middle and Old Age. PLOS Medicine. 2015;12(3):e1001779. 60. 60.Miller KL, Alfaro-Almagro F, Bangerter NK, Thomas DL, Yacoub E, Xu J, et al. Multimodal population brain imaging in the UK Biobank prospective epidemiological study. Nature Neuroscience. 2016;19:1523. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nn.4393&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27643430&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F07%2F2022.06.07.22276081.atom) 61. 61.Alfaro-Almagro F, McCarthy P, Afyouni S, Andersson JLR, Bastiani M, Miller KL, et al. Confound modelling in UK Biobank brain imaging. NeuroImage. 2021;224:117002. 62. 62.Smith DJ, Nicholl BI, Cullen B, Martin D, Ul-Haq Z, Evans J, et al. Prevalence and Characteristics of Probable Major Depression and Bipolar Disorder within UK Biobank: Cross-Sectional Study of 172,751 Participants. PLOS ONE. 2013;8(11):e75362. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0075362&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24282498&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F07%2F2022.06.07.22276081.atom) 63. 63.Kessler RC, Andrews G, Mroczek D, Ustun B, Wittchen HU. The World Health Organization composite international diagnostic interview short-form (CIDI-SF). International journal of methods in psychiatric research. 1998;7(4):171–85. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/mpr.47&link_type=DOI) 64. 64.Spitzer RL, Kroenke K, Williams JBW, Group atPHQPCS. Validation and Utility of a Self-report Version of PRIME-MDThe PHQ Primary Care Study. JAMA. 1999;282(18):1737–44. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jama.282.18.1737&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10568646&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F07%2F2022.06.07.22276081.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000083471400028&link_type=ISI) 65. 65.Alfaro-Almagro F, Jenkinson M, Bangerter NK, Andersson JLR, Griffanti L, Douaud G, et al. Image processing and Quality Control for the first 10,000 brain imaging datasets from UK Biobank. NeuroImage. 2018;166:400–24. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F07%2F2022.06.07.22276081.atom) 66. 66.Liewald DCS. The Liewald-Cox Heatmapper Tool. 67. 67. Martin H. Teicher, M.D., Ph.D.., and, Jacqueline A. Samson, Ph.D.. Childhood Maltreatment and Psychopathology: A Case for Ecophenotypic Variants as Clinically and Neurobiologically Distinct Subtypes. American Journal of Psychiatry. 2013;170(10):1114–33. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1176/appi.ajp.2013.12070957&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23982148&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F07%2F2022.06.07.22276081.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000325383700013&link_type=ISI) 68. 68.Zisook S, Lesser I, Stewart JW, Wisniewski SR, Balasubramani GK, Fava M, et al. Effect of age at onset on the course of major depressive disorder. The American Journal of Psychiatry. 2007;164(10):1539–46. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1176/appi.ajp.2007.06101757&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17898345&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F07%2F2022.06.07.22276081.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000250049600018&link_type=ISI) 69. 69.Shen X, MacSweeney N, Chan SW, Barbu MC, Adams MJ, Romaniuk L, et al. Brain Structural Associations with Depression in a Large Early Adolescent Sample (The ABCD Cohort). Available at SSRN 3719896. 70. 70.Amlien IK, Fjell AM, Tamnes CK, Grydeland H, Krogsrud SK, Chaplin TA, et al. Organizing Principles of Human Cortical Development--Thickness and Area from 4 to 30 Years: Insights from Comparative Primate Neuroanatomy. Cereb Cortex. 2016;26(1):257–67. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cercor/bhu214&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25246511&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F07%2F2022.06.07.22276081.atom) 71. 71.Green C, Stolicyn A, Harris MA, Shen X, Romaniuk L, Barbu MC, et al. Hair glucocorticoids are associated with childhood adversity, depressive symptoms and reduced global and lobar grey matter in Generation Scotland. Translational Psychiatry. 2021;11(1):523. 72. 72.Makris N, Oscar-Berman M, Jaffin SK, Hodge SM, Kennedy DN, Caviness VS, et al. Decreased Volume of the Brain Reward System in Alcoholism. Biological Psychiatry. 2008;64(3):192–202. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.biopsych.2008.01.018&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18374900&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F07%2F2022.06.07.22276081.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000257640700004&link_type=ISI) 73. 73.Lebedeva AK, Westman E, Borza T, Beyer MK, Engedal K, Aarsland D, et al. MRI-Based Classification Models in Prediction of Mild Cognitive Impairment and Dementia in Late-Life Depression. Front Aging Neurosci. 2017;9:13-. 74. 74.Baldwin JR, Reuben A, Newbury JB, Danese A. Agreement Between Prospective and Retrospective Measures of Childhood Maltreatment: A Systematic Review and Meta-analysis. JAMA Psychiatry. 2019;76(6):584–93. 75. 75.Andersen SL, Tomada A, Vincow ES, Valente E, Polcari A, Teicher MH. Preliminary evidence for sensitive periods in the effect of childhood sexual abuse on regional brain development. J Neuropsychiatry Clin Neurosci. 2008;20(3):292–301. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1176/appi.neuropsych.20.3.292&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18806232&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F07%2F2022.06.07.22276081.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000259599500004&link_type=ISI) 76. 76.Lyall DM, Quinn T, Lyall LM, Ward J, Anderson JJ, Smith DJ, et al. Quantifying bias in psychological and physical health in the UK Biobank imaging sub-sample. Brain Communications. 2022.