Smith, Johanna L.

BACKGROUND: Performance and transferability of contemporary polygenic risk scores (PRS) for atherosclerotic cardiovascular disease phenotypes may vary across PRS methods, training data, and trait ascertainment. METHODS: We aimed to investigate the performance and transferability of contemporary PRS for atherosclerotic cardiovascular disease subtypes: coronary heart disease (CHD), abdominal aortic aneurysm (AAA), ischemic stroke (IS), and peripheral artery disease (PAD), using the All of Us Workbench, which consists of a large, diverse cohort with whole-genome sequence data. We also developed and evaluated a multi-trait PRS for each subtype. Performance of PRS for 4 atherosclerotic cardiovascular disease subtypes was compared across genetic similarity groups in 245 388 All of Us participants. Groups genetically similar to European, African, admixed American, and remaining groups (combined as other) were used to assess PRS for CHD, IS, AAA, PAD, and multi-trait. RESULTS: PRS for CHD and AAA performed better than IS and PAD. For CHD, CHD PGS003725 performed the best (hazard ratio per SD increase [95% CI]), across genetic ancestry groups, European, and African (1.72 [1.67–1.78], 1.23 [1.17–1.29]), with CHD PGS004696 being best for admixed American (1.91 [1.70–2.15]), and CHD PGS003356 for other (1.75 [1.58–1.95]). The best performing PRS for AAA was AAA Multi for European, other, and admixed American (1.71 [1.52–1.92], 1.59 [1.07–2.37], 1.50 [0.90–2.52]) and AAA PGS003972 for African (1.39 [1.19–1.63]). For IS, IS Multi performed best for other and European (1.49 [1.17–1.89], 1.33 [1.25–1.42]), and IS PGS000039 performed best in admixed American and African (1.17 [1.07–1.27], 1.09 [1.04–1.15]). For PAD, PAD Multi performed best for all groups (other, 1.51 [1.19–1.92]; European, 1.32 [1.24–1.41]; admixed American, 1.23 [1.05–1.45]; and African, 1.18 [1.04–1.34]). CONCLUSIONS: Multi-trait and multi-ancestry PRS performed better than individual trait and/or single ancestry PRS for each atherosclerotic cardiovascular disease phenotype across ancestrally diverse and admixed individuals, with minimal change including adjustment for conventional risk factors.

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.

BACKGROUND: Predictive performance of polygenic risk scores (PRS) varies across populations. To facilitate equitable clinical use, we developed PRS for coronary heart disease (CHD; PRS CHD ) for 5 genetic ancestry groups. METHODS: We derived ancestry-specific and multi-ancestry PRS CHD based on pruning and thresholding (PRS PT ) and ancestry-based continuous shrinkage priors (PRS CSx ) 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 PRS CHD in 176,988 individuals across 9 diverse cohorts. RESULTS: Multi-ancestry PRS PT and PRS CSx outperformed ancestry-specific PRS PT and PRS CSx across a range of tuning values. Two best-performing multi-ancestry PRS CHD (ie, PRS PTmult and PRS CSxmult ) and 1 ancestry-specific (PRS CSxEUR ) were taken forward for validation. PRS PTmult 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]). PRS CSxmult 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 PRS CHD 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 PRS CHD .

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.

Introduction: Polygenic risk scores (PRS) for coronary heart disease (CHD) are associated with incident CHD events. We assessed whether the strength of this association differs over the life course. Methods: We studied 7,958 genotyped participants of European-ancestry without CHD enrolled in the Atherosclerosis Risk in Communities Study (ARIC). Participants were 45-64 years old at the baseline visit occurring 1987-89 with follow up available through 2015. A genome-wide PRS for CHD composed of 1.7 million variants was computed for each participant. Incident CHD was defined as clinically adjudicated myocardial infarction including silent infarction, coronary revascularization, or fatal coronary event. An extended Cox model was constructed with age as the model time scale. The association of PRS with incident CHD was assumed to be log-linear and further modeled as a function of age assuming differing piecewise constant effect sizes across the intervals ≤55, 55-65, 65-80, and ≥80 years. The regression model also included the following covariates: sex, systolic blood pressure, total and HDL cholesterol, history of diabetes, antihypertensive medication use, and smoking status. Results: In the ARIC cohort (median enrollment age 54, interquartile range [IQR] 49-59, female 52%), over a median of 24.5 years (IQR 14.4-28.1), there were 1,872 incident CHD events. In comparison to average polygenic risk for CHD (50th percentile of PRS distribution), harboring high polygenic risk (95th percentile) was associated with a 2.9-fold increased risk of incident CHD (95% CI 2.12-3.95) up to age 55 years, independent of conventional risk factors. The strength of association was attenuated later in life (Hazard ratio, age 55-65, 2.5 [2.14-2.84]; age 65-80, 1.87 [1.68-2.07]; age ≥80, 1.78 [1.38-2.31]) Conclusions: A genome-wide PRS for CHD is more strongly associated with incident coronary events at younger age. The strength of association declined later in life while remaining significant. The age-dependent association of PRS with CHD may need to be accounted for in risk prediction models.