Genome-wide association study identifies seven novel loci associating with circulating cytokines and cell adhesion molecules in Finns

607 Sliz E, etal. J Med Genet 2019;56:607–616. doi:10.1136/jmedgenet-2018-105965 Original article Genome-wide association study identifies seven novel loci associating with circulating cytokines and cell adhesion molecules inFinns eeva Sliz, 1,2,3 Marita Kalaoja,1,2,3 ari ahola-Olli,4,5 Olli raitakari,5,6 Markus Perola,7,8,9 Veikko Salomaa,7 terho lehtimäki,10 toni Karhu,3,11 Heimo Viinamäki,12 Marko Salmi,13 Kristiina Santalahti,13 Sirpa Jalkanen,13 Jari Jokelainen,2,14 Sirkka Keinänen-Kiukaanniemi,2,14,15 Minna Männikkö,16 Karl-Heinz Herzig,3,11,17,18 Marjo-riitta Järvelin,2,3,19,20 Sylvain Sebert,2,3,21 Johannes Kettunen1,2,3 Complex traits To cite: Sliz e, Kalaoja M, ahola-Olli a, etal. J Med Genet 2019;56:607–616. ►additional material is published online only. to view please visit the journal online (http:// dx. doi. org/ 10. 1136/ jmedgenet2018105965). For numbered affiliations see end of article. Correspondence to Dr Johannes Kettunen, center for life course Health research, Oulun Yliopisto, Oulu 90014, Finland; johannes. kettunen@ oulu. fi SS and JK contributed equally. received 27 December 2018 revised 10 april 2019 accepted 20 april 2019 Published Online First 19 June 2019 © author(s) (or their employer(s)) 2019. re-use permitted under cc BY-nc. no commercial re-use. See rights and permissions. Published by BMJ. AbsTrACT background inflammatory processes contribute to the pathophysiology of multiple chronic conditions. genetic factors play a crucial role in modulating the inflammatory load, but the exact mechanisms are incompletely understood. Objective to assess genetic determinants of 16 circulating cytokines and cell adhesion molecules (inflammatory phenotypes) in Finns. Methods genome-wide associations of the inflammatory phenotypes were studied in northern Finland Birth cohort 1966 (n=5284). a subsequent meta-analysis was completed for 10 phenotypes available in a previous genome-wide association study, adding up to 13 577 individuals in the study. complementary association tests were performed to study the effect of the aBO blood types on soluble adhesion molecule levels. results We identified seven novel and six previously reported genetic associations (p<3.1×10−9). three loci were associated with soluble vascular cell adhesion molecule-1 (sVcaM-1) level, one of which was the ABO locus that has been previously associated with soluble e-selectin (se-selectin) and intercellular adhesion molecule-1 (sicaM-1) levels. Our findings suggest that the blood type B associates primarily with sVcaM-1 level, while the a1 subtype shows a robust effect on seselectin and sicaM-1 levels. the genotypes in the ABO locus associating with higher soluble adhesion molecule levels tend to associate with lower circulating cholesterol levels and lower cardiovascular disease risk. Conclusion the present results extend the knowledge about genetic factors contributing to the inflammatory load. Our findings suggest that two distinct mechanisms contribute to the soluble adhesion molecule levels in the ABO locus and that elevated soluble adhesion molecule levels per se may not increase risk for cardiovascular disease. InTrOduCTIOn It is currently established that inflammatory load may play a role in the aetiology of autoimmune and infectious diseases, but also in a broad range of other diseases, such as chronic cardiometabolic disorders,1 neurodegenerative diseases2 and cancer.3 The risk for these diseases increases with age,4 and due to the world’s ageing population5 their prevalence is likely to expand. Moreover, these diseases often co-occur, which is likely due to shared inflammation-related pathophysiology.6 Inflammation is the body’s physiological response to harmful stimuli involving multiple molecular and cellular interactions attempting to restore disturbances in tissue or systemic homeostasis. Circulating cytokines, growth factors, chemokines and cell adhesion molecules (CAMs) (hereafter inflammatory phenotypes) are fundamental mediators of inflammatory responses. Genes encoding these molecules and their receptors play a crucial role in mediating the related functions. Previous studies have identified loci associating with levels of inflammatory phenotypes,7–9 but the understanding of the exact regulatory mechanisms is still incomplete. To add insights into the genetic mechanisms contributing to the inflammatory load, we performed a genome-wide association study (GWAS) of 16 circulating inflammatory phenotypes in 5284 individuals from Northern Finland Birth Cohort 1966 (NFBC1966) and a subsequent meta-analysis of 10 phenotypes in three other Finnish population cohorts,7 adding up to a total of 13 577 individuals in the study. We report identification of seven novel and replication of six loci associating with levels of the circulating inflammatory markers. MeThOds study populations, genotyping and inflammatory phenotype quantification Northern Finland Birth Cohort 1966 The NFBC1966 comprises 96% of all births during 1966 in the two northernmost provinces in Finland; altogether 12 058 children were live-born into the cohort, and the follow-ups occurred at the ages of 1, 14, 31 and 46 years.10 The data analysed in the present study are from the 31 years’ follow-up when clinical examinations and blood sampling were completed for altogether 6033 individuals, 5284 of whom had body mass index (BMI), inflammatory phenotypes and genotype data available (a maximum number of individuals per inflammatory marker of 5100). Genotyping of the samples was completed of Tampere. Protected by copyright. on October 22, 2019 at Tampere University and University Hospitalhttp://jmg.bmj.com/J Med Genet: first published as 10.1136/jmedgenet-2018-105965 on 19 June 2019. Downloaded from 608 Sliz E, etal. J Med Genet 2019;56:607–616. doi:10.1136/jmedgenet-2018-105965 Complex traits using 370 k Illumina HumanHap arrays (Illumina, California, USA), and subsequent imputation was performed based on the 1000 Genomes reference panel. A total of 16 inflammatory phenotypes were quantified from overnight fasting plasma samples using Bio-Rad’s Bio-Plex 200 system (Bio-Rad Laboratories, California, USA) with Milliplex Human Chemokine/Cytokine and CVD/Cytokine kits (Cat# HCYTOMAG-60K-12 and Cat# SPR349; Millipore, St Charles, Missouri, USA) and Bio-Plex Manager Software V.4.3 as previously described.11 The 16 inflammatory phenotypes studied in the NFBC1966 were interleukin (IL) 1-alpha, IL1-beta (IL1β), IL4, IL6, IL8, IL17, IL1 receptor antagonist (IL1ra), interferon gamma-induced protein 10 (IP10), monocyte chemoattractant protein 1 (MCP1), tumour necrosis factor alpha (TNFα), vascular endothelial growth factor (VEGF), plasminogen activator inhibitor 1, soluble CD40 ligand, soluble E-selectin (sE-selectin), soluble intercellular adhesion molecule-1 (sICAM-1) and soluble vascular cell adhesion molecule 1 (sVCAM-1). GWAS summary statistics from three Finnish population cohorts Meta-analyses were conducted for 10 phenotypes available in a previous GWAS.7 The study included up to 8293 Finnish individuals from the Cardiovascular Risk in Young Finns Study (YFS)12 and FINRISK ( www. thl. fi/ finriski),13 adding up to 13 577 individuals studied in the present meta-analyses. Shortly, YFS is a population-based follow-up study started in 1980 comprising randomly chosen individuals from Finnish cities Helsinki, Kuopio, Tampere, Oulu and Turku. The YFS data included in the previous GWAS are from 2019 individuals who participated in the follow-up in 2007 and who had both inflammatory phenotype and genotype data available. FINRISK is a Finnish population survey conducted every 5 years to monitor chronic diseases and their risk factors. The surveys use independent, random and representative samples from different geographical areas of Finland. The data included in the present meta-analyses were from participants of the 1997 and 2002 surveys. Genotypes were obtained using 670 k Illumina HumanHap arrays (Illumina) and imputed based on 1000 Genomes reference panel. Inflammatory markers were quantified using Bio-Rad’s premixed Bio-Plex Pro Human Cytokine 27-plex Assay and 21-plex Assay, and Bio-Plex 200 reader with Bio-Plex V.6.0 software (Bio-Rad Laboratories) as previously described.14 Samples were serum in YFS, EDTA plasma in FINRISK1997 and heparin plasma in FINRISK2002. statistical analyses GWAS and meta-analysis To allow meta-analysis between the present results and the previous GWAS, the data processing and analysis model were done according to Ahola-Olli et al.7 First, rank-based inverse transformation was applied to normalise the phenotypes. Preceding the analyses, linear regression models were fitted to adjust the transformed inflammatory phenotypes for age, sex, BMI and the 10 first genetic principal components to control for population stratification. The resulting residuals were again normalised with inverse transformation, and the adjusted and transformed residuals were used as phenotypes in the analyses. Genome-wide association tests were performed using snptest V.2.5.1 software.15 Allele effects were estimated using an additive model (-frequentist 1), and the option to centre and scale the phenotypes was disabled (-use_raw_phenotypes). The GWAS results were filtered by including markers with model fit info >0.8 and minor allele count >10. Filtered data were used to perform meta-analyses by METAL software (V.2011-03-25)16 for the 10 phenotypes (IL1β, IL1ra, IL4, IL6, IL8, IL17, IP10, MCP1, TNFα and VEGF) available in the previous GWAS.7 Genomic control correction was enabled (GENOMICCONTROL ON) to account for population stratification and cryptic relatedness. To estimate the heterogeneity of effect sizes between NFBC1966 and the previous GWAS, calculation of heterogeneity statistics based on Cochrane’s Q-test was enabled (ANALYZE HETEROGENEITY). Supplemental genome-wide tests in NFBC1966 Individuals showing symptoms of an acute infection were omitted from the supplemental genome-wide tests performed in the NFBC1966 population. Here, individuals reported having fever at the time of the blood sampling and individuals having C-reactive protein (CRP) level >10 mg/L were excluded. Otherwise the analysis models were as above. Conditional analyses and variance explained To assess whether the identified loci harbour multiple independent association signals, we conducted conditional analyses by further adjusting the models with the locus-specific lead variants. The association tests were repeated within a 2 Mb window around the lead SNP for the phenotypes studied in the NFBC1966 population only. For the meta-analysed phenotypes, we applied a method proposed by Yang et al17 that enables conditional analyses of GWAS summary statistics. NFBC1966 was used as a reference sample to estimate linkage disequilibrium (LD) corrections in these analyses. The proportion of variance explained was calculated using all independent variants using the following formula: Variance explained = ( βx √ 2x MAF ( 1 − MAF ))2 Here β is the variant’s effect estimate on the inflammatory phenotype and MAF denotes minor allele frequency. Complementary association tests on soluble adhesion molecule levels Complementary association tests were conducted to better evaluate the molecular mechanism explaining the two potentially independent association signals with soluble CAM levels in the ABO locus. Here, linear models were repeated within a 2 Mb window and further adjusted for the ABO blood type or rs507666 genotype tagging the A1 subtype.18 In addition, we determined the effect estimates of ABO blood types and ABO blood types stratified by rs507666 genotype on sE-selectin, sICAM-1 and sVCAM-1 levels: the adjusted and transformed CAM concentrations were as outcomes in the linear models and ABO blood types as categorical variables (individuals with blood type A vs non-A, and so on). Corresponding models were fitted for the rs507666-stratified blood types (individuals with blood type A and rs507666 G/G vs others, and so on). Shared genetic influences on inflammatory and cardiovascular phenotypes As previous evidence suggests that elevated concentrations of circulating markers of inflammation increase the risk of cardiovascular diseases (CVD),19 20 we further evaluated how variants in the loci associating with inflammatory phenotypes may relate to other cardiovascular traits. We used the gwas-pw method developed by Pickrell et al21 that estimates whether a locus harbours a genetic variant influencing one of the two phenotypes compared (models 1 and 2), if the same variant influences both phenotypes (model 3), or if separate variants within a locus influence the two phenotypes (model 4). Using the gwas-pw and open-access data provided by CARDIoGRAM,22 MEGASTROKE consortium23 and of Tampere. Protected by copyright. on October 22, 2019 at Tampere University and University Hospitalhttp://jmg.bmj.com/J Med Genet: first published as 10.1136/jmedgenet-2018-105965 on 19 June 2019. Downloaded from 609 Sliz E, etal. J Med Genet 2019;56:607–616. doi:10.1136/jmedgenet-2018-105965 Complex traits Table 1 Basic characteristics of the Northern Finland Birth Cohort 1966 study population Characteristics Total number of individuals 5284 Number of men (%) 2543 (48.1) Age, years 31.1±0.4 Body mass index, kg/m224.4±4.0 Glucose, mmol/L 5.1±0.7 Low-density lipoprotein-cholesterol, mmol/L 3.0±0.9 High-density lipoprotein-cholesterol, mmol/L 1.6±0.4 Systolic blood pressure, mm Hg 124.2±13.6 Diastolic blood pressure, mm Hg 76.8±11.7 Values are mean±SD. Global Lipids Genetics Consortium,24 we evaluated the shared genetic determinants of circulating levels of low-density lipoprotein cholesterol (LDL-C), high-density lipoprotein cholesterol (HDL-C), total cholesterol (TotC) and total triglycerides (TG), as well as risk of coronary artery disease (CAD), ischaemic stroke and the inflammatory phenotypes showing significant genetic associations in the present study. The genome breakpoint data set for individuals with European ancestry provided at https:// bitbucket. org/ nygcresearch/ ldetectdata was used in the gwas-pw analyses to split the genome into approximately independent blocks.25 resulTs The basic characteristics of the NFBC1966 study population are provided in table 1. Inflammatory phenotype distributions are tabulated in online supplementary table S1, and their correlation structure is shown in online supplementary figure S1. Using a threshold of p<3.1×10−9 for statistical significance (standard genome-wide significance level p<5×10−8 corrected for 16 phenotypes tested), we identified seven novel and six previously reported loci associating with one or more of the inflammatory phenotypes. The results are summarised in table 2 and combined Manhattan plots are shown in figure 1. Manhattan plots and Q-Q plots for each inflammatory phenotype are provided in online supplementary figure S2 A–Z. Genomic inflation factor values range between 0.99 and 1.02, suggesting no inflation in the test statistics (online supplementary table S2). Online supplementary table S3 lists the traits associated previously with the loci showing novel associations with inflammatory phenotypes in the present study. Cell adhesion molecules The ABO locus shows large effects on sE-selectin, sICAM-1 and sVCAM-1 levels We observed a novel effect on sVCAM-1 concentration in 9q34.2 near ABO (ABO, alpha 1–3 n-acetylgalactosaminyltransferase and alpha 1–3-galactosyltransferase) in the NFBC1966 population. This locus showed a robust association also with sE-selectin and sICAM-1 concentrations as previously reported.18 26 27 Noteworthy, the lead variant for sE-selectin and sICAM-1 associations (rs2519093) was different from the lead variant for sVCAM-1 association (rs8176746). The former variant is in LD (r2=1 in NFBC1966) with rs507666 tagging the ABO blood type A subtype A1, whereas the latter variant tags the blood type B.18 As the GWAS results suggested two potentially separate association signals with the soluble CAM levels in the ABO locus, we conducted complementary association tests to better evaluate the molecular mechanisms explaining the associations. The association of the rs8176746 with sVCAM-1 concentration was significant when adjusted for the rs507666 indicative of the A1 subtype (p=4.98×10−15). On the contrary, the associations of the rs2519093 with concentrations of sE-selectin and sVCAM-1 were highly significant when adjusted for the ABO blood type (p=3.40×10−123 and p=3.43×10−17, respectively). Overall, these results suggest that the association of the rs8176746 with sVCAM-1 level is independent of the A1 subtype, while the association of the rs2519093 with sE-selectin and sICAM-1 levels is independent of the ABO blood type. Statistical significances were abolished when the rs8176746 association with sVCAM-1 was adjusted for ABO blood type and rs2519093 association with sE-selectin or sICAM-1 was adjusted for rs507666. To further evaluate the related molecular mechanisms, we determined the effect estimates of the ABO blood types and ABO blood types stratified by rs507666 genotype on soluble CAM levels. The blood type A showed negative associations with the levels of all the three CAMs and the effect was the most robust on the sE-selectin level (figure 2, left panel). However, major discrepancies in the effect directions were seen when the analyses were stratified by the rs507666 genotype (figure 2, right panel). Congruent with previous reports,18 26 the present results suggest that the A1 subtype/rs507666 influences sE-selectin or sICAM-1 levels. In contrast, the blood type B seems to attribute predominantly to sVCAM-1 level, while the A1 subtype/rs507666 shows only a modest effect on sVCAM-1. HSP90B1 and ABCA8 loci associate with sVCAM-1 levels We identified two other novel loci for sVCAM-1 (12q23.3 and 17q24.2) in the NFBC1966 population. In chr12 the lead variant rs117238625 is in LD (r2=1 in NFBC1966) with rs117468318 that locates in the 5’ untranslated region (UTR) region of HSP90B1 (heat shock protein 90 kDa beta member 1) and, according to RegulomeDB,28 is likely to affect transcription factor binding. The association signal in chr17 locates near ABCA8 (ATP binding cassette subfamily A member 8) encoding one of the ATP binding cassette transporters. Variations in sialyltransferase encoding genes show an effect on sE-selectin level For sE-selectin level, we identified a novel association in 11q24.2 in the region of ST3GAL4 (ST3 beta-galactoside alpha-2,3-sialyltransferase 4). We identified a suggestive signal with sE-selectin level also in 3q12.1 near ST3GAL6 (ST3 beta-galactoside alpha-2,3-sialyltransferase 6), but the association was not significant after multiple correction (p=1.75×10−08). Both of the sialyltransferase genes have been implicated in the production of functional E-selectin, P-selectin and L-selectin ligands in mice.29 Two independent association signals on sICAM-1 level near ICAM1 We replicated the previously reported association for sICAM-1 level in 19p13.2 near ICAM1 (intracellular adhesion molecule 1).18 30 When the primary association test was conditioned for the lead variant rs117960796, another significant association was detected (rs74428614, p=1.14×10−16) indicative of more than one independent variant contributing to sICAM-1 level in this locus. Vascular endothelial growth factor In the meta-analyses, we identified a novel locus 4p16.2 with a large effect on VEGF (β=−2.38 SD). This locus harbours genes EVC (EvC ciliary complex subunit 1), EVC2 (EvC ciliary complex subunit 2) and STK32B (serine/threonine kinase 32B). In addition, we replicated two previously reported loci associating with of Tampere. Protected by copyright. on October 22, 2019 at Tampere University and University Hospitalhttp://jmg.bmj.com/J Med Genet: first published as 10.1136/jmedgenet-2018-105965 on 19 June 2019. Downloaded from 610 Sliz E, etal. J Med Genet 2019;56:607–616. doi:10.1136/jmedgenet-2018-105965 Complex traits Table 2 Significant loci associating with the circulating inflammatory phenotypes study Marker locus Chr:Position Candidate gene nearest gene(s) Annotation dbsnP reference InFO eA eAF beta P value hetPVal Variance explained Total variance explained NFBC1966 sE-selectin 9q34.2 9:136 141 870 ABO ABO Intronic rs2519093 0.994 T 0.188 −0.903 4.48e-305 NA 0.249 0.258 11q24.2 11:126 266 665 ST3GAL4 ST3GAL4 Intronic rs11220471 0.967 G 0.212 −0.162 7.72e-12 nA 0.009 sICAM-1 9q34.2 9:136 141 870 ABO ABO Intronic rs2519093 0.994 T 0.188 −0.352 7.43e-48 NA 0.038 0.118 19p13.2 19:10 383 403 ICAM1 ICAM1 Intronic rs117960796 0.802 A 0.012 −1.669 8.03e-40 NA 0.066 19p13.2 19:10 497 360 ICAM1 CDC37 Intergenic rs74428614 0.992 A 0.163 0.226 1.14e-16* NA 0.014 sVCAM-1 9q34.2 9:136 131 322 ABO ABO Missense rs8176746 1.000 T 0.129 0.256 5.06e-19 nA 0.015 0.038 12q23.3 12:104 448 391 HSP90B1 GLT8D2 Intronic rs117238625 0.981 A 0.023 0.510 2.90e-14 nA 0.012 17q24.2 17:66 823 805 ABCA8 ABCA8 Intergenic rs112001035 0.883 A 0.060 −0.324 1.04e-13 nA 0.012 Meta-analyses IL1β6p22.1 6:30 017 071 HLA locus Intronic rs6917603 1.000 C 0.251 −0.163 1.76e-12 1.00 0.010 0.015 6p22.1 6:30 013 887 HLA locus Intronic rs9261224 1.000 T 0.035 0.261 1.31e-09 * 1.00 0.005 IP10 4q21.1 4:76 899 176 CXCL10 SAD1 Intronic rs192716315 0.851 C 0.003 1.513 2.71e-13 1.00 0.014 0.014 MCP1 1q23.2 1:159 175 354 ACKR1 ACKR1 Missense rs12075 1.000 A 0.469 0.148 1.43e-33 1.51e-13 0.011 0.011 TNFα13q14.3 13:51 141 997 DLEU1 DLEU1 Intronic rs17074575 0.803 G 0.002 2.131 2.71e-09 1.00 0.018 0.018 VEGF 4p16.2 4:5 636 073 STK32B EVC2 Intronic rs186725382 0.875 A 0.001 −2.380 4.53e-10 1.00 0.011 0.052 6p21.1 6:43 927 050 VEGFA C6orf223 Intergenic rs7767396 1.000 A 0.523 0.284 8.35e-105 1.22e-69 0.040 0.056 9p24.2 9:2 686 273 VLDLR VLDLR, KCNV2 Intergenic rs7030781 0.959 T 0.373 -0.099 1.57e-13 5.34e-04 0.005 Statistical significance is considered at p<3.1×10−9. Novel findings are highlighted with bold font. All positions correspond to human genome build 37. *Indicates associations that are significant after conditioning the analyses on the locus-specific lead variant on the preceding row. EA, effect allele; EAF, effect allele frequency; HLA, human leukocyte antigen; HetPVal, p value of heterogeneity as estimated by Cochrane’s Q-test; IL1β, interleukin 1-beta ; INFO, imputation score in NFBC1966; IP10, interferon gamma-induced protein 10; MCP1, monocyte chemoattractant protein 1; NA, not available; NA, not available; NFBC1966, Northern Finland Birth Cohort 1966; TNFα, tumour necrosis factor alpha; VEGF, vascular endothelial growth factor; dbSNP, single nucleotide polymorphism database; sE-selectin, soluble E-selectin; sICAM-1, soluble intercellular adhesion molecule -1; sVCAM-1, soluble vascular cell adhesion molecule-1. of Tampere. Protected by copyright. on October 22, 2019 at Tampere University and University Hospitalhttp://jmg.bmj.com/J Med Genet: first published as 10.1136/jmedgenet-2018-105965 on 19 June 2019. Downloaded from 611 Sliz E, etal. J Med Genet 2019;56:607–616. doi:10.1136/jmedgenet-2018-105965 Complex traits Figure 1 the combined Manhattan plots for significant associations with inflammatory markers studied in (a) northern Finland Birth cohort 1966 and in (B) meta-analyses with three other Finnish population cohorts. Significance threshold p<3.1×10−9 derives from the standard p value limit for genomewide significance p<5×10−8 corrected for 16 markers examined in the present study. novel association signals are highlighted with red font and replicated loci are marked with black font. se-selectin, soluble e-selectin; il1b, interleukin 1-beta; iP10, interferon gamma-induced protein 10; McP1, monocyte chemoattractant protein 1; sicaM-1, soluble intercellular adhesion molecule-1; sVcaM-1, soluble vascular cell adhesion molecule-1; tnFa, tumour necrosis factor alpha; VegF, vascular endothelial growth factor. of Tampere. Protected by copyright. on October 22, 2019 at Tampere University and University Hospitalhttp://jmg.bmj.com/J Med Genet: first published as 10.1136/jmedgenet-2018-105965 on 19 June 2019. Downloaded from 612 Sliz E, etal. J Med Genet 2019;56:607–616. doi:10.1136/jmedgenet-2018-105965 Complex traits Figure 2 the effects of the aBO blood types and the a1 subtype on soluble adhesion molecule levels. the effects of the aBO blood types on se-selectin, sicaM-1 and sVcaM-1 levels were evaluated in linear models, where adjusted (sex, age, body mass index and the 10 first genetic principal components) and transformed soluble adhesion molecule concentrations were used as outcomes and the aBO blood type served as categorical variable (a vs non-a, and so on). corresponding models were fitted for the aBO blood types stratified by the rs507666-a allele count (0, 1 or 2), where the a allele tags the aBO subtype a1 having enhanced glycosyltransferase activity.18 no individuals were found to have B or O blood type and one or more copies of the rs507666-a allele, and thus it was not possible to perform stratification within these blood types. se-selectin, soluble e-selectin; sicaM-1, soluble intercellular adhesion molecule-1; sVcaM-1, soluble vascular cell adhesion molecule-1. VEGF levels in 6p21.1 near VEGFA (vascular endothelial growth factor A) and in 9p24.2 near VLDLR (very-low-density lipoprotein receptor).7 Proinflammatory cytokines Locus near DLEU1 shows a large effect on TNFa We identified a novel variant with a large effect on TNFα levels (β=2.13 SD) in 13q14.3 near DLEU1 and DLEU7 (deleted in lymphocytic leukaemia 1 and 7) in the meta-analyses. The human leukocyte antigen (HLA) locus shows a small effect on IL1β A novel variant at 6p22.1 in the human leucocyte antigen locus associating with IL1β level was identified in the meta-analyses. In the conditional analyses, we observed two independent association signals at this locus (table 1, online supplementary figure S2J). The same locus and the same lead variant rs6917603 showed also a suggestive effect on IL4 level (online supplementary figure S2L), but the meta-analysed result was not significant after multiple correction (p=5.56×10−09). Chemokines We replicated previously reported loci near CXCL10 (C-X-C motif chemokine ligand 10) and ACKR1 (atypical chemokine receptor 1) associating with IP10 levels and with MCP1 levels, respectively.7 supplemental genome-wide tests in nFbC1966 Altogether 236 individuals having fever or CRP >10 mg/L were excluded from the supplemental genome-wide tests performed in the NFBC1966 population. The results of the supplementary analyses were congruent with the original findings (online supplementary table S4). Comparisons of genetic effects on inflammatory phenotypes versus other traits Elevated circulating concentrations of inflammatory markers increase the risk for CVD.19 20 We used the gwas-pw method21 to evaluate the presence of shared genetic determinants between inflammatory phenotypes showing significant genetic association in the present study and other cardiovascular health-related traits (LDL-C, HDL-C, TotC, TG, CAD risk, ischaemic stroke risk) obtained from open-access sources.22–24 Altogether 56 genomic regions showed robust statistical evidence for containing a genetic variant influencing one or more of the inflammatory phenotypes and at the same time one or more of the other traits studied (model 3 posterior probability greater than 0.99; online supplementary figure S3). The ABO locus was one of the loci harbouring variants influencing multiple traits. In this locus, we observed negative linear relationships between the SNP effects on sE-selectin and sICAM-1 levels and CAD risk, stroke risk, as well as LDL-C, HDL-C and TotC levels (figure 3). The results in the other loci are provided in online supplementary figure S3. dIsCussIOn The present study examines the genetic determinants of 16 circulating inflammatory phenotypes in 5284 individuals from Northern Finland with a subsequent meta-analysis of 10 phenotypes in 3 other Finnish populations, adding up to a total of 13 577 participants. We report seven novel and replication of six previously published genetic associations. We identified a novel association for sVCAM-1 concentration at the ABO locus. This locus was also associated with sE-selectin and sICAM-1 levels as observed previously.18 26 The present GWAS suggested two distinct association signals in the ABO locus for the sE-selectin and sICAM-1 levels versus sVCAM-1 level, and the supplementary tests provided further support for at least two mechanisms contributing to circulating concentrations of CAMs in this locus. The two mechanisms include the blood type A subtype A1, which has a robust lowering effect on sE-selectin and sICAM-1 levels,18 26 and the blood type B which seems to have an increasing effect on sVCAM-1 level. The lowering effect of the A1 subtype on sE-selectin and sICAM-1 could arise from increased glycosyltransferase activity that possibly modifies the shedding of the CAMs from the endothelium and/or their clearance rate from circulation.18 26 The underlying mechanism explaining the association between the blood type B and higher sVCAM-1 concentration remains unknown and warrants research. VCAM-1-mediated adhesion involves interaction with galectin-3, a protein that has a specificity for galactosides.31 As the B antigen holds an additional galactose monomer compared with the A and O antigens, and galectins are known to recognise blood type antigens,32 it raises the speculation that the amount of unbound sVCAM-1 in the of Tampere. Protected by copyright. on October 22, 2019 at Tampere University and University Hospitalhttp://jmg.bmj.com/J Med Genet: first published as 10.1136/jmedgenet-2018-105965 on 19 June 2019. Downloaded from 613 Sliz E, etal. J Med Genet 2019;56:607–616. doi:10.1136/jmedgenet-2018-105965 Complex traits Figure 3 SnP effects on soluble adhesion molecule levels versus other cardiovascular health-related traits in the ABO locus. the Pearson’s r of the genetic effects (Z-scores) were estimated using a set of SnPs that located in the ABO locus (defined as the lD block25 containing ABO gene in the gwas-pw21 analyses) and that were available in both the present study and open-access data sets.22–24 Positive correlations are indicated with red color, negative correlations are indicated with blue color, and correlations with p≥0.05 are left blank. the scatter plot representations as well as correlations in the other loci are shown in online supplementary figure S3. caD, coronary artery disease; HDl-c, high-density lipoprotein cholesterol; lD, linkage disequilibrium; lDl-c, low-density lipoprotein cholesterol; se-selectin; soluble e-selectin; sicaM-1, soluble intercellular adhesion molecule-1; sVcaM-1, soluble vascular cell adhesion molecule-1; tg, total triglycerides; totc, total cholesterol. circulation could be influenced by a possible competitive binding of galectin-3 with sVCAM-1 and the B antigen. To evaluate the shared genetic mechanisms, we compared the correspondence of genetic effects on inflammatory phenotypes versus cardiovascular health-related traits. We observed a negative relationship between the genetic effects on CAM levels and the genetic effects on LDL-C and TotC levels, as well as lower risk for CAD and ischaemic stroke, in the ABO locus. This denotes that the genotypes in the ABO locus associating with higher levels of soluble CAMs tended to associate with lower circulating cholesterol levels as well as lower risk of cardiovascular outcomes. This was unexpected since previous evidence suggests that increased soluble CAM levels are linked with atherosclerosis progression and vascular outcomes.20 33 Possible explanations unravelling the negative correlation advocate that soluble CAMs may compete with leucocyte adhesion to the endothelial molecules or that enhanced ectodomain shedding may contribute to the reduced recruitment of leucocytes to the subendothelial space, thereby promoting cardioprotective effects.34 Our results suggesting a negative relationship between the genetic effects on soluble CAM and circulating cholesterol levels advocate that altered cholesterol metabolism could contribute to the CAD risk associated with the ABO locus; the genetic effects of the same SNPs on LDL-C or TotC show positive correlation with CAD risk. Nevertheless, further studies are warranted to understand the exact mechanisms. Another novel association with sVCAM-1 level was detected in chr12. The lead SNP of this locus is in LD with rs117468318 (r2=1 in NFBC1966) that locates in the 5’UTR of HSP90B1 encoding heat shock protein gp96 and, according to RegulomeDB,28 is likely to affect transcription factor binding, suggesting a possible regulatory mechanism for the detected association. HSP90B1/gp96 is a chaperone that is essential for assembly of 14 of 17 integrin pairs in the haematopoietic system.35 Integrin α4β1 is an important ligand of VCAM-1; if altered transcription of HSP90B1 had a downstream effect on integrin α4β1 level, this could further modify the level of unbound sVCAM-1 in circulation. The third novel locus showing association with sVCAM-1 level was identified in chr17 near ABCA8. The lead SNP rs112001035 is an expression quantitative trait locus (eQTL) for ABCA8 in multiple tissue types.36 If ABCA8 is involved in the regulation of HDL level via interaction with ABCA137 and if plasma HDL levels contribute to VCAM-1 expression,38 then altered expression of the ABCA8 could influence circulating levels of sVCAM-1 by modulating HDL particle concentration. However, this hypothesis is not supported by the fact that the effect of the lead SNP on HDL particle concentration is negligible in a metabolomics GWAS (β=−0.043 SD, p=0.049).39 There is evidence suggesting that ABCA8 may be involved in sphingolipid metabolism,40 and it has been hypothesised that ABCA8 may be involved in the formation of specific membrane domains during ApoA-I lipidation.37 Thus, the association between the ABCA8 and sVCAM-1 level could be related to altered HDL composition possibly contributing to endothelial homeostasis rather than absolute particle concentration. However, more evidence is needed to draw conclusions. We detected a novel effect of rs11220471 in chr11 near ST3GAL4 on sE-selectin levels in the NFBC1966 population. ST3GAL4 encodes a member of the glycosyltransferase 29 family of enzymes involved in protein glycosylation. In mice, St3Gal4 is needed for synthesis of functional selectin ligands.29 The altered levels or structure of selectin ligands due to variation in ST3GAL4 could contribute to the levels of unbound sE-selectin in circulation, providing a biologically rational mechanism for the detected association. In the meta-analyses, we detected a novel large-effect locus for VEGF in chr4 (β=−2.38 SD) near STK32B. Mutations in this locus have been associated previously with coeliac disease,41 CAD42 and Ellis-van Creveld syndrome.43 STK32B may play a role in the hedgehog signalling pathway, which has been implicated in metastasis and angiogenesis in cancer44 and downregulated in coeliac disease.45 The hedgehog signalling has shown to be involved in the regulation of VEGF expression during developmental angiogenesis in avian embryo.46 Thus, previous literature and our results advocate that STK32B may be involved in the regulation of VEGF levels possibly via hedgehog signalling-related mechanism. The other novel findings obtained in meta-analysis include a large-effect locus on TNFα level in chr13 (β=2.13 SD). The locus in 13q14.3 associating with TNFα locates near DLEU1 and DLEU7. This region is recurrently deleted in tumours and haematopoietic malignancies.47 DLEU1 is a part of a transcriptionally coregulated gene cluster that modulates the activity of the nuclear factor kappa B (NF-kB) pathway,48 which is also modulated by TNFα.49 It is largely unknown how the DLEU1 and related DLEU2 regulate NF-kB activity50; our result suggests that TNFα signalling might be involved in this mechanism. At last, we identified a small-effect locus in chr6 harbouring two independent association signals on IL1β and showing suggestive association also on IL4 level. This association signal is in the region coding the human leucocyte antigen proteins, and of Tampere. Protected by copyright. on October 22, 2019 at Tampere University and University Hospitalhttp://jmg.bmj.com/J Med Genet: first published as 10.1136/jmedgenet-2018-105965 on 19 June 2019. Downloaded from 614 Sliz E, etal. J Med Genet 2019;56:607–616. doi:10.1136/jmedgenet-2018-105965 Complex traits further experimental evidence would be needed to identify the exact mechanism how the locus contributes to IL levels. The strengths and limitations of our study should be considered. The sample size of the present study should provide adequate power for detecting genetic associations with circulating markers of systemic inflammation.8 The use of genetically isolated populations, such as inhabitants of Northern Finland, should further enhance the power for locus identification in GWAS settings. We were able to perform meta-analyses only for 10 out of the total of 16 inflammatory phenotypes, and the novel findings are largely based on NFBC1966 population only. Thus, replication of the present findings in other populations would be helpful. In particular, the associations of the novel rare, large-effect variants need to be interpreted with caution until the associations are validated in other populations. The interassay coefficient of variability measures for sE-selectin and VEGF in particular are notably larger than 15%, which is considered to be the limit for acceptable values (online supplementary table S1). However, to our consideration, all the findings identified in the present study locate on genome regions with biologically relevant genes. Furthermore, the replications of the previously reported loci speak for the data adequacy and add confidence to the novel associations. Finally, as we have not included functional experiments in this work, we are limited to previous literature when explaining the potential biological mechanism behind the identified associations. The present results provide novel information on genetic mechanisms influencing levels of inflammatory phenotypes in circulation. The evident role of the ABO locus in the regulation of the soluble CAM levels likely encompasses at least two distinct mechanisms influencing sE-selectin, sICAM-1 and sVCAM-1 levels. Our findings provide evidence that increased soluble CAM concentrations per se may not be a risk factor for cardiovascular outcomes. In particular, genetic variation associating with increased sE-selectin or sICAM-1 levels at the ABO locus seems to contribute to lower cardiovascular risk. Furthermore, genetic effects at the ICAM1 locus providing a direct molecular link to sICAM-1 concentration do not correlate with the genetic effects on CAD risk nor stroke risk. Overall, the present study extends the knowledge about the molecular pathways involved in inflammatory load. Author affiliations 1computational Medicine, Faculty of Medicine, University of Oulu, Oulu, Finland 2center for life course Health research, University of Oulu, Oulu, Finland 3Biocenter Oulu, Oulu, Finland 4Department of internal Medicine, Satakunta central Hospital, Pori, Finland 5research centre of applied and Preventive cardiovascular Medicine, University of turku, turku, Finland 6Department of clinical Physiology and nuclear Medicine, turku University Hospital, turku, Finland 7national institute for Health and Welfare, Helsinki, Finland 8institute for Molecular Medicine Finland, University of Helsinki, Helsinki, Finland 9University of tartu, estonian genome center, tartu, estonia 10Department of clinical chemistry, Fimlab laboratories, and Finnish cardiovascular research center - tampere, Faculty of Medicine and Health technology, tampere University, tampere, Finland 11institute of Biomedicine, University of Oulu, Oulu, Finland 12Department of Psychiatry, University of eastern Finland, and Kuopio University Hospital, Kuopio, Finland 13Medicity research laboratory and institute of Biomedicine, University of turku, turku, Finland 14Unit of general Practice, Oulu University Hospital, Oulu, Finland 15Oulu Deaconess institute/Diapolis Oy research Unit, Oulu, Finland 16northern Finland Birth cohorts, Faculty of Medicine, University of Oulu, Oulu, Finland 17Medical research center (Mrc), University of Oulu, and Oulu University Hospital, Oulu, Finland 18Department of gastroenterology and Metabolism, Poznan University of Medical Sciences, Poznan, Poland 19Department of epidemiology and Biostatistics, Mrc-PHe centre for environment and Health, imperial college london, london, UK 20Unit of Primary care, Oulu University Hospital, Oulu, Finland 21Department of genomics and complex Diseases, School of Public Health, imperial college, london, UK Correction notice this article has been corrected since it was published Online First. Values in table 2 have been amended. Acknowledgements We gratefully acknowledge the contributions of the participants in the northern Finland Birth cohort 1966 study. We also thank all the field workers and laboratory personnel for their efforts. Data on coronary artery disease/myocardial infarction have been contributed by carDiograMplusc4D investigators and have been downloaded from www.carDiOgraMPlUSc4D.Org. Contributors eS, SS and JK conceptualised the study. eS, MK and aa-O performed the statistical analyses. tK, K-HH, MS, KS and SJ performed the quantifications of the inflammatory markers. Or, MP, VS, tl, HV, MS, SJ, JJ, SK-K, MM, K-HH, M-rJ, SS and JK provided funding or other resources to conduct the study. SS and JK supervised the study. eS, SS and JK wrote the original manuscript. all authors contributed to revising the content and approved the final version. Funding this work was supported by the University of Oulu graduate School (eS, MK), the Finnish Foundation for cardiovascular research (VS), Biocenter Oulu (SS), european commission (DynaHealtH – H2020 – 633595; SS), academy of Finland (297338 and 307247; JK) and novo nordisk Foundation (nnF17Oc0026062; JK). nFBc1966 received financial support from University of Oulu (grant no 65354), Oulu University Hospital (grant no 2/97 and 8/97), Ministry of Health and Social affairs (grant no 23/251/97, 160/97 and 190/97), national institute for Health and Welfare, Helsinki (grant no 54121), and regional institute of Occupational Health, Oulu, Finland (grant no 50621 and 54231). the Young Finns Study has been financially supported by the academy of Finland (grants 286284, 134309 (eye), 126925, 121584, 124282, 129378 (Salve), 117787 (gendi) and 41071 (Skidi)); the Social insurance institution of Finland; competitive State research Financing of the expert responsibility area of Kuopio, tampere and turku University Hospitals (grant X51001); Juho Vainio Foundation; Paavo nurmi Foundation; Finnish Foundation for cardiovascular research; Finnish cultural Foundation; Sigrid Juselius Foundation; tampere tuberculosis Foundation; emil aaltonen Foundation; Yrjö Jahnsson Foundation; Signe and ane gyllenberg Foundation; Diabetes research Foundation of Finnish Diabetes association; eU Horizon 2020 (grant 755320 for taXinOMiSiS); european research council (grant 742927 for MUltiePigen project); and tampere University Hospital Supporting Foundation. Competing interests VS has participated in a conference trip sponsored by novo nordisk and received a honorarium from the same source for participating in and advisory board meeting. He also has ongoing research collaboration with Bayer. Provenance and peer review not commissioned; externally peer reviewed. Open access this is an open access article distributed in accordance with the creative commons attribution non commercial (cc BY-nc 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See: http:// creativecommons. org/ licenses/ bync/ 4. 0/. reFerenCes 1 Donath MY. targeting inflammation in the treatment of type 2 diabetes: time to start. Nat Rev Drug Discov 2014;13:465–76. 2 Heneka Mt, Kummer MP, latz e. innate immune activation in neurodegenerative disease. Nat Rev Immunol 2014;14:463–77. 3 coussens lM, Werb Z. inflammation and cancer. Nature 2002;420:860–7. 4 goldberg el, Dixit VD. Drivers of age-related inflammation and strategies for healthspan extension. Immunol Rev 2015;265:63–74. 5 World Health Organization. World report on ageing and health 2015. 6 de craen aJM, Posthuma D, remarque eJ, van den Biggelaar aHJ, Westendorp rgJ, Boomsma Di. Heritability estimates of innate immunity: an extended twin study. Genes Immun 2005;6:167–70. 7 ahola-Olli aV, Würtz P, Havulinna aS, aalto K, Pitkänen n, lehtimäki t, Kähönen M, lyytikäinen l-P, raitoharju e, Seppälä i, Sarin a-P, ripatti S, Palotie a, Perola M, Viikari JS, Jalkanen S, Maksimow M, Salomaa V, Salmi M, Kettunen J, raitakari Ot. genomewide association study identifies 27 loci influencing concentrations of circulating cytokines and growth factors. Am J Hum Genet 2017;100:40–50. 8 Sun BB, Maranville Jc, Peters Je, Stacey D, Staley Jr, Blackshaw J, Burgess S, Jiang t, Paige e, Surendran P, Oliver-Williams c, Kamat Ma, Prins BP, Wilcox SK, Zimmerman eS, chi a, Bansal n, Spain Sl, Wood aM, Morrell nW, Bradley Jr, Janjic n, roberts DJ, Ouwehand WH, todd Ja, Soranzo n, Suhre K, Paul DS, Fox cS, Plenge rM, Danesh J, runz H, Butterworth aS. genomic atlas of the human plasma proteome. Nature 2018;558:73–9. of Tampere. Protected by copyright. on October 22, 2019 at Tampere University and University Hospitalhttp://jmg.bmj.com/J Med Genet: first published as 10.1136/jmedgenet-2018-105965 on 19 June 2019. Downloaded from 615 Sliz E, etal. J Med Genet 2019;56:607–616. doi:10.1136/jmedgenet-2018-105965 Complex traits 9 ligthart S, Vaez a, Võsa U, Stathopoulou Mg, de Vries PS, Prins BP, Van der Most PJ, tanaka t, naderi e, rose lM, Wu Y, Karlsson r, Barbalic M, lin H, Pool r, Zhu g, Macé a, Sidore c, trompet S, Mangino M, Sabater-lleal M, Kemp JP, abbasi a, Kacprowski t, Verweij n, Smith aV, Huang t, Marzi c, Feitosa MF, lohman KK, Kleber Me, Milaneschi Y, Mueller c, Huq M, Vlachopoulou e, lyytikäinen l-P, Oldmeadow c, Deelen J, Perola M, Zhao JH, Feenstra B, amini M, lahti J, Schraut Ke, Fornage M, Suktitipat B, chen W-M, li X, nutile t, Malerba g, luan Jian’an, Bak t, Schork n, Del greco M F, thiering e, Mahajan a, Marioni re, Mihailov e, eriksson J, Ozel aB, Zhang W, nethander M, cheng Y-c, aslibekyan S, ang W, gandin i, Yengo l, Portas l, Kooperberg c, Hofer e, rajan KB, Schurmann c, den Hollander W, ahluwalia tS, Zhao J, Draisma HHM, Ford i, timpson n, teumer a, Huang H, Wahl S, liu Y, Huang J, Uh H-W, geller F, Joshi PK, Yanek lr, trabetti e, lehne B, Vozzi D, Verbanck M, Biino g, Saba Y, Meulenbelt i, O’connell Jr, laakso M, giulianini F, Magnusson PKe, Ballantyne cM, Hottenga JJ, Montgomery gW, rivadineira F, rueedi r, Steri M, Herzig K-H, Stott DJ, Menni c, Frånberg M, St Pourcain B, Felix SB, Pers tH, Bakker SJl, Kraft P, Peters a, Vaidya D, Delgado g, Smit JH, großmann V, Sinisalo J, Seppälä i, Williams Sr, Holliday eg, Moed M, langenberg c, räikkönen K, Ding J, campbell H, Sale MM, chen Y-Di, James al, ruggiero D, Soranzo n, Hartman ca, Smith en, Berenson gS, Fuchsberger c, Hernandez D, tiesler cMt, giedraitis V, liewald D, Fischer K, Mellström D, larsson a, Wang Y, Scott Wr, lorentzon M, Beilby J, ryan Ka, Pennell ce, Vuckovic D, Balkau B, concas MP, Schmidt r, Mendes de leon cF, Bottinger eP, Kloppenburg M, Paternoster l, Boehnke M, Musk aW, Willemsen g, evans DM, Madden PaF, Kähönen M, Kutalik Z, Zoledziewska M, Karhunen V, Kritchevsky SB, Sattar n, lachance g, clarke r, Harris tB, raitakari Ot, attia Jr, van Heemst D, Kajantie e, Sorice r, gambaro g, Scott ra, Hicks aa, Ferrucci l, Standl M, lindgren cM, Starr JM, Karlsson M, lind l, li JZ, chambers Jc, Mori ta, de geus eJcn, Heath ac, Martin ng, auvinen J, Buckley BM, de craen aJM, Waldenberger M, Strauch K, Meitinger t, Scott rJ, Mcevoy M, Beekman M, Bombieri c, ridker PM, Mohlke Kl, Pedersen nl, Morrison ac, Boomsma Di, Whitfield JB, Strachan DP, Hofman a, Vollenweider P, cucca F, Jarvelin M-r, Jukema JW, Spector tD, Hamsten a, Zeller t, Uitterlinden ag, nauck M, gudnason V, Qi l, grallert H, Borecki iB, rotter Ji, März W, Wild PS, lokki M-l, Boyle M, Salomaa V, Melbye M, eriksson Jg, Wilson JF, Penninx BWJH, Becker DM, Worrall BB, gibson g, Krauss rM, ciullo M, Zaza g, Wareham nJ, Oldehinkel aJ, Palmer lJ, Murray SS, Pramstaller PP, Bandinelli S, Heinrich J, ingelsson e, Deary iJ, Mägi r, Vandenput l, van der Harst P, Desch Kc, Kooner JS, Ohlsson c, Hayward c, lehtimäki t, Shuldiner ar, arnett DK, Beilin lJ, robino a, Froguel P, Pirastu M, Jess t, Koenig W, loos rJF, evans Da, Schmidt H, Smith gD, Slagboom Pe, eiriksdottir g, Morris aP, Psaty BM, tracy rP, nolte iM, Boerwinkle e, Visvikis-Siest S, reiner aP, gross M, Bis Jc, Franke l, Franco OH, Benjamin eJ, chasman Di, Dupuis J, Snieder H, Dehghan a, alizadeh BZ, lifelines cohort Study, cHarge inflammation Working group. genome analyses of >200,000 individuals identify 58 loci for chronic inflammation and Highlight Pathways that link inflammation and complex Disorders. Am J Hum Genet 2018;103:691–706. 10 Järvelin Mr, Sovio U, King V, lauren l, Xu B, Mccarthy Mi, Hartikainen al, laitinen J, Zitting P, rantakallio P, elliott P. early life factors and blood pressure at age 31 years in the 1966 northern Finland birth cohort. Hypertension 2004;44:838–46. 11 Saukkonen t, Mutt SJ, Jokelainen J, Saukkonen aM, raza gS, Karhu t, Härkönen P, eckel J, Herzig KH, rajala U, Keinänen-Kiukaanniemi S. adipokines and inflammatory markers in elderly subjects with high risk of type 2 diabetes and cardiovascular disease. Sci Rep 2018;8:1–8. 12 raitakari Ot, Juonala M, rönnemaa t, Keltikangas-Järvinen l, räsänen l, Pietikäinen M, Hutri-Kähönen n, taittonen l, Jokinen e, Marniemi J, Jula a, telama r, Kähönen M, lehtimäki t, Åkerblom HK, Viikari JSa. cohort profile: the cardiovascular risk in young finns study. Int J Epidemiol 2008;37:1220–6. 13 Borodulin K, tolonen H, Jousilahti P, Jula a, Juolevi a, Koskinen S, Kuulasmaa K, laatikainen t, Männistö S, Peltonen M, Perola M, Puska P, Salomaa V, Sundvall J, Virtanen SM, Vartiainen e. cohort profile: the national FinriSK study. Int J Epidemiol 2017;47:696–696i. 14 Santalahti K, Maksimow M, airola a, Pahikkala t, Hutri-Kähönen n, Jalkanen S, raitakari Ot, Salmi M. circulating cytokines predict the development of insulin resistance in a prospective Finnish population cohort. J Clin Endocrinol Metab 2016;101:3361–9. 15 Marchini J, Howie B, Myers S, McVean g, Donnelly P. a new multipoint method for genome-wide association studies by imputation of genotypes. Nat Genet 2007;39:906–13. 16 Willer cJ, li Y, abecasis gr. Metal: fast and efficient meta-analysis of genomewide association scans. Bioinformatics 2010;26:2190–1. 17 Yang J, Ferreira t, Morris aP, Medland Se, Madden PaF, Heath ac, Martin ng, Montgomery gW, Weedon Mn, loos rJ, Frayling tM, Mccarthy Mi, Hirschhorn Jn, goddard Me, Visscher PM. conditional and joint multiple-SnP analysis of gWaS summary statistics identifies additional variants influencing complex traits. Nat Genet 2012;44:369–75. 18 Paré g, chasman Di, Kellogg M, rYl Z, rifai n, Badola S, Miletich JP, ridker PM. novel association of aBO histo-blood group antigen with soluble icaM-1: results of a genome-wide association study of 6,578 women. PLoS Genet 2008;4:e1000118. 19 ridker H, Buring r. c-reactive protein and other markers of inflammation in the prediction of cardiovascular disease in women. Ital Heart J Suppl 2000;1:1066–7. 20 Hwang SJ, Ballantyne cM, Sharrett ar, Smith lc, Davis ce, gotto aM, Boerwinkle e. circulating adhesion molecules VcaM-1, icaM-1, and e-selectin in carotid atherosclerosis and incident coronary heart disease cases: the atherosclerosis risk in communities (aric) study. Circulation 1997;96:4219–25. 21 Pickrell JK, Berisa t, liu JZ, Segurel l, tung JY, Hinds D. Detection and interpretation of shared genetic influences on 42 human traits 2017;48:709–17. 22 nikpay M, goel a, Won H-H, Hall lM, Willenborg c, Kanoni S, Saleheen D, Kyriakou t, nelson cP, Hopewell Jc, Webb tr, Zeng l, Dehghan a, alver M, armasu SM, auro K, Bjonnes a, chasman Di, chen S, Ford i, Franceschini n, gieger c, grace c, gustafsson S, Huang J, Hwang S-J, Kim YK, Kleber Me, lau KW, lu X, lu Y, lyytikäinen l-P, Mihailov e, Morrison ac, Pervjakova n, Qu l, rose lM, Salfati e, Saxena r, Scholz M, Smith aV, tikkanen e, Uitterlinden a, Yang X, Zhang W, Zhao W, de andrade M, de Vries PS, van Zuydam nr, anand SS, Bertram l, Beutner F, Dedoussis g, Frossard P, gauguier D, goodall aH, gottesman O, Haber M, Han B-g, Huang J, Jalilzadeh S, Kessler t, König ir, lannfelt l, lieb W, lind l, lindgren cM, lokki M-l, Magnusson PK, Mallick nH, Mehra n, Meitinger t, Memon F-U-r, Morris aP, nieminen MS, Pedersen nl, Peters a, rallidis lS, rasheed a, Samuel M, Shah SH, Sinisalo J, Stirrups Ke, trompet S, Wang l, Zaman KS, ardissino D, Boerwinkle e, Borecki iB, Bottinger eP, Buring Je, chambers Jc, collins r, cupples la, Danesh J, Demuth i, elosua r, epstein Se, esko t, Feitosa MF, Franco OH, Franzosi Mg, granger cB, gu D, gudnason V, Hall aS, Hamsten a, Harris tB, Hazen Sl, Hengstenberg c, Hofman a, ingelsson e, iribarren c, Jukema JW, Karhunen PJ, Kim B-J, Kooner JS, Kullo iJ, lehtimäki t, loos rJF, Melander O, Metspalu a, März W, Palmer cn, Perola M, Quertermous t, rader DJ, ridker PM, ripatti S, roberts r, Salomaa V, Sanghera DK, Schwartz SM, Seedorf U, Stewart aF, Stott DJ, thiery J, Zalloua Pa, O’Donnell cJ, reilly MP, assimes tl, thompson Jr, erdmann J, clarke r, Watkins H, Kathiresan S, McPherson r, Deloukas P, Schunkert H, Samani nJ. a comprehensive 1,000 genomes-based genome-wide association meta-analysis of coronary artery disease. Nat Genet 2015;47:1121–30. 23 Malik r, chauhan g, traylor M, Sargurupremraj M, Okada Y, Mishra a, rutten-Jacobs l, giese aK, Van Der laan SW, gretarsdottir S, anderson cD, chong M, adams HHH, ago t, almgren P, amouyel P, ay H, Bartz tM, Benavente Or, Bevan S, Boncoraglio gB, Brown rD, Butterworth aS, carrera c, carty cl, chasman Di, chen WM, cole JW, correa a, cotlarciuc i, cruchaga c, Danesh J, De Bakker PiW, Destefano al, Den Hoed M, Duan Q, engelter St, Falcone gJ, gottesman rF, grewal rP, gudnason V, gustafsson S, Haessler J, Harris tB, Hassan a, Havulinna aS, Heckbert Sr, Holliday eg, Howard g, Hsu Fc, Hyacinth Hi, ikram Ma, ingelsson e, irvin Mr, Jian X, Jiménezconde J, Johnson Ja, Jukema JW, Kanai M, Keene Kl, Kissela BM, Kleindorfer DO, Kooperberg c, Kubo M, lange la, langefeld cD, langenberg c, launer lJ, lee JM, lemmens r, leys D, lewis cM, lin WY, lindgren ag, lorentzen e, Magnusson PK, Maguire J, Manichaikul a, Mcardle PF, Meschia JF, Mitchell BD, Mosley tH, nalls Ma, ninomiya t, O’Donnell MJ, Psaty BM, Pulit Sl, rannikmäe K, reiner aP, rexrode KM, rice K, rich SS, ridker PM, rost nS, rothwell PM, rotter Ji, rundek t, Sacco rl, Sakaue S, Sale MM, Salomaa V, Sapkota Br, Schmidt r, Schmidt cO, Schminke U, Sharma P, Slowik a, Sudlow clM, tanislav c, tatlisumak t, taylor KD, thijs VnS, thorleifsson g, thorsteinsdottir U, tiedt S, trompet S, tzourio c, Van Duijn cM, Walters M, Wareham nJ, Wassertheil-Smoller S, Wilson Jg, Wiggins Kl, Yang Q, Yusuf S, Bis Jc, Pastinen t, ruusalepp a, Schadt ee, Koplev S, Björkegren JlM, codoni V, civelek M, Smith nl, trégouët Da, christophersen ie, roselli c, lubitz Sa, ellinor Pt, tai eS, Kooner JS, Kato n, He J, Van Der Harst P, elliott P, chambers Jc, takeuchi F, Johnson aD, Sanghera DK, Melander O, Jern c, Strbian D, Fernandez-cadenas i, longstreth Wt, rolfs a, Hata J, Woo D, rosand J, Pare g, Hopewell Jc, Saleheen D, Stefansson K, Worrall BB, Kittner SJ, Seshadri S, Fornage M, Markus HS, Howson JMM, Kamatani Y, Debette S, Dichgans M, amin n, aparicio HS, arnett DK, attia J, Beiser aS, Berr c, Buring Je, Bustamante M, caso V, cheng Yc, choi SH, chowhan a, cullell n, Dartigues JF, Delavaran H, Delgado P, Dörr M, engström g, Ford i, gurpreet WS, Hamsten a, Heitsch l, Hozawa a, ibanez l, ilinca a, ingelsson M, iwasaki M, Jackson rD, Jood K, Jousilahti P, Kaffashian S, Kalra l, Kamouchi M, Kitazono t, Kjartansson O, Kloss M, Koudstaal PJ, Krupinski J, labovitz Dl, laurie cc, levi cr, li l, lind l, lindgren cM, lioutas V, liu YM, lopez Ol, Makoto H, Martinez-Majander n, Matsuda K, Minegishi n, Montaner J, Morris aP, Muiño e, Müller-nurasyid M, norrving B, Ogishima S, Parati ea, Peddareddygari lr, Pedersen nl, Pera J, Perola M, Pezzini a, Pileggi S, rabionet r, riba-llena i, ribasés M, romero Jr, roquer J, rudd ag, Sarin aP, Sarju r, Sarnowski c, Sasaki M, Satizabal cl, Satoh M, Sattar n, Sawada n, Sibolt g, Á S, Smith a, Sobue K, Soriano-tárraga c, Stanne t, Stine Oc, Stott DJ, Strauch K, takai t, tanaka H, tanno K, teumer a, tomppo l, torres-aguila nP, touze e, tsugane S, Uitterlinden ag, Valdimarsson eM, Van Der lee SJ, Völzke H, Wakai K, Weir D, Williams Sr, Wolfe cDa, Wong Q, Xu H, Yamaji t. Multiancestry genome-wide association study of 520,000 subjects identifies 32 loci associated with stroke and stroke subtypes. Nat Genet 2018;50:524–37. 24 Willer cJ, Schmidt eM, Sengupta S, Peloso gM, gustafsson S, Kanoni S, ganna a, chen J, Buchkovich Ml, Mora S, Beckmann JS, Bragg-gresham Jl, chang H-Y, Demirkan a, Den Hertog HM, Do r, Donnelly la, ehret gB, esko t, Feitosa MF, Ferreira t, Fischer K, Fontanillas P, Fraser rM, Freitag DF, gurdasani D, Heikkilä K, Hyppönen e, isaacs a, Jackson aU, Johansson a, Johnson t, Kaakinen M, Kettunen J, Kleber Me, li X, luan J, lyytikäinen l-P, Magnusson PKe, Mangino M, Mihailov e, Montasser Me, Müller-nurasyid M, nolte iM, O’connell Jr, Palmer cD, Perola M, Petersen a-K, Sanna S, Saxena r, Shah S, Shungin D, Sidore c, Song c, Strawbridge rJ, Surakka i, tanaka t, teslovich tM, thorleifsson g, Van den Herik eg, Voight BF, Volcik Ka, Waite ll, Wong a, Wu Y, Zhang W, absher D, asiki g, Barroso i, Been lF, Bolton Jl, Bonnycastle ll, Brambilla P, Burnett MS, cesana g, Dimitriou M, Doney aSF, Döring of Tampere. Protected by copyright. on October 22, 2019 at Tampere University and University Hospitalhttp://jmg.bmj.com/J Med Genet: first published as 10.1136/jmedgenet-2018-105965 on 19 June 2019. Downloaded from