Pleiotropy between neuroticism and physical and mental health: findings from 108 038 men and women in UK Biobank

Abstract

People with higher levels of neuroticism have an increased risk of several types of mental disorder. Higher neuroticism has also been associated, less consistently, with increased risk of various physical health outcomes. We hypothesised that these associations may, in part, be due to shared genetic influences. We tested for pleiotropy between neuroticism and 17 mental and physical diseases or health traits using linkage disequilibrium regression and polygenic profile scoring. Genetic correlations were derived between neuroticism scores in 108 038 people in the UK Biobank and health-related measures from 14 large genome-wide association studies (GWASs). Summary information for the 17 GWASs was used to create polygenic risk scores for the health-related measures in the UK Biobank participants. Associations between the health-related polygenic scores and neuroticism were examined using regression, adjusting for age, sex, genotyping batch, genotyping array, assessment centre and population stratification. Genetic correlations were identified between neuroticism and anorexia nervosa (r _g =0.17), major depressive disorder (r_g =0.66) and schizophrenia (r_g =0.21). Polygenic risk for several health-related measures were associated with neuroticism, in a positive direction in the case of bipolar disorder, borderline personality, major depressive disorder, negative affect, neuroticism (Genetics of Personality Consortium), schizophrenia, coronary artery disease, and smoking (β between 0.009-0.043), and in a negative direction in the case of body mass index (β=-0.0095). A high level of pleiotropy exists between neuroticism and some measures of mental and physical health, particularly major depressive disorder and schizophrenia.

OPEN

Citation: Transl Psychiatry (2016) 6, e791; doi:10.1038/tp.2016.56

http://www.nature.com/tp

Web End =www.nature.com/tp

ORIGINAL ARTICLE

Pleiotropy between neuroticism and physical and mental health: ndings from 108038 men and women in UK Biobank

CR Gale1,2,3,9, SP Hagenaars1,2,4,9, G Davies1,2, WD Hill1,2, DCM Liewald1,2, B Cullen5, BW Penninx6, International Consortium for Blood Pressure GWAS, CHARGE Consortium Aging and Longevity Group, DI Boomsma7, J Pell5, AM McIntosh4, DJ Smith5, IJ Deary1,2,10 and SE Harris1,8,10

People with higher levels of neuroticism have an increased risk of several types of mental disorder. Higher neuroticism has also been associated, less consistently, with increased risk of various physical health outcomes. We hypothesised that these associations may, in part, be due to shared genetic inuences. We tested for pleiotropy between neuroticism and 17 mental and physical diseases or health traits using linkage disequilibrium regression and polygenic prole scoring. Genetic correlations were derived between neuroticism scores in 108 038 people in the UK Biobank and health-related measures from 14 large genome-wide association studies (GWASs). Summary information for the 17 GWASs was used to create polygenic risk scores for the health-related measures in the UK Biobank participants. Associations between the health-related polygenic scores and neuroticism were examined using regression, adjusting for age, sex, genotyping batch, genotyping array, assessment centre and population stratication. Genetic correlations were identied between neuroticism and anorexia nervosa (rg = 0.17), major depressive disorder (rg = 0.66) and schizophrenia (rg = 0.21). Polygenic risk for several health-related measures were associated with neuroticism, in a positive direction in the case of bipolar disorder, borderline personality, major depressive disorder, negative affect, neuroticism (Genetics of Personality Consortium), schizophrenia, coronary artery disease, and smoking ( between 0.0090.043), and in a negative direction in the case of body mass index ( = 0.0095). A high level of pleiotropy exists between neuroticism and

some measures of mental and physical health, particularly major depressive disorder and schizophrenia.

Translational Psychiatry (2016) 6, e791; doi:10.1038/tp.2016.56; published online 26 April 2016

INTRODUCTIONThere is considerable evidence that the personality trait of neuroticism1which describes stable individual differences in the tendency to experience negative emotionshas profound signicance for public health.2 People who are higher in neuroticism have an increased risk of developing Axis I psychopathology, especially the common mental disorders such as mood, anxiety, somatoform and substance use disorders, and also schizophrenia, bipolar disorder and attention-decit hyperactivity disorder (ADHD).36 Higher neuroticism is associated with increased likelihood of diagnosis with Axis II personality disorders7 and with greater comorbidity between internalising disorders (such as major depression, generalised anxiety, panic disorders and phobias) and externalising disorders (such as alcohol and drug dependence, antisocial personality and conduct disorders).8 There is evidence that higher neuroticism is linked with risk of developing Alzheimers disease.9 People who are higher in neuroticism tend to make greater use of mental health services, regardless of whether they have a mental disorder,10 perhaps because they are more likely to perceive a need for care.11 The

estimated economic costs of neuroticism in terms of health-care use and absenteeism are massive.12 Much of these costs relate to reported chronic somatic conditions.12 Higher neuroticism has been linked with increased somatic complaints,13,14 with percep

tion of poorer health,15,16 with future somatic multi-morbidity, as assessed by a count of reported chronic conditions,17 and with

increased likelihood of reporting a range of physical health problems.18

Evidence that neuroticism is predictive of objectively assessed physical health is still relatively sparse and ndings to date are often mixed. For example, whereas some prospective studies have found that higher neuroticism increases mortality from all causes19 or coronary heart disease,20 and is predictive of raised blood pressure21 or body mass index (BMI),22 others have found no such association.2326 In a pooled analysis of data from ve cohorts, there was no consistent association between neuroticism and incidence of type 2 diabetes: higher neuroticism was linked with increased risk in one cohort, but not in others.27

Part of the explanation for associations between neuroticism and these various mental and physical health outcomes may be

1Centre for Cognitive Ageing and Cognitive Epidemiology, University of Edinburgh, Edinburgh, UK; 2Department of Psychology, University of Edinburgh, Edinburgh, UK; 3MRC Lifecourse Epidemiology Unit, University of Southampton, Southampton General Hospital, Southampton, UK; 4Division of Psychiatry, University of Edinburgh, Edinburgh, UK;

5Institute of Health and Wellbeing, University of Glasgow, Glasgow, UK; 6Department of Psychiatry, EMGO Institute for Health and Care Research and Neuroscience Campus Amsterdam, VU University Medical Center/GGZ inGeest, Amsterdam, The Netherlands; 7Department of Biological Psychology, VU University Amsterdam, Amsterdam, The Netherlands and 8Medical Genetics Section, University of Edinburgh Centre for Genomic and Experimental Medicine and MRC Institute of Genetics and Molecular Medicine, Western General Hospital, Edinburgh, UK. Correspondence: Professor C Gale, MRC Lifecourse Epidemiology Unit, University of Southampton, Southampton General Hospital, Southampton SO16 6YD, UK.

E-mail: mailto:[email protected]

Web End [email protected]

9These authors contributed equally to this work.

10These authors contributed equally to this work.

Received 6 October 2015; revised 15 January 2016; accepted 5 March 2016

Pleiotropy between neuroticism and health CR Gale et al

2

due to shared genetic inuences. Twin and adoption studies suggest that genetic inuences account for between a third and a half of individual differences in neuroticism.1 Many physical and mental illnesses and health-related measures also show moderate heritability.28 Twin studies have shown that there is considerable overlap between the genetic factors that inuence variations in neuroticism and those that determine risk of depression and other internalising disorders.29,30 It is now possible to test for such

pleiotropy in associations between neuroticism and health outcomes using data from single nucleotide polymorphism (SNP) genotyping in unrelated individuals, making it possible to carry out much larger studies without the assumptions made by twin-study methods. A recent genome-wide association meta-analysis based on data from over 70 000 individuals found that neuroticism is inuenced by many genetic variants of small effect, that is, a polygenic effect, that also inuence the risk of major depressive disorder.31 Whether there is pleiotropy between neuroticism and other mental disorders or with physical health outcomes is unclear.

Testing for pleiotropy using SNP-based genetic data can be carried out in several ways. Linkage disequilibrium (LD) regression calculates genetic correlations between health measures using the summary results of genome-wide association studies (GWASs).32 It

determines how much of the genetic inuence on two traits/ diseases is common to both. Polygenic risk scoring33 uses summary GWAS data for a given illness or health trait to test whether polygenic liability to that illness/trait is associated with phenotypes for that illness/trait (for example, neuroticism scores) or others measured in an independent sample. It allows the amount of variance in one trait/disease attributed to the polygenic score for a second trait/disease to be calculated. Polygenic risk of neuroticism was recently associated with major depressive disorder.31

In the present study we aimed to discover whether shared genetic aetiology explains part of the associations between neuroticism and various physical and mental health outcomes, all of which have been shown to be phenotypically correlated with neuroticism in at least one study. We used data on over 108 000 UK Biobank participants who completed a questionnaire on neuroticism and provided DNA for genome-wide genotyping. Using summary data from GWAS meta-analyses on 17 health-related measures, we tested for neuroticism-health pleiotropy using two complementary methods. First, we used LD score regression to derive genetic correlations between health-related measures and neuroticism. Second, we calculated the associations between polygenic risk scores for health-related measures and the neuroticism phenotype in UK Biobank participants.

MATERIALS AND METHODS Participants

The participants in this study took part in the baseline survey of UK Biobank.34 (http://www.ukbiobank.ac.uk

Web End =http://www.ukbiobank.ac.uk). UK Biobank was set up as a resource for identifying determinants of disease in middle-aged and older people. Between 2006 and 2010, 502 655 community-dwelling people aged between 37 and 73 years and living in the United Kingdom were recruited to the study. They underwent assessments of cognitive and physical functions, mood and personality. They provided blood, urine and saliva samples for future analysis, completed questionnaires about their social backgrounds and lifestyle and agreed to have their health followed longitudinally. UK Biobank received ethical approval from the Research Ethics Committee (REC reference 11/NW/0382).

For the present study, genome-wide genotyping data were available on 112 151 individuals (58 914 female) aged 4073 years (mean age = 56.9 years, s.d. = 7.9) after the quality control process (see below).

ProceduresNeuroticism. Participants completed the Neuroticism scale of the Eysenck

Personality Questionnaire-Revised (EPQ-R) Short Form.35 This scale has been concurrently validated in older people against two of the most

widely used measures of neuroticism, taken from the International Personality Item Pool (IPIP) and the NEO-Five Factor Inventory (NEO-FFI); it correlated 0.84 with the IPIP-Emotional Stability scale and 0.85 with the NEO-FFI Neuroticism scale.36 A previous study found a high genetic correlation (0.91) between the EPQ-R Short Form Neuroticism scale and psychological distress assessed in a non-psychiatric population using the 30-item General Health Questionnaire.37

Genotyping and quality control. In all, 152 729 UK Biobank blood samples were genotyped using either the UK BiLEVE38 array (N = 49 979) or the UK Biobank axiom array (N = 102 750). A full description of the genotyping process is available in the Supplementary Materials. Quality control was performed by Affymetrix (Santa Clara, CA, USA), the Wellcome Trust Centre for Human Genetics and by the present authors; this included removal of participants based on missingness, relatedness, gender mismatch, non-British ancestry and other criteria, and is described in the Supplementary Materials.

Genome-wide association analyses in the UK Biobank sample. Genome-wide association analyses were performed on the neuroticism measure in order to use the summary results for LD regression. Details of the GWAS procedures are provided in the Supplementary Materials. Results from the GWAS are published elsewhere.39

Curation of summary results from GWAS consortia on health-related variables. In order to conduct LD score regression and polygenic prole score analyses between the UK Biobank neuroticism data and the genetic predisposition to mental and physical health outcomes, we gathered summary data from published meta-analyses on 17 health-related measures: 9 relating to mental health (ADHD, Alzheimers disease, anorexia nervosa, bipolar disorder, borderline personality, major depressive disorder, negative affect, neuroticism (from the Genetics of Personality Consortium (GPC)) and schizophrenia) and 8 relating to physical health (systolic and diastolic blood pressure, BMI, coronary artery disease, longevity, rheumatoid arthritis, smoking status and type 2 diabetes). Details of these health-related variables, the consortias websites, key references and number of subjects included in each consortiums GWAS are given in Supplementary Table 1.

Statistical analysisComputing genetic associations between neuroticism and health-related variables. We use two methods to compute genetic associations between neuroticism and the health-related variables, LD score regression and polygenic prole/risk scoring. Each provides a different metric to infer the existence of pleiotropy between pairs of traits. LD score regression was used to derive genetic correlations to determine the degree to which the polygenic architecture of a trait overlaps with that of another. Genetic correlations were also derived between the different physical and mental health measures, and neuroticism based on the results of the GWAS from the GPC by de Moor et al.31 A genetic correlation between the two neuroticism GWAS (UK Biobank and GPC) was also calculated. The polygenic risk score method was used to test the extent to which the polygenic information from GWASs of health-related variables could predict neuroticism in the UK Biobank participants. Both LD score regression and polygenic risk scores are dependent on the traits analysed being highly polygenic in nature, that is, where a large number of variants of small effect contribute towards phenotypic variation.32,40 LD score

regression was performed between the 14 health-related traits from GWAS consortia, as some measures did not pass thresholds to be included in LD score analysis, while the polygenic prole score analyses were performed on the complete set of 17 health-related traits from GWAS consortia.

LD score regression. In order to quantify the extent of pleiotropy between neuroticism, measured in UK Biobank, and the collated health traits, we used LD score regression.32,41 This is a class of techniques that exploits the

correlational structure of the SNPs found across the genome (we provide more details of LD score regression in the Supplementary Materials). Here we use LD score regression to derive genetic correlations between neuroticism and health-related measures using 14 large GWAS consortia data sets that enable pleiotropy of their health-related traits to be quantied with the neuroticism trait in UK Biobank. We followed the data-processing pipeline devised by Bulik-Sullivan et al.32,41 described in more detail in the Supplementary Materials. In order to ensure that the genetic correlation for the Alzheimers disease phenotype was not driven by a

Translational Psychiatry (2016), 1 7

Pleiotropy between neuroticism and health CR Gale et al

3

Table 1. Genetic correlations between neuroticism documented in UK Biobank and mental and physical health-related traits curated from GWAS consortia

Trait Category Traits from GWAS consortia Neuroticism (UKB) (n = 108 038) Neuroticism (GPC) (n = 63 661)

rg s.e. P rg s.e. P

Mental health ADHD 0.060 0.082 0.4681 0.090 0.191 0.6371Alzheimers disease 0.118 0.080 0.1378 0.089 0.177 0.6147 Alzheimers disease (500 kb) 0.091 0.063 0.1514 0.051 0.145 0.7262 Anorexia nervosa 0.174 0.04 2.36 105 0.319 0.082 9.56 105 Bipolar disorder 0.083 0.052 0.1096 0.162 0.101 0.1091 Major depressive disorder 0.659 0.087 2.75 1014 1.073 0.199 7.15 108 Schizophrenia 0.212 0.050 2.39 105 0.198 0.068 0.0033

Physical health Blood pressure: diastolic 0.011 0.078 0.887 0.108 0.091 0.2363 Blood pressure: systolic 0.051 0.060 0.3959 0.008 0.083 0.925 BMI 0.036 0.028 0.1921 0.046 0.055 0.4077 Coronary artery disease 0.104 0.053 0.0489 0.221 0.104 0.0331 Longevity 0.083 0.081 0.3011 0.002 0.163 0.9888 Rheumatoid arthritis 0.0816 0.05 0.103 0.093 0.123 0.4463 Smoking status 0.113 0.05 0.0205 0.124 0.090 0.1657 Type 2 diabetes 0.145 0.07 0.0389 0.044 0.107 0.6837

Abbreviations: ADHD, attention-decit hyperactivity disorder; BMI, body mass index; GPC, Genetics of Personality Consortium; GWAS, genome-wide association study; UKB, UK Biobank. Statistically signicant P-values (after false discovery rate correctionPo0.0033) are shown in bold.

single locus or biased the t of the regression model, a 500-kb region centred on the APOE locus was removed and this phenotype was re-run. This additional model is referred to in Table 1 as Alzheimers disease (500 kb).

Polygenic proling. The UK Biobank genotyping data required recoding from numeric (1, 2) allele coding to standard ACGT format before being used in polygenic prole scoring analyses. This was achieved using a bespoke programme developed by one of the present authors (DCML), details of which are provided in the Supplementary Materials.

Polygenic proles were created for 17 health-related phenotypes (Table 2 and Supplementary Table 2) in all genotyped participants using PRSice.42 This software calculates the sum of alleles associated with the phenotype of interest across many genetic loci, weighted by their effect sizes estimated from a GWAS of that phenotype in an independent sample. Before creating the scores, SNPs with a minor allele frequency o0.01 were removed, and clumping was used to obtain SNPs in linkage equilibrium with an r2o0.25 within a 200-bp window. Multiple scores were then created for each phenotype containing SNPs selected according to the signicance of their association with the phenotype. The GWAS summary data for the 12 health-related phenotypes were used to create ve polygenic proles for each in the UK Biobank participants, at thresholds of Po0.01, Po0.05, Po0.1, Po0.5 and all SNPs. The most predictive threshold will be presented in the main tables of this paper. The full results, including all ve thresholds, can be found in Supplementary Table 2.

Associations between the polygenic proles and neuroticism were examined in linear regression models, adjusting for age at measurement, sex, genotyping batch and array, assessment centre and the rst 10 genetic principal components to adjust for population stratication. We corrected for multiple testing across all polygenic prole scores at all signicance thresholds for associations with neuroticism using the false discovery rate method.43 As neuroticism on average is higher in females and likely declines with age,2 age- (under and over 60 years old) and gender-stratied models were examined.

RESULTSWithin UK Biobank, 108 038 individuals with genotype data completed the Neuroticism scale of the EPQ-R Short Form. Their mean (s.d.) score for neuroticism was 4.02 (3.17).

The genetic correlation between neuroticism measured in the UK Biobank and neuroticism measured in the GPC is 1.0 (s.e. = 0.11).

Table 1 shows the genetic correlations obtained using LD score regression between neuroticism in UK Biobank and neuroticism

from the GPC, and the published GWAS results on the health-related traits. Neuroticism in UK Biobank showed signicant positive genetic correlations with three traits, all related to mental health, namely major depressive disorder (rg = 0.66), schizophrenia (rg = 0.21) and anorexia nervosa (rg = 0.17). The ndings replicated when using Neuroticism from the GPC. There were no signicant genetic correlations between neuroticism and either the other mental health-related traits (ADHD, Alzheimers disease and bipolar disorder) or any of the physical health-related traits (systolic and diastolic blood pressure, BMI, coronary artery disease, type 2 diabetes, smoking status, rheumatoid arthritis or longevity). These results are shown graphically in Figure 1.

Table 2 shows the results of the polygenic risk scoring, using the most predictive threshold of the ve that were created. Higher polygenic risk for six mental health-related traits, bipolar disorder, borderline personality, major depressive disorder, negative affect, neuroticism (calculated using summary data from the GPC) and schizophrenia, was signicantly associated with higher levels of neuroticism in UK Biobank (standardised between 0.017 and0.043). There were no signicant associations between neuroticism and polygenic risk for the other mental health-related traits examined, namely ADHD, Alzheimers disease and anorexia nervosa. Polygenic risk scores for three physical health-related traits were signicantly associated with neuroticism: higher polygenic risk for BMI was associated with lower levels of neuroticism ( = 0.0095, P = 0.0015), higher polygenic risk for coronary artery disease was associated with higher levels of neuroticism ( = 0.011, P = 0.0003) and higher polygenic prole scores for smoking were associated with higher levels of neuroticism ( = 0.17, P = 2.48 107). No signicant associations were found between polygenic risk for the other physical health-related traits (systolic and diastolic blood pressure, type 2 diabetes, rheumatoid arthritis or longevity) and neuroticism. To test whether the false discovery rate signicant association between polygenic risk for coronary artery disease and neuroticism was confounded by individuals diagnosed with cardiovascular disease, 2717 individuals who had had a heart attack and 2468 individuals with angina were removed from the regression analysis. Our estimate of the association between polygenic risk of coronary artery disease and neuroticism was unchanged by this exclusion. The same applied for major depressive disorder (excluding 7494 individuals with a probable diagnosis of major

Translational Psychiatry (2016), 1 7

Pleiotropy between neuroticism and health CR Gale et al

4

Table 2. Associations between polygenic risk scores for mental and physical health-related traits created from GWAS consortia summary data and neuroticism in UK Biobank participants, adjusted for age, sex, assessment centre, genotyping batch and array, and 10 principal components for population stratication

Trait category Traits from GWAS consortia Neuroticism (n = 108 038)

Threshold R2 P

Mental health ADHD 0.05 0.0045 2.06 105 0.1292Alzheimers disease 0.05 0.0065 4.16 105 0.0312 Anorexia nervosa 0.1 0.0054 2.87 105 0.0735 Bipolar disorder 0.5 0.0171 0.0003 1.71 108

Borderline personality 1 0.0150 0.0002 5.43 107

Major depressive disorder 1 0.0357 0.0012 1.23 1031

Negative affect (anxiety) 1 0.00949 8.99 105 1.55 103 Neuroticism (GPC) 0.5 0.0433 0.0020 1.77 1047

Schizophrenia 0.1 0.0359 0.0012 7.88 1032

Physical health Blood pressure: diastolic 0.1 0.0047 2.20 105 0.1167Blood pressure: systolic 0.05 0.0016 2.51 106 0.5969BMI 0.01 0.0095 9.00 105 0.0015 Coronary artery disease 0.1 0.0109 0.00012 0.0003 Longevity 0.05 0.0051 2.55 105 0.0917 Rheumatoid arthritis 0.5 0.0007 5.09 105 0.8266 Smoking status 0.5 0.0167 0.0002 2.48 107

Type 2 diabetes 0.01 0.0031 9.67 106 0.2989

Abbreviations: ADHD, attention-decit hyperactivity disorder; BMI, body mass index; GPC, Genetics of Personality Consortium; GWAS, genome-wide association study. The associations between the polygenic risk scores and neuroticism with the largest effect size (threshold) are presented. Statistically signicant P-values (after false discovery rate correctionPo0.0065) are shown in bold.

Figure 1. Barplot of genetic correlations (s.e.) calculated using linkage disequilibrium score regression between neuroticism in UK Biobank (UKB) and the Genetics of Personality Consortium (GPC), and mental and physical health measures from genome-wide association study consortia. *Po0.0033. ADHD, attention-decit hyperactivity disorder; BMI, body mass index.

depression or bipolar disorder44), where the estimate of the association showed little change. The complete polygenic risk score results, including all ve thresholds, are shown in Supplementary Table 2.

Age- and gender-stratied analyses indicated there were no substantial differences by age or gender for 13 out of 17 traits. Four traits did show a potential age or gender effect. Polygenic risk for anorexia nervosa predicts neuroticism in males, but not in females. Polygenic risk for diastolic blood pressure predicts neuroticism only in females over 60 years of age. The strongest

associations between polygenic risk for negative affect and smoking status, and neuroticism, are in females under 60 years of age. The full results for these stratied analyses can be found in Supplementary Table 3.

DISCUSSIONIn this study of 108 038 men and women from UK Biobank who had been genotyped and assessed for neuroticism, we exploited the summary results of 17 large international GWAS consortia to

Translational Psychiatry (2016), 1 7

Pleiotropy between neuroticism and health CR Gale et al

5

examine whether there is pleiotropy between neuroticism and a range of physical and mental health outcomes using two methods, LD score regression and polygenic prole scoring. Summary results from two separate GWASs of neuroticism showed a genetic correlation of 1.0. Using summary data from both of these studies of neuroticism, the genetic correlations that were calculated with physical and mental health were very similar. As regards the six mental health outcomes that were investigated using both LD score regression and polygenic prole scoring, we found consistent evidence of pleiotropy between neuroticism and both major depression and schizophrenia with these methods, showing that, to a signicant degree, the same genetic variants are responsible for the heritability of each pair of phenotypes and that genetic variants associated with major depression or schizophrenia in GWAS consortia are signicantly predictive of variation in neuroticism in UK Biobank. There was some evidence for pleiotropy between neuroticism, and bipolar disorder, borderline personality, anorexia nervosa and negative affect on the basis of results from polygenic prole scoring or LD score regression, respectively, but the extent of these associations varied according to the method used. In all cases where there was a signicant nding using one method but not the other, the direction of effect was the same using both methods. We found no evidence of pleiotropy between neuroticism and the other mental health outcomes examinedADHD and Alzheimers disease. Of the eight physical health outcomes studied, none showed evidence of pleiotropy with neuroticism on the basis of the genetic correlations obtained from LD score regression, but there was some indication of genetic overlap between neuroticism and coronary artery disease, smoking status and BMI. Higher polygenic risk for coronary artery disease and smoking status was signicantly associated with higher levels of neuroticism, and polygenic risk for higher BMI was associated with lower levels of neuroticism. No associations were found between polygenic risk for the other physical health-related traits (systolic and diastolic blood pressure, rheumatoid arthritis, type 2 diabetes or longevity) and neuroticism.

Previous investigations into pleiotropy between neuroticism and mental health outcomes using polygenic risk proling have found evidence of substantial shared genetic aetiology between neuroticism and major depression.31,45,46 Our observations in the

present, much larger, sample conrm those ndings and for, we believe, the rst time quantify the extent to which the same genetic variants are responsible for the heritability in these two phenotypes, using an additional metric, LD score regression.32 The

relatively high genetic correlation between neuroticism and major depression (rg = 0.66), identied in our study, is similar to the genetic correlation identied in a previous twin study (rg = 0.43).30 Our results also provide the rst evidence, to our knowledge, that the phenotypic correlations found between neuroticism and schizophrenia47,48 are due at least in part to genetic overlap.

The extent of pleiotropy between anorexia nervosa or bipolar disorder and neuroticism has been unclear, though there is some, though limited, evidence to link both disorders phenotypically with neuroticism.5,49 A previous study found no signicant association between polygenic risk scores for neuroticism and bipolar disorder, but the sample size was small.46 In this much

larger sample, higher polygenic risk scores for bipolar disorder were signicantly predictive of higher neuroticism, but results of LD score regression showed little indication of genetic correlation between neuroticism and this disorder. We found a small but highly signicant genetic correlation between neuroticism and anorexia (rg = 0.17), but there was no association between polygenic risk scores for this condition and neuroticism. Although there is now considerable evidence that neuroticism is a risk factor for the development of Alzheimers disease,9 there was no indication in our analyses that shared genes account for this link.

So far as we are aware, there have been no previous investigations of the extent of pleiotropy between neuroticism and physical health outcomes. This may be because, whereas there is considerable evidence for phenotypic associations between higher neuroticism and poorer self-rated health or greater somatic complaints,1318 fewer studies have examined neuroticism as a predictor of objectively measured physical health outcomes, and ndings on such outcomes as coronary heart disease, blood pressure, BMI and all-cause mortality have been inconsistent.1926 Of the eight objectively measured physical health outcomes included in the current study, three showed evidence of a degree of genetic pleiotropy with neuroticism: coronary artery disease, smoking status and BMI. Neither demonstrated any measurable genetic correlation with neuroticism, but higher polygenic risk score for coronary artery disease was associated with higher neuroticism. This is consistent with the nding in pooled data from three cohorts that higher neuroticism was associated with increased mortality from coronary heart disease.20 Higher polygenic risk for BMI was associated with lower neuroticism. The direction of this association was unexpected. Findings on the phenotypic relationships between neuroticism and BMI have produced inconsistent results: one study showed higher neuroticism was associated with higher BMI,22 but another found no association.25

The chief strength of our study is the large sample size that permits powerful, robust tests of genetic association. Second, all the UK Biobank genetic data were processed at the same location and on the same platform. Finally, use of summary data from 17 large international GWAS consortia studies allowed us to perform a comprehensive examination of the degree of pleiotropy between neuroticism and a range of physical and mental health-related phenotypes, and to produce many of the rst estimates of the genetic correlation between neuroticism and these phenotypes.

Our study also has some limitations. First, the GWAS studies we curated to carry out LD score regression and extract polygenic risk scores often involved meta-analyses of results from data sets with considerable heterogeneity in sample size, genome-wide imputation quality and measurement of phenotypes. With larger and more consistent independent data sets, it should be possible to use the polygenic risk scores to predict more variance in neuroticism. Second, we restricted the genotyped samples to individuals of white British ancestry in order to minimise any inuence of population structure. Our results therefore need to be replicated in large samples with different genetic backgrounds as we did not have the power to model data from UK Biobank individuals of other ancestries. In the stratied analysis 13 out of 17 traits did not indicate a gender or age effect. These stratied results need to be interpreted with caution; they were mostly null, and no multiple testing correction was applied due to the high correlations between the different models. The four traits that did indicate potential gender or age effects need to be replicated in an independent sample, before any conclusions can be drawn from them.

We note that genetic correlations reect the amount of the genetic inuence on two traits that is common to both. This is independent of the heritability of either trait. Therefore, it is entirely possible for a trait to have a small proportion of variance accounted for by genetic variants, but to have a high genetic correlation with another trait. The amount of variance explained by each polygenic prole score is small, as would be expected by the fact that not all SNPs were genotyped, and those that were do not necessarily accurately tag the causal genetic variants. The polygenic risk score most predictive threshold varied between health traits, suggesting that the amount of shared genetic aetiology between neuroticism and each of the health traits differs. Although testing multiple thresholds may be deemed to increase the multiple testing problem, it should be noted that the

Translational Psychiatry (2016), 1 7

Pleiotropy between neuroticism and health CR Gale et al

6

SNPs included in each threshold are not independent. Also, at least for major depressive disorder and schizophrenia, pleiotropy was quantied using LD score regression, which involved only a single test per pair of phenotypes. The estimate of neuroticism variance explained by each polygenic prole should be considered as the minimum estimate of the variance explained. Owing to pruning SNPs in LD, the polygenic risk score method makes the assumption of a single causal variant being tagged in each LD block considered. If this assumption is not true for the phenotypes considered, the proportion of variance explained will be underestimated here.

In this large sample from UK Biobank, we aimed to discover whether shared genetic aetiology explained part of the associations between neuroticism and various physical and mental health outcomes. Our ndings suggest that associations between neuroticism and several mental health outcomes including major depression, schizophrenia, bipolar disorder and anorexia nervosa are in part due to shared genetic inuences. We found that polygenic risk scores for coronary artery disease, smoking and BMI were predictive of neuroticism scores. This large-scale mapping of the extent of pleiotropy between neuroticism and physical and mental health outcomes adds to our understanding of the cause of links between this important personality trait and later health.

CONFLICT OF INTEREST

The authors declare no conict of interest.

ACKNOWLEDGMENTS

This research has been conducted using the UK Biobank Resource. The work was undertaken in The University of Edinburgh Centre for Cognitive Ageing and Cognitive Epidemiology, part of the cross council Lifelong Health and Wellbeing Initiative (MR/ K026992/1). We thank Riccardo Marioni for conducting some analyses for this paper. Funding from the BBSRC and the Medical Research Council (MRC) is gratefully acknowledged.

REFERENCES

1 Matthews G, Deary IJ, Whiteman MC. Personality Traits. 3rd edn. Cambridge

University Press: Cambridge, UK, 2009.2 Lahey BB. Public health signicance of neuroticism. Am Psychol 2009; 64: 241256.3 Kotov R, Gamez W, Schmidt F, Watson D. Linking big personality traits to anxiety, depressive, and substance use disorders: a meta-analysis. Psychol Bull 2010; 136: 768821.4 Malouff JM, Thorsteinsson EB, Schutte NS. The relationship between the ve-factor model of personality and symptoms of clinical disorders: a meta-analysis. J Psychopathol Behav Assess 2005; 27: 101114.5 Lonnqvist JE, Verkasalo M, Haukka J, Nyman K, Tiihonen J, Laaksonen I et al. Premorbid personality factors in schizophrenia and bipolar disorder: results from a large cohort study of male conscripts. J Abnorm Psychol 2009; 118: 418423.6 Michielsen M, Comijs HC, Semeijn EJ, Beekman AT, Deeg DJ, Kooij JJ. Attention decit hyperactivity disorder and personality characteristics in older adults in the general Dutch population. Am J Geriatr Psychiatry 2014; 22: 16231632.7 Saulsman LM, Page AC. The ve-factor model and personality disorder empirical literature: a meta-analytic review. Clin Psychol Rev 2004; 23: 10551085.8 Khan AA, Jacobson KC, Gardner CO, Prescott CA, Kendler KS. Personality and comorbidity of common psychiatric disorders. Br J Psychiatry 2005; 186: 190196.9 Terracciano A, Sutin AR, An Y, O'Brien RJ, Ferrucci L, Zonderman AB et al. Personality and risk of Alzheimers disease: new data and meta-analysis. Alzheimers Dement 2014; 10: 179186.10 ten Have M, Oldehinkel A, Vollebergh W, Ormel J. Does neuroticism explain variations in care service use for mental health problems in the general population? Results from the Netherlands Mental Health Survey and Incidence Study (NEMESIS). Soc Psychiatry Psychiatr Epidemiol 2005; 40: 425431.11 Seekles WM, Cuijpers P, van de Ven P, Penninx BW, Verhaak PF, Beekman AT et al. Personality and perceived need for mental health care among primary care patients. J Affect Disord 2012; 136: 666674.12 Cuijpers P, Smit F, Penninx BW, de Graaf R, ten Have M, Beekman AT. Economic costs of neuroticism: a population-based study. Arch Gen Psychiatry 2010; 67: 10861093.

13 Costa PT, Mccrae RR. Neuroticism, somatic complaints, and diseaseis the bark worse than the bite. J Pers 1987; 55: 299316.

14 Jorm AF, Christensen H, Henderson S, Korten AE, Mackinnon AJ, Scott R. Neuroticism and self-reported health in an elderly community sample. Pers Indiv Differ 1993; 15: 515521.

15 Goodwin R, Engstrom G. Personality and the perception of health in the general population. Psychol Med 2002; 32: 325332.

16 Watson D, Pennebaker JW. Health complaints, stress, and distress: exploring the central role of negative affectivity. Psychol Rev 1989; 96: 234254.

17 Neeleman J, Bijl R, Ormel J. Neuroticism, a central link between somatic and psychiatric morbidity: path analysis of prospective data. Psychol Med 2004; 34: 521531.

18 Goodwin RD, Cox BJ, Clara I. Neuroticism and physical disorders among adults in the community: Results from the national comorbidity survey. J Behav Med 2006; 29: 229238.

19 Shipley BA, Weiss A, Der G, Taylor MD, Deary IJ. Neuroticism, extraversion, and mortality in the UK Health and Lifestyle Survey: a 21-year prospective cohort study. Psychosom Med 2007; 69: 923931.

20 Jokela M, Pulkki-Raback L, Elovainio M, Kivimaki M. Personality traits as risk factors for stroke and coronary heart disease mortality: pooled analysis of three cohort studies. J Behav Med 2014; 37: 881889.

21 Turiano NA, Pitzer L, Armour C, Karlamangla A, Ryff CD, Mroczek DK. Personality trait level and change as predictors of health outcomes: ndings from a national study of Americans (MIDUS). J Gerontol B Psychol Sci Soc Sci 2012; 67: 412.

22 Magee CA, Heaven PCL. Big-Five personality factors, obesity and 2-year weight gain in Australian adults. J Res Pers 2011; 45: 332335.

23 Nakaya N, Tsubono Y, Hosokawa T, Hozawa A, Kuriyama S, Fukudo S et al. Personality and mortality from ischemic heart disease and stroke. Clin Exp Hypertens 2005; 27: 297305.

24 Jokela M, Batty GD, Nyberg ST, Virtanen M, Nabi H, Singh-Manoux A et al. Personality and all-cause mortality: individual-participant meta-analysis of 3,947 deaths in 76,150 adults. Am J Epidemiol 2013; 178: 667675.

25 Mottus R, McNeill G, Jia X, Craig LC, Starr JM, Deary IJ. The associations between personality, diet and body mass index in older people. Health Psychol 2013; 32: 353360.

26 Almada SJ, Zonderman AB, Shekelle RB, Dyer AR, Daviglus ML, Costa PT Jr. et al. Neuroticism and cynicism and risk of death in middle-aged men: the Western Electric Study. Psychosom Med 1991; 53: 165175.

27 Jokela M, Elovainio M, Nyberg ST, Tabak AG, Hintsa T, Batty GD et al. Personality and risk of diabetes in adults: pooled analysis of 5 cohort studies. Health Psychol 2014; 33: 16181618.

28 Polderman TJ, Benyamin B, de Leeuw CA, Sullivan PF, van Bochoven A, Visscher PM et al. Meta-analysis of the heritability of human traits based on fty years of twin studies. Nat Genet 2015; 47: 702709.

29 Hettema JM, Neale MC, Myers JM, Prescott CA, Kendler KS. A population-based twin study of the relationship between neuroticism and internalizing disorders. Am J Psychiatry 2006; 163: 857864.

30 Kendler KS, Myers J. The genetic and environmental relationship between major depression and the ve-factor model of personality. Psychol Med 2010; 40: 801806.

31 Genetics of Personality Consortium, de Moor MH, van den Berg SM, Verweij KJ, Krueger RF, Luciano M et al. Meta-analysis of genome-wide association studies for neuroticism, and the polygenic association with major depressive disorder. JAMA Psychiatry 2015; 72: 642650.

32 Bulik-Sullivan BK, Loh PR, Finucane HK, Ripke S, Yang J, Schizophrenia Working Group of the Psychiatric Genomics C, et al. LD score regression distinguishes confounding from polygenicity in genome-wide association studies. Nat Genet 2015; 47: 291295.

33 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: 559575.

34 Sudlow C, Gallacher J, Allen N, Beral V, Burton P, Danesh J et al. UK biobank: an open access resource for identifying the causes of a wide range of complex diseases of middle and old age. PLoS Med 2015; 12: e1001779.

35 Deary IJ, Bedford A. Some origins and evolution of the EPQ-R (short form) neuroticism and extraversion items. Pers Indiv Differ 2011; 50: 12131217.

36 Gow AJ, Whiteman MC, Pattie A, Deary IJ. Goldbergs IPIP Big-Five factor markers: Internal consistency and concurrent validation in Scotland. Pers Indiv Differ 2005; 39: 317329.

37 Ivkovic V, Vitart V, Rudan I, Janicijevic B, Smolej-Narancic N, Skaric-Juric T et al. The Eysenck personality factors: psychometric reliability, heritability and phenotypic and genetic correlations with psychological distress in an isolated Croatian population. Pers Indiv Differ 2007; 42: 123133.

38 Wain LV, Shrine N, Miller S, Jackson VE, Ntalla I, Soler Artigas M. Novel insights into the genetics of smoking behaviour, lung function, and chronic obstructive pulmonary disease (UK BiLEVE): a genetic association study in UK Biobank. Lancet Respir Med 2015; 3: 769781.

Translational Psychiatry (2016), 1 7

Pleiotropy between neuroticism and health CR Gale et al

7

39 Smith DJ, Escott-Price V, Davies G, Bailey MES, Colodro Conde L, Ward J et al. Genome-wise analysis of over 106,000 individuals identies 9 neuroticism-associated loci. Mol Psychiatry 2016 (in press).

40 International Schizophrenia C, Purcell SM, Wray NR, Stone JL, Visscher PM, ODonovan MC et al. Common polygenic variation contributes to risk of schizophrenia and bipolar disorder. Nature 2009; 460: 748752.

41 Bulik-Sullivan B, Finucane HK, Anttila V, Gusev A, Day FR, Perry JR et al. An atlas of genetic correlations across human diseases and traits. bioRxiv 2015; 47: 12361241.

42 Eusden J, Lewis CM, OReilly PF. PRSice: Polygenic Risk Score software. Bioinformatics 2015; 31: 14661468.

43 Benjamini Y, Hochberg Y. Controlling the false discovery ratea practical and powerful approach to multiple testing. J Roy Stat Soc B Met 1995; 57: 289300.

44 Smith DJ, Nicholl BI, Cullen B, Martin D, Ul-Haq Z, Evans J et al. Prevalence and characteristics of probable major depression and bipolar disorder within UK Biobank: cross-sectional study of 172,751 participants. PLoS One 2013; 8: e75362.

45 Luciano M, Huffman JE, Arias-Vasquez A, Vinkhuyzen AA, Middeldorp CM, Giegling I et al. Genome-wide association uncovers shared genetic effects among personality traits and mood states. Am J Med Genet B Neuropsychiatr Genet 2012; 159B: 684695.

46 Middeldorp CM, de Moor MH, McGrath LM, Gordon SD, Blackwood DH, Costa PT et al. The genetic association between personality and major depression or bipolar disorder. A polygenic score analysis using genome-wide association data.Transl Psychiatry 2011; 1: e50.

47 Horan WP, Blanchard JJ, Clark LA, Green MF. Affective traits in schizophrenia and schizotypy. Schizophr Bull 2008; 34: 856874.

48 Van Os J, Jones PB. Neuroticism as a risk factor for schizophrenia. Psychol Med 2001; 31: 11291134.

49 Bulik CM, Sullivan PF, Tozzi F, Furberg H, Lichtenstein P, Pedersen NL. Prevalence, heritability, and prospective risk factors for anorexia nervosa. Arch Gen Psychiatry 2006; 63: 305312.

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articles Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

Web End =http://creativecommons.org/licenses/ http://creativecommons.org/licenses/by/4.0/

Web End =by/4.0/

Supplementary Information accompanies the paper on the Translational Psychiatry website (http://www.nature.com/tp)

Translational Psychiatry (2016), 1 7