The public health impact of depression may be substantially mitigated if adequate attention is directed to effectively understand and treat depression early in development. Depressive disorders are associated with impairment, chronic suffering, and early death, and impact about 16% of the population (Kessler, Avenevoli, & Merikangas, Reference Kessler, Avenevoli and Merikangus2001). Historical trends suggest that depression is on the rise and is the third leading cause of global burden of disease worldwide (Berndt et al., Reference Berndt, Koran, Finkelstein, Gelenberg, Kornstein and Miller2000; World Health Organization, 2008). Depression in adolescence is of particular importance (Zalsman, Brent, & Weersing, Reference Zalsman, Brent and Weersing2006). Not only is depression commonly first evident during adolescence, but an early onset of depression is associated with a poor prognosis (Lewinsohn, Clarke, Seeley, & Rohde, Reference Lewinsohn, Clarke, Seeley and Rohde1994; Weissman et al., Reference Weissman, Wolk, Goldstein, Moreau, Adams and Greenwald1999; Zisook et al., Reference Zisook, Lesser, Stewart, Wisniewski, Balasubramani and Fava2007).
Major depressive disorder (MDD) has been characterized as a multisystemic disorder affecting brain and body (Insel & Charney, Reference Insel and Charney2003). Inclusion of multiple levels of analysis provides an opportunity to examine the interplay across relevant systems. The focus on depression early in development is a priority because adolescents are more sensitive to stress (Compas & Wagner, Reference Compas, Wagner, Gore and Colton1991), and the neurobiological systems relevant to threat detection and stress regulation are continuing to undergo maturational refinement (e.g., Lenroot & Giedd, Reference Lenroot and Giedd2006; Luciana & Collins, Reference Luciana and Collins2012; Romeo & McEwen, Reference Romeo and McEwen2006). Neuroscience research on adolescent MDD to date has identified anomalous functioning in systems involved in responding to threats in the environment, including key brain regions (e.g., Cullen et al., Reference Cullen, Gee, Klimes-Dougan, Gabbay, Hulvershorn and Mueller2009, Reference Cullen, Klimes-Dougan, Muetzel, Mueller, Camchong and Houri2010; Thomas et al., Reference Thomas, Drevets, Dahl, Ryan, Birmaher and Eccard2001; Yang et al., Reference Yang, Simmons, Matthews, Tapert, Frank and Max2010) and the hypothalamic–pituitary–adrenal (HPA) axis (e.g., Klimes-Dougan, Hastings, Granger, Usher, & Zahn-Waxler, Reference Klimes-Dougan, Hastings, Granger, Usher and Zahn-Waxler2001; Rao, Hammen, Ortiz, Chen, & Poland, Reference Rao, Hammen, Ortiz, Chen and Poland2008). However, these approaches are limited by focusing primarily on either the neural or the hormonal aspects of the biological threat response system. Research with adult depression has begun to examine interplay across systems; while these findings may have limited developmental relevance, consideration of multiple levels of analysis provides a useful framework for advancing our understanding of the complex neurobiology that underlies the pathophysiology of depression (e.g., Pruessner et al., Reference Pruessner, Dedovic, Pruessner, Lord, Buss and Collins2010). The current work uses multiple levels of analysis to examine the interplay of systems relevant to threat response.
There are several existing models that highlight the challenges of threat processing for those struggling with depression (Drevets, Reference Drevets1999; Ghashghaei & Barbas, Reference Ghashghaei and Barbas2002; Mayberg, Reference Mayberg1997; Nestler et al., Reference Nestler, Barrot, DiLeone, Eisch, Gold and Monteggia2002; Phillips, Drevets, Rauch, & Lane Reference Phillips, Drevets, Rauch and Lane2003; Price & Drevets, Reference Price and Drevets2010). Fronto-limbic circuitry and HPA axis functioning are two key systems important for threat processing, and preclinical studies have demonstrated clear links between these systems (e.g., Diorio, Viau, & Meaney, Reference Diorio, Viau and Meaney1993; McEwen, Reference McEwen and Gazzaniga1995; Reul & de Kloet, Reference Reul and de Kloet1985; Sullivan & Gratton, Reference Sullivan and Gratton2002). Presumably in certain pathological conditions, excessive limbic activation may lead to overstimulation of the HPA axis, resulting in the release of stress hormones whose cumulative effects include alterations in receptor functioning as well as deleterious long-term consequences for neuronal health (e.g., McEwen, Reference McEwen and Gazzaniga1995; Musselman & Nemeroff, Reference Musselman and Nemeroff1993). There is preliminary evidence that HPA axis normalization can be achieved when treatment is effective (Fisher, Gunnar, Chamberlain, & Reid, Reference Fisher, Gunnar, Chamberlain and Reid2000; Pariante, Kim, Makoff, & Kerwin, Reference Pariante, Kim, Makoff and Kerwin2003). Adolescence may represent a critical window of development in which interventions for depression could be most successful in terms of alternating these stress response patterns (Levine Reference Levine1957), reducing the likelihood of neurodegeneration and gene expression alterations (de Kloet, Reference De Kloet2003; Kaufman & Charney, Reference Kaufman and Charney2001; Meaney & Szyf, Reference Meaney and Szyf2005).
Highly complex neural networks make up the threat response system of the brain and the associated hormonal cascade. Given that this field is still in a preliminary phase of inquiry, we chose to focus on the amygdala, a central node of the threat response network circuitry that is posited to affect hormonal systems. There is evidence of excitatory modulation of the amygdala on the HPA axis (Van de Kar & Blair, Reference Van de Kar and Blair1999). The amygdala and the HPA system are both key components in a response system that detects and orchestrates a regulatory response from distributed networks. The amygdala also affects the glucocorticoid feedback processing that is disrupted in depression (Herman, Flak, & Jankord, Reference Herman, Flak and Jankord2008).
Several studies to date have considered the interplay between amygdala functioning and HPA axis functioning in healthy adults and adults suffering from depression. Drevets et al. (Reference Drevets, Price, Bardgett, Reich, Todd and Raichle2002) found a positive association between glucose metabolism in the amygdala, as assessed by positron emission tomography scans, and cortisol levels during the scan in depressed adults. Cunningham-Bussel et al. (Reference Cunningham-Bussel, Root, Butler, Tuescher, Pan and Epstein2009) reported that the right amygdala blood oxygen level dependent response to visual images of the World Trade Center attack was positively correlated with cortisol values during a scan in healthy adults, and right amygdala activation was correlated with the cortisol pre- versus postscan change scores. Another study showed that exogenous cortisol administration led to increased noradrenergic activation in the amygdala (van Stegeren et al., Reference van Stegeren, Wolf, Everaerd, Scheltens, Barkhof and Rombouts2007). Exogenous cortisol has also been found to “decouple” the amygdala from executive control brain regions (Henckens, van Wingen, Joels, & Fernandez, Reference Henckens, van Wingen, Joels and Fernadez2012). While there is some conflicting evidence (e.g., Holsen et al., Reference Holsen, Lancaster, Klibanski, Whitfield-Babrieli, Cherkerzian and Ubka2013; Lovallo, Robinson, Glahn, & Fox, Reference Lovallo, Robinson, Glahn and Fox2010), overall these studies suggest that elevated amygdala activation is associated with cortisol activation in both healthy and MDD populations.
In addition to amygdala functioning, amygdala volume may also be linked with HPA axis functioning. Some work has suggested links between the volumes of other limbic structures, including the pituitary gland and the hippocampus, to HPA axis functioning in depressed adults (Axelson et al., Reference Axelson, Doraiswamy, Boyko, Rodrigo Escalona, McDonald and Ritchie1992, Reference Axelson, Doraiswamy, McDonald, Boyko, Tupler and Patterson1993; Dedovic et al., Reference Dedovic, Engert, Duchesne, Lue, Andrews and Efanov2010; Treadway et al., Reference Treadway, Grant, Ding, Hollon, Gore and Shelton2009), while other studies have failed to document such an association (e.g., Colla et al., Reference Colla, Kronenberg, Deuschle, Meichel, Hagen and Bohrer2007; Vythilingam et al., Reference Vythilingam, Vermetten, Anderson, Luckenbaugh, Anderson and Snow2004). Among the few studies that have examined amygdala volume more specifically, there has been limited evidence of across-system associations (Kronenberg et al., Reference Kronenberg, Tebartz van Elst, Regaen, Deuschle, Heuser and Colla2009; Pruessner et al., Reference Pruessner, Dedovic, Pruessner, Lord, Buss and Collins2010; Schuhmacher et al., Reference Schuhmacher, Mössner, Jessen, Scheef, Block and Belloche2012). These results may be due to variance in methods used to assess HPA axis functioning, which range from considering cortisol under basal conditions to highly variable stress reactive paradigms (Pressner et al., 2010). These results may also be because depression is sometimes associated with larger amygdala volumes and in other cases smaller amygdala volumes (Hamilton, Siemer, & Gotlib, Reference Hamilton, Siemer and Gotlib2008; Leuner & Shors, Reference Leuner and Shors2013). In a meta-analysis, Hamilton et al. (Reference Hamilton, Siemer and Gotlib2008) noted that depressed adults undergoing antidepressant pharmacotherapy tend to have larger amygdala volumes. Finally, because the adolescent neurobiological stress system is undergoing development, applications of the adult literature may be limited when considering adolescent depression.
Research is needed to begin parsing out the contributions of these important factors in adolescence. Performing analyses geared toward elucidating structural and hormonal interplay may be able to demonstrate new links between stress and later vulnerability for depressive and other disorders. To date, there have been only a handful of studies that have examined the interplay between systems in typically developing adolescents or in adolescents at risk for depression. Thomason, Hamilton, and Gotlib (Reference Thomason, Hamilton and Gotlib2011) found that in a sample of healthy adolescents who completed functional magnetic resonance imaging (fMRI) after a social stress paradigm, higher cortisol stress response was correlated with increased functional connectivity between the salience network and the subgenual anterior cingulate cortex. When examining potential correspondence between early life stress and HPA axis function, Burghy et al. (Reference Burghy, Stodola, Ruttle, Molloy, Armstrong and Oler2012) found that increased cortisol due to higher early life stress in young girls predicted lower functional connectivity between the amygdala and the ventromedial prefrontal cortex 14 years later. Another study, by Liu et al. (Reference Liu, Chaplin, Wang, Sinha, Mayes and Blumberg2012), showed that in a group of adolescents experiencing stressful life conditions, who completed the Trier Social Stress Test (TSST) followed by an fMRI, a greater cortisol response correlated with less activity in the left hippocampus while viewing fearful faces. Critically important questions considering links between brain and HPA functioning in clinically depressed adolescents have yet to be addressed.
The purpose of this study was to conduct multilevel assessments of threat systems in adolescents with and without depression using multiple levels of analysis, and to evaluate the correspondence across systems (neuronal structure and responses, hormonal responses, and behavior). Assessments addressed amygdala volume and functioning. Evaluation of the hormonal stress system included two paradigms to measure stress reactivity and recovery. The TSST was used to evaluate HPA axis functioning in response to a social stressor on a different day than the brain scan. To obtain an HPA measure that was more temporally linked to our neurocircuitry measure, we also assessed cortisol levels before and after the brain scan. The first study aim was to evaluate differences between depressed and well adolescents for each of these stress indices. Based on previous work (e.g., Rosso et al., Reference Rosso, Cintron, Steingard, Renshaw, Young and Yurgelun-Todd2005; Yang et al., Reference Yang, Simmons, Matthews, Tapert, Frank and Max2010), we predicted that adolescents with MDD would exhibit smaller amygdala size, greater amygdala functioning, and a different pattern of hormonal response to stress. The second study aim was to examine the interplay between amygdala and HPA axis indices. We predicted that associations between the neuronal and hormonal systems would be found and to some extent these associations would differ between MDD and healthy adolescents; however, no specific predictions were made with regard to the directions of the associations. We also explored possible differences between medicated and unmedicated depressed adolescents (Aihara et al., Reference Aihara, Ida, Yuuki, Oshima, Kumano and Takahashi2007).
Methods
Participants
Participants were 79 adolescents between the ages of 12 and 19 years old (M age = 15.90, SD = 1.87): 52 adolescents with MDD and 27 healthy controls (HC). Participants were primarily females (75.95%) and most identified themselves as Caucasian (65.8%), followed by African American (8.9%), Hispanic (8.9%), Asian (3.8%), and Native American (1.3%), with the remaining participants self-identifying as “other” (20.3%). Participants were able to select more than one option for race/ethnicity. MDD and HC participants were matched at recruitment on sex, age, and race.
Participants were recruited using a variety of different strategies, primarily through community postings, and from inpatient and outpatient clinical services at the University of Minnesota and the surrounding area. On the first visit, diagnostic interviews were conducted. The TSST and the brain scan were conducted on two separate subsequent visits. The study was approved by the University of Minnesota's Institutional Review Board. All participants completed signed informed consent and/or assent (if under 18), and all participants received monetary compensation for their participation after completing each of the three visits.
Measures
Diagnosis and symptom assessment
The presence or absence of a DSM-VI-TR Axis I disorder(s) was confirmed by a semistructured diagnostic interview. Participants under 18 years of age and a legal guardian completed independent interviews using the Kiddie Schedule for Affective Disorders and Schizophrenia—Present and Lifetime Version (KSADS-PL; Kaufman et al., Reference Kaufman, Birmaher, Brent, Rao, Flynn and Moreci1997). The KSADS-PL interviews were conducted by highly trained clinical psychologists, child psychiatrists, or advance trainees enrolled in graduate clinical psychology doctoral programs under the direct supervision of a clinician. For participants 18 or 19 years of age, a parent interview was not conducted. Many MDD participants (67%) suffered from at least one comorbid mental illness, which consisted most commonly of an anxiety disorder (67%) and/or attention-deficit/hyperactivity disorder (ADHD; 13%). In addition to diagnostic status, clinicians obtained information about the duration of illness (MDD). Upon having conducted the KSADS-PL, clinicians compete a Global Assessment of Functioning score rating (M = 54.35, SD = 8.35) and the Children's Depression Rating Scale (mean T score = 77.29, SD = 6.13; Poznanski & Mokros, Reference Poznanski and Mokros1996) on participants in the MDD group. The Children's Depression Rating Scale is a semistructured interview that assesses 17 symptom areas related to depression, including those that serve as criteria in the DSM-IV. All participants completed the Beck Depression Inventory II (Beck, Steer, & Brown, Reference Beck, Steer and Brown1996) at each study visit, which documented a significant difference for the average Beck Depression Inventory II score across visits between the HC (M = 2.33, SD = 3.40) and the MDD groups (M = 25.54, SD 12.28).
Within the context of the KSAD-PL, information about medication status was obtained. MDD participants were classified based on medication status. One subgroup of 16 participants was receiving medication for depression and another subgroup of 36 participants was not receiving medication for depression. Participants diagnosed with ADHD who were only receiving stimulants for the treatment of ADHD were included in the unmedicated sample provided that they abstained from taking the medication on the day of the brain scan. Two participants did not abstain; both had deviated from the recommended protocol and taken psychostimulants earlier in the day, but by the time of the scan the effects of the medication were likely partially or totally worn off given the typical half-life of approximately 6 hr for these medications.
DSM-IV diagnoses were established through a consensus meeting that incorporated information from the independent parent and child interviews, rating scales, and, if available, medical records. Participants were eligible for the HC group if they had no evidence of an Axis I diagnosis. Participants were eligible for the patient group if they had a primary diagnosis of MDD. Exclusionary diagnoses for both the MDD and HC groups included pervasive developmental disorder, bipolar disorder, schizophrenia, a neurological disorder, or a chronic or serious medical condition. An additional exclusion criteria was evidence of a below average IQ. Participants were excluded if their IQ, as assessed by the Wechsler Abbreviated Scale of Intelligence (Wechsler, Reference Wechsler1999), was lower than 80. Although the average IQ scores for both groups fell in the normal range, as shown on Table 1, the HC group had a significantly higher IQ than the MDD group (with the average IQs of the unmedicated MDD group appearing the lowest). In the eight cases where IQ data was missing that were retained in this sample, all participants were achieving at an average or above rate at school and had not been identified for special services, suggesting that their IQ was likely to be at or above the average range of functioning.
Table 1. Demographic and clinical characteristics of adolescents with MDD and HC participants
Note: MDD, Major depressive disorder; HC, healthy control; Med, medicated; No Med, not medicated; CDRS, Children's Depression Rating Scale; BDI, Beck Depression Inventory; ADHD, attention-deficit/hyperactivity disorder; TSST, Trier Social Stress Test; MRI, magnetic resonance imaging.
*p ≤ .05 when compared to healthy controls. **p ≤ .01 when compared to healthy controls. ***p ≤ .001 when compared to healthy controls.
Assessment of brain structure and function
Participants completed a MRI scan using a Siemens 3 Tesla TIM Trio scanner (Erlangen, Germany) that is housed at the Center for Magnetic Resonance Research at the University of Minnesota. A 5-min scan was administered to acquired to obtain a structural image using a T1-weighted high-resolution magnetization prepared gradient echo sequence (repetition time = 2530 ms, echo time = 3.65 ms, inversion time = 1100 ms, flip angle = 7°, 224 coronal slices, field of view = 256 mm, voxel size 1 × 1 × 1 mm, matrix size = 256 × 256, GRAPPA = 2). In addition to several other scans (e.g., resting), an fMRI scan was administered during which participants completed the subsequently described emotion face-matching task. The fMRI scan consisted of an echo planar imaging sequence, which was used to collect 197 T2-weighted whole-brain functional volumes in the context of the task (34 3.0 mm contiguous interleaved axial slices; aligned to anterior and posterior commissures with –30° tilt, repetition time = 2000 ms, echo time = 28 ms, flip angle = 80°, field of view = 200 mm, voxel size 3.1 × 3.1 × 3.0 mm, matrix = 64 × 64).
Amygdala volume analyses
Volumetric data was processed using the FreeSurfer 5.3.0 software (http://surfer.nmr.mgh.harvard.edu/), including brain extraction and parcellation of data into a standard set of anatomically based regions of white and gray matter. FreeSurfer output was visually inspected; when any errors were identified, they were manually corrected, and the pipeline's remaining steps were repeated. Although many regions of interest were identified using this software, for the purposes of this study, we focused on the volumetric data that was produced for the right and left amygdala.
Emotion face-matching task and data analysis
The emotion face-matching task (Hariri, Tessitore, Mattay, Fera, & Weinberger, Reference Hariri, Tessitore, Mattay, Fera and Weinberger2002) used E-Prime software and was projected on a screen inside the bore of the MRI scanner that the participant could see using a mirror attached to the head coil. The task entailed both affective and control stimuli. For the affective stimuli, Ekman faces (Ekman & Friesen, Reference Ekman and Friesen1976) were used to portray anger and fear, which included six images of each gender and emotion. The control stimuli consisted of circles, horizontal ellipses, and vertical ellipses. Participants were instructed to match the stimuli presented on the top row with one of the two stimuli presented in the bottom row using a button box. Specifically, participants were instructed to match the shape for the control stimuli to faces with emotional expression (fear or anger) for the affective stimuli. The task was presented in 13 24-s counterbalanced blocks (3 fixation, 5 shape, and 5 emotion). This task took approximately 6.5 min to complete in the scanner.
Analysis of fMRI data was conducted using software tools from the FMRIB software library (http://fsl.fmrib.ox.ac.uk) version 4.1.8. Preprocessing steps included motion correction, brain extraction, high-pass temporal filtering, prewhitening, regression of motion parameters, and registration to MNI standard space. A first-level analysis for each data set was conducted to regress the task model onto the fMRI data at each voxel of the brain. We included a covariate of no interest, which allowed us to regress out those volumes in which motion (relative to the preceding volume) exceeded our threshold of 1.5 mm in any direction (half the size of one voxel). We considered two explanatory variables from the block-design task (matching emotion faces and matching neutral shapes) and two contrasts (matching emotion faces minus fixation and matching emotion faces minus matching shapes). Left and right amygdala masks from the Harvard Oxford Subcortical Structural Atlas were then used to extract the average z score for all voxels within these regions from the emotion minus fixation and the emotion minus shape contrasts of the regression results for each participant. These average z score values were then used for group comparison and correlation analyses. The primary index of focus here was on amygdala (right and left) activation for the emotion face minus fixation contrast. Follow-up analyses were also conducted with the index that assesses amygdala functioning for emotion minus shapes.
HPA axis assessments
Participants completed a slightly modified version of the TSST, a task that has been found to reliably elicit a stress response (Kirschbaum, Pirke, & Hellhammer, Reference Kirschbaum, Pirke and Hellhammer1993). Participants were asked to spend 5 min preparing a speech to introduce themselves to a job committee and were informed that another task would follow the speech task. After the preparation period, participants were escorted to another room in front of two unfamiliar evaluators wearing white lab coats, who were trained to remain neutral and to avoid giving reassurance or feedback. Participants were first asked to deliver their speech (5 min) and then asked to do a serial subtraction task with corrective feedback provided by experimenters (5 min). Participants were debriefed immediately following the completion of the TSST.
A total of five salivary samples were collected throughout this visit: the first before speech preparation (0 min), the second immediately following the TSST (15 min), and the final three at approximately 30, 45, and 60 min. Instead of lengthening the visit, researchers eased subject burden by having them complete tasks that don't typically activate stress responding (e.g., Gunnar, Talge, & Herrera, Reference Gunnar, Talge and Herrera2009), including rating scales and computerized neuropsychological tests, during the time the two final samples were collected.
Experiences and expressions of stress were recorded within the context of the TSST. After completion of the TSST, participants were asked to rate a series of questions on a scale of 1 (calm) to 5 (high stress), including, “How stressful was giving the speech (job interview)?” and “How stressful was the subtraction task?” A mean score across these two items was used as the summary score for experienced stress. The two examiners for the TSST independently rated participant's behavior on a scale of 1 (not stressed at all) to 6 (discontinued the procedure because the participant was so stressed) for the following items, “How stressed did the participant appear during the story telling task (job interview)?” and “How stressed did the participant appear during the arithmetic task?” These scores were averaged for the participant. The two examiner ratings were positively correlated r (74) = .41, p = .01. A mean score across examiner ratings was used as the summary score for expressed stress.
Cortisol was also assessed within the context of the MRI scan: there is accumulating evidence that the MRI scan procedure may be considered a moderate stressor, likely due to the movement restrictions, loud noises, and novel task features required, particularly for scan-naive participants (Eatough, Shirtcliff, Hansen, & Pollak, Reference Eatough, Shirtcliff, Hanson and Pollak2009; Peters, Cleare, Papadopoulos, & Fu, Reference Peters, Cleare, Papadopoulos and Fu2011; Tessner, Walker, Hochman, & Hamann, Reference Tessner, Walker, Hochman and Hamann2006). Researchers collected a saliva sample upon arrival to the MRI facility and again immediately after completing the 75-min scan.
For each saliva sample, participants facilitated salivary excretion by chewing Trident Original gum for 20–30 s before spitting out the saliva and gum. Participants then pushed their saliva through a straw and into a 1.5 ml vial. Samples were labeled and stored in a –25 °C freezer until they were shipped to Universität Trier in Trier, Germany, for analysis. Researchers used assay methods consistent with Dressendörfer, Kirschbaum, Rohde, Stahl, and Strasburger (Reference Dressendörfer, Kirschbaum, Rohde, Stahl and Strasburger1992). Summary indices of cortisol values: across the TSST, cortisol was represented by the area under the curve ground (AUCg) and the AUC from the initial pretest sample (AUCi; Pruessner, Kirschbaum, Meinlschmid, & Hellhammer, Reference Pruessner, Kirschbaum, Meinlschmid and Hellhammer2003). Difference scores (MRI CORT) were the primary summary index for cortisol levels linked to the brain scan.
Statistical analysis
Preliminary analyses were conducted using Kruskal–Wallis (continuous characteristics) and Fisher exact tests (categorical characteristics) to determine whether there were any demographic or other characteristic differences across the MDD and HC groups (two group analysis). We also considered how the behavioral and biological indices were correlated. To address the first study aim, general linear models were used to assess whether MDD and HC differed on amygdala structure, amygdala functioning, HPA axis functioning on the TSST, and HPA axis functioning during the MRI. To address the second study aim, a series of general linear models were conducted with each of AUCi, AUCg, and MRI CORT as separate dependent variables, to consider if there were associations with the amygdala (structure or function) either overall or differentially by risk status of the adolescent (a total of 12 analyses with three HPA outcomes and four possible predictors: right and left amygdala volume and right and left amygdala functioning). Differential associations with HPA axis functioning, modeled as the interaction between amygdala (structure or function) and risk status, would represent differences between MDD and HC participants in the interplay between the components of the biological stress system.
For the previously described analyses, we carefully considered appropriate statistical controls. With few exceptions, demographic characteristics were comparable across the groups; however, IQ was significantly different across groups. Although IQ was initially considered as an adjusting variable in all models, it was not retained (any changes in the results were most likely to be due to changes in degrees of freedom in the models and are noted subsequently). Because MDD and HC showed robust group differences on experiences and expressions of stress during the TSST (described below), we included this summary index as an adjusting variable in all models of HPA axis functioning on the TSST (e.g., AUCg). In addition, for all analyses that considered amygdala volume as either dependent or independent variable, intracranial volume was included as an adjusting variable. For all analyses that considered MRI CORT as either dependent or independent variable, time of day of the scan's onset, cortisol levels pre-MRI, and previous experience with a brain scan (yes/no) were included as adjusting variables.
Finally, as exploratory analyses, this same modeling process was repeated using three groups: unmedicated MDD, medicated MDD, and HC participants. Results were reported below when the three group comparisons were significant (p < .05).
Results
Participant characteristics by group are shown in Table 1. Preliminary analyses considered the correlations among behavioral, neuronal, and HPA axis variables of relevance to threat processing (Table 2). Right and left amygdala volumes were positively correlated, as were right and left amygdala response to the emotion face-matching task. There was no evidence that the structure and the function of the amygdala were correlated with each other in this sample. With regard to HPA axis functioning, cortisol levels were correlated between the TSST and the MRI CORT. In addition, right amygdala functioning within the context of the emotion face-matching task was positively correlated with AUCg cortisol functioning during the TSST.
Table 2. Correlations for primary study variables
Note: TSST, Trier Social Stress Test; TSST ratings, z score means of self-reports and experimenter observations; AUCg, area under the curve ground levels of salivary cortisol; AUCi, area under the curve initial levels of salivary cortisol; MRI CORT difference, pre magnetic resonance imaging – post magnetic resonance imaging levels of salivary cortisol.
*p = .05 (two tailed). **p = .01 (two tailed).
Our first study aim was to examine whether MDD and HC adolescents differed on indexes of behavioral and neurobiological functioning (descriptive information is provided on Table 3). Consideration of experiences and expression of stress within the context of the TSST was assessed across participant groups. The summary score of self-reported experiences of stress and experimenter ratings of expressed stress during the TSST were significantly higher for the MDD group than the HC group, F (72) = 20.47, p < .001). In addition, there was a significant difference across the three groups, F (72) = 20.47, p < .001). Post hoc comparisons showed that medicated MDD participants were significantly higher on behavioral stress than were HC (p < .001) and that the unmedicated MDD participants were significantly higher on behavioral stress than were HC (p = .004).
Table 3. Descriptive information about adolescents with MDD and HC participants stress functioning
Note: MDD, Major depressive disorder; HC, healthy control; Med, medicated; No Med, not medicated; HPA, hypothalamic–pituitary–adrenal axis; TSST, Trier Social Stress Test; AUCg, area under the curve ground levels of salivary cortisol; AUCi, area under the curve initial levels of salivary cortisol; MRI, magnetic resonance imaging; MRI CORT difference, pre-MRI – post-MRI levels of salivary cortisol.
There were no significant MDD and HC differences for mean volumes for each of left or right amygdalae, nor were there any significant three-group differences. For example, the MDD and HC group results for amygdala volume were F (1, 77) = 0.00, p = .96 for the left amygdala and F (1, 77) = 1.08, p = .30 for the right amygdala. As predicted, mean right amygdala activation, F (1, 64) = 5.48, p = .02, and left amygdala activation, F (1, 64) = 4.10, p = .05, was significantly different between MDD and HC (Figure 1). This significant difference for right amygdala functioning was still evident when IQ was entered as a control variable. There were no significant three-group differences for right or left amygdala functioning.
Figure 1. Amygdala activation during an emotion face-matching task.
We evaluated unadjusted MDD and HC differences for various indices of HPA axis functioning. There were no significant MDD and HC differences for the TSST for AUCg, F (1, 72) = 0.24, p = .63, or AUCi, F (1, 72) = 1.49, p = .23. Results were largely the same after controlling for important confounders. However, as shown on Figure 2, there were differences in patterns of linear slope the shape of the cortisol response versus time for the patterns of responses between unmedicated MDD and HC, with a high and flat-peaked time trajectory for the HC group for which the cortisol measurements started low and ended low (p value for negative second-order effect in the first four time measurements was .0093), a low and flatter cortisol versus time pattern for the unmedicated MDD group (p = .05 for positive second-order effect, relative to HC), and a combination of those two patterns for the medicated MDD group. MRI CORT did not differ across groups, MDD versus HC: F (1, 67) = 1.90, p = .17. Nor did MRI CORT significantly differ across groups after controlling for time of MRI onset, cortisol at MRI onset, and prior scan experience. No significant three-group differences were found for AUCg, AUCi, or MRI CORT.
Figure 2. Cortisol levels for the Trier Social Stress Test. The cortisol responses in the context of the Trier Social Stress Tests are shown. (Top) The patterns of cortisol for the two groups and (bottom) the patterns for the three groups.
Our second study aim was to examine whether MDD and HC adolescents differed in the interplay between amygdala structure and function with HPA axis functioning, after controlling for the confounders identified in the Statistical Analysis section. There was a significant interaction between right amygdala volume and participant group when predicting TSST AUCg for two-group comparison, F (1, 67) = 9.55, p = .003 (Figure 3), and three-group comparison, F (2, 65) = 5.24, p = .008 (Figure 4). There was also a significant interaction between right amygdala volume and participant group when predicting MRI CORT for the three-group comparison, F (2, 60) = 4.51, p = .01 (Figure 5), but not for the two-group comparison, F (1, 62) = 1.56, p = .21. Group-specific slopes (or pooled-group slopes, as appropriate) of volume with cortisol response are shown in Table 4. There was a significant interaction between left amygdala volume and participant group when predicting TSST AUCg, two group: F (1, 67) = 5.33, p = .02 (this pattern was similar to that found for right amygdala); three group: F (2, 65) = 2.81, p = .06 (this pattern was similar to that found for the right amygdala), but not when predicting MRI CORT AUCg, two group: F (1, 62) = 0.92, p = .34; three group: F (2, 60) = 2.69, p = .08. There were no significant findings for AUCi for any of these analyses.
Figure 3. Different patterns of correspondence between salivary cortisol levels during the Trier Social Stress Test (area under the curve with respect to ground) and the right amygdala volume were found for the depressed (major depressive disorder) and healthy control adolescents.
Figure 4. Different patterns of correspondence between salivary cortisol levels during the Trier Social Stress Test (area under the curve with respect to ground) and the right amygdala volume were found for the unmedicated depressed, the medicated depressed, and the healthy control adolescents.
Figure 5. Different patterns of correspondence between salivary cortisol levels during the MRI and the right amygdala volume were found for the unmedicated depressed, the medicated depressed, and the healthy control adolescents.
Table 4. Interplay regressions adjusted for reported/observed stress during the Trier Social Stress Test and intracranial volume
Note: AUCg, Area under the curve ground levels of salivary cortisol; No med MDD, unmedicated adolescents major depressive disorder; HC, healthy controls; Med MDD, medicated adolescents MDD. Subjects diagnosed with MDD with med MDD and no med MDD did not experience a significant left amygdala volume or right amygdala volume effects, unless it was relative to HC.
a The value for the difference from HC or slope = 0.
Next we examined whether MDD and HC adolescents differed in the interplay between amygdala function and HPA axis functioning after controlling for the confounders listed in Statistical Analysis. There was a significant association of right (p = .03; Figure 6), but not left (p = .09) amygdala functioning with AUCg cortisol response to the TSST for all groups combined. Pooled-group slopes of activation with cortisol response are shown in Table 4 and depicted on Figure 6. However, there were no significant interactions between (either right or left) amygdala functioning and participant group when predicting AUCg, AUCi, and MRI CORT, whether examining the two- or three-group comparisons.
Figure 6. (Color online) Correspondence between salivary cortisol levels during the Trier Social Stress Test (area under the curve with respect to ground) and the right amygdala functioning were found for the whole sample.
Discussion
To our knowledge, this is the first study to consider the interplay of the neural and hormonal systems of key relevance to stress regulation in depressed and well adolescents. Multilevel assessment included structural and functional measurement of the amygdala, hormonal response to stress in two different paradigms, and behavioral measures of stress. This study design allowed us to begin exploring whether biological threat system interplay differs by diagnosis across the sample of primarily female participants. Several important patterns emerged: (a) there were important ways in which the MDD and HC groups differed in responses to threat, (b) amygdala volume predicted HPA functioning differently in MDD and HC adolescents, (c) amygdala functioning predicted HPA functioning across the whole group of adolescent participants, and (d) the findings suggested some possible ways in which medication in MDD adolescents status might be important.
Our first aim was to consider possible differences in the biological stress system by assessment of key limbic structures and associated endocrine responses. Experimental evidence showed that the paradigms were generally effective in activating relevant systems; the amygdala was activated within the context of the emotion face-matching task, and cortisol levels were elevated in response to the TSST. There were also some differences for MDD and HC adolescents for biological responses to threat, specifically elevated responsiveness of the amygdala in the emotion face-matching task was found for the MDD adolescents. This is consistent with a recent study by Yang et al. (Reference Yang, Simmons, Matthews, Tapert, Frank and Max2010) using a nearly identical emotion-matching task in a smaller sample of adolescents with depression. Group differences were also noted in the pattern of HPA axis activation in the context of the TSST, with MDD participants showing a flatter, less reactive pattern than HC adolescents. This pattern of HPA axis response differed from some past studies investigating adolescents with internalizing problems (Klimes-Dougan et al., Reference Klimes-Dougan, Hastings, Granger, Usher and Zahn-Waxler2001) or with MDD (Rao et al., Reference Rao, Hammen, Ortiz, Chen and Poland2008). This flatter pattern may be more representative of MDD females (Stroud, Papandonatos, Williamson, & Dahl, Reference Stroud, Papandonatos, Williamson and Dahl2011), which is consistent with our 80% female sample. Consistent with previous work showing a flat reactivity pattern in adolescents with nonsuicidal self-injury (e.g., Kaess et al., Reference Kaess, Hille, Parzer, Maser-Gluth, Resch and Brunner2011), most of the adolescents in our sample (83%) displayed suicidal thoughts, suicidal behavior, and/or self-injury.
The second aim of this study was to examine the interplay of multiple levels of the biological stress system. An important contribution of this study is the evidence documenting that amygdala volume predicted HPA axis functioning differently for the MDD and HC adolescents. There is tentative evidence that in HC participants, large right and left amygdala volume may be reflective of a highly attuned system that is not only quick to marshal a stress response but also quick to abate arousal once the threat is past. In HC, amygdala volume was negatively correlated with cortisol values in the TSST. A similar pattern was noted in the MRI paradigm in the control group, with a larger amygdala size representing a highly responsive stress system, which exhibits a steep decline in cortisol levels over the course of the scan. Together these indices provide a coherent story of how amygdala volume is consistently related to a pattern of highly responsive HPA functioning under threat conditions in HC adolescents. The findings pertaining to the MDD participants provide a contrast to the interplay patterns of HC participants for cortisol levels within the context of both the TSST and the MRI. For the TSST, MDD participants (particularly unmedicated MDD participants) showed a positive association between amygdala volume and cortisol levels (AUCg). This pattern awaits replication given that research with adults has not shown a similar association across biological stress systems (e.g., Schuhmacher et al., Reference Schuhmacher, Mössner, Jessen, Scheef, Block and Belloche2012; Treadway et al., Reference Treadway, Grant, Ding, Hollon, Gore and Shelton2009).
The results of this study suggest that medication status may be important for further consideration. Because past studies typically evaluate either unmedicated adolescents or medicated adolescents (each having different limitations in internal or external validity), we thought it important to begin to explore possible medication effects for adolescents who continued to be symptomatic. There were some group differences noted based on MDD medication status for biological stress response or the interplay between the amygdala volume or functioning and the HPA axis. Most frequently medicated MDD adolescents tended to be more similar to the HC than the unmedicated MDD group. Medication regimens may influence HPA axis functioning by altering the glucocorticoid receptors and mineralocorticoid receptors as suggested by evidence that long-term antidepressant treatment has been shown to upregulate these receptors in the brain, normalizing the HPA axis cascade (Mason & Pariante, Reference Mason and Pariante2006). However, a recent meta-analysis of the treatment response literature with adults has brought into question the evidence that HPA axis functioning is altered with pharmacotherapy (McKay & Zakzanis, Reference McKay and Zakzanis2010). This remains a question for further inquiry in adolescent samples undergoing randomized control medication trials.
There was no evidence of group differences for the associations between amygdala activation and cortisol responses during threat. Across groups, high amygdala activation during matching of negative emotion faces was predictive of higher overall HPA axis functioning. A similar pattern was found (data not shown), even considering a more pure index of amygdala activation that was based on amygdala activation minus activation of the shape-matching task. Given that high activation of the amygdala and a flat cortisol pattern in response to stress was more characteristic of MDD adolescents, regardless of overall patterns it is possible that the positive association between amygdala and HPA axis activation may have different implications for the MDD and HC groups. It was somewhat surprising that, in this study, amygdala activation was related to cortisol levels that were measured on a different day (the TSST) but not to the more time-linked cortisol levels that were collected before and after the scan (MRI CORT). It is likely that engaging in the emotion face-matching task while undergoing a scan (even for those who had not been previously scanned) was not sufficiently potent to activate the HPA axis. By contrast, others have found that amygdala activation in healthy adults viewing traumatic imagery (pictures of the World Trade Center attack) was positively associated with cortisol reactivity that was collected prior to and following the brain scan (Cunningham-Bussel et al., Reference Cunningham-Bussel, Root, Butler, Tuescher, Pan and Epstein2009).
The approach taken in this study represents an important first step in trying to address the extraordinarily complex task of considering multiple and interacting levels of biological stress systems. We measured threat detection systems in the brain and the body. One strength of this study was that it considered both brain volume and function of the amygdala, a pivotal structure to threat detection. While brain structure is thought to underlie brain functioning (Hebb, Reference Hebb1949), rarely are strong associations between structure and function noted. Future multilevel research is also needed to expand the current work to investigate other key regions implicated in threat processing that have high concentrations of glucocorticoid receptors such as the hippocampus, anterior cingulate cortex, and prefrontal cortex (PFC). Future multilevel work should also expand upon the current finding to include assessment of a broader array of regulatory regions implicated in threat response. Models of emotion processing suggest reciprocal ventral and dorsal systems (Phillips et al., Reference Phillips, Drevets, Rauch and Lane2003). Some have suggested a reciprocal function of the PFC and limbic/hormonal response, although others have suggested increased recruitment of these regions in response to threat. For example, one study found an inverse correlation between rostral anterior cingulate cortex volume and average diurnal basal cortisol levels in depressed adults (Treadway et al., Reference Treadway, Grant, Ding, Hollon, Gore and Shelton2009). Another study found that individuals with the greatest increases in glucose metabolism in the medial PFC (BA 9 and BA 10) in response to stress were likely to have the lowest cortisol AUCg scores for the TSST (Kern et al., Reference Kern, Oakes, Stone, McAuliff, Kirschbaum and Davidson2008). Wang et al. (Reference Wang, Rao, Wetmore, Furlan, Korczykowski and Dinges2005) assessed stress responses using a modified TSST within the scan and salivary cortisol prior to and following the scan. They noted a positive association between right PFC and cortisol levels in healthy adults. Similarly, Jahn et al. (Reference Jahn, Fox, Abercrombie, Shelton, Oakes and Davidson2010) reported a positive association between brain metabolism in a number of ventral medial and limbic regions, including the subgenual PFC and pregenual PFC (BA 25/24), and cortisol under threat conditions in adolescent rhesus monkeys. These latter studies go against the reciprocal theory and suggest that as limbic systems respond, regulatory regions also respond to threat. This area of research is yet in its infancy and the conflicting finding highlight that these complex issues of interplay between systems require further investigation.
Another strength of this study was that two assessments of potential threat were assessed. Even though assessments were conducted on different days, in different settings, and with different experiments, the results indicated moderate correlations between the cortisol levels within the context of the MRI and the TSST paradigm. However, it was somewhat surprising to find that HPA axis activation in the MRI scanner was not related to amygdala functioning. Instead, cortisol levels during the TSST were positively correlated with amygdala functioning within the context of the emotion face-viewing task. In an ideal situation, proceeding with methodologies that allow researchers to temporally map stress responses of the neuroendocrine system onto neuronal activation would be desirable (see Dedovic et al., Reference Dedovic, Engert, Duchesne, Lue, Andrews and Efanov2010). However, the same stressor may be sufficiently potent to activate key limbic regions, but may need to be more intense and protracted to activate the peripheral stress response. In addition, a better understanding of HPA functioning under basal and stress conditions is warranted, for the results may differ considerably across these estimates of HPA axis functioning (e.g., Cunningham-Bussel et al., Reference Cunningham-Bussel, Root, Butler, Tuescher, Pan and Epstein2009; Root et al., Reference Root, Tuescher, Cunningham-Bussel, Pan, Epstein and Silbersweig2009). Additional refinements may include (a) assessing more saliva samples over a longer period of time to better evaluate recovery of the HPA axis, (b) using psychological stressors that elicit a stronger stress response within ethical constraints, (c) including a broader array of hormonal assessment (e.g., ACTH, DHEA, and DHEA-sulfate) given the cross-regulation of a broader hormonal network in adolescence (e.g., Marceau et al., Reference Marceau, Shirtcliff, Hastings, Klimes-Dougan, Zahn-Waxler and Dorn2014), and (d) including fewer task demands (e.g., limiting tasks demands before and after the TSST).
In the future it will be important to assess experiences and expressions of stress more comprehensively, particularly in the context of the MRI. Correspondence between experiences, expressions, or physiological responses to strong emotions/stress is often limited (e.g., Bauer, Quas, & Boyce, Reference Bauer, Quas and Boyce2002; Gross, Reference Gross1998). In this study, the findings show moderate correspondence between experiences and expression of emotion, r (74) = .44, p < .0001. Experience and expression of stress in the TSST was entered as a summary control variable in the primary study analyses because the results of this study are consistent with past research (Holsen et al., Reference Holsen, Lancaster, Klibanski, Whitfield-Babrieli, Cherkerzian and Ubka2013) showing that MDD adolescents experienced and expressed higher levels of stress than did HC adolescents. However, there was minimal evidence that experienced or expressed stress was related to neuronal or hormonal indexes. One exception was based on the question “Did the participant appear to be relieved to have completed the TSST?” (a high score represents “extremely relieved”). Experimenter ratings were significantly correlated with right, r (62) = .34, p = .008, and left, r (62) = .28, p = .02, amygdala activation and may reflect an important aspect threat recovery. It will be fruitful to continue to examine further how observed emotions maps onto brain activity.
In conclusion, the results of this study suggest differential links between amygdala structure and volume and neuroendocrine responses to stress in adolescents with MDD in comparison to HC. Further investigation is warranted to disentangle the processes more pertinent to normative maturation (e.g., functional connectivity across affective and cognitive control regions), the timing of onset and intensity of depression, and change over time with treatment. This is especially important in regard to how these clinical characteristics may be reflected in the adolescent's developing brain and extended physiological networks. This line of work would presumably serve to bolster current conceptualizations of the pathophysiology of depression early in development and may have implications for intervention efforts that take advantage of the plasticity that is still evident in neurobiological systems important during adolescence. Preventive intervention efforts may capitalize on this line of research, for it would be important to know in what ways altering the functioning of one system may influence stress functioning more broadly. Schuhmacher et al. (Reference Schuhmacher, Mössner, Jessen, Scheef, Block and Belloche2012) recently showed that larger amygdala volume prior to treatment was associated with a “normalization” of the hormonal stress response as measured by the dexamethasone /corticotrophin releasing hormone tests. Our team has identified that anomalous HPA axis functioning predicts treatment responses in young children with internalizing problems (Klimes-Dougan, Klingbeil, & August, Reference Klimes-Dougan, Klingbeil and August2009) and adolescents with depression (Gunlicks-Stoessel, Mufson, Cullen, & Klimes-Dougan, Reference Gunlicks-Stoessel, Mufson, Cullen and Klimes-Dougan2013). Considering these broader systems and their interplay is likely to provide a richer description of neurobiological functioning that may further enhance prediction of treatment response. Application of these results may also consider how intervention alters biological threat systems' functioning (Fisher, Stoolmiller, Gunnar, & Burraston, Reference Fisher, Stoolmiller, Gunnar and Burraston2007). Longitudinal research by Frodl et al. (Reference Frodl, Koutsouleris, Bottlender, Born, Jäger and Scupin2008) has found brain morphology changes in gray matter density in brain regions implicated in the threat system for depressed adults, and patients that remitted during the follow-up assessment showed less volume decline than those who failed to remit. There are certainly exciting possibilities for the field, and translation of this neurobiological research will be critical to advance current clinical practice and aid in the development of more effective neurobiologically informed interventions.
