Development and cross-validation of a structural brain signature predicting delay discounting in healthy adults (Study 1)
Individual differences in impulsivity
On average, participants had a fitted log(k) parameter of -5.94 (median log(k)=-5.49, corresponding to k = 0.0041). Discounting rates were characterized by substantial individual differences (SD = 2.00), with log(k) ranging from − 11.92 to -2.16. These individual differences were very stable over a 7-week period as reported previously 14. On the I-8 subscale of urgency trait, participants’ average scores varied between 1 and 5 (mean = 2.72; median = 2.5; SD = 0.84). Log(k) showed a trend for a weak positive correlation with the urgency trait (R = 0.17, p = 0.06, 95%-CI= [-0.009, 0.35]).
Cross-validated predictions of delay discounting - Validity of predicted log(k) in healthy participants
The 10-fold cross-validation procedure revealed a significant accuracy of the brain-based prediction (see Fig. 2A and 2B and Supplementary Fig. 1): the predictions had a mean squared error of 3.45 (permutation test: p = 0.0026), a root mean squared error of 1.86 (permutation test: p = 0.0026), a mean absolute error for predicted log(k) of 1.46 (permutation test: p = 0.0022), and a cross-validated prediction-outcome correlation of R = 0.35 (permutation test: p = 0.0028) (Fig. 2C).
Further, supporting the reliability and conceptual validity of the brain-predicted log(k)’s, we found that brain-based predictions at baseline significantly correlated with (out-of-sample) log(k)’s computed from the ITC task performed seven weeks later (R = 0.34, p < 0.001, 95%-CI= [0.18, 1]) (Fig. 2D). This suggests that a relatively stable part of the between-person variability in delay discounting was explained by individual differences in brain structure. Moreover, higher brain-predicted log(k) values were associated with higher self-reported urgency (R = 0.20, p = 0.037, 95%-CI= [0.01, 0.37]) (Fig. 2E).
Like the actual measures of log(k) (see 14), brain-based predictions of log(k) did not significantly correlate with age (R=-0.11, p = 0.24, 95%-CI= [-0.29, 0.07]), education (R=-0.15, p = 0.10, 95%-CI= [-0.33, 0.03]), income (R=-0.12, p = 0.21, 95%-CI= [-0.30, 0.07]), BMI (R= -0.04, p = 0.66, 95%-CI= [-0.22, 0.14]), and percentage of body fat (R= -0.13, p = 0.18, 95%-CI= [-0.31, 0.06]) (see more details in Supplementary Fig. 2).
Performance of the Structural Impulsivity Signature in a second independent sample of healthy participants (Study 2)
Study 2 tests the predictions of the Structural Impulsivity Signature (SIS) in a second MRI dataset of healthy participants, that has used a different protocol, scanner, different preprocessing pipeline, in a socio-demographically different participant population.
Individual differences in impulsivity
The mean log(k) parameter in Study 2 was − 4.09 (median log(k)=-3.94, corresponding to a k of 0.019). Individual differences in the discounting parameter were less variable (SD = 0.98) as compared to Study 1, with log(k) ranging from − 7.08 to -2.12. Participants had average urgency trait scores (means of positive and negative urgency) varying between 1.00 and 3.01 (mean = 1.76; median = 1.68; SD = 0.48). In Study 2, log(k) had a trend for a negative correlation with urgency (R=-0.14, p = 0.06, 95%-CI= [-0.29, 0.008]). Therefore, in Study 2, the discounting rate does not seem to be related to individual differences in impulsivity.
Brain-based predictions of impulsivity - Validity of predicted log(k) in a second independent sample of healthy participants
For each participant in Study 2, we calculated the predicted individual log(k) as the dot-product between the weight map developed in Study 1 and the individual GMD map. We then tested whether predicted log(k) correlated with observed individual log(k) and with the impulsivity trait of urgency (UPPS subscales). While we did not find a significant link between predicted and observed log(k) in Study 2 (R = 0.06, p = 0.21, 95%-CI= [-0.07, 1]), predicted log(k) was positively associated with urgency scales (R = 0.15, p = 0.047, 95%-CI= [0.002, 0.30], see Fig. 2F), as in Study 1. Thus, the results of Study 2 partially validate the developed structural brain signature as a brain signature of impulsivity.
Validation of the structural brain signature in a clinical sample of bvFTD patients and matched controls (Study 3)
Our last analysis step aimed at further testing the generalizability of the SIS by evaluating its validity in a patient population that is characterized by impulsivity. Study 3 employed a distinct protocol from Studies 1 and 2 (different ITC task, different MRI scanner and parameters), and in a different, older population including dementia patients with substantial structural atrophy. This further allowed us to investigate the clinical relevance of the SIS (1) for classifying patients with bvFTD differently from matched control participants and (2) for predicting the core symptoms of disinhibition and executive deficits in patients with bvFTD 8 .
Differences of impulsivity between bvFTD patients and healthy controls
In line with the core symptoms of this disorder, bvFTD patients presented significantly higher delay discounting (i.e. more impatient or impulsive choices) compared to controls, for both money rewards and food rewards (see 35). They also showed higher inhibition deficit (Hayling-error score; t = 5.71, p < 0.001, Cohen’s d = 1.60, 95%-CI=[0.87, 2.33]) and lower executive performances (FAB score; t=-7.31, p < 0.001, Cohen’s d=-2.00, 95%-CI=[-1.23, -2.77]) compared to controls (see Supplementary table 1).
Brain-based predictions of impulsivity – Validity of predicted log(k) in bvFTD patients
To investigate the predictive validity of our classifier in Study 3, we first tested whether predicted log(k)’s (obtained from the brain pattern applied to each participant’s grey matter density map) were correlated with actual log(k)’s calculated in this study across the whole sample (patients and controls). This analysis showed that the predicted log(k) values were positively correlated with actual log(k) values, for both monetary rewards (R = 0.30, p = 0.03, 95%-CI= [0.03, 1], mean absolute error of 2.08) and for food rewards (R = 0.40, p = 0.006, 95%-CI= [0.15, 1], mean absolute error of 2.65) (see Fig. 3.A and 3.B).
We next tested whether the SIS predictions could distinguish bvFTD patients from controls. As expected, we found that brain-predicted log(k) was significantly higher in bvFTD patients than in controls (t = 3.60, p = 0.0009, Cohen’s d = 1.09, 95%-CI=[0.41, 1.76] – see Fig. 3.C). Notably, brain-predicted log(k) significantly predicted whether a grey matter density map was from a bvFTD patient or from a control participant, with a classification accuracy of 81% (p = 0.002, sensitivity = 87.5%, specificity = 72.2%, - see Fig. 3.D). Interestingly, the actual log(k)’s calculated for monetary and food rewards in this sample revealed slightly lower predictive accuracies and especially lower specificities: 73.7% accuracy for monetary rewards (p = 0.07, sensitivity = 100%, specificity = 37.5%,) and 76.3% for food rewards (p = 0.01, sensitivity = 100%, specificity = 47.1%).
We next investigated the relationship between brain-predicted log(k) and clinical measures of bvFTD core symptoms of disinhibition and executive deficits. Across both the patient and control groups, higher predicted log(k) was associated with higher inhibition deficit (higher Hayling-error score; R = 0.55, p = 0.0002, 95%-CI= [0.30, 0.74]) and higher executive troubles (lower FAB score; R=-0.56, p = 0.0001, 95%-CI= [-0.74, -0.30]). More interestingly, even within the group of bvFTD patients, higher predicted log(k) was associated with higher inhibition deficit (higher Hayling-error score; R = 0.52, p = 0.01, 95%-CI= [0.14, 0.77]) and higher executive troubles (lower FAB score; R=-0.43, p = 0.04, 95%-CI= [-0.71, -0.03]) (see Fig. 3.E and 3.F). Further, we checked that predicted log(k) was still significantly related to lack of inhibition (i.e., higher Hayling-error scores; B = 8.63, p = 0.02, 95%-CI= [1.51, 15.7]) within bvFTD patients even after controlling for executive function deficit; this added result showed that the relationship between brain-based predictions and disinhibition symptom was not only due to shared variance with the severity of dysexecutive syndrome. Together, these findings show that the SIS significantly and accurately classified bvFTD patients from matched controls, and that it tracked the severity of key symptoms in these patients.
Spatial distribution of weights in the structural brain signature (Study 1)
Thresholded pattern of structural brain signature
Bootstrapping results revealed the positive and negative weights that most strongly contributed to GMD-based prediction of individual differences in delay discounting. At a threshold of q = 0.05 FDR-corrected, we found two clusters in which grey matter density positively contributed to discounting differences (which means that higher grey matter density was associated with higher impatience); these clusters were in the left lateral parietal cortex (supramarginal gyrus) and left lateral occipital cortex (superior division). At a threshold of p = 0.001 uncorrected, we found additional clusters contributing positive weights, especially in regions of the valuation system 50 such as the right orbitofrontal cortex (OFC), ventromedial prefrontal cortex (vmPFC) and right ventral striatum.
At q = 0.05 FDR-corrected, there was one cluster in the posterior cingulate cortex (PCC) and adjacent lingual gyrus (including retrosplenial cortex) in which grey matter density contributed negatively to discounting differences (i.e., in which lower grey matter density was associated with higher impatience). At a threshold of p = 0.001 uncorrected, other important regions contributing negative weights were found in the left hippocampus, the right anterior insulae (AI), dorsal anterior cingulate cortex (ACC), and amygdalae. For display purposes, the bootstrapped weight map is displayed in Fig. 3A at a more comprehensive threshold (p = 0.05 uncorrected, see also Supplementary table 2).
Similarity of structural brain signature to meta-analytic maps
When comparing the predictive map of log(k) with meta-analytic uniformity maps 47, we observed that the highest similarities (spatial correlation r’s > 0.1 in absolute value) were with the “Emotions”, “Affect”, “Conflict” and “Imagery” meta-analytic maps (Fig. 4B). These spatial correlations were all negative, indicating that greater grey matter density in areas related to emotions, affect, conflict processing, and imagery contributes to predicting lower delay discounting or more ‘patient’ decision-making (or conversely, lower grey matter density in these areas predicts higher discounting and more impulsive decision-making). The “Emotions”, “Affect”, “Conflict” and “Imagery” meta-analytic maps correspond to overlapping functional networks (see Fig. 4.B). Among the most overlapping regions between these four networks (in red), the AI and dorsal ACC, corresponding to robust negative weights in the brain pattern, are known to be major hubs of the salience network 51.
Spatial distribution of brain regions contributing to higher predicted log(k) in bvFTD (Study 3)
To identify the main brain regions which contributed to differentiate bvFTD patients from controls on the brain-predicted log(k), we contrasted bvFTD patients versus controls in terms of voxel-wise pattern expression of the predictive map of log(k). To this end, for each bvFTD patient and each control participant, we computed an ‘importance map’ as the unsummed matrix dot product between the predictive structural weight map and the individual grey matter density map. Since higher resulting dot product contributes to higher predicted discounting, the importance map shows which brain regions contributed to increase (or decrease) predicted discounting in each individual. We performed a t-test contrasting bvFTD patients and controls (bvFTD > controls) on the resulting importance maps, with a family-wise error (FWE) correction applied to p-values to correct for multiple comparisons across the brain (see Fig. 5.C). This contrast shows the regions in which structural atrophy contributed positively to higher predicted discounting in bvFTD than in controls (regions in red). These included the OFC, anterior insulae, dorsal ACC, striatum, thalamus, amygdalae, hippocampus, and middle temporal regions. These regions corresponded to areas combining the presence of negative weights in the predictive brain pattern (i.e., voxels for which higher GMD predicts lower discounting and more patient decision-making, shown in Fig. 5.B) and the presence of significant grey matter atrophy due to bvFTD pathology (see atrophy pattern in Fig. 5.A). Thus, the contrast shown in Fig. 5.C also maps the regions in which the SIS is the most similar to bvFTD atrophy pattern.