A gene-centric association scan for Coagulation Factor VII levels in European and African Americans: the Candidate Gene Association Resource (CARe) Consortium

doi:10.1093/hmg/ddr264

Hum Mol Genet. 2011 September 1; 20(17): 3525–3534.

Published online 2011 June 15. doi: 10.1093/hmg/ddr264

PMCID: PMC3153310

A gene-centric association scan for Coagulation Factor VII levels in European and African Americans: the Candidate Gene Association Resource (CARe) Consortium

Kira C. Taylor,¹^* Leslie A. Lange,² Delilah Zabaneh,³ Ethan Lange,² Brendan J. Keating,⁵ Weihong Tang,⁶ Nicholas L. Smith,⁷ Joseph A. Delaney,⁸ Meena Kumari,⁹ Aroon Hingorani,¹⁰ Kari E. North,¹ Mika Kivimaki,⁹ Russell P. Tracy,¹¹ Christopher J. O'Donnell,¹² Aaron R. Folsom,⁶ David Green,¹³ Steve E. Humphries,⁴ and Alexander P. Reiner⁷

¹Department of Epidemiology, University of North Carolina at Chapel Hill Gillings School of Global Public Health, Chapel Hill, NC 27514, USA,

²Department of Genetics, University of North Carolina at Chapel Hill School of Medicine, Chapel Hill, NC 27514, USA,

³University College London Genetics Institute, Department of Genetics Environment and Evolution, and

⁴Centre for Cardiovascular Genetics, BHF Laboratories, Institute of Cardiovascular Science, University College London, London WC1E 6JJ, UK,

⁵Center for Applied Genomics, Abramson Research Center, The Children's Hospital of Philadelphia, Philadelphia, PA 19104, USA,

⁶Division of Epidemiology and Community Health, University of Minnesota School of Public Health, Minneapolis, MN 55454, USA,

⁷Department of Epidemiology, University of Washington School of Public Health, Seattle, WA 98195, USA,

⁸Department of Pharmaceutical Outcomes and Pharmacy, University of Florida College of Pharmacy, Gainesville, FL 32610, USA,

⁹Department of Epidemiology and Population Health, London School of Hygiene and Tropical Medicine, London WC1E 7HT, UK,

¹⁰Genetic Epidemiology Group, Department of Epidemiology and Public Health, University College London, London WC1E 6BT, UK,

¹¹Department of Pathology and Laboratory Medicine, University of Vermont College of Medicine, Burlington, VT 05405, USA,

¹²National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, MD 20824, USA, and

¹³Division of Hematology/Oncology, Northwestern University Feinberg School of Medicine, Chicago, IL 60611, USA

^*To whom correspondence should be addressed at: Department of Epidemiology, UNC Gillings School of Global Public Health, 137 E Franklin St., Suite 306, Chapel Hill, NC 27514, USA. Tel: +Phone: 1 4048081718; Fax: +1 9199669800; Email: kirat/at/email.unc.edu

Received February 3, 2011; Revised May 13, 2011; Accepted June 6, 2011.

Abstract

Polymorphisms in several distinct genomic regions, including the F7 gene, were recently associated with factor VII (FVII) levels in European Americans (EAs). The genetic determinants of FVII in African Americans (AAs) are unknown. We used a 50 000 single nucleotide polymorphism (SNP) gene-centric array having dense coverage of over 2 000 candidate genes for cardiovascular disease (CVD) pathways in a community-based sample of 16 324 EA and 3898 AA participants from the Candidate Gene Association Resource (CARe) consortium. Our aim was the discovery of new genomic loci and more detailed characterization of existing loci associated with FVII levels. In EAs, we identified three new loci associated with FVII, of which APOA5 on chromosome 11q23 and HNF4A on chromosome 20q12–13 were replicated in a sample of 4289 participants from the Whitehall II study. We confirmed four previously reported FVII-associated loci (GCKR, MS4A6A, F7 and PROCR) in CARe EA samples. In AAs, the F7 and PROCR regions were significantly associated with FVII. Several of the FVII-associated regions are known to be associated with lipids and other cardiovascular-related traits. At the F7 locus, there was evidence of at least five independently associated SNPs in EAs and three independent signals in AAs. Though the variance in FVII explained by the existing loci is substantial (20% in EA and 10% in AA), larger sample sizes and investigation of lower frequency variants may be required to identify additional FVII-associated loci in EAs and AAs and further clarify the relationship between FVII and other CVD risk factors.

INTRODUCTION

Coagulation Factor VII (FVII) is a vitamin K-dependent serine protease that is activated in response to vascular or tissue injury. Activated FVII (FVIIa) interacts with tissue factor to initiate coagulation (1). The association of FVII with cardiovascular disease (CVD) events is uncertain, with some positive and negative findings (2–6). Commonly, FVII had an association in unadjusted analyses, but this was often removed by adjusting for known CVD risk factors, especially lipids. Polymorphisms in the FVII gene (F7) have been associated with increased risk of coronary heart disease and stroke in some studies but not all (7–10).

FVII levels are highly heritable (estimates range from 53 to 63%) and polymorphisms in F7 account for approximately a third of the variance of FVII levels in European Americans (EAs) (11–13). A recent genome-wide association study (GWAS) in EAs confirmed that polymorphisms in the F7 gene region on chromosome 13q34 are strongly associated with FVII levels and also identified novel associations in four additional genomic regions: GCKR, ADH4, MS4A6A and PROCR (14). Whether any of these or other genetic loci are associated with FVII levels in African Americans (AAs) is unknown.

The purpose of this study was to (i) perform more detailed characterization of FVII-associated loci and fine-mapping of previously discovered associations in EAs; and (ii) identify genes associated with FVII levels in AAs. We analyzed samples from a consortium of community-based prospective cohort studies [the Candidate Gene Association Resource (CARe)] with EA and AA participants genotyped using the custom 50 K single nucleotide polymorphism (SNP) cardiovascular gene-centric ITMAT Broad-CARe (IBC) array. The IBC array has dense coverage of over 2000 CVD candidate genes and includes SNPs selected to ensure that common patterns of genetic variation in populations of European and African descent are well-captured.

RESULTS

Participant characteristics

The characteristics of our EA and AA study populations are reported in Table 1. Prior to exclusion of individuals due to low genotyping call rates, there were a total of 16 352 EAs and 3948 AAs. The mean age at FVII measurement ranged from 30 [Coronary Artery Risk in Young Adults (CARDIA)] to 73 [Cardiovascular Health Study (CHS)] years. Mean FVII levels were higher in the studies which sampled older individuals. The prevalence of cardiovascular risk factors tended to be higher among AAs than EAs. FVII levels did not consistently differ by race, with higher levels found among AAs compared with EAs in CHS and lower levels found among AAs compared with EAs in CARDIA. After filtering individuals with less than 90% genotype call rate, our final meta-analysis sample sizes were 16 324 EA and 3898 AA subjects.

Table 1.

Characteristics of the 16352 EA and 3948 AA study participants^a

FVII association results in EAs

In EAs, 7 discrete regions spanning 12 candidate genes on the IBC array contained SNPs significantly associated with FVII levels (Table 2). Four of the FVII association signals in EAs correspond to loci previously reported in a GWAS of three of the same cohorts (GCKR on chromosome 2p23, MS4A2-MS4A6A on chromosome 11q12, F7–F10 on chromosome 13q34 and PROCR on chromosome 20q11.2) (14). The addition of the CARDIA cohort to the three cohorts used in Smith et al. [Atherosclerosis Risk in Communities (ARIC), CHS and Framingham Heart Study (FHS)] strengthened these associations. A fifth FVII-associated locus reported by Smith et al., which spans two alcohol dehydrogenase (ADH) genes, ADH4 and ADH5, was not represented on the IBC array.

Table 2.

Loci associated with FVII in EAs

Three of the genomic regions significantly associated with FVII in EAs have not previously been reported, including eight SNPs in or near the ZNF259/BUD13/APOA5 region on chromosome 11q23, one SNP upstream of MLXIPL on chromosome 7q11 and a non-synonymous SNP in HNF4A on chromosome 20q12. To further investigate the three newly associated FVII loci, we examined the results in 4289 participants from the Whitehall II study (Table 3). For the top SNP in the ZNF259/BUD13/APOA5 region on chromosome 11q23 (rs2266788) and the HNF4A non-synonymous SNP (rs1800961), the magnitudes and directions of association were similar in the Whitehall II study, replicating the CARe results. Upon meta-analysis of the CARe and the Whitehall II results, both SNPs reached or closely approached the conventional level of ‘genome-wide significance' (< 5×10⁻⁸). The chromosome 7q11 SNP upstream of MLXIPL (rs7777102) was not available in the Whitehall II data set, and therefore could not be validated.

Table 3.

Replication of novel FVII-associated loci in EAs in the Whitehall II Study

FVII association results in AAs

In AA participants, three distinct regions were associated with FVII (Table 4). Two of these regions were concordantly associated with FVII in EAs (F7-F10, PROCR), while one region appeared to be specific to AAs (two SNPs at 10p15 near PRKCQ), although the direction of the association was similar to what was observed in EAs. Of the five remaining loci associated with FVII in EAs, GCKR, MS4A2 and HNF4A showed a similar direction of association in AAs (Table 5). In AAs, effect estimates were smaller and P-values were larger for the top SNP associations noted in EAs. There were substantial observed allele frequency differences between the two populations for several of the top SNPs that would result in differences in statistical power to detect an effect when present. Notably, the PRKCQ SNP (rs602419) had a minor allele frequency of only 0.02 in EAs, compared with 0.32 in AAs. In addition, observed linkage disequilibrium (LD) patterns differed between EAs and AAs in the region of F7 and F10 genes, as shown in Supplementary Material, Figures S1 and S2, respectively, with the extent of LD considerably lower among AAs.

Table 4.

Loci associated with FVII in AAs

Table 5.

Comparison of top SNPs from Tables 2 and 3 in EAs and AAs

To further investigate the F7, PROCR and newly associated PRKCQ region in AA, we assessed the association of the top SNPs in 918 AA participants from the Women's Health Initiative (WHI) (Table 6). The magnitude and direction of associations for the top F7 and PROCR SNPs were similar between WHI and CARe, replicating the CARe AA results. While the novel PRKCQ SNP was not significantly associated with FVII levels in the WHI sample, the direction of association was the same.

Table 6.

Replication of array-wide significant regions identified in AAs in WHI-SHARe AAs (N= 918)

Fine-mapping of FVII-associated loci in EAs and AAs

Regional LD association plots for array-wide significant loci in EAs and AAs are shown in Supplementary Material, Figures S3–12. Based on assessment of regional LD association patterns and results from the conditional regression analyses, there was little evidence supporting multiple independent association signals within the PROCR, PRKCQ, HNF4A, MS4A2, APOA5, MLXIPL regions, and moderate evidence for the GCKR region. In contrast, there was evidence for more than one independent association at an array-wide significant level in the F7 gene region in both EAs and AAs (Supplementary Material, Table S1).

The regional plots for F7 (Supplementary Material, Figures S3 and S4) demonstrate multiple significant loci in the F7–F10 region for EA and AA. Step-wise regression analyses revealed five independent signals in EAs (rs561241, rs555212, rs3093265, rs3211727, rs493833) and three independent signals in AAs (rs561241, rs493833, rs3093230) (Table 7). All signals identified were in the F7 gene, except for one in the F10 gene in EAs. In EA, the minor allele of SNPs rs555212 was associated with higher FVII levels, whereas the minor alleles of SNPs rs561241, rs3093265, rs3211727 and rs493833 were associated with lower FVII levels. Together, these SNPs explained 18% of the variation of normalized FVII levels in EAs and 9% of the variation in AAs.

Table 7.

Fine mapping of the F7/F10 region in EAs and AAs

Observed LD patterns differed between EAs and AAs in the region of the F7 and F10 genes, as shown in Supplementary Material, Figures S1 and S2, with the extent of LD considerably lower among AAs. The lower extent of LD between SNPs in the F7, along with the results of conditional association analysis, could be used to fine-map the FVII association signals in this region. For example, the most significant SNP (rs561241) in EAs, located within the F7 promoter region, is in strong LD (r²> 0.9) with several other SNPs that span the entire F7 gene region, including rs1755685, rs6039, rs7981123 and rs493833. In AAs, the same SNP rs561241 showed LD (r²> 0.6) with only one other SNP, rs1755685, also in the promoter region. These findings suggest that the main signal associated with lower FVII levels is localized to the F7 promoter region. Similarly, the promoter SNP associated with higher FVII levels in EA (rs555212) is in strong LD with at least two other SNPs (rs3093229 and rs3093230) located in the promoter region (r²> 0.98). In AA, there is less LD between rs555212, rs3093229 and rs3093230 (pairwise r² ranging between 0.4 and 0.5). Of these three SNPs, rs3093230 is most strongly associated with higher FVII in AA.

Haplotype-based analysis confirmed the individual F7 SNP association results (Supplementary Material, Table S2). There were three prominent F7 haplotypes in EAs. The haplotype tagged by rs555212 was associated with higher FVII levels, while the haplotype tagged by rs561241 was associated with lower FVII levels. In AA, there was greater haplotype diversity, with five haplotypes having frequencies of at least 5%. There are two haplotypes tagged by the minor allele of rs555212; only the haplotype that additionally contains the rs3093230 minor allele is associated with higher FVII levels in AA.

Total variance explained by known and novel loci

When SNPs representing the five independent signals for the F7/F10 region, the two potentially independent signals from the GKCR region and the other loci identified in Tables 2 and and33 were all included simultaneously in a linear regression model for the ARIC study, they explained a total of 20.2% of the variance in age- and sex-adjusted normalized FVII levels in EAs, and 10.2% in AAs (adjusted R² = 0.202 and 0.102, respectively). Most of this variance is due to variation in the F7 gene, with the five independent regions in the F7 gene explaining 18% of the variance, and the remaining loci from Table 2 explaining the remaining 2% in ARIC EAs.

Comparison of genetic association results for FVII activity versus FVII antigen assays

For each genomic region, associations of the top SNPs with FVII antigen levels closely mirrored those with FVII activity levels. Effect estimates and standard errors were similar, validating previous studies that have meta-analyzed activity and antigen levels (Supplementary Material, Table S3).

DISCUSSION

By using a custom 50 K SNP genotyping array with dense, multi-population tag SNP coverage of >2 000 genes in CVD-related pathways including coagulation, we further characterized regions of the genome associated with FVII levels in EAs and AAs. Our findings can be summarized as follows: (i) the chromosome 13q34 F7–F10 gene region contained multiple variants associated with higher and lower FVII in both EAs and AAs. (ii) The lower extent of LD among AA enabled finer localization of the strongest F7 signals. (iii) We identified two novel genomic regions significantly associated with FVII in EAs, the ZNF259–BUD13–APOA5 gene region on chromosome 11q23 and chromosome 20q12 HNF4A gene (4). Of other candidate regions associated with FVII in EAs, we found evidence for replication of the chromosome 20 PROCR SNP rs867186 (Ser219Gly) and GCKR associations with higher FVII levels in AAs; on the other hand, FVII association with polymorphisms within the MS4A2, APOA5 and HNF4A gene regions were not evident in AAs.

Polymorphisms in the F7 gene could be grouped in into three major haplotypes in European descent populations (5). Compared with the most common or ‘wild-type’ haplotype, the haplotype tagged by the minor allele of rs555212 (and rs3093230) was associated with higher FVII levels while the haplotype tagged by the minor allele at rs561241 was associated with lower levels. These results confirm previous findings from smaller European population studies that SNPs tagging these haplotypes constitute the major determinants of FVII activity in EAs (5,15). Because of the large sample size of CARe and the dense SNP coverage of the IBC genotyping array, we were able to identify at least three additional variants within the F7–F10 genes independently associated FVII levels. Two of these are located within intron 5 of the F7 gene (rs493833 and rs3093265) and one within intron 1 of the F10 gene (rs3093268).The strong LD among the promoter variants in Europeans has previously hampered detailed characterization of the functional variants responsible for reduced FVII levels. Functional analysis of the F7 promoter suggested that the −323ins10 and −122C (rs561241) variants are strongly associated with decreased promoter activity. It has remained uncertain whether the −402A variant (rs510317) or other upstream variants are responsible for the increased F7 promoter activity (16). In our analysis, AA showed greater F7 haplotype diversity, which suggested further localization of the higher-expression F7 promoter region to rs3093230 (or an untyped variant in LD).

Characterization of common cis-acting F7 gene variants may have relevance for athero-thrombotic disease risk. The F7 polymorphisms associated with higher FVII levels have been associated with increased risk of stroke and myocardial infarction (MI) in young women (5,17,18).The relationship of haplotype C FVII-lowering polymorphisms to risk of clinical CVD has been equivocal, with published meta-analyses of MI (10) and stroke (19) showing no significant association. Since our results suggest independent associations of several additional F7 polymorphisms on FVII level, evaluation of these variants in the context of clinical CVD risk, or extent of subclinical atherosclerosis in young adults (20) both in EA and AA populations, may be warranted.

Polymorphisms in LD with the PROCR Ser219Gly variant (rs867186) have been strongly associated with higher soluble endothelial protein C receptor (EPCR) levels, as well as with higher levels of clotting FVII, FVIIa, protein C antigen and several downstream markers of activated coagulation in the extrinsic pathway (21–22) in whites. Binding of FVII/FVIIa to soluble and/or membrane-associated EPCR and subsequent endocytosis of the receptor-ligand complexes may be important for regulating the circulating concentrations of FVII and protein C (23). These findings have potential implications for the role of the PROCR Ser219Gly dimorphism in risk of clinical thrombotic disease (24,25).

FVII levels are correlated with other CVD risk factors, such as lipid levels, body mass index and insulin resistance (6,26). The causal nature and mechanisms of these associations remain to be elucidated. There is a fairly well-described connection between high-fat diet, triglycerides and FVII activity. FVII and other vitamin K-dependent coagulation proteins bind to triglyceride-rich lipid particles (27). Moreover, expression of the F7 gene can be modulated by glucose and insulin levels due to a specific promoter element (28). The association of FVII levels with polymorphisms of the GCKR, APOA5, MLXIPL and HNF4A genes suggests an additional level of shared genetic regulation between FVII synthesis, glucose metabolism and low-grade inflammation. The APOA5 SNP rs6589566 that was discovered in our study was also associated with low density lipoprotein levels in a recent genome-wide association scan (29). Other SNPs in LD with rs6589566 (e.g. rs964184, rs1558861, rs4938303, rs12280753) have been associated with high density lipoprotein (HDL-C) levels and triglyceride levels (30–34).

While the same SNPs were not identified in this study, variants in or near PRKCQ have been found to be associated at a genome-wide significance level with type I diabetes (35–36) and rheumatoid arthritis (37–38). One non-synonymous SNP in MXLIPL was also significantly associated with plasma triglycerides in another recent genome-wide scan (33), although LD information was not available between our SNP (rs7777102) and theirs (rs3812316).

The GCKR gene product inhibits glucokinase in liver and pancreatic islet cells and is considered a susceptibility gene candidate for a form of maturity-onset diabetes of the young. Polymorphisms of GCKR have been associated with a number of metabolic and CVD traits, including HDL, triglycerides, glucose, insulin resistance, renal function, C-reactive protein and protein C levels. GWASs have shown that the GCKR SNP identified in our study, rs1260326, is also associated with triglyceride levels (32), chronic kidney disease (39) and 2 h post-load glucose levels (40).

This study has several notable strengths, including a large sample size for both EA and AA participants. The IBC custom chip enables a more comprehensive study of lower frequency variants in comparison to GWAS. When comparing findings between populations, it is important to note differences in findings for individual loci may be due to different LD patterns, lower statistical power in AAs due to a smaller sample size, difference in minor allele frequency or effect modification by genetic or environmental factors.

In summary, we have confirmed and refined known loci, and identified several new genomic loci associated with FVII levels, including APOA5 and HNF4A in EAs. The chromosome 13q34 F7–F10 gene region contains multiple variants independently associated with higher and lower FVII in both EAs and AAs. Of other candidate regions associated with FVII in EA, we found evidence for replication of the chromosome 20 PROCR gene Ser219Gly and GCKR associations with higher FVII levels in AAs; on the other hand, FVII association with polymorphisms within the MS4A2, APOA5 and HNF4A gene regions were not evident in AAs. Larger samples and investigation of lower frequency variants may be required to identify additional FVII-associated loci in AAs.

MATERIALS AND METHODS

CARe discovery population

The purpose of the CARe (Candidate Gene Association Resource) consortium is to bring together deeply phenotyped prospective cohort studies to increase power for genetic association scans of CVD and other disorders (41). Cohorts included in this analysis of FVII levels are ARIC, CARDIA, the CHS and the FHS. Further details of the participating CARe studies are reported in the Supplementary Materials.

To be included in the current analysis, the participant must have provided informed consent for genetic testing and had FVII antigen levels or activity levels measured. Exclusion criteria included use of the anticoagulant warfarin. In ARIC, CHS and CARDIA (year 7 exam), FVII activity levels were measured in citrated plasma by assaying clotting time using FVII-deficient plasma (see (14), for ARIC and CHS; see (42) for CARDIA). In FHS and CARDIA (year 20 exam), FVII antigen levels were measured in ethylenediaminetetraacetic acid plasma by Enzyme-linked immunosorbent assay (Diagnostica Stago) (14,20). As CARe study participants were examined at multiple time points, covariate values that were included in the analyses were selected from the visit closest to the time of FVII measurement. After excluding study participants based on low genotyping quality control (see further information below), there were a total of 16 352 EA and 3948 AA CARe study participants available for the analysis.

Genotyping

Both AAs and EAs in CARe were genotyped using the 50 000 SNP, cardiovascular gene-centric IBC SNP array, described by Keating et al. (43). The IBC SNP array includes 49 320 SNPs selected across ~2100 candidate loci for CVD. Of the five genomic regions previously reported to be associated with FVII in the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium (GCKR, ADH4, MS4A6A, F7 and PROCR), all except ADH4 are represented on the IBC array. Genotyping was done at the Broad Institute (Cambridge, MA, USA). Criteria for DNA sample exclusion based on genotype data included sex mismatch, discordance among duplicate samples or sample call rate <95%. For each set of duplicates or monozygotic twins, data from the sample with the highest genotyping call rate were retained. SNPs were excluded when monomorphic, the call rate was <95% or Hardy–Weinberg equilibrium (HWE) was P < 10⁻⁵ in EAs. Given the genetic admixture in AAs, there was no HWE filter used for AAs. After these exclusions were applied, data remained on 47 539 working SNPs.

Additional SNPs were imputed from HapMap Phases 1 and 2 using MACH (44–45), as described in detail under Supplemental Materials. The total number of SNPs was ~250 000 after adding the imputed SNPs. We reported the genotyped results from the raw genotype data when a SNP was both genotyped and imputed.

Data analysis

Initial IBC-array-wide association analyses were conducted at the Broad Institute, and subsequent analyses were conducted at UNC-Chapel Hill. Family structure (where applicable) was taken into account during the association analysis, using a linear mixed-effects model (46). To generate a normally distributed outcome variable for use in association analyses, we stratified the data by sex, race and cohort, and modeled FVII levels as a function of age and study site (where applicable). We inverse normally transformed the residuals, and recombined the data from race and sex strata within each cohort before conducting the association analyses. For both genotyped and imputed SNPs, we assumed an additive genetic model. Consistent with this assumption, we used dosage values (i.e. a value between 0.0 and 2.0 calculated as the expected number of copies of the reference allele based on the posterior probability of each of the three possible genotypes) in the regression model implemented in PLINK (for cohorts with unrelated individuals) or the linear mixed-effects model implemented in genome-wide association study in families (for cohorts with related individuals) (46,47). To control for potential population stratification, we included the top 10 principal components, calculated in Eigenstrat, as covariates in the association analysis (48). Principal components were calculated using the cleaned CARe IBC genotype data and the HapMap populations (CEU, YRI, CHB + JPT) as reference (seed) populations.

After obtaining cohort- and race-specific results, we used the METAL software to meta-analyze the data within race (49). We used a fixed-effect Z-score approach when combining results from studies that analyzed FVII activity and antigen levels. Meta-analyses were based on P-values and weighted by sample sizes. We also calculated summary effect estimates and standard errors by meta-analyzing the three studies which measured FVII activity levels. As a supplementary analysis, we compared effect estimates based on FVII activity versus FVII antigen levels.

For the IBC array, the effective number of independent tests was calculated to be 26 482 for AA and 20 544 for EA, based on LD between markers on the array (50). To maintain an overall type 1 error rate of 5%, a statistical threshold of α = 2 × 10⁻⁶ (0.05/25000) was thus used to declare array-wide (experiment-wide) significance.

Haplotype analysis was conducted using the R-based software package HAPLO.STATS for the known haplotypes in the F7 gene represented by four SNPs (rs561241, rs2274030, rs3093265 and rs6046) that tag common patterns of LD in Caucasians (5), and an additional SNP that defined another haplotype in AAs in our data (rs3093230). The association of each haplotype with FVII levels, relative to the most common referent haplotype, was estimated using generalized linear models (haplo.glm). The EA and AA samples had the same referent haplotype. Haplotype analyses were restricted to the ARIC study (the largest sample) in EAs and AAs.

Conditional analyses

For genomic regions which initially revealed more than one significant SNP in the association scan, we conducted conditional regression analyses to estimate the number of independently associated SNPs. Conditional analyses were performed iteratively in a forward stepwise manner, with a meta-analysis done at each stage. We began by conditioning on the most significant SNP in the region (or a closely linked proxy) by including it as a covariate in the linear regression model, in each population. We determined whether any other SNPs in that region remained significant (defined as P < 2 × 10⁻⁶) after meta-analysis. If some results were still significant, we took the most significant of those SNPs, and in the next round, conditioned on the top SNPs from both rounds. This continued until there were no longer any additional significant SNPs. To estimate the cumulative proportion of variance in the sex- and age-adjusted, normalized FVII levels explained by the associated SNPs, we obtained the sum of the partial R² for each of the associated SNPs when all were simultaneously included in the regression model, adjusting for 10 PCs. All conditional analyses were conducted in SAS v. 9.2 (Cary, NC, USA).

EA replication cohort: Whitehall II Study

Replication analyses of novel FVII association findings for EAs were performed in the Whitehall II study. Details of this cohort study are provided in ref. (51). Briefly, 10 308 UK participants (70% men) were recruited between 1985 and 1988 of whom >6156 participants had blood samples for DNA collected in 2002–2004. After filtering based on call rates and removal of ancestry outliers based on principal component analysis, custom genotyping for the IBC array (48 032 SNPs) was available on 4289 individual of white ancestry. Fasting blood samples were measured for FVII activity in 1991–1993 by the method of Brozovic et al. (52). Data analysis of inverse normalized FVII activity levels was performed similarly to CARe using PLINK.

AA replication cohort: WHI

For replication of association findings in AAs, we used genotyped and imputed SNPs available through a GWAS performed in the AA participants from the WHI. Of the 8421 AA women in the WHI-SHARe project who were genotyped using the Affymetrix 6.0 chip, 984 of these women also had FVII antigen and activity levels measured. Following the primary data analysis plan, FVII association analyses were performed, using PLINK for genotyped SNPs, and ProbAbel for imputed SNPs (53).

SUPPLEMENTARY MATERIAL

Supplementary Material is available at HMG online.

Conflict of Interest statement. None declared.

FUNDING

This work was supported by the National Heart, Lung and Blood Institute of the National Institutes of Health (#5T32HL007055-34); the Heart Lung and Blood Institute of the National Institutes of Health (HL36310, HL71862-06); the National Institutes of Health (HHSN268200625226C, N01HC65226-6-0-0); the National Institute of Aging of the National Institutes of Health (AG032136, AG034454, AG13196); the Agency for Health Care Policy Research (HS06516); the Medical Research Council; the Health and Safety Executive; and the John D. and Catherine T. MacArthur Foundation Research Networks on Successful Midlife Development and Socioeconomic Status and Health and the BUPA Foundation. D.Z. is funded by the UCL Genetics Institute and S.E.H. by the British Heart Foundation.

Supplementary Material

Supplementary Data

Click here to view.

REFERENCES

1. Lippi G., Favaloro E.J., Franchini M., Guidi G.C. Milestones and perspectives in coagulation and hemostasis. Semin. Thromb. Hemost. 2009;35:9–22. [PubMed]

2. Meade T.W., Mellows S., Brozovic M., Miller G.J., Chakrabarti R.R., North W.R., Haines A.P., Stirling Y., Imeson J.D., Thompson S.G. Haemostatic function and ischaemic heart disease: principal results of the Northwick Park Heart Study. Lancet. 1986;2:533–537. [PubMed]

3. Wu K.K., Folsom A.R., Heiss G., Davis C.E., Conlan M.G., Barnes R. Association of coagulation factors and inhibitors with carotid artery atherosclerosis. Early results of the Atherosclerosis Risk in Communities (ARIC) Study. Ann. Epidemiol. 1992;2:471–480. [PubMed]

4. Junker R., Heinrich J., Schulte H., van de Loo J., Assmann G. Coagulation factor VII and the risk of coronary heart disease in healthy men. Arterioscler. Thromb. Vasc. Biol. 1997;17:1539–1544. [PubMed]

5. Reiner A.P., Carlson C.S., Rieder M.J., Siscovick D.S., Liu K., Chandler W.L., Green D., Schwartz S.M., Nickerson D.A. Coagulation factor VII gene haplotypes, obesity-related traits, and cardiovascular risk in young women. J. Thromb. Haemost. 2007;5:42–49. [PubMed]

6. Karatela R.A., Sainani G.S. Interrelationships of Factor VII activity and plasma leptin with insulin resistance in coronary heart disease. Atherosclerosis. 2010;209:235–240. [PubMed]

7. Carew J.A., Basso F., Miller G.J., Hawe E., Jackson A.A., Humphries S.E., Bauer K.A. A functional haplotype in the 5' flanking region of the factor VII gene is associated with an increased risk of coronary heart disease. J. Thromb. Haemost. 2003;1:2179–2185. [PubMed]

8. Lindman A.S., Pedersen J.I., Arnesen H., Hjerkinn E.M., Seljeflot I. Association between the factor VII haplotype containing the -402A allele and myocardial infarction in a population of elderly men at high risk for coronary heart disease. Thromb. Haemost. 2005;94:226–227. [PubMed]

9. Bozzini C., Girelli D., Bernardi F., Ferraresi P., Olivieri O., Pinotti M., Martinelli N., Manzato F., Friso S., Villa G., et al. Influence of polymorphisms in the factor VII gene promoter on activated factor VII levels and on the risk of myocardial infarction in advanced coronary atherosclerosis. Thromb. Haemost. 2004;92:541–549. [PubMed]

10. Ye Z., Liu E.H., Higgins J.P., Keavney B.D., Lowe G.D., Collins R., Danesh J. Seven haemostatic gene polymorphisms in coronary disease: meta-analysis of 66,155 cases and 91,307 controls. Lancet. 2006;367:651–658. [PubMed]

11. Souto J.C., Almasy L., Muniz-Diaz E., Soria J. M., Borrell M., Bayen L., Mateo J., Madoz P., Stone W., Blangero J., et al. Functional effects of the ABO locus polymorphism on plasma levels of von Willebrand factor, factor VIII, and activated partial thromboplastin time. Arterioscler. Thromb. Vasc. Biol. 2000;20:2024–2028. [PubMed]

12. de Lange M., Snieder H., Ariens R.A., Spector T.D., Grant P.J. The genetics of haemostasis: a twin study. Lancet. 2001;357:101–105. [PubMed]

13. Soria J.M., Almasy L., Souto J.C., Sabater-Lleal M., Fontcuberta J., Blangero J. The F7 gene and clotting factor VII levels: dissection of a human quantitative trait locus. Hum. Biol. 2005;77:561–575. [PubMed]

14. Smith N.L., Chen M.H., Dehghan A., Strachan D.P., Basu S., Soranzo N., Hayward C., Rudan I., Sabater-Lleal M., Bis J.C., et al. Novel associations of multiple genetic loci with plasma levels of factor VII, factor VIII, and von Willebrand factor: The CHARGE (Cohorts for Heart and Aging Research in Genome Epidemiology) Consortium. Circulation. 2010;121:1382–1392. [PMC free article] [PubMed]

15. Ken-Dror G., Drenos F., Humphries S.E., Talmud P.J., Hingorani A.D., Kivimaki M., Kumari M., Bauer K.A., Morrissey J.H., Ireland H.A. Haplotype and genotype effects of the F7 gene on circulating factor VII, coagulation activation markers and incident coronary heart disease in UK men. J. Thromb. Haemost. 2010;8:2394–2403. [PMC free article] [PubMed]

16. Sabater-Lleal M., Chillon M., Howard T.E., Gil E., Almasy L., Blangero J., Fontcuberta J., Soria J.M. Functional analysis of the genetic variability in the F7 gene promoter. Atherosclerosis. 2007;195:262–268. [PubMed]

17. Funk M., Endler G., Lalouschek W., Hsieh K., Schillinger M., Lang W., Mannhalter C. Factor VII gene haplotypes and risk of ischemic stroke. Clin. Chem. 2006;52:1190–1192. [PubMed]

18. Zakai N., Lange L., Longstreth W.T., Jr, O'Meara E.S., Kelley J.L., Fornage M., Nickerson D., Cushman M., Reiner A.P. Association of coagulation-related and inflammation-related genes and factor VIIc levels with stroke: the Cardiovascular Health Study. J. Thromb. Haemost. 2011;9:267–274. [PMC free article] [PubMed]

19. Maguire J.M., Thakkinstian A., Sturm J., Levi C., Lincz L., Parsons M., Whyte S., Attia J. Polymorphisms in platelet glycoprotein 1balpha and factor VII and risk of ischemic stroke: a meta-analysis. Stroke. 2008;39:1710–1716. [PubMed]

20. Green D., Foiles N., Chan C., Kang J., Schreiner P.J., Liu K. An association between clotting factor VII and carotid intima-media thickness: the CARDIA study. Stroke. 2010;41:1417–1422. [PMC free article] [PubMed]

21. Ireland H.A., Cooper J.A., Drenos F., Acharya J., Mitchell J.P., Bauer K.A., Morrissey J.H., Esnouf M.P., Humphries S.E. FVII, FVIIa, and downstream markers of extrinsic pathway activation differ by EPCR Ser219Gly variant in healthy men. Arterioscler. Thromb. Vasc. Biol. 2009;29:1968–1974. [PMC free article] [PubMed]

22. Reiner A.P., Carty C.L., Jenny N.S., Nievergelt C., Cushman M., Stearns-Kurosawa D.J., Kurosawa S., Kuller L.H., Lange L.A. PROC, PROCR and PROS1 polymorphisms, plasma anticoagulant phenotypes, and risk of cardiovascular disease and mortality in older adults: the Cardiovascular Health Study. J. Thromb. Haemost. 2008;6:1625–1632. [PMC free article] [PubMed]

23. Nayak R.C., Sen P., Ghosh S., Gopalakrishnan R., Esmon C.T., Pendurthi U.R., Rao L.V. Endothelial cell protein C receptor cellular localization and trafficking: potential functional implications. Blood. 2009;114:1974–1986. [PubMed]

24. Biguzzi E., Merati G., Liaw P.C., Bucciarelli P., Oganesyan N., Qu D., Gu J.M., Fetiveau R., Esmon C.T., Mannucci P.M., et al. A 23bp insertion in the endothelial protein C receptor (EPCR) gene impairs EPCR function. Thromb. Haemost. 2001;86:945–948. [PubMed]

25. Ireland H., Konstantoulas C.J., Cooper J.A., Hawe E., Humphries S.E., Mather H., Goodall A.H., Hogwood J., Juhan-Vague I., Yudkin J.S., et al. EPCR Ser219Gly: elevated sEPCR, prothrombin F1+2, risk for coronary heart disease, and increased sEPCR shedding in vitro. Atherosclerosis. 2005;183:283–292. [PubMed]

26. Mitropoulos K.A., Miller G.J., Reeves B.E., Wilkes H.C., Cruickshank J.K. Factor VII coagulant activity is strongly associated with the plasma concentration of large lipoprotein particles in middle-aged men. Atherosclerosis. 1989;76:203–208. [PubMed]

27. Moyer M.P., Tracy R.P., Tracy P.B., van't Veer C., Sparks C.E., Mann K.G. Plasma lipoproteins support prothrombinase and other procoagulant enzymatic complexes. Arterioscler. Thromb. Vasc. Biol. 1998;18:458–465. [PubMed]

28. Jackson A.A., Cronin K.R., Zachariah R., Carew J.A. CCAAT/enhancer-binding protein-beta participates in insulin-responsive expression of the factor VII gene. J. Biol. Chem. 2007;282:31156–31165. [PubMed]

29. Wallace C., Newhouse S.J., Braund P., Zhang F., Tobin M., Falchi M., Ahmadi K., Dobson R.J., Marcano A.C., Hajat C., et al. Genome-wide association study identifies genes for biomarkers of cardiovascular disease: serum urate and dyslipidemia. Am. J. Hum. Genet. 2008;82:139–149. [PubMed]

30. Smith E.N., Chen W., Kahonen M., Kettunen J., Lehtimaki T., Peltonen L., Raitakari O.T., Salem R.M., Schork N.J., Shaw M., et al. Longitudinal genome-wide association of cardiovascular disease risk factors in the bogalusa heart study. PLoS Genet. 2010;6:pii: e1001094. [PMC free article] [PubMed]

31. Waterworth D.M., Ricketts S.L., Song K., Chen L., Zhao J.H., Ripatti S., Aulchenko Y.S., Zhang W., Yuan X., Lim N., et al. Genetic variants influencing circulating lipid levels and risk of coronary artery disease. Arterioscler. Thromb. Vasc. Biol. 2010;30:2264–2276. [PubMed]

32. Kathiresan S., Willer C.J., Peloso G.M., Demissie S., Musunuru K., Schadt E.E., Kaplan L., Bennett D., Li Y., Tanaka T., et al. Common variants at 30 loci contribute to polygenic dyslipidemia. Nat. Genet. 2009;41:56–65. [PMC free article] [PubMed]

33. Kooner J.S., Chambers J.C., Aguilar-Salinas C.A., Hinds D.A., Hyde C.L., Warnes G.R., Gomez Perez F.J., Frazer K.A., Elliott P., Scott J., et al. Genome-wide scan identifies variation in MLXIPL associated with plasma triglycerides. Nat. Genet. 2008;40:149–151. [PubMed]

34. Willer C.J., Sanna S., Jackson A.U., Scuteri A., Bonnycastle L.L., Clarke R., Heath S.C., Timpson N.J., Najjar S.S., Stringham H.M., et al. Newly identified loci that influence lipid concentrations and risk of coronary artery disease. Nat. Genet. 2008;40:161–169. [PubMed]

35. Barrett J.C., Clayton D.G., Concannon P., Akolkar B., Cooper J.D., Erlich H.A., Julier C., Morahan G., Nerup J., Nierras C., et al. Genome-wide association study and meta-analysis find that over 40 loci affect risk of type 1 diabetes. Nat. Genet. 2009;41:703–707. [PMC free article] [PubMed]

36. Cooper J.D., Smyth D.J., Smiles A.M., Plagnol V., Walker N.M., Allen J.E., Downes K., Barrett J.C., Healy B.C., Mychaleckyj J.C., et al. Meta-analysis of genome-wide association study data identifies additional type 1 diabetes risk loci. Nat. Genet. 2008;40:1399–1401. [PMC free article] [PubMed]

37. Stahl E.A., Raychaudhuri S., Remmers E.F., Xie G., Eyre S., Thomson B.P., Li Y., Kurreeman F.A., Zhernakova A., Hinks A., et al. Genome-wide association study meta-analysis identifies seven new rheumatoid arthritis risk loci. Nat. Genet. 2010;42:508–514. [PubMed]

38. Raychaudhuri S., Remmers E.F., Lee A.T., Hackett R., Guiducci C., Burtt N.P., Gianniny L., Korman B.D., Padyukov L., Kurreeman F.A., et al. Common variants at CD40 and other loci confer risk of rheumatoid arthritis. Nat. Genet. 2008;40:1216–1223. [PMC free article] [PubMed]

39. Kottgen A., Pattaro C., Boger C.A., Fuchsberger C., Olden M., Glazer N.L., Parsa A., Gao X., Yang Q., Smith A.V., et al. New loci associated with kidney function and chronic kidney disease. Nat. Genet. 2010;42:376–384. [PMC free article] [PubMed]

40. Saxena R., Hivert M.F., Langenberg C., Tanaka T., Pankow J.S., Vollenweider P., Lyssenko V., Bouatia-Naji N., Dupuis J., Jackson A.U., et al. Genetic variation in GIPR influences the glucose and insulin responses to an oral glucose challenge. Nat. Genet. 2010;42:142–148. [PMC free article] [PubMed]

41. Musunuru K., Lettre G., Young T., Farlow D.N., Pirruccello J.P., Ejebe K.G., Keating B.J., Yang Q., Chen M.H., Lapchyk N., et al. Candidate gene association resource (CARe): design, methods, and proof of concept. Circ. Cardiovasc. Genet. 2010;3:267–275. [PMC free article] [PubMed]

42. Green D., Ruth K.J., Folsom A.R., Liu K. Hemostatic factors in the Coronary Artery Risk Development in Young Adults (CARDIA) Study. Arterioscler. Thromb. 1994;14:686–693. [PubMed]

43. Keating B.J., Tischfield S., Murray S.S., Bhangale T., Price T.S., Glessner J.T., Galver L., Barrett J.C., Grant S.F., Farlow D.N., et al. Concept, design and implementation of a cardiovascular gene-centric 50 k SNP array for large-scale genomic association studies. PLoS ONE. 2008;3:e3583. [PMC free article] [PubMed]

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

45. Li Y., Willer C., Sanna S., Abecasis G. Genotype imputation. Annu. Rev. Genomics Hum. Genet. 2009;10:387–406. [PMC free article] [PubMed]

46. Chen M-.H., Yang Q. GWAF: an R package for genome-wide association analyses with family data. Bioinformatics. 2010;26:580–581. [PMC free article] [PubMed]

47. Purcell S., Neale B., Todd-Brown K., Thomas L., Ferreira M.A., Bender D., Maller J., Sklar P., de Bakker P.I., Daly M.J., et al. PLINK: a toolset for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 2007;81:551–575. [PubMed]

48. Price A.L., Patterson N.J., Plenge R.M., Weinblatt M.E., Shadick N.A., Reich D. Principal components analysis corrects for stratification in genome-wide association studies. Nat. Genet. 2006;38:904–909. [PubMed]

49. Willer C.J., Li Y., Abecasis G.R. METAL: fast and efficient meta-analysis of genomewide association scans. Bioinformatics. 2010;26:2190–2191. [PMC free article] [PubMed]

50. Lo K.S., Wilson J.G., Lange L.A., Folsom A.R., Galarneau G., Ganesh S.K., Grant S.K., Keating B.J., McCarroll S.A., Mohler E.R., 3rd, et al. Genetic association analysis highlights new loci that modulate hematological trait variation in Caucasians and African Americans. Hum. Genet. 2011;129:307–317. [PMC free article] [PubMed]

51. Talmud P.J., Drenos F., Shah S., Shah T., Palmen J., Verzilli C., Gaunt T.R., Pallas J., Lovering R., Li K., et al. Gene-centric association signals for lipids and apolipoproteins identified via the Human CVD BeadChip. Am. J. Hum. Genet. 2009;85:628–642. [PubMed]

52. Brozovic M., Stirling Y., Harricks C., North W.R., Meade T.W. Factor VII in an industrial population. Br. J. Haematol. 1974;28:381–391. [PubMed]

53. Aulchenko Y.S., Struchalin M.V., van Duijn C.M. ProbABEL package for genome-wide association analysis of imputed data. BMC Bioinformatics. 2010;11:134. [PMC free article] [PubMed]

Articles from Human Molecular Genetics are provided here courtesy of
Oxford University Press