Feng, QiPing

Background: Susceptibility to cardiometabolic diseases (CDs) is shaped by genetic, socio-environmental, and lifestyle factors. We evaluated the associations of disease-specific polysocial risk scores (PSSs)— derived from social determinants of health (SDOH) and lifestyle factors—and genetic risk factors with six CDs, across self-identified race/ethnicity (SIRE) groups in the eMERGE IV cohort. Methods: The six CDs included: atrial fibrillation, coronary heart disease, chronic kidney disease, hyperlipidemia, obesity, and type 2 diabetes. Disease-specific PSSs were developed using machine-learning models based on 65 SDOH and lifestyle factors. Participants were grouped based on their SIRE as Asian, Black, Hispanic/Latino, and White. We evaluated PSS distributions across SIRE groups. Genetic risk was assessed through disease-specific polygenic risk scores (PRSs), monogenic etiologies, and family history (FamHx). High PSS and high PRS were defined as the top 5th percentile. Associations were tested using multivariable logistic regression; model performance was assessed using the C-statistic. Results: Among 20178 participants (50±15 years, 68% female, 40% non-White), the most frequently associated PSS components were physical activity, sleep, self-perceived health, and income. PSS distributions varied across SIRE groups, with Whites having the lowest and Black and Hispanic/Latinos the highest scores for most CDs (Figure 1). All PSSs were significantly associated with their respective CDs, independent of confounders (OR per 1 SD increase: 1.19–2.68; Figure 2-A). No significant interactions were observed between PSSs and PRSs or monogenic etiologies. However, significant interactions were noted between PSSs and FamHx for CKD and obesity. PSS and PRS effects on the odds of CDs were independent (Figure 2-B) and additive (Figure 2-C). Incorporating PSSs into models that included clinical and genetic risk factors improved predictive performance of risk prediction models (Figure 3). Conclusion: PSSs, reflecting socio-environmental and lifestyle disadvantage, varied significantly across SIRE groups and were associated with risk of six CDs, independent of and additive to PRS. Given the substantial variation in socio-environmental and lifestyle risks across SIRE groups, integrating these factors into risk prediction models for CDs may help address health disparities and enhance preventive efforts, especially in marginalized populations.

Background: We investigated the impact of genetic risk factors for coronary heart disease (CHD) -polygenic risk score (PRS), familial hypercholesterolemia (FH), and family history (FamHx)- on CHD risk estimates, across the age spectrum, in two diverse cohorts of US adults -eMERGE IV (eIV) and All of Us (AoU). Methods: CHD was defined as myocardial infarction, unstable angina, and coronary revascularization. Self-identified race/ethnicity (SIRE) was used as a population descriptor. We calculated a PRS for CHD (PRS CHD , PGS004698), ascertained FH as the presence of pathogenic/likely pathogenic variants in FH genes, and defined FamHx as early-onset CHD in first-degree relatives. We employed Pooled Cohort Equations (PCE) to estimate the 10y risk of CHD for adults ≥40y and modeled the association of conventional risk factors with CHD in adults <40y. AoU served as the training set, and eIV as the testing set. We analyzed the impact of PRS CHD and FamHx on CHD risk estimates by a) using multivariable logistic regression and Cox proportional hazard models, assessing discrimination and the extent of risk reclassification; and b) net benefit analysis and decision curves to assess the performance of prediction models across actionable thresholds. Results: We analyzed data for 19348 participants from eIV (age 50±15, 68% female, 41% non-White) and 239645 participants from AoU (age 55±17, 61% female, 48% non-White). Genetic risk factors were significantly associated with CHD. PRS CHD performance varied by SIRE groups, while FamHx was consistent. The effects of PRS CHD and FamHx on CHD were independent and additive (Figure 1). In adults ≥40y, incorporating PRS CHD and FamHx into PCE improved discrimination (C-statistic increased from 0.719 to 0.753; P -diff=9.1×10 -3 , Figure 2) and reclassified risk in 19% and 20% of participants at the 7.5% and 10% 10y CHD risk thresholds, respectively. Between the 7.5% and 10% 10y CHD risk thresholds, incorporating PRS CHD and FamHx into the PCE improved the net benefit of the risk prediction models across White, Black, and Hispanic/Latino groups (Figure 3). Conclusion: PRS CHD and FamHx were independently and additively associated with CHD in two large diverse cohorts in the US. Incorporating PRS CHD and FamHx into PCE improved risk discrimination, reclassified risk in a significant portion of participants, and improved net benefit of the PCE across all three major SIRE groups, motivating the addition of these factors to clinical risk calculators.

Apart from ancestry, personal or environmental covariates may contribute to differences in polygenic score (PGS) performance. We analyzed the effects of covariate stratification and interaction on body mass index (BMI) PGS (PGS BMI ) across four cohorts of European (N = 491,111) and African (N = 21,612) ancestry. Stratifying on binary covariates and quintiles for continuous covariates, 18/62 covariates had significant and replicable R 2 differences among strata. Covariates with the largest differences included age, sex, blood lipids, physical activity, and alcohol consumption, with R 2 being nearly double between best- and worst-performing quintiles for certain covariates. Twenty-eight covariates had significant PGS BMI –covariate interaction effects, modifying PGS BMI effects by nearly 20% per standard deviation change. We observed overlap between covariates that had significant R 2 differences among strata and interaction effects – across all covariates, their main effects on BMI were correlated with their maximum R 2 differences and interaction effects (0.56 and 0.58, respectively), suggesting high-PGS BMI individuals have highest R 2 and increase in PGS effect. Using quantile regression, we show the effect of PGS BMI increases as BMI itself increases, and that these differences in effects are directly related to differences in R 2 when stratifying by different covariates. Given significant and replicable evidence for context-specific PGS BMI performance and effects, we investigated ways to increase model performance taking into account nonlinear effects. Machine learning models (neural networks) increased relative model R 2 (mean 23%) across datasets. Finally, creating PGS BMI directly from GxAge genome-wide association studies effects increased relative R 2 by 7.8%. These results demonstrate that certain covariates, especially those most associated with BMI, significantly affect both PGS BMI performance and effects across diverse cohorts and ancestries, and we provide avenues to improve model performance that consider these effects.

Apart from ancestry, personal or environmental covariates may contribute to differences in polygenic score (PGS) performance. We analyzed the effects of covariate stratification and interaction on body mass index (BMI) PGS (PGS BMI ) across four cohorts of European (N = 491,111) and African (N = 21,612) ancestry. Stratifying on binary covariates and quintiles for continuous covariates, 18/62 covariates had significant and replicable R 2 differences among strata. Covariates with the largest differences included age, sex, blood lipids, physical activity, and alcohol consumption, with R 2 being nearly double between best- and worst-performing quintiles for certain covariates. Twenty-eight covariates had significant PGS BMI –covariate interaction effects, modifying PGS BMI effects by nearly 20% per standard deviation change. We observed overlap between covariates that had significant R 2 differences among strata and interaction effects – across all covariates, their main effects on BMI were correlated with their maximum R 2 differences and interaction effects (0.56 and 0.58, respectively), suggesting high-PGS BMI individuals have highest R 2 and increase in PGS effect. Using quantile regression, we show the effect of PGS BMI increases as BMI itself increases, and that these differences in effects are directly related to differences in R 2 when stratifying by different covariates. Given significant and replicable evidence for context-specific PGS BMI performance and effects, we investigated ways to increase model performance taking into account nonlinear effects. Machine learning models (neural networks) increased relative model R 2 (mean 23%) across datasets. Finally, creating PGS BMI directly from GxAge genome-wide association studies effects increased relative R 2 by 7.8%. These results demonstrate that certain covariates, especially those most associated with BMI, significantly affect both PGS BMI performance and effects across diverse cohorts and ancestries, and we provide avenues to improve model performance that consider these effects.

Background: Lp(a) (lipoprotein [a]) levels are higher in individuals of African ancestry (AA) than in individuals of European ancestry (EA). We examined associations of genetically predicted Lp(a) levels with (1) atherosclerotic cardiovascular disease subtypes: coronary heart disease, cerebrovascular disease, peripheral artery disease, and abdominal aortic aneurysm and (2) nonatherosclerotic cardiovascular disease phenotypes, stratified by ancestry. Methods: We performed (1) Mendelian randomization analyses for previously reported cardiovascular associations and (2) Mendelian randomization–phenome-wide association analyses for novel associations. Analyses were stratified by ancestry in electronic Medical Records and Genomics, United Kingdom Biobank, and Million Veteran Program cohorts separately and in a combined cohort of 804 507 EA and 103 580 AA participants. Results: In Mendelian randomization analyses using the combined cohort, a 1-SD genetic increase in Lp(a) level was associated with atherosclerotic cardiovascular disease subtypes in EA—odds ratio and 95% CI for coronary heart disease 1.28 (1.16–1.41); cerebrovascular disease 1.14 (1.07–1.21); peripheral artery disease 1.22 (1.11–1.34); abdominal aortic aneurysm 1.28 (1.17–1.40); in AA, the effect estimate was lower than in EA and nonsignificant for coronary heart disease 1.11 (0.99–1.24) and cerebrovascular disease 1.06 (0.99–1.14) but similar for peripheral artery disease 1.16 (1.01–1.33) and abdominal aortic aneurysm 1.34 (1.11–1.62). In EA, a 1-SD genetic increase in Lp(a) level was associated with aortic valve disorders 1.34 (1.10–1.62), mitral valve disorders 1.18 (1.09–1.27), congestive heart failure 1.12 (1.05–1.19), and chronic kidney disease 1.07 (1.01–1.14). In AA, no significant associations were noted for aortic valve disorders 1.08 (0.94–1.25), mitral valve disorders 1.02 (0.89–1.16), congestive heart failure 1.02 (0.95–1.10), or chronic kidney disease 1.05 (0.99–1.12). Mendelian randomization–phenome-wide association analyses identified novel associations in EA with arterial thromboembolic disease, nonaortic aneurysmal disease, atrial fibrillation, cardiac conduction disorders, and hypertension. Conclusions: Many cardiovascular associations of genetically increased Lp(a) that were significant in EA were not significant in AA. Lp(a) was associated with atherosclerotic cardiovascular disease in four major arterial beds in EA but only with peripheral artery disease and abdominal aortic aneurysm in AA. Additionally, novel cardiovascular associations were detected in EA.

Objectives Betaine‐homocysteine methyltransferase (BHMT) catalyzes the remethylation of homocysteine. BHMT in the human liver is responsible for up to 25% of homocysteine metabolism. Knowledge of the relationship between genetic polymorphisms and BHMT phenotypes could enhance our understanding of the role of this reaction in drug metabolism. Methods Levels of BHMT enzyme activity and immunoreactive protein were measured in 268 human liver biopsies. Twelve tag BHMT SNPs from our own gene resequencing and HapMap data were used to genotype DNA from these samples. Genotype‐phenotype association analysis was then performed. BHMT enzyme activity averaged 0.077 ± 0.028 units/μg, with a 12 fold variation. Results BHMT protein levels averaged 6.9 ± 3.3 units/μg, with a 50‐fold variation. Three SNPs (rs41272270, rs16876512, and rs6875201), located 28kb upstream, in the 5′‐UTR and in intron 1 of BHMT , respectively, were significantly correlated with both phenotypes (p<10 −8 for enzyme activity and p<10 −5 for protein level). rs41272270 was associated with altered protein binding during electrophoresis mobility shift assay. Conclusions Three tag SNPs were significantly associated with human hepatic BHMT enzyme activity and immunoreactive protein. These observations may help us understand the regulation of BHMT in the human liver and its role in drug metabolism. Supported by NIH U01 GM61388, R01 GM28157, R01 CA132780.

Background: The pathophysiology of lipoprotein (a) [Lp(a)] remains obscure. We conducted a phenome wide analysis study (PheWAS) to identify pleiotropic effects of LPA variants. Methods: Among 51,843 European-ancestry participants (age 58±16 years, 54% female) of the electronic MEdical Records and Genomics (eMERGE) network we assessed whether LPA variants exhibited pleiotropic affects through an electronic health record-based PheWAS. We adjusted significant associations by presence of atherosclerotic cardiovascular disease (ASCVD; defined as a composite of coronary heart disease, cerebrovascular disease, and peripheral artery disease) to determine whether those associations were independent of ASCVD. Results: In PheWAS analysis, we tested 864 phecodes with 36 independent LPA variants. We identified 21 significant associations with non-ASCVD phenotypes including heart failure, aortic valve stenosis, and chronic kidney disease of which 14 were replicated across independent cohorts. The strength of associations between the LPA variant rs10455872 and both heart failure and aortic valve stenosis related phecodes were significantly attenuated after adjustment for ASCVD. However, the association with chronic kidney disease phecode remained independent following adjustment for ASCVD without any attenuation in the estimated odds ratio (Table). Conclusions: LPA genetic variants were associated with aortic valve stenosis, heart failure and chronic kidney disease. The observed association of LPA variant rs10455872 with aortic stenosis and heart failure appear to be partially mediated by ASCVD, while the association with chronic kidney disease appears to be independent of ASCVD. Additional studies are needed to elucidate potential mechanisms by which Lp(a) may contribute to the pathogenesis of non-ASCVD disease phenotypes.

AbstractCarotid artery atherosclerotic disease (CAAD) is a risk factor for stroke. We used a genome‐wide association (GWAS) approach to discover genetic variants associated with CAAD in participants in the electronic Medical Records and Genomics (eMERGE) Network. We identified adult CAAD cases with unilateral or bilateral carotid artery stenosis and controls without evidence of stenosis from electronic health records at eight eMERGE sites. We performed GWAS with a model adjusting for age, sex, study site, and genetic principal components of ancestry. In eMERGE we found 1793 CAAD cases and 17,958 controls. Two loci reached genome‐wide significance, on chr6 in LPA (rs10455872, odds ratio [OR] (95% confidence interval [CI]) = 1.50 (1.30–1.73), p = 2.1 × 10−8) and on chr7, an intergenic single nucleotide variant (SNV; rs6952610, OR (95% CI) = 1.25 (1.16–1.36), p = 4.3 × 10−8). The chr7 association remained significant in the presence of the LPA SNV as a covariate. The LPA SNV was also associated with coronary heart disease (CHD; 4199 cases and 11,679 controls) in this study (OR (95% CI) = 1.27 (1.13–1.43), p = 5 × 10−5) but the chr7 SNV was not (OR (95% CI) = 1.03 (0.97–1.09), p = .37). Both variants replicated in UK Biobank. Elevated lipoprotein(a) concentrations ([Lp(a)]) and LPA variants associated with elevated [Lp(a)] have previously been associated with CAAD and CHD, including rs10455872. With electronic health record phenotypes in eMERGE and UKB, we replicated a previously known association and identified a novel locus associated with CAAD.

Background: Coronary heart disease (CHD) is a leading cause of death globally. Although therapy with statins decreases circulating levels of low-density lipoprotein cholesterol and the incidence of CHD, additional events occur despite statin therapy in some individuals. The genetic determinants of this residual cardiovascular risk remain unknown. Methods: We performed a 2-stage genome-wide association study of CHD events during statin therapy. We first identified 3099 cases who experienced CHD events (defined as acute myocardial infarction or the need for coronary revascularization) during statin therapy and 7681 controls without CHD events during comparable intensity and duration of statin therapy from 4 sites in the Electronic Medical Records and Genomics Network. We then sought replication of candidate variants in another 160 cases and 1112 controls from a fifth Electronic Medical Records and Genomics site, which joined the network after the initial genome-wide association study. Finally, we performed a phenome-wide association study for other traits linked to the most significant locus. Results: The meta-analysis identified 7 single nucleotide polymorphisms at a genome-wide level of significance within the LPA/PLG locus associated with CHD events on statin treatment. The most significant association was for an intronic single nucleotide polymorphism within LPA/PLG (rs10455872; minor allele frequency, 0.069; odds ratio, 1.58; 95% confidence interval, 1.35–1.86; P =2.6×10 − 10 ). In the replication cohort, rs10455872 was also associated with CHD events (odds ratio, 1.71; 95% confidence interval, 1.14–2.57; P =0.009). The association of this single nucleotide polymorphism with CHD events was independent of statin-induced change in low-density lipoprotein cholesterol (odds ratio, 1.62; 95% confidence interval, 1.17–2.24; P =0.004) and persisted in individuals with low-density lipoprotein cholesterol ≤70 mg/dL (odds ratio, 2.43; 95% confidence interval, 1.18–4.75; P =0.015). A phenome-wide association study supported the effect of this region on coronary heart disease and did not identify noncardiovascular phenotypes. Conclusions: Genetic variations at the LPA locus are associated with CHD events during statin therapy independently of the extent of low-density lipoprotein cholesterol lowering. This finding provides support for exploring strategies targeting circulating concentrations of lipoprotein(a) to reduce CHD events in patients receiving statins.