Genome-Wide Investigation of Exogenous Female Hormones, Genetic Variation, and Venous Thromboembolism Risk.

BACKGROUND: Increased risk of venous thromboembolism (VTE) is a life-threatening side effect for users of oral contraceptives (OCs) or hormone therapy (HT). OBJECTIVES: To investigate the potential for genetic predisposition to VTE in OC or HT users, we conducted a gene-by-environment (GxE) case-only meta-analysis of genome-wide association studies (GWAS). METHODS: Use or non-use of OCs (7 studies) or HT (8 studies) at the time of the VTE event was determined by pharmacy records or self-report. A synergy index (SI) was modeled for each variant in each study and estimated supra-multiplicative GxE interaction. The SI parameters were first meta-analyzed across OC and HT studies, and subsequently meta-analyzed to obtain an overall estimate. The primary analysis was agnostic GWAS and interrogated all imputed genotypes using a p-value threshold of <5.0x10-8; secondary analyses were candidate-based. RESULTS: The VTE case-only OC meta-analysis included 2,895 OC users and 6,607 non-users; the case-only HT meta-analysis included 2,434 HT users and 12,793 non-users. In primary GWAS meta-analyses, no variant reached genome-wide significance, but the smallest p-value approached statistical significance: rs9386463 (p = 5.03x10-8). We tested associations for 138 candidate variants and identified 2 that exceeded statistical significance (0.05/138=3.62x10-4): F5 rs6025 (p = 1.87x10-5, SI = 1.29: previously observed) and F11 rs2036914 (p = 2.0x10-4, SI = 0.91; new observation). CONCLUSIONS: The candidate-variant approach to identify supra-multiplicative associations between genetic variation and OC-HT use identified a new association with common genetic variation in F11 while the agnostic interrogations did not yield new discoveries.

Methylation patterns associated with C-reactive protein in racially and ethnically diverse populations.

Systemic low-grade inflammation is a feature of chronic disease. C-reactive protein (CRP) is a common biomarker of inflammation and used as an indicator of disease risk; however, the role of inflammation in disease is not completely understood. Methylation is an epigenetic modification in the DNA which plays a pivotal role in gene expression. In this study we evaluated differential DNA methylation patterns associated with blood CRP level to elucidate biological pathways and genetic regulatory mechanisms to improve the understanding of chronic inflammation. The racially and ethnically diverse participants in this study were included as 50% White, 41% Black or African American, 7% Hispanic or Latino/a, and 2% Native Hawaiian, Asian American, American Indian, or Alaska Native (total n = 13,433) individuals. We replicated 113 CpG sites from 87 unique loci, of which five were novel (CADM3, NALCN, NLRC5, ZNF792, and cg03282312), across a discovery set of 1,150 CpG sites associated with CRP level (p < 1.2E-7). The downstream pathways affected by DNA methylation included the identification of IFI16 and IRF7 CpG-gene transcript pairs which contributed to the innate immune response gene enrichment pathway along with NLRC5, NOD2, and AIM2. Gene enrichment analysis also identified the nuclear factor-kappaB transcription pathway. Using two-sample Mendelian randomization (MR) we inferred methylation at three CpG sites as causal for CRP levels using both White and Black or African American MR instrument variables. Overall, we identified novel CpG sites and gene transcripts that could be valuable in understanding the specific cellular processes and pathogenic mechanisms involved in inflammation.

Time-to-Event Genome-Wide Association Study for Incident Cardiovascular Disease in People With Type 2 Diabetes.

OBJECTIVE: To identify genetic risk factors for incident cardiovascular disease (CVD) among people with type 2 diabetes (T2D). RESEARCH DESIGN AND METHODS: We conducted a multiancestry time-to-event genome-wide association study for incident CVD among people with T2D. We also tested 204 known coronary artery disease (CAD) variants for association with incident CVD. RESULTS: Among 49,230 participants with T2D, 8,956 had incident CVD events (event rate 18.2%). We identified three novel genetic loci for incident CVD: rs147138607 (near CACNA1E/ZNF648, hazard ratio [HR] 1.23, P = 3.6 × 10-9), rs77142250 (near HS3ST1, HR 1.89, P = 9.9 × 10-9), and rs335407 (near TFB1M/NOX3, HR 1.25, P = 1.5 × 10-8). Among 204 known CAD loci, 5 were associated with incident CVD in T2D (multiple comparison-adjusted P < 0.00024, 0.05/204). A standardized polygenic score of these 204 variants was associated with incident CVD with HR 1.14 (P = 1.0 × 10-16). CONCLUSIONS: The data point to novel and known genomic regions associated with incident CVD among individuals with T2D.

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 polygenic risk score of atrial fibrillation improves prediction of lifetime risk for heart failure.

AIMS: Heart failure (HF) has shared genetic architecture with its risk factors: atrial fibrillation (AF), body mass index (BMI), coronary heart disease (CHD), systolic blood pressure (SBP), and type 2 diabetes (T2D). We aim to assess the association and risk prediction performance of risk-factor polygenic risk scores (PRSs) for incident HF and its subtypes in bi-racial populations. METHODS AND RESULTS: Five PRSs were constructed for AF, BMI, CHD, SBP, and T2D in White participants of the Atherosclerosis Risk in Communities (ARIC) study. The associations between PRSs and incident HF and its subtypes were assessed using Cox models, and the risk prediction performance of PRSs was assessed using C statistics. Replication was performed in the ARIC study Black and Cardiovascular Health Study (CHS) White participants. In 8624 ARIC study Whites, 1922 (31% cumulative incidence) HF cases developed over 30 years of follow-up. PRSs of AF, BMI, and CHD were associated with incident HF (P < 0.001), where PRSAF showed the strongest association [hazard ratio (HR): 1.47, 95% confidence interval (CI): 1.41-1.53]. Only the addition of PRSAF to the ARIC study HF risk equation improved C statistics for 10 year risk prediction from 0.812 to 0.829 (∆C: 0.017, 95% CI: 0.009-0.026). The PRSAF was associated with both incident HF with reduced ejection fraction (HR: 1.43, 95% CI: 1.27-1.60) and incident HF with preserved ejection fraction (HR: 1.46, 95% CI: 1.33-1.62). The associations between PRSAF and incident HF and its subtypes, as well as the improved risk prediction, were replicated in the ARIC study Blacks and the CHS Whites (P < 0.050). Protein analyses revealed that N-terminal pro-brain natriuretic peptide and other 98 proteins were associated with PRSAF. CONCLUSIONS: The PRSAF was associated with incident HF and its subtypes and had significant incremental value over an established HF risk prediction equation.

Multi-trait analysis characterizes the genetics of thyroid function and identifies causal associations with clinical implications.

To date only a fraction of the genetic footprint of thyroid function has been clarified. We report a genome-wide association study meta-analysis of thyroid function in up to 271,040 individuals of European ancestry, including reference range thyrotropin (TSH), free thyroxine (FT4), free and total triiodothyronine (T3), proxies for metabolism (T3/FT4 ratio) as well as dichotomized high and low TSH levels. We revealed 259 independent significant associations for TSH (61% novel), 85 for FT4 (67% novel), and 62 novel signals for the T3 related traits. The loci explained 14.1%, 6.0%, 9.5% and 1.1% of the total variation in TSH, FT4, total T3 and free T3 concentrations, respectively. Genetic correlations indicate that TSH associated loci reflect the thyroid function determined by free T3, whereas the FT4 associations represent the thyroid hormone metabolism. Polygenic risk score and Mendelian randomization analyses showed the effects of genetically determined variation in thyroid function on various clinical outcomes, including cardiovascular risk factors and diseases, autoimmune diseases, and cancer. In conclusion, our results improve the understanding of thyroid hormone physiology and highlight the pleiotropic effects of thyroid function on various diseases.

Multi-ancestry polygenic risk scores for venous thromboembolism.

Venous thromboembolism (VTE) is a significant contributor to morbidity and mortality, with large disparities in incidence rates between Black and White Americans. Polygenic risk scores (PRSs) limited to variants discovered in genome-wide association studies in European-ancestry samples can identify European-ancestry individuals at high risk of VTE. However, there is limited evidence on whether high-dimensional PRS constructed using more sophisticated methods and more diverse training data can enhance the predictive ability and their utility across diverse populations. We developed PRSs for VTE using summary statistics from the International Network against Venous Thrombosis (INVENT) consortium GWAS meta-analyses of European- (71,771 cases and 1,059,740 controls) and African-ancestry samples (7,482 cases and 129,975 controls). We used LDpred2 and PRSCSx to construct ancestry-specific and multi-ancestry PRSs and evaluated their performance in an independent European- (6,261 cases and 88,238 controls) and African-ancestry sample (1,385 cases and 12,569 controls). Multi-ancestry PRSs with weights tuned in European- and African-ancestry samples, respectively, outperformed ancestry-specific PRSs in European- (PRSCSXEUR: AUC=0.61 (0.60, 0.61), PRSCSX_combinedEUR: AUC=0.61 (0.60, 0.62)) and African-ancestry test samples (PRSCSXAFR: AUC=0.58 (0.57, 0.6), PRSCSX_combined AFR: AUC=0.59 (0.57, 0.60)). The highest fifth percentile of the best-performing PRS was associated with 1.9-fold and 1.68-fold increased risk for VTE among European- and African-ancestry subjects, respectively, relative to those in the middle stratum. These findings suggest that the multi-ancestry PRS may be used to identify individuals at highest risk for VTE and provide guidance for the most effective treatment strategy across diverse populations.

Polygenic Risk Score Associates With Atherosclerotic Plaque Characteristics at Autopsy.

BACKGROUND: Polygenic risk scores (PRSs) for coronary artery disease (CAD) potentially improve cardiovascular risk prediction. However, their relationship with histopathologic features of CAD has never been examined systematically. METHODS: From 4327 subjects referred to CVPath by the State of Maryland Office Chief Medical Examiner for sudden death between 1994 and 2015, 2455 cases were randomly selected for genotyping. We generated PRS from 291 known CAD risk loci. Detailed histopathologic examination of the coronary arteries was performed in all subjects. The primary study outcome measurements were histopathologic plaque features determining severity of atherosclerosis, including %stenosis, calcification, thin-cap fibroatheromas, and thrombotic CAD. RESULTS: After exclusion of cases with insufficient DNA sample quality or with missing data, 954 cases (mean age, 48.8±14.7 years; 75.7% men) remained in the final study cohort. Subjects in the highest PRS quintile exhibited more severe atherosclerosis compared with subjects in the lowest quintile, with greater %stenosis (80.3%±27.0% versus 50.4%±38.7%; adjusted P<0.001) and a higher frequency of calcification (69.6% versus 35.8%; adjusted P=0.004) and thin-cap fibroatheroma (26.7% versus 9.5%; adjusted P=0.007). Even after adjustment for traditional CAD risk factors, subjects within the highest PRS quintile had higher odds of severe atherosclerosis (ie, ≥75% stenosis; adjusted odds ratio, 3.77 [95% CI, 2.10-6.78]; P<0.001) and plaque rupture (adjusted odds ratio, 4.05 [95% CI, 2.26-7.24]; P<0.001). Moreover, subjects within the highest quintile had higher odds of CAD-associated cause of death, especially among those aged ≤50 years (adjusted odds ratio, 4.08 [95% CI, 2.01-8.30]; P<0.001). No statistically significant associations were observed with plaque erosion after adjusting for covariates. CONCLUSIONS: This is the first autopsy study investigating associations between PRS and atherosclerosis severity at the histopathologic level in subjects with sudden death. Our pathological analysis suggests PRS correlates with plaque burden and features of advanced atherosclerosis and may be useful as a method for CAD risk stratification, especially in younger subjects.

Proteomic prediction of incident heart failure and its main subtypes.

AIM: To examine the ability of serum proteins in predicting future heart failure (HF) events, including HF with reduced or preserved ejection fraction (HFrEF or HFpEF), in relation to event time, and with or without considering established HF-associated clinical variables. METHODS AND RESULTS: In the prospective population-based Age, Gene/Environment Susceptibility Reykjavik Study (AGES-RS), 440 individuals developed HF after their first visit with a median follow-up of 5.45 years. Among them, 167 were diagnosed with HFrEF and 188 with HFpEF. A least absolute shrinkage and selection operator regression model with non-parametric bootstrap were used to select predictors from an analysis of 4782 serum proteins, and several pre-established clinical parameters linked to HF. A subset of 8-10 distinct or overlapping serum proteins predicted different future HF outcomes, and C-statistics were used to assess discrimination, revealing proteins combined with a C-index of 0.80 for all incident HF, 0.78 and 0.80 for incident HFpEF or HFrEF, respectively. In the AGES-RS, protein panels alone encompassed the risk contained in the clinical information and improved the performance characteristics of prediction models based on N-terminal pro-B-type natriuretic peptide and clinical risk factors. Finally, the protein predictors performed particularly well close to the time of an HF event, an outcome that was replicated in the Cardiovascular Health Study. CONCLUSION: A small number of circulating proteins accurately predicted future HF in the AGES-RS cohort of older adults, and they alone encompass the risk information found in a collection of clinical data. Incident HF events were predicted up to 8 years, with predictor performance significantly improving for events occurring less than 1 year ahead, a finding replicated in an external cohort study.

A statistical framework for powerful multi-trait rare variant analysis in large-scale whole-genome sequencing studies.

Large-scale whole-genome sequencing (WGS) studies have improved our understanding of the contributions of coding and noncoding rare variants to complex human traits. Leveraging association effect sizes across multiple traits in WGS rare variant association analysis can improve statistical power over single-trait analysis, and also detect pleiotropic genes and regions. Existing multi-trait methods have limited ability to perform rare variant analysis of large-scale WGS data. We propose MultiSTAAR, a statistical framework and computationally-scalable analytical pipeline for functionally-informed multi-trait rare variant analysis in large-scale WGS studies. MultiSTAAR accounts for relatedness, population structure and correlation among phenotypes by jointly analyzing multiple traits, and further empowers rare variant association analysis by incorporating multiple functional annotations. We applied MultiSTAAR to jointly analyze three lipid traits (low-density lipoprotein cholesterol, high-density lipoprotein cholesterol and triglycerides) in 61,861 multi-ethnic samples from the Trans-Omics for Precision Medicine (TOPMed) Program. We discovered new associations with lipid traits missed by single-trait analysis, including rare variants within an enhancer of NIPSNAP3A and an intergenic region on chromosome 1.

Type 2 Diabetes Modifies the Association of CAD Genomic Risk Variants With Subclinical Atherosclerosis.

BACKGROUND: Individuals with type 2 diabetes (T2D) have an increased risk of coronary artery disease (CAD), but questions remain about the underlying pathology. Identifying which CAD loci are modified by T2D in the development of subclinical atherosclerosis (coronary artery calcification [CAC], carotid intima-media thickness, or carotid plaque) may improve our understanding of the mechanisms leading to the increased CAD in T2D. METHODS: We compared the common and rare variant associations of known CAD loci from the literature on CAC, carotid intima-media thickness, and carotid plaque in up to 29 670 participants, including up to 24 157 normoglycemic controls and 5513 T2D cases leveraging whole-genome sequencing data from the Trans-Omics for Precision Medicine program. We included first-order T2D interaction terms in each model to determine whether CAD loci were modified by T2D. The genetic main and interaction effects were assessed using a joint test to determine whether a CAD variant, or gene-based rare variant set, was associated with the respective subclinical atherosclerosis measures and then further determined whether these loci had a significant interaction test. RESULTS: Using a Bonferroni-corrected significance threshold of P<1.6×10-4, we identified 3 genes (ATP1B1, ARVCF, and LIPG) associated with CAC and 2 genes (ABCG8 and EIF2B2) associated with carotid intima-media thickness and carotid plaque, respectively, through gene-based rare variant set analysis. Both ATP1B1 and ARVCF also had significantly different associations for CAC in T2D cases versus controls. No significant interaction tests were identified through the candidate single-variant analysis. CONCLUSIONS: These results highlight T2D as an important modifier of rare variant associations in CAD loci with CAC.

Identification of circulating proteins associated with general cognitive function among middle-aged and older adults.

Identifying circulating proteins associated with cognitive function may point to biomarkers and molecular process of cognitive impairment. Few studies have investigated the association between circulating proteins and cognitive function. We identify 246 protein measures quantified by the SomaScan assay as associated with cognitive function (p < 4.9E-5, n up to 7289). Of these, 45 were replicated using SomaScan data, and three were replicated using Olink data at Bonferroni-corrected significance. Enrichment analysis linked the proteins associated with general cognitive function to cell signaling pathways and synapse architecture. Mendelian randomization analysis implicated higher levels of NECTIN2, a protein mediating viral entry into neuronal cells, with higher Alzheimer's disease (AD) risk (p = 2.5E-26). Levels of 14 other protein measures were implicated as consequences of AD susceptibility (p < 2.0E-4). Proteins implicated as causes or consequences of AD susceptibility may provide new insight into the potential relationship between immunity and AD susceptibility as well as potential therapeutic targets.

The functional impact of rare variation across the regulatory cascade.

Each human genome has tens of thousands of rare genetic variants; however, identifying impactful rare variants remains a major challenge. We demonstrate how use of personal multi-omics can enable identification of impactful rare variants by using the Multi-Ethnic Study of Atherosclerosis, which included several hundred individuals, with whole-genome sequencing, transcriptomes, methylomes, and proteomes collected across two time points, 10 years apart. We evaluated each multi-omics phenotype's ability to separately and jointly inform functional rare variation. By combining expression and protein data, we observed rare stop variants 62 times and rare frameshift variants 216 times as frequently as controls, compared to 13-27 times as frequently for expression or protein effects alone. We extended a Bayesian hierarchical model, "Watershed," to prioritize specific rare variants underlying multi-omics signals across the regulatory cascade. With this approach, we identified rare variants that exhibited large effect sizes on multiple complex traits including height, schizophrenia, and Alzheimer's disease.

Mosaic chromosomal alterations in blood across ancestries using whole-genome sequencing.

Megabase-scale mosaic chromosomal alterations (mCAs) in blood are prognostic markers for a host of human diseases. Here, to gain a better understanding of mCA rates in genetically diverse populations, we analyzed whole-genome sequencing data from 67,390 individuals from the National Heart, Lung, and Blood Institute Trans-Omics for Precision Medicine program. We observed higher sensitivity with whole-genome sequencing data, compared with array-based data, in uncovering mCAs at low mutant cell fractions and found that individuals of European ancestry have the highest rates of autosomal mCAs and the lowest rates of chromosome X mCAs, compared with individuals of African or Hispanic ancestry. Although further studies in diverse populations will be needed to replicate our findings, we report three loci associated with loss of chromosome X, associations between autosomal mCAs and rare variants in DCPS, ADM17, PPP1R16B and TET2 and ancestry-specific variants in ATM and MPL with mCAs in cis.

WHOLE GENOME SEQUENCING ANALYSIS OF BODY MASS INDEX IDENTIFIES NOVEL AFRICAN ANCESTRY-SPECIFIC RISK ALLELE.

Obesity is a major public health crisis associated with high mortality rates. Previous genome-wide association studies (GWAS) investigating body mass index (BMI) have largely relied on imputed data from European individuals. This study leveraged whole-genome sequencing (WGS) data from 88,873 participants from the Trans-Omics for Precision Medicine (TOPMed) Program, of which 51% were of non-European population groups. We discovered 18 BMI-associated signals (P < 5 × 10-9). Notably, we identified and replicated a novel low frequency single nucleotide polymorphism (SNP) in MTMR3 that was common in individuals of African descent. Using a diverse study population, we further identified two novel secondary signals in known BMI loci and pinpointed two likely causal variants in the POC5 and DMD loci. Our work demonstrates the benefits of combining WGS and diverse cohorts in expanding current catalog of variants and genes confer risk for obesity, bringing us one step closer to personalized medicine.

Time-to-Event Genome-Wide Association Study for Incident Cardiovascular Disease in People with Type 2 Diabetes Mellitus.

BACKGROUND: Type 2 diabetes mellitus (T2D) confers a two- to three-fold increased risk of cardiovascular disease (CVD). However, the mechanisms underlying increased CVD risk among people with T2D are only partially understood. We hypothesized that a genetic association study among people with T2D at risk for developing incident cardiovascular complications could provide insights into molecular genetic aspects underlying CVD. METHODS: From 16 studies of the Cohorts for Heart & Aging Research in Genomic Epidemiology (CHARGE) Consortium, we conducted a multi-ancestry time-to-event genome-wide association study (GWAS) for incident CVD among people with T2D using Cox proportional hazards models. Incident CVD was defined based on a composite of coronary artery disease (CAD), stroke, and cardiovascular death that occurred at least one year after the diagnosis of T2D. Cohort-level estimated effect sizes were combined using inverse variance weighted fixed effects meta-analysis. We also tested 204 known CAD variants for association with incident CVD among patients with T2D. RESULTS: A total of 49,230 participants with T2D were included in the analyses (31,118 European ancestries and 18,112 non-European ancestries) which consisted of 8,956 incident CVD cases over a range of mean follow-up duration between 3.2 and 33.7 years (event rate 18.2%). We identified three novel, distinct genetic loci for incident CVD among individuals with T2D that reached the threshold for genome-wide significance (P<5.0×10-8): rs147138607 (intergenic variant between CACNA1E and ZNF648) with a hazard ratio (HR) 1.23, 95% confidence interval (CI) 1.15 - 1.32, P=3.6×10-9, rs11444867 (intergenic variant near HS3ST1) with HR 1.89, 95% CI 1.52 - 2.35, P=9.9×10-9, and rs335407 (intergenic variant between TFB1M and NOX3) HR 1.25, 95% CI 1.16 - 1.35, P=1.5×10-8. Among 204 known CAD loci, 32 were associated with incident CVD in people with T2D with P<0.05, and 5 were significant after Bonferroni correction (P<0.00024, 0.05/204). A polygenic score of these 204 variants was significantly associated with incident CVD with HR 1.14 (95% CI 1.12 - 1.16) per 1 standard deviation increase (P=1.0×10-16). CONCLUSIONS: The data point to novel and known genomic regions associated with incident CVD among individuals with T2D.

Genetic insights into resting heart rate and its role in cardiovascular disease.

Resting heart rate is associated with cardiovascular diseases and mortality in observational and Mendelian randomization studies. The aims of this study are to extend the number of resting heart rate associated genetic variants and to obtain further insights in resting heart rate biology and its clinical consequences. A genome-wide meta-analysis of 100 studies in up to 835,465 individuals reveals 493 independent genetic variants in 352 loci, including 68 genetic variants outside previously identified resting heart rate associated loci. We prioritize 670 genes and in silico annotations point to their enrichment in cardiomyocytes and provide insights in their ECG signature. Two-sample Mendelian randomization analyses indicate that higher genetically predicted resting heart rate increases risk of dilated cardiomyopathy, but decreases risk of developing atrial fibrillation, ischemic stroke, and cardio-embolic stroke. We do not find evidence for a linear or non-linear genetic association between resting heart rate and all-cause mortality in contrast to our previous Mendelian randomization study. Systematic alteration of key differences between the current and previous Mendelian randomization study indicates that the most likely cause of the discrepancy between these studies arises from false positive findings in previous one-sample MR analyses caused by weak-instrument bias at lower P-value thresholds. The results extend our understanding of resting heart rate biology and give additional insights in its role in cardiovascular disease development.

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.

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.

Polygenic Risk Score Associates with Atherosclerotic Plaque Characteristics at Autopsy.

Background: Polygenic risk scores (PRS) for coronary artery disease (CAD) potentially improve cardiovascular risk prediction. However, their relationship with histopathologic features of CAD has never been examined systematically. Methods: From 4,327 subjects referred to CVPath by the State of Maryland Office Chief Medical Examiner (OCME) for sudden death between 1994 and 2015, 2,455 cases were randomly selected for genotyping. We generated PRS from 291 known CAD risk loci. Detailed histopathologic examination of the coronary arteries was performed in all subjects. The primary study outcome measurements were histopathologic plaque features determining severity of atherosclerosis, including %stenosis, calcification, thin-cap fibroatheromas (TCFA), and thrombotic CAD. Results: After exclusion of cases with insufficient DNA sample quality or with missing data, 954 cases (mean age 48.8±14.7; 75.7% men) remained in the final study cohort. Subjects in the highest PRS quintile exhibited more severe atherosclerosis compared to subjects in the lowest quintile, with greater %stenosis (80.3%±27.0% vs. 50.4%±38.7%; adjusted p<0.001) and a higher frequency of calcification (69.6% vs. 35.8%; adjusted p=0.004) and TCFAs (26.7% vs. 9.5%; adjusted p=0.007). Even after adjustment for traditional CAD risk factors subjects within the highest PRS quintile had higher odds of severe atherosclerosis (i.e., ≥75% stenosis; adjusted OR 3.77; 95%CI 2.10-6.78; p<0.001) and plaque rupture (adjusted OR 4.05; 95%CI 2.26-7.24; p<0.001). Moreover, subjects within the highest quintile had higher odds of CAD-associated cause of death, especially among those aged 50 years and younger (adjusted OR 4.08; 95%CI 2.01-8.30; p<0.001). No associations were observed with plaque erosion. Conclusions: This is the first autopsy study investigating associations between PRS and atherosclerosis severity at the histopathologic level in subjects with sudden death. Our pathological analysis suggests PRS correlates with plaque burden and features of advanced atherosclerosis and may be useful as a method for CAD risk stratification, especially in younger subjects. Highlights: In this autopsy study including 954 subjects within the CVPath Sudden Death Registry, high PRS correlated with plaque burden and atherosclerosis severity.The PRS showed differential associations with plaque rupture and plaque erosion, suggesting different etiologies to these two causes of thrombotic CAD.PRS may be useful for risk stratification, particularly in the young. Further examination of individual risk loci and their association with plaque morphology may help understand molecular mechanisms of atherosclerosis, potentially revealing new therapy targets of CAD. Graphic Abstract: A polygenic risk score, generated from 291 known CAD risk loci, was assessed in 954 subjects within the CVPath Sudden Death Registry. Histopathologic examination of the coronary arteries was performed in all subjects. Subjects in the highest PRS quintile exhibited more severe atherosclerosis as compared to subjects in the lowest quintile, with a greater plaque burden, more calcification, and a higher frequency of plaque rupture.