Introduction
Major depressive disorder (MDD) is a severe mental health problem affecting millions worldwide. Despite a variety of treatment modalities, depression is prone to become a difficult-to-treat and chronic psychiatric illness, often referred to as treatment-resistant depression (TRD) (Akil et al., Reference Akil, Gordon, Hen, Javitch, Mayberg, McEwen and Nestler2018). Although still an issue of debate, personality features may not only add to the vulnerability to develop depression (Newton-Howes, Tyrer, & Johnson, Reference Newton-Howes, Tyrer and Johnson2006), it may drive treatment outcome and clinical non-response (Cloninger, Svrakic, & Przybeck, Reference Cloninger, Svrakic and Przybeck2006). Depressed patients with a comorbid personality disorder may respond less well to treatment, regardless of the chosen intervention (Rosenbluth, Macqueen, McIntyre, Beaulieu, & Schaffer, Reference Rosenbluth, Macqueen, McIntyre, Beaulieu and Schaffer2012). Furthermore, a long history of depression, with multiple episodes of resistance to treatment may induce personality changes within patients, which can be stable over time (Teraishi et al., Reference Teraishi, Hori, Sasayama, Matsuo, Ogawa, Ishida and Kunugi2015; Zaninotto et al., Reference Zaninotto, Solmi, Toffanin, Veronese, Cloninger and Correll2016), and this might also influence the results concerning the role of personality traits in relation to clinical outcome.
In pharmacotherapy research, personality dimensions have been evaluated to examine their ability to guide or even predict clinical outcome in MDD (Kaneda, Yasui-Furukori, Nakagami, Sato, & Kaneko, Reference Kaneda, Yasui-Furukori, Nakagami, Sato and Kaneko2011; Mulder, Reference Mulder2002; Mulder, Joyce, & Luty, Reference Mulder, Joyce and Luty2003). For instance, concerning the Temperament and Character Inventory (TCI; Cloninger, Przybeck, Svrakic, & Wetzel, Reference Cloninger, Przybeck, Svrakic and Wetzel1994), higher scores on the temperament dimension Harm Avoidance (HA) has been related to an increased occurrence of mood and anxiety disorders (Hansenne et al., Reference Hansenne, Reggers, Pinto, Kjiri, Ajamier and Ansseau1999; Hirano et al., Reference Hirano, Sato, Narita, Kusunoki, Ozaki, Kimura and Uehara2002; Nery et al., Reference Nery, Hatch, Nicoletti, Monkul, Najt, Matsuo and Soares2009; Spittlehouse et al., Reference Spittlehouse, Pearson, Luty, Mulder, Carter, McKenzie and Joyce2010). HA refers to the inhibition and cessation of behavior and is theoretically associated with the serotonergic system. Depressed patients displaying higher scores on Novelty Seeking (NS) – associated with the dopaminergic system and implying the activation and initiation of behavior – have been, however not consistently, linked with a more favorable antidepressant (AD) treatment outcome (Mulder et al., Reference Mulder, Joyce and Luty2003). Lower scores on the temperament dimension Reward Dependence (RD) – thought to be related to noradrenergic functioning and described as the maintenance and continuation of ongoing behavior – could be more at risk for TRD (Takahashi et al., Reference Takahashi, Shirayama, Muneoka, Suzuki, Sato and Hashimoto2013). The temperament dimension Persistence, formerly part of the RD personality construct but is currently regarded as a separate temperament dimension, has not been related to any monoaminergic neurotransmitter (NT) system yet (Cloninger et al., Reference Cloninger, Przybeck, Svrakic and Wetzel1994). Importantly, the brainstem is home to a group of modulatory NT pathways, such as those arising from the raphe nuclei (serotonergic), ventral tegmental area (VTA) (dopaminergic), and locus coeruleus (LC) (noradrenergic), which form a modulatory network that coordinates interactions between other networks in the brain via the hypothalamus and thalamus providing further input in the frontal, cingulate, and other regions of the cerebral cortex (Venkatraman, Edlow, & Immordino-Yang, Reference Venkatraman, Edlow and Immordino-Yang2017).
In a recent large cohort study conducted by Balestri et al. (Reference Balestri, Porcelli, Souery, Kasper, Dikeos, Ferentinos and Serretti2019), these authors showed that MDD patients who did not remit after adequate trials of AD continued to display high scores on HA. They also found that clinical non-response was associated with high scores on HA and low scores on RD, but most importantly they showed that lower scores on RD were found to be related to TRD. Unfortunately, no neurobiological markers were included. 18Fluorodeoxyglucose positron emission tomography (18FDG PET) is a well-known neuroimaging method which has been shown to be very insightful to detect brain metabolic disturbances in the depressed state (Mayberg, Reference Mayberg2009; Sublette et al., Reference Sublette, Milak, Galfalvy, Oquendo, Malone and Mann2013). Several 18FDG PET studies have examined the relationship between the TCI and regional brain glucose metabolism in healthy subjects (Hakamata et al., Reference Hakamata, Iwase, Iwata, Kobayashi, Tamaki, Nishio and Inada2009; Park, Park, & Kim, Reference Park, Park and Kim2016; Youn et al., Reference Youn, Lyoo, Kim, Park, Ha, Lee and Kwon2002). However, these studies employed univariate analytical techniques to identify brain–personality relationships that may mask within-group heterogeneity. Critically important for translational neuroscience, machine learning-based approaches may be better placed to predict personality traits using neuroimaging signatures at the individual subject level (Jiang et al., Reference Jiang, Calhoun, Zuo, Lin, Li, Fan and Sui2018).
Because it is largely unclear how regional glucose metabolism (CMRglc) may be implicated in the treatment-resistant depressed state and no studies have examined yet its influence on the prediction of temperament dimensions, we collected a group of well-defined TRD patients and we expected that metabolic activity in particular parts of the brain stem (and their respective NT systems) would significantly modulate the different temperament scores in TRD. We aimed to find CMRglc-based biomarkers for performing personalized prediction of temperament in TRD.
Methods
Participants
The study was approved by the ethics committee (UZBrussel) and all participants gave written informed consent. Forty AD-free depressed patients were included by using the Dutch version of the Mini-International Neuropsychiatric Interview (MINI; Sheehan et al., Reference Sheehan, Lecrubier, Sheehan, Amorim, Janavs, Weiller and Dunbar1998) excluding bipolar and psychotic depression. Subjects with substance abuse/dependence were also not included. To evaluate depression severity, patients were assessed with the 17-item Hamilton Depression Rating Scale (HDRS; Hamilton, Reference Hamilton1967). As described by Rush, Thase, & Dube (Reference Rush, Thase and Dube2003), all included MDD patients were considered at least stage III treatment-resistant: they had had a minimum of two unsuccessful treatment trials with a serotonin reuptake inhibitor/noradrenaline and/or serotonin reuptake inhibitor and one failed clinical trial with a tricyclic antidepressant. To avoid influences of concomitant AD medication, all patients went through a medication washout and they were AD-free for at least 2 weeks before the 18FDG PET scan. Right handedness was assessed with the Van Strien questionnaire (Van Strien & Van Beek, Reference Van Strien and Van Beek2000). This study was part of a larger project investigating several neurobiological and neurocognitive markers in TRD.
Assessment
Before performing the 18FDG PET scan, all participants were assessed with the HDRS (Hamilton, Reference Hamilton1967) and with the Dutch version of the TCI (de la Rie, Duijsens, & Cloninger, Reference de la Rie, Duijsens and Cloninger1998). The TCI contains three character dimensions: Cooperativeness, Self-Directedness, and Self-Transcendence, which are the complex result of learning processes, emotional maturation during one's entire life and comprise socio-cultural factors, and four temperament dimensions – which are thought to be largely hereditary, manifesting early in life and remaining largely stable throughout a lifetime – barely influenced by external factors (Cloninger, Reference Cloninger1987; Cloninger, Svrakic, & Przybeck, Reference Cloninger, Svrakic and Przybeck1993). Because only HA, NS, and RD are thought to be related to specific neurobiological (NTs) systems (Cloninger et al., Reference Cloninger, Przybeck, Svrakic and Wetzel1994), we specifically focused on these three temperament dimensions. Given that gender differences in MDD are well documented and a proximally 17.1% of women and 9.4% of men worldwide will experience a depressive episode during their lifetime (Alonso, Lepine, & Committee, Reference Alonso and Lepine2007), and given that differences have been reported on temperament dimensions, in the depressed as well as in the non-depressed state (Miettunen, Veijola, Lauronen, Kantojarvi, & Joukamaa, Reference Miettunen, Veijola, Lauronen, Kantojarvi and Joukamaa2007), we incorporated gender (but also age) in all our analyses.
All behavioral data were analyzed with SPSS 25 (Statistical Package for the Social Sciences, IBM, Chicago, USA). The significance level was set at p ⩽ 0.05, two-tailed, for all analyses (see Table 1).
Table 1. Demographics
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20221123115408919-0671:S0033291720005425:S0033291720005425_tab1.png?pub-status=live)
MDD, major depressive disorder; HDRS, Hamilton Depression Rating Scale; NS, novelty seeking; HA, harm avoidance; RD, reward dependence; P, persistence; SD, self-directedness; CO, cooperativeness; ST, self-transcendence.
Ratio's, means and standard deviations are provided.
a Temperament dimensions focused on for this study.
18FDG-PET brain imaging
Aligned with our previous CMRglc research (Baeken et al., Reference Baeken, De Raedt, Van Hove, Clerinx, De Mey and Bossuyt2009, Reference Baeken, Marinazzo, Everaert, Wu, Van Hove, Audenaert and De Raedt2015; Baeken, Wu, & De Raedt, Reference Baeken, Wu and De Raedt2018; van Heeringen, Wu, Vervaet, Vanderhasselt, & Baeken, Reference van Heeringen, Wu, Vervaet, Vanderhasselt and Baeken2017), static 18FDG-PET scan imaging was performed by a Siemens Ecat Accel scanner 30 min after intravenous tracer administration (222 MBq 18FDG). All 18FDG-PET scans were performed between 9:00 and 11:00 h and all participants had to lie supine with their eyes closed. Timings of producing 18FDG, injecting, and scanning were in accordance with the European Association of Nuclear Medicine Neuroimaging guidelines (Varrone et al., Reference Varrone, Asenbaum, Vander Borght, Booij, Nobili and Nagren2009). Emission data were obtained in 3D mode over 10 min. For transmission, germanium-68 sources (3 × 185 MBq; decay corrected) were used and data were acquired in 2D mode over 3 min. The emission data were reconstructed iteratively (OSEM 10 iterations, 32 subsets) applying a post-reconstruction filter (Gaussian filter 6 mm full-width-at-half-maximum). For the transmission, we used filtered back reconstruction and all data were subsequently segmented into regions with similar attenuation factors. To obtain attenuation correction factors for each line of response, these segmented images were forward projected. The pixel size was 2.57 mm transaxial and 3.38 mm axial. Each scan was projected onto a normalized brain template with the SPM12 (Statistical Parametric Mapping, Welcome Department of Imaging Neuroscience, University College of London, UK) software. The voxel-wise metabolic activity in the brain stem was extracted for individual prediction using the Desikan–Killiany Atlas (Desikan et al., Reference Desikan, Segonne, Fischl, Quinn, Dickerson, Blacker and Killiany2006).
Individualized prediction
We employed a sparse multiple kernel learning (MKL) with elastic-net regularization to select the minimal set of voxels that best predicted temperament scores, implementing in the SpicyMKL toolbox (version 3; Suzuki & Tomioka, Reference Suzuki and Tomioka2011). A nested cross-validation (CV) was used to perform hyperparameter optimization (10-fold CV for the inner loop, leave-one-out for the outer loop). Age, gender, and mean global CMRglc value were combined in one single linear kernel. For each voxel in the brainstem, a linear kernel was computed. A voxel with higher kernel weight (unsigned) can be interpretative as a higher contribution to the predictive model. The Pearson correlation coefficient (r) and root mean squared error (rMSE) were calculated between the actual and the predicted temperament scores for the overall predictive performance. Permutation tests were carried out to assess the statistical significance of the correlation coefficient and rMSE (randomly shuffle temperament scores 1000 times).
Results
For one 56-year-old male patient, the 18FDG-PET scan showed severe artifacts and was not suitable for analysis. Consequently, this patient was removed from the analysis.
18FDG-PET brain imaging results
As shown in Fig. 1, the predicted RD scores significantly correlated with the actual RD scores (r = 0.425, p = 0.004; rMSE = 3.418, p = 0.015). Contributing to the predictions of RD, the MKL regression identified metabolic activity in the pontine nuclei (PO) and LC, the medial reticular formation (MRF), the dorsal/median raphe nuclei, the VTA, red nucleus, and substantia nigra [Fig. 2, only those voxels that survived at CV procedure are reported, the spatial location is identified by AAL3 (Rolls, Huang, Lin, Feng, & Joliot, Reference Rolls, Huang, Lin, Feng and Joliot2020) and the Harvard Ascending Arousal Network Atlas (Edlow et al., Reference Edlow, Takahashi, Wu, Benner, Dai, Bu and Folkerth2012)]. In addition, the age, gender, and mean global CMRglc value also contribute to MKL prediction (4.64% of total weight). The voxel-wise correlations indicated that the metabolic activity in these regions was positively correlated with RD scores. In addition, brain stem CMRglc was not predictive for HA (r = 0.164, p = 0.115; rMSE = 8.34, p = 0.152) and NS (r = −0.122, p = 0.516; rMSE = 6.26, p = 0.533) scores.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20221123115408919-0671:S0033291720005425:S0033291720005425_fig1.png?pub-status=live)
Fig. 1. Scatter plot showing actual and predicted reward dependence (RD) scores (r = 0.425, p = 0.004; rMSE = 3.418, p = 0.015). The size of scatter point is proportional to the mean global CMRglc value.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20221123115408919-0671:S0033291720005425:S0033291720005425_fig2.png?pub-status=live)
Fig. 2. Three-dimensional overview of voxel-wise kernel weight in the brainstem for reward dependence (RD) prediction, with cluster size >10 voxels. The color represents the relative informativeness of different voxels. PO, pontine nuclei; LC, locus coeruleus; MRF, medial reticular formation; DRN, dorsal raphe nuclei; MRN, median raphe nuclei; VTA, the ventral tegmental area.
Discussion
This is the first study applying the MKL algorithm to predict the individual differences in temperament scores based on brainstem metabolic activity in TRD. The MKL model was able to predict the individual scores on RD from the brainstem metabolic activity, but brainstem CMRglc was not predictive for the individual scores on HA and NS. More precisely, we found that individual differences in RD are closely linked to metabolic increases in the raphe nuclei, the PO/LC, and the VTA, respectively, suggesting the implication of the serotonergic, noradrenergic, and dopaminergic system. These findings provide additional support that, at least in our case in the treatment-resistant depressed state, this temperament dimension shows significant increased CMRglc in brainstem regions associated with different monoaminergic NT systems (Samochowiec et al., Reference Samochowiec, Rybakowski, Czerski, Zakrzewska, Stepien, Pelka-Wysiecka and Hauser2001). Indeed, research not always has found the association between temperaments and their respective theoretically NT systems, and it remains an open question whether specific individual traits can be attributed to only one monoaminergic system (Robbins, Reference Robbins2018).
According to Cloningers' theory of personality, the temperament dimension RD has been linked with noradrenergic activity (Cloninger et al., Reference Cloninger, Przybeck, Svrakic and Wetzel1994). Norepinephrine, in addition to its effects on sensory processing and waking behavior, it contributes to cognition, flexibility, and mnemonic processes, and dense noradrenergic innervations in the prefrontal cortical and hippocampal areas are thought to modulate these functions (Borodovitsyna, Flamini, & Chandler, Reference Borodovitsyna, Flamini and Chandler2017). The personality temperament RD has been related to the formation of conditioned signals of reward (Gerra et al., Reference Gerra, Zaimovic, Timpano, Zambelli, Delsignore and Brambilla2000) and may be implicated to maintain the behavior previously associated with reward (Bond, Reference Bond2001). Subjects with lower scores on RD are described as less sentimental and less socially attached; they generally are less dependent on the approval of others, and those with higher scores than average are clinically characterized as eager to help and please others (Hansenne et al., Reference Hansenne, Reggers, Pinto, Kjiri, Ajamier and Ansseau1999). Given that RD also reflects the maintenance of socially rewarded behavior, our 18FDG-PET findings in the LC and the pontine noradrenergic nuclei suggests that in the treatment-resistant state, this could be mediated by the noradrenergic system, reflecting behavioral flexibility in social contacts. Besides that pontine microbleeds in stroke may lead to depression (Tang et al., Reference Tang, Chen, Liang, Chu, Mok, Ungvari and Wong2013, Reference Tang, Liu, Chen, Abrigo, Chu, Mok and Wong2014), the PO are the largest of the precerebellar nuclei, providing the principal input to the cerebellum (Kratochwil, Maheshwari, & Rijli, Reference Kratochwil, Maheshwari and Rijli2017; Schwarz & Thier, Reference Schwarz and Thier1999). Not surprisingly that current evidence also points to an important role of the cerebellum in the processing of social behaviors (Adamaszek et al., Reference Adamaszek, D'Agata, Ferrucci, Habas, Keulen, Kirkby and Verhoeven2017).
It is of interest to note that metabolic activity in serotonergic and dopaminergic-related brainstem structures was also predictive for RD scores when clinically depressed. Serotonin is considered as a fundamental mediator of emotional, motivational, and cognitive aspects of reward representation (Kranz, Kasper, & Lanzenberger, Reference Kranz, Kasper and Lanzenberger2010). The increased CMRglc in the raphe nuclei is suggestive for the involvement of the serotonergic system in this temperament construct in the desire for social relationships (Gerretsen et al., Reference Gerretsen, Graff-Guerrero, Menon, Pollock, Kapur, Vasdev and Mamo2010). It encodes reward probability, reward value and cost, and how beneficial the current environmental context may be (Luo, Li, & Zhong, Reference Luo, Li and Zhong2016). Besides that the dorsal raphe nucleus (DRN), the largest serotonergic nucleus in the brainstem, strongly moderates reinforcement learning, its neuronal activity has been associated with reward size (Nagai et al., Reference Nagai, Takayama, Nishitani, Andoh, Koda, Shirakawa and Kaneko2020) and reward waiting (Luo, Zhou, & Liu, Reference Luo, Zhou and Liu2015). Together with the medial raphe nuclei (MRN), the DRN modulates dopaminergic activity in the VTA, a region which has been primarily implicated in reward-motivated behavior, similar to the red nucleus and the substantia nigra (Browne et al., Reference Browne, Abela, Chu, Li, Ji, Lambe and Fletcher2019; Cha et al., Reference Cha, Carlson, Dedora, Greenberg, Proudfit and Mujica-Parodi2014). Of interest, it is the interaction between the serotonergic and the dopaminergic system that controls the behavioral response to reward (or punishment) (Macoveanu, Reference Macoveanu2014), potentially mediated by excitatory glutamatergic neurons (Wang et al., Reference Wang, Zhang, Qi, Wang, Cachope, Mejias-Aponte and Morales2019).
Although we only have focused on brainstem CMRglc, the positive relationship between RD scores and several specific brainstem areas associated with specific monoaminergic NT involvement suggest that, at least in the TRD state, the temperament RD construct goes beyond the noradrenergic system. This assumption is partly supported by the predictive metabolic activity of RD in the MRF, a medial brainstem network part of the emotional motor system mediated by multiple NT systems, also moderating emotional arousal and maintenance of consciousness (Venkatraman et al., Reference Venkatraman, Edlow and Immordino-Yang2017). Because we only measured CMRglc and no specific monoaminergic system, we can only assume that the higher metabolic activity predictive for higher scores on RD may represent higher order reward-related behavior and is mediated by noradrenaline, serotonin, and dopamine.
Our study had some major advantages: all TRD patients were AD-free on the time of scanning, and we implemented an MKL approach to identify CMRglc in brainstem regions predictive to determine temperament at the individual level. This overcomes the limitations of using univariate analysis which possesses limited sensitivity to identify subtle and spatially distributed effects, which could be important in complex personality dimensions with complex behavioral outcomes, involving multiple brain stem areas with a variety of functions (Fernandes et al., Reference Fernandes, Portugal, Alves, Arruda-Sanchez, Rao, Volchan and Mourao-Miranda2017; Jiang et al., Reference Jiang, Calhoun, Zuo, Lin, Li, Fan and Sui2018). MKL avoids the traditional in-sample correlation approach, prone to overfitting and not always allowing for out-of-sample generalizations, and not readily applicable to use on the individualized level. Notwithstanding, our CMRglc results cannot be generalized to the non-depressed state as it has been reported that RD scores may differ between resistant and non-resistant depressed samples (Balestri et al., Reference Balestri, Porcelli, Souery, Kasper, Dikeos, Ferentinos and Serretti2019). Of course, because we used 18FDG PET as a tracer and no specific ligand to quantify any of the potentially involved NT systems, our brainstem CMRglc findings cannot specifically be attributed to NT processes. Given the a priori hypothesis of brainstem involvement, we cannot make claims on CMRglc predictions related to temperaments in other parts of the brain. Also, given the relatively poor resolution of 18FDG PET and that each scan was projected onto a normalized brain template and not to the individual anatomical scan, all result interpretations should be done cautiously. Lastly, the TCI was merely used as a screening for personality temperaments and it was not intended for an actual diagnosis of an underlying personality disorder.
In conclusion, our MKL model showed that brainstem CMRglc was predictive only for the temperament dimension RD. The predictive value of brainstem metabolic activity on RD scores indicates that this temperament dimension could be mediated by different monoaminergic systems. Our findings expand our neurobiological insight that in the treatment-resistant depressed state, potentially other NT systems besides the noradrenergic system may contribute to with the expressed desire for social relationships (Gerretsen et al., Reference Gerretsen, Graff-Guerrero, Menon, Pollock, Kapur, Vasdev and Mamo2010). Machine learning paradigms may have good potential to further elaborate on the neurobiological mechanisms of personality. Future brain imaging studies using specific NT ligands are needed not only to investigate whether our findings on temperament also hold in the healthy and non-treatment-resistant depressed state.
Financial support
This research was supported by the National Natural Science Foundation of China (Grant No. 61876156), and a grant from the Scientific Fund W. Gepts UZBrussel. This work was also supported by the ‘Bijzonder Onderzoeksfonds’ [grant number BOF 16/GOA/017] and the ‘Rode Neuzen’ funding for scientific research [grant number G0F4617N].
Conflict of interest
None.
Ethical standards
The authors assert that all procedures contributing to this work comply with the ethical standards of the relevant national and institutional committees on human experimentation and with the Helsinki Declaration of 1975, as revised in 2008.