A genetic association study of circulating coagulation factor VIII and von Willebrand factor levels

ABSTRACT

Coagulation factor VIII (FVIII) and its carrier protein von Willebrand factor (VWF) are critical to coagulation and platelet aggregation. We leveraged whole-genome sequence data from the Trans-Omics for Precision Medicine (TOPMed) program along with TOPMed-based imputation of genotypes in additional samples to identify genetic associations with circulating FVIII and VWF levels in a single-variant meta-analysis, including up to 45 289 participants. Gene-based aggregate tests were implemented in TOPMed. We identified 3 candidate causal genes and tested their functional effect on FVIII release from human liver endothelial cells (HLECs) and VWF release from human umbilical vein endothelial cells. Mendelian randomization was also performed to provide evidence for causal associations of FVIII and VWF with thrombotic outcomes. We identified associations (P < 5 × 10-9) at 7 new loci for FVIII (ST3GAL4, CLEC4M, B3GNT2, ASGR1, F12, KNG1, and TREM1/NCR2) and 1 for VWF (B3GNT2). VWF, ABO, and STAB2 were associated with FVIII and VWF in gene-based analyses. Multiphenotype analysis of FVIII and VWF identified another 3 new loci, including PDIA3. Silencing of B3GNT2 and the previously reported CD36 gene decreased release of FVIII by HLECs, whereas silencing of B3GNT2, CD36, and PDIA3 decreased release of VWF by HVECs. Mendelian randomization supports causal association of higher FVIII and VWF with increased risk of thrombotic outcomes. Seven new loci were identified for FVIII and 1 for VWF, with evidence supporting causal associations of FVIII and VWF with thrombotic outcomes. B3GNT2, CD36, and PDIA3 modulate the release of FVIII and/or VWF in vitro.

A multitrait genetic study of hemostatic factors and hemorrhagic transformation after stroke treatment

BACKGROUND

Thrombolytic recombinant tissue plasminogen activator (r-tPA) treatment is the only pharmacologic intervention available in the ischemic stroke acute phase. This treatment is associated with an increased risk of intracerebral hemorrhages, known as hemorrhagic transformations (HTs), which worsen the patient's prognosis.

OBJECTIVES

To investigate the association between genetically determined natural hemostatic factors' levels and increased risk of HT after r-tPA treatment.

METHODS

Using data from genome-wide association studies on the risk of HT after r-tPA treatment and data on 7 hemostatic factors (factor [F]VII, FVIII, von Willebrand factor [VWF], FXI, fibrinogen, plasminogen activator inhibitor-1, and tissue plasminogen activator), we performed local and global genetic correlation estimation multitrait analyses and colocalization and 2-sample Mendelian randomization analyses between hemostatic factors and HT.

RESULTS

Local correlations identified a genomic region on chromosome 16 with shared covariance: fibrinogen-HT, P = 2.45 × 10-11. Multitrait analysis between fibrinogen-HT revealed 3 loci that simultaneously regulate circulating levels of fibrinogen and risk of HT: rs56026866 (PLXND1), P = 8.80 × 10-10; rs1421067 (CHD9), P = 1.81 × 10-14; and rs34780449, near ROBO1 gene, P = 1.64 × 10-8. Multitrait analysis between VWF-HT showed a novel common association regulating VWF and risk of HT after r-tPA at rs10942300 (ZNF366), P = 1.81 × 10-14. Mendelian randomization analysis did not find significant causal associations, although a nominal association was observed for FXI-HT (inverse-variance weighted estimate [SE], 0.07 [-0.29 to 0.00]; odds ratio, 0.87; 95% CI, 0.75-1.00; raw P = .05).

CONCLUSION

We identified 4 shared loci between hemostatic factors and HT after r-tPA treatment, suggesting common regulatory mechanisms between fibrinogen and VWF levels and HT. Further research to determine a possible mediating effect of fibrinogen on HT risk is needed.

Genome-wide association meta-analysis identifies risk loci for abdominal aortic aneurysm and highlights PCSK9 as a therapeutic target

Abdominal aortic aneurysm (AAA) is a common disease with substantial heritability. In this study, we performed a genome-wide association meta-analysis from 14 discovery cohorts and uncovered 141 independent associations, including 97 previously unreported loci. A polygenic risk score derived from meta-analysis explained AAA risk beyond clinical risk factors. Genes at AAA risk loci indicate involvement of lipid metabolism, vascular development and remodeling, extracellular matrix dysregulation and inflammation as key mechanisms in AAA pathogenesis. These genes also indicate overlap between the development of AAA and other monogenic aortopathies, particularly via transforming growth factor β signaling. Motivated by the strong evidence for the role of lipid metabolism in AAA, we used Mendelian randomization to establish the central role of nonhigh-density lipoprotein cholesterol in AAA and identified the opportunity for repurposing of proprotein convertase, subtilisin/kexin-type 9 (PCSK9) inhibitors. This was supported by a study demonstrating that PCSK9 loss of function prevented the development of AAA in a preclinical mouse model.

Polygenic risk of major depressive disorder as a risk factor for venous thromboembolism

Major depressive disorder (MDD), bipolar disorder (BD), and schizophrenia (SCZ) are associated with an increased risk of cardiovascular diseases, including venous thromboembolism (VTE). The reasons for this are complex and include obesity, smoking, and use of hormones and psychotropic medications. Genetic studies have increasingly provided evidence of the shared genetic risk of psychiatric and cardiometabolic illnesses. This study aimed to determine whether a genetic predisposition to MDD, BD, or SCZ is associated with an increased risk of VTE. Genetic correlations using the largest genome-wide genetic meta-analyses summary statistics for MDD, BD, and SCZ (Psychiatric Genetics Consortium) and a recent genome-wide genetic meta-analysis of VTE (INVENT Consortium) demonstrated a positive association between VTE and MDD but not BD or SCZ. The same summary statistics were used to construct polygenic risk scores for MDD, BD, and SCZ in UK Biobank participants of self-reported White British ancestry. These were assessed for impact on self-reported VTE risk (10 786 cases, 285 124 controls), using logistic regression, in sex-specific and sex-combined analyses. We identified significant positive associations between polygenic risk for MDD and the risk of VTE in men, women, and sex-combined analyses, independent of the known risk factors. Secondary analyses demonstrated that this association was not driven by those with lifetime experience of mental illness. Meta-analyses of individual data from 6 additional independent cohorts replicated the sex-combined association. This report provides evidence for shared biological mechanisms leading to MDD and VTE and suggests that, in the absence of genetic data, a family history of MDD might be considered when assessing the risk of VTE.

Antithrombin, Protein C, and Protein S: Genome and Transcriptome-Wide Association Studies Identify 7 Novel Loci Regulating Plasma Levels

BACKGROUND

Antithrombin, PC (protein C), and PS (protein S) are circulating natural anticoagulant proteins that regulate hemostasis and of which partial deficiencies are causes of venous thromboembolism. Previous genetic association studies involving antithrombin, PC, and PS were limited by modest sample sizes or by being restricted to candidate genes. In the setting of the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium, we meta-analyzed across ancestries the results from 10 genome-wide association studies of plasma levels of antithrombin, PC, PS free, and PS total.

METHODS

Study participants were of European and African ancestries, and genotype data were imputed to TOPMed, a dense multiancestry reference panel. Each of the 10 studies conducted a genome-wide association studies for each phenotype and summary results were meta-analyzed, stratified by ancestry. Analysis of antithrombin included 25 243 European ancestry and 2688 African ancestry participants, PC analysis included 16 597 European ancestry and 2688 African ancestry participants, PSF and PST analysis included 4113 and 6409 European ancestry participants. We also conducted transcriptome-wide association analyses and multiphenotype analysis to discover additional associations. Novel genome-wide association studies and transcriptome-wide association analyses findings were validated by in vitro functional experiments. Mendelian randomization was performed to assess the causal relationship between these proteins and cardiovascular outcomes.

RESULTS

Genome-wide association studies meta-analyses identified 4 newly associated loci: 3 with antithrombin levels (GCKR, BAZ1B, and HP-TXNL4B) and 1 with PS levels (ORM1-ORM2). transcriptome-wide association analyses identified 3 newly associated genes: 1 with antithrombin level (FCGRT), 1 with PC (GOLM2), and 1 with PS (MYL7). In addition, we replicated 7 independent loci reported in previous studies. Functional experiments provided evidence for the involvement of GCKR, SNX17, and HP genes in antithrombin regulation.

CONCLUSIONS

The use of larger sample sizes, diverse populations, and a denser imputation reference panel allowed the detection of 7 novel genomic loci associated with plasma antithrombin, PC, and PS levels.

Whole genome analysis of plasma fibrinogen reveals population-differentiated genetic regulators with putative liver roles

Genetic studies have identified numerous regions associated with plasma fibrinogen levels in Europeans, yet missing heritability and limited inclusion of non-Europeans necessitates further studies with improved power and sensitivity. Compared with array-based genotyping, whole genome sequencing (WGS) data provides better coverage of the genome and better representation of non-European variants. To better understand the genetic landscape regulating plasma fibrinogen levels, we meta-analyzed WGS data from the NHLBI's Trans-Omics for Precision Medicine (TOPMed) program (n=32,572), with array-based genotype data from the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium (n=131,340) imputed to the TOPMed or Haplotype Reference Consortium panel. We identified 18 loci that have not been identified in prior genetic studies of fibrinogen. Of these, four are driven by common variants of small effect with reported MAF at least 10% higher in African populations. Three ( SERPINA1, ZFP36L2 , and TLR10) signals contain predicted deleterious missense variants. Two loci, SOCS3 and HPN , each harbor two conditionally distinct, non-coding variants. The gene region encoding the protein chain subunits ( FGG;FGB;FGA ), contains 7 distinct signals, including one novel signal driven by rs28577061, a variant common (MAF=0.180) in African reference panels but extremely rare (MAF=0.008) in Europeans. Through phenome-wide association studies in the VA Million Veteran Program, we found associations between fibrinogen polygenic risk scores and thrombotic and inflammatory disease phenotypes, including an association with gout. Our findings demonstrate the utility of WGS to augment genetic discovery in diverse populations and offer new insights for putative mechanisms of fibrinogen regulation.

Key Points

Largest and most diverse genetic study of plasma fibrinogen identifies 54 regions (18 novel), housing 69 conditionally distinct variants (20 novel).Sufficient power achieved to identify signal driven by African population variant.Links to (1) liver enzyme, blood cell and lipid genetic signals, (2) liver regulatory elements, and (3) thrombotic and inflammatory disease.

DNA methylation analysis is used to identify novel genetic loci associated with circulating fibrinogen levels in blood

BACKGROUND

Fibrinogen plays an essential role in blood coagulation and inflammation. Circulating fibrinogen levels may be determined based on interindividual differences in DNA methylation at cytosine-phosphate-guanine (CpG) sites and vice versa.

OBJECTIVES

To perform an EWAS to examine an association between blood DNA methylation levels and circulating fibrinogen levels to better understand its biological and pathophysiological actions.

METHODS

We performed an epigenome-wide association study of circulating fibrinogen levels in 18 037 White, Black, American Indian, and Hispanic participants, representing 14 studies from the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium. Circulating leukocyte DNA methylation was measured using the Illumina 450K array in 12 904 participants and using the EPIC array in 5133 participants. In each study, an epigenome-wide association study of fibrinogen was performed using linear mixed models adjusted for potential confounders. Study-specific results were combined using array-specific meta-analysis, followed by cross-replication of epigenome-wide significant associations. We compared models with and without CRP adjustment to examine the role of inflammation.

RESULTS

We identified 208 and 87 significant CpG sites associated with fibrinogen levels from the 450K (p < 1.03 × 10-7) and EPIC arrays (p < 5.78 × 10-8), respectively. There were 78 associations from the 450K array that replicated in the EPIC array and 26 vice versa. After accounting for overlapping sites, there were 83 replicated CpG sites located in 61 loci, of which only 4 have been previously reported for fibrinogen. The examples of genes located near these CpG sites were SOCS3 and AIM2, which are involved in inflammatory pathways. The associations of all 83 replicated CpG sites were attenuated after CRP adjustment, although many remained significant.

CONCLUSION

We identified 83 CpG sites associated with circulating fibrinogen levels. These associations are partially driven by inflammatory pathways shared by both fibrinogen and CRP.

Implicating genes, pleiotropy, and sexual dimorphism at blood lipid loci through multi-ancestry meta-analysis

BACKGROUND

Genetic variants within nearly 1000 loci are known to contribute to modulation of blood lipid levels. However, the biological pathways underlying these associations are frequently unknown, limiting understanding of these findings and hindering downstream translational efforts such as drug target discovery.

RESULTS

To expand our understanding of the underlying biological pathways and mechanisms controlling blood lipid levels, we leverage a large multi-ancestry meta-analysis (N = 1,654,960) of blood lipids to prioritize putative causal genes for 2286 lipid associations using six gene prediction approaches. Using phenome-wide association (PheWAS) scans, we identify relationships of genetically predicted lipid levels to other diseases and conditions. We confirm known pleiotropic associations with cardiovascular phenotypes and determine novel associations, notably with cholelithiasis risk. We perform sex-stratified GWAS meta-analysis of lipid levels and show that 3-5% of autosomal lipid-associated loci demonstrate sex-biased effects. Finally, we report 21 novel lipid loci identified on the X chromosome. Many of the sex-biased autosomal and X chromosome lipid loci show pleiotropic associations with sex hormones, emphasizing the role of hormone regulation in lipid metabolism.

CONCLUSIONS

Taken together, our findings provide insights into the biological mechanisms through which associated variants lead to altered lipid levels and potentially cardiovascular disease risk.

Cross-Ancestry Investigation of Venous Thromboembolism Genomic Predictors

BACKGROUND

Venous thromboembolism (VTE) is a life-threatening vascular event with environmental and genetic determinants. Recent VTE genome-wide association studies (GWAS) meta-analyses involved nearly 30 000 VTE cases and identified up to 40 genetic loci associated with VTE risk, including loci not previously suspected to play a role in hemostasis. The aim of our research was to expand discovery of new genetic loci associated with VTE by using cross-ancestry genomic resources.

METHODS

We present new cross-ancestry meta-analyzed GWAS results involving up to 81 669 VTE cases from 30 studies, with replication of novel loci in independent populations and loci characterization through in silico genomic interrogations.

RESULTS

In our genetic discovery effort that included 55 330 participants with VTE (47 822 European, 6320 African, and 1188 Hispanic ancestry), we identified 48 novel associations, of which 34 were replicated after correction for multiple testing. In our combined discovery-replication analysis (81 669 VTE participants) and ancestry-stratified meta-analyses (European, African, and Hispanic), we identified another 44 novel associations, which are new candidate VTE-associated loci requiring replication. In total, across all GWAS meta-analyses, we identified 135 independent genomic loci significantly associated with VTE risk. A genetic risk score of the significantly associated loci in Europeans identified a 6-fold increase in risk for those in the top 1% of scores compared with those with average scores. We also identified 31 novel transcript associations in transcriptome-wide association studies and 8 novel candidate genes with protein quantitative-trait locus Mendelian randomization analyses. In silico interrogations of hemostasis and hematology traits and a large phenome-wide association analysis of the 135 GWAS loci provided insights to biological pathways contributing to VTE, with some loci contributing to VTE through well-characterized coagulation pathways and others providing new data on the role of hematology traits, particularly platelet function. Many of the replicated loci are outside of known or currently hypothesized pathways to thrombosis.

CONCLUSIONS

Our cross-ancestry GWAS meta-analyses identified new loci associated with VTE. These findings highlight new pathways to thrombosis and provide novel molecules that may be useful in the development of improved antithrombosis treatments.

A saturated map of common genetic variants associated with human height

Common single-nucleotide polymorphisms (SNPs) are predicted to collectively explain 40-50% of phenotypic variation in human height, but identifying the specific variants and associated regions requires huge sample sizes1. Here, using data from a genome-wide association study of 5.4 million individuals of diverse ancestries, we show that 12,111 independent SNPs that are significantly associated with height account for nearly all of the common SNP-based heritability. These SNPs are clustered within 7,209 non-overlapping genomic segments with a mean size of around 90 kb, covering about 21% of the genome. The density of independent associations varies across the genome and the regions of increased density are enriched for biologically relevant genes. In out-of-sample estimation and prediction, the 12,111 SNPs (or all SNPs in the HapMap 3 panel2) account for 40% (45%) of phenotypic variance in populations of European ancestry but only around 10-20% (14-24%) in populations of other ancestries. Effect sizes, associated regions and gene prioritization are similar across ancestries, indicating that reduced prediction accuracy is likely to be explained by linkage disequilibrium and differences in allele frequency within associated regions. Finally, we show that the relevant biological pathways are detectable with smaller sample sizes than are needed to implicate causal genes and variants. Overall, this study provides a comprehensive map of specific genomic regions that contain the vast majority of common height-associated variants. Although this map is saturated for populations of European ancestry, further research is needed to achieve equivalent saturation in other ancestries.

A multi-layer functional genomic analysis to understand noncoding genetic variation in lipids

A major challenge of genome-wide association studies (GWASs) is to translate phenotypic associations into biological insights. Here, we integrate a large GWAS on blood lipids involving 1.6 million individuals from five ancestries with a wide array of functional genomic datasets to discover regulatory mechanisms underlying lipid associations. We first prioritize lipid-associated genes with expression quantitative trait locus (eQTL) colocalizations and then add chromatin interaction data to narrow the search for functional genes. Polygenic enrichment analysis across 697 annotations from a host of tissues and cell types confirms the central role of the liver in lipid levels and highlights the selective enrichment of adipose-specific chromatin marks in high-density lipoprotein cholesterol and triglycerides. Overlapping transcription factor (TF) binding sites with lipid-associated loci identifies TFs relevant in lipid biology. In addition, we present an integrative framework to prioritize causal variants at GWAS loci, producing a comprehensive list of candidate causal genes and variants with multiple layers of functional evidence. We highlight two of the prioritized genes, CREBRF and RRBP1, which show convergent evidence across functional datasets supporting their roles in lipid biology.

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.

Integrated GWAS and Gene Expression Suggest ORM1 as a Potential Regulator of Plasma Levels of Cell-Free DNA and Thrombosis Risk

Plasma cell-free DNA (cfDNA) is a surrogate marker of neutrophil extracellular traps (NETs) that contribute to immunothrombosis. There is growing interest about the mechanisms underlying NET formation and elevated cfDNA, but little is known about the factors involved. We aimed to identify genes involved in the regulation of cfDNA levels using data from the Genetic Analysis of Idiopathic Thrombophilia (GAIT-2) Project.

Imputed genotypes, whole blood RNA-Seq data, and plasma cfDNA quantification were available for 935 of the GAIT-2 participants from 35 families with idiopathic thrombophilia. We performed heritability and GWAS analysis for cfDNA. The heritability of cfDNA was 0.26 ( p  = 3.7 × 10 ⁻⁶ ), while the GWAS identified a significant association (rs1687391, p  = 3.55 × 10 ⁻¹⁰ ) near the ORM1 gene, on chromosome 9. An eQTL (expression quantitative trait loci) analysis revealed a significant association between the lead GWAS variant and the expression of ORM1 in whole blood ( p  = 6.14 × 10 ⁻⁹ ). Additionally, ORM1 expression correlated with levels of cfDNA ( p  = 4.38 × 10 ⁻⁴ ). Finally, genetic correlation analysis between cfDNA and thrombosis identified a suggestive association ( ρ_g  = 0.43, p  = 0.089).

All in all, we show evidence of the role of ORM1 in regulating cfDNA levels in plasma, which might contribute to the susceptibility to thrombosis through mechanisms of immunothrombosis.

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.

Organ-on-a-chip technology: a novel approach to investigate cardiovascular diseases

The development of organs-on-chip (OoC) has revolutionized in vitro cell-culture experiments by allowing a better mimicry of human physiology and pathophysiology that has consequently led researchers to gain more meaningful insights into disease mechanisms. Several models of hearts-on-chips and vessels-on-chips have been demonstrated to recapitulate fundamental aspects of the human cardiovascular system in the recent past. These 2D and 3D systems include synchronized beating cardiomyocytes in hearts-on-chips and vessels-on-chips with layer-based structures and the inclusion of physiological and pathological shear stress conditions. The opportunities to discover novel targets and to perform drug testing with chip-based platforms have substantially enhanced, thanks to the utilization of patient-derived cells and precise control of their microenvironment. These organ models will provide an important asset for future approaches to personalized cardiovascular medicine and improved patient care. However, certain technical and biological challenges remain, making the global utilization of OoCs to tackle unanswered questions in cardiovascular science still rather challenging. This review article aims to introduce and summarize published work on hearts- and vessels-on chips but also to provide an outlook and perspective on how these advanced in vitro systems can be used to tailor disease models with patient-specific characteristics.

The power of genetic diversity in genome-wide association studies of lipids

Increased blood lipid levels are heritable risk factors of cardiovascular disease with varied prevalence worldwide owing to different dietary patterns and medication use1. Despite advances in prevention and treatment, in particular through reducing low-density lipoprotein cholesterol levels2, heart disease remains the leading cause of death worldwide3. Genome-wideassociation studies (GWAS) of blood lipid levels have led to important biological and clinical insights, as well as new drug targets, for cardiovascular disease. However, most previous GWAS4-23 have been conducted in European ancestry populations and may have missed genetic variants that contribute to lipid-level variation in other ancestry groups. These include differences in allele frequencies, effect sizes and linkage-disequilibrium patterns24. Here we conduct a multi-ancestry, genome-wide genetic discovery meta-analysis of lipid levels in approximately 1.65 million individuals, including 350,000 of non-European ancestries. We quantify the gain in studying non-European ancestries and provide evidence to support the expansion of recruitment of additional ancestries, even with relatively small sample sizes. We find that increasing diversity rather than studying additional individuals of European ancestry results in substantial improvements in fine-mapping functional variants and portability of polygenic prediction (evaluated in approximately 295,000 individuals from 7 ancestry groupings). Modest gains in the number of discovered loci and ancestry-specific variants were also achieved. As GWAS expand emphasis beyond the identification of genes and fundamental biology towards the use of genetic variants for preventive and precision medicine25, we anticipate that increased diversity of participants will lead to more accurate and equitable26 application of polygenic scores in clinical practice.

Effects of Thyroid Function on Hemostasis, Coagulation, and Fibrinolysis: A Mendelian Randomization Study

Background: Untreated hypothyroidism is associated with acquired von Willebrand syndrome, and hyperthyroidism is associated with increased thrombosis risk. However, the causal effects of thyroid function on hemostasis, coagulation, and fibrinolysis are unknown. Methods: In a two-sample Mendelian randomization (MR) study with genome-wide association variants, we assessed causality of genetically predicted hypothyroidism (N = 134,641), normal-range thyrotropin (TSH; N = 54,288) and free thyroxine (fT4) (N = 49,269), hyperthyroidism (N = 51,823), and thyroid peroxidase antibody positivity (N = 25,821) on coagulation (activated partial thromboplastin time, von Willebrand factor [VWF], factor VIII [FVIII], prothrombin time, factor VII, fibrinogen) and fibrinolysis (D-dimer, tissue plasminogen activator [TPA], plasminogen activator inhibitor-1) from the CHARGE Hemostasis Consortium (N = 2583-120,246). Inverse-variance-weighted random effects were the main MR analysis followed by sensitivity analyses. Two-sided p < 0.05 was nominally significant, and p < 0.0011[ = 0.05/(5 exposures × 9 outcomes)] was Bonferroni significant for the main MR analysis. Results: Genetically increased TSH was associated with decreased VWF [β(SE) = -0.020(0.006), p = 0.001] and with decreased fibrinogen [β(SE) = -0.008(0.002), p = 0.001]. Genetically increased fT4 was associated with increased VWF [β(SE) = 0.028(0.011), p = 0.012]. Genetically predicted hyperthyroidism was associated with increased VWF [β(SE) = 0.012(0.004), p = 0.006] and increased FVIII [β(SE) = 0.013(0.005), p = 0.007]. Genetically predicted hypothyroidism and hyperthyroidism were associated with decreased TPA [β(SE) = -0.009(0.024), p = 0.024] and increased TPA [β(SE) = 0.022(0.008), p = 0.008], respectively. MR sensitivity analyses showed similar direction but lower precision. Other coagulation and fibrinolytic factors were inconclusive. Conclusions: In the largest genetic studies currently available, genetically increased TSH and fT4 may be associated with decreased and increased synthesis of VWF, respectively. Since Bonferroni correction may be too conservative given the correlation between the analyzed traits, we cannot reject nominal associations of thyroid traits with coagulation or fibrinolytic factors.

FGL1 as a modulator of plasma D-dimer levels: Exome-wide marker analysis of plasma tPA, PAI-1, and D-dimer

BACKGROUND

Use of targeted exome-arrays with common, rare variants and functionally enriched variation has led to discovery of new genes contributing to population variation in risk factors. Plasminogen activator-inhibitor 1 (PAI-1), tissue plasminogen activator (tPA), and the plasma product D-dimer are important components of the fibrinolytic system. There have been few large-scale genome-wide or exome-wide studies of PAI-1, tPA, and D-dimer.

OBJECTIVES

We sought to discover new genetic loci contributing to variation in these traits using an exome-array approach.

METHODS

Cohort-level analyses and fixed effects meta-analyses of PAI-1 (n = 15 603), tPA (n = 6876,) and D-dimer (n = 19 306) from 12 cohorts of European ancestry with diverse study design were conducted, including single-variant analyses and gene-based burden testing.

RESULTS

Five variants located in NME7, FGL1, and the fibrinogen locus, all associated with D-dimer levels, achieved genome-wide significance (P < 5 × 10-8 ). Replication was sought for these 5 variants, as well as 45 well-imputed variants with P < 1 × 10-4 in the discovery using an independent cohort. Replication was observed for three out of the five significant associations, including a novel and uncommon (0.013 allele frequency) coding variant p.Trp256Leu in FGL1 (fibrinogen-like-1) with increased plasma D-dimer levels. Additionally, a candidate-gene approach revealed a suggestive association for a coding variant (rs143202684-C) in SERPINB2, and suggestive associations with consistent effect in the replication analysis include an intronic variant (rs11057830-A) in SCARB1 associated with increased D-dimer levels.

CONCLUSION

This work provides new evidence for a role of FGL1 in hemostasis.

Association between ABO haplotypes and the risk of venous thrombosis: impact on disease risk estimation

Genetic risk score (GRS) analysis is a popular approach to derive individual risk prediction models for complex diseases. In venous thrombosis (VT), such type of analysis shall integrate information at the ABO blood group locus, which is one of the major susceptibility loci. However, there is no consensus about which single nucleotide polymorphisms (SNPs) must be investigated when properly assessing association between ABO locus and VT risk. Using comprehensive haplotype analyses of ABO blood group tagging SNPs in 5425 cases and 8445 controls from 6 studies, we demonstrate that using only rs8176719 (tagging O1) to correctly assess the impact of ABO locus on VT risk is suboptimal, because 5% of rs8176719-delG carriers do not have an increased risk of developing VT. Instead, we recommend the use of 4 SNPs, rs2519093 (tagging A1), rs1053878 (A2), rs8176743 (B), and rs41302905 (O2), when assessing the impact of ABO locus on VT risk to avoid any risk misestimation. Compared with the O1 haplotype, the A2 haplotype is associated with a modest increase in VT risk (odds ratio, ∼1.2), the A1 and B haplotypes are associated with an ∼1.8-fold increased risk, whereas the O2 haplotype tends to be slightly protective (odds ratio, ∼0.80). In addition, although the A1 and B blood groups are associated with increased von Willebrand factor and factor VIII plasma levels, only the A1 blood group is associated with ICAM levels, but in an opposite direction, leaving additional avenues to be explored to fully understand the spectrum of biological effects mediated by ABO locus on cardiovascular traits.

Genome-wide association study of circulating interleukin 6 levels identifies novel loci

Interleukin 6 (IL-6) is a multifunctional cytokine with both pro- and anti-inflammatory properties with a heritability estimate of up to 61%. The circulating levels of IL-6 in blood have been associated with an increased risk of complex disease pathogenesis. We conducted a two-staged, discovery and replication meta genome-wide association study (GWAS) of circulating serum IL-6 levels comprising up to 67 428 (ndiscovery = 52 654 and nreplication = 14 774) individuals of European ancestry. The inverse variance fixed effects based discovery meta-analysis, followed by replication led to the identification of two independent loci, IL1F10/IL1RN rs6734238 on chromosome (Chr) 2q14, (Pcombined = 1.8 × 10-11), HLA-DRB1/DRB5 rs660895 on Chr6p21 (Pcombined = 1.5 × 10-10) in the combined meta-analyses of all samples. We also replicated the IL6R rs4537545 locus on Chr1q21 (Pcombined = 1.2 × 10-122). Our study identifies novel loci for circulating IL-6 levels uncovering new immunological and inflammatory pathways that may influence IL-6 pathobiology.

Association Between Hemostatic Profile and Migraine: A Mendelian Randomization Analysis

OBJECTIVE

To assess support for a causal relationship between hemostatic measures and migraine susceptibility using genetic instrumental analysis.

METHODS

Two-sample Mendelian randomization instrumental analyses leveraging available genome-wide association study (GWAS) summary statistics were applied to hemostatic measures as potentially causal for migraine and its subtypes, migraine with aura (MA) and migraine without aura (MO). Twelve blood-based measures of hemostasis were examined, including plasma level or activity of 8 hemostatic factors and 2 fibrinopeptides together with 2 hemostasis clinical tests.

RESULTS

There were significant instrumental effects between increased coagulation factor VIII activity (FVIII; odds ratio [95% confidence interval] 1.05 [1.03, 1.08]/SD, p = 6.08 × 10^-05), von Willebrand factor level (vWF; 1.05 [1.03, 1.08]/SD, p = 2.25 × 10^-06), and phosphorylated fibrinopeptide A level (1.13 [1.07, 1.19]/SD, p = 5.44 × 10^-06) with migraine susceptibility. When extended to migraine subtypes, FVIII, vWF, and phosphorylated fibrinopeptide A showed slightly stronger effects with MA than overall migraine. Fibrinogen level was inversely linked with MA (0.76 [0.64, 0.91]/SD, p = 2.32 × 10^-03) but not overall migraine. None of the hemostatic factors was linked with MO. In sensitivity analysis, effects for fibrinogen and phosphorylated fibrinopeptide A were robust, whereas independent effects of FVIII and vWF could not be distinguished, and FVIII associations were potentially affected by pleiotropy at the ABO locus. Causal effects from migraine to the hemostatic measures were not supported in reverse Mendelian randomization. However, MA was not included due to lack of instruments.

CONCLUSIONS

The findings support potential causality of increased FVIII, vWF, and phosphorylated fibrinopeptide A and decreased fibrinogen in migraine susceptibility, especially for MA, potentially revealing etiologic relationships between hemostasis and migraine.

Expression of microRNAs in human platelet-poor plasma: analysis of the factors affecting their expression and association with proximal genetic variants

To translate circulating microRNAs (miRNAs) into the clinic, a deeper understanding of the factors affecting their expression is needed. In this study, we explored the features affecting the expression of miRNAs and their genetic regulation using the expression data of 103 miRNAs obtained by qPCR in the platelet-poor plasma of 104 subjects. The principal components (PCs) of the expression of the miRNAs were associated with technical and biological features (e.g., synthetic controls or sex) and with blood cell counts. Also, the associations with proximal genetics variants were analysed. We found that haemolysis marker (dCt hsa-miR-23a-3p-hsa-miR-451a) was correlated strongly (β = 0.84, p = 2.07x10-29) with the second PC, which explained 10.1% of the overall variability. Thus, we identified haemolysis as a source of variability for miRNA expression even in mild hemolyzed samples (haemolysis marker dCt <5). In addition to hsa-miR-23a-3p and hsa-miR-451a, the miRNAs most stable and most susceptible to haemolysis were identified. Then, we discovered that the expression of miRNAs in platelet-poor plasma was not biased by any blood cell count, and thus, our results supported their role as biomarkers of tissue-specific conditions. Finally, we identified 1,323 genetic variants that corresponded to 158 miRNA expression quantitative trait loci for 14 miRNAs (FDR <0.2), which were enriched in promoter regions (p = 0.03). This enrichment corresponded to a range of specific tissues (e.g., breast or fat) although not to blood tissue, supporting the concept that the expression of circulating miRNAs is under the genetic control of different tissues.

Mendelian Randomization Analysis of Hemostatic Factors and Their Contribution to Peripheral Artery Disease-Brief Report

BACKGROUND AND OBJECTIVE

Peripheral artery disease (PAD) is the third most common form of atherosclerotic vascular disease and is characterized by significant functional disability and increased cardiovascular mortality. Recent genetic data support a role for a procoagulation protein variant, the factor V Leiden mutation, in PAD. The role of other hemostatic factors in PAD remains unknown. We evaluated the role of hemostatic factors in PAD using Mendelian randomization. Approach and Results: Two-sample Mendelian randomization to evaluate the roles of FVII (factor VII), FVIII (factor VIII), FXI (factor XI), VWF (von Willebrand factor), and fibrinogen in PAD was performed using summary statistics from GWAS for hemostatic factors performed within the Cohorts for Heart and Aging Research in the Genome Epidemiology Consortium and from GWAS performed for PAD within the Million Veteran Program. Genetically determined FVIII and VWF, but not FVII, FXI, or fibrinogen, were associated with PAD in Mendelian randomization experiments (FVIII: odds ratio, 1.41 [95% CI, 1.23-1.62], P=6.0×10^-7, VWF: odds ratio, 1.28 [95% CI, 1.07-1.52], P=0.0073). In single variant sensitivity analysis, the ABO locus was the strongest genetic instrument for both FVIII and VWF.

CONCLUSIONS

Our results suggest a role for hemostasis, and by extension, thrombosis in PAD. Further study is warranted to determine whether VWF and FVIII independently affect the biology of PAD.

The TAGA Study: A Study of Factors Determining Aortic Diameter in Families at High Risk of Abdominal Aortic Aneurysm Reveal Two New Candidate Genes

A variety of disorders are known to be related with aortic geometry, among them abdominal aortic aneurysm (AAA). This work aims to present the main determinants of abdominal aortic diameter in a new cohort of families at high risk of AAA. The Triple-A Genomic Analysis (TAGA) study comprises 407 individuals related in 12 families. Each family was collected through a proband with AAA. We calculated heritability and genetic correlations between abdominal aortic diameter and clinical parameters. A genome-wide linkage scan was performed based on 4.6 million variants. A predictive model was calculated with conditional forest. Heritability of the abdominal aortic diameter was 34%. Old age, male sex, higher height, weight, creatinine levels in serum, and better lung capacity were the best predictors of aortic diameter. Linkage analyses suggested the implication of Epidermal Growth Factor Receptor (EGFR) and Betacellulin (BTC) genes with aortic diameter. This is the first study to evaluate genetic components of variation of the aortic diameter in a population of AAA high-risk individuals. These results reveal EGFR, a gene that had been previously implicated in AAA, as a determinant of aortic diameter variation in healthy genetically enriched individuals, and might indicate that a common genetic background could determine the diameter of the aorta and future risk of AAA.

PCSK6 Is a Key Protease in the Control of Smooth Muscle Cell Function in Vascular Remodeling

RATIONALE

PCSKs (Proprotein convertase subtilisins/kexins) are a protease family with unknown functions in vasculature. Previously, we demonstrated PCSK6 upregulation in human atherosclerotic plaques associated with smooth muscle cells (SMCs), inflammation, extracellular matrix remodeling, and mitogens.

OBJECTIVE

Here, we applied a systems biology approach to gain deeper insights into the PCSK6 role in normal and diseased vessel wall.

METHODS AND RESULTS

Genetic analyses revealed association of intronic PCSK6 variant rs1531817 with maximum internal carotid intima-media thickness progression in high-cardiovascular risk subjects. This variant was linked with PCSK6 mRNA expression in healthy aortas and plaques but also with overall plaque SMA+ cell content and pericyte fraction. Increased PCSK6 expression was found in several independent human cohorts comparing atherosclerotic lesions versus healthy arteries, using transcriptomic and proteomic datasets. By immunohistochemistry, PCSK6 was localized to fibrous cap SMA+ cells and neovessels in plaques. In human, rat, and mouse intimal hyperplasia, PCSK6 was expressed by proliferating SMA+ cells and upregulated after 5 days in rat carotid balloon injury model, with positive correlation to PDGFB (platelet-derived growth factor subunit B) and MMP (matrix metalloprotease) 2/MMP14. Here, PCSK6 was shown to colocalize and cointeract with MMP2/MMP14 by in situ proximity ligation assay. Microarrays of carotid arteries from Pcsk6^-/- versus control mice revealed suppression of contractile SMC markers, extracellular matrix remodeling enzymes, and cytokines/receptors. Pcsk6^-/- mice showed reduced intimal hyperplasia response upon carotid ligation in vivo, accompanied by decreased MMP14 activation and impaired SMC outgrowth from aortic rings ex vivo. PCSK6 silencing in human SMCs in vitro leads to downregulation of contractile markers and increase in MMP2 expression. Conversely, PCSK6 overexpression increased PDGFBB (platelet-derived growth factor BB)-induced cell proliferation and particularly migration.

CONCLUSIONS

PCSK6 is a novel protease that induces SMC migration in response to PDGFB, mechanistically via modulation of contractile markers and MMP14 activation. This study establishes PCSK6 as a key regulator of SMC function in vascular remodeling. Visual Overview: An online visual overview is available for this article.

Genome-Wide Association Transethnic Meta-Analyses Identifies Novel Associations Regulating Coagulation Factor VIII and von Willebrand Factor Plasma Levels

BACKGROUND

Factor VIII (FVIII) and its carrier protein von Willebrand factor (VWF) are associated with risk of arterial and venous thrombosis and with hemorrhagic disorders. We aimed to identify and functionally test novel genetic associations regulating plasma FVIII and VWF.

METHODS

We meta-analyzed genome-wide association results from 46 354 individuals of European, African, East Asian, and Hispanic ancestry. All studies performed linear regression analysis using an additive genetic model and associated ≈35 million imputed variants with natural log-transformed phenotype levels. In vitro gene silencing in cultured endothelial cells was performed for candidate genes to provide additional evidence on association and function. Two-sample Mendelian randomization analyses were applied to test the causal role of FVIII and VWF plasma levels on the risk of arterial and venous thrombotic events.

RESULTS

We identified 13 novel genome-wide significant ( P≤2.5×10^-8) associations, 7 with FVIII levels ( FCHO2/TMEM171/TNPO1, HLA, SOX17/RP1, LINC00583/NFIB, RAB5C-KAT2A, RPL3/TAB1/SYNGR1, and ARSA) and 11 with VWF levels ( PDHB/PXK/KCTD6, SLC39A8, FCHO2/TMEM171/TNPO1, HLA, GIMAP7/GIMAP4, OR13C5/NIPSNAP, DAB2IP, C2CD4B, RAB5C-KAT2A, TAB1/SYNGR1, and ARSA), beyond 10 previously reported associations with these phenotypes. Functional validation provided further evidence of association for all loci on VWF except ARSA and DAB2IP. Mendelian randomization suggested causal effects of plasma FVIII activity levels on venous thrombosis and coronary artery disease risk and plasma VWF levels on ischemic stroke risk.

CONCLUSIONS

The meta-analysis identified 13 novel genetic loci regulating FVIII and VWF plasma levels, 10 of which we validated functionally. We provide some evidence for a causal role of these proteins in thrombotic events.

Associations of autozygosity with a broad range of human phenotypes

In many species, the offspring of related parents suffer reduced reproductive success, a phenomenon known as inbreeding depression. In humans, the importance of this effect has remained unclear, partly because reproduction between close relatives is both rare and frequently associated with confounding social factors. Here, using genomic inbreeding coefficients (FROH) for >1.4 million individuals, we show that FROH is significantly associated (p < 0.0005) with apparently deleterious changes in 32 out of 100 traits analysed. These changes are associated with runs of homozygosity (ROH), but not with common variant homozygosity, suggesting that genetic variants associated with inbreeding depression are predominantly rare. The effect on fertility is striking: FROH equivalent to the offspring of first cousins is associated with a 55% decrease [95% CI 44-66%] in the odds of having children. Finally, the effects of FROH are confirmed within full-sibling pairs, where the variation in FROH is independent of all environmental confounding.

Influence of ABO Locus on PFA-100 Collagen-ADP Closure Time Is Not Totally Dependent on the Von Willebrand Factor. Results of a GWAS on GAIT-2 Project Phenotypes

(1) Background: In a previous study, we found that two phenotypes related to platelet reactivity, measured with the PFA-100 system, were highly heritable. The aim of the present study was to identify genetic determinants that influence the variability of these phenotypes: closure time of collagen-ADP (Col-ADP) and of collagen-epinephrine (Col-Epi). (2) Methods: As part of the GAIT-2 (Genetic Analysis of Idiopathic Thrombophilia (2) Project, 935 individuals from 35 large Spanish families were studied. A genome-wide association study (GWAS) with ≈ 10 M single nucleotide polymorphisms (SNPs) was carried out with Col-ADP and Col-Epi phenotypes. (3) Results: The study yielded significant genetic signals that mapped to the ABO locus. After adjusting both phenotypes for the ABO genotype, these signals disappeared. After adjusting for von Willebrand factor (VWF) or for coagulation factor VIII (FVIII), the significant signals disappeared totally for Col-Epi phenotype but only partially for Col-ADP phenotype. (4) Conclusion: Our results suggest that the ABO locus exerts the main genetic influence on PFA-100 phenotypes. However, while the effect of the ABO locus on Col-Epi phenotype is mediated through VWF and/or FVIII, the effect of the ABO locus on Col-ADP phenotype is partly produced through VWF and/or FVIII, and partly through other mechanisms.

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.

Author Correction: Discovery and refinement of genetic loci associated with cardiometabolic risk using dense imputation maps

In the version of the article published, the surname of author Aaron Isaacs is misspelled as Issacs.

Genome Analyses of >200,000 Individuals Identify 58 Loci for Chronic Inflammation and Highlight Pathways that Link Inflammation and Complex Disorders

C-reactive protein (CRP) is a sensitive biomarker of chronic low-grade inflammation and is associated with multiple complex diseases. The genetic determinants of chronic inflammation remain largely unknown, and the causal role of CRP in several clinical outcomes is debated. We performed two genome-wide association studies (GWASs), on HapMap and 1000 Genomes imputed data, of circulating amounts of CRP by using data from 88 studies comprising 204,402 European individuals. Additionally, we performed in silico functional analyses and Mendelian randomization analyses with several clinical outcomes. The GWAS meta-analyses of CRP revealed 58 distinct genetic loci (p < 5 × 10-8). After adjustment for body mass index in the regression analysis, the associations at all except three loci remained. The lead variants at the distinct loci explained up to 7.0% of the variance in circulating amounts of CRP. We identified 66 gene sets that were organized in two substantially correlated clusters, one mainly composed of immune pathways and the other characterized by metabolic pathways in the liver. Mendelian randomization analyses revealed a causal protective effect of CRP on schizophrenia and a risk-increasing effect on bipolar disorder. Our findings provide further insights into the biology of inflammation and could lead to interventions for treating inflammation and its clinical consequences.

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.

Abstract 23: Identification of NCF4 as a Novel Regulator in Arterial Remodeling and Advanced Atherosclerosis

Arterial wall remodeling is a central multifactorial process in the development and progression of cardiovascular diseases. We employed an approach aimed at observing genetic variants associated with the progression of carotid intima-media thickness (cIMT) in order to identify novel pathways effecting vessel remodeling. This was achieved by conducting gene-centric analysis of 400,000 variants in 3,042 subjects with repeated cIMT measurements.

Rs16997464 on chr22 intergenic between neutrophil cytosolic factor-4 ( NCF4) and colony stimulating factor 2 receptor beta ( CSF2RB ) was associated with cIMT progression at array-wide significance (p <4.5x10 -7 ). The potential causative genes within this locus were investigated using a human vascular and non-vascular tissue biobank. Expression of 9 genes near rs16997464, were analyzed with the most significant association being with NCF4 in aortic adventitia. The effect of the variant on the function of the NCF4 gene product was further analyzed by comparing the oxidative burst capacity of neutrophils from subjects with different rs16997464 genotypes. We observed that neutrophils homozygous for the minor T allele, associated with slower cIMT progression, produced more extracellular ROS than neutrophils homozygous for the G allele, indicating a functional effect of rs16997464 on the NCF4 gene product p40 phox , a component of the NADPH oxidase 2 complex (NOX2).

In parallel, we investigated if the chr22 locus also influenced the cellular composition of the atherosclerotic plaque, by utilizing data from the Athero-Express Biobank. Here we found that the minor T allele associated with a higher smooth muscle cell (SMC) content in the plaque. Finally, using a partial ligation model in mice where ncf4 is mutated, resulting in a reduced but not absent NOX2-associated ROS formation, we observed a reduced neointima formation in the ncf4 -mutated strain compared with wild-type littermates.

Thus, this study identified rs16997464 in the NCF4-CSF2RB locus as a novel genetic determinant of cIMT progression, and provides evidence suggesting that NCF4 is involved in SMC proliferation and alteration of vessel wall pathophysiology.

Mapping of 79 loci for 83 plasma protein biomarkers in cardiovascular disease

Recent advances in highly multiplexed immunoassays have allowed systematic large-scale measurement of hundreds of plasma proteins in large cohort studies. In combination with genotyping, such studies offer the prospect to 1) identify mechanisms involved with regulation of protein expression in plasma, and 2) determine whether the plasma proteins are likely to be causally implicated in disease. We report here the results of genome-wide association (GWA) studies of 83 proteins considered relevant to cardiovascular disease (CVD), measured in 3,394 individuals with multiple CVD risk factors. We identified 79 genome-wide significant (p<5e-8) association signals, 55 of which replicated at P<0.0007 in separate validation studies (n = 2,639 individuals). Using automated text mining, manual curation, and network-based methods incorporating information on expression quantitative trait loci (eQTL), we propose plausible causal mechanisms for 25 trans-acting loci, including a potential post-translational regulation of stem cell factor by matrix metalloproteinase 9 and receptor-ligand pairs such as RANK-RANK ligand. Using public GWA study data, we further evaluate all 79 loci for their causal effect on coronary artery disease, and highlight several potentially causal associations. Overall, a majority of the plasma proteins studied showed evidence of regulation at the genetic level. Our results enable future studies of the causal architecture of human disease, which in turn should aid discovery of new drug targets.

Several proteins that circulate in blood have been linked to cardiovascular disease through the use of classic epidemiology and correlation studies. If individuals with higher risk of disease have higher levels of a protein, the protein may be associated with disease. However, this does not necessarily mean that the protein causes disease; it may merely be an innocent bystander or a consequence of the disease process. To establish whether a protein causes disease, a genetic approach, insensitive to reverse causation, can be used. Instead of correlating the levels of the protein itself, gene variants that regulate the protein levels are used in the analysis. This approach requires prior knowledge of which genetic variants are linked to individual proteins. Therefore we completed a map of how common genetic variants affect the blood concentration levels of 83 proteins that have been implicated in cardiovascular disease. By using this map of cause-to-effect findings, we gained insights into the regulation of a majority of the proteins under study and how they relate to risk of coronary artery disease. This study provides a map of genetic regulation of important cardiovascular plasma proteins, insights into their upstream regulatory environment, as well as novel leads for cardiovascular drug development.

Comparison of HapMap and 1000 Genomes Reference Panels in a Large-Scale Genome-Wide Association Study

An increasing number of genome-wide association (GWA) studies are now using the higher resolution 1000 Genomes Project reference panel (1000G) for imputation, with the expectation that 1000G imputation will lead to the discovery of additional associated loci when compared to HapMap imputation. In order to assess the improvement of 1000G over HapMap imputation in identifying associated loci, we compared the results of GWA studies of circulating fibrinogen based on the two reference panels. Using both HapMap and 1000G imputation we performed a meta-analysis of 22 studies comprising the same 91,953 individuals. We identified six additional signals using 1000G imputation, while 29 loci were associated using both HapMap and 1000G imputation. One locus identified using HapMap imputation was not significant using 1000G imputation. The genome-wide significance threshold of 5×10-8 is based on the number of independent statistical tests using HapMap imputation, and 1000G imputation may lead to further independent tests that should be corrected for. When using a stricter Bonferroni correction for the 1000G GWA study (P-value < 2.5×10-8), the number of loci significant only using HapMap imputation increased to 4 while the number of loci significant only using 1000G decreased to 5. In conclusion, 1000G imputation enabled the identification of 20% more loci than HapMap imputation, although the advantage of 1000G imputation became less clear when a stricter Bonferroni correction was used. More generally, our results provide insights that are applicable to the implementation of other dense reference panels that are under development.

Phenotypic Modulation of Smooth Muscle Cells in Atherosclerosis Is Associated With Downregulation of LMOD1, SYNPO2, PDLIM7, PLN, and SYNM

OBJECTIVE

Key augmented processes in atherosclerosis have been identified, whereas less is known about downregulated pathways. Here, we applied a systems biology approach to examine suppressed molecular signatures, with the hypothesis that they may provide insight into mechanisms contributing to plaque stability.

APPROACH AND RESULTS

Muscle contraction, muscle development, and actin cytoskeleton were the most downregulated pathways (false discovery rate=6.99e-21, 1.66e-6, 2.54e-10, respectively) in microarrays from human carotid plaques (n=177) versus healthy arteries (n=15). In addition to typical smooth muscle cell (SMC) markers, these pathways also encompassed cytoskeleton-related genes previously not associated with atherosclerosis. SYNPO2, SYNM, LMOD1, PDLIM7, and PLN expression positively correlated to typical SMC markers in plaques (Pearson r>0.6, P<0.0001) and in rat intimal hyperplasia (r>0.8, P<0.0001). By immunohistochemistry, the proteins were expressed in SMCs in normal vessels, but largely absent in human plaques and intimal hyperplasia. Subcellularly, most proteins localized to the cytoskeleton in cultured SMCs and were regulated by active enhancer histone modification H3K27ac by chromatin immunoprecipitation-sequencing. Functionally, the genes were downregulated by PDGFB (platelet-derived growth factor beta) and IFNg (interferron gamma), exposure to shear flow stress, and oxLDL (oxidized low-density lipoprotein) loading. Genetic variants in PDLIM7, PLN, and SYNPO2 loci associated with progression of carotid intima-media thickness in high-risk subjects without symptoms of cardiovascular disease (n=3378). By eQTL (expression quantitative trait locus), rs11746443 also associated with PDLIM7 expression in plaques. Mechanistically, silencing of PDLIM7 in vitro led to downregulation of SMC markers and disruption of the actin cytoskeleton, decreased cell spreading, and increased proliferation.

CONCLUSIONS

We identified a panel of genes that reflect the altered phenotype of SMCs in vascular disease and could be early sensitive markers of SMC dedifferentiation.

Analysis with the exome array identifies multiple new independent variants in lipid loci

It has been hypothesized that low frequency (1-5% minor allele frequency (MAF)) and rare (<1% MAF) variants with large effect sizes may contribute to the missing heritability in complex traits. Here, we report an association analysis of lipid traits (total cholesterol, LDL-cholesterol, HDL-cholesterol triglycerides) in up to 27 312 individuals with a comprehensive set of low frequency coding variants (ExomeChip), combined with conditional analysis in the known lipid loci. No new locus reached genome-wide significance. However, we found a new lead variant in 26 known lipid association regions of which 16 were >1000-fold more significant than the previous sentinel variant and not in close LD (six had MAF <5%). Furthermore, conditional analysis revealed multiple independent signals (ranging from 1 to 5) in a third of the 98 lipid loci tested, including rare variants. Addition of our novel associations resulted in between 1.5- and 2.5-fold increase in the proportion of heritability explained for the different lipid traits. Our findings suggest that rare coding variants contribute to the genetic architecture of lipid traits.

Abstract 173: Proprotein Convertase Subtilisin/Kexin Type 6 is a Key Protease in the Control of Smooth Muscle Cell Function in Vascular Disease

Proprotein convertases (PCSKs) process matrix metalloproteases (MMPs) and cytokines. Apart from PCSK9, the role of these enzymes in vascular disease is largely unknown. Previously, we demonstrated upregulation of PCSK6 in carotid atherosclerosis, primarily localized to smooth muscle cells (SMCs) and positively correlated to inflammation, extracellular matrix remodeling and cytokines. Here, we extended these findings to determine the role of PCSK6 in vascular development and disease. Increased expression of PCSK6 in vascular disease was validated by microarrays from two non-overlapping cohorts of carotid plaques vs. non-atherosclerotic arteries (n=50 patients and n=32 patients, p<0.0001), as well as abdominal (AAA, n=14, p<0.0001) and thoracic aortic aneurysms (TAA, n=244, p=0.012). By eQTL, variants in the PCSK6 gene were found to influence it’s expression in both plaques and aneurysms. Among these, rs6598465 also showed association with maximum progression of carotid intima-media thickness in high-risk coronary artery disease subjects (n=3388, p=0.037). By IHC, PCSK6 localized mainly to SMCs in the fibrous cap and neovessels in atherosclerotic, AAA and TAA tissues. In mouse-, rat-, and human intimal hyperplasia, PCSK6 was expressed in proliferating SMCs. By microarrays, after rat carotid balloon injury there was an early downregulation of PCSK6 followed by an upregulation in later phases during SMC activation, as well as positive correlation to PDGFB and IGF1 (Spearman r>0.7, p<0.0001) and to MMP2 and MMP14 (r>0.5, p<0.0001). In zebrafish embryos, PCSK6 localized to heart and vasculature and its ablation caused defective peripheral vascular patterning with cerebral and myocardial hemorrhage. PCSK6 -/- mice did not present an obvious vascular phenotype but showed reduced intimal hyperplasia compared to wild-type mice after carotid artery ligation (p=0.015). In vitro, PCSK6 overexpression markedly increased SMC migration upon PDGFBB stimulation (p<0.0001). The present study establishes PCSK6 as a key modulator of SMC function in vascular disease and demonstrates a functional link between PCSK6 expression and SMC migration in vascular remodeling.

Genome-wide association analysis of self-reported events in 6135 individuals and 252 827 controls identifies 8 loci associated with thrombosis

Thrombotic diseases are among the leading causes of morbidity and mortality in the world. To add insights into the genetic regulation of thrombotic disease, we conducted a genome-wide association study (GWAS) of 6135 self-reported blood clots events and 252 827 controls of European ancestry belonging to the 23andMe cohort of research participants. Eight loci exceeded genome-wide significance. Among the genome-wide significant results, our study replicated previously known venous thromboembolism (VTE) loci near the F5, FGA-FGG, F11, F2, PROCR and ABO genes, and the more recently discovered locus near SLC44A2 In addition, our study reports for the first time a genome-wide significant association between rs114209171, located upstream of the F8 structural gene, and thrombosis risk. Analyses of expression profiles and expression quantitative trait loci across different tissues suggested SLC44A2, ILF3 and AP1M2 as the three most plausible candidate genes for the chromosome 19 locus, our only genome-wide significant thrombosis-related locus that does not harbor likely coagulation-related genes. In addition, we present data showing that this locus also acts as a novel risk factor for stroke and coronary artery disease (CAD). In conclusion, our study reveals novel common genetic risk factors for VTE, stroke and CAD and provides evidence that self-reported data on blood clots used in a GWAS yield results that are comparable with those obtained using clinically diagnosed VTE. This observation opens up the potential for larger meta-analyses, which will enable elucidation of the genetics of thrombotic diseases, and serves as an example for the genetic study of other diseases.

A meta-analysis of 120 246 individuals identifies 18 new loci for fibrinogen concentration

Genome-wide association studies have previously identified 23 genetic loci associated with circulating fibrinogen concentration. These studies used HapMap imputation and did not examine the X-chromosome. 1000 Genomes imputation provides better coverage of uncommon variants, and includes indels. We conducted a genome-wide association analysis of 34 studies imputed to the 1000 Genomes Project reference panel and including ∼120 000 participants of European ancestry (95 806 participants with data on the X-chromosome). Approximately 10.7 million single-nucleotide polymorphisms and 1.2 million indels were examined. We identified 41 genome-wide significant fibrinogen loci; of which, 18 were newly identified. There were no genome-wide significant signals on the X-chromosome. The lead variants of five significant loci were indels. We further identified six additional independent signals, including three rare variants, at two previously characterized loci: FGB and IRF1. Together the 41 loci explain 3% of the variance in plasma fibrinogen concentration.

Sex-Specific Effects of Adiponectin on Carotid Intima-Media Thickness and Incident Cardiovascular Disease

BACKGROUND

Plasma adiponectin levels have previously been inversely associated with carotid intima-media thickness (IMT), a marker of subclinical atherosclerosis. In this study, we used a sex-stratified Mendelian randomization approach to investigate whether adiponectin has a causal protective influence on IMT.

METHODS AND RESULTS

Baseline plasma adiponectin concentration was tested for association with baseline IMT, IMT progression over 30 months, and occurrence of cardiovascular events within 3 years in 3430 participants (women, n=1777; men, n=1653) with high cardiovascular risk but no prevalent disease. Plasma adiponectin levels were inversely associated with baseline mean bifurcation IMT after adjustment for established risk factors (β=-0.018, P<0.001) in men but not in women (β=-0.006, P=0.185; P for interaction=0.061). Adiponectin levels were inversely associated with progression of mean common carotid IMT in men (β=-0.0022, P=0.047), whereas no association was seen in women (0.0007, P=0.475; P for interaction=0.018). Moreover, we observed that adiponectin levels were inversely associated with coronary events in women (hazard ratio 0.57, 95% CI 0.37 to 0.87) but not in men (hazard ratio 0.82, 95% CI 0.54 to 1.25). A gene score of adiponectin-raising alleles in 6 loci, reported recently in a large multi-ethnic meta-analysis, was inversely associated with baseline mean bifurcation IMT in men (β=-0.0008, P=0.004) but not in women (β=-0.0003, P=0.522; P for interaction=0.007).

CONCLUSIONS

This report provides some evidence for adiponectin protecting against atherosclerosis, with effects being confined to men; however, compared with established cardiovascular risk factors, the effect of plasma adiponectin was modest. Further investigation involving mechanistic studies is warranted.

Abstract 367: Pcsk6 Is a Key Protease Modulating Smooth Muscle Cell Activation in Vascular Remodeling and Plaque Vulnerability

Proprotein convertases (PCSKs) are conserved among species and involved in processing of MMPs and growth factors, but poorly characterised in atherosclerosis. Previously we demonstrated upregulation of PCSK6 in a large cohort of plaques from symptomatic vs. asymptomatic patients. This protease localized to smooth muscle cells (SMCs) and showed a positive correlation to markers of inflammation, extracellular matrix remodeling and cytokines in plaques. Here we aimed at elucidating its role in vascular development and disease.

In zebrafish embryos PCSK6 localized to heart and vasculature and its ablation caused defective peripheral vascular patterning with cerebral and myocardial hemorrhage. Increased expression of PCSK6 in vascular pathologies was validated by microarrays from carotid plaques vs. undiseased arteries (n=32 patients, p<0.0001), abdominal (AAA, n=14, p<0.0001) and thoracic aortic aneurysms (TAA, n=244, p=0.012). By immunohistochemistry, PCSK6 localized mainly to SMCs in the fibrous cap and neo-vessels in plaques, AAA and TAA tissues. Investigation of mouse carotid ligation, rat artery balloon injury and human restenosis tissues with pronounced intimal hyperplasia confirmed induction of PCSK6 in proliferating SMCs. By microarrays following progression of rat carotid injury, PCSK6 was downregulated in early phases defined by inflammatory pathways, while upregulated in later phases with SMCs activation and localized in SMCs. Expression of PCSK6 in this model was positively correlated to PDGFB and IGF1 (Spearman r>0.7, p<0.0001). PCSK6 overexpression in SMCs in vitro increased their migration in a wound-healing assay, especially upon PDGFBB stimulation. By eQTL analyses several polymorphisms in the PCSK6 gene were found to influence its expression in plaques and aneurysm tissue. Among these, rs6598465 showed association with maximum progression of arterial intima-media thickness in high-risk coronary artery disease subjects (n=3400, p=0.037).

We identified a functional link between elevated expression of PCSK6 and vascular remodeling characterized by SMC activation. The present study establishes PCSK6 as one of the key modulators of pathological processes in relation to plaque vulnerability.

Abstract 150: Identification of SYNPO2, SYNM, LMOD1, PDLIM7 and PLN as Novel Markers of Smooth Muscle Cells in Atherosclerosis

Smooth muscle cell (SMC) activation is a hallmark of vascular remodeling associated with downregulation of SMC specific genes and altered cell function. Our aim was to identify a novel set of genes silenced following SMC activation in vascular disease.

We compared microarrays from carotid plaques (n=177) vs. undiseased arteries and symptomatic vs. asymptomatic patients and found that ‘muscle contraction’, ‘muscle differentiation’ and ‘cytoskeleton based migration’ were the most significantly downregulated pathways (p=6.06e-23; 5.22e-11; 4.62e-11 respectively). Genes clustered in these categories were known SMC markers together with a novel set of genes functionally coupled to actin cytoskeleton: SYNPO2, SYNM, LMOD1, PDLIM7 and PLN. Transcription factors regulating the expression of these genes showed enrichment of SRF, MYOD1 and MYOGENIN, known to control muscle differentiation. Downregulation of these genes was validated in independent microarrays from carotid (n=21) and coronary plaques (n=38). Their expression was positively correlated to typical SMC markers both in human plaques (Pearson r>0.8, p<0.0001) and rat carotid intimal hyperplasia (r>0.9, p<0.0001). In microarrays from rat carotid balloon injury, these genes were downregulated in the early phases of tissue remodeling but gradually reappeared in the mature intima. By immunohistochemistry, the proteins localized to SMCs in normal human vessels, but were mostly absent in plaques and restenosis tissues. In vitro, SYNPO2 and PDLIM7 localized to actin filaments, SYNM to cortical cytoskeleton, PLN to nucleus while LMOD1 was rapidly lost in cultured human and primary rat SMCs. Polymorphisms in some of these genes associated with carotid intima-media thickness in high-risk coronary artery disease subjects (n=3400). By eQTL analyses, rs11746443 influenced PDLIM7 expression in plaques. Interestingly, this polymorphism constitutes the binding site for HEY1 transcription factor implicated in vascular development.

We identified a panel of novel SMC proteins, which are lost in vascular disease and reflect the altered phenotype of SMCs in vascular remodeling. Our results indicate that these could be early and sensitive markers of the SMCs dedifferentiation in atherosclerosis.

Predicting venous thrombosis in women using a combination of genetic markers and clinical risk factors

BACKGROUND

Family history of venous thromboembolism (VTE) has been suggested to be more useful in risk assessment than thrombophilia testing.

OBJECTIVES

We investigated established genetic susceptibility variants for association with VTE and evaluated a genetic risk score in isolation and combined with known trigger factors, including family history of VTE.

PATIENTS/METHOD

A total of 18 single nucleotide polymorphisms (SNPs) selected from the literature were genotyped in 2835 women participating in a Swedish nationwide case-control study (the ThromboEmbolism Hormone Study [TEHS]). Association with VTE was assessed by odds ratios (ORs) with 95% confidence interval (CI) using logistic regression. Clinical and genetic predictors that contributed significantly to the fit of the logistic regression model were included in the prediction models. SNP-SNP interactions were investigated and incorporated into the models if found significant. Risk scores were evaluated by calculating the area under the receiver-operating characteristics curve (AUC).

RESULTS

Seven SNPs (F5 rs6025, F2 rs1799963, ABO rs514659, FGG rs2066865, F11 rs2289252, PROC rs1799810 and KNG1 rs710446) with four SNP-SNP interactions contributed to the genetic risk score for VTE, with an AUC of 0.66 (95% CI, 0.64-0.68). After adding clinical risk factors, which included family history of VTE, the AUC reached 0.84 (95% CI, 0.82-0.85). The goodness of fit of the genetic and combined scores improved when significant SNP-SNP interaction terms were included.

CONCLUSION

Prediction of VTE in high-risk individuals was more accurate when a combination of clinical and genetic predictors with SNP-SNP interactions was included in a risk score.

Low-frequency and rare exome chip variants associate with fasting glucose and type 2 diabetes susceptibility

Fasting glucose and insulin are intermediate traits for type 2 diabetes. Here we explore the role of coding variation on these traits by analysis of variants on the HumanExome BeadChip in 60,564 non-diabetic individuals and in 16,491 T2D cases and 81,877 controls. We identify a novel association of a low-frequency nonsynonymous SNV in GLP1R (A316T; rs10305492; MAF=1.4%) with lower FG (β=-0.09±0.01 mmol l(-1), P=3.4 × 10(-12)), T2D risk (OR[95%CI]=0.86[0.76-0.96], P=0.010), early insulin secretion (β=-0.07±0.035 pmolinsulin mmolglucose(-1), P=0.048), but higher 2-h glucose (β=0.16±0.05 mmol l(-1), P=4.3 × 10(-4)). We identify a gene-based association with FG at G6PC2 (pSKAT=6.8 × 10(-6)) driven by four rare protein-coding SNVs (H177Y, Y207S, R283X and S324P). We identify rs651007 (MAF=20%) in the first intron of ABO at the putative promoter of an antisense lncRNA, associating with higher FG (β=0.02±0.004 mmol l(-1), P=1.3 × 10(-8)). Our approach identifies novel coding variant associations and extends the allelic spectrum of variation underlying diabetes-related quantitative traits and T2D susceptibility.