Hostname: page-component-745bb68f8f-mzp66 Total loading time: 0 Render date: 2025-02-11T04:27:29.419Z Has data issue: false hasContentIssue false

The cortisol awakening response predicts major depression: predictive stability over a 4-year follow-up and effect of depression history

Published online by Cambridge University Press:  01 June 2012

S. Vrshek-Schallhorn*
Affiliation:
Department of Psychology, Northwestern University, Evanston, IL, USA
L. D. Doane
Affiliation:
Department of Psychology, Arizona State University, Tempe, AZ, USA
S. Mineka
Affiliation:
Department of Psychology, Northwestern University, Evanston, IL, USA
R. E. Zinbarg
Affiliation:
Department of Psychology, Northwestern University, Evanston, IL, USA The Family Institute at Northwestern University, Evanston, IL, USA
M. G. Craske
Affiliation:
Department of Psychology, University of California, Los Angeles (UCLA), Los Angeles, CA, USA
E. K. Adam
Affiliation:
School of Education and Social Policy, Northwestern University, Evanston, IL, USA Cells to Society Center, Institute for Policy Research, Northwestern University, Evanston, IL, USA
*
*Address for correspondence: S. Vrshek-Schallhorn, Ph.D., Department of Psychology, Northwestern University, 2029 Sheridan Road, Evanston, IL 60208, USA. (Email: suzanne.schallhorn@gmail.com)
Rights & Permissions [Opens in a new window]

Abstract

Background

The cortisol awakening response (CAR) has been shown to predict major depressive episodes (MDEs) over a 1-year period. It is unknown whether this effect: (a) is stable over longer periods of time; (b) is independent of prospective stressful life events; and (c) differentially predicts first onsets or recurrences of MDEs.

Method

A total of 270 older adolescents (mean age 17.06 years at cortisol measurement) from the larger prospective Northwestern-UCLA Youth Emotion Project completed baseline diagnostic and life stress interviews, questionnaires, and a 3-day cortisol sampling protocol measuring the CAR and diurnal rhythm, as well as up to four annual follow-up interviews of diagnoses and life stress.

Results

Non-proportional person-month survival analyses revealed that higher levels of the baseline CAR significantly predict MDEs for 2.5 years following cortisol measurement. However, the strength of prediction of depressive episodes significantly decays over time, with the CAR no longer significantly predicting MDEs after 2.5 years. Elevations in the CAR did not significantly increase vulnerability to prospective major stressful life events. They did, however, predict MDE recurrences more strongly than first onsets.

Conclusions

These results suggest that a high CAR represents a time-limited risk factor for onsets of MDEs, which increases risk for depression independently of future major stressful life events. Possible explanations for the stronger effect of the CAR for predicting MDE recurrences than first onsets are discussed.

Type
Original Articles
Copyright
Copyright © Cambridge University Press 2012

Introduction

Dysregulation of the hypothalamic–pituitary–adrenal (HPA) axis, and particularly the hormone cortisol, has been implicated in the etiology of major depressive disorder by a large body of animal research and cross-sectional research in humans (Chrousos & Gold, Reference Chrousos and Gold1992; Ehlert et al. Reference Ehlert, Gaab and Heinrichs2001). Several prospective studies of HPA functioning and depression have shown that cortisol levels during the morning (rather than average or total cortisol levels throughout the day), around the peak of cortisol's diurnal rhythm, prospectively predict depressive episodes or symptoms (Harris et al. Reference Harris, Borsanyi, Messari, Stanford, Brown, Cleary, Shiers and Herbert2000; Halligan et al. Reference Halligan, Herbert, Goodyer and Murray2007; Goodyer et al. Reference Goodyer, Croudace, Dudbridge, Ban and Herbert2010). One aspect of morning cortisol that has received particular attention is the cortisol awakening response (CAR) – the pronounced (on average, 50–60%) increase in cortisol level that peaks approximately 30 to 40 min after waking (Pruessner et al. Reference Pruessner, Wolf, Hellhammer, Buske-Kirschbaum, Von Auer, Jobst, Kaspers and Kirschbaum1997). The CAR is frequently operationalized as the difference between the level of salivary cortisol sampled immediately after waking and the level sampled 30 to 40 min after waking (Clow et al. Reference Clow, Hucklebridge, Stalder, Evans and Thorn2010). Experimental research indicates that the CAR represents a true response to waking, rather than an underlying diurnal rhythm (Wilhelm et al. Reference Wilhelm, Born, Kudielka, Schlotz and Wüst2007). Though the precise function of the CAR has not been fully characterized, current theory suggests that it serves to mobilize resources to meet the perceived demands of the upcoming day (e.g. Adam et al. Reference Adam, Hawkley, Kudielka and Cacioppo2006; Fries et al. Reference Fries, Dettenborn and Kirschbaum2009).

The CAR is thought to be an important indicator of risk for major depressive episodes (MDE) in part because it captures a large and rapid increase on top of already elevated waking cortisol levels (Adam et al. Reference Adam, Doane, Zinbarg, Mineka, Craske and Griffith2010), potentially reaching physiological levels capable of activating relatively low-affinity glucocorticoid receptors implicated in the etiology of major depressive disorder (e.g. Holsboer, Reference Holsboer2000). Recently, we showed that the CAR prospectively predicted the occurrence of MDEs across a 1-year follow-up period (odds ratio 3.0, p = 0.04; Adam et al. Reference Adam, Doane, Zinbarg, Mineka, Craske and Griffith2010). Other basal cortisol indices, namely the diurnal slope and total waking cortisol (area under the curve with respect to ground; AUCg), were not related to MDE onsets during this follow-up period. Participants in this study were a subset of those from the Northwestern-UCLA Youth Emotion Project (YEP), a two-site longitudinal study of risk for emotional disorders (Zinbarg et al. Reference Zinbarg, Mineka, Craske, Griffith, Sutton, Rose, Nazarian, Mor and Waters2010). Since that time, additional follow-up waves have been conducted and the number of participants in the cortisol subsample of the YEP reporting any MDE over the follow-up period has increased from 18 to 42, permitting additional analyses over a longer time period.

Thus, using data from the same study over an extended 4-year follow-up period, the present study examined three questions. First, does the CAR predict MDE onsets beyond 1 year (what we refer to as a question of predictive stability)? Several prospective studies suggest that morning or average daytime cortisol predicts depressive symptoms over 2 to 3 years (Halligan et al. Reference Halligan, Herbert, Goodyer and Murray2007; Shirtcliff & Essex, Reference Shirtcliff and Essex2008); to our knowledge, no one has examined the CAR as a predictor of MDE onset over similar time periods, nor has anyone explicitly examined the predictive stability of the CAR. We will use the CAR in relation to MDE onsets over an extended period of time by using a CAR × time interaction variable to detect decay in predictive power. Knowledge regarding the stability of this effect could inform our understanding of the CAR, and has practical significance for the design of prospective research studies. The robustness of the CAR effect and the CAR × time interaction will also be considered. In particular, models will examine whether the CAR and its interaction with time predict MDE above and beyond other known predictors of MDE, including the personality trait neuroticism (e.g. Kendler et al. Reference Kendler, Kuhn and Prescott2004) and a history of MDE or anxiety disorders at the time of cortisol measurement (e.g. Boland & Keller, Reference Boland, Keller, Gotlib and Hammen2009). Models will also consider whether the CAR is a significant predictor of MDE above and beyond the effects of possible confounding variables, including fatigue and negative affect at the time of cortisol sampling (e.g. Adam et al. Reference Adam, Hawkley, Kudielka and Cacioppo2006), and chronic and episodic life stress in the year prior to cortisol sampling (Adam & Gunnar, Reference Adam and Gunnar2001; Gerritsen et al. Reference Gerritsen, Geerlings, Beekman, Deeg, Penninx and Comijs2010).

Second, do elevations in the CAR increase vulnerability to prospective major stressful life events? Major stressful life events are known to precipitate MDEs, often within the month after the occurrence of an event (Kendler et al. Reference Kendler, Karkowski and Prescott1999). Our research has shown that the CAR's ability to predict depression onset was independent of an index of moderate to severe episodic stressful events occurring in approximately the year prior to cortisol measurement (Adam et al. Reference Adam, Doane, Zinbarg, Mineka, Craske and Griffith2010). What is unknown is whether elevations in the CAR associated with increased risk for MDEs represent enhanced vulnerability to future stressful events, or instead whether such CAR elevations increase risk for MDEs independently of these stressful events. To our knowledge, prior research has demonstrated unique effects of elevated morning cortisol and of stressful life events in the month prior to depressive onset (Goodyer et al. Reference Goodyer, Tamplin, Herbert and Altham2000), but research has not tested whether morning cortisol peaks (or an elevated CAR) significantly interact with stressful life events in predicting MDEs.

Third, does the CAR predict first onsets and recurrences of MDE differentially? While our prior report after a 1-year follow-up presented descriptive data supporting that both first onsets and recurrences of MDE were predicted by elevations in the CAR (Adam et al. Reference Adam, Doane, Zinbarg, Mineka, Craske and Griffith2010), an increased number of MDE cases in the present report permits a formal statistical examination of whether the CAR interacts with MDE history to predict new onsets and recurrences differentially, where one prospective MDE per person is considered. It might be that CARs of similar magnitude are more predictive of MDE risk for first onsets or recurrences. To our knowledge, the present study is the first to examine each of these questions, which are intended to clarify the relationship of the CAR to onsets of depressive episodes.

Method

Participants

Detailed information regarding recruitment and demographics of the larger YEP sample is reported by Zinbarg et al. (Reference Zinbarg, Mineka, Craske, Griffith, Sutton, Rose, Nazarian, Mor and Waters2010). In brief, over 1900 high school juniors were screened for neuroticism using the Revised Eysenck Personality Questionnaire Neuroticism Scale (Eysenck et al. Reference Eysenck, Eysenck and Barrett1985). Those scoring high on this measure were oversampled, comprising 59% of the final group of 627 (432 females) who gave consent to participate in a longitudinal study and who completed the baseline interviews and questionnaires described below. Participants were asked to repeat all these measures annually.

Of these participants, 491 individuals were randomly selected for invitation to participate in a cortisol sampling protocol; a total of 344 participants (250 females) completed this protocol, 70.1% of those invited. Participants were excluded from analyses for one or more of the following reasons: corticosteroid medication (n = 13); insufficient cortisol data (n = 29); no follow-up interviews after the baseline interview (n = 18); current major depression at baseline (n = 17); current or past post-traumatic stress disorder at baseline (n = 2), which is associated with lower cortisol levels (for a meta-analysis, see Meewisse et al. Reference Meewisse, Reitsma, De Vries, Gersons and Olff2007); or at any time during the study being diagnosed with dysthymic disorder (n = 8), clinically significant psychotic symptoms (n = 3) or bipolar disorder I or II (n = 7). Individuals with these three diagnoses were excluded a priori from analyses in order to yield a cleaner comparison group. After exclusions, a total of 270 participants remained, yielding 75 males and 195 females in the current analyses. The sample comprises more females than males both because females were more likely than males to agree to participate if invited, and because we oversampled for high levels of neuroticism, a trait on which females score higher on average (Costa et al. Reference Costa, Terracciano and McCrae2001). Participants were Caucasian (49%), African American (10%), Hispanic or Latino (18%), Asian (4%), or multiple/other races and ethnicities (19%). Hollingshead's socio-economic status index (Hollingshead, Reference Hollingshead1975) was determined from information gathered at the baseline interview.

Procedures

Assessment of psychopathology

The presence or absence of clinically significant MDEs was ascertained at each assessment using the Structured Clinical Interview for DSM-IV, non-patient edition (SCID; First et al. Reference First, Spitzer, Gibbon and Williams2001). All diagnoses reported here were clinically significant manifestations. The baseline SCID assessed lifetime diagnoses, and four subsequent annual follow-up SCIDs assessed diagnoses of psychopathology occurring during the interim since the participant's previous assessment. Interviewers were blind to the results of previous assessments. Inter-rater reliability was assessed for approximately 10% of all SCIDs administered in the larger study. Across the five assessment periods, κ for individual interviewers' diagnoses of MDEs, adjusted due to departure from equiprobability (R. Zinbarg et al. unpublished observations), ranged from 0.82 to 0.94 (mean = 0.89, s.d. = 0.05).

Life stress assessment

Baseline chronic stress, which is used as a covariate in models to rule out the possibility that stress accounts for CAR's effects on MDE, was measured using the Life Stress Interview (LSI; Hammen et al. Reference Hammen, Adrian, Gordon, Burge, Jaenicke and Hiroto1987). Chronic stress over the past year was measured across 10 life domains: best friend relationship; social circle; romantic relationships; family relationships; academics; work; finances; neighborhood conditions; physical health; and family's health. Ratings were assigned by the interviewer for each domain on a scale from 1 to 5 in half-point increments, with 1 representing the best possible circumstances, and 5 representing the worst possible circumstances. Scores in each domain were standardized, then averaged to create a chronic stress composite score. Inter-rater reliability (intra-class correlations) for the baseline chronic stress composite score was 0.87 within site and 0.84 across sites.

Episodic life stress was also assessed at annual interviews using the LSI (Hammen, Reference Hammen1991). Episodic events reported by participants in the LSI were presented by the interviewer to a team of independent raters blind to the participant's diagnoses and response to the event. Event severity ratings were assigned by the consensus of the independent rating team, ranging from 1 (no significant threat or negative implications) to 5 (a very severe event, maximal negative impact or threat) in half-point increments. Inter-rater reliability (intra-class correlations) cross-site for the five interviews examined ranged from 0.69 to 0.76 (mean = 0.72, s.d. = 0.03). Here, events are deemed to be ‘major’ if assigned a severity rating of 2.5 or greater, reflecting moderately to severely stressful events. To covary episodic life stress prior to cortisol, the severity scores of all major events in the year prior were summed to create a single episodic stress score. To examine prospective stressful events following cortisol, the occurrence of a major event in a given month was coded dichotomously (absent = 0, present = 1) in ‘person-month’ datasets, described in more detail in the ‘Analytic plan’ section. If a follow-up interview was missed, the subsequent SCID covered the period since the prior interview, but the LSI covered only the past 12 months due to concerns about reliability of recalling stressful life events occurring over longer durations. As such, when prospective stressful life events are incorporated into statistical models, months for which LSI data were not available were excluded from analyses.

Personality assessment

As previously described (Zinbarg et al. Reference Zinbarg, Mineka, Craske, Griffith, Sutton, Rose, Nazarian, Mor and Waters2010), the personality trait neuroticism was measured using four separate scales measuring neuroticism or one of its facets: the Revised Eysenck Personality Questionnaire Neuroticism Scale (Eysenck et al. Reference Eysenck, Eysenck and Barrett1985), the International Personality Item Pool (IPIP) NEO Revised Personality Inventory (Goldberg, Reference Goldberg, Mervielde, Deary, De Fruyt and Ostendorf1999), the Behavioral Inhibition Scale from the Behavioral Inhibition Scale/Behavioral Activation Scale (Carver & White, Reference Carver and White1994) and the Big Five Mini-Markers Neuroticism Scale (Saucier, Reference Saucier1994). These were standardized and then averaged to yield a composite neuroticism score (Zinbarg et al. Reference Zinbarg, Mineka, Craske, Griffith, Sutton, Rose, Nazarian, Mor and Waters2010).

Cortisol assessment

The cortisol collection protocol used here has been described in detail elsewhere (Adam et al. Reference Adam, Doane, Zinbarg, Mineka, Craske and Griffith2010; Doane et al. Reference Doane, Adam, Mineka, Zinbarg, Craske and Griffithin press). In brief, cortisol was collected by passive drool on three consecutive weekdays, six times per day: immediately upon waking, 40 min after waking, three semi-random times throughout the day spaced to avoid mealtimes and signaled by a programmed watch beep, and, finally, at bedtime. Prior to each sample, participants completed a diary report including self-reported mood, which yields the measures of negative affect and fatigue used here. Participants were instructed not to eat, drink or brush their teeth in the 30 min prior to the expected morning and bedtime sample collection times. Several health variables associated with HPA axis functioning were measured by questionnaire or in daily diaries: asthma, hormonal contraceptives, time of waking and bedtime on cortisol measurement days, and daily nicotine use.

Once returned to the laboratory, samples were stored at −20 °C and later shipped to Trier, Germany, where they were assayed in duplicate using time-resolved fluorescent-detection immunoassay (Dressendörfer et al. Reference Dressendörfer, Kirschbaum, Rohde, Stahl and Strasburger1992). Intra-assay variation ranged from 4.0% to 6.7%, while inter-assay variation ranged from 7.1% to 9.0%. Raw cortisol data were logarithmically transformed due to positive skew prior to analyses; raw salivary cortisol values are presented in graphs for ease of interpretation. In addition to cortisol levels at waking, 40 min after waking, and bedtime, three indices of diurnal cortisol were calculated for use in analyses: the CAR (value 40 min after waking minus value at waking), the diurnal rhythm slope (a regression line fit to all measures between waking and bedtime except the measure taken at 40 min after waking), and the average cortisol level (the area under the curve of all six cortisol measures with respect to ground, AUCg, divided by total time awake). Indices represent the mean from the 3 days of cortisol collection.

Analytic plan

Cox regression models (continuous time survival analyses; Cox, Reference Cox1972) examining the primary research questions utilized a person-month dataset, permitting unique values for stressful life events and MDE onsets for each month. In order to examine the predictive stability of the CAR, non-proportional hazards models were employed. Unlike proportional hazards models, such models do not assume that the hazard associated with a particular predictor is constant across time, permitting the examination of interactions with time (Singer & Willett, Reference Singer and Willett2003). Interactions of the CAR × time were calculated by multiplying median-centered time in months by the CAR value. All Cox models yield hazard ratios (HRs) and their 95% confidence intervals for each predictor. Continuous predictors are standardized to facilitate interpretation of HRs: a one-unit increase in a predictor corresponds to one standard deviation.

Stressful life events and MDE onsets were dated to the nearest month after cortisol measurement based on dates collected during the SCID and LSI interviews. Whether life events indicated as occurring in the same month as MDE onsets in fact occurred before the MDE onsets (i.e. temporal precedence) was carefully scrutinized. When the order of the event and the MDE onset was indeterminate, or when the event clearly occurred after the MDE onset, that event was excluded from consideration. A maximum of one MDE onset per participant was included in the model, whether that MDE represented a first onset or a recurrence of MDE; that is, once a participant experienced an MDE onset after cortisol measurement, their remaining months were not included. Participants without such diagnoses remained in the model until the month corresponding to their last follow-up SCID assessment or a maximum of 54 months.

In models for the prediction of MDEs, unless otherwise specified below, the following sets of variables are entered as four blocks into Cox regression models: demographic variables (i.e. gender, SES, African American race, Hispanic ethnicity), health variables (i.e. hours of sleep, time of waking, nicotine use, asthma, hormonal contraceptives), cortisol variables (CAR, slope, and AUCg), and finally the interaction of CAR × time. For all analyses, p values ⩽0.05 are considered statistically significant.

Results

Descriptive statistics for all variables are presented in Table 1. Zero-order correlations for all variables are presented in Table 2. Participants who were included in analyses did not differ significantly from those who were excluded in: gender [χ2(1) = 0.79, n.s.]; minority race/ethnicity [χ2(1) = 3.29, n.s.]; or baseline Hollingshead socio-economic status level [F(1, 334) = 0.66, n.s.]. A total of 42 MDEs were observed (35 in females, seven in males; 27 first onsets, 15 recurrences) after participation in the cortisol protocol, with time to onset ranging from less than 1 month to 53 months. A total of 10 723 person-months were available for analysis of MDE onsets for the primary model.

Table 1. Descriptive statistics of raw data for cortisol levels and covariates

s.d., Standard deviation; MDE, major depressive episode; SES, Hollingshead socio-economic status index; AUCg, area under the curve with respect to ground, or total waking cortisol; CAR, cortisol awakening response; slope, rate of cortisol decline from wake-up to bedtime.

a Cortisol levels at times 1–6 are given in μg/dl and represent average cortisol levels over the 3 days of sample collection at the following times: sample immediately upon waking (S1); sample 40 min after waking (S2); sample approximately 3 h post-waking (S3); sample approximately 8 h post-waking (S4); sample approximately 12 h post-waking (S5); and sample immediately before bedtime (S6). All analyses were conducted using natural log-transformed values.

Table 2. Zero-order Pearson correlations

SES, Hollingshead socio-economic status index; MDE, major depressive episode; AUCg, area under the curve with respect to ground, or total waking cortisol; CAR, cortisol awakening response.

a Variable names are condensed to conserve space; refer to Table 1 for full variables names.

b ‘MDE history’ (variable 12) refers to whether or not (yes = 1, no = 0) a participant was diagnosed with having experienced at least one clinically significant MDE prior to cortisol measurement.

* Statistically significant correlation (p < 0.05).

The CAR predicting MDEs

In the final primary model, the CAR was a significant prospective predictor of MDEs; at the median of time (21 months), any one standard deviation increase in the CAR is associated with approximately double the risk for MDE onset. The CAR × time interaction was also significant; the predictive power of the CAR degraded (about 3% per month) with successive months (Table 3, model 1). There was no evidence for a decline in the underlying risk of depressive onsets per person-month (hazard function) over time. To estimate how long the CAR predicts MDEs at a statistically significant level, the time intercept was re-centered to reflect each 6-month time point since the cortisol assessment, providing estimates of effect size at each time point. As implied by the linear decline in prediction over time, maximal effects were seen shortly after cortisol measurement, with a one standard deviation increase in the CAR corresponding to a roughly 3.5-fold increase in risk for MDE onset (Fig. 1). Results also indicate that the CAR significantly predicts MDEs in our sample up to 30 months, or 2.5 years, but not significantly so at 36 months and beyond.

Fig. 1. Model-estimated hazard ratios illustrate the observed cortisol awakening response (CAR) × time interaction. Values are hazard ratios for z-scored baseline ln-transformed CAR in 6-month intervals after baseline. MDE, Major depressive episode. * p < 0.05, ** p < 0.01, n.s., non-significant (p > 0.05).

Table 3. Results of several Cox regression models for MDE onsetsa

MDE, Major depressive episode; s.e., standard error; CI, confidence interval; CAR, cortisol awakening response; AUCg, area under the curve with respect to ground, or total waking cortisol.

a All models include CAR, slope, AUCg and the CAR × time interaction as well as additional covariates for gender, Hollingshead socio-economic status index, African American race, Hispanic ethnicity, tobacco use, asthma, time of waking, hours of sleep, and hormonal contraceptive use. Continuous variables were entered in the model as z scores to aid interpretation; dichotomous variables (i.e. gender) were entered untransformed.

To examine the robustness of the CAR main effect and CAR × time interaction following the addition of other known predictors of depression and potential confounding variables, four separate models were tested with the respective addition of: (a) baseline neuroticism (Table 3, model 2); (b) anxiety disorder at baseline and history of an MDE; (c) baseline episodic and chronic life stress; and (d) negative affect and fatigue at the time of cortisol participation. Although in each case the additional block made a significant contribution to predicting onsets of MDE [all block changes in χ2 (df ⩽ 2) ⩾13.60, all p's <0.05], the CAR and its interaction with time remained significant predictors of MDEs with very similar HRs as in the primary model. Including time interactions for the additional covariates did not alter the findings.

Does the CAR interact with prospective major stressful life events?

When person-months without LSI data were excluded, a total of 9645 person-months were available, including 39 MDE onsets. While a main effect of major stressful life events occurring in the past month robustly predicted MDE onsets, the CAR × major stressful life events interaction was not significant (Table 3, model 3). The CAR and its interaction with time remained significant predictors of MDEs independently of major stressful life events.

Does the CAR predict first onsets and recurrences of MDE differentially?

To evaluate whether the CAR similarly predicted both first onsets and recurrences of an MDE when one prospective MDE per person was considered, an interaction term was created between MDE history at the time of the cortisol sampling (no = 0, yes = 1) and the CAR. When this MDE history × CAR interaction was added, the simple main effect of CAR, indicating its relationship to first onsets of MDE, remained significant, as did the CAR × time interaction. The MDE history × CAR interaction effect was also significant (Table 3, model 4). Whereas elevated CARs significantly predicted first onsets, they predicted recurrences with greater strength. Diurnal rhythms of raw salivary cortisol levels are presented by diagnostic group (never depressed, first-onset MDE, recurrent MDE, or history of MDE only) in Fig. 2. This interaction effect was not a matter of differing CAR magnitudes in the first-onset and recurrence groups, as CAR magnitudes did not significantly differ as a function of MDE history [F(1, 40) = 0.22, n.s.]. Further, the number of months between cortisol measurement and MDE onset did not differ significantly between first onsets (mean = 22.19, s.d. = 12.64) and recurrences (mean = 21.13, s.d. = 16.59) [F(1, 40) = 0.053, p = 0.82].

Fig. 2. Baseline 3-day averages of raw salivary cortisol level (μg/dl) across the day for never-depressed participants, participants who experienced a first onset of major depressive episode (MDE) within the 4-year follow-up, participants who experienced a recurrence of MDE within the 4-year follow-up, and participants with a history of MDE but no recurrence during the 4-year follow-up. Analyses were conducted with natural logarithmically transformed values.

Discussion

Using data from 270 participants from the larger Youth Emotion Project over four annual follow-ups, this is the first study to explicitly test the extent to which the CAR predicts MDE onsets consistently over an extended period of time. Results showed that at 6 months after cortisol sampling, any one-unit increase in the baseline standardized CAR was associated with a greater than 3-fold increase in risk per month for MDE onset. This risk then decays, remaining a statistically significant predictor in our sample for approximately 2.5 years (out of 4 years), but not significantly so after such time. These effects persisted despite covariance of other known predictors of MDEs and possible confounding variables. Thus, the CAR appears to serve as a time-limited risk factor for depression. Studies interested in predicting onsets of MDE over a several-year period may benefit from including a baseline CAR measure, but longitudinal studies using the CAR as a predictor of MDE onsets over more extended periods of time may benefit from using repeated measures of the CAR. This is also the first study to examine whether the CAR predicts MDE onsets in interaction with prospectively occurring major stressful life events. The results indicate that both the CAR and major events predict MDEs uniquely and independently of one another.

In addition, the results showed that when a single prospective MDE onset per person is considered, the CAR significantly predicts both first onsets and recurrences of depression, but that the CAR predicts recurrences with greater strength than it predicts first onsets. Both individuals experiencing first onsets and those experiencing recurrences have CARs of similar magnitude (Fig. 2), but the CAR appears particularly important in predicting whether someone with a history of depression will experience a recurrence. Because the CAR prospectively predicts first onsets, these data suggest that an elevated CAR should not be construed solely as a ‘scar’ of depression (Cowen, Reference Cowen2010), i.e. a marker that is only present following an MDE, which may be indicative of damage caused by the episode. The observed effect is consistent with several possible explanations. Individuals with a history of MDE who have recurrences during the follow-up period experienced earlier first onsets than those whose first onsets were observed during the follow-up period. Earlier age of first onset is associated with a worse course of illness (Boland & Keller, Reference Boland, Keller, Gotlib and Hammen2009) and may account for the observed effect. Alternatively, Post's (Reference Post1992) kindling theory has suggested that within-person change occurs as a result of a first MDE, which may erode resilience and increase an individual's propensity for later episodes. While the present data do not examine within-person change in the relation of the CAR to MDE, a stronger relationship between the CAR and recurrences of MDE compared with first onsets of MDE may be consistent with such within-person change.

Taken together, these findings suggest that the CAR and not other aspects of the diurnal cortisol rhythm predict future MDE onsets. Current theory posits that the CAR functions to prepare the individual to rise to the expected demands of the upcoming day (Adam et al. Reference Adam, Hawkley, Kudielka and Cacioppo2006; Fries et al. Reference Fries, Dettenborn and Kirschbaum2009); thus, the presence of the CAR can be seen as adaptive. However, frequent large CAR values may indicate excessive perceived demands and/or insufficient perceived resources (e.g. competencies and social support) to meet those demands. In such instances, large CAR increases may be indicative of struggle above and beyond anything captured by measuring mood. By contrast, flattened diurnal cortisol slopes may serve as a marker of the day's levels of negative emotions such as tension and anger, leading to elevated evening cortisol levels (Adam et al. Reference Adam, Hawkley, Kudielka and Cacioppo2006). Flattened diurnal slopes are also suggested to serve as a ‘scar’ marker for prior depression (Doane et al. in press), suggesting they may indicate allostatic load associated with accommodating recent psychopathology, among other stressors.

Although the present data do not speak to the biological mechanism by which the CAR is related to MDE risk, several possibilities exist. As noted earlier, the CAR reflects both the highest and the most rapidly increasing cortisol levels of the day. Over time, an elevated CAR may contribute to more intense occupation and eventual down-regulation of glucocorticoid receptors, a pattern which has been implicated in depression (Holsboer, Reference Holsboer2000). Further, the CAR appears to be regulated by a range of factors, including extra-pituitary efferent pathways from the suprachiasmatic nucleus to the adrenal influencing sensitivity to adrenocorticotropic hormone (ACTH), and a role for hippocampal functioning, perhaps particularly pre-awakening (Clow et al. Reference Clow, Hucklebridge, Stalder, Evans and Thorn2010). Such pathways may uniquely influence the CAR, and not the diurnal slope or evening levels, leading to the differential importance of the CAR for prediction of MDE onsets. Finally, it is possible that the CAR does not play a causal role in precipitating MDE onsets, but rather serves as a proxy marker for risk conferred by unknown but correlated neural processes.

Limitations

Despite strengths of the present study, including the use of structured diagnostic and life stress interviews, the longitudinal design, the naturalistic measurement of the CAR over a 3-day period, and inclusion of potentially confounding covariates, this study has several limitations. First, due to concerns about feasibility with 17-year-olds and participant burden, only two samples were collected to define the CAR (0 and 40 min after waking), rather than the original approach, which relied on samples taken every 10–15 min for 1 h after waking (Pruessner et al. Reference Pruessner, Wolf, Hellhammer, Buske-Kirschbaum, Von Auer, Jobst, Kaspers and Kirschbaum1997). Second, although the examination of stressful life events in person-month format represents an advancement over other approaches, we have examined major stressful life events in an aggregate fashion. It may be that CAR increases vulnerability to certain types of major stressful events, but not all major events. Third, we cannot rule out a developmental interpretation whereby the CAR's interaction with time indicates that the CAR was particularly important during the age range closest to the CAR measurement (i.e. roughly ages 17 to 19 years), rather than a decay in predictive ability. Fourth, we would have liked to examine whether the CAR's relation to prospective MDE onsets differed between males and females, i.e. a gender × CAR interaction; however, the number of onsets among males remains too low to address this question (seven cases). Finally, given that the present study oversampled for participants with high levels of neuroticism (59%), it may be important to replicate these findings among a representative community sample or a sample selected based on a different risk factor, such as family history of depression.

Conclusions

This is the first study to demonstrate that the CAR predicts depression onsets for 2.5 years of a 4-year follow-up, initially serving as a strong predictor but ‘decaying’ in its predictive ability with the passage of time. That is, the CAR appears to function as a time-limited risk factor for depression. Further, the CAR does not significantly enhance vulnerability to major stressful life events. Third, while the CAR predicts both first onsets and recurrences of depression when a maximum of one prospective MDE onset per person is considered, it predicts recurrences with greater strength than first onsets. Research which examines factors contributing to individual differences in the CAR, and within-person change in the CAR over time may be particularly important for understanding the role of the HPA axis in the etiology and course of major depressive disorder.

Acknowledgements

The present study was supported by a two-site grant from the National Institute of Mental Health to S.M. and R.E.Z. (no. R01-MH065652) and to M.G.C. (no. R01-MH065651), by a William T. Grant Foundation Scholars Award to E.K.A., a Faculty Fellowship (E.K.A.) and Graduate Fellowship (L.D.D.) from the Institute for Policy Research at Northwestern University, and a Postdoctoral NRSA from the National Institute of Mental Health to S.V.-S. (no. F32-MH091955). The content is the responsibility of the authors and does not necessarily represent the official views of the National Institute of Mental Health.

Declaration of Interest

None.

References

Adam, E, Doane, L, Zinbarg, R, Mineka, S, Craske, M, Griffith, J (2010). Prospective prediction of major depressive disorder from cortisol awakening responses in adolescence. Psychoneuroendocrinology 35, 921931.CrossRefGoogle ScholarPubMed
Adam, E, Gunnar, M (2001). Relationship functioning and home and work demands predict individual differences in diurnal cortisol patterns in women. Psychoneuroendocrinology 26, 189208.CrossRefGoogle ScholarPubMed
Adam, E, Hawkley, L, Kudielka, B, Cacioppo, J (2006). Day-to-day dynamics of experience–cortisol associations in a population-based sample of older adults. Proceedings of the National Academy of Sciences USA 103, 1705817063.CrossRefGoogle Scholar
Boland, R, Keller, M (2009). Course and outcome of depression. In Handbook of Depression (ed. Gotlib, I. and Hammen, C.), pp. 2344. Guilford Press: New York.Google Scholar
Carver, CS, White, TL (1994). Behavioral inhibition, behavioral activation, and affective responses to impending reward and punishment: the BIS/BAS scales. Journal of Personality and Social Psychology 67, 319333.CrossRefGoogle Scholar
Chrousos, G, Gold, P (1992). The concepts of stress and stress system disorders: overview of physical and behavioral homeostasis. JAMA 267, 12441252.CrossRefGoogle ScholarPubMed
Clow, A, Hucklebridge, F, Stalder, T, Evans, P, Thorn, L (2010). The cortisol awakening response: more than a measure of HPA axis function. Neuroscience and Biobehavioral Reviews 35, 97103.CrossRefGoogle ScholarPubMed
Costa, PT, Terracciano, A, McCrae, RR (2001). Gender differences in personality traits across cultures: robust and surprising findings. Journal of Personality and Social Psychology 81, 322331.CrossRefGoogle ScholarPubMed
Cowen, P (2010). Not fade away: the HPA axis and depression. Psychological Medicine 40, 14.CrossRefGoogle ScholarPubMed
Cox, D (1972). Regression models and life-tables. Journal of the Royal Statistical Society Series B (Methodological) 34, 187220.CrossRefGoogle Scholar
Doane, LD, Adam, EK, Mineka, S, Zinbarg, R, Craske, M, Griffith, J (in press). Are flatter diurnal cortisol rhythms associated with major depression and anxiety disorders in late adolescence? The role of life stress and daily negative emotion. Development and Psychopathology.Google Scholar
Dressendörfer, R, Kirschbaum, C, Rohde, W, Stahl, F, Strasburger, C (1992). Synthesis of a cortisol–biotin conjugate and evaluation as a tracer in an immunoassay for salivary cortisol measurement. Journal of Steroid Biochemistry and Molecular Biology 43, 683692.CrossRefGoogle Scholar
Ehlert, U, Gaab, J, Heinrichs, M (2001). Psychoneuroendocrinological contributions to the etiology of depression, posttraumatic stress disorder, and stress-related bodily disorders: the role of the hypothalamus–pituitary–adrenal axis. Biological Psychology 57, 141152.CrossRefGoogle Scholar
Eysenck, SBG, Eysenck, HJ, Barrett, P (1985). A revised version of the psychoticism scale. Personality and Individual Differences 6, 2129.CrossRefGoogle Scholar
First, M, Spitzer, R, Gibbon, M, Williams, J (2001). Structured Clinical Interview for DSM-IV-TR Axis I Disorders-Non-Patient Edition. New York State Psychiatric Institute, Biometrics Research Department: New York.Google Scholar
Fries, E, Dettenborn, L, Kirschbaum, C (2009). The cortisol awakening response (CAR): facts and future directions. International Journal of Psychophysiology 72, 6773.CrossRefGoogle ScholarPubMed
Gerritsen, L, Geerlings, M, Beekman, A, Deeg, D, Penninx, B, Comijs, H (2010). Early and late life events and salivary cortisol in older persons. Psychological Medicine 40, 15691578.CrossRefGoogle ScholarPubMed
Goldberg, L (1999). A broad-bandwidth, public domain, personality inventory measuring the lower-level facets of several five-factor models. In Personality Psychology in Europe (ed. Mervielde, I., Deary, I., De Fruyt, F. and Ostendorf, F.), pp. 728. Tilburg University Press: Tilburg, The Netherlands.Google Scholar
Goodyer, I, Croudace, T, Dudbridge, F, Ban, M, Herbert, J (2010). Polymorphisms in BDNF (Val66Met) and 5-HTTLPR, morning cortisol and subsequent depression in at-risk adolescents. British Journal of Psychiatry 197, 365371.CrossRefGoogle ScholarPubMed
Goodyer, I, Tamplin, A, Herbert, J, Altham, P (2000). Recent life events, cortisol, dehydroepiandrosterone and the onset of major depression in high-risk adolescents. British Journal of Psychiatry 177, 499504.CrossRefGoogle ScholarPubMed
Halligan, S, Herbert, J, Goodyer, I, Murray, L (2007). Disturbances in morning cortisol secretion in association with maternal postnatal depression predict subsequent depressive symptomatology in adolescents. Biological Psychiatry 62, 4046.CrossRefGoogle ScholarPubMed
Hammen, C (1991). Generation of stress in the course of unipolar depression. Journal of Abnormal Psychology 100, 555561.CrossRefGoogle ScholarPubMed
Hammen, C, Adrian, C, Gordon, D, Burge, D, Jaenicke, C, Hiroto, D (1987). Children of depressed mothers: maternal strain and symptom predictors of dysfunction. Journal of Abnormal Psychology 96, 190198.CrossRefGoogle ScholarPubMed
Harris, T, Borsanyi, S, Messari, S, Stanford, K, Brown, G, Cleary, S, Shiers, H, Herbert, J (2000). Morning cortisol as a risk factor for subsequent major depressive disorder in adult women. British Journal of Psychiatry 177, 505510.CrossRefGoogle ScholarPubMed
Hollingshead, A (1975). Four Factor Index of Social Status. Yale University: New Haven, CT.Google Scholar
Holsboer, F (2000). The corticosteroid receptor hypothesis of depression. Neuropsychopharmacology 23, 477501.CrossRefGoogle ScholarPubMed
Kendler, K, Karkowski, L, Prescott, C (1999). Causal relationship between stressful life events and the onset of major depression. American Journal of Psychiatry 156, 837841.CrossRefGoogle ScholarPubMed
Kendler, K, Kuhn, J, Prescott, C (2004). The interrelationship of neuroticism, sex, and stressful life events in the prediction of episodes of major depression. American Journal of Psychiatry 161, 631636.CrossRefGoogle ScholarPubMed
Meewisse, ML, Reitsma, JB, De Vries, GJ, Gersons, BPR, Olff, M (2007). Cortisol and post-traumatic stress disorder in adults: systematic review and meta-analysis. British Journal of Psychiatry 191, 387392.CrossRefGoogle ScholarPubMed
Post, R (1992). Transduction of psychosocial stress into the neurobiology of recurrent affective disorder. American Journal of Psychiatry 149, 9991010.Google ScholarPubMed
Pruessner, J, Wolf, O, Hellhammer, D, Buske-Kirschbaum, A, Von Auer, K, Jobst, S, Kaspers, F, Kirschbaum, C (1997). Free cortisol levels after awakening: a reliable biological marker for the assessment of adrenocortical activity. Life Sciences 61, 25392549.CrossRefGoogle ScholarPubMed
Saucier, G (1994). Mini-markers: a brief version of Goldberg's unipolar big-five markers. Journal of Personality Assessment 63, 506516.CrossRefGoogle ScholarPubMed
Shirtcliff, EA, Essex, MJ (2008). Concurrent and longitudinal associations of basal and diurnal cortisol with mental health symptoms in early adolescence. Developmental Psychobiology 50, 690703.CrossRefGoogle ScholarPubMed
Singer, J, Willett, J (2003). Applied Longitudinal Data Analysis: Modeling Change and Event Occurrence. Oxford University Press: New York.CrossRefGoogle Scholar
Wilhelm, I, Born, J, Kudielka, B, Schlotz, W, Wüst, S (2007). Is the cortisol awakening rise a response to awakening? Psychoneuroendocrinology 32, 358366.CrossRefGoogle Scholar
Zinbarg, R, Mineka, S, Craske, M, Griffith, J, Sutton, J, Rose, R, Nazarian, M, Mor, N, Waters, A (2010). The Northwestern-UCLA Youth Emotion Project: associations of cognitive vulnerabilities, neuroticism and gender with past diagnoses of emotional disorders in adolescents. Behaviour Research and Therapy 48, 347358.CrossRefGoogle ScholarPubMed
Figure 0

Table 1. Descriptive statistics of raw data for cortisol levels and covariates

Figure 1

Table 2. Zero-order Pearson correlations

Figure 2

Fig. 1. Model-estimated hazard ratios illustrate the observed cortisol awakening response (CAR) × time interaction. Values are hazard ratios for z-scored baseline ln-transformed CAR in 6-month intervals after baseline. MDE, Major depressive episode. * p < 0.05, ** p < 0.01, n.s., non-significant (p > 0.05).

Figure 3

Table 3. Results of several Cox regression models for MDE onsetsa

Figure 4

Fig. 2. Baseline 3-day averages of raw salivary cortisol level (μg/dl) across the day for never-depressed participants, participants who experienced a first onset of major depressive episode (MDE) within the 4-year follow-up, participants who experienced a recurrence of MDE within the 4-year follow-up, and participants with a history of MDE but no recurrence during the 4-year follow-up. Analyses were conducted with natural logarithmically transformed values.