Dissociable roles of reward prediction error in the contrasting mood dynamics of depression and anxiety
Figures
Experimental protocol.
(A) Study outline. This study examines how depression and anxiety influence mood fluctuations. The first experiment assesses the bifactor structure that disentangles shared and specific components of depression and anxiety in the psychometric dataset (N=901). The second experiment tests associations between depression- and anxiety-specific traits and RPE-induced mood fluctuations using the questionnaire battery and a gambling task with momentary mood ratings in the laboratory dataset (N=44). The third and fourth experiments replicate the second experiment in online samples (online dataset 1, N=747; online dataset 2, N=235). The fifth experiment tests whether these findings generalize to a clinical dataset of patients with affective disorders (N=61). (B) Mapping between the three factors and 128 items in the bifactor model. (C) Factor loadings of items on the general factor, anxiety-specific factor, and depression-specific factor. (D) Orthogonality among the general, anxiety-specific, and depression-specific factor scores across datasets with complete questionnaire data. (E) Mean correlations between factor scores and questionnaire scores. Overall, the general factor showed high correlations with all questionnaires. The depression-specific factor correlated most strongly with TAIdep and MASQad, whereas the anxiety-specific factor correlated most strongly with the remaining questionnaires. These correlational results supported the bifactor structure of anxiety and depression. (F) Gambling task design. On each trial, participants were asked to choose between a certain option and a gamble option. Once an option was selected, the corresponding outcome was displayed in the center of the screen. The cumulative score was always shown in the upper-right corner. Every two or three trials, participants were asked to complete a self-paced rating of “How happy are you at this moment?” on a slider scale ranging from 0 (very unhappy) to 100 (very happy). (G) Temporal dynamics of happiness ratings for representative individuals with high (top 25%) and low (top 75%) mood variation, healthy datasets, and the clinical dataset. (H) Momentary mood model. (I) Results of momentary mood model for each dataset. MASQaa, the subscale of anxious arousal in the Mood and Anxiety Symptoms Questionnaire; TAIanx, the subscale of anxiety in the Trait Anxiety Inventory; CESD, Center for Epidemiologic Studies Depression Scale; BDI, Beck Depression Inventory; BFIn, the subscale of neuroticism in the Big Five Inventory; PSWQ, Penn State Worry Questionnaire; MASQad, the subscale of anhedonic depression in the Mood and Anxiety Symptoms Questionnaire; TAIdep, the subscale of depression in the Trait Anxiety Inventory; CR, certain reward; EV, expected value; RPE, reward prediction error.
Gambling rate in each dataset.
All datasets showed the same gambling rate pattern, as well as in line with literature.
Model-agnostic mood results for all datasets.
All datasets showed the same initial mood overall before the task, mean mood, and mood variation across the task.
Replication of the classic mood effects.
(A) Happier after winning than losing. (B) Mood drift across the task.
Parameter recovery for the winning model (M1) for the healthy and clinical datasets.
Parameter recovery analysis showed high correlations between the matched simulated and real parameters and low correlations between the mismatched simulated and real parameters.
Replication of Rutledge et al., 2017’s findings using BDI (N=1026).
Depression symptom measured by BDI was negatively correlated with the baseline mood parameter. Data was combined from laboratory dataset, online dataset 1, and online dataset 2.
Correlations of factor scores from the winning model (the bifactor model) with questionnaire scores for each dataset.
Dissociable associations of depression and anxiety with mood fluctuations.
Correlations of depression- and anxiety-specific factor scores with mood variation and RPE-related mood sensitivity () for the laboratory dataset (A, E), online dataset 1 (B, F), online dataset 2 (C, G), and the combined dataset (N=1026; D, H). (I) Reduced RPE-related mood sensitivity statistically mediated the association between depression-specific scores and reduced mood fluctuations. (J) Increased RPE-related mood sensitivity statistically mediated the association between anxiety-specific scores and increased mood fluctuations. (K) Among the mood sensitivity parameters, the depression-specific factor was selectively associated with decreased RPE-related mood sensitivity. (L) Among the mood sensitivity parameters, the anxiety-specific factor was selectively associated with increased RPE-related mood sensitivity. Regression coefficients are shown as bootstrapped mean ± SE. CR, certain reward; EV, expected value; RPE, reward prediction error. *p<0.05.
Clinical validation of reduced RPE-related mood sensitivity in depression.
(A, B) Correlations of depression and anxiety factor score with mood variation and mood parameter of RPE () for the clinical dataset. (C) The mediation model among depression, , and mood variation in the clinical population. The regression coefficients were represented by mean ± SE, which were estimated by bootstrap. RPE, reward prediction error; *p<0.05.
Consistent negative associations of the common factor with mood baseline and gain-domain risk attitude across datasets.
(A) Negative association between the common factor and the baseline mood parameter. (B) Negative association between the common factor and gain-domain risk attitude. Regression coefficients are shown as mean ± SE. *p<0.05.
Tables
Basic demographic details.
| Psychometric dataset (N=901) | Laboratory dataset (N=44) | Online dataset 1 (N=747) | Online dataset 2 (N=235) | Clinical dataset (N=116) | |
|---|---|---|---|---|---|
| Gender (female) | 676 | 18 | 500 | 143 | 89 |
| Age | 22.04±2.10 | 20.05±1.70 | 20.90±2.41 | 21.67±2.52 | 16.40±3.86 |
| MASQaa | 24.70±8.13 | 22.98±6.14 | 24.16±7.60 | 23.92±7.62 | 43.51±10.51 |
| TAIanx | 16.06±4.86 | 15.00±4.95 | 15.69±4.97 | 15.99±5.34 | 23.35±5.44 |
| CESD | 35.28±10.72 | 33.30±10.73 | 34.82±10.74 | 35.06±11.41 | 55.64±13.04 |
| BDI | 11.38±9.35 | 7.91±8.11 | 10.04±8.94 | 10.26±9.18 | 22.66±17.91 |
| BFIn | 33.99±8.82 | 30.52±8.26 | 32.82±9.37 | 32.79±9.81 | 43.75±9.31 |
| PSWQ | 48.28±11.61 | 46.89±12.52 | 48.25±12.15 | 48.36±13.29 | 59.18±12.87 |
| MASQad | 61.49±14.00 | 60.23±16.63 | 58.07±15.30 | 59.57±15.71 | 77.16±15.39 |
| TAIdep | 29.39±5.47 | 27.32±6.16 | 27.61±6.27 | 27.85±6.34 | 35.97±4.91 |
| Data collection periods | 2021.6.22–2021.6.28 | 2022.6.2–2022.7.5 | 2022.4.24–2022.5.22 | 2022.7.22–2022.7.26 | 2023.3.16–2023.8.4 |
-
Descriptive data are presented as mean ± SD.
-
MASQaa, the subscale of anxious arousal in the Mood and Anxiety Symptoms Questionnaire; TAIanx, the subscale of anxiety in the Trait Anxiety Inventory; CESD, Center for Epidemiologic Studies Depression Scale; BDI, Beck Depression Inventory; BFIn, the subscale of neuroticism in the Big Five Inventory; PSWQ, Penn State Worry Questionnaire; MASQad, the subscale of anhedonic depression in the Mood and Anxiety Symptoms Questionnaire; TAIdep, the subscale of depression in the Trait Anxiety Inventory.
Exploratory analysis for the common factor.
| Models | Healthy datasets (n=1026) | Clinical dataset (n=116) | ||
|---|---|---|---|---|
| BIC | Results | BIC | Results | |
| β0 ~ common | –830.23 | b=–0.03, t=–6.54, p<0.001, 95% CI = [-0.042–0.022] | –113.74 | b=–0.03, t=–2.05, p=0.043, 95% CI = [-0.052–0.001] |
| β0 ~ common + common2 | –823.74 | –110.91 | ||
| αgain ~ common | 551.76 | b=–0.03, t=–3.08, p=0.002, 95% CI = [-0.049–0.011] | 101.7 | b=–0.11, t=–3.24, p=0.002, 95% CI = [-0.172–0.041] |
| αgain ~ common + common2 | 553.73 | 106.12 | ||
-
Note: For healthy datasets, we added datasets as random variables. β0, the mood baseline parameter; αgain, risk attitude for gain.
Psychometric model comparisons.
| Models | CFI | TLI | RMSEA | SRMR |
|---|---|---|---|---|
| Bifactor model | 0.989 | 0.989 | 0.046 | 0.061 |
| Three-factor model | 0.794 | 0.789 | 0.200 | 0.205 |
| High-order model | 0.981 | 0.981 | 0.060 | 0.072 |
Mood model comparisons.
| Model number | # of parameters | Δ BIC | |||
|---|---|---|---|---|---|
| Laboratory dataset N=44 | Online dataset 1 N=747 | Online dataset 2 N=235 | Clinical dataset N=116 | ||
| Model #1 | 6 | 0 | 0 | 0 | 0 |
| Model #2 | 5 | 43.99 | 918.62 | 228.00 | 41.57 |
| Model #3 | 8 | 122.13 | 1815.15 | 586.31 | 549.73 |
| Model #4 | 7 | 155.22 | 2662.46 | 814.28 | 424.10 |
| Model #5 | 7 | 82.24 | 709.23 | 204.32 | 151.56 |
| Model #6 | 6 | 47.76 | 2514.71 | 743.17 | 601.16 |
-
Δ BIC, Bayesian information criterion relative to the winning model (Model #1).
Objective vs. subjective happiness models.
| # of parameters | Δ BIC | ||||
|---|---|---|---|---|---|
| Healthy dataset N=1026 | Clinical dataset N=116 | Vanhasbroeck et al., dataset N=49 | Rutledge’s dataset N=46,204 | ||
| Model #1 | 6 | 0 | 0 | 0 | 0 |
| Model #6 | 5 | 3305.64 | 601.16 | 118.00 | 4997.29 |
| Model #1 vs 6 EP | - | 1.00 | 1.00 | 1.00 | 1.00 |
-
Δ BIC, Bayesian information criterion relative to the winning model (Model #1); EP, exceedance probability; EP greater than 0.95 was considered significant.
Time effect on happiness.
| Model number | # of parameters | Δ BIC | |||
|---|---|---|---|---|---|
| Healthy dataset N=1026 | Clinical dataset N=116 | Vanhasbroeck et al., dataset N=49 | Rutledge’s dataset N=46,204 | ||
| Model #1 | 6 | 0 | 0 | 0 | 0 |
| Model #7 | 6 | 1883.39 | 84.88 | 98.15 | 4152.61 |
| Model #8 | 6 | 1487.79 | 66.37 | 61.73 | 2295.40 |
| Model #9 | 6 | 1064.22 | 59.53 | 74.40 | 1458.34 |
-
Δ BIC, Bayesian information criterion relative to the winning model (Model #1).
Contributions to model fit.
| Model number | # of parameters | Mean R2 | |||
|---|---|---|---|---|---|
| Healthy dataset N=1026 | Clinical dataset N=116 | Vanhasbroeck et al., dataset N=49 | Rutledge’s dataset N=46,204 | ||
| M1: CR +EV + RPE | 6 | 0.68 | 0.47 | 0.58 | 0.69 |
| M10: No CR | 5 | 0.62 | 0.42 | 0.52 | 0.60 |
| M11: No EV | 5 | 0.59 | 0.40 | 0.54 | 0.60 |
| M12: No RPE | 5 | 0.42 | 0.29 | 0.34 | 0.50 |
Interaction between different wins and risk preference for gain on z-scored happiness.
| Datasets | Trial types | |
|---|---|---|
| Certain trials | Gamble trials (better outcomes; nonzero) | |
| Healthy dataset N=1026 | beta = –0.00; t=–0.07; p=0.943 95% CI = [–0.01, 0.01] | beta = –0.00; t=–0.24; p=0.807 95% CI = [–0.00, 0.00] |
| Clinical dataset N=116 | beta = –0.01; t=–0.64; p=0.524 95% CI = [–0.04, 0.02] | beta = –0.00; t=–2.13; p=0.034 95% CI = [-0.01,–0.00] |
| Vanhasbroeck et al., dataset N=49 | beta = 1.41; t=0.51; p=0.609 95% CI = [–4.00, 6.82] | beta = 0.38; t=0.55; p=0.581 95% CI = [–0.96, 1.71] |
| Rutledge’s dataset N=46,204 | beta = 0.10; t=0.84; p=0.402 95% CI = [–0.14, 0.34] | beta = 0.03; t=0.82; p=0.412 95% CI = [–0.04, 0.10] |
Correlations of the depression-specific factor with mood variation and βRPE.
| Laboratory dataset | Online dataset 1 | Online dataset 2 | Clinical dataset | |
|---|---|---|---|---|
| Mood variation | r=–0.309 p=0.041 | r=–0.191 p<0.001 | r=–0.135 p=0.039 | r=–0.239 p=0.009 |
| βRPE | r=–0.352 p=0.019 | r=–0.142 p<0.001 | r=–0.189 p=0.004 | r=–0.216 p=0.020 |
| Mediation | a×b = –0.306 95% CI: [-0.547,–0.065] p=0.015 | a×b = –0.096 95% CI: [-0.145,–0.047] p<0.001 | a×b = –0.126 95% CI: [-0.212,–0.040] p=0.004 | a×b = –0.141 95% CI: [-0.261,–0.038] p=0.021 |
Clinical characteristics of patients with affective disorders.
| Clinical variables | N=116 |
|---|---|
| Diagnosis (MDD/AD/MDD and AD/others) | 51/22/32/11 |
| Illness duration (months) | 19.57±19.22 |
| Medications (yes) | |
| SSRI | 96 |
| Antipsychotics | 57 |
| BZDs | 38 |
| Mood stabilizer | 14 |
-
Note that mood stabilizer refers to Lithium in this dataset.
Choice model comparisons.
| Model number | # of parameters | Δ BIC | |||
|---|---|---|---|---|---|
| Laboratory dataset N=44 | Online dataset 1 N=747 | Online dataset 2 N=235 | Clinical dataset N=116 | ||
| cM1 | 1 | 320.92 | 10959.49 | 3818.36 | 2368.27 |
| cM2 | 4 | 62.70 | 3655.22 | 1030.07 | 620.85 |
| cM3 | 6 | 0 | 0 | 0 | 0 |
-
Δ BIC, Bayesian information criterion relative to the winning model (cM3).