Personality Polygenes, Positive Affect, and Life Satisfaction

HHS Public Access Author manuscript Author Manuscript Twin Res Hum Genet. Author manuscript; available in PMC 2017 October 01. Published in final edited form as: Twin Res Hum Genet. 2016 October ; 19(5): 407–417. doi:10.1017/thg.2016.65. Personality Polygenes, Positive Affect, and Life Satisfaction Author Manuscript Alexander Weiss1, Bart M. L. Baselmans2,3, Edith Hofer4,5, Jingyun Yang6,7, Aysu Okbay8,9,10, Penelope A. Lind11, Mike B. Miller12, Ilja M. Nolte13, Wei Zhao14, Saskia P. Hagenaars1,15, Jouke-Jan Hottenga2,3, Lindsay K. Matteson12, Harold Snieder13, Jessica D. Faul16, Catharina A. Hartman17, Patricia A. Boyle6,18, Henning Tiemeier9,19,20, Miriam A. Mosing21,22, Alison Pattie1, Gail Davies1, David C. Liewald1, Reinhold Schmidt4, Philip L. De Jager23,24,25, Andrew C. Heath26, Markus Jokela27, John M. Starr28, Albertine J. Oldehinkel17, Magnus Johannesson29, David Cesarini30,31, Albert Hofman9, Sarah E. Harris1,32, Jennifer A. Smith14, Liisa Keltikangas-Järvinen33, Laura Pulkki-Råback33,34, Helena Schmidt35,36, Jacqui Smith37, William G. Iacono12, Matt McGue12, David A. Bennett6,7, Nancy L. Pedersen22, Patrik K. E. Magnusson22, Ian J. Deary1, Nicholas G. Martin11, Dorret I. Boomsma2,3,38, Meike Bartels2,3,38,†, and Michelle Luciano1,† 1Centre Author Manuscript for Cognitive Ageing and Cognitive Epidemiology, Department of Psychology, School of Philosophy, Psychology and Language Sciences, The University of Edinburgh, Edinburgh, UK 2Department of Biological Psychology, Netherlands Twin Register, VU University Amsterdam, Amsterdam, the Netherlands 3EMGO+ Institute for Health and Care Research, VU University Medical Centre, Amsterdam, the Netherlands 4Clinical Division of Neurogeriatrics, Department of Neurology, Medical University Graz, Austria 5Institute of Medical Informatics, Statistics and Documentation, Medical University Graz, Austria 6Rush Alzheimer’s Disease Center, Rush University Medical Center, Chicago, IL, USA 7Department of Neurological Sciences, Rush University Medical Center, Chicago, IL, USA 8Department of Applied Economics, Erasmus School of Economics, Erasmus University Rotterdam, Rotterdam, the Netherlands 9Department of Epidemiology, Erasmus Medical Center, Rotterdam, the Netherlands 10Erasmus University Rotterdam Institute for Behavior and Biology, Rotterdam, the Netherlands 11Quantitative Genetics, QIMR Berghofer Institute of Medical Research, Brisbane, Queensland, Australia 12Department of Psychology, University of Minnesota, USA 13Department of Epidemiology, University of Groningen, Groningen, the Netherlands 14Survey Research Center, Institute for Social Research, Author Manuscript This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/ licenses/by/4.0/), which permits unrestricted re-use, distribution, and reproduction in any medium, provided the original work is properly cited. Address for correspondence: Michelle Luciano, Centre for Cognitive Ageing and Cognitive Epidemiology, Department of Psychology, School of Philosophy, Psychology and Language Sciences, The University of Edinburgh, 7 George Square, EH8 9JZ Edinburgh, UK. michelle.luciano@ed.ac.uk. †Shared authorship. *See online supplementary Lifelines author list. Conflict of Interest None. Ethical Standards The authors assert that all procedures contributing to this work comply with the ethical standards of the relevant national and institutional committees on human experimentation and with the Helsinki Declaration of 1975, as revised in 2008. Supplementary Material To view supplementary material for this article, please visit http://dx.doi.org/10.1017/thg.2016.65. Weiss et al. Page 2 Author Manuscript Author Manuscript Author Manuscript University of Michigan, Ann Arbor, MI, USA 15Division of Psychiatry, University of Edinburgh, Royal Edinburgh Hospital, Edinburgh, UK 16Department of Epidemiology, School of Public Health, University of Michigan, Ann Arbor, MI, USA 17Interdisciplinary Center Psychopathology and Emotion regulation, University Medical Center, University of Groningen, Groningen, the Netherlands 18Department of Behavioral Sciences, Rush University Medical Center, Chicago, Illinois, USA 19Department of Psychiatry, Erasmus Medical Center, Rotterdam, the Netherlands 20Department of Child and Adolescent Psychiatry, Erasmus Medical Center, Rotterdam, the Netherlands 21Department of Neuroscience, Karolinska Institutet, Stockholm, Sweden 22Department of Medical Epidemiology and Biostatistics, Karolinska Institutet, Stockholm, Sweden 23Program in Translational NeuroPsychiatric Genomics, Institute for the Neurosciences, Departments of Neurology and Psychiatry, Brigham and Women’s Hospital, Boston, MA, USA 24Harvard Medical School, Boston, MA, USA 25Program in Medical and Population Genetics, Broad Institute, Cambridge, MA, USA 26Division of Biology and Biomedical Sciences, Washington University, MO, USA 27Institute of Behavioural Sciences, University of Helsinki, Finland 28Geriatric Medicine Unit, Western General Hospital, Edinburgh, and Centre for Cognitive Ageing and Cognitive Epidemiology, University of Edinburgh, UK 29Department of Economics, Stockholm School of Economics, Stockholm, Sweden 30Department of Economics, New York University, New York, USA 31Research Institute for Industrial Economics, Stockholm, Sweden 32Medical Genetics Section, University of Edinburgh Centre for Genomic and Experimental Medicine and MRC Institute of Genetics and Molecular Medicine, Western General Hospital, Crewe Road, Edinburgh, UK 33IBS, Unit of Personality, Work and Health, Institute of Behavioural Sciences, University of Helsinki, Finland 34Helsinki Collegium for Advanced Studies, University of Helsinki, Finland 35Department of Neurology, Medical University Graz, Austria 36Institute of Molecular Biology and Biochemistry, Centre for Molecular Medicine, Medical University of Graz, Graz, Austria 37Department of Psychology, University of Michigan, Ann Arbor, MI, USA 38Neuroscience Campus Amsterdam, Amsterdam, the Netherlands Abstract Author Manuscript Approximately half of the variation in wellbeing measures overlaps with variation in personality traits. Studies of non-human primate pedigrees and human twins suggest that this is due to common genetic influences. We tested whether personality polygenic scores for the NEO FiveFactor Inventory (NEO-FFI) domains and for item response theory (IRT) derived extraversion and neuroticism scores predict variance in wellbeing measures. Polygenic scores were based on published genome-wide association (GWA) results in over 17,000 individuals for the NEO-FFI and in over 63,000 for the IRT extraversion and neuroticism traits. The NEO-FFI polygenic scores were used to predict life satisfaction in 7 cohorts, positive affect in 12 cohorts, and general wellbeing in 1 cohort (maximal N = 46,508). Meta-analysis of these results showed no significant association between NEO-FFI personality polygenic scores and the wellbeing measures. IRT extraversion and neuroticism polygenic scores were used to predict life satisfaction and positive affect in almost 37,000 individuals from UK Biobank. Significant positive associations (effect sizes <0.05%) were observed between the extraversion polygenic score and wellbeing measures, and a negative association was observed between the polygenic neuroticism score and life satisfaction. Furthermore, using GWA data, genetic correlations of −0.49 and −0.55 were Twin Res Hum Genet. Author manuscript; available in PMC 2017 October 01. Weiss et al. Page 3 Author Manuscript estimated between neuroticism with life satisfaction and positive affect, respectively. The moderate genetic correlation between neuroticism and wellbeing is in line with twin research showing that genetic influences on wellbeing are also shared with other independent personality domains. Keywords wellbeing; genetics; polygenic prediction; happiness; genetic correlation Author Manuscript Happiness is a desirable state that is universally pursued. It is also linked to personality traits, such as those of the Five-Factor Model (Adams et al., 2012; DeNeve & Cooper, 1998). Individuals who score lower on neuroticism and higher on extraversion, agreeableness, and conscientiousness report being happier and more satisfied with their lives (meta-analytic correlations ranged 0.17–0.22; DeNeve & Cooper, 1998). Genetic influences account for approximately 40% of variation in wellbeing (Bartels, 2015), which is comparable to the heritability estimates for personality traits (Bouchard & Loehlin, 2001). Genetic analysis has shown that although unique, non-additive genetic effects were found for happiness and general quality of life (Bartels & Boomsma, 2009), a common additive genetic factor influences different well-being measures (i.e., general quality of life, present quality of life, life satisfaction, and subjective happiness/positive affect). Author Manuscript Evidence for shared genetic variance between personality and wellbeing comes from biometric genetic studies of great ape pedigrees (Adams et al., 2012; Weiss et al., 2002). It also comes from studies of human twins and siblings. Using a three-item wellbeing measure (present and general life satisfaction, control over one’s life), Weiss et al. (2008) showed that a general personality additive genetic factor explained 2.2% of the variance in wellbeing. Additional genetic contributions to wellbeing were via independent factors that influenced neuroticism (5.3% of variance), extraversion (13%), and conscientiousness (0.8%). Hahn et al. (2013) confirmed the absence of unique genes influencing a multidimensional measure of life satisfaction in their extended twin study, additionally showing shared non-additive genetic variance between neuroticism and life satisfaction. Author Manuscript A complementary test of the hypothesis that common genes underlie variation in personality and happiness is to use molecular data, such as single nucleotide polymorphisms (SNPs). In a recent large study (N ≈ 300K), a polygenic score constructed from a genome-wide association (GWA) meta-analysis on subjective wellbeing explained ~0.7% of the variance in neuroticism and ∼0.4% of the variance in extraversion (Okbay et al., 2016). Applying bivariate linkage disequilibrium score regression (Bulik-Sullivan et al., 2015) to the GWA summary statistics for wellbeing and neuroticism resulted in a SNP-based genetic correlation of −0.75 (SE = 0.034; Okbay et al., 2016). This genetic correlation represents the correlation of common, additive genetic effects between the two traits. Whereas the variance in a trait explained by polygenic scores is typically low, methods to infer the expected SNPderived variance from polygenic scores show agreement with their empirical and simulationbased estimates (Dudbridge, 2013; Wray et al., 2014). Twin Res Hum Genet. Author manuscript; available in PMC 2017 October 01. Weiss et al. Page 4 Author Manuscript To provide greater support for a genetic association between personality and wellbeing, our aim here is to predict phenotypic scores for wellbeing and its subcomponents of life satisfaction and positive affect by using information about SNP effects on neuroticism, extraversion, openness, agreeableness, and conscientiousness. We used a method involving polygenic prediction models that enabled us to test whether genes influencing one trait influence another trait (for a review, see Wray et al., 2014). In this method, GWA results of a trait are used to create a polygenic score representing the sum of the effects of individual SNPs on that trait in an independent sample. This score is then used to predict the trait of interest. Polygenic prediction models do not require family designs, enabling the use of a large number of population-based studies with wellbeing and genotyping data. Author Manuscript We furthermore established genetic correlations between neuroticism and wellbeing measures by using a bivariate restricted maximum likelihood (REML) estimation (Lee et al., 2012) that has not previously been applied to these traits. This method uses genome-wide SNP data to calculate a genetic relationship matrix between unrelated individuals which within a REML framework allows estimation of the heritability due to all SNPs. This extends to the bivariate case from which genetic correlations can be ascertained. Author Manuscript We created polygenic scores using GWA results for the NEO Five-Factor Inventory (NEOFFI; de Moor et al., 2012) and for extraversion and neuroticism from item response theory (IRT) analyses of varying personality scales (de Moor et al., 2015; van den Berg et al., 2015). Whereas the NEO-FFI GWA meta-analysis comprised a smaller total sample size (N = 17,375) than the IRT extraversion and neuroticism GWA meta-analyses (N ~ 63,000), importantly, it measures all five personality domains, and polygenic prediction based on these results has been successful for extraversion (predicting bipolar disorder) and neuroticism (predicting major depressive disorder; Middeldorp et al., 2011). We used unitweighted tests to determine whether the polygenic score of any personality domain was associated with phenotypic variance in life satisfaction, positive affect, and wellbeing. For the NEO-FFI GWA results, polygenic prediction was tested in 14 cohorts that were independent of the GWA, and for the IRT extraversion and neuroticism GWA results, polygenic prediction was tested in the UK Biobank, which was independent of the GWA meta-analyses. To establish genetic correlations between neuroticism and wellbeing using bivariate REML, we used a large cohort of unrelated individuals with genome-wide data and measurements on all the traits of interest. Methods Participants Author Manuscript NEO-FFI polygenic prediction in 14 cohorts—Cohorts were drawn from a GWA study meta-analysis of wellbeing conducted by the Social Sciences Genetic Association Consortium (SSGAC; http://www.thessgac.org), with the proviso that none of the cohorts were part of the GWAS meta-analysis of the NEO-FFI (de Moor et al., 2012); personality data were not required for analysis. Participants were (or were ancestors of) white Europeans. Thirteen cohorts with positive affect (n ranged 351–11,971) and seven cohorts with life satisfaction (n ranged 351–9,938) were available (five cohorts had positive affect and life satisfaction measures) for inclusion in our meta-analysis. An additional cohort (n = Twin Res Hum Genet. Author manuscript; available in PMC 2017 October 01. Weiss et al. Page 5 Author Manuscript 6,960) had a measure of general wellbeing that was analyzed separately. Individual cohort descriptions, including the scales and/or items used to measure wellbeing, are provided in the Supplementary information. The relevant institutional ethics review boards approved the individual studies. DNA was extracted using standard protocols. Genotyping procedures are summarized in Supplementary Table S1. Cohorts used HapMap II imputed data or, if unavailable, observed genotypes for analysis. Imputed data were preferred because the GWAS personality results were based on HapMap II data, thus ensuring that all SNPs would be matched to those available in the GWAS. One cohort used 1000G imputation but removed SNPs that were not available in HapMap II. Author Manuscript Author Manuscript Extraversion and neuroticism polygenic prediction in UK Biobank—Five of the SSGAC cohorts participated in the IRT extraversion and neuroticism GWA studies (de Moor et al., 2015; van den Berg et al., 2015); therefore, another independent cohort was sought for this prediction analysis. Participants were drawn from the baseline survey of the UK Biobank (http://www.ukbiobank.ac.uk), a resource established for investigating factors influencing disease in middle and older age. These measures (including questionnaire and biological samples) were collected between 2006 and 2010 on 502,655 British community residing individuals, a subset of whom were used in the present study. Positive affect was measured by the item ‘In general how happy are you?’ on a six-point scale (extremely happy, very happy, moderately happy, moderately unhappy, very unhappy, and extremely unhappy). General life satisfaction was surveyed across family relationships, financial situation, friendship, health, and work/job domains on the same six-point scale. Responses on these items demonstrated positive manifold and were best described by a single factor that explained 37% of variance. An averaged life satisfaction score was used to account for missing data where a person was currently unemployed (n = 11,679), did not know (n ranged 97–380), or preferred not to answer (n ranged 30–170). Neuroticism was measured by 12 items from the Eysenck Personality Questionnaire Revised (Eysenck & Eysenck, 1991). Wellbeing data were available for 36,737 (positive affect) and 36,911 (life satisfaction) individuals with genome-wide genotyping data. These data were skewed in the direction of lower positive affect/life satisfaction, but no ceiling effect was present. Ages ranged between 40 and 70 years (mean age = 57.31 years, SD = 7.92). Author Manuscript DNA was obtained via blood samples and genotyping performed with either the UK BiLEVE array or the UK Biobank axiom array. Standard quality control procedures were followed, including checks for gender mismatch and non-British ancestry. Further description can be found in Hagenaars et al. (2016). Polygenic scores were created on the observed genotypes. UK Biobank received ethical approval from the Research Ethics Committee (REC reference 11/NW/0382). Statistical Analysis NEO-FFI polygenic prediction in 14 cohorts—Five sets of polygenic scores representing the personality domains of neuroticism, extraversion, openness, agreeableness, and conscientiousness were estimated using SNP association results from the largest GWA Twin Res Hum Genet. Author manuscript; available in PMC 2017 October 01. Weiss et al. Page 6 Author Manuscript meta-analysis of NEO-FFI domains to date (de Moor et al., 2012). This GWA study included 10 discovery samples (N = 17,375). None of the cohorts – except NTR – in the present study were part of this personality GWA. For their analyses, NTR removed the participants who were part of the personality GWA meta-analysis. Author Manuscript Personality polygenic scores were estimated in each cohort using five probability thresholds for choosing SNPs to include in the score. These were based on the significance value for each SNP from the GWAS meta-analysis: p < .01, p < .05, p < .1, p < .5, and p < 1. Polygenic scores were formed by summing the meta-analytic effect size coefficients (betas) weighted by the number of copies (0/1/2) of the effect allele carried by the individual across all SNPs within the threshold. For imputed data, best guess genotypes were used but excluding SNPs with an imputation quality estimated r2 less than 0.80. Before score calculation, SNPs with a minor allele frequency <0.05 and Hardy–Weinberg Equilibrium test <p × 10−7 were removed. SNPs were then pruned for linkage disequilibrium using an r2 cutoff of 0.25 within a 200-SNP sliding window, following Purcell et al. (2009). Missing SNPs for an individual were imputed dependent on the observed allele frequency in the cohort. Polygenic scores were calculated using PLINK (Purcell et al., 2007). Supplementary Table S2 shows the number of SNPs included in the calculation of the polygenic score at each threshold for all the cohorts. Author Manuscript Author Manuscript To predict phenotypic wellbeing scores from the polygenic personality scores, regression analysis was used. The dependent measures (positive affect, life satisfaction, and general wellbeing) were residualized on age, age squared (if significant), sex, population stratification components, and number of non-missing SNPs contributing to each individual’s score (where observed genotypes were used or where sparse genotyping led to poorer imputation quality). Standardized residual scores were then used as the dependent variable. A series of univariate regression analyses using each of the five polygenic personality scores as predictors was run for each polygenic score threshold (i.e., 25 tests). For the MCTFR cohort, a feasible generalized least squares regression was used to account for familial correlations. For NTR, a generalized estimating equating model was used to account for family structure. A meta-analysis of the standardized regression coefficients from the regression models for life satisfaction and positive affect was performed assuming random effects in R (MAc package; http://cran.r-project.org/web/packages/MAc/ index.html). This produced an overall effect size and standard error. A false discovery rate correction (Benjamini and Hochberg method) to an alpha level of 0.05 was applied to each of the meta-analyses and to the analysis of general wellbeing. Cohort estimate heterogeneity was assessed by Cochran’s Q, which uses the sum of squared deviations of each study’s effect size from the meta-analytic estimate to determine significance. A supplementary meta-analysis was performed on combined life satisfaction, positive affect and general wellbeing measures to obtain a maximal sample size (∼10,000 more individuals than the positive affect analysis). Where a cohort had two measures, the measure with the larger sample was chosen. Extraversion and neuroticism polygenic prediction in UK Biobank—Five polygenic scores were calculated for extraversion and neuroticism based on the significance value for each SNP from the largest respective GWA meta-analysis of these traits (p < .01, p Twin Res Hum Genet. Author manuscript; available in PMC 2017 October 01. Weiss et al. Page 7 Author Manuscript < .05, p < .1, p < .5, and p < 1 (de Moor et al., 2015; van den Berg et al., 2015). Both GWA studies were based on the same 29 meta-analysis samples that included 63,030 individuals for extraversion and 63,661 individuals for neuroticism. Because there was variation in the personality scale used across samples, an IRT procedure was used to harmonize the personality traits prior to GWA (van den Berg et al., 2014). Polygenic scores (as described in the previous section) for extraversion and neuroticism were created using PRSice software (Euesden et al., 2015) at the five SNP inclusion levels. Before calculating the scores, exclusions were made of SNPs with low minor allele frequency (<0.01) and of SNPs in linkage disequilibrium (r2 > 0.25) using a clumping method within a 250 kb window. A lower minor allele frequency level exclusion was set for this sample due to its much larger size than the samples comprising the meta-analysis described above; and given the increased reliability of individual effects from the larger GWA meta-analysis, the clumping procedure, which preferentially selects SNPs showing the greatest association, was preferred. For extraversion, the polygenic scores were the composite of 4,271, 18,606, 34,981, 143,525, and 238,487 SNPs for respective p < .01, p < .05, p < .1, p < .5, and p < 1 inclusion thresholds. For neuroticism, the polygenic scores were the composite of 4,266, 18,427, 34,700, 143,520, and 205,751 SNPs for respective p < .01, p < .05, p < .1, p < .5, and p < 1 inclusion thresholds. The regression models for polygenic extraversion and neuroticism scores predicting wellbeing included additional independent variables: age at survey, sex, genotyping batch and array, assessment center, and the first 10 genetic principal components (to correct for population stratification). FDR correction was applied to these analyses. Author Manuscript Author Manuscript Genetic correlations between neuroticism and wellbeing in UK Biobank— Given the large size of UK biobank and the availability of neuroticism and two wellbeing measures, genetic correlations were derived using SNP-based methods (bivariate REML; Lee et al., 2012). This method uses a standard bivariate linear model in which random polygenic effects are fitted and the variance covariance matrix conditioned by a genomic similarity relationship matrix that is estimated from genome-wide SNP information. The program GCTA (Yang et al., 2011) was used for this analysis on unrelated individuals only (individuals with a genetic similarity >0.025 were removed) to remove potential confounding from environmental influences. Observed genotypes were used excluding SNPs with a minor allele frequency less than 0.01. All phenotypes were regressed for the effects of age, sex, assessment center, genetic batch, genetic array, and 10 population stratification components; resulting residual scores were used in the GCTA analysis. Results NEO-FFI Polygenic Prediction in 14 Cohorts Author Manuscript Meta-analysis results for univariate regression models where personality polygenic scores predict life satisfaction and positive affect can be found in Tables 1 and 2, respectively. These tables display the regression beta, standard error and p value for each personality domain at each of the polygenic score inclusion thresholds (i.e., p < .01, p < .05, p < .1, p < . 5, and p < 1). Twin Res Hum Genet. Author manuscript; available in PMC 2017 October 01. Weiss et al. Page 8 Author Manuscript No tests were significant for life satisfaction or positive affect at the false discovery rate corrected alpha (q = 0.002). For positive affect, heterogeneity between cohorts was observed for all neuroticism polygenic scores, four of the extraversion polygenic score estimates and three of the agreeableness polygenic scores (see Supplementary Table S3, for individual cohort betas). The correlations between personality polygenic scores and wellbeing (and corresponding p values) are shown in Table 3. In this analysis, no correlations surpassed the FDR corrected significance level. Results from the meta-analysis in which all measures were combined are presented in Supplementary Table S4. No regression coefficients differed significantly from zero and there was significant heterogeneity between cohort estimates for five tests (neuroticism at SNP inclusion p < .01, extraversion at SNP inclusion p < .5 and p < 1, and agreeableness at SNP inclusion p < 0.5 and p < 1). IRT Extraversion and Neuroticism Polygenic Prediction in UK Biobank Author Manuscript The significance value and amount of variance explained by the polygenic extraversion and neuroticism scores in predicting life satisfaction and positive affect are shown in Figure 1. The FDR significance level was 0.0325. Extraversion polygenic scores significantly predicted both wellbeing measures (at all SNP inclusion thresholds for positive affect and at three thresholds for life satisfaction), whereas neuroticism polygenic scores significantly predicted only life satisfaction (at all thresholds). In all models, polygenic scores at the more liberal SNP inclusion thresholds explained more variance than the more restrictive SNP inclusion sets. The direction of the effect was as predicted with polygenic neuroticism scores negatively related to life satisfaction and extraversion positively related to measures of wellbeing. The amount of variance explained was extremely small, not exceeding 0.04%. Author Manuscript Genetic correlations between neuroticism and wellbeing in UK Biobank—For the analysis of neuroticism and positive affect, 30,367 individuals were included. SNP-based heritabilities of 0.15 (SE = 0.02) and 0.08 (SE = 0.02) were estimated for respective neuroticism and positive affect measures with a genetic correlation of −0.55 (SE = 0.09). The analysis of neuroticism and life satisfaction (N = 30,494) gave a heritability of 0.13 (SE = 0.02) for life satisfaction and a genetic correlation of −0.49 (SE = 0.07) with neuroticism. Discussion Author Manuscript These results build upon biometric research showing that common genes influence personality and happiness. The polygenic prediction based on the larger GWA of IRT-based extraversion and neuroticism showed significant association with wellbeing measures at a corrected false discovery rate. The personality polygenic prediction of wellbeing based on the smaller GWA of personality was non-significant for all five NEO-FFI domains. In the NEO-FFI meta-analysis heterogeneity was evident in, at most, four cohorts, suggesting that there were few differences owing to study specific factors (e.g., variation in measurement instrument). Because the meta-analysis and UK Biobank prediction samples were of comparable size (and resulting power), the limiting factor then for these analyses was the difference in power between the GWA studies of the NEO-FFI traits and IRT-based extraversion and neuroticism, on which the polygenic scores were based. In our test of the genetic correlation between neuroticism and wellbeing measures using genetic relationships Twin Res Hum Genet. Author manuscript; available in PMC 2017 October 01. Weiss et al. Page 9 Author Manuscript based on genome-wide SNP data, we found a moderate degree of genetic overlap for both positive affect and life satisfaction. Author Manuscript The finding in UK Biobank that extraversion polygenic scores predicted both life satisfaction and positive affect (measures showing a 0.62 phenotypic correlation in our sample) but that neuroticism polygenic scores predicted only life satisfaction was unexpected given that the combined measure of happiness and satisfaction with life used in the recent GWA of wellbeing significantly predicted neuroticism and extraversion (Okbay et al., 2016). Our finding in UK Biobank of similar-sized genetic correlations between neuroticism with positive affect and life satisfaction would also predict that polygenic neuroticism should relate to positive affect. The null finding might point to type 2 error rather than an interpretation that positive and negative affect are not opposite poles of the same dimension (e.g., Russell & Carroll, 1999). It is likely that the SNP-based genetic correlation between extraversion and positive affect will be stronger than for neuroticism, but we were unable to test this here because no other personality traits were collected in UK Biobank. Author Manuscript Author Manuscript The amount of variance in the wellbeing measures explained by the polygenic scores was extremely small, less than half a percent. But given that polygenic neuroticism only predicts 0.66% of variance in neuroticism itself (de Moor et al., 2015), our finding is not unexpected. As GWA meta-analysis studies of personality get larger, this effect size should increase; this is demonstrated by the superior reverse prediction of personality from polygenic wellbeing (Okbay et al., 2016). However, given the low estimated SNP-based heritabilities for neuroticism and wellbeing (<0.15 in our study), the limit for variance explained by a polygenic measure will necessarily be small. Twin and family studies show that heritabilities for personality and wellbeing are at least double that of the SNP-based estimates, which only consider the genetic variation due to common variants. Therefore, further gains in prediction might be achieved by investigating rare and/or structural genetic variants. There are no rare variant studies on personality, but in the only study (Power & Pluess, 2015) to estimate the heritability of all the Five-Factor Model domains using genome-wide SNP data (N = 5,011), only neuroticism and openness showed significant genetic influences, suggesting that rare variants might be important. With regard to structural variants, preliminary investigations do not show an effect of large copy number variants on personality (Luciano, MacLeod et al., 2012). Additionally, by using an additive composite of personality SNP effects we may have restricted the prediction of wellbeing. Extended twin studies show non-additive genetic effects for extraversion, neuroticism, and conscientiousness (Hahn et al., 2013; Keller et al., 2005), and measures of wellbeing (Bartels & Boomsma, 2009; Hahn et al., 2016). Further studies are therefore needed to confirm whether different personality traits share greater additive or non-additive genetic variance with wellbeing. Our study confirms that improvements in polygenic score prediction results from larger meta-analysis GWA studies of the predictor trait. However, it should be noted that Middeldorp et al. (2011) used a subsample (N = 13,835) of de Moor et al.’s (2012) NEO-FFI GWA study to create polygenic personality scores that predicted major depressive disorder (from neuroticism) and bipolar disorder (from extraversion). Moreover, Luciano, Huffman et al. (2012) predicted depressive symptoms from polygenic neuroticism using a GWA sample Twin Res Hum Genet. Author manuscript; available in PMC 2017 October 01. Weiss et al. Page 10 Author Manuscript Author Manuscript that was even smaller. Accepting that their results were not type 1 errors, one must ask why we failed to predict wellbeing here. One possibility is that the genetic correlations between neuroticism and extraversion are stronger with major depressive disorder (∼0.72; Middeldorp et al., 2005) and bipolar disorder (0.44; Hare et al., 2012) than the genetic correlations between personality and wellbeing (0.20–0.66; Weiss et al., 2008). These estimates, however, are based on twin studies where the similarity across all types of genetic variation is considered. The polygenic scores focus only on common variants, so genetic correlations based on these are more relevant. Using GWA results to estimate genetic correlations, neuroticism showed the same absolute correlation (0.75) with wellbeing (combined positive affect and life satisfaction) and depression (Okbay et al., 2016), although in our bivariate SNP-based method using raw genotypes, genetic correlations between neuroticism and separate positive affect and life satisfaction measures were lower (−0.55 and −0.49). Genetic correlations between extraversion and wellbeing using genome-wide SNP data will be informative. It may well be that personality has stronger genetic links with mental illness than wellbeing. That wellbeing is influenced predominantly by environmental factors unrelated to personality (Weiss et al., 2008) might also limit polygenic prediction. Author Manuscript Using the largest GWA studies to date of extraversion and neuroticism (independent of the UK Biobank sample) we confirmed that polygenic effects for these personality domains influenced wellbeing. Prediction tended to be better when using all SNP data rather than limiting prediction to a smaller number of SNPs with larger effects on personality. This suggests that many genes of very small effect are important for extraversion, neuroticism, and wellbeing. Although neuroticism has captured the interest of many researchers in cognitive psychology and psychiatry, our study also shows an important role of extraversion in mental wellbeing. We expect that genes influencing agreeableness, conscientiousness and openness will also have some role in explaining wellbeing, but our analysis could not reliably address this. Supplementary Material Refer to Web version on PubMed Central for supplementary material. Acknowledgments Author Manuscript We thank the study participants and the support staff involved in participant recruitment, phenotype collection, genotyping, data management and processing. Individual study acknowledgements can be found in the online Supplementary material. ASPS: The research reported in this article was funded by the Austrian Science Fond (FWF) grant numbers P20545-P05 and P13180. The Medical University of Graz supports the databank of the ASPS. HRS: HRS is supported by the National Institute on Aging (NIA U01AG009740). The genotyping was funded separately by the National Institute on Aging (RC2 AG036495, RC4 AG039029). LBC1921: Phenotype collection was supported by the UK’s Biotechnology and Biological Sciences Research Council (BBSRC), the Royal Society, and the Chief Scientist Office of the Scottish Government. Genotyping was funded by the BBSRC. The work was undertaken by the University of Edinburgh Centre for Cognitive Ageing and Cognitive Epidemiology, part of the cross council Lifelong Health and Wellbeing Initiative (MR/K026992/1). Funding from the BBSRC and Medical Research Council (MRC) is gratefully acknowledged. LifeLines Cohort Study: The LifeLines Cohort Study, and generation and management of GWAS genotype data for the LifeLines Cohort Study is supported by the Netherlands Organization of Scientific Research NWO (grant 175.010.2007.006), the Economic Structure Enhancing Fund (FES) of the Dutch government, the Ministry of Economic Affairs, the Ministry of Education, Culture and Science, the Ministry for Health, Welfare and Sports, the Northern Netherlands Collaboration of Provinces (SNN), the Province of Groningen, University Medical Center Groningen, the University of Groningen, Dutch Kidney Foundation and Dutch Diabetes Research Foundation. NTR: This work was supported by the Netherlands Organization for Scientific Research (NWO: MagW/ZonMW grants 904-61-090, Twin Res Hum Genet. Author manuscript; available in PMC 2017 October 01. Weiss et al. Page 11 Author Manuscript Author Manuscript Author Manuscript 985-10-002,904-61-193,480-04-004, 400-05-717, NOW-bilareral agreement 463-06-001, NWO-VENI 451-04-034, Addiction-31160008 Middelgroot-911-09-032, Spinozapremie 56-464-14192), Biobanking and Biomolecular Resources Research Infrastructure (BBMRI; NL, 184.021.007), the VU University’s Institute for Health and Care Research (EMGO+) and Neuroscience Campus Amsterdam (NCA), the European Science Council (ERC Advanced, 230374), the Avera Institute for Human Genetics, Sioux Falls, South Dakota (USA) and the National Institutes of Health (NIH, R01D0042157-01A). Part of the genotyping was funded by the Genetic Association Information Network (GAIN) of the Foundation for the US National Institutes of Health (NIMH, MH081802) and by the Grand Opportunity grants 1RC2MH089951-01 and 1RC2 MH089995-01 from the NIMH. Part of the analyses (for RS & TRAILS, too) were carried out on the Genetic Cluster Computer (http:// www.geneticcluster.org), which is financially supported by the Netherlands Scientific Organization (NWO 480-05-003), the Dutch Brain Foundation, and the Department of Psychology and Education of the VU University Amsterdam. M. Bartels is/was financially supported by a senior fellowship of the (EMGO+) Institute for Health and Care and a VU University Research Chair position. MCTFR: Research reported in this publication was supported by the National Institutes of Health under award numbers R37DA005147, R01AA009367, R01AA011886, R01DA013240, R01MH066140, and U01DA024417. QIMRB: Funding was provided by the Australian NHMRC (241944, 339462, 389927, 389875, 389891, 389892, 389938, 442915, 442981, 496739, 552485, 552498, 613608), the FP-5 GenomEUtwin Project (QLG2-CT-2002-01254), and the U.S. NIH grants (AA07535, AA10248, AA13320, AA13321, AA13326, AA14041, MH66206). RS: The generation and management of GWAS genotype data for the Rotterdam Study is supported by the Netherlands Organisation of Scientific Research NWO Investments (nr. 175.010.2005.011, 911-03-012). This study is funded by the Research Institute for Diseases in the Elderly (014-93-015; RIDE2), the Netherlands Genomics Initiative (NGI)/Netherlands Organisation for Scientific Research (NWO) project no. 050-060-810. The Rotterdam Study is funded by Erasmus Medical Center and Erasmus University, Rotterdam, Netherlands Organization for the Health Research and Development (ZonMw), the Research Institute for Diseases in the Elderly (RIDE), the Ministry of Education, Culture and Science, the Ministry for Health, Welfare and Sports, the European Commission (DG XII), and the Municipality of Rotterdam. STR: The Jan Wallander and Tom Hedelius Foundation (P2012-0002:1), the Ragnar Söderberg Foundation (E9/11), The Swedish Research Council (421-2013-1061), the Ministry for Higher Education, The Swedish Research Council (M-2205-1112), GenomEUtwin (EU/QLRT-2001-01254; QLG2-CT-2002-01254), NIH DK U01-066134, the Swedish Foundation for Strategic Research (SSF). The Rush Memory and Aging Project is supported by NIA grants R01AG15819 and R01AG17917, and the Translational Genomics Research Institute. TRAILS has been financially supported by various grants from the Netherlands Organization for Scientific Research NWO (Medical Research Council program grant GB-MW 940-38-011; ZonMW Brainpower grant 100-001-004; ZonMw Risk Behavior and Dependence grants 60-60600-97-118; ZonMw Culture and Health grant 261-98-710; Social Sciences Council medium-sized investment grants GB-MaGW 480-01-006 and GB-MaGW 480-07-001; Social Sciences Council project grants GB-MaGW 452-04-314 and GB-MaGW 452-06-004; NWO large-sized investment grant 175.010.2003.005; NWO Longitudinal Survey and Panel Funding 481-08-013 and 481-11-001), the Dutch Ministry of Justice (WODC), the European Science Foundation (EuroSTRESS project FP-006), Biobanking and Biomolecular Resources Research Infrastructure BBMRI-NL (CP 32), and the participating universities. YFS: Has been financially supported by the Academy of Finland (grants 126925, 121584, 124282, 129378 (Salve), 117787 (Gendi), 41071 (Skidi), and 265869 (Mind)), the Social Insurance Institution of Finland, Kuopio, Tampere and Turku University Hospital Medical Funds (grant 9N035 for Dr. Lehtimäki), Juho Vainio Foundation, Paavo Nurmi Foundation, Finnish Foundation of Cardiovascular Research and Finnish Cultural Foundation, Tampere Tuberculosis Foundation and Emil Aaltonen Foundation (for Dr. Lehtimäki). References Author Manuscript Adams MJ, King JE, Weiss A. The majority of genetic variation in orangutan personality and subjective well-being is nonadditive. Behavior Genetics. 2012; 42:866–866. Bartels M. Genetics of wellbeing and its components satisfaction with life, happiness, and quality of life: A review and meta-analysis of heritability studies. Behavior Genetics. 2015; 45:137–156. [PubMed: 25715755] Bartels M, Boomsma DI. Born to be happy? The etiology of subjective well-being. Behavior Genetics. 2009; 39:605–615. [PubMed: 19728071] Bouchard TJ Jr, Loehlin JC. Genes, evolution, and personality. Behavior Genetics. 2001; 31:243–273. [PubMed: 11699599] Bulik-Sullivan B, Finucane HK, Anttila V, Gusev A, Day FR, Loh PR, Neale BM. An atlas of genetic correlations across human diseases and traits. Nature Genetics. 2015; 47:1236–1241. [PubMed: 26414676] de Moor MH, Costa PT, Terracciano A, Krueger RF, de Geus EJ, Toshiko T, Boomsma DI. Metaanalysis of genome-wide association studies for personality. Molecular Psychiatry. 2012; 17:337– 349. [PubMed: 21173776] Twin Res Hum Genet. Author manuscript; available in PMC 2017 October 01. Weiss et al. Page 12 Author Manuscript Author Manuscript Author Manuscript Author Manuscript de Moor MH, van den Berg SM, Verweij KJ, Krueger RF, Luciano M, Arias Vasquez A, Boomsma DI. Meta-analysis of genome-wide association studies for neuroticism, and the polygenic association with major depressive disorder. JAMA Psychiatry. 2015; 72:642–650. [PubMed: 25993607] DeNeve KM, Cooper H. The happy personality: A meta-analysis of 137 personality traits and subjective well-being. Psychological Bulletin. 1998; 124:197–229. [PubMed: 9747186] Dudbridge F. Power and predictive accuracy of polygenic risk scores. PLoS Genetics. 2013; 9:e1003348. [PubMed: 23555274] Euesden J, Lewis CM, O’Reilly PF. PRSice: Polygenic risk score software. Bioinformatics. 2015; 31:1466–1468. [PubMed: 25550326] Eysenck, HJ.; Eysenck, SBG. Eysenck Personality Questionnaire – Revised (EPQ-R). London: L Hodder & Stoughton; 1991. Hagenaars SP, Harris SE, Davies G, Hill WD, Liewald DC, Ritchie SJ, Deary IJ. Shared genetic aetiology between cognitive functions and physical and mental health in UK Biobank (N = 112 151) and 24 GWAS consortia. Molecular Psychiatry. 2016 Advance online publication. Hahn E, Gottschling J, König CJ, Spinath FM. The heritability of job satisfaction reconsidered: Only unique environmental influences beyond personality. Journal of Business and Psychology. 2016 Advance online publication. Hahn E, Johnson W, Spinath FM. Beyond the heritability of life satisfaction – The roles of personality and twin-specific influences. Journal of Research in Personality. 2013; 47:757–767. Hare E, Contreras J, Raventos H, Flores D, Jerez A, Nicolini H, Escamilla M. Genetic structure of personality factors and bipolar disorder in families segregating bipolar disorder. Journal of Affective Disorders. 2012; 136:1027–1033. [PubMed: 21903278] Keller MC, Coventry WL, Heath AC, Martin NG. Widespread evidence for non-additive genetic variation in Cloninger’s and Eysenck’s personality dimensions using a twin plus sibling design. Behavior Genetics. 2005; 35:707–721. [PubMed: 16273321] Lee SH, Yang J, Goddard ME, Visscher PM, Wray NR. Estimation of pleiotropy between complex diseases using single-nucleotide polymorphism-derived genomic relationships and restricted maximum likelihood. Bioinformatics. 2012; 28:2540–2542. [PubMed: 22843982] Luciano M, Huffman JE, Arias-Vasquez A, Vinkhuyzen AA, Middeldorp CM, Giegling I, Deary IJ. Genome-wide association uncovers shared genetic effects among personality traits and mood states. American Journal of Medical Genetics Part B: Neuropsychiatric Genetics. 2012; 159B:684– 695. Luciano M, MacLeod AK, Payton A, Davies G, Ke XY, Tenesa A, Deary IJ. Effects of gene copy number variants on personality and mood in ageing cohorts. Personality and Individual Differences. 2012; 53:393–397. Middeldorp CM, Cath DC, Van Dyck R, Boomsma DI. The co-morbidity of anxiety and depression in the perspective of genetic epidemiology. A review of twin and family studies. Psychological Medicine. 2005; 35:611–624. [PubMed: 15918338] Middeldorp CM, de Moor MH, McGrath LM, Gordon SD, Blackwood DH, Costa PT, Boomsma DI. The genetic association between personality and major depression or bipolar disorder. A polygenic score analysis using genome-wide association data. Translational Psychiatry. 2011; 1:e50. [PubMed: 22833196] Okbay A, Baselmans BM, De Neve JE, Turley P, Nivard MG, Fontana MA, Cesarini D. Genetic variants associated with subjective well-being, depressive symptoms, and neuroticism identified through genome-wide analyses. Nature Genetics. 2016; 48:624–633. [PubMed: 27089181] Power RA, Pluess M. Heritability estimates of the big five personality traits based on common genetic variants. Translational Psychiatry. 2015; 5:e604. [PubMed: 26171985] Purcell SM, Wray NR, Stone JL, Visscher PM, O’Donovan MC, Sullivan PF, Sklar P. Common polygenic variation contributes to risk of schizophrenia and bipolar disorder. Nature. 2009; 460:748–752. [PubMed: 19571811] Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, Sham PC. PLINK: A tool set for whole-genome association and population-based linkage analyses. American Journal of Human Genetics. 2007; 81:559–575. [PubMed: 17701901] Twin Res Hum Genet. Author manuscript; available in PMC 2017 October 01. Weiss et al. Page 13 Author Manuscript Author Manuscript Russell JA, Carroll JM. On the bipolarity of positive and negative affect. Psychological Bulletin. 1999; 125:3–30. [PubMed: 9990843] van den Berg SM, de Moor MH, McGue M, Pettersson E, Terracciano A, Verweij KJ, Boomsma DI. Harmonization of neuroticism and extraversion phenotypes across inventories and cohorts in the genetics of personality consortium: An application of item response theory. Behavior Genetics. 2014; 44:295–313. [PubMed: 24828478] van den Berg SM, de Moor MM, Verweij KH, Krueger R, Luciano M, Arias Vasquez A, Boomsma DI. Meta-analysis of genome-wide association studies for extraversion: Findings from the genetics of personality consortium. Behavior Genetics. 2015; 46:170–182. [PubMed: 26362575] Weiss A, Bates TC, Luciano M. Happiness is a personal(ity) thing – The genetics of personality and well-being in a representative sample. Psychological Science. 2008; 19:205–210. [PubMed: 18315789] Weiss A, King JE, Enns RM. Subjective well-being is heritable and genetically correlated with dominance in chimpanzees (Pan troglodytes). Journal of Personality and Social Psychology. 2002; 83:1141–1149. [PubMed: 12416918] Wray NR, Lee SH, Mehta D, Vinkhuyzen AA, Dudbridge F, Middeldorp CM. Research review: Polygenic methods and their application to psychiatric traits. Journal of Child Psychology and Psychiatry. 2014; 55:1068–1087. [PubMed: 25132410] Yang JA, Lee SH, Goddard ME, Visscher PM. GCTA: A tool for genome-wide complex trait analysis. American Journal of Human Genetics. 2011; 88:76–82. [PubMed: 21167468] Author Manuscript Author Manuscript Twin Res Hum Genet. Author manuscript; available in PMC 2017 October 01. Weiss et al. Page 14 Author Manuscript Author Manuscript Author Manuscript FIGURE 1. Author Manuscript Neuroticism and extraversion polygenic scores at five SNP inclusion thresholds (x-axis) predicting life satisfaction and positive affect in UK Biobank. Amount of variance explained by the polygenic scores is depicted on the y-axis and the significance value of the polygenic predictor is displayed on the bars. Twin Res Hum Genet. Author manuscript; available in PMC 2017 October 01. Author Manuscript Author Manuscript Author Manuscript Author Manuscript TABLE 1 Twin Res Hum Genet. Author manuscript; available in PMC 2017 October 01. p < .01 p < .05 p < .1 p < .5 p<1 Neuroticism 0.01 (0.007) 0.01 (0.010) 0.016 (0.011) −0.001 (0.007) 0.002 (0.007) p = .16 p = .27 p = .15 p = .70 p = .66 Extraversion 0.014 (0.010) 0.012 (0.007) 0.015 (0.007) 0.012 (0.007) 0.009 (0.007) p = .15 p = .08 p = .031 p = .09 p = .21 Openness −0.012 (0.010) −0.014 (0.01) −0.014 (0.009) −0.012 (0.007) −0.012 (0.007) p = .25 p = .14 p = .12 p = .08 p = .09 Agreeableness −0.008 (0.007) 0 (0.008) 0.002 (0.007) 0.004 (0.008) 0.006 (0.009) p = .28 p = .75 p = .62 p = .56 p = .46 Conscientiousness 0.004 (0.007) 0.01 (0.007) 0.002 (0.007) 0.017 (0.007) 0.015 (0.007) p = .51 p = .16 p = .63 p = .021 p = .042 Weiss et al. Meta-Analysis Results (Regression Coefficient, Standard Error, p Value) for Univariate Analyses of Personality Polygenic Scores (at Five SNP Inclusion Thresholds) Predicting Life Satisfaction (Total N = 19,270) Note: False discovery rate q = 0.002. Page 15 Author Manuscript Author Manuscript Author Manuscript Author Manuscript TABLE 2 p < .01 p < .05 p < .1 p < .5 p<1 −0.006 (0.011)b −0.007 (0.013)b −0.01 (0.014)b −0.019 (0.016)b −0.013 (0.016)b p = .52 p = .51 p = .43 p = .22 p = .37 0.001 (0.005) 0.012 (0.008)a 0.015 (0.009)a 0.02 (0.010)b 0.019 (0.010)b p = .68 p = .10 p = .08 p = .048 p = .047 −0.006 (0.005) −0.001 (0.005) 0.000 (0.005) −0.004 (0.005) −0.003 (0.005) p = .17 p = .66 p = .73 p = .39 p = .45 Agreeableness 0.012 (0.006) 0.02 (0.007)a 0.02 (0.007) 0.020 (0.009)b 0.021 (0.009)a p = .033 p = .006 p = .004 p = .029 p = .019 Conscientiousness 0.004 (0.005) 0.005 (0.005) 0.003 (0.005) 0.002 (0.005) 0.000 (0.005) p = .38 p = .24 p = .50 p = .60 p = .74 Neuroticism Twin Res Hum Genet. Author manuscript; available in PMC 2017 October 01. Extraversion Openness Weiss et al. Meta-Analysis Results (Regression Coefficient, Standard Error, p Value) for Univariate Analyses of Personality Polygenic Scores (at Five SNP Inclusion Thresholds) Predicting Positive Affect (Total N = 46,508) Note: False discovery rate q = 0.002. a Significant heterogeneity p < .05. b significant heterogeneity p < .001. Page 16 Author Manuscript Author Manuscript Author Manuscript Author Manuscript TABLE 3 Twin Res Hum Genet. Author manuscript; available in PMC 2017 October 01. p < .01 p < .05 p < .1 p < .5 p<1 Neuroticism −0.009 −0.017 −0.018 −0.023 −0.026 p = .42 p = .15 p = .13 p = .05 p = .03 Extraversion 0.011 0.022 0.02 0.015 0.015 p = .35 p = .06 p = .09 p = .21 p = .22 Openness 0.011 0.002 0.003 0.01 0.01 p = .34 p = .86 p = .82 p = .39 p = .42 Agreeableness −0.015 0.002 0.008 0.008 0.007 p = .22 p = .87 p = .53 p = .49 p = .55 Conscientiousness −0.011 −0.004 0.012 0.007 0.001 p = .34 p = .73 p = .31 p = .54 p = .69 Weiss et al. Correlation and p Value for Univariate Analyses of Personality Polygenic Scores (at Five SNP Inclusion Thresholds) Predicting General Wellbeing in the MCTFR (N = 6,960) Note: False discovery rate q = 0.002. Page 17