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

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.

Characterization of caffeine response regulatory variants in vascular endothelial cells

Genetic variants in gene regulatory sequences can modify gene expression and mediate the molecular response to environmental stimuli. In addition, genotype-environment interactions (GxE) contribute to complex traits such as cardiovascular disease. Caffeine is the most widely consumed stimulant and is known to produce a vascular response. To investigate GxE for caffeine, we treated vascular endothelial cells with caffeine and used a massively parallel reporter assay to measure allelic effects on gene regulation for over 43,000 genetic variants. We identified 665 variants with allelic effects on gene regulation and 6 variants that regulate the gene expression response to caffeine (GxE, false discovery rate [FDR] < 5%). When overlapping our GxE results with expression quantitative trait loci colocalized with coronary artery disease and hypertension, we dissected their regulatory mechanisms and showed a modulatory role for caffeine. Our results demonstrate that massively parallel reporter assay is a powerful approach to identify and molecularly characterize GxE in the specific context of caffeine consumption.

Genetic architecture of cardiac dynamic flow volumes

Cardiac blood flow is a critical determinant of human health. However, the definition of its genetic architecture is limited by the technical challenge of capturing dynamic flow volumes from cardiac imaging at scale. We present DeepFlow, a deep-learning system to extract cardiac flow and volumes from phase-contrast cardiac magnetic resonance imaging. A mixed-linear model applied to 37,653 individuals from the UK Biobank reveals genome-wide significant associations across cardiac dynamic flow volumes spanning from aortic forward velocity to aortic regurgitation fraction. Mendelian randomization reveals a causal role for aortic root size in aortic valve regurgitation. Among the most significant contributing variants, localizing genes (near ELN, PRDM6 and ADAMTS7) are implicated in connective tissue and blood pressure pathways. Here we show that DeepFlow cardiac flow phenotyping at scale, combined with genotyping data, reinforces the contribution of connective tissue genes, blood pressure and root size to aortic valve function.

Non-Random Enrichment of Single-Nucleotide Polymorphisms Associated with Clopidogrel Resistance within Risk Loci Linked to the Severity of Underlying Cardiovascular Diseases: The Role of Admixture

Cardiovascular disease (CVD) is one of the leading causes of death in Puerto Rico, where clopidogrel is commonly prescribed to prevent ischemic events. Genetic contributors to both a poor clopidogrel response and the severity of CVD have been identified mainly in Europeans. However, the non-random enrichment of single-nucleotide polymorphisms (SNPs) associated with clopidogrel resistance within risk loci linked to underlying CVDs, and the role of admixture, have yet to be tested. This study aimed to assess the possible interaction between genetic biomarkers linked to CVDs and those associated with clopidogrel resistance among admixed Caribbean Hispanics. We identified 50 SNPs significantly associated with CVDs in previous genome-wide association studies (GWASs). These SNPs were combined with another ten SNPs related to clopidogrel resistance in Caribbean Hispanics. We developed Python scripts to determine whether SNPs related to CVDs are in close proximity to those associated with the clopidogrel response. The average and individual local ancestry (LAI) within each locus were inferred, and 60 random SNPs with their corresponding LAIs were generated for enrichment estimation purposes. Our results showed no CVD-linked SNPs in close proximity to those associated with the clopidogrel response among Caribbean Hispanics. Consequently, no genetic loci with a dual predictive role for the risk of CVD severity and clopidogrel resistance were found in this population. Native American ancestry was the most enriched within the risk loci linked to CVDs in this population. The non-random enrichment of disease susceptibility loci with drug-response SNPs is a new frontier in Precision Medicine that needs further attention.

Intracellular tPA-PAI-1 interaction determines VLDL assembly in hepatocytes

Apolipoprotein B (apoB)-lipoproteins initiate and promote atherosclerotic cardiovascular disease. Plasma tissue plasminogen activator (tPA) activity is negatively associated with atherogenic apoB-lipoprotein cholesterol levels in humans, but the mechanisms are unknown. We found that tPA, partially through the lysine-binding site on its Kringle 2 domain, binds to the N terminus of apoB, blocking the interaction between apoB and microsomal triglyceride transfer protein (MTP) in hepatocytes, thereby reducing very-low-density lipoprotein (VLDL) assembly and plasma apoB-lipoprotein cholesterol levels. Plasminogen activator inhibitor 1 (PAI-1) sequesters tPA away from apoB and increases VLDL assembly. Humans with PAI-1 deficiency have smaller VLDL particles and lower plasma levels of apoB-lipoprotein cholesterol. These results suggest a mechanism that fine-tunes VLDL assembly by intracellular interactions among tPA, PAI-1, and apoB in hepatocytes.

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.

Inflammation compared to low-density lipoprotein cholesterol: two different causes of atherosclerotic cardiovascular disease

PURPOSE OF REVIEW

Inflammation is gaining attention as a target for prevention of atherosclerotic cardiovascular disease (ASCVD). The purpose of this review is to compare the evidence for inflammation with the evidence for low-density lipoprotein (LDL) cholesterol in ASCVD.

RECENT FINDINGS

Evidence from human genetic studies and randomized controlled trials implicate the inflammatory pathway from the inflammasome through interleukin (IL)-1 to IL-6 as a cause of ASCVD. Higher levels of IL-6 may lead to proportionally increased risk of ASCVD, and randomized controlled trials of IL-6 inhibitors are underway. The causal evidence for LDL cholesterol in ASCVD is overwhelming and recent important findings instead revolve around development of improved LDL cholesterol lowering therapy through RNA and DNA based therapeutics. Even though some lipid-lowering therapies lower IL-6, the IL-6 inflammatory pathway and LDL cholesterol are two separate causes of ASCVD.

SUMMARY

IL-6 mediated inflammation most likely causes ASCVD, in parallel with LDL cholesterol. However, fewer individuals in the general population are exposed to high IL-6 than high LDL cholesterol. For inflammation, future research should focus on improving efficacy and safety of anti-inflammatory therapy, and for LDL cholesterol, future research should focus on wider and more effective implementation of LDL cholesterol lowering therapy.

Association between PAI-1 Polymorphisms and Ischemic Stroke in a South Korean Case-Control Cohort

Stroke is the second leading cause of death in the world. Approximately 80% of strokes are ischemic in origin. Many risk factors have been linked to stroke, including an increased level of plasminogen activator inhibitor-1 (PAI-1). PAI-1 levels increase and remain elevated in blood during the acute phase of ischemic stroke, which can impair fibrinolytic activity, leading to coronary artery disease and arterial thrombotic disorders. Here, we present a case-control study of 574 stroke patients and 425 controls seen for routine health examination or treatment for nonspecific dizziness, nonorganic headache, or anxiety for positive family history of stroke at the Bundang Medical Center in South Korea. Polymorphisms in PAI-1 were identified by polymerase chain reaction/restriction fragment length polymorphism analysis using genomic DNA. Specifically, three variations (-675 4G>5G, 10692T>C, and 12068G>A) were linked to a higher overall prevalence of stroke as well as a higher prevalence of certain stroke subtypes. Haplotype analyses also revealed combinations of these variations (-844G>A, -675 4G>5G, 43G>A, 9785A>G, 10692T>C, 11053T>G, and 12068G>A) that were significantly associated with a higher prevalence of ischemic stroke. To the best of our knowledge, this is the first strong evidence that polymorphic sites in PAI-1 promoter and 3'-UTR regions are associated with higher ischemic stroke risk. Furthermore, the PAI-1 genotypes and haplotypes identified here have potential as clinical biomarkers of ischemic stroke and could improve the prognosis and future management of stroke patients.

Microvascular and macrovascular complications of type 2 diabetes mellitus: Exome wide association analyses

BACKGROUND

Type 2 diabetes mellitus (T2DM) is a chronic, metabolic disorder in which concomitant insulin resistance and β-cell impairment lead to hyperglycemia, influenced by genetic and environmental factors. T2DM is associated with long-term complications that have contributed to the burden of morbidity and mortality worldwide. The objective of this manuscript is to conduct an Exome-Wide Association Study (EWAS) on T2DM Emirati individuals to improve our understanding on diabetes-related complications to improve early diagnostic methods and treatment strategies.

METHODS

This cross-sectional study recruited 310 Emirati participants that were stratified according to their medically diagnosed diabetes-related complications: diabetic retinopathy, diabetic neuropathy, diabetic nephropathy, and cardiovascular complications. The Illumina's Infinium Exome-24 array was used and 39,840 SNPs remained for analysis after quality control.

FINDINGS

The analysis revealed the associations of various genes with each complication category: 1) diabetic retinopathy was associated to SHANK3 gene in locus 22q13.33 (SNP rs9616915; p=5.18 x10^-4), ZSCAN5A gene in locus 19q13.43 (SNP rs7252603; p=7.55 x10^-4), and DCP1B gene in locus 12p13.33 (SNPs rs715146, rs1044950, rs113147414, rs34730825; p=7.62 x10^-4); 2) diabetic neuropathy was associated to ADH4 gene in locus 4q23 (SNP rs4148883; p=1.23 x10^-4), SLC11A1 gene in locus 2q35 (SNP rs17235409; p=1.85 x10^-4), and MATN4 gene in locus 20q13.12 (SNP rs2072788; p=2.68 x10^-4); 3) diabetic nephropathy was associated to PPP1R3A gene in locus 7q31.1 (SNP rs1799999; p=1.91 x10^-4), ZNF136 gene in locus 19p13.2 (SNP rs140861589; p=2.80 x10^-4), and HSPA12B gene in locus 20p13 (SNP rs6076550; p=2.86 x10^-4); and 4) cardiovascular complications was associated to PCNT gene in locus 21q22.3 (SNPs rs7279204, rs6518289, rs2839227, rs2839223; p=2.18 x10^-4,3.04 x10^-4,4.51 x10^-4,5.22 x10^-4 respectively), SEPT14 gene in locus 7p11.2 (SNP rs146350220; p=2.77 x10^-4), and WDR73 gene in locus 15q25.2 (SNP rs72750868; p=4.47 x10^-4).

INTERPRETATION

We have identified susceptibility loci associated with each category of T2DM-related complications in the Emirati population. Given that only 16% of the markers from the Illumina's Infinium Exome chip passed quality control assessment, this demonstrates that multiple variants were, either, monomorphic in the Arab population or were not genotyped due to the use of a Euro-centric EWAS array that limits the possibility of including targeted ethnic-specific SNPs. Our results suggest the alarming possibility that lack of representation in reference panels could inhibit discovery of functionally important loci associated to T2DM complications. Further effort must be conducted to improve the representation of diverse populations in genotyping and sequencing studies.