Introduction
According to the Centers for Disease Control and Prevention, childhood maltreatment (CM) is defined as ‘any act or series of acts of commission or omission by a parent or other caregiver that results in harm, potential for harm, or threat of harm to a child’ (Leeb et al., Reference Leeb, Paulozzi, Melanson, Simon and Arias2008). CM may be physical, sexual or emotional and may result in inadequate environmental input (e.g. deprivation or neglect) or excessive harmful input (Sheridan and McLaughlin, Reference Sheridan and McLaughlin2014). About one-quarter of all adults have encountered CM in their life (Fu et al., Reference Fu, Feng, Qin, Wang, Wu, Cai, Lan and Yang2018) and this statistic may even be higher as a history of childhood adversity is likely under-reported. In fact, a recent meta-analysis has found that over half of children globally had experienced violence in just the past year alone (Hillis et al., Reference Hillis, Mercy, Amobi and Kress2016). Given the prevalence of CM, this is especially alarming as CM is strongly associated with a wide range of adverse consequences, not only causing suffering in the immediate aftermath, but also long-term detrimental effects to mental and physical health. For example, children with a history of CM are more prone to smoking and obesity, as well as of being perpetrators and victims of violence (WHO, 2016, November). Also, both prospective and retrospective reports of maltreatment were found to be associated with adult psychiatric disorders in a recent study, though the strongest associations were found when maltreatment was retrospectively self-reported (Newbury et al., Reference Newbury, Arseneault, Moffitt, Caspi, Danese, Baldwin and Fisher2018). In that study, it was also shown that young adults who recall being maltreated have a particularly elevated risk for psychopathology. Notably, CM is one of the strongest factors in the development of major depressive disorder (MDD) (Bernet and Stein, Reference Bernet and Stein1999), the leading cause of disability worldwide according to the World Health Organization, with increasing rates over the past decade (WHO, 2016, November).
As both CM and MDD have a high incidence in the general population, the interplay between these two phenomena is important to investigate, both for prevention and treatment. Depressed patients with CM, for example, respond more poorly to antidepressant treatment than those without CM (Nanni et al., Reference Nanni, Uher and Danese2012). CM and MDD may be causally linked, as MDD is a disorder characterized by pathological responses to stress (Frodl et al., Reference Frodl, Koutsouleris, Bottlender, Born, Jager, Scupin, Reiser, Moller and Meisenzahl2008). In experimental studies, chronic social stress induces glucocorticoid-mediated pyramidal dendrite retraction in the hippocampus and changes in dendrite arborization in the prefrontal cortex (PFC) (Woolley et al., Reference Woolley, Gould, Frankfurt and McEwen1990; Magarinos et al., Reference Magarinos, McEwen, Flugge and Fuchs1996; Wellman, Reference Wellman2001; Kole et al., Reference Kole, Czeh and Fuchs2004), which may be associated with the behavioral manifestations of stress-related disorders like MDD (Macqueen and Frodl, Reference Macqueen and Frodl2010). Therefore, based on the extant literature, one hypothesis is that CM in humans acts as a chronic stressor contributing to changes of brain structure and function, which in turn may increase vulnerability to psychiatric disorders such as MDD. Supporting this theory, CM was found to be associated with reduced brain volumes in the amygdala, PFC, and cerebellum (Frodl et al., Reference Frodl, Reinhold, Koutsouleris, Reiser and Meisenzahl2010; Edmiston et al., Reference Edmiston, Wang, Mazure, Guiney, Sinha, Mayes and Blumberg2011; Dannlowski et al., Reference Dannlowski, Stuhrmann, Beutelmann, Zwanzger, Lenzen, Grotegerd, Domschke, Hohoff, Ohrmann, Bauer, Lindner, Postert, Konrad, Arolt, Heindel, Suslow and Kugel2012; Teicher et al., Reference Teicher, Samson, Anderson and Ohashi2016) – regions also reported to be affected in MDD (Wise et al., Reference Wise, Radua, Via, Cardoner, Abe, Adams, Amico, Cheng, Cole, de Azevedo Marques Perico, Dickstein, Farrow, Frodl, Wagner, Gotlib, Gruber, Ham, Job, Kempton, Kim, Koolschijn, Malhi, Mataix-Cols, McIntosh, Nugent, O'Brien, Pezzoli, Phillips, Sachdev, Salvadore, Selvaraj, Stanfield, Thomas, van Tol, van der Wee, Veltman, Young, Fu, Cleare and Arnone2016; Schmaal et al., Reference Schmaal, Hibar, Samann, Hall, Baune, Jahanshad, Cheung, van Erp, Bos, Ikram, Vernooij, Niessen, Tiemeier, Hofman, Wittfeld, Grabe, Janowitz, Bulow, Selonke, Volzke, Grotegerd, Dannlowski, Arolt, Opel, Heindel, Kugel, Hoehn, Czisch, Couvy-Duchesne, Renteria, Strike, Wright, Mills, de Zubicaray, McMahon, Medland, Martin, Gillespie, Goya-Maldonado, Gruber, Kramer, Hatton, Lagopoulos, Hickie, Frodl, Carballedo, Frey, van Velzen, Penninx, van Tol, van der Wee, Davey, Harrison, Mwangi, Cao, Soares, Veer, Walter, Schoepf, Zurowski, Konrad, Schramm, Normann, Schnell, Sacchet, Gotlib, MacQueen, Godlewska, Nickson, McIntosh, Papmeyer, Whalley, Hall, Sussmann, Li, Walter, Aftanas, Brack, Bokhan, Thompson and Veltman2017). A correlation between CM and medial prefrontal gray matter volume was also detected irrespective of diagnosis with MDD or anxiety (van Harmelen et al., Reference van Harmelen, van Tol, van der Wee, Veltman, Aleman, Spinhoven, van Buchem, Zitman, Penninx and Elzinga2010). Finally, in an ENIGMA-MDD mega-analysis focusing on subcortical structures, CM was found to be associated with lower caudate volume in females. Those alterations were more strongly associated with emotional and physical neglect than with other forms of CM (Frodl et al., Reference Frodl, Janowitz, Schmaal, Tozzi, Dobrowolny, Stein, Veltman, Wittfeld, van Erp, Jahanshad, Block, Hegenscheid, Volzke, Lagopoulos, Hatton, Hickie, Frey, Carballedo, Brooks, Vuletic, Uhlmann, Veer, Walter, Schnell, Grotegerd, Arolt, Kugel, Schramm, Konrad, Zurowski, Baune, van der Wee, van Tol, Penninx, Thompson, Hibar, Dannlowski and Grabe2017).
Research in animals and humans also suggests important distinctions between types of CM on brain structure. Specifically, researchers theorized that experiences characterized by deprivation (e.g. emotional and physical neglect) compared with experiences characterized by threat (e.g. emotional abuse and physical violence) lead to different effects on neuronal development (McLaughlin et al., Reference McLaughlin, Sheridan and Lambert2014a). A community study in 287 adolescents showed that exposure to threat and violence was associated with automatic emotion regulation deficits, but not cognitive control disturbances. In contrast, exposure to poverty was associated with worse cognitive control, but no deficits automatic emotion regulation. On the other hand, both violence and poverty predicted poor inhibition in an emotional context (Lambert et al., Reference Lambert, King, Monahan and McLaughlin2017). Interestingly, children exposed to severe deprivation in the form of institutional rearing exhibited widespread cortical thinning in the superior and inferior parietal cortex (McLaughlin et al., Reference McLaughlin, Sheridan, Winter, Fox, Zeanah and Nelson2014b), and children exposed to neglect often have deficits in language abilities (Farah et al., Reference Farah, Shera, Savage, Betancourt, Giannetta, Brodsky, Malmud and Hurt2006). Individuals with a history of deprivation showed smaller gray matter volumes compared with individuals with a history of abuse in the fusiform gyrus and the middle occipital gyrus (Everaerd et al., Reference Everaerd, Klumpers, Zwiers, Guadalupe, Franke, van Oostrom, Schene, Fernandez and Tendolkar2016). Therefore, exploring effects from different types of CM on brain structure was an important goal of the current study.
In this mega-analysis, we first aimed to investigate the association between CM severity and cortical brain structure in MDD patients and healthy subjects. We hypothesized that more severe CM would be related to lower cortical thickness and surface area, especially of the OFC, ACC, medial PFC, and insula – regions affected in adult MDD (Fischl et al., Reference Fischl, Salat, Busa, Albert, Dieterich, Haselgrove, van der Kouwe, Killiany, Kennedy, Klaveness, Montillo, Makris, Rosen and Dale2002) and involved in emotion regulation (Desikan et al., Reference Desikan, Segonne, Fischl, Quinn, Dickerson, Blacker, Buckner, Dale, Maguire, Hyman, Albert and Killiany2006). We also hypothesized that MDD patients with a more severe history of CM would show smaller cortical brain measures than healthy controls (HC) with a similar history of CM. Prior studies detected effects of CM on dorsomedial PFC volume irrespective of diagnosis, but did not fully consider the severity of CM. Besides severity, we also investigated the relationship between different types of CM and brain structure. Furthermore, we hypothesized changes to be more prominent in females than males (Frodl et al., Reference Frodl, Janowitz, Schmaal, Tozzi, Dobrowolny, Stein, Veltman, Wittfeld, van Erp, Jahanshad, Block, Hegenscheid, Volzke, Lagopoulos, Hatton, Hickie, Frey, Carballedo, Brooks, Vuletic, Uhlmann, Veer, Walter, Schnell, Grotegerd, Arolt, Kugel, Schramm, Konrad, Zurowski, Baune, van der Wee, van Tol, Penninx, Thompson, Hibar, Dannlowski and Grabe2017) and thus investigated the interactions between sex and CM on brain structure. Finally, given the large sample size and wide age range, we aimed to explore the interactive effects of CM with age on brain structure.
Methods and materials
Samples
In the current study, 12 international sites participating in the ENIGMA MDD Workgroup with information on CM agreed to participate in the Childhood Adversity Subgroup. Detailed demographics and clinical characteristics for each sample may be found in eTables 1 and 2. Most studies used SCID-1, CIDI, or another form of a standardized interview (eTable 3). Exclusion criteria for study enrollment are given in eTable 3. In total, we analyzed data from 3872 participants: 1284 patients with a lifetime history of MDD and 2588 HC. All participating sites obtained approval from local institutional review boards and ethics committees. In addition, this mega-analysis was approved by the ethics board of the medical faculty of the Otto von Guericke University Magdeburg, Germany. All study participants provided written consent at their local site. In the case of adolescent participants, parent/legal guardian provided written consent and the adolescent provided written assent.
Assessment
Severity of CM (CM-severity) was measured across all sites with the Childhood Trauma Questionnaire (CTQ) (Bernstein et al., Reference Bernstein, Fink, Handelsman, Foote, Lovejoy, Wenzel, Sapareto and Ruggiero1994). The short form of the CTQ is a standardized self-report instrument consisting of 28 items containing five major subscales of CM. Each one also features a cut-off to determine the presence of emotional (⩾12), physical (⩾10), and sexual abuse (⩾8) or emotional (⩾15) and physical neglect (⩾10). Three additional items provide information on responders' tendencies toward minimization and denial. For our analyses, we assessed CM in two ways. First, based on a score above the cut-off for at least one of the abuse or neglect subscales, we divided our participants into four groups (CM-type): no CM, neglect (no abuse), abuse (no neglect), abuse + neglect. In a second analysis, we explored the effect of CTQ total sum score as a continuous variable (CM-severity).
Severity of depressive symptoms at the time of scanning was measured in some sites with the Hamilton Depression Questionnaire (HDRS-17), in others, the Beck Depression Inventory (BDI-II) or Inventory of Depressive Symptomatology-Self Report (IDS-SR) was used. Age of onset and antidepressant medication use at the time of scan were also recorded in 11 and 12 sites, respectively.
Image processing and analysis
Participants all underwent structural T1-weighted MRI brain scans locally at each site, where scans were analyzed using the fully-automated and validated segmentation software FreeSurfer (version 5.0 or higher) (Fischl et al., Reference Fischl, Salat, Busa, Albert, Dieterich, Haselgrove, van der Kouwe, Killiany, Kennedy, Klaveness, Montillo, Makris, Rosen and Dale2002). Image acquisition parameters and software descriptions for each sample are given in eTable 4. Deep brain structure volumes were extracted and visually inspected for segmentation accuracy. Parcellations for cortical thickness and surface area of 68 (34 left and 34 right) regions based on the Desikan–Killiany atlas (Desikan et al., Reference Desikan, Segonne, Fischl, Quinn, Dickerson, Blacker, Buckner, Dale, Maguire, Hyman, Albert and Killiany2006) and left and right hemisphere measures were derived and visually inspected for accuracy following a protocol designed to facilitate harmonized image analysis across multiple sites (http://enigma.ini.usc.edu/protocols/imaging-protocols/). Association between CM and subcortical measures were previously published (Frodl et al., Reference Frodl, Janowitz, Schmaal, Tozzi, Dobrowolny, Stein, Veltman, Wittfeld, van Erp, Jahanshad, Block, Hegenscheid, Volzke, Lagopoulos, Hatton, Hickie, Frey, Carballedo, Brooks, Vuletic, Uhlmann, Veer, Walter, Schnell, Grotegerd, Arolt, Kugel, Schramm, Konrad, Zurowski, Baune, van der Wee, van Tol, Penninx, Thompson, Hibar, Dannlowski and Grabe2017).
Statistical framework of mega-analysis
Statistical analyses were performed using SPSS statistics (version 24).
We performed ANOVAs, or Kruskal–Wallis tests as appropriate, to compare at scan, age at MDD onset, clinical severity of depression, and CTQ scores between groups and cohorts. The χ2 tests were used to analyze differences between frequencies of males and females.
Then, we built generalized estimating equation (GEE) models with thickness or surface area of each region as the dependent variable. Our models had a linear scale response. All participants were included, irrespective of diagnosis. The independent between-subject variable CM was defined in two ways: as the factor CM-type (0 = no CM, 1 = neglect, 2 = abuse, 3 = neglect + abuse) or as CM-severity (continuous: total CTQ). Each of these variables was included in separate models. In all models, we included the between-subjects factors diagnosis (factor: 1 = patients, 0 = HC), sex (factor: 0 = males, 1 = females), and the within-subject factor hemisphere (left, right). Age (continuous), neuroimaging cohort (factor), and total intracranial volume (continuous) were used as between-subject covariates. FreeSurfer version and scanner type were comprised in the factor neuroimaging cohort. As we did not expect CM-severity effects to be lateralized, hemisphere was only included as a main effect. Our prior research showed differential effects of CM-severity in predicting the volume of subcortical structures depending on sex and MDD diagnosis (Frodl et al., Reference Frodl, Janowitz, Schmaal, Tozzi, Dobrowolny, Stein, Veltman, Wittfeld, van Erp, Jahanshad, Block, Hegenscheid, Volzke, Lagopoulos, Hatton, Hickie, Frey, Carballedo, Brooks, Vuletic, Uhlmann, Veer, Walter, Schnell, Grotegerd, Arolt, Kugel, Schramm, Konrad, Zurowski, Baune, van der Wee, van Tol, Penninx, Thompson, Hibar, Dannlowski and Grabe2017). Therefore, we explored in our models all possible interactions between CM, sex, and diagnosis for both surface and thickness. To assess the effect of CM across all brain areas, we first ran analyses on the total thickness and surface across all regions respectively, adding region as a within-subject factor. Then, we repeated the process for each region individually. Finally, we explored the interaction between age and severity of CM while keeping all other terms in the model as the main effects.
In all models, Wald χ2 tests were used to assess the significance of each term. To account for multiple tests (34 regions), a false discovery rate (FDR) correction was computed on the resulting p-values. Findings were considered significant if p FDR<0.05. Any significant interactive effects resulting from the models described above were followed up with post-hoc testing.
Investigation of clinical confounds
A subset of our MDD cohort (N = 966) had more detailed clinical information and allowed us to explore additional potential confounding effects. Therefore, we investigated if the thickness or surface of all regions was significantly predicted by clinical severity (continuous: BDI total score, since HAM-D was available only for a minor subset of participants), recurrence (factor: 0 = first episode, 1 = recurrent episode), current antidepressant use (factor: 0 = no, 1 = yes), age of depression onset (continuous). Information about clinical remission was not provided by all sites; therefore, for this analysis, we defined it as current BDI⩽12 (Riedel et al., Reference Riedel, Möller, Obermeier, Schennach-Wolff, Bauer, Adli, Kronmüller, Nickel, Brieger, Laux, Bender, Heuser, Zeiler, Gaebel and Seemüller2010). We then built GEE models that featured these measures as predictors (main effects), together with CM (severity or type), age, sex, site, hemisphere, and total intracranial volume. Detailed information on this subsample is presented in eTable 14.
Results
Demographics
For details on participant's demographics and clinical features, see Table 1. Overall, data significantly differed between centers with respect to sex, age, CTQ scores, and clinical features (see eTable 5). Frequency of co-occurrence between abuse and neglect is 12.7%. CM severity is influenced by abuse and neglect to a similar extent (β = 0.41 and β = 0.42, respectively).
CTQ, Childhood Trauma Questionnaire; ICV, total intracranial volume; BDI#, Beck Depression Inventory; HDRS-17#, Hamilton Depression Rating Scale.
Shown are mean values±standard deviation. For CTQ subscales, number of subjects above the cut-off are given in brackets.
a Mann Whitney U test used.
b From sites that used these ratings.
Cortical thickness
Main effects of CM-severity
A summary of all significant findings is reported in Table 2. For an overview of the results of the models run on each region, see eTable 6.
Wald χ2 and p values of CM severity and type are shown for the regions where they were significant. For all effects see Supplemental Tables. Effects are coefficients for the model term CM severity or the estimates of the indicated contrast for CM type. CM, childhood maltreatment; FDR, false discovery rate.
We detected a significant main effect showing an inverse relation between CM-severity and thickness of the banks of the superior temporal sulcus (Wald χ2 = 14.583, p FDR = 0.033, B = −0.001, Fig. 1). A significant main effect of CM was also present on the thickness of the supramarginal gyrus (Wald χ2 = 8.889, p FDR = 0.049, B = −0.001, Fig. 1).
CM-severity and sex interaction
When considering all regions, the interaction between CM-severity and sex was significant (Wald χ2 = 5.220, p = 0.022). Dividing the data by sex, post-hoc analyses showed a significant negative effect of CM-severity on cortical thickness in females (Wald χ2 = 4.861, p = 0.027, B = −0.000649), but not in males (Wald-χ2 = 1.287, p = 0.257, B = −0.000136).
When running models for each region separately, we found a significant interaction between CM-severity and sex on the cortical thickness of the rostral anterior cingulate cortex (Wald χ2 = 13.556, p FDR = 0.008). Post-hoc analysis revealed a significant positive effect of CM-severity on the cortical thickness of this region in males (Wald χ2 = 14.426, p < 0.001, B = 0.002, Fig. 1) but not in females (Wald χ2 = 3.174, p = 0.075, B = −0.0006).
CM-severity and age interaction
When considering all regions, a significant interaction between age and severity of CM was detected (Wald χ2 = 11.105, p = 0.001, B = −0.000035).
Models ran for each region separately indicated that this interaction between age and severity of CM was significant across all participants in the rostral anterior cingulate, isthmus of the cingulate, posterior cingulate, lateral orbitofrontal gyrus, parahippocampal gyrus, inferior frontal gyrus (IFG) pars opercularis, IFG pars triangularis, superior frontal gyrus, banks of the superior temporal sulcus, cuneus, fusiform gyrus, insula, precentral gyrus, precuneus, supramarginal gyrus, and transverse temporal gyrus (see Table 3, eTable 7, Fig. 2).
Wald χ2 and p values of CM severity × Age and CM type × Age are shown for the regions where they were significant. For all effects see Supplemental Tables. Effects are coefficients for the model term CM severity × Age or those for Age in the Neglect + Abuse group v. the no-CM group. CM, childhood maltreatment; FDR, false discovery rate.
Main effects of CM-type
We found a significant main effect of CM type (eTable 10) in the banks of the superior temporal sulcus (Wald χ2 = 19.888, p FDR = 0.006), inferior parietal lobe (Wald χ2 = 15.273, p FDR = 0.023), middle temporal lobe (Wald χ2 = 12.123, p FDR = 0.048), precuneus (Wald χ2 = 15.325, p FDR = 0.023), and supramarginal gyrus (Wald χ2 = 13.990, p FDR = 0.026). In all cases, the neglect + abuse group had lower mean thickness values compared to the no CM group (all p < 0.01, Fig. 1) and there was no difference between the abuse only as well as neglect only CM types and the no CM group.
CM-type and age interaction
The interaction between age and type of CM was significant across all participants for most regions (see Table 3, eTable 11). In all cases, the effects of age were more negative in the neglect + abuse group compared to the CM group (all p < 0.05).
Cortical surface area
A summary of all significant findings is reported in Table 2. For an overview of the results of the models run on each region, see eTable 8.
Main effects of CM-severity
Across all regions, a negative main effect of CM-severity on the cortical surface area was observed (Wald χ2 = 4.413, p = 0.036, B = −0.414). When running separate models for each region (eTable 8), we detected a significant inverse main effect of CM-severity on the surface area of the middle temporal gyrus (Wald χ2 = 12.368, p FDR = 0.015, B = −1.504, Fig. 3).
CM-type, diagnosis, and sex interaction
We found a significant interaction between CM type, diagnosis, and sex (eTable 12) in the caudal anterior cingulate (Wald χ2 = 17.807, p FDR < 0.001). Post-hoc testing revealed that, in depressed males, those having suffered from either abuse or neglect had a lower average cortical surface area of the caudal anterior cingulate cortex than those who had no history of CM (p = 0.003 and p = 0.017, respectively, Fig. 3).
CM-severity/type and age interaction
We found no effects of the interaction between CM-severity or CM type and age in predicting cortical surface area (eTable 9).
Investigation of clinical confounds
Our post-hoc investigation in a subset of patients with detailed information showed no significant effects of clinical variables on thickness or surface (all p > 0.05). See eTables 15–18 for the model effects.
Discussion
This study represents the largest effort worldwide to investigate the association between CM and cortical brain structure in a large sample of MDD patients and healthy subjects. We found that CM has a subtle but widespread association with cortical thickness and surface area, which is likely influenced by sex and age.
Two procedures of describing CM were used. First, the dimensional measure of CM severity allowed for continuous analysis. Second, the categorical classification of CM into no CM, only neglect, only abuse, and both abuse and neglect, allowed for an analysis of the type of CM. It should be highlighted that participants exposed to both neglect and abuse also had higher total CM values.
Severity of CM was associated with lower mean cortical surface area regardless of the region across all MDD patients and HC. In women, higher severity of CM was also associated with thinner thickness across all regions. These findings are consistent with prior research showing widespread effects of severity of CM on the brain irrespective of psychopathological status (Chaney et al., Reference Chaney, Carballedo, Amico, Fagan, Skokauskas, Meaney and Frodl2014).
Regardless of diagnosis, when individual regions were investigated, this effect of CM severity survived correction for multiple testing in temporal and temporo-parietal regions. Specifically, participants with higher CM severity had significantly thinner cortex in the banks of the superior temporal sulcus and the supramarginal gyrus. When considering the type of abuse, high severity of CM, represented by concurrent childhood neglect and abuse, was once again associated with reduced thickness in these two areas and additionally in the precuneus, middle temporal lobe, and inferior parietal cortex. Moreover, participants with higher severity of CM also showed a smaller surface area of the middle temporal gyrus. The magnitude of these negative effects on thickness and surface area pointed toward a reduction around 0.001 mm of thickness and around 0.4–1.5 mm2 of surface with each one-point increase in CTQ score depending on the brain region. For example, an increase of 100 points in the CTQ scale would be associated with a 0.1 mm (or 4%) decrease of cortical thickness in the banks of the superior temporal sulcus. The middle temporal lobe is thought to be essential for our ability to understand actions and semantic associations (Davey et al., Reference Davey, Thompson, Hallam, Karapanagiotidis, Murphy, De Caso, Krieger-Redwood, Bernhardt, Smallwood and Jefferies2016). One possibility is that CM may lead to difficulties in semantic retrieval through alterations in regions of temporal cortex and the default mode network. Indeed, other studies have also demonstrated that individuals with higher severity of CM showed reduced cortical surface on the left middle temporal area and lingual gyrus (Kelly et al., Reference Kelly, Viding, Wallace, Schaer, De Brito, Robustelli and McCrory2013). In contrast, in a study of adolescents and young adults exposed to CM, increased cortical volume was observed in the left inferior and middle temporal gyri relative to HC (Lim et al., Reference Lim, Hart, Mehta, Worker, Simmons, Mirza and Rubia2018). In the present study, the other regions we report showing an impact from CM type are located in the temporo-parietal area and around the temporo-parietal junction; both of these regions play a role in the theory of mind processing which is important during daily social interactions (Saxe and Kanwisher, Reference Saxe and Kanwisher2003). Deficits in these areas might suggest a disadvantage for subjects with a history of CM, in particular those with increased severity and more types of CM. No significant main effects of severity of CM were detected for other regions we hypothesized to be vulnerable. For example, prior studies found a significant main effect of CM in different regions such as fronto-limbic areas, visual cortex, and cerebellum (Kelly et al., Reference Kelly, Viding, Wallace, Schaer, De Brito, Robustelli and McCrory2013; Yang et al., Reference Yang, Cheng, Mo, Bai, Shen, Liu, Li, Jiang, Chen, Lu, Sun and Xu2017). This might be due to smaller sample sizes and more homogeneity in prior studies: in our analysis, which features a larger sample size of N = 3872, we detected an overall effect of severity of CM on the whole cortex with some prominence in the temporal and temporo-parietal regions.
Another interesting finding was that males, but not females, with a more severe history of CM, regardless of diagnosis, showed distinctly thicker rostral anterior cingulate cortices. These results suggest sex differences in the effects of CM on the structure of this region (Fallucca et al., Reference Fallucca, MacMaster, Haddad, Easter, Dick, May, Stanley, Rix and Rosenberg2011; Canu et al., Reference Canu, Kostic, Agosta, Munjiza, Ferraro, Pesic, Copetti, Peljto, Lecic Tosevski and Filippi2015). The anterior cingulate cortex is involved in emotional and inhibitory processes (Garavan et al., Reference Garavan, Hester, Murphy, Fassbender and Kelly2006; Steele et al., Reference Steele, Aharoni, Munro, Calhoun, Nyalakanti, Stevens, Pearlson and Kiehl2013). Thus, males seem to be particularly sensitive to CM with regards to the thickness in a region relevant to emotion regulation and might show a reactive increase of thickness. Whether this thickness increase of the rostral ACC is adaptive cannot be effectively addressed in the current sample, as longitudinal data and resiliency measures were not available for the bulk of the cohort. In this context, it is interesting that male patients with a history of neglect and abuse had significantly smaller surface areas of the caudal ACC compared to those without CM, pointing toward a negative effect of CM in the caudal ACC in participants who developed MDD. In healthy subjects, such an effect of CM was not seen in the caudal ACC and we could tentatively ascribe to resilience (Feder et al., Reference Feder, Nestler and Charney2009).
Childhood maltreatment by age interaction
A novel finding detected in our secondary analysis was that CM severity and age interacted regardless of diagnosis to predict the thickness of several cortical regions involved in emotional processing, such as portions of the cingulate, orbitofrontal, insular, dorsolateral prefrontal, and medial prefrontal cortices. In these areas, older people including patients and controls with higher CM severity had lower cortical thickness. The orbitofrontal and dorsomedial prefrontal cortices can be seen as separate networks interacting closely with limbic structures, but also showing cortico-cortical interconnections with each other (Ongur and Price, Reference Ongur and Price2000; Phillips et al., Reference Phillips, Ladouceur and Drevets2008). These regions allow the brain to process emotionally salient information and help with the regulation of emotional behavior (Phillips et al., Reference Phillips, Ladouceur and Drevets2008). The insula is closely interconnected with the orbitofrontal cortex and is involved in emotion and executive processing as well as working memory (Levens and Phelps, Reference Levens and Phelps2010). The cingulate cortex is also well known to have cognitive and emotional functions: its dorsal parts are involved in emotion evaluation, whereas the ventral parts and the dorsomedial PFC are involved in emotion regulation (Etkin et al., Reference Etkin, Egner and Kalisch2011). Overall, these results are consistent with previous studies showing that CM-severity impacts regions involved in emotion regulation, including the insula (Teicher et al., Reference Teicher, Anderson, Ohashi and Polcari2014).
Also, our cross-sectional data suggest that cortical thickness might decrease more rapidly with age in individuals with a more severe history of CM, although this still needs to be confirmed by longitudinal analyses. It will be critical for future studies to assess the effects of abuse across multiple time points and to consider such abuse in relation to the age of participants.
Diagnosis and CM interactions
In the present study, no effect of diagnosis was detected and we also did not find a significant interaction between diagnosis and CM. It is possible that we could not replicate the main effect of MDD diagnosis on OFC and ACC thickness because of the smaller sample size of the current study (overall N = 3872) compared to the sample size in our previous ENIGMA MDD meta-analysis that focused on the effects of MDD (overall N = 10 105) (Schmaal et al., Reference Schmaal, Hibar, Samann, Hall, Baune, Jahanshad, Cheung, van Erp, Bos, Ikram, Vernooij, Niessen, Tiemeier, Hofman, Wittfeld, Grabe, Janowitz, Bulow, Selonke, Volzke, Grotegerd, Dannlowski, Arolt, Opel, Heindel, Kugel, Hoehn, Czisch, Couvy-Duchesne, Renteria, Strike, Wright, Mills, de Zubicaray, McMahon, Medland, Martin, Gillespie, Goya-Maldonado, Gruber, Kramer, Hatton, Lagopoulos, Hickie, Frodl, Carballedo, Frey, van Velzen, Penninx, van Tol, van der Wee, Davey, Harrison, Mwangi, Cao, Soares, Veer, Walter, Schoepf, Zurowski, Konrad, Schramm, Normann, Schnell, Sacchet, Gotlib, MacQueen, Godlewska, Nickson, McIntosh, Papmeyer, Whalley, Hall, Sussmann, Li, Walter, Aftanas, Brack, Bokhan, Thompson and Veltman2017). However, since this is the largest joint mega-analysis concerning cortical thickness and CM to date, this null finding could also suggest that the effects of MDD commonly reported in studies and meta-analyses may be the result of the interaction of several underlying variables. For example, different effects of MDD depending on age and onset were earlier highlighted in Schmaal et al., (Reference Schmaal, Hibar, Samann, Hall, Baune, Jahanshad, Cheung, van Erp, Bos, Ikram, Vernooij, Niessen, Tiemeier, Hofman, Wittfeld, Grabe, Janowitz, Bulow, Selonke, Volzke, Grotegerd, Dannlowski, Arolt, Opel, Heindel, Kugel, Hoehn, Czisch, Couvy-Duchesne, Renteria, Strike, Wright, Mills, de Zubicaray, McMahon, Medland, Martin, Gillespie, Goya-Maldonado, Gruber, Kramer, Hatton, Lagopoulos, Hickie, Frodl, Carballedo, Frey, van Velzen, Penninx, van Tol, van der Wee, Davey, Harrison, Mwangi, Cao, Soares, Veer, Walter, Schoepf, Zurowski, Konrad, Schramm, Normann, Schnell, Sacchet, Gotlib, MacQueen, Godlewska, Nickson, McIntosh, Papmeyer, Whalley, Hall, Sussmann, Li, Walter, Aftanas, Brack, Bokhan, Thompson and Veltman2017). Our findings suggest that CM could be another factor that accounts in part for structural differences between depressed patients and HC.
Strengths and limitations
A major strength of the study is the large sample size with a relevant control sample allowing inclusion of all 34 left and right cortical brain regions in the analyses. However, a larger sample from different sites also limits the common information collected, since not all sites used the same assessments. Not all patients were drug free and, further, the history of antidepressants use as well as duration, type, and dosage of antidepressant treatment were not assessed completely during the lifetime, so we cannot rule out effects of antidepressant exposure influencing our results. Additionally, psychopathology was assessed with different ratings, so that it was not possible to use depression severity as a single covariate in the analysis. Overall, the datasets from the different samples included in the mega-analysis were significantly heterogeneous regarding demographics and clinical features. This is a common limitation of multi-site analyses and we accounted for this effect by adding site as a covariate to all of our models. It is worth noting that we were able to analyze the influence of clinical confounds in a subset of our MDD cohort, where we did not find any significant effect of recurrence, antidepressant medication, remission, severity or age of depression onset in predicting cortical thickness or surface. However, measures of socioeconomic status and education have been shown to play a role in brain structure (Ritchie et al., Reference Ritchie, Tucker-Drob, Cox, Dickie, Del, Corley, Royle, Redmond, Munoz Maniega, Pattie, Aribisala, Taylor, Clarke, Gow, Starr, Bastin, Wardlaw and Deary2017), but were unfortunately unavailable in our sample. Finally, we considered hemisphere as a within-subject effect in our dataset. Our hypothesis was that CM would affect anatomically distinct regions differently rather than be selective for a specific region on a specific hemisphere. Therefore, we believed that including all possible interactions between regions and hemisphere would lead to an unnecessarily complex model. However, it is possible that besides the bilateral effects we report, subtler lateralized effects of CM might exist in specific areas.
Even if our investigation features the broad variation of ‘real life’ clinical populations, future studies are needed to confirm our findings in carefully controlled datasets. Here, we explored the effect of CM in a sample of healthy participants and patients with MDD. Because it is not clear how the severity or type of CM may affect the development of structural brain measures, it will be important to consider the onset and timing of CM in future (Ho et al., Reference Ho, Dennis, Thompson and Gotlib2018). In addition, future longitudinal data are required to establish whether cortical thickness might decrease more rapidly with age in individuals with a more severe history of CM, as our current cross-sectional data may suggest. For this analysis, while it was possible to use extracted cortical measures from specific regions of interest, it was not possible to retrospectively analyze the original MRI datasets to perform a whole-cortex analysis with FreeSurfer. A surface-based analysis across the entire cortex may afford more sensitivity in detecting effects of CM and thus could be a future step.
Conclusions
The results of our study support the idea that CM-severity appears to affect the structure of temporal and parietal regions in particular. Thus, there are effects in the default mode network and in regions involved in theory of mind as well as perception. Interestingly, CM may interact with the effect of age on cortical thickness in these regions and others involved in emotion regulation. Thus, future longitudinal studies should investigate if subjects with a history of CM may be more prone to cortical thinning during aging or if CM results in changes that mimic aging.
Supplementary material
The supplementary material for this article can be found at https://doi.org/10.1017/S003329171900093X
Acknowledgements
The ENIGMA-Major Depressive Disorder working group gratefully acknowledges support from the NIH BD2 K award, U54 EB020403. The infrastructure for the NESDA study (http://www.nesda.nl) is funded through the Geestkracht program of the Netherlands Organisation for Health Research and Development (Zon-Mw, grant number 10-000-1002) and is supported by participating universities (VU University Medical Center, GGZ inGeest, Arkin, Leiden University Medical Center, GGZ Rivierduinen, University Medical Center Groningen) and mental health care organizations, see http://www.nesda.nl. Lianne Schmaal is supported by The Netherlands Brain Foundation Grant number F2014(1)-24 and the Neuroscience Campus Amsterdam Scientific Excellence grant (IPB-SE-15-PSYCH-Schmaal).
Financial support
The Münster Neuroimaging Cohort was funded by the German Research Foundation (SFB-TRR58, Projects C09 and Z02 to UD) and the Interdisciplinary Center for Clinical Research (IZKF) of the medical faculty of Münster (grant Dan3/012/17 to UD). FOR2107-Marburg was funded by the German Research Foundation (DFG, grant FOR2107 KI588/14-1, KI588/15-1, KI588/14-2, KI588/15-2 to TK, KR3822/5-1 to AK). FOR2107-Münster was funded by the German Research Foundation (DFG, grant FOR2107 HA 7070/2-2 to TH; DA1151/5-1 and DA1151/5-2 to UD). ENIGMA was supported in part by a Consortium grant (U54 EB020403 to PMT, DPH, NJ) from the NIH Institutes contributing to the Big Data to Knowledge (BD2 K) Initiative, and by R01 MH116147 and R56 AG058854. Moreover, an NIH BD2 K award, U54EB020403 to T.G.M. van Erp. London: This paper represents independent research [part] funded by the National Institute for Health Research (NIHR) Biomedical Research Centre at South London and Maudsley NHS Foundation Trust and King's College London. The views expressed are those of the author(s) and not necessarily those of the NHS, the NIHR or the Department of Health and Social Care. SHIP: The Study of Health in Pomerania (SHIP) is supported by the German Federal Ministry of Education and Research (grants 01ZZ9603, 01ZZ0103, and 01ZZ0403) the Ministry of Cultural Affairs as well as the Social Ministry of the Federal State of Mecklenburg-West Pomerania. MRI scans were supported by Siemens Healthcare, Erlangen, Germany. SHIP-LEGEND was supported by the German Research Foundation (GR1912/5-1). Dublin: The study was supported by a Science Foundation Ireland (SFI) Stokes Professorship Grant to Thomas Frodl. San Francisco: This work was supported by the Brain and Behavior Research Foundation (formerly NARSAD) to T.T.Y.; the National Institute of Mental Health (R01MH085734 to T.T.Y.; K01MH097978 to K.Z.L.; K01MH117442 to T.C.H.) and by the American Foundation for Suicide Prevention (PDF-1-064-13) to T.C.H.
Conflict of interest
Nic van der Wee has consulted for Lilly, Wyeth, Servier, Pfizer, and GlaxoSmithKline. Thomas Frodl has received fees for presentations for educational programs with Servier, Janssen, and Lundbeck. Hans Grabe has consulted for Fresenius Medical Care. Carsten Konrad received fees for an educational program from Aristo Pharma, Janssen-Cilag, Lilly, MagVenture, Servier, and Trommsdorff as well as travel support and speakers honoraria from Aristo Pharma, Janssen-Cilag, Lundbeck, Neuraxpharm, and Servier. Theo G.M. van Erp has consulted for Roche Pharmaceuticals, Ltd., and has a contract with Otsuka Phamaceutical Co., Ltd (OPCJ). Knut Schnell has consulted for Roche Pharmaceuticals and Servier Pharmaceuticals.