|Home | About | Journals | Submit | Contact Us | Français|
Studies of individuals who do not meet criteria for major depressive disorder (MDD) but with subclinical levels of depressive symptoms may aid in the identification of neurofunctional abnormalities that possibly precede and predict the development of MDD. The purpose of this study was to evaluate relations between subclinical levels of depressive symptoms and neural activation patterns during tasks previously shown to differentiate individuals with and without MDD.
Functional magnetic resonance imaging (fMRI) was used to assess neural activations during active emotion regulation, a resting state scan, and reward processing. Participants were twelve females with a range of depressive symptoms who did not meet criteria for MDD.
Increased depressive symptom severity predicted (1) decreased left midfrontal gyrus activation during reappraisal of sad stimuli; (2) increased right midfrontal gyrus activation during distraction from sad stimuli; (3) increased functional connectivity between a precuneus seed region and left orbitofrontal cortex during a resting state scan; and (4) increased paracingulate activation during non-win outcomes during a reward-processing task.
These pilot data shed light on relations between subclinical levels of depressive symptoms in the absence of a formal MDD diagnosis and neural activation patterns. Future studies will be needed to test the utility of these activation patterns for predicting MDD onset in at-risk samples.
Despite a large body of research addressing impaired neural functioning in major depressive disorder (MDD), relatively little is known about the neurofunctional characteristics of individuals without MDD but with subclinical levels of depressive symptoms. Studies of individuals who do not meet criteria for MDD but with subclinical levels of depressive symptoms may yield a number of insights: first, they may suggest potential neurobiologic markers of those at risk for MDD [1,2]; second, they may shed light on hereditary and environmental influences on depressive temperament ; and third, they may suggest avenues of inquiry regarding the neurobiology of MDD risk and resilience, and thereby inform treatment and preventative intervention approaches [4,5]. The goal of the present pilot study was to investigate linkages between regional brain activation patterns and subclinical depressive symptoms while participants were engaged in emotion regulation, resting state, and reward processing paradigms that have previously been shown to differentiate MDD and nondepressed samples. Below we outline the rationale for selecting each of these three paradigms.
Deficits in emotion regulation (ER) are thought to be central to the core features of MDD . Human and preclinical studies suggest that the prefrontal cortex influences activity in the limbic system in a top-down manner [7,8], and numerous investigations have demonstrated abnormal recruitment of prefrontal cortical regions while individuals with MDD actively regulate responses to emotional stimuli [9,10]. For example, Johnstone and colleagues  reported that when instructed to reappraise negatively valenced images, individuals with MDD demonstrated impaired prefrontal cortical inhibition of limbic regions as well as a divergent pattern of covariation between ventromedial prefrontal cortex and amygdala activity. Similarly, Beauregard and colleagues  reported right anterior cingulate and insular cortex hyperactivation while individuals with MDD down-regulated their emotional responses to sad images. Finally, Kanske and colleagues  reported anterior cingulate and lateral orbitofrontal cortex hyperactivation during emotion regulation in individuals with remitted MDD , supporting the framework that altered neural mechanisms of emotion regulation may indeed represent not only a state marker of MDD illness, but potentially a trait marker of MDD risk as well.
Second, we investigated resting state functional brain connectivity. A mostly midline network of correlated intrinsic brain activity is active during rest and deactive during goal-directed tasks [13,14]. This network has been called both the “baseline state” and the “default mode network (DMN)” [15,16]. The DMN regulates self-referential activities, including evaluating the salience of internal and external cues, remembering the past, and planning the future [13,17]. Individuals with MDD are characterized by altered patterns of DMN activity during emotion processing tasks as well as at rest. Sheline and colleagues  reported increased activity in MDD within DMN regions (i.e., ventromedial prefrontal cortex, anterior cingulate, and lateral parietal and temporal cortices), suggesting that MDD is characterized by a failure to down-regulate DMN activity during emotion processing. Another study by the same research group  reported increased intrinsic connectivity in MDD between a precuneus seed region and dorsolateral prefrontal cortex (the so-called “dorsal nexus” for overconnected brain regions in MDD) (see also  and ). Given the centrality of the DMN for self-referential activities, abnormal DMN activity in MDD has been suggested to reflect excessive self-focus accompanied by a decreased ability to attend to cognitive tasks .
Finally, we measured frontostriatal responses during a reward processing task because multiple studies have reported altered neural mechanisms of reward processing in MDD (for a review, see ). For example, Smoski and colleagues  reported that individuals with MDD demonstrated orbital frontal cortex hyperactivity during reward selection that predicted depressive symptom severity and caudate nucleus hypoactivation during reward anticipation. Pizzagalli and colleagues  found that individuals with MDD showed decreased putamen responses during reward anticipation and decreased nucleus accumbens and caudate nucleus responses to monetary gains. Finally, Knutson and colleagues  reported increased anterior cingulate activation during reward anticipation in MDD. Taken together, these studies highlight not only that MDD is characterized by altered neural activation patterns during reward processing, but that patterns of brain activation differences in MDD are contingent on the temporal phase of the reward response. For this reason, we used a task that allowed for the dissociation of responses during reward selection, reward anticipation, and reward outcome phases of reward responding.
In summary, the purpose of the present study was to extend the literature addressing altered neural mechanisms of emotion regulation, DMN connectivity, and reward processing in MDD by examining the covariation between levels of depressive symptoms and patterns of brain activation in individuals with subclinical levels of depressive symptoms but without MDD diagnoses. Predictions were informed by the extant MDD literature reviewed above [6,19,22]: we hypothesized relations between subclinical depressive symptom severity and (1) lateral prefrontal cortex activity during ER; (2) precuneus-dorsolateral prefrontal cortex connectivity during the resting state scan; (3) lateral prefrontal cortex activity during reward selection; (4) striatal activity during reward anticipation; and (5) medial prefrontal cortex activity during reward outcomes. Given the higher rates of MDD in women  and because gender moderates responses to emotional images [26,27], this pilot study examined only females, a strategy that is consistent with other neuroimaging studies examining MDD , depression risk  and depressive traits .
Participants consented to protocols approved by the local Human Investigations Committees at both the University of North Carolina Chapel Hill and Duke University Medical Center and were recruited from a database maintained at the Duke-UNC Brain Imaging and Analysis Center and by flyers posted in campus and medical center locations. To identify individuals with moderate but not severe depression symptoms, interested respondents first completed a modified the Beck Depression Inventory-II [BDI,  via a secure web survey. Next, potential participants completed an in-person screening session that included the full BDI and administration of the Structured Clinical Interview for DSM-IV, Patient Version, with Psychotic Screen (SCID) . The final sample was recruited based on a stratification strategy to recruit equal numbers of individuals with low (i.e., 0–9), medium (i.e., 10–18), and high (i.e., 19+) BDI scores. Exclusion criteria included: 1) current anxiety, substance abuse, or mood disorder beyond dysthymia (as assessed by the SCID), 2) history of mania, 3) currently taking psychoactive medication, 4) magnetic resonance imaging contraindicated (e.g., metal in body), 5) history of neurological injury or disease, and 6) pregnancy. The final sample was twelve females (mean age: 23.3±4.3; 3 African American; 3 Asian, and six Caucasian) who participated in the scan session.
Participants completed a training session prior to the scan day at which they were taught strategies to distract themselves and reappraise their interpretations while viewing sad stimuli. For the distract condition, participants were trained to think of something positive that was unrelated to what was depicted in the images, such as a fond memory, family member, or an anticipated event. For the reappraise condition, participants were trained to either give the story suggested in the picture a happy ending or alternate meaning, or to tell themselves that no one close to them was affected by the depicted situation. As a control strategy, participants were instructed to experience the image as they normally would without trying to increase or decrease their emotional response (the “view” condition).
On the day of the scan, participants reviewed task instructions and completed the BDI again. During scanning, participants completed the three paradigms described below in a single scan session. Due to technical problems, three participants returned to complete the emotion regulation task at a separate imaging session, at which time the BDI was readministered. Analyses included only BDI scores obtained on the day that each task was completed. The mean BDI scores in final analyses were 11.33 (SD=10.21, Range=0 to 26, Skewness=0.10) on the day of the resting state and reward processing scans, and 9.92 (SD=9.02, Range=0 to 24, Skewness=0.29) on the day of the ER scan.
The ER task consisted of four 4-minute 24-second runs with nine sad trials and three neutral trials each (see the top of Figure1). The stimuli used in the task are described elsewhere [33,34]. During each 22-second trial, participants first fixated for 6seconds and then saw a sad or neutral image for 11seconds with an auditory instruction presented 3–6seconds after image onset. Sad images were presented on 75% of trials and were accompanied by an instruction to “reappraise,” “distract,” or “view” the image. Neutral images, on the other 25% of trials, were accompanied by an instruction to “view” the image. Following image offset, participants reported their ER success on a 4-point scale anchored by “not successful” and “very successful,” or they indicated they felt “no emotion” (trials with “no emotion” were not included in analyses).
Participants were instructed to close their eyes and relax but not fall asleep for the duration of the 6-minute resting state scan.
The WOF is described in greater detail elsewhere . Briefly, three 12-minute runs with 46 trials each were composed of selection, anticipation, and outcome phases (see the bottom of Figure1). During the selection phase (3seconds), participants viewed a segmented wheel of fortune and indicated their selection. During the anticipation phase (3.5-7.5seconds), participants rated how sure they were of winning on a scale of 1 to 6 (anchored by “unsure” and “sure”). During the outcome phase (4seconds), participants viewed the amount won (“$0” for not winning), the total cumulative amount won in the run, and rated how they felt on a pictorial Likert scale. Responses were contrasted with a control condition that lacked monetary stimuli but required active responding.
Scanning was performed on a GE 3 Tesla Signa Excite HD scanner with 50-mT/m gradients (General Electric, Waukesha, Wisconsin, USA). Head movement was restricted by foam cushions. An eight-channel head coil was used for parallel imaging. A high resolution T1-weighted anatomical image was acquired using a 3D fast SPGR pulse sequence (68 slices, TR=500ms; TE=20ms; FOV=24cm; image matrix=2562; voxel size=0.9375 × 0.9375 × 1.9mm) aligned in the near axial plane defined by the anterior and posterior commissures. Whole brain functional images consisting of 30 slices parallel to the AC-PC plane were collected using a BOLD-sensitive gradient-echo sequence with spiral-in k-space sampling and SENSE encoding to take advantage of the 8-channel coil (TE: 30ms; FOV: 25.6cm; isotropic voxel size: 4mm3; SENSE factor=2). Functional imaging sessions began with four discarded RF excitations to allow for steady state equilibrium. TR’s were 1500ms for the ER run and 2000ms for the WoF and resting state runs. Images were presented using ePrime (Psychology Software Tools Inc., Pittsburgh, PA) for the WOF task, and with CIGAL  for the ER task.
Functional data were preprocessed using FEAT 5.92 within FSL 4.0.4 (Oxford Centre for Functional Magnetic Resonance Imaging of the Brain (FMRIB), Oxford University, U.K.). Preprocessing was applied in the following steps: (i) brain extraction for non-brain removal on the T1-weighted image , (ii) motion correction using MCFLIRT , (iii) spatial smoothing using a Gaussian kernel of FWHM 5mm, (iv) mean-based intensity normalization of all volumes by the same factor, and (v) high-pass filtering (60s) . Functional images were co-registered to the brain-extracted T1-weighted image in native space, and the anatomical image was normalized to Montreal Neurological Institute (MNI) standard stereotaxic space. The same transformation matrices used for anatomical-to-standard transformations were then used for functional-to-standard space transformations of co-registered functional images. All registrations were carried out using an intermodal registration tool [37,39]. Voxel-wise temporal autocorrelation was estimated and corrected using FMRIB's Improved Linear Model .
For ER and WOF data, onset times of events were used to model a signal response containing a regressor for each response type, which was convolved with a double-gamma function to model the predicted hemodynamic response. A general linear model (GLM) approach generated whole brain images of parameter estimates and variances, representing average signal change from baseline. The same high-pass filtering applied to the functional data was applied to the GLM. Group-wise activation images were calculated by a mixed effects analysis using Bayesian estimation techniques, FMRIB Local Analysis of Mixed Effects [FILM, . Group-level models, with mean-centered BDI scores as regressors, were created by combining contrast maps from individual subjects. Finally, scatterplots of correlations between brain activation and BDI score are presented for illustrative purposes .
Analysis of resting state data used the following seed-based connectivity approach: (1) raw functional scans were temporally band-pass filtered (0.1<f<.08); (2) as in , the mean signal intensity timecourses from voxels inside a spherical ROI (radius=10mm) in the precuneus (MNI coordinates 0, -56, 30) were extracted for all participants in native space; (3) average signal intensity timecourses derived from the previous step were then used as unconvolved regressors in a GLM that included motion parameters, resulting in clusters significantly correlated with the timecourse of the seed region; (4) a group-level model was created by combining contrast maps from individual participants.
Given the importance of protecting against Type II errors in this exploratory pilot study, activation maps were thresholded at Z>2.3 with a spatial extent of at least ten voxels, consistent with published recommendations [43,44] and our prior MDD studies [25,45-47].
There were no differences in self-reported ER success across the four ER conditions (Neutral-View, Sad-View, Sad-Distract, and Sad-Reappraise), p’s>0.80. There was a significant inverse relation between BDI scores and success ratings only in the Sad-Reappraise condition, r=−0.71, p<.01, reflecting poorer self-reported ER in those with higher BDI scores.
The top of Figure2 depicts results during Sad-Reappraise>Sad-View (see Table1). Areas that showed activation included the paracingulate gyrus (Z Max=4.36), right middle frontal gyrus (Z Max=3.07), and bilateral temporal gyrus (Z Max=4.14 for left and 3.81 for right). BDI scores were inversely related to activation magnitudes in the left middle frontal gyrus (Z Max=2.98) and right orbitofrontal cortex (Z Max=2.82) (see Table2).
The bottom of Figure2 depicts group-average results during Sad-Distract>Sad-View (see also Table3). Areas that showed activation included the cingulate gyrus (Z Max=4.54), right supramarginal gyrus (Z Max=3.63), and right lateral occipital cortex (Z Max=3.39). BDI scores were directly related to activation magnitudes in the right middle frontal gyrus (Z Max=2.69), right precentral gyrus (Z Max=3.09), and left superior temporal gyrus (Z Max=2.87) (see Table4).
Brain regions that covaried with the precuneus seed are depicted in (Figure3 and Table5). Intrinsic functional connectivity was observed between the precuneus seed and left postcentral gyrus (Z Max=5.19), right orbital frontal cortex (Z Max=3.61), and bilateral temporal pole (Z Max=4.38 for right and 3.43 for left). BDI scores directly predicted the strength of functional connectivity between the precuneus seed and left orbitofrontal cortex (Z Max=4.42), left thalamus (Z Max=3.58), right occipital cortex (Z Max=3.96), bilateral parietal lobule (Z Max=5.81 for right and 3.91 for left), bilateral postcentral gyrus (Z Max=4.71 for right and 3.89 for left), and bilateral precentral gyrus (Z Max=4.55 for right and 4.59 for left) (see Table6).
Average (SD) valence ratings on a scale of 1–6 after non-win trials was 2.87 (1.19) (1 depicted a very sad face and 6 depicted a neutral face). Average (SD) valence after win trials was 4.94 (0.86) (1 depicted a neutral face and 6 depicted a very happy face). BDI scores were not significantly correlated with valence ratings after win trials, r=−0.44, p=0.15 or after non-win trials, r=−0.03, p =0.93.
The top of Figure4 depicts results during reward selection (see also Table7). Areas that showed activation included dorsal paracingulate gyrus (Z Max=3.45), anterior cingulate gyrus (Z Max=2.97), and right frontal orbital cortex (Z Max=3.15). BDI scores did not predict activation magnitudes in any region during this phase of the task.
The bottom of Figure4 depicts results during reward anticipation (see also Table8). Areas that showed activation included bilateral nucleus accumbens (Z Max=3.44 for right and 2.79 for left,), right putamen (Z Max=2.67), and right thalamus (Z Max=2.8). BDI scores did not predict activation magnitudes in any regions during this phase of the task.
The top of Figure5 depicts results during win outcomes (see also Table9). Areas that showed activation included the left putamen (Z Max=3.1), right thalamus (Z Max =2.89), and right amygdala (Z max=3.02). BDI scores did not predict activation magnitudes in any canonical reward regions during win outcomes. The bottom of Figure5 depicts average results during non-win outcomes (see also Table10). Areas that showed activation included the left paracingulate gyrus (Z max=3.61), left putamen (Z Max=2.80), and left thalamus (Z Max=2.96). BDI scores were directly related to activation during non-win outcomes in the left paracingulate gyrus (Z Max=2.62) and the right caudate nucleus (Z Max=2.96) (See Table11).
The goal of the present study was to use fMRI to investigate linkages between patterns of brain activation and connectivity and subclinical levels of depressive symptoms. Emotion regulation, resting state, and reward processing paradigms were used that have been shown to differentiate MDD and control samples and to be related to core features of MDD. During the ER task, participants with higher BDI scores reported significantly poorer self-reported ER success. This pattern is consistent with the findings of Beauregard and colleagues , who found that individuals with MDD reported relatively greater difficulty down-regulating feelings of sadness. Consistent with reports from nonclinical samples , group-level activation during reappraisal of sad stimuli was evident in the right middle frontal gyrus and bilateral temporal gyrus. Depressive symptom severity was related to activation magnitudes in right middle frontal gyrus during emotion regulation, though the laterality and direction of this finding was contingent on the type of ER strategy: left middle frontal gyrus activation was inversely related to activation during reappraisal, whereas right middle frontal gyrus activation was directly related to activation during distraction. These findings are consistent with neuroimaging data indicating that MDD is characterized by increased right middle frontal activation during reappraisal . More broadly, these results suggest that brain activation during ER is related to the magnitude of subclinical levels of depressive symptoms.
During the resting state scan, intrinsic connectivity was observed between the precuneus seed region and the temporal cortex and anterior cingulate, similar to findings previously reported in nonclinical samples . There was a significant positive correlation between depression symptom severity and the strength of functional connectivity between the precuneus seed region and the thalamus and left orbital frontal cortex. Given that MDD is characterized by increased DMN connectivity both during emotional processing paradigms and at rest [18,19], the present findings suggest that increased DMN connectivity may characterize individuals with increasing levels of subclinical depressive symptoms, possibly reflecting heightened self-referential activities and ruminative cognitions.
During all phases of the reward processing task, participants demonstrated frontostriatal brain activation patterns similar to published reports from nonclinical samples, including dorsal paracingulate cortex activation during reward selection, striatal activation during reward anticipation, and medial prefrontal activation during reward outcomes [25,35]. Depressive symptom severity was not related to activation magnitudes in canonical reward regions during reward selection, reward anticipation, or win outcomes. However, depressive symptom severity was directly related to paracingulate gyrus activation during non-win outcomes. This pattern of findings is noteworthy when considered alongside the results of Pizzagalli and colleagues , who found greater differences between MDD and control samples during reward outcomes than during reward anticipation. There is also evidence that pediatric MDD is characterized by altered neural mechanisms of reward anticipation following reward outcomes , suggesting that neural responses to reward outcomes may impact future experiences with rewards. More generally, this pattern of results suggests that altered neural correlates of non-win reward outcome may be a more sensitive marker of subclinical depressive symptoms than neural responses during reward selection, reward anticipation or win outcomes.
The present study has a number of methodological limitations that should be considered when interpreting results. The relatively small sample and inclusion of only females suggests the need for replication in larger and more heterogeneous samples. Monthly variability in mood symptoms tied to menstrual phase and use of hormone contraceptives were not assessed, nor was participant handedness 1. Finally, brain imaging analyses were conducted at uncorrected thresholds to protect against Type II errors in this pilot study, and future studies with larger samples will allow for replication with corrected statistical thresholds.
Despite these limitations, this study suggests that subclinical depressive symptom severity is selectively related to regional brain activation patterns in key emotion regulation and reward processing regions and to DMN connectivity implicated in MDD. Although the three domains addressed here (i.e., emotion regulation, resting state connectivity, and reward processing) were assessed via separate tasks, there is clear conceptual overlap between these three domains. For example, prefrontal cortical recruitment is anomalous in MDD during regulation of rewarding stimuli  and DMN activity in MDD is abnormal during emotion regulation . In this regard, the present findings may represent a conservative estimate of the effects of subclinical depressive symptoms on neural responses during tasks that simultaneously address these constructs.
Given that individuals with subclinical levels of depressive symptoms are at heightened risk for developing MDD [1,2], these results suggest candidate neural markers that may reflect trait vulnerability to MDD. Because MDD represents a highly prevalent, chronic, and costly disorder, identification of at risk individuals may provide critical windows for preventative interventions that represents a public health priority . Future research using longitudinal and high-risk designs will be needed to more fully evaluate etiologic questions about MDD, MDD risk, and MDD resilience.
aWe thank an anonymous reviewer for raising these points.
BDI: Beck Depression Inventory-II; BOLD: Blood Oxygen Level Dependent; DMN: Default Mode Network; ER: Emotion Regulation; fMRI: Functional magnetic resonance imaging; GLM: General linear model; MDD: Major Depressive Disorder; MNI: Montreal Neurological Institute; WOF: Wheel of Fortune; SCID: Structured Clinical Interview for DSM-IV Patient Version.
The authors declare that they have no competing interests.
JNF, MJS, and GSD conceived of and conducted the study and drafted the manuscript. RVK, BF, JM, JB, and CP participated in study design and data analysis and helped to draft the manuscript. All authors read and approved the final manuscript.
The pre-publication history for this paper can be accessed here:
This research was supported by a Duke-UNC Brain Imaging and Analysis Center Post-Baccalaureate fellowship to J. Felder and by a UNC Research Council Award to G. Dichter. M. Smoski was supported by NIMH T32-MH070448, a NARSAD Young Investigator award, and a career development award from Duke University Medical Center, NICHD K12 HD043446. G. Dichter was supported by NARSAD Young Investigator awards and by K23 MH081285. F. McClernon was supported by K23DA017261. The funding sources had no further role in study design, in the collection, analysis and interpretation of data; in the writing of the report; and in the decision to submit the paper for publication.
The authors would like to thank Todd Harshbarger and Syam Gadde for assistance with image analysis, Prue Cuper for SCID supervision, Dr. Nan-kuei Chen and Dr. James Voyvodic for helpful discussions regarding analysis of resting state data, MRI technologists Susan Music, Natalie Goutkin, and Luke Pool for assistance with data acquisition, and Dr. Allen Song for assistance with various aspects of this research. Assistance for this study was also provided by the Neuroimaging Core of the UNC Carolina Institute for Developmental Disabilities (P30 HD03110).
Publication of this article was funded by the University of Colorado Boulder Libraries Open Access Fund.