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

Whole-exome sequencing of 14 389 individuals from the ESP and CHARGE consortia identifies novel rare variation associated with hemostatic factors

Plasma levels of fibrinogen, coagulation factors VII and VIII and von Willebrand factor (vWF) are four intermediate phenotypes that are heritable and have been associated with the risk of clinical thrombotic events. To identify rare and low-frequency variants associated with these hemostatic factors, we conducted whole-exome sequencing in 10 860 individuals of European ancestry (EA) and 3529 African Americans (AAs) from the Cohorts for Heart and Aging Research in Genomic Epidemiology Consortium and the National Heart, Lung and Blood Institute's Exome Sequencing Project. Gene-based tests demonstrated significant associations with rare variation (minor allele frequency < 5%) in fibrinogen gamma chain (FGG) (with fibrinogen, P = 9.1 × 10-13), coagulation factor VII (F7) (with factor VII, P = 1.3 × 10-72; seven novel variants) and VWF (with factor VIII and vWF; P = 3.2 × 10-14; one novel variant). These eight novel rare variant associations were independent of the known common variants at these loci and tended to have much larger effect sizes. In addition, one of the rare novel variants in F7 was significantly associated with an increased risk of venous thromboembolism in AAs (Ile200Ser; rs141219108; P = 4.2 × 10-5). After restricting gene-based analyses to only loss-of-function variants, a novel significant association was detected and replicated between factor VIII levels and a stop-gain mutation exclusive to AAs (rs3211938) in CD36 molecule (CD36). This variant has previously been linked to dyslipidemia but not with the levels of a hemostatic factor. These efforts represent the largest integration of whole-exome sequence data from two national projects to identify genetic variation associated with plasma hemostatic factors.

Multi-phenotype analyses of hemostatic traits with cardiovascular events reveal novel genetic associations

BACKGROUND

Multi-phenotype analysis of genetically correlated phenotypes can increase the statistical power to detect loci associated with multiple traits, leading to the discovery of novel loci. This is the first study to date to comprehensively analyze the shared genetic effects within different hemostatic traits, and between these and their associated disease outcomes.

OBJECTIVES

To discover novel genetic associations by combining summary data of correlated hemostatic traits and disease events.

METHODS

Summary statistics from genome wide-association studies (GWAS) from seven hemostatic traits (factor VII [FVII], factor VIII [FVIII], von Willebrand factor [VWF] factor XI [FXI], fibrinogen, tissue plasminogen activator [tPA], plasminogen activator inhibitor 1 [PAI-1]) and three major cardiovascular (CV) events (venous thromboembolism [VTE], coronary artery disease [CAD], ischemic stroke [IS]), were combined in 27 multi-trait combinations using metaUSAT. Genetic correlations between phenotypes were calculated using Linkage Disequilibrium Score Regression (LDSC). Newly associated loci were investigated for colocalization. We considered a significance threshold of 1.85 × 10-9 obtained after applying Bonferroni correction for the number of multi-trait combinations performed (n = 27).

RESULTS

Across the 27 multi-trait analyses, we found 4 novel pleiotropic loci (XXYLT1, KNG1, SUGP1/MAU2, TBL2/MLXIPL) that were not significant in the original individual datasets, were not described in previous GWAS for the individual traits, and that presented a common associated variant between the studied phenotypes.

CONCLUSIONS

The discovery of four novel loci contributes to the understanding of the relationship between hemostasis and CV events and elucidate common genetic factors between these traits.

Elucidating mechanisms of genetic cross-disease associations at the PROCR vascular disease locus

Many individual genetic risk loci have been associated with multiple common human diseases. However, the molecular basis of this pleiotropy often remains unclear. We present an integrative approach to reveal the molecular mechanism underlying the PROCR locus, associated with lower coronary artery disease (CAD) risk but higher venous thromboembolism (VTE) risk. We identify PROCR-p.Ser219Gly as the likely causal variant at the locus and protein C as a causal factor. Using genetic analyses, human recall-by-genotype and in vitro experimentation, we demonstrate that PROCR-219Gly increases plasma levels of (activated) protein C through endothelial protein C receptor (EPCR) ectodomain shedding in endothelial cells, attenuating leukocyte-endothelial cell adhesion and vascular inflammation. We also associate PROCR-219Gly with an increased pro-thrombotic state via coagulation factor VII, a ligand of EPCR. Our study, which links PROCR-219Gly to CAD through anti-inflammatory mechanisms and to VTE through pro-thrombotic mechanisms, provides a framework to reveal the mechanisms underlying similar cross-phenotype associations.

Biochemical, molecular and clinical aspects of coagulation factor VII and its role in hemostasis and thrombosis

Activated factor VII (FVIIa), the first protease of clotting, expresses its physiological procoagulant potential only after complexing with tissue factor (TF) exposed to blood. Deep knowledge of the FVIIa-TF complex and F7 gene helps to understand the Janus-faced clinical findings associated to low or elevated FVII activity (FVIIc). Congenital FVII deficiency, the most frequent among the recessively inherited bleeding disorders, is caused by heterogeneous mutations in the F7 gene. Complete FVII deficiency causes perinatal lethality. A wide range of bleeding symptoms, from life-threatening intracranial hemorrhage to mild mucosal bleeding, is observed in patients with apparently modest differences in FVIIc levels. Though clinically relevant FVIIc threshold levels are still uncertain, effective management, including prophylaxis, has been devised, substantially improving the quality of life of patients. The exposure of TF in diseased arteries fostered investigation on the role of FVII in cardiovascular disease. FVIIc levels were found to be predictors of cardiovascular death and to be markedly associated to F7 gene variation. These genotype-phenotype relationships are among the most extensively investigated in humans. Genome-wide analyses extended association to numerous loci that, together with F7, explain >50% of FVII level plasma variance. However, the ability of F7 variation to predict thrombosis was not consistently evidenced in the numerous population studies. Main aims of this review are to highlight i) the biological and clinical information that distinguishes FVII deficiency from the other clotting disorders and ii) the impact exerted by genetically predicted FVII level variation on bleeding as well as on the thrombotic states.

A genome-wide association study identifies new loci for factor VII and implicates factor VII in ischemic stroke etiology

Factor VII (FVII) is an important component of the coagulation cascade. Few genetic loci regulating FVII activity and/or levels have been discovered to date. We conducted a meta-analysis of 9 genome-wide association studies of plasma FVII levels (7 FVII activity and 2 FVII antigen) among 27 495 participants of European and African ancestry. Each study performed ancestry-specific association analyses. Inverse variance weighted meta-analysis was performed within each ancestry group and then combined for a trans-ancestry meta-analysis. Our primary analysis included the 7 studies that measured FVII activity, and a secondary analysis included all 9 studies. We provided functional genomic validation for newly identified significant loci by silencing candidate genes in a human liver cell line (HuH7) using small-interfering RNA and then measuring F7 messenger RNA and FVII protein expression. Lastly, we used meta-analysis results to perform Mendelian randomization analysis to estimate the causal effect of FVII activity on coronary artery disease, ischemic stroke (IS), and venous thromboembolism. We identified 2 novel (REEP3 and JAZF1-AS1) and 6 known loci associated with FVII activity, explaining 19.0% of the phenotypic variance. Adding FVII antigen data to the meta-analysis did not result in the discovery of further loci. Silencing REEP3 in HuH7 cells upregulated FVII, whereas silencing JAZF1 downregulated FVII. Mendelian randomization analyses suggest that FVII activity has a positive causal effect on the risk of IS. Variants at REEP3 and JAZF1 contribute to FVII activity by regulating F7 expression levels. FVII activity appears to contribute to the etiology of IS in the general population.

Pleiotropic effects of n-6 and n-3 fatty acid-related genetic variants on circulating hemostatic variables

INTRODUCTION

Data from epidemiological studies and clinical trials suggest an influence of dietary and circulating polyunsaturated fatty acids (PUFAs) on the hemostasis profile. Genome-wide association studies (GWAS) have identified single nucleotide polymorphisms (SNPs) related to plasma PUFAs levels. We aimed to investigate whether the SNPs related to plasma PUFAs levels were also associated with plasma levels of hemostatic variables.

MATERIALS AND METHODS

We tested the associations between 9 PUFA-related SNPs and 6 hemostatic variables in 9035 European Americans (EAs) and 2702 African Americans (AAs) in the Atherosclerosis Risk in Communities (ARIC) Study. We then conducted a replication study by looking-up our novel observed associations in three published GWAS for hemostatic factors in different EA populations.

RESULTS

We observed a novel linoleic acid-related locus at the JMJD1C region associated with factor VII activity (FVIIc): rs10740118 and rs1935, Beta (p) = -1.31 (1 × 10-3) and 1.37 (5 × 10-4) in EAs, respectively, and - 1.24 (5 × 10-4) and 1.28 (3 × 10-4) in meta-analysis of EAs and AAs of ARIC. This novel association was replicated in two of three independent EA populations (p = 0.01 and 0.03 in meta-analyses). We confirmed previously reported associations at the docosapentaenoic acid-related GCKR locus with protein C and FVIIc and at JMJD1C with fibrinogen. Adjustment for plasma PUFAs did not abolish the associations between these loci and hemostatic variables.

CONCLUSIONS

Our study identified a novel association for FVIIc at JMJD1C, a histone demethylase that plays a role in DNA repair and possibly transcription regulation and RNA processing.

Associations of activated coagulation factor VII and factor VIIa-antithrombin levels with genome-wide polymorphisms and cardiovascular disease risk

ESSENTIALS

Essentials A fraction of coagulation factor VII circulates in blood as an activated protease (FVIIa). We evaluated FVIIa and FVIIa-antithrombin (FVIIa-AT) levels in the Cardiovascular Health Study. Polymorphisms in the F7 and PROCR loci were associated with FVIIa and FVIIa-AT levels. FVIIa may be an ischemic stroke risk factor in older adults and FVIIa-AT may assess mortality risk.

SUMMARY

Background A fraction of coagulation factor (F) VII circulates as an active protease (FVIIa). FVIIa also circulates as an inactivated complex with antithrombin (FVIIa-AT). Objective Evaluate associations of FVIIa and FVIIa-AT with genome-wide single nucleotide polymorphisms (SNPs) and incident coronary heart disease, ischemic stroke and mortality. Patients/Methods We measured FVIIa and FVIIa-AT in 3486 Cardiovascular Health Study (CHS) participants. We performed a genome-wide association scan for FVIIa and FVIIa-AT in European-Americans (n = 2410) and examined associations of FVII phenotypes with incident cardiovascular disease. Results In European-Americans, the most significant SNP for FVIIa and FVIIa-AT was rs1755685 in the F7 promoter region on chromosome 13 (FVIIa, β = -25.9 mU mL-1 per minor allele; FVIIa-AT, β = -26.6 pm per minor allele). Phenotypes were also associated with rs867186 located in PROCR on chromosome 20 (FVIIa, β = 7.8 mU mL-1 per minor allele; FVIIa-AT, β = 9.9 per minor allele). Adjusted for risk factors, a one standard deviation higher FVIIa was associated with increased risk of ischemic stroke (hazard ratio [HR], 1.12; 95% confidence interval [CI], 1.01, 1.23). Higher FVIIa-AT was associated with mortality from all causes (HR, 1.08; 95% CI, 1.03, 1.12). Among European-American CHS participants the rs1755685 minor allele was associated with lower ischemic stroke (HR, 0.69; 95% CI, 0.54, 0.88), but this association was not replicated in a larger multi-cohort analysis. Conclusions The results support the importance of the F7 and PROCR loci in variation in circulating FVIIa and FVIIa-AT. The findings suggest FVIIa is a risk factor for ischemic stroke in older adults, whereas higher FVIIa-AT may reflect mortality risk.

Genome-wide association study with additional genetic and post-transcriptional analyses reveals novel regulators of plasma factor XI levels

Coagulation factor XI (FXI) has become increasingly interesting for its role in pathogenesis of thrombosis. While elevated plasma levels of FXI have been associated with venous thromboembolism and ischemic stroke, its deficiency is associated with mild bleeding. We aimed to determine novel genetic and post-transcriptional plasma FXI regulators.We performed a genome-wide association study (GWAS) for plasma FXI levels, using novel data imputed to the 1000 Genomes reference panel. Individual GWAS analyses, including a total of 16,169 European individuals from the ARIC, GHS, MARTHA and PROCARDIS studies, were meta-analysed and further replicated in 2,045 individuals from the F5L family, GAIT2 and MEGA studies. Additional association with activated partial thromboplastin time (aPTT) was tested for the top SNPs. In addition, a study on the effect of miRNA on FXI regulation was performed using in silico prediction tools and in vitro luciferase assays.Three loci showed robust, replicating association with circulating FXI levels: KNG1 (rs710446, P-value = 2.07 × 10-302), F11 (rs4253417, P-value  = 2.86 × 10-193), and a novel association in GCKR (rs780094, P-value  = 3.56 ×10-09), here for the first time implicated in FXI regulation. The two first SNPs (rs710446 and rs4253417) also associated with aPTT. Conditional and haplotype analyses demonstrated a complex association signal, with additional novel SNPs modulating plasma FXI levels in both the F11 and KNG1 loci. Finally, eight miRNAs were predicted to bind F11 mRNA. Over-expression of either miR-145 or miR-181 significantly reduced the luciferase activity in cells transfected with a plasmid containing FXI-3'UTR.These results should open the door to new therapeutic targets for thrombosis prevention.

Modulation of Glucokinase Regulatory Protein: A Double-Edged Sword?

The continuous search for drugs targeting type 2 diabetes mellitus (T2DM) has led to the identification of small molecules that disrupt the binding between glucokinase and glucokinase regulatory protein (GKRP). Although mice studies are encouraging, it will take years before these disruptors can be introduced to T2DM patients. Recently, genome-wide association studies (GWASs) have shown that variants in the gene encoding GKRP protect against T2DM and kidney disease but predispose to gout, nonalcoholic fatty liver disease, and dyslipidemia. These genetic data, together with previous experience with systemic and hepatospecific glucokinase activators, provide insight into the anticipated efficacy and safety of small-molecule disruptors in humans. Interestingly, they suggest that the opposite--enhanced GKRP-glucokinase binding--could be beneficial in selected patients.

Exome Genotyping Identifies Pleiotropic Variants Associated with Red Blood Cell Traits

Red blood cell (RBC) traits are important heritable clinical biomarkers and modifiers of disease severity. To identify coding genetic variants associated with these traits, we conducted meta-analyses of seven RBC phenotypes in 130,273 multi-ethnic individuals from studies genotyped on an exome array. After conditional analyses and replication in 27,480 independent individuals, we identified 16 new RBC variants. We found low-frequency missense variants in MAP1A (rs55707100, minor allele frequency [MAF] = 3.3%, p = 2 × 10(-10) for hemoglobin [HGB]) and HNF4A (rs1800961, MAF = 2.4%, p < 3 × 10(-8) for hematocrit [HCT] and HGB). In African Americans, we identified a nonsense variant in CD36 associated with higher RBC distribution width (rs3211938, MAF = 8.7%, p = 7 × 10(-11)) and showed that it is associated with lower CD36 expression and strong allelic imbalance in ex vivo differentiated human erythroblasts. We also identified a rare missense variant in ALAS2 (rs201062903, MAF = 0.2%) associated with lower mean corpuscular volume and mean corpuscular hemoglobin (p < 8 × 10(-9)). Mendelian mutations in ALAS2 are a cause of sideroblastic anemia and erythropoietic protoporphyria. Gene-based testing highlighted three rare missense variants in PKLR, a gene mutated in Mendelian non-spherocytic hemolytic anemia, associated with HGB and HCT (SKAT p < 8 × 10(-7)). These rare, low-frequency, and common RBC variants showed pleiotropy, being also associated with platelet, white blood cell, and lipid traits. Our association results and functional annotation suggest the involvement of new genes in human erythropoiesis. We also confirm that rare and low-frequency variants play a role in the architecture of complex human traits, although their phenotypic effect is generally smaller than originally anticipated.

Rare and low-frequency variants and their association with plasma levels of fibrinogen, FVII, FVIII, and vWF

Fibrinogen, coagulation factor VII (FVII), and factor VIII (FVIII) and its carrier von Willebrand factor (vWF) play key roles in hemostasis. Previously identified common variants explain only a small fraction of the trait heritabilities, and additional variations may be explained by associations with rarer variants with larger effects. The aim of this study was to identify low-frequency (minor allele frequency [MAF] ≥0.01 and <0.05) and rare (MAF <0.01) variants that influence plasma concentrations of these 4 hemostatic factors by meta-analyzing exome chip data from up to 76,000 participants of 4 ancestries. We identified 12 novel associations of low-frequency (n = 2) and rare (n = 10) variants across the fibrinogen, FVII, FVIII, and vWF traits that were independent of previously identified associations. Novel loci were found within previously reported genes and had effect sizes much larger than and independent of previously identified common variants. In addition, associations at KCNT1, HID1, and KATNB1 identified new candidate genes related to hemostasis for follow-up replication and functional genomic analysis. Newly identified low-frequency and rare-variant associations accounted for modest amounts of trait variance and therefore are unlikely to increase predicted trait heritability but provide new information for understanding individual variation in hemostasis pathways.

Amino acid change in the carbohydrate response element binding protein is associated with lower triglycerides and myocardial infarction incidence depending on level of adherence to the Mediterranean diet in the PREDIMED trial

BACKGROUND

A variant (rs3812316, C771G, and Gln241His) in the MLXIPL (Max-like protein X interacting protein-like) gene encoding the carbohydrate response element binding protein has been associated with lower triglycerides. However, its association with cardiovascular diseases and gene-diet interactions modulating these traits are unknown.

METHODS AND RESULTS

We studied 7166 participants in the PREvención with DIeta MEDiterránea trial testing a Mediterranean diet (MedDiet) intervention versus a control diet for cardiovascular prevention, with a median follow-up of 4.8 years. Diet, lipids, MLXIPL polymorphisms, and cardiovascular events were assessed. Data were analyzed at baseline and longitudinally. We used multivariable-adjusted Cox regression to estimate hazard ratios for cardiovascular outcomes. The MLXIPL-rs3812316 was associated with lower baseline triglycerides (P=5.5×10(-5)) and lower hypertriglyceridemia (odds ratio, 0.73; 95% confidence interval [CI], 0.63-0.85; P=1.4×10(-6) in G-carriers versus CC). This association was modulated by baseline adherence to MedDiet. When adherence to MedDiet was high, the protection was stronger (odds ratio, 0.63; 95% CI, 0.51-0.77; P=8.6×10(-6)) than when adherence to MedDiet was low (odds ratio, 0.88; 95% CI, 0.70-1.09; P=0.219). Throughout the follow-up, both the MLXIPL-rs3812316 (P=3.8×10(-6)) and the MedDiet intervention (P=0.030) were significantly associated with decreased triglycerides. Likewise in G-carriers MedDiet intervention was associated with greater total cardiovascular risk reduction and specifically for myocardial infarction. In the MedDiet, but not in the control group, we observed lower myocardial infarction incidence in G-carriers versus CC (hazard ratios, 0.34; 95% CI, 0.12-0.93; P=0.036 and 0.90; 95% CI, 0.35-2.33; P=0.830, respectively).

CONCLUSIONS

Our novel results suggest that MedDiet enhances the triglyceride-lowering effect of the MLXIPL-rs3812316 variant and strengthens its protective effect on myocardial infarction incidence.

CLINICAL TRIAL REGISTRATION

URL: www.controlled-trials.com. Unique Identifier: ISRCTN35739639.

Gene-centric association signals for haemostasis and thrombosis traits identified with the HumanCVD BeadChip

Coagulation phenotypes show strong intercorrelations, affect cardiovascular disease risk and are influenced by genetic variants. The objective of this study was to search for novel genetic variants influencing the following coagulation phenotypes: factor VII levels, fibrinogen levels, plasma viscosity and platelet count. We genotyped the British Women's Heart and Health Study (n=3,445) and the Whitehall II study (n=5,059) using the Illumina HumanCVD BeadArray to investigate genetic associations and pleiotropy. In addition to previously reported associations (SH2B3, F7/F10, PROCR, GCKR, FGA/FGB/FGG, IL5), we identified novel associations at GRK5 (rs10128498, p=1.30x10(-6)), GCKR (rs1260326, p=1.63x10(-6)), ZNF259-APOA5 (rs651821, p=7.17x10(-6)) with plasma viscosity; and at CSF1 (rs333948, p=8.88x10(-6)) with platelet count. A pleiotropic effect was identified in GCKR which associated with factor VII (p=2.16x10(-7)) and plasma viscosity (p=1.63x10(-6)), and, to a lesser extent, ZNF259-APOA5 which also associated with factor VII and fibrinogen (p<1.00x10-²) and plasma viscosity (p<1.00x10(-5)). Triglyceride associated variants were overrepresented in factor VII and plasma viscosity associations. Adjusting for triglyceride levels resulted in attenuation of associations at the GCKR and ZNF259-APOA5 loci. In addition to confirming previously reported associations, we identified four single nucleotide polymorphisms (SNPs) associated with plasma viscosity and platelet count and found evidence of pleiotropic effects with SNPs in GCKR and ZNF259-APOA5. These triglyceride-associated, pleiotropic SNPs suggest a possible causal role for triglycerides in coagulation.

Associations between a polymorphism in the pleiotropic GCKR and Age-related phenotypes: the HALCyon programme

BACKGROUND

The glucokinase regulatory protein encoded by GCKR plays an important role in glucose metabolism and a single nucleotide polymorphism (SNP) rs1260326 (P446L) in the gene has been associated with several age-related biomarkers, including triglycerides, glucose, insulin and apolipoproteins. However, associations between SNPs in the gene and other ageing phenotypes such as cognitive and physical capability have not been reported.

METHODS

As part of the Healthy Ageing across the Life Course (HALCyon) collaborative research programme, men and women from five UK cohorts aged between 44 and 90+ years were genotyped for rs1260326. Meta-analysis was used to pool within-study genotypic associations between the SNP and several age-related phenotypes, including body mass index (BMI), blood lipid levels, lung function, and cognitive and physical capability.

RESULTS

We confirm the associations between the minor allele of the SNP and higher triglycerides and lower glucose levels. We also observed a triglyceride-independent association between the minor allele and lower BMI (pooled beta on z-score= -0.04, p-value=0.0001, n=16,251). Furthermore, there was some evidence for gene-environment interactions, including physical activity attenuating the effects on triglycerides. However, no associations were observed with measures of cognitive and physical capability.

CONCLUSION

Findings from middle-aged to older adults confirm associations between rs1260326 GCKR and triglycerides and glucose, suggest possible gene-environment interactions, but do not provide evidence that its relevance extends to cognitive and physical capability.