Hostname: page-component-78c5997874-j824f Total loading time: 0 Render date: 2024-11-13T07:03:49.990Z Has data issue: false hasContentIssue false

Genetic risk scores associated with temperament clusters in Finnish depression patients

Published online by Cambridge University Press:  04 September 2023

Simo-Pekko Salminen*
Affiliation:
Department of Psychiatry, Faculty of Medicine and Health Technology, Tampere University, Tampere, Finland
Anssi Solismaa
Affiliation:
Department of Psychiatry, Faculty of Medicine and Health Technology, Tampere University, Tampere, Finland Department of Psychiatry, Tampere University Hospital, Pirkanmaa Hospital District, Tampere, Finland
Leo-Pekka Lyytikäinen
Affiliation:
Department of Clinical Chemistry, Fimlab Laboratories and Finnish Cardiovascular Research Center, Faculty of Medicine and Health Technology, Tampere University, Tampere, Finland
Vesa Paavonen
Affiliation:
Department of Psychiatry, Faculty of Medicine and Health Technology, Tampere University, Tampere, Finland
Nina Mononen
Affiliation:
Department of Clinical Chemistry, Fimlab Laboratories and Finnish Cardiovascular Research Center, Faculty of Medicine and Health Technology, Tampere University, Tampere, Finland
Terho Lehtimäki
Affiliation:
Department of Clinical Chemistry, Fimlab Laboratories and Finnish Cardiovascular Research Center, Faculty of Medicine and Health Technology, Tampere University, Tampere, Finland
Esa Leinonen
Affiliation:
Department of Psychiatry, Faculty of Medicine and Health Technology, Tampere University, Tampere, Finland Department of Psychiatry, Tampere University Hospital, Pirkanmaa Hospital District, Tampere, Finland
Olli Kampman
Affiliation:
Department of Psychiatry, Faculty of Medicine and Health Technology, Tampere University, Tampere, Finland Department of Psychiatry, Tampere University Hospital, Pirkanmaa Hospital District, Tampere, Finland Department of Clinical Sciences (Psychiatry), Umeå University, and Västerbotten Welfare Region, Umeå, Sweden Department of Clinical Sciences (Psychiatry), University of Turku, Turku, Finland Wellbeing Services County of Ostrobothnia, Vaasa, Finland
*
Corresponding author: Simo-Pekko Salminen; Email: [email protected]
Rights & Permissions [Opens in a new window]

Abstract

Objective:

Cloninger’s temperament dimensions have been studied widely in relation to genetics. In this study, we examined Cloninger’s temperament dimensions grouped with cluster analyses and their association with single nucleotide polymorphisms (SNPs). This study included 212 genotyped Finnish patients from the Ostrobothnia Depression Study.

Methods:

The temperament clusters were analysed at baseline and at six weeks from the beginning of the depression intervention study. We selected depression-related catecholamine and serotonin genes based on a literature search, and 59 SNPs from ten different genes were analysed. The associations of single SNPs with temperament clusters were studied. Using the selected genes, genetic risk score (GRS) analyses were conducted considering appropriate confounding factors.

Results:

No single SNP had a significant association with the temperament clusters. Associations between GRSs and temperament clusters were observed in multivariate models that were significant after permutation analyses. Two SNPs from the DRD3 gene, two SNPs from the SLC6A2 gene, one SNP from the SLC6A4 gene, and one SNP from the HTR2A gene associated with the HHA/LRD/LP (high harm avoidance, low reward dependence, low persistence) cluster at baseline. Two SNPs from the HTR2A gene were associated with the HHA/LRD/LP cluster at six weeks. Two SNPs from the HTR2A gene and two SNPs from the COMT gene were associated with the HP (high persistence) cluster at six weeks.

Conclusion:

GRSs seem to associate with an individual’s temperament profile, which can be observed in the clusters used. Further research needs to be conducted on these types of clusters and their clinical applicability.

Type
Original Article
Copyright
© The Author(s), 2023. Published by Cambridge University Press on behalf of Scandinavian College of Neuropsychopharmacology

Significant outcomes

  • Genetic risk scores in a multivariate model associated with a clustered model of Cloninger’s temperament dimensions. These clusters have been shown to be associated with the comorbidity of anxiety disorders in depressed patients and treatment response.

Limitations

  • No predetermined widespread criteria for clustering multiple Cloninger’s temperament dimensions exist.

  • This study had a limited number of participants, which might have resulted in the inability to detect subtle clinical differences.

Introduction

In previous research, Cloninger has proposed seven dimensions, all with several subscales, which could be used to describe and study the temperament and character of an individual (Cloninger, Reference Cloninger1986; Cloninger et al., Reference Cloninger, Svrakic and Przybeck1993). He proposed that especially four of these dimensions, novelty seeking (NS; initiation or activation of appetitive behaviour in response to novelty), harm avoidance (HA; inhibition of behaviour in response to signs of punishment), reward dependence (RD; maintenance of behaviour in response to cues of social reward), and persistence (P; maintenance of behaviour despite frustration) are mainly genetically defined, and are linked to specific neurotransmitters: NS to the dopaminergic system, HA to the serotonergic system, and RD to the noradrenergic system. Cloninger also created a questionnaire, which could be used to grade an individual’s temperament and character traits on these dimensions (Cloninger et al., Reference Cloninger, Przybeck, Svrakic Dragan and Wetzel Richard1994). This questionnaire is called the Temperament and Character Inventory (TCI), which was later revised to create the Temperament and Character Inventory-Revised (TCI-R), which has 240 items, each with a scale of 5 possible answers. This questionnaire is the primary tool used when assigning an individual’s temperament and character traits to Cloninger’s model.

An association between Cloninger’s temperament and character traits and depression (Kampman & Poutanen, Reference Kampman and Poutanen2011; Zaninotto et al., Reference Zaninotto, Solmi, Toffanin, Veronese, Cloninger and Correll2016) and anxiety (Kampman et al., Reference Kampman, Viikki, Jarventausta and Leinonen2014) has been established in meta-analyses. Research has also been conducted about the relationship between the genes defining the catecholamine and serotonin systems and Cloninger’s temperament dimensions (Suzuki et al., Reference Suzuki, Matsumoto, Ishii, Oshino, Goto and Otani2007; Lee et al., Reference Lee, Yang, Lee, Kim, Joe, Jung, Choi and Ham2008; Alfimova et al., Reference Alfimova, Monakhov, Golimbet, Korovaitseva and Lyashenko2010; Tsuchimine et al., Reference Tsuchimine, Yasui-Furukori, Sasaki, Kaneda, Sugawara, Yoshida and Kaneko2012; Tsuchimine et al., Reference Tsuchimine, Taniguchi, Sugawara, Kaneda and Yasui-Furukori2013; Glavina Jelas et al., Reference Glavina Jelas, Devic and Karlovic2018). At this point in time, the results of these studies are inconclusive, and further studies need to be conducted.

Paavonen et al. (Reference Paavonen, Luoto, Lassila, Leinonen and Kampman2018) have proposed that instead of describing an individual’s risk for specific psychiatric disorders based on single TCI-R dimensions, it would be more sensible to create a cluster defined by multiple TCI dimensions to assess the risk of these conditions. Paavonen et al. proposed the following clusters: HNS, HHA/LRD/LP, LNS, and LHA/HP, where H stands for high level of a dimension, L stands for low level of a dimension, and non-specified dimensions being at an intermediate level. Based on their findings, these clusters have an association with specific forms of anxiety in depressed patients: the HNS cluster with panic disorder and agoraphobia, the HHA/LRD/LP cluster with social anxiety disorder (SAD), and the LNS cluster with generalised anxiety disorder. The LHA/HP cluster was associated with a lower risk of anxiety in general. In an earlier study, different temperament clusters have been associated with the severity of depression, and together with vegetative symptoms, predicted how well individuals responded to antidepressants (Paavonen et al., Reference Paavonen, Kampman, Illi, Viikki, Setälä-Soikkeli and Leinonen2014).

The aim of this study is to investigate the association between single nucleotide polymorphisms in the catecholamine and serotonin systems and the temperament clusters in Finnish population who were depressed. In this study, we investigated if specific single nucleotide polymorphisms (SNPs) related to the catecholamine and serotonin systems that had a connection with dimensions in the TCI-R questionnaire based on previous studies are associated with the temperament clusters proposed by Paavonen et al. (Reference Paavonen, Luoto, Lassila, Leinonen and Kampman2018).

Material and methods

Patients

The study included 242 patients from the Finnish region of Southern Ostrobothnia, who were referred to psychiatric specialised healthcare units for symptoms of depression, anxiety, suicidality, insomnia, and substance abuse during 2009–2013. The participants of the study were screened before inclusion in regards to alcohol use with Alcohol Use Disorders Identification test (AUDIT) (Bohn et al., Reference Bohn, Babor and Kranzler1995), and depressive symptoms using Beck Depression Inventory (BDI) version 1A (Beck et al., Reference Beck, Guth, Steer and Ball1997). The inclusion criterion was a score of at least 17 on the BDI scale, which is an indication of moderate depression. Patients whose symptoms were present due to organic conditions of the brain or had psychotic symptoms (ICD-10 F2* diagnoses) were excluded.

At baseline, the participants were assessed with the Mini International Neuropsychiatric Interview 5.0 (MINI) (Sheehan et al., Reference Sheehan, Lecrubier, Sheehan, Amorim, Janavs, Weiller, Hergueta, Baker and Dunbar1998) and the Montgomery-Åsberg Depression Rating Scale (MADRS) (Montgomery & Asberg, Reference Montgomery and Asberg1979), and the TCI-R. There were some dropouts after screening, a total of 228 patients (94% of all recruited patients) were examined with MADRS at baseline. A psychiatric assessment was also conducted on the participants at baseline, where their medication was assessed and optimised. The fluoxetine equivalent dose was calculated (Hansen et al., Reference Hansen, Moore, Dusetzina, Leinwand, Gartlehner and Gaynes2009) to make comparing the participants’ medication possible. At 6 weeks, the MADRS score and TCI-R questionnaire were repeated.

Demographic and clinical assessments were conducted at the beginning of the study. In this study, depressive disorders were considered as primary, anxiety disorders as secondary, and other diagnoses as tertiary, and all diagnoses were assigned based on clinical standards. For their main diagnoses, 181 (88.7%) patients had MDD, eight (3.9%) had dysthymic disorder, eleven (5.4%) had an anxiety disorder), three (1.5%) had self-destructiveness, and one (0.5%) had alcohol use disorder. 67% of patients with a mood disorder as their primary diagnosis had comorbid anxiety disorders. 65.7% of patients reported having at least one episode of MDD in their history. The patients were also divided into two groups based on their AUDIT questionnaire scores. Patients with a score of <11 were considered to have no alcohol use problems (non-AUP), and patients with a score of ≥11 were considered to have alcohol use problems (AUP). After some patient dropouts, these criteria created two groups: 127 non-AUP patients (mean age 38.5 ± 12.8 years), and 89 AUP patients (mean age 38.7 ± 11.4 years). The MINI was conducted on 82 AUP patients (data missing in 7 cases), and 61 (74.4%) of them were diagnosed with current alcohol use disorder based on the DSM-IV criteria, and eight (9.8%) of them with other current substance use disorders. The specific demographics of each temperament cluster are presented in Table 1.

Table 1. Demographics of each temperament cluster including age, sex, smoking status, MADRS and AUDIT scores at the beginning of the study and antidepressant doses in fluoxetine equivalents

  • MADRS score at the time of creating the corresponding cluster (either at the beginning of the study or 6 weeks).

  • AUDIT score at the beginning of the study.

Further specifics of the study protocol and patient details can be found in previous publications of the Ostrobothnia Depression Study (ClinicalTrials.gov identifier NCT02520271), Paavonen et al. (Reference Paavonen, Luoto, Koivukangas, Lassila, Leinonen and Kampman2016). The local Human Subjects Review Committee approved gathering the data used in this study, and it was carried out in accordance with the latest version of the Declaration of Helsinki. Since individuals with psychosis or organic brain diseases were excluded, all participants were able to give their informed written consent. For further information on the ethical aspects of this study, visit ClinicalTrials.org, identifier NCT02520271.

DNA extraction and genotyping

Blood samples for genotyping were gathered from a total of 212 participants. DNA was extracted from peripheral blood leukocytes using QIAamp DNA Blood Mini kit and an automated biorobot M48 extraction (versus, Hilden, Germany). Genotyping was performed using Illumina Infinium HumanCoreExome-12 DNA Analysis Beadchip version 1.0, according to the manufacturer’s recommendation at Helmholtz Zentrum, München, Germany. The following quality control filters were applied: GenCall score <0.15, GenTrain score <0.20, sample and SNP call rate <0.95, Hardy–Weinberg equation p-value <10−6, excess heterozygosity, cryptic relatedness (pi-hat >0.2), gender check and multidimensional scaling.

Literature search and selection of genes

A search of PubMed was carried out on the 18th of December 2018, to explore which SNPs had been studied in relation to Cloninger’s temperament dimensions. A search that included the keywords ‘temperament’, ‘harm avoidance’ or ‘temperament and character inventory’, and had a mention of ‘genetics’, as well as ‘5-HT’, ‘norepinephrine’, ‘noradrenaline’, or ‘dopamine’ were included in the search. With this search, 217 studies were included. All abstracts from these studies from 1997 onwards were read, from which specific SNPs and genes were recorded. SNPs and genes with only one research paper written on them were excluded. From this list, all the genes that were related to the dopaminergic, noradrenergic, or serotonergic synaptic cleft or the immediate metabolism of the neurotransmitters from the synaptic cleft were included in the next step. From these genes, and from a 5000 base pair distance backwards from the beginning of these genes, all SNPs included in the HumanCoreExome-12 beadchip were included in the study. Uncommon SNPs with a mean allele frequency (MAF) of <5% in the population of this study were excluded.

The genes that were included in the study were the catechol-O-methyltransferase (COMT), the dopamine receptor D2 (DRD2), the dopamine receptor D3 (DRD3), the dopamine receptor D4 (DRD4), the 5-hydroxytryptamine receptor 2A (HTR2A), the 5-hydroxytryptamine receptor 2B (HTR2B), the 5-hydroxytryptamine receptor 2C (HTR2C), the monoamine oxidase A, the solute carrier family 6 member 2 (SLC6A2), the solute carrier family 6 member 3 (SLC6A3), and the solute carrier family 6 member 4 (SLC6A4). The chip used did not have any SNPs with a MAF of over 5% in our sample in the DRD4 gene; thus, the DRD4 gene was excluded from the study. A total of 59 SNPs from the selected genes were available for the analyses from the used genotyping array. A combination of three models was created from the SNPs: the recessive, additive, and dominant model, as described in a study by Solismaa et al. (Reference Solismaa, Kampman, Lyytikäinen, Seppälä, Viikki, Mononen, Lehtimäki and Leinonen2017). Other types of gene polymorphisms, such as insertions and deletions, were not available with the genotyping methods used in this study.

Statistical methods

Confounding factors were controlled in the analyses. Sex, age, and smoking were controlled for all the analyses. The MADRS score at the beginning of the study was also considered for all the clusters created at baseline, since the TCI dimensions used to create the temperament clusters may be subject to confoundment from major depressive disorder (MDD) or the other psychiatric conditions the participants had. At six weeks after re-assessing the optimal medication, the TCI dimensions likely reflect the individual’s temperament and character with less confoundment as the baseline TCI dimensions (Hansenne & Bianchi, Reference Hansenne and Bianchi2009; Kampman & Poutanen, Reference Kampman and Poutanen2011). Thus, taking the MADRS score at six weeks into account was not done automatically for the analyses.

Other possible confounding factors that were explored were the fluoxetine equivalent dose of antidepressants, the MADRS score at the sixth week of the study, the AUDIT score at baseline, if the participant’s MDD is recurrent or not and BMI. The correlation of each temperament cluster and these possible confounding factors were analysed with the Welch two-sample t-test (continuous confounding factors) or Pearson’s chi-square test (dichotomic confounding factors). All the confounding factors with an association with the temperament clusters with a p-value of <0.05 were used as explanatory variables in the analyses. Different models with different explanatory factors were explored. In these explorations, confounding factors with a p-value of <0.05 were included in all analyses, while other possible factors were explored in different combinations. Before performing multivariate analyses, interactions between confounding factors were explored with Spearman’s rank-order correlation. Correlations of r > 0.5 were considered significant.

The SNPs were analysed in relation to the temperament clusters, which were based on the studies of Paavonen et al. Reference Paavonen, Kampman, Illi, Viikki, Setälä-Soikkeli and Leinonen(2014, Reference Paavonen, Luoto, Lassila, Leinonen and Kampman2018). In these studies, only clusters at baseline were used. In our study, we also included clusters which were formed at six weeks from the beginning of the study, created with similar methods as the clusters at baseline. At six weeks, the clustering algorithm formed slightly different clusters: HNS/HP (corresponding with HNS at baseline), HHA/LRD/LP (corresponding with HHA/LRD/LP at baseline), LNS/HHA/HRD (corresponding with LNS at baseline), and HP (corresponding with LHA/HP at baseline).

The clusters at both baseline and at six weeks were analysed, for a total of eight different cluster analyses. Eight different dichotomised variables were formed: the cluster of interest was compared against the other clusters at the same timepoint. The association of SNPs with the temperament clusters was analysed with logistic regression models. Appropriate confounding factors were used in the logistic regression analyses as explanatory factors. The p-values of this analysis were adjusted with the false discovery rate (FDR) method to counteract the multiple comparison problem when analysing a large number of SNPs with their different modes of inheritance. Genetic risk scores (GRS) were calculated for each temperament cluster separately. The final GRS consists of the weighted sum of risk alleles related to each cluster.

To form the GRSs, SNPs with a non-adjusted p-value of ≤0.05 were inserted into the same model as the confounding factors, and analysed with a stepwise logistic regression analysis, using Akaike information criteria (Akaike, Reference Akaike1981). This model was used to select SNPs that together best explain the response variable. This model also takes into account the problem of multiple SNPs having a high linkage disequilibrium (LD) by only including the single best SNP of these for the model. The score for each SNP was calculated by taking the beta estimate values obtained from the initial logistic regression analysis and using it to weight the individual SNPs in the GRS. The GRS was used in the final logistic regression analysis, which again included the confounding factors as explanatory variables.

When combining the genetic effects of individual SNPs that were not statistically significant after FDR adjustment, there is a risk of overfitting the data and causing type I errors. To avoid this, it is important to create an estimate of the null distribution of the GRS test statistics (Hu et al., Reference Hu, Muise, Xing, Brumell, Silverberg and Xu2013). Thus, a permutation test was conducted. Patient ID numbers were sampled from the phenotype data, mixing the genotype and phenotype data. This data was then analysed with identical methods as the original data. This sampling was repeated at least 999 times for each cluster and each combination of confounding factors, and up to 9999 times for each cluster with possible significant findings. The p-values from all these permutations were collected. If the p-value from the actual data ended up being among the smallest 5% of all the collected p-values, the results were deemed statistically significant. Lastly, LD analyses were conducted on all the SNPs in the study, with a window of 10,000 kbp and a threshold of r 2 = 0.8.

Plink version v1.90b5.4 (updated on the 10th of April 2018) and version 2.00a3 (updated on the 23rd of September 2020) were used to extract the SNPs from the gene data and to calculate the SNPs’ LD. Statistical analyses were performed with R (version 1.2.5042, R Core Team (2020). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. URL https://www.R-project.org/).

Results

Exploration and correlation of confounding factors in significant clusters

At baseline, patients in the HHA/LRD/LP cluster had a higher MADRS score (mean 24.5 ± 5.39 (HHA/LRD/LP) versus mean 22.6 ± 6.60 (rest of the sample), p = 0.033) and had higher antidepressant doses (fluoxetine equivalent doses: mean 32.4 ± 19.2 SD (HHA/LRD/LP) versus mean 25.8 ± 20.5 SD (rest of the sample), p = 0.031) compared to the rest of the sample. Patients in the HHA/LRD/LP cluster at baseline were also younger than in the rest of the sample (mean 34.8 ± 11.4 SD (HHA/LRD/LP) versus mean 40.0 ± 12.2 SD (rest of the sample), p = 0.005).

At six weeks, patients in the HP cluster had lower MADRS scores (mean 14.3 ± 7.75 SD (HP) versus mean 18.2 ± 7.81 SD (rest of the sample), p = 0.0027), and tended to be older (mean 41.5 ± 12.6 SD (HP) versus mean 37.6 ± 12.1 SD (rest of the sample), p = 0.055) than the rest of the sample. At six weeks, there were less smokers in the HP cluster when compared to the rest of the sample (40.8% (HP) versus 59.8% (rest of the sample), χ2 = 3.98, p = 0.046). At six weeks, the HHA/LRD/LP cluster had a slightly higher MADRS score than the rest of the sample (mean 19.0 ± 7.23 SD (HHA/LRD/LP) versus mean 16.2 ± 8.15 SD (rest of the sample), p = 0.034). The other analysed clusters also had their own combination of significant confounding factors, but as they had no significant results in SNP analyses, so the group differences aren’t reported here.

No two confounding factors had an association of 0.5 or higher with each other in Spearman’s rank-order correlation analyses, thus none of the factors were dropped from multivariate models based on this calculation. The HHA/LRD/LP cluster at six weeks was found to have both the MADRS score at baseline and at six weeks as significant confounding factors based on the exploration of these factors. To avoid bias, only the MADRS score at six weeks was chosen as an explanatory variable out of these two variables. The confounding factors considered in the clusters with significant results can be found in Table 2.

Table 2. Logistic regression models with different temperament clusters as dependent variables, explanatory variables used in the model and results of the permutation analysis

*= Additive coding, ** = Recessive coding, *** = Dominant coding.

Results of SNP analyses

No single SNP had a significant adjusted p-value in the models where appropriate confounding factors were included as explanatory variables. In the GRS analyses, significant associations were observed in the HHA/LRD/LP cluster at baseline, in the HP cluster at six weeks, and in the HHA/LRD/LP cluster at six weeks. In the model used, when explanatory factors were considered, six SNPs had an association with the HHA/LRD/LP cluster at baseline, four SNPs had an association with the HP cluster at six weeks, and two SNPs had an association with the HHA/LRD/LP cluster at six weeks. All significant findings are summed up in Table 2.

When exploring the combination of explanatory factors in the analyses, the combination of significant SNPs in the model remained constant in the HHA/LRD/LP and HP clusters at six weeks. In the HHA/LRD/LP cluster at baseline, taking out sex as an explanatory factor affected the combination of these SNPs. Compared to the analysis with sex as an explanatory variable, this latter set of SNPs did not include rs582385 in the HTR2A gene and did include rs582854, rs4942582 and rs7330636 in the HTR2A gene. We conducted a post hoc exploration about the linkage of these SNPs. The largest LD was between rs582385 and rs582854 in the HTR2A gene, with an LD of 0.28, and rs7330636 had no linkages in the analysis. In the other clusters analysed, we observed no significant results based on the p-values from the permutations conducted. In the LD analyses, rs6280 in the DRD3 gene was linked to rs324029 in the same gene (LD r 2 = 0.84), and rs223939 in the COMT gene was linked to rs4818 in the same gene (LD r 2 = 0.96). Detailed results about the SNPs and LDs can be found in Table 3.

Table 3. Linkage disequilibrium (LD) with other SNPs in the sample and annotations of SNPs in the GRSs found to be statistically significant after permutations

(a) A literature search on the significantly associated SNPs in this study can be found in the supplementary material.

(b) Significant SNPs in the analysis without sex as a confounding factor.

(c) Also as an Intron in gene HTR2A-AS1.

(d) Also nearGene-5 in gene MIR4761.

Discussion

The main result of this study is that significant associations were observed between the combination model of SNPs from the DRD3, SLC6A2, SLC6A4 and HTR2A genes and the HHA/LRD/LP cluster at baseline, the combination models of the SNPs from the HTR2A gene and the HHA/LRD/LP cluster at six weeks, and the combination models of the SNPs from the HTR2A and COMT genes and the HP clusters at six weeks. Based on this finding, it seems that polymorphisms in the DRD3, SL6A2, SLC6A3, SLC6A4, HTR2A, and COMT genes may alter the catecholaminergic and serotonergic transmission processes significantly, which influences an individual’s temperament. In the study by Paavonen et al. (Reference Paavonen, Luoto, Lassila, Leinonen and Kampman2018), the HHA/LRD/LP cluster at baseline was associated with SAD in depressed patients. In the same study, the LHA/HP cluster at baseline, which corresponds with the HP cluster at six weeks, was associated with the lowest prevalence of anxiety disorders in general. Paavonen et al. (Reference Paavonen, Kampman, Illi, Viikki, Setälä-Soikkeli and Leinonen2014) also observed these temperament clusters to predict treatment response in depression when including high vegetative symptoms in the analysis, though this effect was only observed in post-treatment. Based on these findings, the risk of anxiety disorders during depressive episodes associated with these temperament profiles could be mediated by certain genotype clusters. By genotyping the SNPs associated with these clusters in depressed patients, we could estimate the risk of specific anxiety disorders in these patients. Further research needs to be conducted to repeat these results and estimates the clinical importance of these types of findings before widespread genotyping could be considered in MDD patients. However, having a genetic model to help with diagnostics and treatment would be beneficial in clinical practice.

The results of the HHA/LRD/LP cluster at baseline changed depending on whether sex was controlled as a confounding factor. In the initial exploration of confounding factors, sex did not associate with the cluster. The SNPs in question varied only in the HTR2A gene. A possible explanation for this is that there is a considerable difference in the how the SNPs in the HTR2A gene behave between the sexes. We speculated that this finding might be due to a linkage between these SNPs, but the post hoc analysis only revealed a minor linkage between some of these SNPs, thus not explaining the result. The reason for the difference in these analyses remains unclear. One speculation would be that these genes are expressed differently between the two sexes. In a meta-analysis by Miettunen et al. (Reference Miettunen, Veijola, Lauronen, Kantojärvi and Joukamaa2007), sex was observed to have an effect on the HA and RD dimensions, which could also at least partially explain this finding, as the HHA/LRD/LP cluster includes the extremities of both of these dimensions.

The participants in this study were depressed, which can affect the dimensions in the TCI-R questionnaire. Previous studies have shown that HA and P seem to change during the acute phase of MDD compared to the remission phase (Hansenne & Bianchi, Reference Hansenne and Bianchi2009; Kampman & Poutanen, Reference Kampman and Poutanen2011). In this study, we assumed that each individual has a ‘trait’ temperament, as well as a ’state’ temperament profile. The trait temperament is what the individual’s genetics define based on Cloninger’s theory, and the state temperament is the temperament that can be quantified but is influenced by psychiatric and external conditions, such as acute depression. Since this study focused on genetics, the aim was to get as close to the trait temperament profiles as possible. Especially during the acute phase of MDD, which corresponds to the baseline of this study, the state temperament recorded by the TCI-R questionnaire is probably different from individual’s genetic trait temperament. Therefore, the MADRS score was used as an explanatory variable to try to diminish the effect of MDD symptoms on the temperament dimensions. Also, a depressive episode has shown to leave a permanent ’scar’ on an individual’s TCI profile (Kampman & Poutanen, Reference Kampman and Poutanen2011). Striving to study the trait temperament of individuals with psychiatric disorders is somewhat difficult, since the exact effect of psychiatric disorders on an individual’s temperament dimensions is difficult to establish without prior data about the individual’s temperament. However, despite the variety of known risk factors for MDD, we cannot reliably predict future episodes in individuals who end up with MDD. In a clinical setting, patients are assessed after symptoms of psychiatric disorders begin, thus the temperament dimensions collected would probably be affected by these disorders. Thus, either the state temperament or trait temperament altered by MDD might be more relevant in a clinical setting than the pure trait temperaments of healthy individuals.

Clustering the TCI dimensions has not been done widely in previous studies, and no predetermined criteria for these clusters have been established. Despite this, creating a more complex model of an individual’s temperament might be a sensible approach to studying the behaviour, trajectories, and genetic influences in general, as the nature of human temperament is complex. A problem with the clustering used in this study is the applicability of these findings in a clinical environment. SNPs have been widely studied in previous literature separately, as well as in combinations. A strength of this study is using the SNPs as combination, since finding the psychiatric relevance of a single SNP would need a large pool of high-quality data, and even then, the effect and clinical relevance of individual SNPs would most likely be low.

Based on the preliminary literature search, a considerable number of studies were found which have investigated the influence of genetics on temperament, though many of the studies focused on other genetic polymorphisms than SNPs, and some used other methods than the TCI-R to quantify an individual’s temperament. For example, a substantial amount of research has been conducted on the effects of the short and long alleles of SLC6A4 in regard to temperament and psychiatric conditions (Ma et al., Reference Ma, Li, Wang, Shi, Sun, Sheng, Zhang, Zhang, Rao and Han2014; Borkowska et al., Reference Borkowska, Bieliński, Szczęsny, Szwed, Tomaszewska, Kałwa, Lesiewska, Junik, Gołębiewski, Sikora, Tretyn, Akiskal and Akiskal2015). Some SNPs, for example rs2020936 in the SLC6A4 combined with another SNP in the same gene, have been found to have an association with the short allele of SLC6A4 (Wray et al., Reference Wray, James, Gordon, Dumenil, Ryan, Coventry, Statham, Pergadia, Madden, Heath, Montgomery and Martin2009), though this seems to be more of an exception than a rule. As another example, studies have been conducted on the variable number of tandem repeat mutations in the DRD4 gene (Kang et al., Reference Kang, Namkoong and Kim2008; Oniszczenko & Dragan, Reference Oniszczenko and Dragan2012). Due to the available data, this study only included SNPs, so validating and studying other types of gene polymorphisms was not possible.

A limitation of this study is the small number of participants. This may result in the inability to detect subtle clinical differences between individuals with different SNPs. A strength of this study is the inclusion of both inpatients and outpatients, thus including a wide range and severity of symptoms. Even though the participants’ medication was optimised, compliance was only assessed with patient reports.

In general, the possibility of using genetic tests in clinical applications, for example in the treatment of depression, currently has its limitations. In a recent review by Lewis & Vassos (Reference Lewis and Vassos2020), the authors estimated an individual’s risk of depression to be 2% from genetics and 98% from unaccounted variation from unmodelled genetic and environmental factors. When comparing polygenic risk scores, they observed a 2.5-fold increase in risk of depression in the highest decile compared to the lowest decile. At this point in time, the genetic impact on an individuals’ risk of depression seems to be small, but measurable. The clinical significance of measuring an individual’s genome in the context of depression is most likely small but might be an additional tool in the future if genotyping of patients becomes more widespread in a clinical setting, especially if an individual’s genome data would be available before their psychiatric symptoms arise.

In conclusion, associations between combinations of SNPs and HHA/LRD/LP and HP clusters were observed in the multivariate models used in this study. Further research needs to be conducted to validate these findings and to explore the possibility of using these temperament clusters more widely. In addition, researching other types of gene polymorphisms in relation to a cluster model would be beneficial, as well as studying the temperament clusters of individuals with no history of psychiatric disorders.

Supplementary material

The supplementary material for this article can be found at https://doi.org/10.1017/neu.2023.33.

Acknowledgements

We would like to express our sincere gratitude to the research nurses Susanna Hotakainen, Marja Koivumäki and Kati Huhtala, and MD, PhD Antti Koivukangas, MD, PhD Antero Lassila, and MD Kaisa Luoto for their important work on the ODS project.

Author contribution

Simo-Pekko Salminen, Anssi Solismaa, and Olli Kampman were responsible for planning the research. Olli Kampman and Esa Leinonen participated in organising and designing the ODS study and collecting the data. Leo-Pekka Lyytikäinen, Nina Mononen, and Terho Lehtimäki performed the laboratory analyses and genotyping. Simo-Pekko Salminen, Anssi Solismaa, and Leo-Pekka Lyytikäinen conducted the statistical analyses. Simo-Pekko Salminen, Anssi Solismaa, Vesa Paavonen, Terho Lehtimäki, and Olli Kampman contributed to the interpretation of the data. Simo-Pekko Salminen and Anssi Solismaa drafted the manuscript. Vesa Paavonen, Terho Lehtimäki, Leo-Pekka Lyytikäinen, and Olli Kampman edited the manuscript. All authors have reviewed and approved the manuscript.

Financial support

This research was funded by Southern Ostrobothnia Hospital District research fund (EVO1114).

Competing interests

The authors of this study have no conflicts of interest.

References

Akaike, H (1981) Likelihood of a model and information criteria. Journal of Econometrics 16(1), 314.CrossRefGoogle Scholar
Alfimova, MV, Monakhov, MV, Golimbet, VE, Korovaitseva, GI and Lyashenko, GL (2010) Analysis of associations between 5-HTT, 5-HTR2A, and GABRA6 gene polymorphisms and health-associated personality traits. Bulletin of Experimental Biology and Medicine 149(4), 434436.CrossRefGoogle ScholarPubMed
Beck, AT, Guth, D, Steer, RA and Ball, R (1997) Screening for major depression disorders in medical inpatients with the Beck Depression Inventory for Primary Care. Behaviour Research and Therapy 35(8), 785791.CrossRefGoogle ScholarPubMed
Bohn, MJ, Babor, TF and Kranzler, HR (1995) The Alcohol Use Disorders Identification Test (AUDIT): validation of a screening instrument for use in medical settings. Journal of Studies on Alcohol 56(4), 423432.CrossRefGoogle ScholarPubMed
Borkowska, A, Bieliński, M, Szczęsny, W, Szwed, K, Tomaszewska, M, Kałwa, A, Lesiewska, N, Junik, R, Gołębiewski, M, Sikora, M, Tretyn, A, Akiskal, K and Akiskal, H (2015) Effect of the 5-HTTLPR polymorphism on affective temperament, depression and body mass index in obesity. Journal of Affective Disorders 184, 193197.CrossRefGoogle ScholarPubMed
Cloninger, CR (1986) A unified biosocial theory of personality and its role in the development of anxiety states. Psychiatric Developments 4(3), 167226.Google ScholarPubMed
Cloninger, CR, Svrakic, DM and Przybeck, TR (1993) A psychobiological model of temperament and character. Archives of General Psychiatry 50(12), 975990.CrossRefGoogle ScholarPubMed
Cloninger, CR, Przybeck, TR, Svrakic Dragan, M and Wetzel Richard, D (1994) The Temperament and Character Inventory (TCI): A Guide to its Development and Use. St Louis, Missouri: Center for Psychobiology of Personality, Washington University. Google Scholar
Glavina Jelas, I, Devic, I and Karlovic, D (2018) Cloninger’s temperament and character dimensions and dopaminergic genes: DAT1 VNTR and COMT Val158Met polymorphisms. Psychiatria Danubina 30, 4756.CrossRefGoogle ScholarPubMed
Hansen, RA, Moore, CG, Dusetzina, SB, Leinwand, BI, Gartlehner, G and Gaynes, BN (2009) Controlling for drug dose in systematic review and meta-analysis: a case study of the effect of antidepressant dose. Medical Decision Making 29(1), 91103.CrossRefGoogle ScholarPubMed
Hansenne, M and Bianchi, J (2009) Emotional intelligence and personality in major depression: trait versus state effects. Psychiatry Research 166(1), 6368.CrossRefGoogle ScholarPubMed
Hu, P, Muise, AM, Xing, X, Brumell, JH, Silverberg, MS and Xu, W (2013) Association between a multi-locus genetic risk score and inflammatory bowel disease. Bioinformatics and Biology Insights 7, 143152.CrossRefGoogle ScholarPubMed
Kampman, O and Poutanen, O (2011) Can onset and recovery in depression be predicted by temperament? A systematic review and meta-analysis. Journal of Affective Disorders 135(1-3), 2027.CrossRefGoogle ScholarPubMed
Kampman, O, Viikki, M, Jarventausta, K and Leinonen, E (2014) Meta-analysis of anxiety disorders and temperament. Neuropsychobiology 69(3), 175186.CrossRefGoogle ScholarPubMed
Kang, JI, Namkoong, K and Kim, SJ (2008) The association of 5-HTTLPR and DRD4 VNTR polymorphisms with affective temperamental traits in healthy volunteers. Journal of Affective Disorders 109(1-2), 157163.CrossRefGoogle ScholarPubMed
Lee, BC, Yang, JW, Lee, SH, Kim, SH, Joe, SH, Jung, IK, Choi, IG and Ham, BJ (2008) An interaction between the norepinephrine transporter and monoamine oxidase A polymorphisms, and novelty-seeking personality traits in Korean females. Progress in Neuro-Psychopharmacology and Biological Psychiatry 32(1), 238242.CrossRefGoogle ScholarPubMed
Lewis, CM and Vassos, E (2020) Polygenic risk scores: from research tools to clinical instruments. Genome Medicine 12(1), 44.CrossRefGoogle ScholarPubMed
Ma, Y, Li, B, Wang, C, Shi, Z, Sun, Y, Sheng, F, Zhang, Y, Zhang, W, Rao, Y and Han, S (2014) 5-HTTLPR polymorphism modulates neural mechanisms of negative self-reflection. Cerebral Cortex 24(9), 24212429.CrossRefGoogle ScholarPubMed
Miettunen, J, Veijola, J, Lauronen, E, Kantojärvi, L and Joukamaa, M (2007) Sex differences in Cloninger’s temperament dimensions--a meta-analysis. Comprehensive Psychiatry 48(2), 161169.CrossRefGoogle ScholarPubMed
Montgomery, SA and Asberg, M (1979) A new depression scale designed to be sensitive to change. British Journal of Psychiatry 134(4), 382389.CrossRefGoogle ScholarPubMed
Oniszczenko, W and Dragan, WŁ. (2012) Association between temperament in terms of the Regulative Theory of Temperament and DRD4 and DAT1 gene polymorphisms. Comprehensive Psychiatry 53(6), 789796.CrossRefGoogle ScholarPubMed
Paavonen, V, Kampman, O, Illi, A, Viikki, M, Setälä-Soikkeli, E and Leinonen, E (2014) A cluster model of temperament as an indicator of antidepressant response and symptom severity in major depression. Psychiatry Investigation 11(1), 1823.CrossRefGoogle ScholarPubMed
Paavonen, V, Luoto, K, Koivukangas, A, Lassila, A, Leinonen, E and Kampman, O (2016) Temperament and character profiles associated with depression and treatment response in patients with or without comorbid substance abuse. Psychiatry Research 245, 250258.CrossRefGoogle ScholarPubMed
Paavonen, V, Luoto, K, Lassila, A, Leinonen, E and Kampman, O (2018) Temperament clusters associate with anxiety disorder comorbidity in depression. Journal of Affective Disorders 236, 252258.CrossRefGoogle ScholarPubMed
Sheehan, DV, Lecrubier, Y, Sheehan, KH, Amorim, P, Janavs, J, Weiller, E, Hergueta, T, Baker, R and Dunbar, GC (1998) The Mini-International Neuropsychiatric Interview (M.I.N.I.): the development and validation of a structured diagnostic psychiatric interview for DSM-IV and ICD-10. Journal of Clinical Psychiatry 59, 2257.Google ScholarPubMed
Solismaa, A, Kampman, O, Lyytikäinen, LP, Seppälä, N, Viikki, M, Mononen, N, Lehtimäki, T and Leinonen, E (2017) Histaminergic gene polymorphisms associated with sedation in clozapine-treated patients. European Neuropsychopharmacology 27(5), 442449.CrossRefGoogle ScholarPubMed
Suzuki, A, Matsumoto, Y, Ishii, G, Oshino, S, Goto, K and Otani, K (2007) No association between the -3081A/T polymorphism in the norepinephrine transporter gene promoter and personality traits in healthy subjects. Neuroscience Letters 425(3), 192194.CrossRefGoogle ScholarPubMed
Tsuchimine, S, Taniguchi, T, Sugawara, N, Kaneda, A and Yasui-Furukori, N (2013) No association between a polymorphism in the serotonin receptor 2B (HTR2B) gene and personality traits in healthy Japanese subjects. Neuropsychobiology 68(1), 5962.CrossRefGoogle ScholarPubMed
Tsuchimine, S, Yasui-Furukori, N, Sasaki, K, Kaneda, A, Sugawara, N, Yoshida, S and Kaneko, S (2012) Association between the dopamine D2 receptor (DRD2) polymorphism and the personality traits of healthy Japanese participants. Progress in Neuro-Psychopharmacology and Biological Psychiatry 38(2), 190193.CrossRefGoogle ScholarPubMed
Wray, NR, James, MR, Gordon, SD, Dumenil, T, Ryan, L, Coventry, WL, Statham, DJ, Pergadia, ML, Madden, PA, Heath, AC, Montgomery, GW and Martin, NG (2009) Accurate, large-scale genotyping of 5HTTLPR and flanking single nucleotide polymorphisms in an association study of depression, anxiety, and personality measures. Biological Psychiatry 66(5), 468476.CrossRefGoogle Scholar
Zaninotto, L, Solmi, M, Toffanin, T, Veronese, N, Cloninger, CR and Correll, CU (2016) A meta-analysis of temperament and character dimensions in patients with mood disorders: comparison to healthy controls and unaffected siblings. Journal of Affective Disorders 194, 8497.CrossRefGoogle ScholarPubMed
Figure 0

Table 1. Demographics of each temperament cluster including age, sex, smoking status, MADRS and AUDIT scores at the beginning of the study and antidepressant doses in fluoxetine equivalents

Figure 1

Table 2. Logistic regression models with different temperament clusters as dependent variables, explanatory variables used in the model and results of the permutation analysis

Figure 2

Table 3. Linkage disequilibrium (LD) with other SNPs in the sample and annotations of SNPs in the GRSs found to be statistically significant after permutations

Supplementary material: File

Salminen et al. supplementary material

Salminen et al. supplementary material

Download Salminen et al. supplementary material(File)
File 36 KB