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.

Genome-wide analysis in over 1 million individuals of European ancestry yields improved polygenic risk scores for blood pressure traits.

Hypertension affects more than one billion people worldwide. Here we identify 113 novel loci, reporting a total of 2,103 independent genetic signals (P < 5 × 10-8) from the largest single-stage blood pressure (BP) genome-wide association study to date (n = 1,028,980 European individuals). These associations explain more than 60% of single nucleotide polymorphism-based BP heritability. Comparing top versus bottom deciles of polygenic risk scores (PRSs) reveals clinically meaningful differences in BP (16.9 mmHg systolic BP, 95% CI, 15.5-18.2 mmHg, P = 2.22 × 10-126) and more than a sevenfold higher odds of hypertension risk (odds ratio, 7.33; 95% CI, 5.54-9.70; P = 4.13 × 10-44) in an independent dataset. Adding PRS into hypertension-prediction models increased the area under the receiver operating characteristic curve (AUROC) from 0.791 (95% CI, 0.781-0.801) to 0.826 (95% CI, 0.817-0.836, ∆AUROC, 0.035, P = 1.98 × 10-34). We compare the 2,103 loci results in non-European ancestries and show significant PRS associations in a large African-American sample. Secondary analyses implicate 500 genes previously unreported for BP. Our study highlights the role of increasingly large genomic studies for precision health research.

An approach to identify gene-environment interactions and reveal new biological insight in complex traits.

There is a long-standing debate about the magnitude of the contribution of gene-environment interactions to phenotypic variations of complex traits owing to the low statistical power and few reported interactions to date. To address this issue, the Gene-Lifestyle Interactions Working Group within the Cohorts for Heart and Aging Research in Genetic Epidemiology Consortium has been spearheading efforts to investigate G × E in large and diverse samples through meta-analysis. Here, we present a powerful new approach to screen for interactions across the genome, an approach that shares substantial similarity to the Mendelian randomization framework. We identify and confirm 5 loci (6 independent signals) interacted with either cigarette smoking or alcohol consumption for serum lipids, and empirically demonstrate that interaction and mediation are the major contributors to genetic effect size heterogeneity across populations. The estimated lower bound of the interaction and environmentally mediated heritability is significant (P < 0.02) for low-density lipoprotein cholesterol and triglycerides in Cross-Population data. Our study improves the understanding of the genetic architecture and environmental contributions to complex traits.

Genome-Wide Interaction Analysis With DASH Diet Score Identified Novel Loci for Systolic Blood Pressure.

BACKGROUND: The Dietary Approaches to Stop Hypertension (DASH) diet score lowers blood pressure (BP). We examined interactions between genotype and the DASH diet score in relation to systolic BP. METHODS: We analyzed up to 9 420 585 single nucleotide polymorphisms in up to 127 282 individuals of 6 population groups (91% of European population) from the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium (n=35 660) and UK Biobank (n=91 622) and performed European population-specific and cross-population meta-analyses. RESULTS: We identified 3 loci in European-specific analyses and an additional 4 loci in cross-population analyses at Pinteraction<5e-8. We observed a consistent interaction between rs117878928 at 15q25.1 (minor allele frequency, 0.03) and the DASH diet score (Pinteraction=4e-8; P for heterogeneity, 0.35) in European population, where the interaction effect size was 0.42±0.09 mm Hg (Pinteraction=9.4e-7) and 0.20±0.06 mm Hg (Pinteraction=0.001) in Cohorts for Heart and Aging Research in Genomic Epidemiology and the UK Biobank, respectively. The 1 Mb region surrounding rs117878928 was enriched with cis-expression quantitative trait loci (eQTL) variants (P=4e-273) and cis-DNA methylation quantitative trait loci variants (P=1e-300). Although the closest gene for rs117878928 is MTHFS, the highest narrow sense heritability accounted by single nucleotide polymorphisms potentially interacting with the DASH diet score in this locus was for gene ST20 at 15q25.1. CONCLUSIONS: We demonstrated gene-DASH diet score interaction effects on systolic BP in several loci. Studies with larger diverse populations are needed to validate our findings.

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.

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.

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.

Genome-Wide Interaction Analysis with DASH Diet Score Identified Novel Loci for Systolic Blood Pressure.

Objective: We examined interactions between genotype and a Dietary Approaches to Stop Hypertension (DASH) diet score in relation to systolic blood pressure (SBP). Methods: We analyzed up to 9,420,585 biallelic imputed single nucleotide polymorphisms (SNPs) in up to 127,282 individuals of six population groups (91% of European population) from the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium (CHARGE; n=35,660) and UK Biobank (n=91,622) and performed European population-specific and cross-population meta-analyses. Results: We identified three loci in European-specific analyses and an additional four loci in cross-population analyses at P for interaction < 5e-8. We observed a consistent interaction between rs117878928 at 15q25.1 (minor allele frequency = 0.03) and the DASH diet score (P for interaction = 4e-8; P for heterogeneity = 0.35) in European population, where the interaction effect size was 0.42±0.09 mm Hg (P for interaction = 9.4e-7) and 0.20±0.06 mm Hg (P for interaction = 0.001) in CHARGE and the UK Biobank, respectively. The 1 Mb region surrounding rs117878928 was enriched with cis-expression quantitative trait loci (eQTL) variants (P = 4e-273) and cis-DNA methylation quantitative trait loci (mQTL) variants (P = 1e-300). While the closest gene for rs117878928 is MTHFS, the highest narrow sense heritability accounted by SNPs potentially interacting with the DASH diet score in this locus was for gene ST20 at 15q25.1. Conclusion: We demonstrated gene-DASH diet score interaction effects on SBP in several loci. Studies with larger diverse populations are needed to validate our findings.

A new Approach to Identify Gene-Environment Interactions and Reveal New Biological Insight in Complex traits.

There is a long-standing debate about the magnitude of the contribution of gene-environment interactions to phenotypic variations of complex traits owing to the low statistical power and few reported interactions to date. To address this issue, the CHARGE Gene-Lifestyle Interactions Working Group has been spearheading efforts to investigate G×E in large and diverse samples through meta-analysis. Here, we present a powerful new approach to screen for interactions across the genome, an approach that shares substantial similarity to the Mendelian randomization framework. We identified and confirmed 5 loci (6 independent signals) interacting with either cigarette smoking or alcohol consumption for serum lipids, and empirically demonstrated that interaction and mediation are the major contributors to genetic effect size heterogeneity across populations. The estimated lower bound of the interaction and environmentally mediated contribution ranges from 1.76% to 14.05% of SNP heritability of serum lipids in Cross-Population data. Our study improves the understanding of the genetic architecture and environmental contributions to complex traits.

Multi-tissue epigenetic analysis identifies distinct associations underlying insulin resistance and Alzheimer's disease at CPT1A locus.

BACKGROUND: Insulin resistance (IR) is a major risk factor for Alzheimer's disease (AD) dementia. The mechanisms by which IR predisposes to AD are not well-understood. Epigenetic studies may help identify molecular signatures of IR associated with AD, thus improving our understanding of the biological and regulatory mechanisms linking IR and AD. METHODS: We conducted an epigenome-wide association study of IR, quantified using the homeostatic model assessment of IR (HOMA-IR) and adjusted for body mass index, in 3,167 participants from the Framingham Heart Study (FHS) without type 2 diabetes at the time of blood draw used for methylation measurement. We identified DNA methylation markers associated with IR at the genome-wide level accounting for multiple testing (P < 1.1 × 10-7) and evaluated their association with neurological traits in participants from the FHS (N = 3040) and the Religious Orders Study/Memory and Aging Project (ROSMAP, N = 707). DNA methylation profiles were measured in blood (FHS) or dorsolateral prefrontal cortex (ROSMAP) using the Illumina HumanMethylation450 BeadChip. Linear regressions (ROSMAP) or mixed-effects models accounting for familial relatedness (FHS) adjusted for age, sex, cohort, self-reported race, batch, and cell type proportions were used to assess associations between DNA methylation and neurological traits accounting for multiple testing. RESULTS: We confirmed the strong association of blood DNA methylation with IR at three loci (cg17901584-DHCR24, cg17058475-CPT1A, cg00574958-CPT1A, and cg06500161-ABCG1). In FHS, higher levels of blood DNA methylation at cg00574958 and cg17058475 were both associated with lower IR (P = 2.4 × 10-11 and P = 9.0 × 10-8), larger total brain volumes (P = 0.03 and P = 9.7 × 10-4), and smaller log lateral ventricular volumes (P = 0.07 and P = 0.03). In ROSMAP, higher levels of brain DNA methylation at the same two CPT1A markers were associated with greater risk of cognitive impairment (P = 0.005 and P = 0.02) and higher AD-related indices (CERAD score: P = 5 × 10-4 and 0.001; Braak stage: P = 0.004 and P = 0.01). CONCLUSIONS: Our results suggest potentially distinct epigenetic regulatory mechanisms between peripheral blood and dorsolateral prefrontal cortex tissues underlying IR and AD at CPT1A locus.

Rare variants in long non-coding RNAs are associated with blood lipid levels in the TOPMed whole-genome sequencing study.

Long non-coding RNAs (lncRNAs) are known to perform important regulatory functions in lipid metabolism. Large-scale whole-genome sequencing (WGS) studies and new statistical methods for variant set tests now provide an opportunity to assess more associations between rare variants in lncRNA genes and complex traits across the genome. In this study, we used high-coverage WGS from 66,329 participants of diverse ancestries with measurement of blood lipids and lipoproteins (LDL-C, HDL-C, TC, and TG) in the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) program to investigate the role of lncRNAs in lipid variability. We aggregated rare variants for 165,375 lncRNA genes based on their genomic locations and conducted rare-variant aggregate association tests using the STAAR (variant-set test for association using annotation information) framework. We performed STAAR conditional analysis adjusting for common variants in known lipid GWAS loci and rare-coding variants in nearby protein-coding genes. Our analyses revealed 83 rare lncRNA variant sets significantly associated with blood lipid levels, all of which were located in known lipid GWAS loci (in a ±500-kb window of a Global Lipids Genetics Consortium index variant). Notably, 61 out of 83 signals (73%) were conditionally independent of common regulatory variation and rare protein-coding variation at the same loci. We replicated 34 out of 61 (56%) conditionally independent associations using the independent UK Biobank WGS data. Our results expand the genetic architecture of blood lipids to rare variants in lncRNAs.

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.

Carriers of rare damaging CCR2 genetic variants are at lower risk of atherosclerotic disease.

Background: The CCL2/CCR2 axis governs monocyte trafficking and recruitment to atherosclerotic lesions. Human genetic analyses and population-based studies support an association between circulating CCL2 levels and atherosclerosis. Still, it remains unknown whether pharmacological targeting of CCR2, the main CCL2 receptor, would provide protection against human atherosclerotic disease. Methods: In whole-exome sequencing data from 454,775 UK Biobank participants (40-69 years), we identified predicted loss-of-function (LoF) or damaging missense (REVEL score >0.5) variants within the CCR2 gene. We prioritized variants associated with lower monocyte count (p<0.05) and tested associations with vascular risk factors and risk of atherosclerotic disease over a mean follow-up of 14 years. The results were replicated in a pooled cohort of three independent datasets (TOPMed, deCODE and Penn Medicine BioBank; total n=441,445) and the effect of the most frequent damaging variant was experimentally validated. Results: A total of 45 predicted LoF or damaging missense variants were identified in the CCR2 gene, 4 of which were also significantly associated with lower monocyte count, but not with other white blood cell counts. Heterozygous carriers of these variants were at a lower risk of a combined atherosclerosis outcome, showed a lower burden of atherosclerosis across four vascular beds, and were at a lower lifetime risk of coronary artery disease and myocardial infarction. There was no evidence of association with vascular risk factors including LDL-cholesterol, blood pressure, glycemic status, or C-reactive protein. Using a cAMP assay, we found that cells transfected with the most frequent CCR2 damaging variant (3:46358273:T:A, M249K, 547 carriers, frequency: 0.14%) show a decrease in signaling in response to CCL2. The associations of the M249K variant with myocardial infarction were consistent across cohorts (ORUKB: 0.62 95%CI: 0.39-0.96; ORexternal: 0.64 95%CI: 0.34-1.19; ORpooled: 0.64 95%CI: 0.450.90). In a phenome-wide association study, we found no evidence for higher risk of common infections or mortality among carriers of damaging CCR2 variants. Conclusions: Heterozygous carriers of damaging CCR2 variants have a lower burden of atherosclerosis and lower lifetime risk of myocardial infarction. In conjunction with previous evidence from experimental and epidemiological studies, our findings highlight the translational potential of CCR2-targeting as an atheroprotective approach.

Rare variants in long non-coding RNAs are associated with blood lipid levels in the TOPMed Whole Genome Sequencing Study.

Long non-coding RNAs (lncRNAs) are known to perform important regulatory functions. Large-scale whole genome sequencing (WGS) studies and new statistical methods for variant set tests now provide an opportunity to assess the associations between rare variants in lncRNA genes and complex traits across the genome. In this study, we used high-coverage WGS from 66,329 participants of diverse ancestries with blood lipid levels (LDL-C, HDL-C, TC, and TG) in the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) program to investigate the role of lncRNAs in lipid variability. We aggregated rare variants for 165,375 lncRNA genes based on their genomic locations and conducted rare variant aggregate association tests using the STAAR (variant-Set Test for Association using Annotation infoRmation) framework. We performed STAAR conditional analysis adjusting for common variants in known lipid GWAS loci and rare coding variants in nearby protein coding genes. Our analyses revealed 83 rare lncRNA variant sets significantly associated with blood lipid levels, all of which were located in known lipid GWAS loci (in a ±500 kb window of a Global Lipids Genetics Consortium index variant). Notably, 61 out of 83 signals (73%) were conditionally independent of common regulatory variations and rare protein coding variations at the same loci. We replicated 34 out of 61 (56%) conditionally independent associations using the independent UK Biobank WGS data. Our results expand the genetic architecture of blood lipids to rare variants in lncRNA, implicating new therapeutic opportunities.

Evaluating the use of blood pressure polygenic risk scores across race/ethnic background groups.

We assess performance and limitations of polygenic risk scores (PRSs) for multiple blood pressure (BP) phenotypes in diverse population groups. We compare "clumping-and-thresholding" (PRSice2) and LD-based (LDPred2) methods to construct PRSs from each of multiple GWAS, as well as multi-PRS approaches that sum PRSs with and without weights, including PRS-CSx. We use datasets from the MGB Biobank, TOPMed study, UK biobank, and from All of Us to train, assess, and validate PRSs in groups defined by self-reported race/ethnic background (Asian, Black, Hispanic/Latino, and White). For both SBP and DBP, the PRS-CSx based PRS, constructed as a weighted sum of PRSs developed from multiple independent GWAS, perform best across all race/ethnic backgrounds. Stratified analysis in All of Us shows that PRSs are better predictive of BP in females compared to males, individuals without obesity, and middle-aged (40-60 years) compared to older and younger individuals.

Coronary heart disease and ischemic stroke polygenic risk scores and atherosclerotic cardiovascular disease in a diverse, population-based cohort study.

The predictive ability of coronary heart disease (CHD) and ischemic stroke (IS) polygenic risk scores (PRS) have been evaluated individually, but whether they predict the combined outcome of atherosclerotic cardiovascular disease (ASCVD) remains insufficiently researched. It is also unclear whether associations of the CHD and IS PRS with ASCVD are independent of subclinical atherosclerosis measures. 7,286 White and 2,016 Black participants from the population-based Atherosclerosis Risk in Communities study who were free of cardiovascular disease and type 2 diabetes at baseline were included. We computed previously validated CHD and IS PRS consisting of 1,745,179 and 3,225,583 genetic variants, respectively. Cox proportional hazards models were used to test the association between each PRS and ASCVD, adjusting for traditional risk factors, ankle-brachial index, carotid intima media thickness, and carotid plaque. The hazard ratios (HR) for the CHD and IS PRS were significant with HR of 1.50 (95% CI: 1.36-1.66) and 1.31 (95% CI: 1.18-1.45) respectively for the risk of incident ASCVD per standard deviation increase in CHD and IS PRS among White participants after adjusting for traditional risk factors. The HR for the CHD PRS was not significant with an HR of 0.95 (95% CI: 0.79-1.13) for the risk of incident ASCVD in Black participants. The HR for the IS PRS was significant with an HR of 1.26 (95%CI: 1.05-1.51) for the risk of incident ASCVD in Black participants. The association of the CHD and IS PRS with ASCVD was not attenuated in White participants after adjustment for ankle-brachial index, carotid intima media thickness, and carotid plaque. The CHD and IS PRS do not cross-predict well, and predict better the outcome for which they were created than the composite ASCVD outcome. Thus, the use of the composite outcome of ASCVD may not be ideal for genetic risk prediction.

Whole-Genome Sequencing Analysis of Human Metabolome in Multi-Ethnic Populations.

Circulating metabolite levels may reflect the state of the human organism in health and disease, however, the genetic architecture of metabolites is not fully understood. We have performed a whole-genome sequencing association analysis of both common and rare variants in up to 11,840 multi-ethnic participants from five studies with up to 1666 circulating metabolites. We have discovered 1985 novel variant-metabolite associations, and validated 761 locus-metabolite associations reported previously. Seventy-nine novel variant-metabolite associations have been replicated, including three genetic loci located on the X chromosome that have demonstrated its involvement in metabolic regulation. Gene-based analysis have provided further support for seven metabolite-replicated loci pairs and their biologically plausible genes. Among those novel replicated variant-metabolite pairs, follow-up analyses have revealed that 26 metabolites have colocalized with 21 tissues, seven metabolite-disease outcome associations have been putatively causal, and 7 metabolites might be regulated by plasma protein levels. Our results have depicted the genetic contribution to circulating metabolite levels, providing additional insights into understanding human disease.

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.