A functional genomic framework to elucidate novel causal metabolic dysfunction-associated fatty liver disease genes.

BACKGROUND AIMS: Metabolic dysfunction-associated fatty liver disease (MASLD) is the most prevalent chronic liver pathology in western countries, with serious public health consequences. Efforts to identify causal genes for MASLD have been hampered by the relative paucity of human data from gold-standard magnetic resonance quantification of hepatic fat. To overcome insufficient sample size, genome-wide association studies using MASLD surrogate phenotypes have been used, but only a small number of loci have been identified to date. In this study, we combined GWAS of MASLD composite surrogate phenotypes with genetic colocalization studies followed by functional in vitro screens to identify bona fide causal genes for MASLD. APPROACH RESULTS: We used the UK Biobank to explore the associations of our novel MASLD score, and genetic colocalization to prioritize putative causal genes for in vitro validation. We created a functional genomic framework to study MASLD genes in vitro using CRISPRi. Our data identify VKORC1, TNKS, LYPLAL1 and GPAM as regulators of lipid accumulation in hepatocytes and suggest the involvement of VKORC1 in the lipid storage related to the development of MASLD. CONCLUSIONS: Complementary genetic and genomic approaches are useful for the identification of MASLD genes. Our data supports VKORC1 as a bona fide MASLD gene. We have established a functional genomic framework to study at scale putative novel MASLD genes from human genetic association studies.

Genetically predicted lipoprotein(a) associates with coronary artery plaque severity independent of low-density lipoprotein cholesterol.

AIMS: Elevated Lipoprotein(a) [Lp(a)] is a causal risk factor for atherosclerotic cardiovascular disease, but the mechanisms of risk are debated. Studies have found inconsistent associations between Lp(a) and measurements of atherosclerosis. We aimed to assess the relationship between Lp(a), low-density lipoprotein cholesterol (LDL-C) and coronary artery plaque severity. METHODS: The study population consisted of participants of the Million Veteran Program who have undergone an invasive angiogram. The primary exposure was genetically predicted Lp(a), estimated by a polygenic score. Genetically predicted LDL-C was also assessed for comparison. The primary outcome was coronary artery plaque severity, categorized as normal, non-obstructive disease, 1-vessel disease, 2-vessel disease, and 3-vessel or left main disease. RESULTS: Among 18,927 adults of genetically inferred European ancestry and 4,039 adults of genetically inferred African ancestry, we observed consistent associations between genetically predicted Lp(a) and obstructive coronary plaque, with effect sizes trending upward for increasingly severe categories of disease. Associations were independent of risk factors, clinically measured LDL-C and genetically predicted LDL-C. However, we did not find strong or consistent evidence for an association between genetically predicted Lp(a) and risk for non-obstructive plaque. CONCLUSIONS: Genetically predicted Lp(a) is positively associated with coronary plaque severity independent of LDL-C, consistent with Lp(a) promoting atherogenesis. However, the effects of Lp(a) may be greater for progression of plaque to obstructive disease than for the initial development of non-obstructive plaque. A limitation of this study is that Lp(a) was estimated using genetic markers and could not be directly assayed, nor could apo(a) isoform size.

This study assessed the association between genetic propensity towards higher lipoprotein(a) [Lp(a)] in the blood and the severity of coronary artery plaque seen on clinical angiograms, independent of other factors, including low-density lipoprotein cholesterol (LDL-C). The study was conducted in a large U.S. population using data from the Million Veteran Program. Genetically predicted high Lp(a) was associated with obstructive coronary plaque, but it was not associated with non-obstructive coronary plaque. This association was independent of LDL-C, and the association was greater for more severe forms of disease.The mechanisms of association between Lp(a) and cardiovascular events are debated. Prior studies have shown that Lp(a) does not associate with early markers of atherosclerosis. Our analyses support the idea that Lp(a) plays less of a role in early plaque initiation but plays a significant role in the progression of plaque towards more severe disease, independent of LDL-C.

Mendelian randomization analyses clarify the effects of height on cardiovascular diseases.

An inverse correlation between stature and risk of coronary artery disease (CAD) has been observed in several epidemiologic studies, and recent Mendelian randomization (MR) experiments have suggested causal association. However, the extent to which the effect estimated by MR can be explained by cardiovascular, anthropometric, lung function, and lifestyle-related risk factors is unclear, with a recent report suggesting that lung function traits could fully explain the height-CAD effect. To clarify this relationship, we utilized a well-powered set of genetic instruments for human stature, comprising >1,800 genetic variants for height and CAD. In univariable analysis, we confirmed that a one standard deviation decrease in height (~6.5 cm) was associated with a 12.0% increase in the risk of CAD, consistent with previous reports. In multivariable analysis accounting for effects from up to 12 established risk factors, we observed a >3-fold attenuation in the causal effect of height on CAD susceptibility (3.7%, p = 0.02). However, multivariable analyses demonstrated independent effects of height on other cardiovascular traits beyond CAD, consistent with epidemiologic associations and univariable MR experiments. In contrast with published reports, we observed minimal effects of lung function traits on CAD risk in our analyses, indicating that these traits are unlikely to explain the residual association between height and CAD risk. In sum, these results suggest the impact of height on CAD risk beyond previously established cardiovascular risk factors is minimal and not explained by lung function measures.

Increased BMI associated with decreased breastfeeding initiation in Million Veteran Program participants.

Background: Breastfeeding has been associated with maternal and infant health benefits but has been inversely associated with body mass index (BMI) prepartum. Breastfeeding and BMI are both linked to socioeconomic factors. Methods: Data from parous female participants with available breastfeeding information from the Million Veteran Program cohort was included. BMI at enrollment and earliest BMI available were extracted, and polygenic scores (PGS) for BMI were calculated. We modeled breastfeeding for one month or more as a function of BMI at enrollment; earliest BMI where available pre-pregnancy; and PGS for BMI. We conducted Mendelian randomization for breastfeeding initiation using PGS as an instrumental variable. Results: A higher BMI predicted a lower likelihood of breastfeeding for one month or more in all analyses. A +5 kg/m 2 BMI pre-pregnancy was associated with a 24% reduced odds of breastfeeding, and a +5 kg/m 2 genetically predicted BMI was associated with a 17% reduced odds of breastfeeding. Conclusions: BMI predicts a lower likelihood of breastfeeding for one month or longer. Given the high success of breastfeeding initiation regardless of BMI in supportive environments as well as potential health benefits, patients with elevated BMI may benefit from additional postpartum breastfeeding support.

Digital Footprints of Obesity Treatment: GLP-1 Receptor Agonists and the Health Equity Divide.

Our research investigates the societal implications of access to glucagon-like peptide-1 (GLP-1) agonists, particularly in light of recent clinical trials demonstrating the efficacy of semaglutide in reducing cardiovascular mortality. A decade-long analysis of Google Trends indicates a significant increase in searches for GLP-1 agonists, primarily in North America. This trend contrasts with the global prevalence of obesity. Given the high cost of GLP-1 agonists, a critical question arises: Will this disparity in medication accessibility exacerbate the global health equity gap in obesity treatment? This viewpoint explores strategies to address the health equity gap exacerbated by this emerging medication. Because GLP-1 agonists hold the potential to become a cornerstone in obesity treatment, ensuring equitable access is a pressing public health concern.

Diversity and scale: Genetic architecture of 2068 traits in the VA Million Veteran Program.

One of the justifiable criticisms of human genetic studies is the underrepresentation of participants from diverse populations. Lack of inclusion must be addressed at-scale to identify causal disease factors and understand the genetic causes of health disparities. We present genome-wide associations for 2068 traits from 635,969 participants in the Department of Veterans Affairs Million Veteran Program, a longitudinal study of diverse United States Veterans. Systematic analysis revealed 13,672 genomic risk loci; 1608 were only significant after including non-European populations. Fine-mapping identified causal variants at 6318 signals across 613 traits. One-third (n = 2069) were identified in participants from non-European populations. This reveals a broadly similar genetic architecture across populations, highlights genetic insights gained from underrepresented groups, and presents an extensive atlas of genetic associations.

Design and pilot results from the Million Veteran Program Return Of Actionable Results (MVP-ROAR) Study.

BACKGROUND: As a mega-biobank linked to a national healthcare system, the Million Veteran Program (MVP) can directly improve the health care of participants. To determine the feasibility and outcomes of returning medically actionable genetic results to MVP participants, the program launched the MVP Return of Actionable Results (MVP-ROAR) Study, with familial hypercholesterolemia (FH) as an exemplar actionable condition. METHODS: The MVP-ROAR Study consists of a completed single-arm pilot phase and an ongoing randomized clinical trial (RCT), in which MVP participants are recontacted and invited to receive clinical confirmatory gene sequencing testing and a telegenetic counseling intervention. The primary outcome of the RCT is 6-month change in low-density lipoprotein cholesterol (LDL-C) between participants receiving results at baseline and those receiving results after 6 months. RESULTS: The pilot developed processes to identify and recontact participants nationally with probable pathogenic variants in low-density lipoprotein receptor (LDLR) on the MVP genotype array, invite them to clinical confirmatory gene sequencing, and deliver a telegenetic counseling intervention. Among participants in the pilot phase, 8 (100%) had active statin prescriptions after 6 months. Results were shared with 16 first-degree family members. Six-month ΔLDL-C (low-density lipoprotein cholesterol) after the genetic counseling intervention was -37 mg/dL (95% CI: -12 to -61; P = .03). The ongoing RCT will determine between-arm differences in this primary outcome. CONCLUSION: While underscoring the importance of clinical confirmation of research results, the pilot phase of the MVP-ROAR Study marks a turning point in MVP and demonstrates the feasibility of returning genetic results to participants and their providers. The ongoing RCT will contribute to understanding how such a program might improve patient health care and outcomes. CLINICAL TRIAL REGISTRATION: ClinicalTrials.gov ID NCT04178122.

Identifying therapeutic targets for cancer among 2074 circulating proteins and risk of nine cancers.

Circulating proteins can reveal key pathways to cancer and identify therapeutic targets for cancer prevention. We investigate 2,074 circulating proteins and risk of nine common cancers (bladder, breast, endometrium, head and neck, lung, ovary, pancreas, kidney, and malignant non-melanoma) using cis protein Mendelian randomisation and colocalization. We conduct additional analyses to identify adverse side-effects of altering risk proteins and map cancer risk proteins to drug targets. Here we find 40 proteins associated with common cancers, such as PLAUR and risk of breast cancer [odds ratio per standard deviation increment: 2.27, 1.88-2.74], and with high-mortality cancers, such as CTRB1 and pancreatic cancer [0.79, 0.73-0.85]. We also identify potential adverse effects of protein-altering interventions to reduce cancer risk, such as hypertension. Additionally, we report 18 proteins associated with cancer risk that map to existing drugs and 15 that are not currently under clinical investigation. In sum, we identify protein-cancer links that improve our understanding of cancer aetiology. We also demonstrate that the wider consequence of any protein-altering intervention on well-being and morbidity is required to interpret any utility of proteins as potential future targets for therapeutic prevention.

Development and utility of a clinical research informatics application for participant recruitment and workflow management for a return of results pilot trial in familial hypercholesterolemia in the Million Veteran Program.

Objective: The development of clinical research informatics tools and workflow processes associated with re-engaging biobank participants has become necessary as genomic repositories increasingly consider the return of actionable research results. Materials and Methods: Here we describe the development and utility of an informatics application for participant recruitment and enrollment management for the Veterans Affairs Million Veteran Program Return Of Actionable Results Study, a randomized controlled pilot trial returning individual genetic results associated with familial hypercholesterolemia. Results: The application is developed in Python-Flask and was placed into production in November 2021. The application includes modules for chart review, medication reconciliation, participant contact and biospecimen logging, survey recording, randomization, and documentation of genetic counseling and result disclosure. Three primary users, a genetic counselor and two research coordinators, and 326 Veteran participants have been integrated into the system as of February 23, 2023. The application has successfully handled 3367 task requests involving greater than 95 000 structured data points. Specifically, application users have recorded 326 chart reviews, 867 recruitment telephone calls, 158 telephone-based surveys, and 61 return of results genetic counseling sessions, among other available study tasks. Conclusion: The development of usable, customizable, and secure informatics tools will become increasingly important as large genomic repositories begin to return research results at scale. Our work provides a proof-of-concept for developing and using such tools to aid in managing the return of results process within a national biobank.

A functional genomic framework to elucidate novel causal non-alcoholic fatty liver disease genes.

Background & Aims: Non-alcoholic fatty liver disease (NAFLD) is the most prevalent chronic liver pathology in western countries, with serious public health consequences. Efforts to identify causal genes for NAFLD have been hampered by the relative paucity of human data from gold-standard magnetic resonance quantification of hepatic fat. To overcome insufficient sample size, genome-wide association studies using NAFLD surrogate phenotypes have been used, but only a small number of loci have been identified to date. In this study, we combined GWAS of NAFLD composite surrogate phenotypes with genetic colocalization studies followed by functional in vitro screens to identify bona fide causal genes for NAFLD. Approach & Results: We used the UK Biobank to explore the associations of our novel NAFLD score, and genetic colocalization to prioritize putative causal genes for in vitro validation. We created a functional genomic framework to study NAFLD genes in vitro using CRISPRi. Our data identify VKORC1, TNKS, LYPLAL1 and GPAM as regulators of lipid accumulation in hepatocytes and suggest the involvement of VKORC1 in the lipid storage related to the development of NAFLD. Conclusions: Complementary genetic and genomic approaches are useful for the identification of NAFLD genes. Our data supports VKORC1 as a bona fide NAFLD gene. We have established a functional genomic framework to study at scale putative novel NAFLD genes from human genetic association studies.

Multi-Ancestry Polygenic Risk Score for Coronary Heart Disease Based on an Ancestrally Diverse Genome-Wide Association Study and Population-Specific Optimization.

BACKGROUND: Predictive performance of polygenic risk scores (PRS) varies across populations. To facilitate equitable clinical use, we developed PRS for coronary heart disease (CHD; PRSCHD) for 5 genetic ancestry groups. METHODS: We derived ancestry-specific and multi-ancestry PRSCHD based on pruning and thresholding (PRSPT) and ancestry-based continuous shrinkage priors (PRSCSx) applied to summary statistics from the largest multi-ancestry genome-wide association study meta-analysis for CHD to date, including 1.1 million participants from 5 major genetic ancestry groups. Following training and optimization in the Million Veteran Program, we evaluated the best-performing PRSCHD in 176,988 individuals across 9 diverse cohorts. RESULTS: Multi-ancestry PRSPT and PRSCSx outperformed ancestry-specific PRSPT and PRSCSx across a range of tuning values. Two best-performing multi-ancestry PRSCHD (ie, PRSPTmult and PRSCSxmult) and 1 ancestry-specific (PRSCSxEUR) were taken forward for validation. PRSPTmult demonstrated the strongest association with CHD in individuals of South Asian ancestry and European ancestry (odds ratio per 1 SD [95% CI, 2.75 [2.41-3.14], 1.65 [1.59-1.72]), followed by East Asian ancestry (1.56 [1.50-1.61]), Hispanic/Latino ancestry (1.38 [1.24-1.54]), and African ancestry (1.16 [1.11-1.21]). PRSCSxmult showed the strongest associations in South Asian ancestry (2.67 [2.38-3.00]) and European ancestry (1.65 [1.59-1.71]), lower in East Asian ancestry (1.59 [1.54-1.64]), Hispanic/Latino ancestry (1.51 [1.35-1.69]), and the lowest in African ancestry (1.20 [1.15-1.26]). CONCLUSIONS: The use of summary statistics from a large multi-ancestry genome-wide meta-analysis improved the performance of PRSCHD in most ancestry groups compared with single-ancestry methods. Despite the use of one of the largest and most diverse sets of training and validation cohorts to date, improvement of predictive performance was limited in African ancestry. This highlights the need for larger genome-wide association study datasets of underrepresented populations to enhance the performance of PRSCHD.

Diet Quality and Epigenetic Aging in the Women's Health Initiative.

BACKGROUND: Higher diet quality scores are associated with a lower risk for many chronic diseases and all-cause mortality; however, it is unclear if diet quality is associated with aging biology. OBJECTIVE: This study aimed to examine the association between diet quality and a measure of biological aging known as epigenetic aging. DESIGN: A cross-sectional data analysis was used to examine the association between three diet quality scores based on self-reported food frequency questionnaire data and five measures of epigenetic aging based on DNA methylation (DNAm) data from peripheral blood. PARTICIPANTS/SETTING: This study included 4,500 postmenopausal women recruited from multiple sites across the United States (1993-98), aged 50 to 79 years, with food frequency questionnaire and DNAm data available from the Women's Health Initiative baseline visit. MAIN OUTCOME MEASURES: Five established epigenetic aging measures were generated from HumanMethylation450 Beadchip DNAm data, including AgeAccelHannum, AgeAccelHorvath, AgeAccelPheno, AgeAccelGrim, and DunedinPACE. STATISTICAL ANALYSES PERFORMED: Linear mixed models were used to test for associations between three diet quality scores (Healthy Eating Index, Dietary Approaches to Stop Hypertension, and alternate Mediterranean diet scores) and epigenetic aging measures, adjusted for age, race and ethnicity, education, tobacco smoking, physical activity, Women's Health Initiative substudy from which DNAm data were obtained, and DNAm-based estimates of leukocyte proportions. RESULTS: Healthy Eating Index, Dietary Approaches to Stop Hypertension, and alternate Mediterranean diet scores were all inversely associated with AgeAccelPheno, AgeAccelGrim, and DunedinPACE (P < 0.05), with the largest effects with DunedinPACE. A one standard deviation increment in diet quality scores was associated with a decrement (ß ± SE) in DunedinPACE z score of -0.097 ± 0.014 (P = 9.70 x 10-13) for Healthy Eating Index, -0.107 ± 0.014 (P = 1.53 x 10-14) for Dietary Approaches to Stop Hypertension, and -0.068 ± 0.013 (P = 2.31 x 10-07) for the alternate Mediterranean diet. CONCLUSIONS: In postmenopausal women, diet quality scores were inversely associated with DNAm-based measures of biological aging, particularly DunedinPACE.

Cardiorespiratory fitness and risk of heart failure with preserved ejection fraction.

AIMS: Preventive strategies for heart failure with preserved ejection fraction (HFpEF) include pharmacotherapies and lifestyle modifications. However, the association between cardiorespiratory fitness (CRF) assessed objectively by a standardized exercise treadmill test (ETT) and the risk of HFpEF has not been evaluated. Thus, we evaluated the association between CRF and HFpEF incidence. METHODS AND RESULTS: We assessed CRF in US Veterans (624 551 men; mean age 61.2 ± 9.7 years and 43 179 women; mean age 55.0 ± 8.9 years) by a standardized ETT performed between 1999 and 2020 across US Veterans Affairs Medical Centers. All had no evidence of heart failure or myocardial infarction prior to completion of the ETT. We assigned participants to one of five age- and gender-specific CRF categories (quintiles) based on peak metabolic equivalents (METs) achieved during the ETT and four categories based on CRF changes in those with two ETT evaluations (n = 139 434) ≥1.0 year apart. During a median follow-up of 10.1 years (interquartile range 6.0-14.3 years), providing 6 879 229 person-years, there were 16 493 HFpEF events with an average annual rate of 2.4 events per 1000 person-years. The adjusted risk of HFpEF decreased across CRF categories as CRF increased, independent of comorbidities. For fit individuals (≥10.5 METs) the hazard ratio (HR) was 0.48 (95% confidence interval [CI] 0.46-0.51) compared with least fit (≤4.9 METs; referent). Being unfit carried the highest risk (HR 2.88, 95% CI 2.67-3.11) of any other comorbidity. The risk of unfit individuals who became fit was 37% lower (HR 0.63, 95% CI 0.57-0.71), compared to those who remained unfit. CONCLUSIONS: Higher CRF levels are independently associated with lower HFpEF in a dose-response manner. Changes in CRF reflected proportional changes in HFpEF risk, suggesting that the HFpEF risk was modulated by CRF.

Mendelian randomization analyses suggest a causal role for circulating GIP and IL-1RA levels in homeostatic model assessment-derived measures of ß-cell function and insulin sensitivity in Africans without type 2 diabetes.

BACKGROUND: In vitro and in vivo studies have shown that certain cytokines and hormones may play a role in the development and progression of type 2 diabetes (T2D). However, studies on their role in T2D in humans are scarce. We evaluated associations between 11 circulating cytokines and hormones with T2D among a population of sub-Saharan Africans and tested for causal relationships using Mendelian randomization (MR) analyses. METHODS: We used logistic regression analysis adjusted for age, sex, body mass index, and recruitment country to regress levels of 11 cytokines and hormones (adipsin, leptin, visfatin, PAI-1, GIP, GLP-1, ghrelin, resistin, IL-6, IL-10, IL-1RA) on T2D among Ghanaians, Nigerians, and Kenyans from the Africa America Diabetes Mellitus study including 2276 individuals with T2D and 2790 non-T2D individuals. Similar linear regression models were fitted with homeostatic modelling assessments of insulin sensitivity (HOMA-S) and ß-cell function (HOMA-B) as dependent variables among non-T2D individuals (n = 2790). We used 35 genetic variants previously associated with at least one of these 11 cytokines and hormones among non-T2D individuals as instrumental variables in univariable and multivariable MR analyses. Statistical significance was set at 0.0045 (0.05/11 cytokines and hormones). RESULTS: Circulating GIP and IL-1RA levels were associated with T2D. Nine of the 11 cytokines and hormones (exceptions GLP-1 and IL-6) were associated with HOMA-S, HOMA-B, or both among non-T2D individuals. Two-stage least squares MR analysis provided evidence for a causal effect of GIP and IL-RA on HOMA-S and HOMA-B in multivariable analyses (GIP ~ HOMA-S ß = - 0.67, P-value = 1.88 × 10-6 and HOMA-B ß = 0.59, P-value = 1.88 × 10-5; IL-1RA ~ HOMA-S ß = - 0.51, P-value = 8.49 × 10-5 and HOMA-B ß = 0.48, P-value = 5.71 × 10-4). IL-RA was partly mediated via BMI (30-34%), but GIP was not. Inverse variance weighted MR analysis provided evidence for a causal effect of adipsin on T2D (multivariable OR = 1.83, P-value = 9.79 × 10-6), though these associations were not consistent in all sensitivity analyses. CONCLUSIONS: The findings of this comprehensive MR analysis indicate that circulating GIP and IL-1RA levels are causal for reduced insulin sensitivity and increased ß-cell function. GIP's effect being independent of BMI suggests that circulating levels of GIP could be a promising early biomarker for T2D risk. Our MR analyses do not provide conclusive evidence for a causal role of other circulating cytokines in T2D among sub-Saharan Africans.

Plasma Protein Profiling of Incident Cardiovascular Diseases: A Multisample Evaluation.

BACKGROUND: Proteomic profiling could potentially disclose new pathophysiological pathways for cardiovascular diseases (CVD) and improve prediction at the individual level. We therefore aimed to study the plasma protein profile associated with the incidence of different CVDs. METHODS: Plasma levels of 245 proteins suspected to be linked to CVD or metabolism were measured in 4 Swedish prospective population-based cohorts (SIMPLER [Swedish Infrastructure for Medical Population-Based Life-Course and Environmental Research], ULSAM (Uppsala Longitudinal Study of Adult Men), EpiHealth, and POEM [Prospective Investigation of Obesity, Energy Production, and Metabolism]) comprising 11 869 individuals, free of CVD diagnoses at baseline. Our primary CVD outcome was defined by a combined end point that included either incident myocardial infarction, stroke, or heart failure. RESULTS: Using a discovery/validation approach, 42 proteins were associated with our primary composite end point occurring in 1163 subjects. In separate meta-analyses for each of the 3 CVD outcomes, 49 proteins were related to myocardial infarction, 34 to ischemic stroke, and 109 to heart failure. Thirteen proteins were related to all 3 outcomes. Of those, urokinase plasminogen activator surface receptor, adrenomedullin, and KIM-1 (kidney injury molecule 1) were also related to several markers of subclinical CVD in Prospective Investigation of Obesity, Energy production and Metabolism, reflecting myocardial or arterial pathologies. In prediction analysis, a lasso selection of 11 proteins in ULSAM improved the discrimination of CVD by 3.3% (P<0.0001) in SIMPLER when added to traditional risk factors. CONCLUSIONS: Protein profiling in multiple samples disclosed several new proteins to be associated with subsequent myocardial infarction, stroke, and heart failure, suggesting common pathophysiological pathways for these diseases. KIM-1, urokinase plasminogen activator surface receptor, and adrenomedullin were novel early markers of CVD. A selection of 11 proteins improved the discrimination of CVD.

CYP2C19 Polymorphisms and Clinical Outcomes Following Percutaneous Coronary Intervention (PCI) in the Million Veterans Program.

Background: CYP2C19 loss-of-function (LOF) alleles decrease the antiplatelet effect of clopidogrel following percutaneous coronary intervention (PCI) in patients presenting with acute coronary syndrome (ACS). The impact of genotype in stable ischemic heart disease (SIHD) is unclear. Objectives: Determine the association of CYP2C19 genotype with major adverse cardiac events (MACE) after PCI for ACS or SIHD. Methods: Million Veterans Program (MVP) participants age <65 years with a PCI documented in the VA Clinical Assessment, Reporting and Tracking (CART) Program between 1/1/2009 to 9/30/2017, treated with clopidogrel were included. Time to MACE defined as the composite of all-cause death, stroke or myocardial infarction within 12 months following PCI. Results: Among 4,461 Veterans (mean age 59.1 ± 5.1 years, 18% Black); 44% had ACS, 56% had SIHD and 29% carried a CYP2C19 LOF allele. 301 patients (6.7%) experienced MACE while being treated with clopidogrel, 155 (7.9%) in the ACS group and 146 (5.9%) in the SIHD group. Overall, MACE was not significantly different between LOF carriers vs. noncarriers (adjusted hazard ratio [HR] 1.18, confidence interval [95%CI] 0.97-1.45, p=0.096). Among patients presenting with ACS, MACE risk in LOF carriers versus non-carriers was numerically higher (HR 1.30, 95%CI 0.98-1.73, p=0.067). There was no difference in MACE risk in patients with SIHD (HR 1.09, 95%CI 0.82-1.44; p=0.565). Conclusions: CYP2C19 LOF carriers presenting with ACS treated with clopidogrel following PCI experienced a numerically greater elevated risk of MACE events. CYP2C19 LOF genotype is not associated with MACE among patients presenting with SIHD.

CXCL12 regulates coronary artery dominance in diverse populations and links development to disease.

Mammalian cardiac muscle is supplied with blood by right and left coronary arteries that form branches covering both ventricles of the heart. Whether branches of the right or left coronary arteries wrap around to the inferior side of the left ventricle is variable in humans and termed right or left dominance. Coronary dominance is likely a heritable trait, but its genetic architecture has never been explored. Here, we present the first large-scale multi-ancestry genome-wide association study of dominance in 61,043 participants of the VA Million Veteran Program, including over 10,300 Africans and 4,400 Admixed Americans. Dominance was moderately heritable with ten loci reaching genome wide significance. The most significant mapped to the chemokine CXCL12 in both Europeans and Africans. Whole-organ imaging of human fetal hearts revealed that dominance is established during development in locations where CXCL12 is expressed. In mice, dominance involved the septal coronary artery, and its patterning was altered with Cxcl12 deficiency. Finally, we linked human dominance patterns with coronary artery disease through colocalization, genome-wide genetic correlation and Mendelian Randomization analyses. Together, our data supports CXCL12 as a primary determinant of coronary artery dominance in humans of diverse backgrounds and suggests that developmental patterning of arteries may influence one's susceptibility to ischemic heart disease.

Assessing efficiency of fine-mapping obesity-associated variants through leveraging ancestry architecture and functional annotation using PAGE and UKBB cohorts.

Inadequate representation of non-European ancestry populations in genome-wide association studies (GWAS) has limited opportunities to isolate functional variants. Fine-mapping in multi-ancestry populations should improve the efficiency of prioritizing variants for functional interrogation. To evaluate this hypothesis, we leveraged ancestry architecture to perform comparative GWAS and fine-mapping of obesity-related phenotypes in European ancestry populations from the UK Biobank (UKBB) and multi-ancestry samples from the Population Architecture for Genetic Epidemiology (PAGE) consortium with comparable sample sizes. In the investigated regions with genome-wide significant associations for obesity-related traits, fine-mapping in our ancestrally diverse sample led to 95% and 99% credible sets (CS) with fewer variants than in the European ancestry sample. Lead fine-mapped variants in PAGE regions had higher average coding scores, and higher average posterior probabilities for causality compared to UKBB. Importantly, 99% CS in PAGE loci contained strong expression quantitative trait loci (eQTLs) in adipose tissues or harbored more variants in tighter linkage disequilibrium (LD) with eQTLs. Leveraging ancestrally diverse populations with heterogeneous ancestry architectures, coupled with functional annotation, increased fine-mapping efficiency and performance, and reduced the set of candidate variants for consideration for future functional studies. Significant overlap in genetic causal variants across populations suggests generalizability of genetic mechanisms underpinning obesity-related traits across populations.

Proteomic analysis of 92 circulating proteins and their effects in cardiometabolic diseases.

BACKGROUND: Human plasma contains a wide variety of circulating proteins. These proteins can be important clinical biomarkers in disease and also possible drug targets. Large scale genomics studies of circulating proteins can identify genetic variants that lead to relative protein abundance. METHODS: We conducted a meta-analysis on genome-wide association studies of autosomal chromosomes in 22,997 individuals of primarily European ancestry across 12 cohorts to identify protein quantitative trait loci (pQTL) for 92 cardiometabolic associated plasma proteins. RESULTS: We identified 503 (337 cis and 166 trans) conditionally independent pQTLs, including several novel variants not reported in the literature. We conducted a sex-stratified analysis and found that 118 (23.5%) of pQTLs demonstrated heterogeneity between sexes. The direction of effect was preserved but there were differences in effect size and significance. Additionally, we annotate trans-pQTLs with nearest genes and report plausible biological relationships. Using Mendelian randomization, we identified causal associations for 18 proteins across 19 phenotypes, of which 10 have additional genetic colocalization evidence. We highlight proteins associated with a constellation of cardiometabolic traits including angiopoietin-related protein 7 (ANGPTL7) and Semaphorin 3F (SEMA3F). CONCLUSION: Through large-scale analysis of protein quantitative trait loci, we provide a comprehensive overview of common variants associated with plasma proteins. We highlight possible biological relationships which may serve as a basis for further investigation into possible causal roles in cardiometabolic diseases.

A Multi-Ancestry Polygenic Risk Score for Coronary Heart Disease Based on an Ancestrally Diverse Genome-Wide Association Study and Population-Specific Optimization.

Background: Predictive performance of polygenic risk scores (PRS) varies across populations. To facilitate equitable clinical use, we developed PRS for coronary heart disease (PRSCHD) for 5 genetic ancestry groups. Methods: We derived ancestry-specific and multi-ancestry PRSCHD based on pruning and thresholding (PRSP+T) and continuous shrinkage priors (PRSCSx) applied on summary statistics from the largest multi-ancestry genome-wide meta-analysis for CHD to date, including 1.1 million participants from 5 continental populations. Following training and optimization of PRSCHD in the Million Veteran Program, we evaluated predictive performance of the best performing PRSCHD in 176,988 individuals across 9 cohorts of diverse genetic ancestry. Results: Multi-ancestry PRSP+T outperformed ancestry specific PRSP+T across a range of tuning values. In training stage, for all ancestry groups, PRSCSx performed better than PRSP+T and multi-ancestry PRS outperformed ancestry-specific PRS. In independent validation cohorts, the selected multi-ancestry PRSP+T demonstrated the strongest association with CHD in individuals of South Asian (SAS) and European (EUR) ancestry (OR per 1SD[95% CI]; 2.75[2.41-3.14], 1.65[1.59-1.72]), followed by East Asian (EAS) (1.56[1.50-1.61]), Hispanic/Latino (HIS) (1.38[1.24-1.54]), and weakest in African (AFR) ancestry (1.16[1.11-1.21]). The selected multi-ancestry PRSCSx showed stronger associacion with CHD in comparison within each ancestry group where the association was strongest in SAS (2.67[2.38-3.00]) and EUR (1.65[1.59-1.71]), progressively decreasing in EAS (1.59[1.54-1.64]), HIS (1.51[1.35-1.69]), and lowest in AFR (1.20[1.15-1.26]). Conclusions: Utilizing diverse summary statistics from a large multi-ancestry genome-wide meta-analysis led to improved performance of PRSCHD in most ancestry groups compared to single-ancestry methods. Improvement of predictive performance was limited, specifically in AFR and HIS, despite use of one of the largest and most diverse set of training and validation cohorts to date. This highlights the need for larger GWAS datasets of AFR and HIS individuals to enhance performance of PRSCHD.