Hostility in adolescents and adults: a genome-wide association study of the Young Finns

Abstract

Hostility is a multidimensional personality trait with changing expression over the life course. We performed a genome-wide association study (GWAS) of the components of hostility in a population-based sample of Finnish men and women for whom a total of 2.5 million single-nucleotide polymorphisms (SNPs) were available through direct or in silico genotyping. Hostility dimensions (anger, cynicism and paranoia) were assessed at four time points over a 15-year interval (age range 15-30 years at phase 1 and 30-45 years at phase 4) in 982-1780 participants depending on the hostility measure. Few promising areas from chromosome 14 at 99 cM (top SNPs rs3783337, rs7158754, rs3783332, rs2181102, rs7159195, rs11160570, rs941898, P values <3.9 × 10 ^-8 with nearest gene Enah/Vasp-like (EVL)) were found suggestively to be related to paranoia and from chromosome 7 at 86 cM (top SNPs rs802047, rs802028, rs802030, rs802026, rs802036, rs802025, rs802024, rs802032, rs802049, rs802051, P values <6.9 × 10 ^-7 with nearest gene CROT (carnitine O-octanoyltransferase)) to cynicism, respectively. Some shared suggestive genetic influence for both paranoia and cynicism was also found from chromosome 17 at 2.8 cM (SNPs rs12936442, rs894664, rs6502671, rs7216028) and chromosome 22 at 43 cM (SNPs rs7510759, rs7510924, rs7290560), with nearest genes RAP1 GTPase activating protein 2 (RAP1GAP2) and KIAA1644, respectively. These suggestive associations did not replicate across all measurement times, which warrants further study on these SNPs in other populations.

Citation: Transl Psychiatry (2011) 1, e11, doi:10.1038/tp.2011.13

& 2011 Macmillan Publishers Limited All rights reserved 2158-3188/11 http://www.nature.com/tp

Web End =www.nature.com/tp

Hostility in adolescents and adults: a genome-wide association study of the Young Finns

P Merjonen1, L Keltikangas-Jarvinen1, M Jokela1, I Seppala2,3, L-P Lyytikanen2,3, L Pulkki-Rback1,4, M Kivimaki1, M Elovainio1,5, J Kettunen6,7, S Ripatti5,7, M Kahnen8, J Viikari9, A Palotie6,7, L Peltonen6,7,{, OT Raitakari10 and T Lehtimaki2,3

Hostility is a multidimensional personality trait with changing expression over the life course. We performed a genome-wide association study (GWAS) of the components of hostility in a population-based sample of Finnish men and women for whom a total of 2.5 million single-nucleotide polymorphisms (SNPs) were available through direct or in silico genotyping. Hostility dimensions (anger, cynicism and paranoia) were assessed at four time points over a 15-year interval (age range 1530 years at phase 1 and 3045 years at phase 4) in 9821780 participants depending on the hostility measure. Few promising areas from chromosome 14 at 99 cM (top SNPs rs3783337, rs7158754, rs3783332, rs2181102, rs7159195, rs11160570, rs941898, P values o3.9 10 8 with nearest gene Enah/Vasp-like (EVL)) were found suggestively to be related to paranoia and from chromosome 7

at 86 cM (top SNPs rs802047, rs802028, rs802030, rs802026, rs802036, rs802025, rs802024, rs802032, rs802049, rs802051, P values o6.9 10 7 with nearest gene CROT (carnitine O-octanoyltransferase)) to cynicism, respectively. Some shared

suggestive genetic inuence for both paranoia and cynicism was also found from chromosome 17 at 2.8 cM (SNPs rs12936442, rs894664, rs6502671, rs7216028) and chromosome 22 at 43 cM (SNPs rs7510759, rs7510924, rs7290560), with nearest genes RAP1 GTPase activating protein 2 (RAP1GAP2) and KIAA1644, respectively. These suggestive associations did not replicate across all measurement times, which warrants further study on these SNPs in other populations.

Translational Psychiatry (2011) 1, e11; doi:http://dx.doi.org/10.1038/tp.2011.13

Web End =10.1038/tp.2011.13 ; published online 21 June 2011

Introduction

Hostility is a personality trait characterizing how trustworthy individuals perceive other people and how they handle these feelings toward others. The cognitive component of hostility characterizes cynical and distrustful attitudes, which is the primary reference of the term hostility,1 whereas the affective component reects feelings of irritability and anger.The behavioral component covers expression of these attitudes and feelings as either expressing them out, that is, aggression, or as suppressing or repressing them. Hostility traits have been found to be related to various social and health problems, such as criminality and violence,2,3 isolation and relationship aggression,4 depression,5 cardiovascular diseases6 and all-cause mortality risk,7 although the ndings are not entirely consistent.8 Identifying the origins of hostility may help to understand the developmental paths related to hostility and to develop effective preventions to reduce problems related to hostile behaviors.

Both genetic and environmental factors are involved in the development of hostility,9 with heritability estimated to be B3050%.1012 However, the molecular nature of the genetic

background and the specic regions of the genome that underlie hostility remain mainly unknown. To our knowledge, only one genome scan study of hostility has been published to date.13 That study covered 387 autosomal short-tandem-repeat polymorphisms and did not nd signicant linkage with hostility.13 In the present study, we report a large-scale genome-wide association study (GWAS) of hostility where over 2.5 million single-nucleotide polymorphisms (SNPs) were analyzed, thereby mapping the most potentially signicant areas of the genome regarding hostility for further inspection and providing preliminary evidence of the genetic basis of hostility.

As cognitive, affective and behavioral components of hostility may vary in their etiology and have different genetic backgrounds,14 we used three different scales of hostility, each of which was measured four times over a 15-year time span extending from adolescence and young adulthood (age 1530 years) into adulthood (age 3045 years) in a Caucasian Finnish population. It has been argued that personality is still transient and amenable to environmental effects in young adulthood, but between 30 and 50 years of age, it begins to stabilize and genetic effects become more prominent.11 Thus,

1IBS, Unit of Personality, Work and Health Psychology, University of Helsinki, Helsinki, Finland; 2Department of Clinical Chemistry, Centre for Laboratory Medicine, Tampere University Hospital, Tampere, Finland; 3Medical School, University of Tampere, Tampere, Finland; 4Finnish Institute of Occupational Health, Helsinki, Finland;

5THL, Helsinki, Finland; 6Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge, UK; 7FIMM, Institute for Molecular Medicine Finland, University of Helsinki, Helsinki, Finland; 8Department of Clinical Physiology, Tampere University Hospital, Tampere, Finland; 9Department of Medicine, University of Turku and Turku University Hospital, Turku, Finland and 10Department of Clinical Physiology, Turku University Hospital and Research Centre of Applied and Preventive Cardiovascular Medicine, University of Turku, Turku, FinlandCorrespondence: Professor L Keltikangas-Jarvinen, IBS, Unit of Personality, Work and Health Psychology, University of Helsinki, P.O. Box 9, FIN-00014 University of Helsinki, Helsinki, Finland.

E-mail: mailto:[email protected]

Web End =liisa.keltikangas-jarvinen@helsinki.

{Deceased.

Keywords: hostility; GWAS; personality; gene; development

Received 6 May 2011; accepted 10 May 2011

Genome-wide association study of hostility P Merjonen et al

2

an additional aim of the present study was to test whether genetic effects underlying hostility are stable across different ages or whether they gain importance with advanced age in adulthood.

Materials and methods

Population and study design. Participants were from the population-based prospective Young Finns (YF) cohort study, which started in 1980 with 3596 boys and girls from different geographical areas of Finland.15 The genome of the participants was genotyped in 2009 and personality tests assessing hostility were administered in four follow-up phases in 1992, 1997, 2001 and 2007. At the baseline of the present study (1992), participants were 15, 18, 21, 24, 27 and 30 years old, and they were followed up for 15 years until they were 30, 33, 36, 39, 42 and 45 years, respectively. The nal study sample with complete measurements consisted of 9821781 men and women depending on the measure of hostility (n for anger between the four measurement phases ranges between 1619 and 1776, n for cynicism between 1622 and 1781 and n for paranoia between 1622 and 1780).

Measures of hostility components. We assessed three aspects of hostility, that is, cynicism and paranoia, which both reect the cognitive component, and anger, which represents the affective component. Cynicism was measured with a seven-item self-completion cynicism scale derived from the Minnesota Multiphasic Personality Inventory (for example, It is safer to trust nobody).16,17 Paranoia was assessed with the six-item self-completion paranoid ideation subscale of the Symptom Checklist-90R (for example, Others do not give me proper credit for my achievements).18 Anger was assessed with a seven-item Irritability Scale of the Buss-Durkee Hostility Inventory (for example, I lose my temper easily but get over it quickly).19 Detailed description of the scales has been published in previous papers.9,20 Response format for all scales was on a ve-point scale, ranging from totally disagree (1) to totally agree (5), and the mean of each scale was calculated for only those who had responded to at least 50% of the items on the scale. In addition, for each scale we calculated the mean score over the four measurements in 1992, 1997, 2001 and 2007 to capture a more stable trait of hostility. Reliability for the four-measurement mean score was high, with Cronbachs a being 0.85, 0.84 and 0.82 for cynicism, paranoia and anger scales, respectively.

Genotyping and quality control of YF study. The genome-wide SNP genotyping of YF study was done by a custom Illumina BeadChip (San Diego, CA, USA) containing 670 000 SNPs and copy-number variant probes from 2442 YF participants (1123 males and 1319 females). The custom 670K chip shares 562 643 SNPs in common with the Illumina Human610 BeadChip. Genotypes were called using Illuminas clustering algorithm (Illuminus).21 A total of 2556 samples were genotyped. After initial clustering, we removed 2 subjects for poor call rates (call rate o0.90), and 54

samples failed subsequent quality control (that is, duplicated samples, heterozygosity, low call rate or custom SNP ngerprint genotype discrepancy). The following lters were applied to the remaining data: minor allele frequency 0.01, genotyped call rate (GENO) 0.05, MIND 0.05 and Hardy Weinberg equilibrium 1 10 6. Of 2500 individuals, 3 were

removed for low genotyping (MIND 40.05), 11 766 markers were excluded based on HardyWeinberg equilibrium test (Pp1 10 6), 7746 SNPs failed missingness test (GENO

40.05), 34 596 SNPs failed frequency test (minor allele frequency o0.01) and 1 individual failed gender check. None were removed by subsequent heterozygosity check. In that point, there were 546 770 SNPs and 2496 individuals who were utilized to generate an identity-by-descent matrix le in PLINK.22 There were 51 pairs of individuals with pi-hat 40.2, and thus these individuals were removed because of possible relatedness. One of the pair was removed using greater missingness as criteria. After nal frequency and genotyping running, there was 546 677 SNPs available from a sample of 2442 YF subjects. Genotype imputation was performed for the YF SNP data using MACH23 with the HapMap (phase II, release 22 CEU, NCBI build 36, dbSNP 126) haplotypes as reference.

Statistics of GWAS. Quasi-continuous mean variables of hostility subscales were BoxCox transformed. Residuals were obtained using linear regression model in which hostility variables were adjusted for sex and age in order to control the most obvious environmental factors related to hostility. Residuals were standardized (mean 0, s.d. 1) and their distributions were conrmed to be very close to normal by visual Q-Q plot analysis. We also veried that the estimates for the b-coefcients from the GWAS are not driven by few outliers by plotting leverage vs standardized residuals plots for the residuals.

We have an 80% power of identifying SNPs that explain at least 4% of the variability with sample size of 985 (mean of four measurement). For the four measurements within each hostility scale the sample sizes are increased to

B1780 subjects. These analyses were powered to detect the effects of common variants down to 2.1% of explained variability.

Tests for additive genetic effects were carried out on a linear scale using linear regression. Genotypes were coded as 0, 1 or 2 when the SNP was genotyped and by dosage (scale 02) when imputed. These tests were performed to assess association of SNPs with the standardized residuals using PLINK22 for the genotyped data. ProbABEL24 was used to t the model, taking account of the genotype uncertainty at imputed SNPs. The P-values were combined from the analysis by favoring genotyped tests over imputed ones. The Q-Q and Manhattan plots were drawn for the analysis of the results. The P-value for genome-wide signicance was set at Po9 10 8, corresponding to a target a of 0.05 with a

Bonferroni correction for 550 000 million independent tests with direct genotyping. Cynicism was normally distributed, whereas the distributions of paranoia and anger were slightly positively skewed. Thus BoxCox transformations were used for all the outcomes.

Translational Psychiatry

Genome-wide association study of hostility P Merjonen et al

3

Results

As shown in Table 1, the average age of the genotyped sample is 37.56 (s.d. 5.03). The bivariate correlations

between hostility measures are shown in Table 2. The stability of the measures (rs range 0.450.69) as well as their bivariate correlations (rs range 0.380.77) are moderate (all Ps o0.001). Cynicism and paranoia correlate higher with each other than with anger. Younger participants scored higher on the three hostility measures (r 0.12, Po0.001, r 0.08,

P 0.01, and r 0.05, P 0.123 for mean cynicism,

paranoia and anger, respectively). Females scored higher on anger (r 0.20, Po0.001) and males on cynicism

(r 0.18, Po0.001) and paranoia (r 0.09, Po0.01). All the

subsequent models were therefore adjusted for sex and age.

We tested 2 577 640 SNPs for association with the three hostility scales measured in four different time points. The top SNPs derived from SNPs with P-values p1 10 5 are

presented in the Tables 3 and 4. Table 3 shows the top SNPs when hostility is measured as a mean score of four measurement phases ((phase 1 phase 2 phase 3 phase

4/4)). Chromosome 14 at 99 cM (SNPs rs3783337, rs7158754, rs3783332, rs2181102, rs7159195, rs11160570, rs941898) predicted suggestively the mean paranoia during the 15 years at the genome-wide statistical signicance (Po9 108, Table 2). However, this suggestive association

did not replicate at each single measurement point over time (Table 4). Table 4 shows the top SNPs when the most signicant associations are selected, irrespective of measurement phase (selected from phase 1, phase 2, phase 3 or phase 4). The most signicant SNP suggestively associated with anger was found on chromosome 17 at 11 cM SNP rs11656526 (P-value o9 10 8, Table 4) for anger mea

sured in 1992. Also, loci on chromosome 6 at 6.7 cM seemed promising when predicting anger in 2007, which shows the most reliable results for anger according to Q-Q plot analyses. However, these suggestive associations did not replicate in other measurement years, and hence the stability of these associations was weak.

Table 5 shows those SNPs that replicate at different measurement phases or at different scales. The most systematic replicating evidence for suggestive genetic effects was found for cynicism, although the signicance levels (Po1 105) did not reach the Bonferroni corrected genome-

wide signicance level. Promising SNPs suggestively predicting cynicism were found on chromosome 7 at 86 cM (rs802047, rs802028, rs802030, rs802026, rs802036, rs802025, rs802024, rs802032, rs802049, rs802051), which replicated on two different measurements of cynicism (1992 and 1997) as well as the rst four of the SNPs on above on the mean of all four measurements of cynicism (Table 5). In addition, SNPs in chromosome 22 at 43 cM (rs7510759, rs7510924) were associated with cynicism in 1997 and the mean of all four measurements of cynicism. The genetic background of different components of hostility appears to be largely distinct from each other, although a group of SNPs from chromosome 17 at 2.8 cM

Table 1 Characteristics of the study group

Variable n % Mean (s.d.)

Range

Sex 2443Male 1123 46.0 Female 1320 54.0

Age in 2007 2443 37.56 (5.03) 3045

HostilityMean of 1992, 1997, 2001 and 2007

Anger 987 2.50 (0.58) 1.004.25 Cynicism 986 2.71 (0.59) 1.114.54 Paranoia 985 2.30 (0.53) 1.004.62 1992

Anger 1776 2.52 (0.77) 1.005.00 Cynicism 1781 2.82 (0.67) 1.005.00 Paranoia 1780 2.35 (0.64) 1.004.67 1997

Anger 1619 2.61 (0.75) 1.005.00 Cynicism 1622 2.89 (0.72) 1.004.86 Paranoia 1622 2.46 (0.64) 1.004.50 2001

Anger 1750 2.51 (0.71) 1.005.00 Cynicism 1740 2.70 (0.70) 1.005.00 Paranoia 1739 2.31 (0.64) 1.005.00 2007

Anger 1738 2.42 (0.68) 1.004.86 Cynicism 1737 2.51 (0.71) 1.004.71 Paranoia 1737 2.14 (0.64) 1.004.83

Table 2 Correlations between hostility measures (n ranges between 983 and 2443 from mean cynicismmean paranoia correlation with agesex correlation)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16

1. Sex2. Age in 2007 0.00

3. Mean anger 0.20*** 0.05

4. Mean cynicism 0.18*** 0.12*** 0.44***

5. Mean paranoia 0.09** 0.08** 0.57*** 0.77***

6. Anger in 1992 0.18*** 0.14*** 0.77*** 0.36*** 0.47***

7. Cynicism in 1992 0.15*** 0.16*** 0.33*** 0.76*** 0.56*** 0.41***

8. Paranoia in 1992 0.02 0.16*** 0.48*** 0.59*** 0.78*** 0.56*** 0.63***

9. Anger in 1997 0.19*** 0.02 0.83*** 0.36*** 0.47*** 0.52*** 0.25*** 0.35***

10. Cynicism in 1997 0.08** 0.11*** 0.37*** 0.85*** 0.65*** 0.27*** 0.56*** 0.45*** 0.38***

11. Paranoia in 1997 0.03 0.06* 0.47*** 0.66*** 0.84*** 0.36*** 0.44*** 0.56*** 0.49*** 0.69***

12. Anger in 2001 0.18*** 0.01 0.84*** 0.35*** 0.45*** 0.49*** 0.21*** 0.30*** 0.62*** 0.30*** 0.37***

13. Cynicism in 2001 0.11*** 0.03 0.40*** 0.87*** 0.67*** 0.30*** 0.52*** 0.44*** 0.30*** 0.66*** 0.53*** 0.38***

14. Paranoia in 2001 0.06* 0.00 0.47*** 0.66*** 0.85*** 0.32*** 0.39*** 0.50*** 0.36*** 0.51*** 0.62*** 0.46*** 0.72***15. Anger in 2007 0.12*** 0.07** 0.79*** 0.38*** 0.46*** 0.45*** 0.20*** 0.32*** 0.50*** 0.26*** 0.34*** 0.62*** 0.33*** 0.37***

16. Cynicism in 2007 0.17*** 0.04 0.37*** 0.84*** 0.68*** 0.25*** 0.46*** 0.42*** 0.28*** 0.60*** 0.54*** 0.30*** 0.69*** 0.56*** 0.39***

17. Paranoia in 2007 0.10*** 0.02 0.45*** 0.63*** 0.83*** 0.28*** 0.34*** 0.49*** 0.33*** 0.44*** 0.59*** 0.33*** 0.54*** 0.65*** 0.46*** 0.73***

*Po0.05; **Po0.01; ***Po0.001.

Sex: 1 male, 0 female.

Translational Psychiatry

Genome-wide association study of hostility

P Merjonen et al

4

Table 3 Genetic markers showing top 10 SNPs within mean of four measurement years in each hostility scale

CHR SNP Base pair Minor allele

MAF n b s.e. P-value r2 Closest gene

Anger

2 rs2882650 6 517 472 C 0.34 986 0.22 0.05 4.3 10 6 0.02

2 rs10929436 6 518 881 T 0.34 986 0.22 0.05 4.3 10 6 0.02

2 rs4668497 6 517 422 T 0.34 986 0.22 0.05 4.4 10 6 0.02

2 rs7593230 6 519 359 T 0.34 986 0.22 0.05 4.4 10 6 0.02

4 rs4859315 32 867 764 C 0.06 986 0.83 0.18 4.6 10 6 0.02

8 rs17648656 30 973 921 T 0.47 986 0.21 0.05 4.6 10 6 0.02 PURG

8 rs11776713 30 981 149 T 0.47 986 0.21 0.05 4.6 10 6 0.02 PURG

8 rs11779521 30 983 843 T 0.47 986 0.21 0.05 4.6 10 6 0.02 PURG

8 rs11775287 30 983 881 C 0.47 986 0.21 0.05 4.7 10 6 0.02 PURG

2 rs10929438 6 522 878 A 0.34 986 0.21 0.05 4.7 10 6 0.02 Cynicism

7 rs802047 86 795 721 C 0.12 985 0.35 0.07 2.9 10 7 0.03

20 rs2426192 48 695 861 A 0.29 985 0.24 0.05 6.3 10 7 0.02 FAM65C

14 rs1884535 94 806 200 A 0.05 983 0.53 0.11 6.3 10 7 0.02 CLMN

20 rs2245361 48 695 563 C 0.29 985 0.24 0.05 6.5 10 7 0.02 FAM65C

22 rs8136107 35 697 254 A 0.10 985 0.41 0.08 8.3 10 7 0.02

8 rs7833231 4 669 830 G 0.30 982 0.24 0.05 9.8 10 7 0.02 CSMD1

22 rs16997638 35 734 113 C 0.11 985 0.34 0.07 1.6 10 6 0.02 TST

7 rs802030 86 831 487 G 0.10 985 0.38 0.08 2.0 10 6 0.02 CROT

7 rs802028 86 829 611 T 0.10 985 0.38 0.08 2.0 10 6 0.02 CROT

10 rs10510007 116 626 711 G 0.33 983 0.22 0.05 2.1 10 6 0.02 FAM160B1 Paranoia14 rs3783337 99 665 031 T 0.17 984 0.34 0.06 3.5 10 8 0.03 EVL

14 rs7158754 99 653 102 A 0.17 984 0.34 0.06 3.5 10 8 0.03 EVL

14 rs3783332 99 656 510 A 0.17 984 0.34 0.06 3.5 10 8 0.03 EVL

14 rs2181102 99 653 702 G 0.17 984 0.34 0.06 3.5 10 8 0.03 EVL

14 rs7159195 99 653 083 G 0.17 984 0.34 0.06 3.5 10 8 0.03 EVL

14 rs11160570 99 651 389 T 0.17 984 0.34 0.06 3.5 10 8 0.03 EVL

14 rs941898 99 669 190 G 0.17 984 0.34 0.06 3.8 10 8 0.03 EVL

14 rs941900 99 673 152 C 0.19 984 0.28 0.06 1.3 10 6 0.02 EVL

22 rs7510759 43 038 359 A 0.16 984 0.50 0.1 1.6 10 6 0.02 KIAA1644

22 rs7510924 43 039 988 T 0.16 984 0.50 0.1 1.6 10 6 0.02 KIAA1644

Abbreviations: CHR, chromosome; MAF, minor allele frequency; SNP, single-nucleotide polymorphism.

Bold values Po9 10 8.

(rs12936442, rs894664, rs6502671, rs7216028) and from chromosome 22 at 43 cM (rs7510759, rs7510924, rs7290560) and at 36 cM (rs8136107) were suggestively associated with both cynicism and paranoia. Replications of the genetic linkage between different measurement of hostility and different measurement years are presented in Table 5.

Discussion

Our study reports results of a large-scale GWA analysis of hostility, with hostility measured in four follow-ups across 15 years of time span with three different scales. Although only few associations achieved genome-wide signicance, many associations approached signicance. We attempted to capture more reliable ndings of the genotype over time by using the mean of hostility levels in the four time points as the outcome. Most of the suggestive associations did not replicate across measurement times, which undermines the robustness of the single signicant associations. These suggestive associations should therefore be interpreted with appropriate caution. The inconsistent ndings resemble those from many previous GWA studies of personality, most of which have not yet found robust evidence for specic candidate genes.13,2527

The strongest associations were found for mean score of paranoia with a number of closely linked SNPs in chromosome

14 at 99 cM, although this suggestive association had limited replicability over time. Chromosome 14 atB100 cM has been previously linked to neuroticism and anxiety27,28 and at

103 cM to bipolar disorder.29 The present study thus adds evidence that this region may include genetic markers or determinants for general anxiety and distrust (that is, paranoia) as well as for susceptibility to psychiatric diagnoses involving distrust against others. The nding that the mean paranoia for four different time points had signicant genetic linkage, but single measurements of paranoia did not, may imply that paranoia as a stable trait has wider genetic basis, but high distrust in one point in time may depend more on transient environmental factors and be more prone to uctuate. The closest gene for the found paranoia linked SNPs is EVL gene in chromosome 14, which is proposed as a possible candidate gene for colorectal cancer.30

Another signicant genomic region found in the current study is in chromosome 17 at 2.8 cM, which was suggestively linked with both paranoia and cynicism in the most recent measurement when the participants were at age 3045 years. The closest gene for this region is RAP1GAP2, which affects GTPase-activating protein, has a role in regulating the platelet aggregation, and is expressed especially in heart, testis and blood leukocytes, and also in stomach, pancreas and intestines, and slightly in brain.31 Thus, this might

Translational Psychiatry

Genome-wide association study of hostility P Merjonen et al

5

Table 4 Genetic markers showing top 10 SNPs within each hostility scale across four measurements

CHR SNP Base pair Minor allele

MAF n b s.e. P-value r2 Closest gene

Year

Anger

17 rs11656526 11 289 530 T 0.04 1775 0.58 0.11 7.3 10 8 0.02 SHISA6 1992

6 rs17647258 67 181 498 A 0.14 1737 0.26 0.05 1.3 10 7 0.02 2007

6 rs9445708 67 181 730 A 0.14 1737 0.26 0.05 1.3 10 7 0.02 2007

6 rs9445711 67 182 360 C 0.14 1737 0.26 0.05 1.3 10 7 0.02 2007

6 rs10223593 67 183 069 G 0.14 1737 0.26 0.05 1.3 10 7 0.02 2007

6 rs10223721 67 183 177 A 0.14 1737 0.26 0.05 1.3 10 7 0.02 2007

6 rs17647306 67 183 568 G 0.14 1737 0.26 0.05 1.3 10 7 0.02 2007

6 rs10223766 67 183 733 T 0.14 1737 0.26 0.05 1.3 10 7 0.02 2007

6 rs10223625 67 183 675 G 0.14 1737 0.26 0.05 1.3 10 7 0.02 2007

6 rs10223661 67 184 072 C 0.14 1737 0.26 0.05 1.3 10 7 0.02 2007 Cynicism7 rs802047 86 795 721 C 0.12 1780 0.28 0.05 5.1 10 8 0.02 1992

7 rs802047 86 795 721 C 0.12 1621 0.27 0.05 2.6 10 7 0.02 1997

7 rs802028 86 829 611 T 0.10 1780 0.30 0.06 2.7 10 7 0.01 CROT 1992

7 rs802030 86 831 487 G 0.10 1780 0.30 0.06 2.7 10 7 0.01 CROT 1992

7 rs802026 86 826 975 A 0.10 1780 0.28 0.05 2.9 10 7 0.01 CROT 1992

7 rs802026 86 826 975 A 0.10 1621 0.28 0.06 4.7 10 7 0.02 CROT 1997

7 rs802028 86 829 611 T 0.10 1621 0.30 0.06 6.6 10 7 0.02 CROT 1997

14 rs1884535 94 806 200 A 0.05 1737 0.40 0.08 6.7 10 7 0.01 CLMN 2001

7 rs802030 86 831 487 G 0.10 1621 0.30 0.06 6.8 10 7 0.02 CROT 1997

9 rs17320021 113 380 854 G 0.05 1779 0.38 0.08 8.2 10 7 0.01 LTB4DH 1992 Paranoia12 rs10506598 69 302 206 G 0.28 1734 0.19 0.04 2.3 10 7 0.02 PTPRB 2007

14 rs2281515 92 473 316 T 0.34 1738 0.21 0.04 7.5 10 7 0.01 ITPK1 2001

13 rs9592675 69 535 509 T 0.42 1779 0.16 0.03 1.1 10 6 0.01 KLHL1 1992

19 rs11671165 43 962 864 G 0.23 1621 0.21 0.04 1.2 10 6 0.01 1997

13 rs9317872 69 538 646 C 0.42 1777 0.16 0.03 1.2 10 6 0.01 KLHL1 1992

13 rs12853326 69 533 917 G 0.42 1779 0.16 0.03 1.9 10 6 0.01 KLHL1 1992

20 rs348790 58 487 910 C 0.40 1779 0.17 0.03 2.0 10 6 0.01 1992

20 rs17724512 6 612 177 A 0.07 1779 0.33 0.07 2.0 10 6 0.01 1992

18 rs10514232 72 987 524 G 0.15 1736 0.23 0.05 2.2 10 6 0.01 2007

13 rs12871523 69 532 587 A 0.41 1779 0.16 0.03 2.2 10 6 0.01 KLHL1 1992

Abbreviations: CHR, chromosome; MAF, minor allele frequency; SNP, single-nucleotide polymorphism.

Bold values Po9 10 8.

Table 5 SNPs replicating in different years or different hostility scales

CHR SNP BP Minor allele (A1)

MAF P-value Replication Closest gene

7 rs802047 86 795 721 C 0.12 o3 10 7 Cynicism in 1992, 1997 and mean cynicism score

7 rs802028 86 829 611 T 0.10 o2 10 6 Cynicism in 1992, 1997 and mean cynicism score CROT

7 rs802030 86 831 487 G 0.10 o2 10 6 Cynicism in 1992, 1997 and mean cynicism score CROT

7 rs802026 86 826 975 A 0.10 o4 10 6 Cynicism in 1992, 1997 and mean cynicism score CROT

7 rs802036 86 815 830 G 0.09 o7 10 6 Cynicism in 1992 and 1997 CROT

7 rs802025 86 824 568 T 0.07 o1 10 5 Cynicism in 1992 and 1997 CROT

7 rs802024 86 823 655 T 0.07 o1 10 5 Cynicism in 1992 and 1997 CROT

7 rs802032 86 801 186 A 0.07 o1 10 5 Cynicism in 1992 and 1997

7 rs802049 86 797 791 T 0.07 o1 10 5 Cynicism in 1992 and 1997

7 rs802051 86 798 396 T 0.07 o1 10 5 Cynicism in 1992 and 1997

17 rs12936442 2 879 859 A 0.10 o6 10 6 Cynicism and paranoia in 2007 RAP1GAP2

17 rs894664 2 857 234 A 0.13 o8 10 6 Cynicism and paranoia in 2007 RAP1GAP2

17 rs6502671 2 852 848 A 0.13 o7 10 6 Cynicism and paranoia in 2007 RAP1GAP2

17 rs7216028 2 880 423 T 0.11 o8 10 6 Cynicism and paranoia in 2007 RAP1GAP2

22 rs7510759 43 038 359 A 0.16 o5 10 6 Cynicism in 1997 and mean cynicism and mean paranoia score KIAA1644

22 rs7510924 43 039 988 T 0.16 o5 10 6 Cynicism in 1997 and mean cynicism and mean paranoia score KIAA1644

22 rs8136107 35 697 254 A 0.10 o5 10 6 Mean cynicism and paranoia score

22 rs7290560 43 036 573 A 0.15 o6 10 6 Mean cynicism and paranoia score KIAA1644

Abbreviations: BP, base pair; CHR, chromosome; MAF, minor allele frequency; SNP, single-nucleotide polymorphism.

also be a possible link between hostility and health problems. Both cynicism and paranoia mean scores were also associated with a region in chromosome 22 at 36 cM and at 43 cM for which the closest gene is KIAA1644. Neither RAP1GAP2 nor KIAA1644 have previously been

linked to personality traits, although chromosome 22 at 36 cM has been linked to bipolar disorder and schizophrenia.32

Although the results for cynicism did not reach the

Bonferroni corrected statistical signicance level, there were

Translational Psychiatry

Genome-wide association study of hostility P Merjonen et al

6

many marginally signicant associations. Especially, areas on chromosome 7 at 86 cM were related to cynicism in 1992, 1997 and mean of cynicism measurements. The nearest gene for this region is CROT (carnitine O-octanoyltransferase) that affects fatty acid functioning in a cell level and is expressed at least in mice almost everywhere in the body, especially in liver and intestines, and also slightly in heart and brain.33

The observed suggestive associations may have some theoretical implications. Theoretically, cynicism is assumed to develop more in response to environmental experiences, which may explain the less signicant relation to genetic background. However, it may be that there are multiple overlapping genetic effects and gene gene

and gene environment interactions that prevent SNPs to

reach the Bonferroni corrected signicance level. Same locations in chromosome 17 at 2.8 cM and chromosome 22 at 36 and 43 cM were associated with both cynicism and paranoia, which may imply shared genetic background with these hostility dimensions. Such hostile attitudes might be seen as core of the hostility construct.1,7,9 Anger, on the

other hand, is theoretically a separate construct having its developmental roots in temperament-like characteristics.34

Our nding that anger did not share similar genetic background with cynicism or paranoia implies that the consideration of anger as a separate construct seems justied also from the genetic perspective. The phenotypic and genotypic differences behind hostility measures may thus in part explain the mixed ndings between hostility and cardiovascular health.8

Measuring complex personality traits, like hostility, involves challenges of accurate measurement of the phenotype. Measurement error due to imperfect assessment of the phenotype reduces the ability to capture stabile phenotype over time (testretest correlations) and introduces time-specic variance in the measures. The lack of adjustment for relevant environmental factors inuencing personality development may partly explain why the GWAS ndings of personality traits rarely replicate in different time points or in different samples. This is not a unique problem to our study, as previous studies with well-established personality scales, for example, Temperament and Character Inventory25,35 and

Big Five,26,36,37 have rarely found consistent associations

with GWAS.

In summary, this GWAS showed preliminary evidence for specic regions possibly related to hostility. The suggestive associations were small in magnitude and did not replicate across all measurement times, and thus they warrant further study in other populations. Single SNPs are likely to have small and thereby variable effects on personality traits, and many real effects may be lost in plenty of associations because of insufcient statistical power and measurement imprecision related to the identication of the phenotype. Accumulating evidence from several cohorts should provide more accurate and reliable data on the genetic background of hostility and other personality traits.

Conict of interest

The authors declare no conict of interest.

Acknowledgements. We thank Ville Aalto and Irina Lisinen for their valuable help in data management. This study was supported by the Academy of Finland (Grant No. 123621 (LP-R)), the University of Helsinki Research Funds (Project No. 2106012 (LP-R)), the Turku University Hospital research grants (OTR), Tampere University Hospital Medical Fund (TL), Emil Aaltonen Foundation (TL), the Signe and Ane Gyllenberg Foundation (LK-J), the Alfred Kordelin Foundation (PM), the Otto A Malm Foundation (PM), the Jenny and Antti Wihuri Foundation (PM) and the National Doctoral Programme of Psychology in Finland (PM). MiKi is supported by the National Heart Lung and Blood Institute (R01HL36310), US, NIH: National Institute on Aging (R01AG013196; R01AG034454), US, NIH; the Academy of Finland, EU New OSH ERA research program, and the BUPA Foundation Specialist Research Grant, UK. The Young Finns study is also nancially supported by the Academy of Finland (Grant No. 77841, 210283, 117832 and 121584), the Social Insurance Institution of Finland, Juho Vainio Foundation, Finnish Foundation of Cardiovascular Research and Finnish Cultural Foundation.

Disclaimer. The funding sources had neither been involved in the submission of the manuscript nor in the decision to publish the data.

1. Greenglass ER, Julkunen J. Construct validity and sex differences in Cook-Medley hostility. Person Individ Diff 1989; 10: 209218.

2. Soyka M, Graz C, Bottlender R, Dirschedl P, Schoech H. Clinical correlates of later violence and criminal offences in schizophrenia. Schizophr Res 2007; 94: 8998.

3. Abu-Akel A, Abushualeh K. Theory of mind in violent and nonviolent patients with paranoid schizophrenia. Schizophr Res 2004; 69: 4653.

4. Vandervoort D. Quality of social support in mental and physical health. Curr Psychol 1999; 18: 205222.

5. Heponiemi T, Elovainio M, Kivimaki M, Pulkki L, Puttonen S, Keltikangas-Jarvinen L. The longitudinal effects of social support and hostility on depressive tendencies. Soc Sci Med 2006; 63: 13741382.

6. Chida Y, Steptoe A. The association of anger and hostility with future coronary heart disease: a meta-analytic review of prospective evidence. J Am Coll Cardiol 2009; 53: 936946.

7. Miller TQ, Smith TW, Turner CW, Guijarro ML, Hallet AJ. A meta-analytic review of research on hostility and physical health. Psychol Bull 1996; 119: 322348.

8. Rozanski A, Blumenthal JA, Kaplan J. Impact of psychological factors on the pathogenesis of cardiovascular disease and implications for therapy. Circulation 1999; 99: 21922217.

9. Merjonen P, Pulkki-Rback L, Lipsanen J, Lehtimaki T, Rontu R, Viikari J et al. Development of adulthood hostile attitudes: childhood environment and serotonin receptor gene interactions. Pers Relationship 2011; 18: 184197.

10. Hur YM. Nonadditive genetic effects on hostility in South Korean adolescent and young adult twins. Twin Res Hum Genet 2006; 9: 637641.

11. Rebollo I, Boomsma DI. Genetic analysis of anger: genetic dominance or competitive sibling interaction. Behav Genet 2006; 36: 216228.

12. Coccaro EF, Bergeman CS, Kavoussi RJ, Seroczynski AD. Heritability of aggression and irritability: a twin study of the Buss-Durkee aggression scales in adult male subjects. Biol Psychiatry 1997; 41: 273284.

13. Knox SS, Wilk JB, Zhang Y, Weidner G, Ellison RC. A genome scan for hostility: the National Heart, Lung, and Blood Institute Family Heart Study. Mol Psychiatry 2004; 9: 124126.

14. Gustavsson JP, Pedersen NL, Asberg M, Schalling D. Exploration into the sources of individual differences in aggression-, hostility- and anger-related (AHA) personality traits. Pers Individ Diff 1996; 21: 10671071.

15. Raitakari OT, Juonala M, Rnnemaa T, Keltikangas-Jarvinen L, Rasanen L, Pietikainen M et al. Cohort prole: the Cardiovascular Risk in Young Finns Study. Int J Epidemiol 2008; 37: 12201226.

16. Comrey AL. A factor analysis of items on the MMPI paranoia scale. Educ Psychol Meas 1958; 18: 99107.

17. Comrey AL. A factor analysis of items on the MMPI hysteria scale. Educ Psychol Meas 1957; 17: 586592.

18. Derogatis LR, Cleary PA. Conrmation of the dimensional structure of the SCL-90: a study in construct validation. J Clin Psychol 1977; 33: 981989.

19. Buss AH, Durkee A. An inventory for assessing different kinds of hostility. J Consult Psychol 1957; 21: 343349.

20. Merjonen P, Pulkki-Rback L, Puttonen S, Keskivaara P, Juonala M, Telama R et al. Anger is associated with subclinical atherosclerosis in low SES but not in higher SES men and women. The Cardiovascular Risk in Young Finns Study. J Behav Med 2008; 31: 3544.

21. Teo YY, Inouye M, Small KS, Gwilliam R, Deloukas P, Kwiatkowski DP et al. A genotype calling algorithm for the Illumina BeadArray platform. Bioinformatics 2007; 23: 27412746.

22. 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.

Translational Psychiatry

Genome-wide association study of hostility P Merjonen et al

7

23. Li Y, Abecasis GR. Mach 1.0: rapid haplotype reconstruction and missing genotype inference. Am J Hum Genet 2006; S79: 2290.

24. Aulchenko YS, Ripke S, Isaacs A, van Duijn CM. GenABEL: an R library for genome-wide association analysis. Bioinformatics 2007; 23: 12941296.

25. Verweij KJH, Zietsch BP, Medland SE, Gordon SD, Benyamin B, Nyholt DR et al. A genome-wide association study of Cloningers temperament scales: implications for the evolutionary genetics of personality. Biol Psychol 2010; 85: 306317.

26. Terracciano A, Sanna S, Uda M, Deiana B, Usala G, Busonero F et al. Genome-wide association scan for ve major dimensions of personality. Mol Psychiatry 2010; 15: 647656.

27. Wray NR, Middeldorp CM, Birley AJ, Gordon SD, Sullivan PF, Visscher PM et al. Genome-wide linkage analysis of multiple measures of neuroticism of 2 large cohorts from Australia and the Netherlands. Arch Gen Psychiatry 2008; 65: 649658.

28. Middeldorp CM, Hottenga JJ, Slagboom PE, Sullivan PF, de Geus EJ, Posthuma D et al. Linkage on chromosome 14 in a genome-wide linkage study of a broad anxiety phenotype. Mol Psychiatry 2008; 13: 8489.

29. Burton PR, Clayton DG, Cardon LR, Craddock N, Deloukas P, Duncanson A et al. Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls. Nature 2007; 447: 661678.

30. Sjblom T, Jones S, Wood LD, Parsons DW, Lin J, Barber TD et al. The consensus coding sequences of human breast and colorectal cancers. Science 2006; 314: 268274.

31. Schultess J, Danielewski O, Smolenski AP. Rap1GAP2 is a new GTPase-activating protein of Rap1 expressed in human platelets. Blood 2005; 105: 31853192.

32. Badner JA, Gershon ES. Meta-analysis of whole-genome linkage scans of bipolar disorder and schizophrenia. Mol Psychiatry 2002; 7: 405411.

33. Westin MA, Hunt MC, Alexson SE. Short- and medium-chain carnitine acyltransferases and acyl-CoA thioesterases in mouse provide complementary systems for transport of beta-oxidation products out of peroxisomes. Cell Mol Life Sci 2008; 65: 982990.

34. Nigg JT. Temperament and developmental psychopathology. J Child Psychol Psychiatry 2006; 47: 395422.

35. Cloninger CR, Van Eerdewegh P, Goate A, Edenberg HJ, Blangero J, Hesselbrock V et al. Anxiety proneness linked to epistatic loci in genome scan of human personality traits. Am J Med Genet 1998; 81: 313317.

36. Terracciano A, Tanaka T, Sutin AR, Sanna S, Deiana B, Lai S et al. Genome-wide association scan of trait depression. Biol Psychiatry 2010; 68: 811817.

37. de Moor MHM, Costa PT, Terracciano A, Krueger RF, de Geus EJC, Toshiko T et al. Meta-analysis of genome-wide association studies for personality. Mol Psychiatry; e-pub ahead of print 21 December 2010; doi:10.1038/mp.2010.128 (in press).

Translational Psychiatry is an open-access journal published by Nature Publishing Group. This work is licensed under the Creative Commons Attribution-Noncommercial-No Derivative Works 3.0 Unported License. To view a copy of this license, visit http://creativecommons.org/licenses/by-nc-nd/3.0/

Translational Psychiatry