Introduction
Experiences of childhood maltreatment (CM) are one of the main contributors to mental illness (Brown, Harris, & Craig, Reference Brown, Harris and Craig2019; Hughes et al., Reference Hughes, Bellis, Hardcastle, Sethi, Butchart, Mikton and Dunne2017). However, CM is non-specifically associated with psychiatric disorders, i.e. several types of CM can increase vulnerability for a specific disorder in different patients (Vachon, Krueger, Rogosch, & Cicchetti, Reference Vachon, Krueger, Rogosch and Cicchetti2015). CM has been associated with early onset of psychiatric illness, increased symptom severity and comorbidity, and poor clinical outcomes characterized by requiring higher medication dosages, increased suicidal behavior, and more and longer hospitalizations (Lippard & Nemeroff, Reference Lippard and Nemeroff2020). Furthermore, factors such as time of exposure, chronicity, and severity of childhood abuse or neglect play a role in clinical outcomes. Studies indicate a dose–response relationship between multiplicity of exposure, severity or frequency, and risk of mental disorders (Anda et al., Reference Anda, Felitti, Bremner, Walker, Whitfield, Perry and Giles2006).
CM is associated with dysregulation of stress-mediating systems, thereby increasing the risk of mental and physical health problems. Specifically, disruptions in hypothalamic–pituitary–adrenal (HPA)-axis regulation have been studied as a potential mediator of this association (Koss & Gunnar, Reference Koss and Gunnar2018; Kudielka & Wüst, Reference Kudielka and Wüst2010). The HPA-axis is one of the main stress response systems; cortisol, its final effector, released in direct response to acute stressors, triggers a wide range of actions by regulating gene transcription and epigenetic modifications in several brain areas (Provençal, Arloth, Cattaneo, Anacker, & Cattane, Reference Provençal, Arloth, Cattaneo, Anacker and Cattane2019). Furthermore, HPA-axis maintains a diurnal rhythm, with the highest cortisol levels in the morning which decrease progressively during the day until reaching the lowest at midnight. Since the HPA-axis continues to mature during early stages of life, environmental factors such as early-life stress may induce long-lasting changes in its functioning, resulting in the emergence of different disorders (Tarullo & Gunnar, Reference Tarullo and Gunnar2006). However, findings regarding alterations in the patterns of cortisol associated with early-life stress have been inconsistent (Fogelman & Canli, Reference Fogelman and Canli2018).
A recent meta-analysis focusing on CM and diurnal HPA-axis activity in children and adults reported no overall effect on diurnal cortisol slope (Bernard, Frost, Bennett, & Lindhiem, Reference Bernard, Frost, Bennett and Lindhiem2017). However, a moderate association was found between CM and blunted awakening cortisol concentrations when considering only sufferers of CM who were referred from child welfare system agencies. In contrast, another recent meta-analysis showed that CM affects HPA-axis reactivity during stressful situations as evidenced by a flattened cortisol pattern during an acute psychosocial stress task in children and adults who faced early-life adversities (Bunea, Szentágotai-t, & Miu, Reference Bunea, Szentágotai-t and Miu2017). Interestingly, the effects were more pronounced in studies focused on adults and CM. These findings suggest a pattern of blunted cortisol response during the peak and recovery phases of acute stress, and overall hypocortisolism in individuals exposed to CM. However, some studies report hypercortisolism in subjects exposed to early-life stress, childhood trauma, or insensitive interactions with caregivers (Hunter, Minnis, & Wilson, Reference Hunter, Minnis and Wilson2011). Besides, it has been suggested that distinct patterns of cortisol responses may be partially explained by CM severity and frequency (Ouellet-Morin et al., Reference Ouellet-Morin, Robitaille, Langevin, Cantave, Brendgen and Lupien2019), pubertal stage (King et al., Reference King, Colich, LeMoult, Humphreys, Ordaz, Price and Gotlib2017), or sex (Trickett, Gordis, Peckins, & Susman, Reference Trickett, Gordis, Peckins and Susman2014).
Notably, HPA-axis dysregulation, both hyperactivity and hypoactivity, has been associated with different psychiatric disorders and other disease outcomes (Turner et al., Reference Turner, Smyth, Hall, Torres, Hussein, Jayasinghe and Clow2020). Although infancy is a sensitive period for HPA-axis regulation, this system remains plastic and it can be recalibrated during specific ontogenic periods, if the environmental conditions improve. In fact, recent studies support puberty as a key recalibration period to trigger shifts in HPA-axis functioning in postinstitutionalized children (DePasquale, Donzella, & Gunnar, Reference DePasquale, Donzella and Gunnar2019).
Thus, the main aim of the current research was to establish the proximal effects of CM on HPA-axis regulation and anxiety perception in children and adolescents, under basal conditions and in response to a psychosocial stressor, as compared with youth without CM. In addition, the differential impact of the severity and frequency of the CM experiences was also analyzed to better dissect the relationship between CM and HPA-axis dysfunction. Finally, anxiety perception was assessed throughout the experimental stress paradigm to verify that all participants underwent a subjective experience of stress (regardless of their CM history); thus, the potential differences in stress perception with regard to CM can be disentangled from actual differences in HPA-axis functioning. Complementarily, anxiety trait was also assessed in relationship with basal diurnal cortisol output. Specifically, we hypothesized that exposure to CM would be associated with blunted HPA-axis functioning and higher anxiety perception. Moreover, more severe and frequent exposure to CM would be associated with greater dysregulation of the HPA-axis following a dose–response relationship.
Methods
The EPI-Young-Stress project is a multi-center study which aims to evaluate HPA-axis functioning, associated epigenetic signatures, and immunological biomarkers involved in the association between CM and youth mental disorders. The research was conducted at the University of Barcelona and six child and adolescent psychiatry departments in Spain: Hospital Benito Menni, Hospital Clínic Barcelona, Hospital Gregorio Marañón, Hospital Puerta de Hierro, Hospital Santiago Apóstol, and Day Hospital Orienta Gavà.
The study was approved by the Ethical Review Board of each participating hospital and university. Families were explicitly informed about the voluntary nature of the study, their rights, and the procedures, risks, and potential benefits involved. Written consent was required from all parents or legal guardians; the children provided written assent after the nature of the procedure had been fully explained.
Participants
A total of 187 children and adolescents aged 7–17 years participated in this study. Children without psychopathology were recruited from advertisements, primary healthcare centers, schools, and other community facilities. Children with current psychopathology were recruited from the above-mentioned hospitals (inpatient clinics, partial hospitalization programs, and outpatient clinics) (see Table 1). Recruitment lasted from April 2016 to March 2020. Exclusion criteria for all participants included diagnosis of autism spectrum disorder, eating disorder with body mass index (BMI)<18.5, intellectual disability (IQ < 70), current drug dependence, non-fluency in Spanish, extreme premature birth (<1500 g at birth), head injury with loss of consciousness, and severe neurological or other pathological conditions likely to affect HPA-axis functioning (such as cancer or autoimmune diseases).
ADHD, attention-deficit/hyperactivity disorder; BMI, body mass index; CGAS, Children's Global Assessment Scale, rating from 1 to 100 with higher ratings indicating better functioning in a wide range of activities; CM, childhood maltreatment (CM group refers to subjects with a confirmed or suspected history of CM); GFES, The Global Family Environment Scale, ranging from 1 to 90, with higher scores indicating a better family environment; SES, socioeconomic status, raw scores range from 8 to 66, with higher scores reflecting higher SES; WHR, waist-to-hip ratio.
a Student's t test.
b χ2 test.
c Other ethnicities included Latin American (66%), Maghrebin (16%), sub-Saharan (9%), and others (9%).
dThis analysis was conducted with 183 subjects.
eThis analysis was only conducted with the 116 subjects with a current psychiatric diagnosis.
fThis analysis was only conducted with the 108 female subjects.
gThis analysis was conducted with 176 subjects.
hThis analysis was conducted within 171 subjects.
p values: *p < 0.05, **p < 0.01, and ***p < 0.001. d = Cohen's effect size.
Procedures
Sociodemographic and clinical measures
The interview package included basic demographic information including socioeconomic status (SES) based on the Hollingshead Four-Factor Index of SES (Hollingshead, Reference Hollingshead1975). Pubertal development was assessed using the Tanner staging questionnaire (Morris & Udry, Reference Morris and Udry1980) and participants were classified as either children (Tanner stages 1–3) or adolescents (Tanner stages 4–5). The Global Family Environment Scale (GFES) was used to measure the quality of the family environment (Rey et al., Reference Rey, Singh, Hung, Dossetor, Newman, Plapp and Bird1997). Additionally, ethnicity, BMI and waist-to-hip ratio were recorded.
Both participants and their parents directly recounted the youth's medical history. Psychopathology was assessed using the Spanish version of the Schedule for Affective Disorders and Schizophrenia for School-Age Children Present and Lifetime Version DSM-5 (K-SADS-PL-5) (de la Peña et al., Reference de la Peña, Villavicencio, Palacio, Félix, Larraguibel, Viola and Ulloa2018). Information was completed whenever possible using medical records. Final diagnoses were established by consensus, and based on DSM-5 criteria (APA: American Psychiatric Association, 2013), primary psychiatric diagnoses were later classified into dimensions to better characterize the sample (see Table 1). The global level of functioning was measured by the Children's Global Assessment Scale (CGAS) (Shaffer et al., Reference Shaffer, Gould, Brasic, Ambrosini, Fisher, Bird and Aluwahlia1983). The use of psychiatric medication was dichotomized as absence/presence, since there were no differences in cortisol levels according to the different drugs (data available upon request). Current illegal drug use was classified into four frequency groups: never, less than once a month, once a month or more, and daily use (Forti et al., Reference Forti, Quattrone, Freeman, Tripoli, Gayer-anderson and Quigley2019).
Childhood maltreatment assessment
All participants and their parents/legal guardians were interviewed separately, face to face, by one trained psychologist or psychiatrist. They were assessed by means of an exhaustive interview focused on the identification of signs of child vulnerability, adverse experiences, and family interactions, based on the criteria of the instrument ‘Tool for assessing the severity of situations in which children are vulnerable’ (TASSCV), which has been validated by professionals working in child and adolescent care units (see online Supplementary material) (CARM, 2012). Additionally, adolescents older than 12 were assessed for history of CM via the short version of the Childhood Trauma Questionnaire (CTQ-SF) (Bernstein et al., Reference Bernstein, Stein, Newcomb, Walker, Pogge, Ahluvalia and Zule2003) and the Childhood Experience of Care and Abuse Questionnaire (CECA-Q2) (Kaess et al., Reference Kaess, Parzer, Mattern, Resch, Bifulco and Brunner2011). Children under 12 years answered an adapted hetero-administered ad hoc questionnaire (see online Supplementary material). Afterwards, the clinicians completed a table summarizing the different forms of CM effected by caregivers or other adults (not by peers), being TASSCV the main measure of CM used in the primary analyses, while the other measures (CTQ-SF, CECA-Q2, ad hoc questionnaire, and reports from social services or teachers) were used as an additional source of information for the clinicians. The exhaustive participants' evaluation during the recruitment process allowed for clinicians to enrich their praxis. In addition, after the interviews of this study, a referral system of urgent appointment was implemented for those subjects who requested it, activating the usual protocols that guarantee the children's protection rights. Following the TASSCV criteria, each CM type was coded as either: (i) absent, (ii) suspected (if significant signs of neglect or abuse emerged during the evaluation), or (iii) confirmed (with clear evidence from social services or family). Severity and frequency of different types of CM were rated on a four-point Likert scale according to TASSCV criteria. CM severity was coded according to the characteristics of the experience suffered as low (1), moderate (2), severe (3), or very severe (4); while frequency was coded as whether CM had occurred once (1), sometimes (2), often (3), or frequently (4). Five types of CM were considered in the following analysis: physical neglect, emotional neglect, physical abuse, emotional abuse, and sexual abuse.
HPA-axis functioning
Four saliva samples were collected during a normal day with the aim to assess HPA-axis diurnal functioning (basal condition), specifically, on waking up (B1), 30 min after waking (B2), before lunch (B3), and before bedtime (B4). On a different day, in order to explore HPA-axis reactivity during acute psychosocial stress, the Trier Social Stress Test for children (TSST-C), a validated protocol that reliably induces HPA-axis activation, was applied (Buske-Kirschbaum et al., Reference Buske-Kirschbaum, Jobst, Wustmans, Kirschbaum, Rauh and Hellhammer1997). Briefly, upon arrival at the lab, the participants waited in a quiet room for 30 min before entering the examination room, where a panel of judges awaited. During the 20 min of the stress situation, the participants had to perform a speaking and an arithmetic task following instructions from the judges while being videotaped. After the stress task, the participants returned to the first room for 30 min (see online Supplementary material for a more detailed description of the procedure). Five saliva samples were collected during this procedure: 30 min before the stressor (T1), immediately before the stressor (T2), immediately after the stressor (T3), 15 min after the stressor (T4), and 30 min after the stressor (T5) (see Fig. 1). All participants were scheduled at 16:00 h to control for diurnal cortisol variability. Previously, further instructions were given to the participants to avoid factors that have been reported to influence cortisol levels (details in the online Supplementary material). Details about collection time of each salivary cortisol sample are available in Table 2.
AUCg, area under the curve with respect to ground (indicating the total cortisol output); AUCi, area under the curve with respect to increase (reflecting cortisol changes over time); CM, childhood maltreatment (CM group refers to the subjects with a confirmed or suspected history of CM based on TASSCV criteria); STAI-State (PC), percentile scores of state anxiety inventory scale (for adolescents 16–17 years old) and state anxiety inventory for children scale (for participants under 15); STAI-Trait (PC), percentile scores of anxiety trait inventory scale (for adolescents 16–17 years old) and anxiety trait inventory for children scale (for participants under 15); TSST-C, Trier Social Stress Test for children.
Diurnal salivary cortisol was measured at: B1, immediately after awakening; B2, 30 min after waking; B3, before lunch; B4, before bedtime. Mean time for saliva sample collection: 08:52 ± 1:27 (6:00–12:00) (B1); 09:24 ± 1:26 (6:30–12:59) (B2); 14:19 ± 0:53 (12:15–16:40) (B3); and 22:37 ± 01:16 [20:00–2:50(+1day)] (B4). Saliva samples for cortisol measurement during TSST-C were collected at: T1, 30 min before stressor; T2, immediately before stressor; T3, immediately after stressor; T4, 15 min after stressor; T5, 30 min after stressor. Mean time for saliva sample collection during the TSST-C procedure: 16:04 ± 0:11 (15:13–17:15) (T1); 16:33 ± 0:12 (15:42–17:45) (T2); 16:53 ± 0:13 (15:59–18:00) (T3); 17:08 ± 0:13 (16:08-18:16) (T4); and 17:23 ± 0:13 (16:30–18:30)(T5).
Dichotomous CM refers to the analysis comparing youth without CM with youth exposed to any type of CM. Severity of CM refers to the analysis comparing youth without CM, youth exposed to low CM, and youth exposed to moderate/severe CM. Frequency of CM refers to the analysis comparing youth without CM, youth exposed to CM once/sometimes, and youth exposed to CM often/frequently.
a Mixed-effects model (for single measurements) and ANOVA (for AUCg and AUCi). The analyses include the following covariates: clinical status, sex, pubertal stage, psychopharmacological treatment, illegal drugs use, oral contraceptive use, corticosteroid medication, ethnicity, SES, and BMI [additionally adjusting by the time of the first cortisol sample collection (B1) for diurnal analysis]. Values in superscript (B4, T3, T4, T5) indicate the samples with a significant difference in the simple effects test in the context of mixed-effect model.
p values: *p < 0.05, **p < 0.01, and ***p < 0.001. p ⩽ 0.006 [as the Bonferroni-corrected level of significance for multiple testing (0.05/9 = 0.006)].
Saliva samples were collected using Salivette® tubes (Sarstedt, Inc., Newton, NC, USA) for diurnal cortisol assessment and with Saliva Bio Oral Swabs (SOS) (Salimetrics, LLC, State College, PA, USA) for TSST-C cortisol reactivity. The subjects were asked to chew a swab for 1 min and then transfer it directly from their mouth to the tube. They were instructed to store their saliva samples for diurnal cortisol assessment in a freezer until they could be delivered to the research center, where samples were stored at −20 °C. The saliva samples collected during the TSST-C were directly stored at the research center. Details of salivary cortisol determination procedures are explained in the online Supplementary material.
Anxiety trait and anxiety perception during acute stress
The subscale trait of the State-Trait Anxiety Inventory (STAI) was used to evaluate general proneness to anxious behavior (STAI-Trait for children, for subjects 15 years old and under; STAI-Trait, for adolescents 16–17 years old) (Spielberger, Reference Spielberger1973). During the TSST-C, the perceived emotional arousal was assessed via the STAI-State for children scale (for children 15 and under) and the STAI-State subscale (for adolescents 16–17 years old) (Spielberger, Reference Spielberger1973). Participants answered the STAI-State questionnaire: 30 min before the stressor (T1), immediately before the stressor (T2), immediately after the stressor (T3), and 30 min after the stressor (T5) (see Fig. 1).
Statistical analysis
All statistical analyses were performed using SPSS 26 for Windows (IBM, Chicago, Illinois, USA). Descriptive statistics were analyzed by Student's t test for continuous variables and a χ2 test for categorical variables. Cortisol data were log-transformed to reduce skewness. The presence of any type of suspected or confirmed history of CM was included in downstream analysis as a dichotomic variable. The effects of both (i) CM severity (classified as: absent, low, or moderate/severe) and (ii) the frequency of CM (classified as: never, once/sometimes, or often/frequently) were also tested through independent analyses. Sensitivity analysis was conducted to explore the effects of CM when considering only subjects with a confirmed history of CM (with clear evidence from social service reports or family), aggregating those with suspected history of CM together with those without CM (see online Supplementary material).
To examine the effect of CM in diurnal cortisol slopes and changes in cortisol and anxiety perception across the TSST-C, mixed-effects models with a random effect of intercept and a random slope of time, to account for within-subject correlations, were used. Interaction with time was considered the main effect of interest of the model. Time factor had four categories (time-points) for diurnal cortisol and anxiety perception during TSST-C, and five categories for cortisol during TSST-C. In addition, simple effects tests were performed to evaluate the specific time point interaction between groups. Additionally, the overall cortisol secretion during a normal day and throughout the experimental protocol was summarized applying: (i) the area under the curve with respect to ground (AUCg) to explore the total hormonal output, and (ii) the area under the curve with respect to increase (AUCi) to reflect hormonal changes over time (Pruessner, Kirschbaum, Meinlschmidt, & Hellhammer, Reference Pruessner, Kirschbaum, Meinlschmidt and Hellhammer2003). Differences in AUCg, AUCi, and STAI-Trait scores between CM groups were tested by ANOVA. All the analyses were adjusted for the following covariates, as previously described to influence cortisol output during the TSST (Allen, Kennedy, Cryan, Dinan, & Clarke, Reference Allen, Kennedy, Cryan, Dinan and Clarke2014; De Punder, Heim, & Entringer, Reference De Punder, Heim and Entringer2019; Lê-scherban et al., Reference Lê-scherban, Brenner, Hicken, Needhamm, Seeman, Sloan and Roux2018; Marceau & Abel, Reference Marceau and Abel2018): clinical status, sex, pubertal stage, psychopharmacological treatment, illegal drugs use, oral contraceptive use, corticosteroid medication, ethnicity, SES, and BMI. In the diurnal cortisol analyses, the time of first cortisol sample (B1) collection was also included as a covariate. Specifically, in the ANOVA analysis, in order to study the direct effect of clinical status, sex, and pubertal stage on cortisol and anxiety, as well as their potential interactions with CM, these variables were included as inter-subject factors. To correct for the testing of three different CM variables (presence/absence of CM, CM severity, and CM frequency) and three different cortisol summary measures (mixed model, AUCg, and AUCi), in Table 2, a Bonferroni correction was applied by dividing the original α level (p < 0.05) by 9 (3 × 3), and obtained a new significance level of p < 0.006. Spearman's non-parametric correlation was calculated separately in participants without CM and those with a history of CM, to explore the relationship between anxiety perception and salivary cortisol during basal conditions and during the TSST-C.
Results
Attrition and descriptive analysis
Nine subjects had no information available on diurnal cortisol levels, so they were not included in the diurnal cortisol analysis. Three participants had no information available on cortisol and anxiety perception during the TSST-C, so they were not included in the corresponding analysis. Sixteen subjects were excluded from the analysis due to missing information on covariates such as BMI or SES. All the excluded participants due to missing BMI or SES values were diagnosed with a current psychiatric disorder. There were no significant differences in either sociodemographic factors or cortisol values when comparing the participants excluded and subjects with psychiatric diagnostic included in the analysis; however, the excluded participants exhibited significantly higher CGAS than those included (t = 2.360, p = 0.020).
A brief summary of the sociodemographic and anthropometric variables, by CM history, is provided in Table 1. Significant group differences according to CM exposure were observed with regard to age, ethnicity, SES, illegal drug use, CGAS, GFES, BMI, current psychiatric disorder, and type of clinical care unit. Mean cortisol values by CM group measures at each diurnal and TSST-C timepoint, AUCg and AUCi values, and STAI-Trait and STAI-State scores are summarized in Table 1.
Childhood maltreatment and diurnal salivary cortisol
As expected, cortisol levels fluctuated significantly throughout the day, following a circadian rhythm (F = 218.307, p < 0.001). No global interaction between time and CM was detected (F = 1.467, p = 0.225), reflecting a similar cortisol diurnal trajectory in both groups (see Table 2), also evidenced by AUCi levels, F (1,160) = 3.040, p = 0.083, ηp 2 = 0.021. However, the simple effects analysis in the context of mixed-effect model revealed a significant time point-specific interaction at B4 (before bedtime) between CM groups (F = 4.678, p = 0.032). Although cortisol levels consistently decreased from lunchtime to bedtime in both groups, this was less pronounced in the CM group, leading to a higher total hormonal output over the whole day, as evidenced by a higher AUCg, F (1,160) = 12.244, p = 0.001, ηp 2 = 0.079 (see Table 2 and Fig. 2). No significant interactions have been reported between CM and clinical status, pubertal stage, or sex. The effect of clinical status, pubertal stage, and sex on diurnal cortisol levels is reported in the online Supplementary material. Similar results were observed in the diurnal cortisol response when considering only subjects with confirmed CM (see online Supplementary material).
Neither the frequency nor the severity of CM was associated with diurnal cortisol slope during the day, F severity = 1.214, p = 0.300; F frequency = 1.085, p = 0.372, reflecting a similar cortisol diurnal trajectory between groups, also evidenced by AUCi, F severity(2,160) = 1.716, p = 0.184, ηp 2 = 0.024; F frequency(2,160) = 1.276, p = 0.282, ηp 2 = 0.018. However, the simple effect analysis revealed a significant interaction at B4 (before bedtime); participants exposed to moderate/severe CM experiences or often/frequently exposed to CM showed higher cortisol levels before bedtime when compared with subjects without CM (p severity = 0.020; p frequency = 0.048). The AUCg levels suggested a dose–response relationship between CM severity/frequency and total cortisol output during the day, F severity(2,160) = 6.349, p = 0.002, ηp 2 = 0.084; F frequency(2,160) = 6.068, p = 0.003, ηp 2 = 0.081. As expected, these results were even more significant when dichotomizing the sample according to the severity/frequency of CM as either: (1) no/low exposure or (2) moderate/severe exposure (see online Supplementary material).
Childhood maltreatment and salivary cortisol response during acute psychosocial stress (TSST-C)
Cortisol levels during the TSST-C significantly differed as a function of time (F = 8.953, p < 0.001), indicating the validity of this procedure to stimulate cortisol secretion in our cohort. A significant interaction between CM and time was identified (F = 4.530, p = 0.002), indicating a different trajectory of cortisol levels during the protocol between groups of CM. Specifically, the simple effects analysis in the context of mixed-effect model revealed a significant time point-specific interaction when comparing cortisol levels at T3 (immediately after the stressful situation) (F = 4.993; p = 0.027), at T4 (15 min after the stressful situation finished) (F = 10.404, p = 0.001), and at T5 (30 min after the stressful situation finished) (F = 4.561, p = 0.034). While in individuals without CM the cortisol levels increased after acute stress, there were no changes in cortisol concentration in subjects with CM (see Fig. 3a and Table 2). In line with this, participants with CM showed lower levels of AUCi than those without CM, F (1,165) = 4.779, p = 0.030, ηp 2 = 0.031, reflecting fewer hormonal changes over time. In contrast, CM was not associated with a global difference in cortisol levels throughout the entire TSST-C procedure (F = 3.015, p = 0.084), as also indicated by AUCg, F (1,165) = 0.091, p = 0.763, ηp 2 = 0.001. Similar results were observed in cortisol response during TSST-C when considering only subjects with a confirmed history of CM (see online Supplementary material). Sex, pubertal stage, and clinical status did not interact with CM, and none of these variables explained a different response pattern during the TSST-C. However, significant differences were observed in the overall cortisol levels according to pubertal stage and clinical status. Adolescents showed higher levels of cortisol (AUCg) when compared with children, and subjects with a current psychiatric diagnosis reported lower levels of cortisol (AUCg) when compared with healthy participants (further details in the online Supplementary material).
When the severity and frequency of CM were analyzed, significant interactions were again identified between CM and time (F severity = 2.773, p = 0.006; F frequency = 2.665, p = 0.008). Specifically, the simple effects analysis revealed a significant time point-specific interaction when comparing cortisol levels at T3 (immediately after the stressful situation) (p severity = 0.012; p frequency = 0.026), at T4 (15 min after the stressful situation finished) (p severity = 0.001; p frequency = 0.001), and at T5 (30 min after the stressful situation finished) (p severity = 0.033; p frequency = 0.023). While subjects without CM showed an increase in cortisol levels after the stressor, those exposed to moderate/severe or often/frequent CM were characterized by a blunted response, suggesting a dose–response relationship between CM severity/frequency and cortisol fluctuation during the TSST-C (see Fig. 3c and e). In this vein, participants exposed to moderate/severe and often/frequent CM displayed significantly lower values of AUCi than those without CM or exposed to low severity/frequency of CM, F severity(2,165) = 3.921, p = 0.022, ηp 2 = 0.052; F frequency(2,165) = 3.194, p = 0.044, ηp 2 = 0.042 (see Table 2). As expected, these results were even more significant when a new dichotomization was performed for severity/frequency of CM as either: (1) none or low and (2) moderate or severe exposure (for details see online Supplementary material). No significant differences in overall cortisol levels during the protocol were observed between severity/frequency groups of CM (F severity = 1.736, p = 0.179; F frequency = 1.839, p = 0.162), also evidenced by AUCg, F severity(2, 165) = 0.057, p = 0.945, ηp 2 = 0.001; F frequency(2, 165) = 0.074, p = 0.929, ηp 2 = 0.001.
Childhood maltreatment, anxiety trait, and anxiety perception during acute psychosocial stress (TSST-C)
Participants with CM exhibited significantly higher levels of anxiety trait than those without CM, F (1,160) = 9.129, p = 0.003, ηp 2 = 0.060. The severity and frequency of CM were also associated with anxiety trait, F severity(2,160) = 5.109, p = 0.007, ηp 2 = 0.062; F frequency(2,160) = 4.102, p = 0.019, ηp 2 = 0.056, with the lowest anxiety trait levels exhibited by subjects none exposed to CM (see Table 2). No significant correlation between anxiety trait and overall diurnal cortisol levels was found (see online Supplementary material).
As seen in Fig. 3b, the TSST-C consistently increased perceived anxiety after acute stress in all the subjects (F = 34.544, p < 0.001). However, there were no interactions between time and CM (F = 1.742, p = 0.160), reflecting similar trajectories of perceived anxiety during the acute psychosocial stress in both subjects with and those without CM. Furthermore, those with CM showed higher overall perceived anxiety during the entire procedure than subjects without CM (F = 23.836, p < 0.001). Moreover, in youth without CM, anxiety perception during the TSST-C was negatively correlated with cortisol levels, but not in youth exposed to CM (see online Supplementary material). Subjects exposed to both low and high severity/frequency of CM showed higher overall levels of anxiety during the whole protocol than subjects without CM, F severity = 11.112, p < 0.001; F frequency = 12.142, p < 0.001 (see Fig. 3c and d). However, there were no differences between groups in the magnitude of the increase of perceived anxiety after the acute stressor, F severity = 1.670, p = 0.131; F frequency = 1.240, p = 0.287, ηp 2 = 0.022, with all groups exhibiting the same trajectory (see Table 2). Similar results were obtained when considering only subjects with a confirmed history of CM (for details see online Supplementary material).
Discussion
The present study elucidated how the proximal CM in children and adolescents impacts on HPA-axis functioning and on anxiety perception. In summary, youth exposed to CM, regardless of the presence of a current psychopathology, showed (i) a basal disruption of the HPA-axis circadian rhythm with increased daily cortisol levels, (ii) reduced HPA-axis reactivity during an acute psychosocial stress, and (iii) increased anxiety perception as a trait and during the whole psychosocial stress episode. Interestingly, all the subjects exposed to CM experienced heightened anxiety but only those exposed to more severe or frequent CM exhibited significant HPA-axis dysregulation. To the best of our knowledge, this is the first study to date to report the impact of CM severity measured as the gravity of the experiences suffered, rather than as the accumulation of different types of CM (e.g. pinch with momentary redness considered as low physical abuse, v. physical aggression that needs medical intervention considered as very severe).
Our results suggest that subjects who have suffered CM have higher overall diurnal cortisol levels. Specifically, the participants with CM were characterized by a blunted decline of cortisol levels from lunchtime to bedtime, compared with those without CM. This alteration of the circadian cortisol rhythm is consistent with the presence of hypercortisolism, as evidenced by higher AUCg scores in the group exposed to CM, especially those exposed to more severe and frequent CM. This may indicate a desynchrony trend in this intrinsic biological process, which has been described as a risk factor for rising mental health symptoms. Our findings are accordant with other studies focused on CM, which have reported both a blunted decline in HPA-axis activity throughout the day (Bernard, Zwerling, & Dozier, Reference Bernard, Zwerling and Dozier2015) and higher overall cortisol output (Cicchetti & Rogosch, Reference Cicchetti and Rogosch2001). Our results could help to elucidate the co-occurrence of hypercortisolism and a flattened diurnal cortisol response, as high diurnal cortisol levels may be explained by an atypical diurnal decline. Similar findings have been reported in adults exposed to childhood adversities, suggesting the persistence of a less pronounced diurnal cortisol slope (Kuras et al., Reference Kuras, Assaf, Thoma, Gianferante, Hanlin, Chen and Rohleder2017). This HPA-axis dysregulation has important implications for other biological functions, as immune system (e.g. compromising the release of pro and anti-inflammatory substances) ultimately contributes to the increased risk of chronic disease later in life.
Although a recent meta-analysis (Bernard et al., Reference Bernard, Frost, Bennett and Lindhiem2017) reported no overall effect of CM on the diurnal cortisol slope, the authors also discussed the impact of many confounders. For example, age may influence the association between CM and cortisol rhythms; whereas cortisol levels could be elevated soon after the onset of a stressor (hypercortisolism), they could decrease over time, reflecting a pattern of hypocortisolism in adulthood (Miller, Chen, & Zhou, Reference Miller, Chen and Zhou2007). Although we did not observe this interaction between CM and pubertal stage, diurnal cortisol levels showed to be higher in adolescents when compared to children. Furthermore, our findings suggest that CM is associated with biological alterations also in youth without psychiatric disorders. In this regard, different approaches suggest that resilient subjects, who were exposed to CM but are asymptomatic, may present a particular neurobiological adaptive response, as brain connectivity changes to compensate for the alterations caused by abuse (Ohashi et al., Reference Ohashi, Anderson, Bolger, Khan, McGreenery and Teicher2019).
Secondly, regarding the HPA-axis response to acute psychosocial stress, consistent with the extant literature (Bunea et al., Reference Bunea, Szentágotai-t and Miu2017), children and adolescents exposed to CM exhibited a blunted cortisol response during the TSST-C, compared with those without CM. While previous literature supports that the blunted cortisol response is better observed in adult populations (while arguing that smaller effects are seen in children and adolescents due to HPA-axis hyperactivation following immediate trauma), an early hypoactivation is already observed in our sample, as has been reported previously (MacMillan et al., Reference MacMillan, Georgiades, Duku, Shea, Steiner, Niec and Schmidt2009). Although subjects exposed to CM remained hyporeactive under acute stress, in terms of HPA-axis activity, they experienced a significant increase in perceived anxiety, equivalent to that experienced by those not exposed to CM. This reveals a clear dissociation between anxiety perception and the physiological response to stressful situations in young people with CM, which might impair their ability to manage appropriately and cope with everyday emotionally negative situations (Liu et al., Reference Liu, Chaplin, Wang, Sinha, Mayes and Blumberg2012). Notably, emotion regulation deficits have been suggested as a key pathway linking CM with psychopathology (Dvir, Ford, Hill, & Frazier, Reference Dvir, Ford, Hill and Frazier2014; Hart et al., Reference Hart, Lim, Mehta, Simmons, Mirza and Rubia2018). Further studies are required to explore which biomarkers other than cortisol might be linked with heightened anxiety in subjects exposed to CM (Quidé et al., Reference Quidé, Bortolasci, Spolding, Kidnapillai, Watkeys, Cohen-Woods and Green2019). Our results further suggest that, although participants with a current psychopathology tended to have lower cortisol levels in general, the HPA-axis alterations in subjects exposed to CM were present in both subjects with and without a current psychopathology. Contrary to some previous findings, in our sample neither pubertal stage (Gunnar, Wewerka, Frenn, Long, & Griggs, Reference Gunnar, Wewerka, Frenn, Long and Griggs2009) nor sex (Trickett et al., Reference Trickett, Gordis, Peckins and Susman2014) interacted with CM to predict HPA-axis reactivity during the TSST-C.
Furthermore, besides the impact of CM on HPA-axis activity and increased levels of anxiety (trait and state), our findings also showed that the severity and frequency of CM play a key role, thereby supporting a dose–response relationship (Anda et al., Reference Anda, Felitti, Bremner, Walker, Whitfield, Perry and Giles2006). Thus, in line with Trickett et al. (Reference Trickett, Gordis, Peckins and Susman2014), subjects exposed to more severe or more frequent forms of CM manifested the most subdued HPA-axis responses under basal conditions and in response to psychosocial stress; notably, Trickett et al. considered severity as the accumulation of different types of CM, rather than according to the specific characteristics of the experiences suffered. These findings warn of the deleterious impact that milder forms of CM may have once they become chronic. This is important as children who experience mild CM are often not detected or receive less clinical and social care (Humphreys, Reference Humphreys2020). Furthermore, it seems that these children start showing higher levels of perceived anxiety before there is a marked biological dysregulation, offering a window of opportunity for early detection and intervention. Hence, the use of accurate child screening instruments at subclinical stages should be generalized, since most children are only identified once they already have severe psychiatric symptomatology (Bailhache, Leroy, Pillet, & Salmi, Reference Bailhache, Leroy, Pillet and Salmi2013). Moreover, since dysfunction in neurobiological systems negatively impacts treatment outcomes, youth with CM may also require specific treatment adapted to their condition (Tyrka, Burgers, Philip, Price, & Carpenter, Reference Tyrka, Burgers, Philip, Price and Carpenter2013).
The methodology used in the present study includes a wide range of CM experiences reported from different sources, since there is often a substantial gap between subjects identified in informant-based studies and self-report assessments (Baldwin, Reuben, Newbury, & Danese, Reference Baldwin, Reuben, Newbury and Danese2019). Thus, our findings suggest that participants with a suspected history of CM identified by clinicians show the same HPA-axis dysfunctions as subjects with a confirmed history of CM. Likewise, given that CM studies may lack sensitivity when the experiences are not qualitatively assessed (via the severity and frequency of exposure), key information may be lost and findings distorted. This highlights the need for specific training of clinicians in child psychiatric and pediatric services, so CM assessment can be routinely implemented, despite the time and effort required to perform such complex assessments (Zeanah & Edm, Reference Zeanah and Edm2018).
Although prior evidence suggests that exposure to CM during middle childhood has the greatest effects on emotional dysregulation (Dunn, Nishimi, Gomez, Powers, & Bradley, Reference Dunn, Nishimi, Gomez, Powers and Bradley2018), it is difficult to pinpoint the exact developmental period when HPA-axis functioning is disrupted. Future research should incorporate more detailed information about the timing and proximity of CM to delineate vulnerable periods (Andersen & Teicher, Reference Andersen and Teicher2008). It would be interesting to study the clinical course of the children to identify possible risk and protective factors for the future onset of psychopathology. A more dimensional approach focused on symptom dimensions might reveal varying patterns of adrenocortical regulation (Cicchetti & Rogosch, Reference Cicchetti and Rogosch2001). It is important to note that CM is not a phenomenon that can be studied in isolation, since both its causes and consequences are systematic and there are many factors that must be taken into account in order to fully understand it.
The blunted reactivity observed in our study supports plausible habituation, i.e. chronic exposure to stress may be linked with an adaptive desensitization to new stressors over time (Murali & Chen, Reference Murali and Chen2005). These latent neurobiological alterations could drive an increased vulnerability to psychopathology during childhood and adolescence (Busso et al., Reference Busso, Mclaughlin, Brueck, Peverill, Gold and Sheridan2017), which may persist, leading to the onset of a wide range of psychiatric conditions in adulthood (Kudielka & Wüst, Reference Kudielka and Wüst2010). Other factors with the potential to moderate the consequences of CM should also be taken into account, such as the type of CM suffered, the relationship with the abuser, social support received, and coexistence of other types of trauma such as bullying (Arseneault, Reference Arseneault2018), domestic violence (Osofsky, Reference Osofsky2018), or recent stressful life events (March-Llanes, Marqués-Feixa, Mezquita, Fañanás, & Moya-Higueras, Reference March-Llanes, Marqués-Feixa, Mezquita, Fañanás and Moya-Higueras2017).
One of the limitations of the current study is the methodology used for assessing the presence and characteristics of CM exposure. Widely used questionnaires such as the CTQ cannot be administered to children younger than 12 years; indeed, there is no validated questionnaire to assess the presence of CM in the 7–17 years range. The main reason behind this is that younger children have a limited understanding of their own exposure, since they are still cognitively immature. Thus, any assessment of CM in this vulnerable population needs to be adjusted to maximize the reliable information that can be captured from the different informants (not only the child) and, at the same time, to minimize the trauma that the interview itself can represent to a victimized child. Thus, use of TASSCV allows the proper assessment of children and adolescents exposed from milder to severe forms of CM, which would have otherwise not been identified. Unfortunately, use of TASSCV requires a longer time for a proper assessment together with the gathering of information from multiple informants, which might make it more challenging to use than simply relying on short self-administered questionnaires such as the CTQ or considering only the most severe children already detected by social services. Since most of the sample was recruited in psychiatric units, there is an unusually high proportion of ADHD cases in the non-CM group; thus, our findings might not be generalizable to other populations. At the same time, the majority of CM-exposed subjects suffered from some sort of psychiatric condition, while most participants non-exposed to CM had no psychopathological history. Further research including a higher proportion of subjects exposed to CM with no psychiatric symptomatology (i.e. resilient) is required to disentangle the role of CM in the development of HPA-axis disturbances and whether the later precede the onset of psychiatric disorders.
Conclusions
CM affects multiple domains of life such as intimate relationships, violence and criminal offending, employment, drug abuse, and physical and mental health (Hughes et al., Reference Hughes, Bellis, Hardcastle, Sethi, Butchart, Mikton and Dunne2017). It is a serious global health problem with staggering long-term economic costs (Thielen et al., Reference Thielen, ten Have, de Graaf, Cuijpers, Beekman, Evers and Smit2016). This study is intended to raise awareness of the biological and clinical repercussions of CM during or proximally to exposure, encouraging clinicians to ask patients about CM history and to respond accordingly, seeking therapeutic alternatives to manage acute stress better. Children exposed to CM and attended in child protection units, child psychiatric, or pediatric units are still at a sensitive period of neurological, cognitive, social, and emotional development, during which high-quality interventions can make an important difference and shift the balance between risk and protective factors (Chinitz, Guzman, Amstutz, Kohchi, & Alkon, Reference Chinitz, Guzman, Amstutz, Kohchi and Alkon2017). Thus, family psychotherapeutic interventions have the potential to normalize HPA-axis function if implemented promptly (Gunnar, DePasquale, Reid, Donzella, & Miller, Reference Gunnar, DePasquale, Reid, Donzella and Miller2019).
Supplementary material
The supplementary material for this article can be found at https://doi.org/10.1017/S003329172100249X
Acknowledgements
We are indebted to all the participants and their families for taking part in a study with such deep emotional involvement. We also thank all the nurses and lab technicians who made this research possible, especially Anna Valldeperas and José Luís Monteserín.Thank you Sergi Valero for your statistical support.
Financial support
Supported by the Spanish Ministry of Economy and Competitiveness, Instituto de Salud Carlos III through the University of Barcelona multicenter project (PI15/00097)-PI L. Fañanas, Hospital Universitario Araba (PI15/00793)-PI I. Zorrilla, Hospital Gregorio Marañón (PI15/00723)-PI M. Rapado and Hospital Clinic (PI15/00685)-PI S. Romero, co-financed by the European Regional Development Fund project, ‘A way of making Europe’. We thank the Network Centre for Biomedical Research in Mental Health (CIBERSAM). This study was facilitated by a pre-doctoral research grant from the Catalonian authorities to Laia Marques-Feixa (AGAUR-FI_B100023-2018). Supported by a NARSAD Distinguished Investigator Grant awarded to Professor Lourdes Fañanas (26887) and by the Comissionat per a Universitats i Recerca del DIUE, of the Generalitat de Catalunya regional authorities (2017SGR1577). Dr Rapado-Castro was supported by a Ramon y Cajal Research Fellowship (RYC-2017-23144), Spanish Ministry of Science, Innovation and Universities; an NARSAD independent investigator grant (24628) from the Brain & Behaviour Research Foundation; Fundación Familia Alonso and Fundación Alicia Koplowitz.
Conflict of interest
In the last 24 months, Hilario Blasco-Fontecilla received lecture fees from AB-Biotics, Praxis Pharmaceuticals, and Shire. The remaining authors do not have any conflict of interest regarding the publication of this manuscript.
Ethical standards
The authors assert that all procedures contributing to this work comply with the ethical standards of the relevant national and institutional committees on human experimentation and with the Helsinki Declaration of 1975, as revised in 2008.
EPI-Young Stress GROUP
María Martín9, Eulalia Anglada12, Pilar Santamarina2,3,4, Patricia Rubio6, Iría Méndez2,3,4, Dolores Moreno2,6, María Ramos2,6, Jorge Vidal2,6, Juanjo Carballo2,6, Elena Font2,3,4, Lydia Gayubo11, Laura Colino11, María Rodrigo-Yanguas11, Maddi Laborde8, Jaume March-Llanes5.