Original Research ARTICLE
Association of the Polygenic Scores for Personality Traits and Response to Selective Serotonin Reuptake Inhibitors in Patients with Major Depressive Disorder
- 1Discipline of Psychiatry, School of Medicine, University of Adelaide, Adelaide, SA, Australia
- 2Northern Adelaide Local Health Network, Mental Health Services, Adelaide, SA, Australia
- 3Epidemiology Branch, Division of Intramural Population Health Research, National Institute of Child Health and Human Development, National Institutes of Health, Bethesda, MD, United States
- 4HSL Institute for Aging Research, Harvard Medical School, Boston, MA, United States
- 5Program for Quantitative Genomics, Harvard School of Public Health, Boston, MA, United States
- 6Broad Institute of MIT and Harvard, Cambridge, MA, United States
- 7Biomedical Data Science, Stanford University, Stanford, CA, United States
- 8Department of Health Sciences Research, Mayo Clinic, Rochester, NY, United States
- 9Department of Bioengineering, Stanford University, Stanford, CA, United States
- 10Department of Psychiatry and Psychotherapy, University of Muenster, Muenster, Germany
- 11Department of Clinical Pharmacology, University Göttingen, Göttingen, Germany
- 12Department of Psychiatry, Taipei Medical University-Shuangho Hospital, New Taipei City, Taiwan
- 13Department of Psychiatry and Psychotherapy, Faculty of Medicine, University of Freiburg, Freiburg, Germany
- 14Department of Psychiatry and Psychology, Mayo Clinic, Rochester, NY, United States
- 15Department of Psychiatry, Taipei Veterans General Hospital, Taipei, Taiwan
- 16Division of Psychiatry, School of Medicine, National Yang-Ming University, Taipei, Taiwan
- 17Department of Psychiatry, Faculty of Medicine and Life Sciences, University of Tampere, Tampere, Finland
- 18Department of Molecular Pharmacology and Experimental Therapeutics, Mayo Clinic Rochester, Rochester, MN, United States
- 19Department of Psychiatry, Seinäjoki Hospital District, Seinäjoki, Finland
- 20Department of Neuropsychiatry, Kansai Medical University, Osaka, Japan
- 21Department of Psychiatry, Tampere University Hospital, Tampere, Finland
- 22RIKEN Center for Integrative Medical Sciences, Kanagawa, Japan
- 23Department of Pharmacy, Hyogo University of Health Sciences, Hyogo, Japan
- 24Center for Neuropsychiatric Research, National Health Research Institutes, Miaoli, Taiwan
- 25Center for Medical Genetics Research, Rajanukul Institute, Department of Mental Health, Ministry of Public Health Bangkok, Bangkok, Thailand
- 26Research Division Federal Institute for Drugs and Medical Devices, Bonn, Germany
Studies reported a strong genetic correlation between the Big Five personality traits and major depressive disorder (MDD). Moreover, personality traits are thought to be associated with response to antidepressants treatment that might partly be mediated by genetic factors. In this study, we examined whether polygenic scores (PGSs) derived from the Big Five personality traits predict treatment response and remission in patients with MDD who were prescribed selective serotonin reuptake inhibitors (SSRIs). In addition, we performed meta-analyses of genome-wide association studies (GWASs) on these traits to identify genetic variants underpinning the cross-trait polygenic association. The PGS analysis was performed using data from two cohorts: the Pharmacogenomics Research Network Antidepressant Medication Pharmacogenomic Study (PGRN-AMPS, n = 529) and the International SSRI Pharmacogenomics Consortium (ISPC, n = 865). The cross-trait GWAS meta-analyses were conducted by combining GWAS summary statistics on SSRIs treatment outcome and on the personality traits. The results showed that the PGS for openness and neuroticism were associated with SSRIs treatment outcomes at p < 0.05 across PT thresholds in both cohorts. A significant association was also found between the PGS for conscientiousness and SSRIs treatment response in the PGRN-AMPS sample. In the cross-trait GWAS meta-analyses, we identified eight loci associated with (a) SSRIs response and conscientiousness near YEATS4 gene and (b) SSRI remission and neuroticism eight loci near PRAG1, MSRA, XKR6, ELAVL2, PLXNC1, PLEKHM1, and BRUNOL4 genes. An assessment of a polygenic load for personality traits may assist in conjunction with clinical data to predict whether MDD patients might respond favorably to SSRIs.
A major depressive disorder (MDD) is the most common and disabling mental health diseases worldwide (1, 2) with a lifetime prevalence of ~12% (3). Studies estimated a 61.6 million years of life lived with disability caused by MDD accounting for 2.5% of the total disability-adjusted life years and for 8.1% of the total years lived with disability resulted from all diseases (2, 4).
Selective serotonin reuptake inhibitors (SSRIs) are commonly used as the first-line pharmacological treatment for MDD (5). However, treatment efficacy with SSRIs varies widely between individual patients and is inadequate in many cases. Clinical response rates range from 48 to 64% (6, 7) and reported remission rates are as low as 23.5% (7, 8). To improve this situation, an investigation of the biological and psychosocial factors that drive heterogeneity in treatment outcomes is necessary.
There is growing evidence from genetic studies that antidepressant treatment response is substantially influenced by genes (7, 9–17). A study involving nearly 3,000 MDD patients estimated that genetic factors explain 42% of the differences in the level of treatment response (18). A number of genes and single nucleotide polymorphisms (SNPs) that could influence antidepressant treatment outcomes have been reported, including polymorphisms within the COMT (9), HTR2A (10), HTR1A (11), CNR1 (11), SLC6A4 (12), NPY (13), MAOA (14), and IL1B (15) genes. A pharmacogenomic study on SSRIs response by the International SSRIs Pharmacogenomics Consortium (ISPC) identified several SNPs with suggestive association after 4 weeks of treatment, including the neuregulin-1 gene, which is involved in many aspects of brain development, such as neuronal maturation (7).
In addition to genetic factors, multiple demographic, clinical, and psychological predictors of SSRI response in MDD have been identified, collectively explaining 5–15% of the variance in treatment outcomes (19–23). Among the psychological predictors, personality traits defined by the Five-Factor Model of Personality (“Big Five”: extraversion, agreeableness, conscientiousness, neuroticism, openness) (24) have previously been reported to influence antidepressant treatment response and remission (25–29). Of these, neuroticism is a frequently reported predisposing factor for depression and was shown to negatively affect antidepressants treatment response (30, 31). In a recent study, MDD patients resistant to antidepressants were more likely to report high clinical scores for neuroticism, but low scores for openness, conscientiousness, and extraversion (26). In a large study of patients with MDD (n = 8,229), pre-existing personality dysfunction was associated with poor response to antidepressants (27). Further, some studies have suggested that SSRIs have a direct positive impact on scores for neuroticism or extraversion in MDD patients, and that part of the antidepressant effect might be explained through these adjustments (28, 29, 32, 33). Moreover, shared genes are thought to play a key role in the association between personality factors and MDD (34). For example, studies have estimated the genetic correlation between MDD and neuroticism at 55–75% (35, 36). However, no previous work has directly addressed the question whether there is a genetic relationship between the Big Five personality traits and SSRI treatment response and remission in MDD. It has been shown that the genetic architecture of personality traits is highly polygenic, in which several genes of small effect contribute to the overall phenotype (35, 37). Thus, a polygenic score (PGS) analysis approach proposed by the schizophrenia consortium (38) and later applied in several studies (16, 39), is potentially powerful to investigate the genetic influence of each of the Big Five personality traits on antidepressant treatment outcomes. A PGS for each of the Big five personality traits quantifies the combined effects of genetic variants across the whole genome, computed as a weighted summation of effect sizes obtained from genome-wide association studies (GWASs). A successful multi-trait polygenic model may assist for an early screening of diseases risk, clinical diagnosis, and the prediction of treatment response and prognosis (38, 39).
Implicitly, one could also interpret a polygenic association as a biological relationship partly explained by the role of shared genes and common molecular mechanisms. With this in mind, we conducted GWAS meta-analyses by combining GWAS summary statistics on the Big Five personality traits and SSRIs treatment outcome to identify shared genes involved in the cross-trait association.
Materials and Methods
The characteristics of the clinical and genetic data, as well as the sources of the GWAS summary statistics used in our analysis are described below.
Pharmacogenomics Research Network Antidepressant Medication Pharmacogenomic Study (PGRN-AMPS)
The PGRN-AMPS is a clinical trial on the response to escitalopram or citalopram of 529 MDD patients over 8 weeks of treatment. The baseline and follow-up assessment of depression severity were performed using the 16-item Quick Inventory of Depressive Symptomatology (QIDS-C16) (40).
The ISPC is an International Consortium established to discover genes that are responsible for SSRIs treatment response in patients with MDD. For our study, we used data from 865 MDD patients recruited in the USA, Germany, Thailand, Taiwan, and Japan who received SSRI treatment. The 17-item Hamilton Depression Rating Scale was used as a measurement tool to assess and follow-up the treatment progress (7).
Genotyping and Quality Control
The genotype and clinical data for the PGRN-AMPS were available via a controlled access system at the database of Genotypes and Phenotypes: dbGaP1 and the ISPC data were obtained from the ISPC consortium (7).
For the genotype data of both samples, we implemented quality control (QC) steps using PLINK (41) and samples with low genotype rates <95%, sex inconsistencies (X-chromosome heterozygosity), and genetically related individuals were excluded. We also excluded SNPs that had poor genotyping rate <95%, an ambiguity (A/T and C/G SNPs), a minor allele frequency (MAF ≤ 1%), or showed deviation from Hardy–Weinberg Equilibrium (p < 10–6).
After excluding the low-frequency SNPs (MAF < 10%), poor-quality variants (imputation INFO <0.9 and indels), the imputed dosages were converted to best guess genotypes. The subsequent PGS analyses were performed using the best guess genotypes.
GWAS Summary Statistics Data
The PGSs were calculated using the approach previously described by the International Schizophrenia Consortium (38). This method requires an estimated effect size for each SNP to compute weighted PGS. The effect estimates (betas) for this study were the summary statistics obtained from previously published GWASs on extraversion, openness, agreeableness, conscientiousness (37), and on neuroticism (35). The data were publicly available for download at http://www.tweelingenregister.org/GPC/ and http://www.thessgac.org/data, respectively. The effect size estimates for each SNP—quantified as beta was extracted from the download file and used to compute weighted PGS in the PGRN-AMPS and ISPC cohorts.
Definition of SSRI Treatment Outcomes
Treatment response and remission to SSRIs were defined after 4 weeks of treatment follow-up of MDD patients in both cohorts. In addition, PGS associations were evaluated at 8 weeks in PGRN-AMPS. While treatment response was determined as a ≥50% reduction from baseline in the HRSD-17 or QIDS-C16 total scores, SSRI treatment remission was defined as achieving a HRSD-17 score ≤7 or a QIDS-C16 score ≤5 at 4 or 8 weeks of treatment.
PGS Computation and Association Analyses
The PGSs were computed for each of the Big Five personality traits using imputed genetic data weighted by GWAS summary statistics of the respective personality traits, separately for the two cohorts: PGRN-AMPS (n = 529) and ISPC (n = 865) (Table 1; Figure 1). First, quality-controlled SNPs were clumped for linkage disequilibrium (LD) using genome-wide association p-value informed clumping with r2 = 0.1 in a 250-kb window to create an independent SNP-set using PLINK software run on Linux. Next, weighted PGSs were calculated for each individual at a range of p-value thresholds (PT) as a weighted sum of allele dosages (0, 1, or 2). The PT refers to the p-values associated with the effect size of each of the SNPs, as listed in the GWAS summary statistics (35, 37). The weighting was performed by multiplying the dosage of each effect increasing allele by its effect size derived from the GWAS summary statistics (β-coefficient), then divided by the total number of SNPs in each threshold. The PGS was computed at a range of PT (<1 × 10–2, <5 × 10–2, <0.1, <0.2, <0.3, <0.4, <0.5, and <1.0) separately for each of the two cohorts. Performing the PGS at different PT provides a range of alternative scores to choose the most significantly associated (optimal) PGS that will be used for prediction modeling. At each PT, a logistic regression modeling was applied to response/remission to SSRIs (dependent variables) using the PGS for each of the Big Five personality traits as the main predictor variable and adjusting for common covariates, such as age, sex, and cohort-specific covariates including four principal components in the PGRN-AMPS and “study sites” in the ISPC. A statistically significant association between the PGSs for the Big Five personality traits and response/remission to SSRIs was determined at p < 0.05, across the PT in both study samples. The prediction accuracy, the percentage of variance explained, Nagelkerke R2, by the PGSs were calculated as the Nagelkerke R2 of the full model with PGS and covariates minus the Nagelkerke R2 of the model with only covariates. To determine the effect of high or low polygenic load on treatment outcomes, the study subjects were grouped into PGS quartiles (Q1–Q4) at the optimal PT. Then, we estimated the odds of treatment response/remission to SSRIs for MDD patients within the group with a high polygenic load for the Big Five personality traits (Q2, Q3, Q4) compared to patients in the lowest PGS quartile (Q1).
Table 1. Baseline characteristics of major depressive disorder patients and their treatment outcomes with selective serotonin reuptake inhibitors after 4 weeks of follow-up.
Figure 1. The bar graphs (A–C) show the association of the PGSs for the Big Five personality traits with SSRIs response or remission at different p-value thresholds (PT) after 4 weeks of treatment in the PGRN-AMPS (n = 529) and ISPC (n = 865) samples. The y-axis (Nagelkerke’s Pseudo-R2) refers to the percentage of variance in SSRIs treatment response/remission accounted for the PGSs of the Big Five personality traits at a particular PT in each sample. On the x-axis, plotted from left to right, are the GWAS PT for personality traits used to group the SNPs for the PGSs. The *sign on the top of each bar signify the statistical significance of the PGS association as *p < 0.05, **p < 0.01, ***p < 0.001. Abbreviations: PGRN-AMPS, the Pharmacogenomics Research Network Antidepressant Medication Pharmacogenomic Study; ISPC, the International SSRI Pharmacogenomics Consortium study; SNP, single nucleotide polymorphism; PGS, polygenic score; SSRIs, selective serotonin reuptake inhibitors.
Cross-Trait Meta-Analyses of GWASs
In the cross-trait meta-analyses, we applied the O’Brien’s (OB) method and the direct Linear Combination of dependent test statistics (dLC) approach (39, 44, 45), which are implemented in the C++ eLX package. Briefly, the OB method and the dLC approach help to combine GWAS effect estimates of genome-wide SNPs, obtained from univariate GWASs and generated two test statistics and associated p-values—one for the OB method and one for the dLC method. More details can be found elsewhere (44, 45). The eLX package is available at https://sites.google.com/site/multivariateyihsianghsu/.
Here, GWAS on personality traits that have shown a significant association in the PGS analysis were combined with GWAS on SSRIs treatment outcome. The GWAS summary statistics on SSRIs treatment response (7) were combined with those on (i) conscientiousness (34) and (ii) openness personality (34). Similarly, the GWAS summary statistics on SSRIs treatment remission (7) was meta-analyzed with (i) openness personality (34) and (ii) neuroticism (35).
Statistical significance was determined based on the smaller of the OB or the dLC p-values. A significant association was determined if (1) the p-value for the cross-trait meta-analysis reached genome-wide significance (p < 5 × 10−8) and (2) the univariate GWAS effects were at least nominally significant (p < 0.05). For each cross-trait meta-analysis, only one lead SNP per locus was reported. Nearby SNPs in LD (r2 > 0.1) with the lead SNP were considered dependent and belonging to the same locus.
Patient Characteristics and Treatment Outcomes
In this study, we analyzed data from 1,394 MDD patients who had SSRI treatment divided into PGRN-AMPS (n = 529) and ISPC (n = 865) samples. The average age of the patients was 42.2 years and the majority of them (64.3%) were females (Table 1).
Of all patients, 622 (46.8%) were classified as treatment responders with a slight variation across the study samples 44.4% in the PGRN-AMPS and 48.1% in the ISPC. Remission rates were 27.6 and 26.1% in the PGRN-AMPS and ISPC samples, respectively. The rate of remission combined across the two studies was 26.7% (Table 1).
Association of the PGS for the Big Five Personality Traits with SSRIs Treatment Outcomes
Polygenic scores were computed for each of the Big Five personality traits, and we investigated their association with two SSRI treatment outcomes—response and remission, after 4 weeks (PGRN-AMPS and ISPC) and 8 weeks (PGRN-AMPS) of treatment.
After 4 weeks of treatment, genetic predisposition to openness, conscientiousness, and neuroticism were associated with SSRIs treatment response and/or remission at p < 0.05 across PT thresholds, in at least one of the two assessed cohorts (Figures 1A–C). Genetic loading for openness was associated with response and remission in both cohorts (Figure 1A1,2). An elevated PGS for conscientiousness was associated with treatment response, but not remission, in the PGRN-AMPS sample only (Figure 1B). A PGS association for neuroticism with remission, but not treatment response, was shown in both cohorts (Figure 1C). The PGSs for extraversion and agreeableness were associated with neither response nor remission.
We also assessed the level of observed variation in SSRI treatment outcomes accounted for by these personality traits, and found that personality traits at the most significant thresholds explained a considerable amount of variance in treatment outcomes. For example, the PGS for openness accounted for ~1.5% of the observed variation in SSRIs treatment response and ~2.8% of the variance in remission. The PGS for neuroticism explained ~1.5% of the variance in remission and the PGS for conscientiousness contributed to ~1.5% of the variability in SSRI treatment response.
The status of treatment response and remission for patients in personality trait quartiles (Q2–Q4) was compared with those in the lowest personality trait PGS quartile (Q1) (Figure 2). Our analysis revealed that MDD patients with a high polygenic load for openness personality had initially poorer remission and response rates at 4 weeks of treatment, with Q4 versus Q1 odds ratios (ORs) ranging from 0.30 [ISPC: 95%CI, 0.15–0.59] to 0.52 [PGRN-AMPS: 95%CI, 0.29–0.90] (Figure 2A1,2, green and brown graphs). After longer treatment duration, we observed a reverse effect. Here, a higher polygenic load for openness was associated with a better SSRIs treatment response at 8 weeks in the PGRN-AMPS, with OR of 1.58 [95%CI, 1.10–2.90] (Figure 2A1,2, blue graphs).
Figure 2. The line plots represent the ORs for favorable response or remission to selective serotonin reuptake inhibitors treatment in patients with MDD with a high personality traits polygenic load (Q2, Q3, and Q4) compared to patients with the lowest polygenic load (Q1), estimated at the most significant p-value thresholds. The quartile-based polygenic scores analyses were performed using data at 4 weeks in the ISPC and at 4 and 8 weeks in the PGRN-AMPS. A polygenic loading for openness personality trait was initially associated with poor response and remission to selective serotonin reuptake inhibitors (SSRIs) in the first 4 weeks of treatment (ISPC, PGRN-AMPS at 4 weeks). After a longer (8 weeks) treatment follow-up, the genetic loading for openness had shown a favorable effect to SSRIs response and remission (PGRN-AMPS at 8 weeks). The polygenic loading for conscientiousness personality was favorably associated with response to SSRIs treatment. However, a polygenic loading for neuroticism personality had shown a negative impact on SSRIs remission. The ORs are reported on the lines and the *sign indicates the statistical significance of the ORs as *p < 0.05, **p < 0.01, ***p < 0.001. Abbreviations: PGRN-AMPS, the Pharmacogenomics Research Network Antidepressant Medication Pharmacogenomic Study; ISPC, the International SSRI Pharmacogenomics Consortium study. OR, odds ratio; Q1, quartiles 1; Q2, quartiles 2; Q3, quartiles 3; Q4, quartiles 4; MDD, major depressive disorder.
Major depressive disorder patients with a higher polygenic load for conscientiousness personality had 1.95 [95% CI, 1.13–3.36] times better SSRIs treatment response compared to those patients in the lowest PGS, although this association was only significant in the PGRN-AMPS sample at 4 weeks of treatment (Figure 2B).
Conversely, MDD patients with a higher polygenic load for neuroticism personality had poorer treatment outcomes with SSRIs. After 4 weeks of treatment, patients in Q4 based on the PGS for neurotic personality had about 50% lower odds of remission compared to patients in Q1 with OR ranging from 0.50 [PGRN-AMPS: 95%CI, 0.28–0.90] to 0.54 [ISPC: 95%CI, 0.33–0.89] (Figure 2C). Constantly, results after 8 weeks of treatment showed a trend inverse association between the PGS for neurotic personality and SSRIs treatment remission, although this was not statistically significant (Figure 2C).
To assess the potential effect of false-positive findings, the association p-values were corrected for multiple testing at each PT for SSRIs treatment response and remission using the Benjamini and Hochberg (BH) method. Each of the p-values was adjusted assuming a conventionally accepted level of 5% false discovery rate (FDR) (46). After FDR adjustment, the associations of the PGS for openness personality with SSRIs treatment response remained statistically significant (in the ISPC sample: FDR adjusted p-value = 0.02 at PT < 1 × 10−2) and with remission (in the PGRN-AMPS sample: FDR adjusted p-value = 0.04 at PT < 5 × 10−2). The PGSs for conscientiousness and neuroticism were not associated with SSRIs treatment outcome after implementing the FDR adjusted p-value <0.05.
Cross-Trait Meta-Analyses of GWASs
For personality traits that showed a significantly associated PGS, cross-trait GWAS meta-analyses was performed by combining summary GWAS data on SSRIs treatment outcomes and personality traits. Table 2 and Figure 3 summarize the cross-trait meta-analyses findings, including the list of genetic loci and nearest genes that are potentially overlapping between the traits. At a p-value of <5 × 10−8, we identified eight genetic loci located within or near to protein-coding genes with possible overlapping effects on SSRIs treatment outcomes and personality traits. We found (i) one locus associated with conscientiousness and SSRI response near the YEATS4 gene (Table 2; Figure 3A) and (ii) seven loci associated with remission and neuroticism located at or near PRAG1, MSRA, XKR6, ELAVL2, PLXNC1, PLEKHM1, and BRUNOL4 genes (Table 2; Figure 3B). From the meta-analyses of SSRIs treatment outcomes with openness personality, we identified only suggestive evidence at significance p < 1 × 10−6 (Table 2).
Table 2. Significant loci resulting from the cross-trait meta-analyses of genome-wide association studies (GWASs) on selective serotonin reuptake inhibitors (SSRIs) treatment response/remission and GWAS on the Big Five personality traits at univariate GWAS p-value <5 × 10−2 and Cross-trait meta-analysis p-value <5 × 10−8.
Figure 3. The Manhattan plots show the results of the cross-trait meta-analysis of genome-wide association studies (GWASs) on selective serotonin reuptake inhibitors treatment outcomes (response or remission) with GWASs on: (A) conscientiousness personality trait; (B) neuroticism personality, highlighting the loci that showed genome-wide significance (orange), and the nearest genes. The −log10 (cross-trait p-value) is plotted against the physical position of each SNP on each chromosome. The threshold for genome-wide significance (cross-trait p-value <5 × 10−8) is indicated by the red dotted horizontal line.
In this study, we analyzed data from 1,394 MDD patients who had been treated with SSRIs and assessed whether it is possible to predict antidepressants treatment outcomes—response and remission, using PGS for the Big Five personality traits. To further validate the PGS association findings and provide additional evidence, cross-trait meta-analyses of GWASs on SSRIs treatment outcomes versus GWASs on the Big Five personality traits were performed. Our findings from both analyses found complementary evidence that the association of the Big Five personality traits with SSRIs treatment outcomes is partly genetic.
Among the Big Five personality traits, the PGS for openness, conscientiousness, and neuroticism were significantly associated with SSRI treatment outcomes in patients with MDD. A high polygenic load for openness predicted poorer odds of response and remission to SSRIs after 4 weeks of treatment. However, after 8 weeks of treatment, the odds of response and remission was reversed and high loading for openness was associated with favorable outcomes. Patients with a high polygenic load for conscientiousness had a better odd of response to SSRIs after 4 weeks of treatment, but were neither more nor less likely to have good outcomes after 8 weeks. In contrast, patients who possessed a higher polygenic load for neuroticism risk genetic variants responded poorer to SSRIs treatment at both time points.
The discrepancy between short-term and intermediate-term treatment outcomes in patients with high polygenic loading for openness was unexpected in the context of the previous literature (26, 27), and raises the question whether statements about personality impact on SSRI treatment outcomes can be reliably reached on the basis of assessments conducted within the first month. While longitudinal studies of treatment outcomes in MDD suggest that treatment response within the first month occurs for a majority of patients who will eventually remit (47), they also indicate that there is a considerable proportion of patients who achieve response and remission after much longer treatment periods (48, 49). In this context, our finding raises the possibility that the different Big Five personality traits could have differential effects on early- versus delayed responses to treatment in MDD.
Moreover, the inconsistences in the direction of the relationship between the Big Five personality traits and response to long-term versus short-term treatment to SSRIs might be explained by a psychological theory (50–52). Studies suggested that antidepressants have a primary effect on emotional processing, providing a platform for long-term cognitive and psychological recovery (50), and the clinical effects of antidepressant treatment may be mediated by early changes in emotional processing (51, 52).
In our data, consistency between the outcome parameters—treatment response and remission was variable. Only the PGS for openness showed a significant association with both treatment response and remission. The PGS for conscientiousness was associated with better treatment response, but not with remission. The PGS for neuroticism predicted lower odds of treatment remission, but not poorer treatment response. At face value, these findings suggest that openness and neuroticism could play more important roles in predicting ultimate remission from depressive episodes, whereas conscientiousness might drive early treatment effects rather than longer term outcomes. However, another explanation is that our cohorts might have been underpowered to detect more consistent effects, or that some of the observed associations were chance findings, perhaps driven by multiple testing. Indeed, only the associations of the PGS for openness personality with SSRIs treatment response remained statistically significant after FDR adjustment. Therefore, future genetic studies with higher patients’ numbers are required to confirm our findings.
In all, our genetic findings are in line with previous clinical investigations of the influence of personality characteristics on antidepressant treatment response in MDD. A study in Japan revealed as depressed patients who were resistant to treatment had a higher neuroticism score and lower scores for openness, conscientiousness, and extraversion than patients who remitted and healthy controls (26). In another study, higher clinical scores for openness at baseline were associated with improved treatment response to antidepressants, whereas a higher score for neuroticism was associated with poor treatment outcomes (53). More generally, poor treatment response was associated with personality dysfunction in a large sample study of more than 8,000 antidepressant-treated adults with MDD (27). Similarly, a meta-analysis of 34 clinical studies concluded that MDD patients with a comorbid personality disorder had double the risk of overall poor clinical and treatment outcomes, compared to patients no co-occurring personality disorder (54).
Additionally, previous studies have shown genetic correlations between Big Five personality traits and psychiatric disorders and the PGS for neuroticism was significantly associated with MDD (55).
Since the PGS association reflects a shared genetic etiology, we applied cross-trait GWAS meta-analyses by combining summary statistics on SSRI treatment outcomes with personality traits, and identified eight overlapping genetic loci. The YEATS4 gene locus was associated with treatment response to SSRIs and conscientiousness. Previously, a gene expression analysis in depressed patients further replicated in mice found lower levels of YEATS4 in depressed patients compared to healthy controls. Moreover, the expression level of this gene was correlated with the dose of imipramine (a tricyclic antidepressant) (56).
The second gene locus (rs144733372) in PLEKHM, which was found in the cross-trait meta-analysis of neuroticism and SSRIs treatment remission, is highly linked (LD: r2 > 0.8) with several other SNPs located within the CRHR1 gene. The CRHR1 gene encodes a G-protein coupled receptor that binds with the neuropeptides of the corticotrophin-releasing hormone family, a major regulator of the hypothalamic–pituitary–adrenal pathway (57). Functional gene polymorphisms in the CRHR1 gene have been associated with SSRIs treatment response (58), and it moderates the association of maltreatment with neuroticism (59). Corticotrophin-releasing hormone signaling has previously been implicated in mood disorders and treatment response to antidepressants (60).
Another gene showing shared associations with SSRI treatment response and neurotic personality is MSRA, which has shown the highest levels of expression in brain tissue (61). Previous studies reported that genetic variants within the MSRA gene could be associated with schizophrenia, bipolar disorder (62, 63), executive cognitive function (64), fluid intelligence (63), and self-reported irritable temperament (65).
Further, loci within the PRAG1 and PLXNC1 genes have shown overlapping influence on SSRI treatment and neuroticism personality. A genetic polymorphism rs706895C/T within the FYN gene belonging to the same family of genes (tyrosine protein kinase family) was significantly associated with personality traits (66). SNPs within the plexin family gene PLXNA2 have previously been implicated in neuroticism, depression, and psychological distress (67).
Overall, these findings lend further weight to our PGS analyses and reinforce the idea that certain gene polymorphisms have a dual impact on personality structure and antidepressant treatment outcomes in MDD. Studying the individual mechanism of each significant genetic locus in relation to antidepressants in the future studies might lead to novel insights in the molecular underpinnings of these drugs. In conclusion, our study provides evidence in the potential ability of the PGS for the Big Five personality traits to elucidate shared biological mechanisms and to predict SSRI treatment outcomes. Whether these PGSs could be applied to everyday clinical practice in the future relies on their ability to stratify MDD patients into categories of good treatment responders versus nonresponders. Further research is required to determine if this is the case. However, the small effect sizes found in our study give rise to cautious interpretation. In our view, their full clinical value likely lies in their contribution to multi-variable models that also comprise clinical and environmental factors influencing medication response.
AA, KS, and BB developed the study concept and design, performed the statistical analysis, and drafted the manuscript. Other authors contributed data, resources and were involved in critical revision of the manuscript, obtained funding, and contributed to cohort and genetic data and study supervision.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The authors are grateful to the study subjects who participated in the studies, and we appreciate the contributions of research staffs who helped in the patient recruitment and data collection for the studies. The authors also would like to thank the National Institutes of Health (NIH), USA for making the PGRN-AMPS accessible to us. The complete clinical data for the ISPC is available at http://www.pharmgkb.org/downloads/. The PGRN-AMPS data were obtained through controlled access distributed from the NIH in the dbGaP (https://www.ncbi.nlm.nih.gov/gap). The analysis of this study was carried out using the high-performance computational capabilities of the University of Adelaide, Phoenix supercomputer (https://www.adelaide.edu.au/phoenix/).
AA received a Postgraduate Research Scholarship support from the University of Adelaide through the Adelaide Scholarship International program. Funding support for the PGRN-AMPS was provided by the National Institute of General Medical Sciences, National Institutes of Health, through the PGRN grant to Principal Investigators RW and LW (U19 GM61388). Dr. D. Mrazek served as the Principal Investigator for the PGRN-AMPS study within the Mayo Clinic PGRN program. The main sources of funding for the ISPC study are presented in the earlier publication (7).
1. Whiteford HA, Degenhardt L, Rehm J, Baxter AJ, Ferrari AJ, Erskine HE, et al. Global burden of disease attributable to mental and substance use disorders: findings from the global burden of disease study 2010. Lancet (2013) 382(9904):1575–86. doi:10.1016/S0140-6736(13)61611-6
2. Global Burden of Disease Study 2013 Collaborators. Global, regional, and national incidence, prevalence, and years lived with disability for 301 acute and chronic diseases and injuries in 188 countries, 1990-2013: a systematic analysis for the global burden of disease study 2013. Lancet (2015) 386(9995):743–800. doi:10.1016/S0140-6736(15)60692-4
3. Andrade L, Caraveo-Anduaga JJ, Berglund P, Bijl RV, De Graaf R, Vollebergh W, et al. The epidemiology of major depressive episodes: results from the International Consortium of Psychiatric Epidemiology (ICPE) surveys. Int J Methods Psychiatr Res (2003) 12(1):3–21. doi:10.1002/mpr.138
4. Ferrari AJ, Charlson FJ, Norman RE, Patten SB, Freedman G, Murray CJ, et al. Burden of depressive disorders by country, sex, age, and year: findings from the global burden of disease study 2010. PLoS Med (2013) 10(11):e1001547. doi:10.1371/journal.pmed.1001547
5. Hieronymus F, Emilsson JF, Nilsson S, Eriksson E. Consistent superiority of selective serotonin reuptake inhibitors over placebo in reducing depressed mood in patients with major depression. Mol Psychiatry (2016) 21(4):523–30. doi:10.1038/mp.2015.53
7. Biernacka JM, Sangkuhl K, Jenkins G, Whaley RM, Barman P, Batzler A, et al. The International SSRI Pharmacogenomics Consortium (ISPC): a genome-wide association study of antidepressant treatment response. Transl Psychiatry (2015) 5:e553. doi:10.1038/tp.2015.47
8. Papakostas GI, Fava M, Thase ME. Treatment of SSRI-resistant depression: a meta-analysis comparing within- versus across-class switches. Biol Psychiatry (2008) 63(7):699–704. doi:10.1016/j.biopsych.2007.08.010
9. Baune BT, Hohoff C, Berger K, Neumann A, Mortensen S, Roehrs T, et al. Association of the COMT val158met variant with antidepressant treatment response in major depression. Neuropsychopharmacology (2008) 33(4):924–32. doi:10.1038/sj.npp.1301462
10. Horstmann S, Lucae S, Menke A, Hennings JM, Ising M, Roeske D, et al. Polymorphisms in GRIK4, HTR2A, and FKBP5 show interactive effects in predicting remission to antidepressant treatment. Neuropsychopharmacology (2010) 35(3):727–40. doi:10.1038/npp.2009.180
11. Domschke K, Dannlowski U, Ohrmann P, Lawford B, Bauer J, Kugel H, et al. Cannabinoid receptor 1 (CNR1) gene: impact on antidepressant treatment response and emotion processing in major depression. Eur Neuropsychopharmacol (2008) 18(10):751–9. doi:10.1016/j.euroneuro.2008.05.003
12. Baune BT, Hohoff C, Mortensen LS, Deckert J, Arolt V, Domschke K. Serotonin transporter polymorphism (5-HTTLPR) association with melancholic depression: a female specific effect? Depress Anxiety (2008) 25(11):920–5. doi:10.1002/da.20433
13. Baffa A, Hohoff C, Baune BT, Muller-Tidow C, Tidow N, Freitag C, et al. Norepinephrine and serotonin transporter genes: impact on treatment response in depression. Neuropsychobiology (2010) 62(2):121–31. doi:10.1159/000317285
14. Domschke K, Hohoff C, Mortensen LS, Roehrs T, Deckert J, Arolt V, et al. Monoamine oxidase A variant influences antidepressant treatment response in female patients with major depression. Prog Neuropsychopharmacol Biol Psychiatry (2008) 32(1):224–8. doi:10.1016/j.pnpbp.2007.08.011
15. Baune BT, Dannlowski U, Domschke K, Janssen DG, Jordan MA, Ohrmann P, et al. The interleukin 1 beta (IL1B) gene is associated with failure to achieve remission and impaired emotion processing in major depression. Biol Psychiatry (2010) 67(6):543–9. doi:10.1016/j.biopsych.2009.11.004
16. Amare AT, Schubert KO, Baune BT. Pharmacogenomics in the treatment of mood disorders: strategies and opportunities for personalized psychiatry. EPMA J (2017) 8(3):211–27. doi:10.1007/s13167-017-0112-8
17. Direk N, Williams S, Smith JA, Ripke S, Air T, Amare AT, et al. An analysis of two genome-wide association meta-analyses identifies a new locus for broad depression phenotype. Biol Psychiatry (2017) 82(5):322–9. doi:10.1016/j.biopsych.2016.11.013
18. Tansey KE, Guipponi M, Hu X, Domenici E, Lewis G, Malafosse A, et al. Contribution of common genetic variants to antidepressant response. Biol Psychiatry (2013) 73(7):679–82. doi:10.1016/j.biopsych.2012.10.030
19. Iniesta R, Malki K, Maier W, Rietschel M, Mors O, Hauser J, et al. Combining clinical variables to optimize prediction of antidepressant treatment outcomes. J Psychiatr Res (2016) 78:94–102. doi:10.1016/j.jpsychires.2016.03.016
20. Smagula SF, Butters MA, Anderson SJ, Lenze EJ, Dew MA, Mulsant BH, et al. Antidepressant response trajectories and associated clinical prognostic factors among older adults. JAMA Psychiatry (2015) 72(10):1021–8. doi:10.1001/jamapsychiatry.2015.1324
22. Caudill MM, Hunter AM, Cook IA, Leuchter AF. The antidepressant treatment response index as a predictor of reboxetine treatment outcome in major depressive disorder. Clin EEG Neurosci (2015) 46(4):277–84. doi:10.1177/1550059414532443
23. De Carlo V, Calati R, Serretti A. Socio-demographic and clinical predictors of non-response/non-remission in treatment resistant depressed patients: a systematic review. Psychiatry Res (2016) 240:421–30. doi:10.1016/j.psychres.2016.04.034
25. Quilty LC, De Fruyt F, Rolland JP, Kennedy SH, Rouillon PF, Bagby RM. Dimensional personality traits and treatment outcome in patients with major depressive disorder. J Affect Disord (2008) 108(3):241–50. doi:10.1016/j.jad.2007.10.022
26. Takahashi M, Shirayama Y, Muneoka K, Suzuki M, Sato K, Hashimoto K. Low openness on the revised NEO personality inventory as a risk factor for treatment-resistant depression. PLoS One (2013) 8(9):e71964. doi:10.1371/journal.pone.0071964
27. Gorwood P, Rouillon F, Even C, Falissard B, Corruble E, Moran P. Treatment response in major depression: effects of personality dysfunction and prior depression. Br J Psychiatry (2010) 196(2):139–42. doi:10.1192/bjp.bp.109.067058
29. Tang TZ, DeRubeis RJ, Hollon SD, Amsterdam J, Shelton R, Schalet B. Personality change during depression treatment: a placebo-controlled trial. Arch Gen Psychiatry (2009) 66(12):1322–30. doi:10.1001/archgenpsychiatry.2009.166
32. Du L, Bakish D, Ravindran AV, Hrdina PD. Does fluoxetine influence major depression by modifying five-factor personality traits? J Affect Disord (2002) 71(1–3):235–41. doi:10.1016/S0165-0327(01)00370-6
33. Jylhä P, Ketokivi M, Mantere O, Melartin T, Holma M, Rytsälä H, et al. Do antidepressants change personality?—A five-year observational study. J Affect Disord (2012) 142(1):200–7. doi:10.1016/j.jad.2012.04.026
34. Lo M-T, Hinds DA, Tung JY, Franz C, Fan C-C, Wang Y, et al. Genome-wide analyses for personality traits identify six genomic loci and show correlations with psychiatric disorders. Nat Genet (2017) 49(1):152–6. doi:10.1038/ng.3736
35. Okbay A, Baselmans BML, De Neve J-E, Turley P, Nivard MG, Fontana MA, et al. Genetic variants associated with subjective well-being, depressive symptoms, and neuroticism identified through genome-wide analyses. Nat Genet (2016) 48(6):624–33. doi:10.1038/ng.3552
36. Kendler KS, Neale MC, Kessler RC, Heath AC, Eaves LJ. A longitudinal twin study of personality and major depression in women. Arch Gen Psychiatry (1993) 50(11):853–62. doi:10.1001/archpsyc.1993.01820230009001
37. de Moor MH, Costa PT, Terracciano A, Krueger RF, de Geus EJ, Toshiko T, et al. Meta-analysis of genome-wide association studies for personality. Mol Psychiatry (2012) 17(3):337–49. doi:10.1038/mp.2010.128
38. International Schizophrenia Consortium, Purcell SM, Wray NR, Stone JL, Visscher PM, O’Donovan MC, et al. Common polygenic variation contributes to risk of schizophrenia that overlaps with bipolar disorder. Nature (2009) 460(7256):748–52. doi:10.1038/nature08185
39. Amare AT, Schubert KO, Hou L, Clark SR, Papiol S, Heilbronner U, et al. Association of polygenic score for schizophrenia and HLA antigen and inflammation genes with response to lithium in bipolar affective disorder: a genome-wide association study. JAMA Psychiatry (2018) 75(1):65–74. doi:10.1001/jamapsychiatry.2017.3433
40. Mrazek DA, Biernacka JM, McAlpine DE, Benitez J, Karpyak VM, Williams MD, et al. Treatment outcomes of depression: the pharmacogenomic research network antidepressant medication pharmacogenomic study. J Clin Psychopharmacol (2014) 34(3):313–7. doi:10.1097/JCP.0000000000000099
41. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet (2007) 81(3):559–75. doi:10.1086/519795
43. Fava M, Rush AJ, Trivedi MH, Nierenberg AA, Thase ME, Sackeim HA, et al. Background and rationale for the sequenced treatment alternatives to relieve depression (STAR*D) study. Psychiatr Clin North Am (2003) 26(2):457–94,x. doi:10.1016/S0193-953X(02)00107-7
47. Henkel V, Seemuller F, Obermeier M, Adli M, Bauer M, Mundt C, et al. Does early improvement triggered by antidepressants predict response/remission? Analysis of data from a naturalistic study on a large sample of inpatients with major depression. J Affect Disord (2009) 115(3):439–49. doi:10.1016/j.jad.2008.10.011
48. Seemuller F, Riedel M, Obermeier M, Bauer M, Adli M, Kronmuller K, et al. Outcomes of 1014 naturalistically treated inpatients with major depressive episode. Eur Neuropsychopharmacol (2010) 20(5):346–55. doi:10.1016/j.euroneuro.2009.11.011
49. Uher R, Mors O, Rietschel M, Rajewska-Rager A, Petrovic A, Zobel A, et al. Early and delayed onset of response to antidepressants in individual trajectories of change during treatment of major depression: a secondary analysis of data from the Genome-Based Therapeutic Drugs for Depression (GENDEP) study. J Clin Psychiatry (2011) 72(11):1478–84. doi:10.4088/JCP.10m06419
50. Harmer CJ, Goodwin GM, Cowen PJ. Why do antidepressants take so long to work? A cognitive neuropsychological model of antidepressant drug action. Br J Psychiatry (2009) 195(2):102–8. doi:10.1192/bjp.bp.108.051193
51. Godlewska BR, Norbury R, Selvaraj S, Cowen PJ, Harmer CJ. Short-term SSRI treatment normalises amygdala hyperactivity in depressed patients. Psychol Med (2012) 42(12):2609–17. doi:10.1017/S0033291712000591
52. Di Simplicio M, Norbury R, Reinecke A, Harmer CJ. Paradoxical effects of short-term antidepressant treatment in fMRI emotional processing models in volunteers with high neuroticism. Psychol Med (2014) 44(2):241–52. doi:10.1017/S0033291713000731
53. Bagby RM, Quilty LC, Segal ZV, McBride CC, Kennedy SH, Costa PT. Personality and differential treatment response in major depression: a randomized controlled trial comparing cognitive-behavioural therapy and pharmacotherapy. Can J Psychiatry (2008) 53(6):361–70. doi:10.1177/070674370805300605
55. de Moor MH, van den Berg SM, Verweij KJ, Krueger RF, Luciano M, Arias Vasquez A, et al. Meta-analysis of genome-wide association studies for neuroticism, and the polygenic association with major depressive disorder. JAMA Psychiatry (2015) 72(7):642–50. doi:10.1001/jamapsychiatry.2015.0554
56. Yamagata H, Uchida S, Matsuo K, Harada K, Kobayashi A, Nakashima M, et al. Identification of commonly altered genes between in major depressive disorder and a mouse model of depression. Sci Rep (2017) 7(1):3044. doi:10.1038/s41598-017-03291-x
57. O’Leary NA, Wright MW, Brister JR, Ciufo S, Haddad D, McVeigh R, et al. Reference sequence (RefSeq) database at NCBI: current status, taxonomic expansion, and functional annotation. Nucleic Acids Res (2016) 44(D1):D733–45. doi:10.1093/nar/gkv1189
58. Dong C, Wong ML, Licinio J. Sequence variations of ABCB1, SLC6A2, SLC6A3, SLC6A4, CREB1, CRHR1 and NTRK2: association with major depression and antidepressant response in Mexican-Americans. Mol Psychiatry (2009) 14(12):1105–18. doi:10.1038/mp.2009.92
59. DeYoung CG, Cicchetti D, Rogosch FA. Moderation of the association between childhood maltreatment and neuroticism by the corticotropin-releasing hormone receptor 1 gene. J Child Psychol Psychiatry (2011) 52(8):898–906. doi:10.1111/j.1469-7610.2011.02404.x
60. Amare AT, Schubert KO, Klingler-Hoffmann M, Cohen-Woods S, Baune BT. The genetic overlap between mood disorders and cardiometabolic diseases: a systematic review of genome wide and candidate gene studies. Transl Psychiatry (2017) 7:e1007. doi:10.1038/tp.2016.261
61. Pruitt KD, Brown GR, Hiatt SM, Thibaud-Nissen F, Astashyn A, Ermolaeva O, et al. RefSeq: an update on mammalian reference sequences. Nucleic Acids Res (2014) 42(Database issue):D756–63. doi:10.1093/nar/gkt1114
62. Walss-Bass C, Soto-Bernardini MC, Johnson-Pais T, Leach RJ, Ontiveros A, Nicolini H, et al. Methionine sulfoxide reductase: a novel schizophrenia candidate gene. Am J Med Genet B Neuropsychiatr Genet (2009) 150b(2):219–25. doi:10.1002/ajmg.b.30791
63. Ma X, Deng W, Liu X, Li M, Chen Z, He Z, et al. A genome-wide association study for quantitative traits in schizophrenia in China. Genes Brain Behav (2011) 10(7):734–9. doi:10.1111/j.1601-183X.2011.00712.x
64. Ni P, Ma X, Lin Y, Lao G, Hao X, Guan L, et al. Methionine sulfoxide reductase A (MsrA) associated with bipolar I disorder and executive functions in A Han Chinese population. J Affect Disord (2015) 184:235–8. doi:10.1016/j.jad.2015.06.004
65. Day FR, Helgason H, Chasman DI, Rose LM, Loh P-R, Scott RA, et al. Physical and neuro-behavioural determinants of reproductive onset and success. Nat Genet (2016) 48(6):617–23. doi:10.1038/ng.3551
66. Li J, Ma H, Deng S, Wu L, Huang Y, Zhu G. Fyn polymorphisms are associated with distinct personality traits in healthy Chinese-Han subjects. J Mol Neurosci (2011) 44(1):1–5. doi:10.1007/s12031-010-9485-z
Keywords: pharmacogenomics, polygenic score, personality traits, major depression, antidepressants, selective serotonin reuptake inhibitors
Citation: Amare AT, Schubert KO, Tekola-Ayele F, Hsu Y-H, Sangkuhl K, Jenkins G, Whaley RM, Barman P, Batzler A, Altman RB, Arolt V, Brockmöller J, Chen C-H, Domschke K, Hall-Flavin DK, Hong C-J, Illi A, Ji Y, Kampman O, Kinoshita T, Leinonen E, Liou Y-J, Mushiroda T, Nonen S, Skime MK, Wang L, Kato M, Liu Y-L, Praphanphoj V, Stingl JC, Bobo WV, Tsai S-J, Kubo M, Klein TE, Weinshilboum RM, Biernacka JM and Baune BT (2018) Association of the Polygenic Scores for Personality Traits and Response to Selective Serotonin Reuptake Inhibitors in Patients with Major Depressive Disorder. Front. Psychiatry 9:65. doi: 10.3389/fpsyt.2018.00065
Received: 08 December 2017; Accepted: 19 February 2018;
Published: 06 March 2018
Edited by:Stefan Borgwardt, University of Basel, Switzerland
Reviewed by:Ju Wang, Tianjin Medical University, China
Kurt Leroy Hoffman, Autonomous University of Tlaxcala, Mexico
Copyright: © 2018 Amare, Schubert, Tekola-Ayele, Hsu, Sangkuhl, Jenkins, Whaley, Barman, Batzler, Altman, Arolt, Brockmöller, Chen, Domschke, Hall-Flavin, Hong, Illi, Ji, Kampman, Kinoshita, Leinonen, Liou, Mushiroda, Nonen, Skime, Wang, Kato, Liu, Praphanphoj, Stingl, Bobo, Tsai, Kubo, Klein, Weinshilboum, Biernacka and Baune. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Bernhard T. Baune, email@example.com