Pathway-Specific Polygenic Risk Scores Identify Obstructive Sleep Apnea–Related Pathways Differentially Moderating Genetic Susceptibility to Coronary Artery Disease

Abstract

Background: Obstructive sleep apnea (OSA) and its features, such as chronic intermittent hypoxia, may differentially affect specific molecular pathways and processes in the pathogenesis of coronary artery disease (CAD) and influence the subsequent risk and severity of CAD events. In particular, competing adverse (eg, inflammatory) and protective (eg, increased coronary collateral blood flow) mechanisms may operate, but remain poorly understood. We hypothesize that common genetic variation in selected molecular pathways influences the likelihood of CAD events differently in individuals with and without OSA, in a pathway-dependent manner. Methods: We selected a cross-sectional sample of 471 877 participants from the UK Biobank, with 4974 ascertained to have OSA, 25 988 to have CAD, and 711 to have both. We calculated pathway-specific polygenic risk scores for CAD, based on 6.6 million common variants evaluated in the CARDIoGRAMplusC4D genome-wide association study (Coronary ARtery DIsease Genome wide Replication and Meta-analysis [CARDIoGRAM] plus The Coronary Artery Disease [C4D] Genetics), annotated to specific genes and pathways using functional genomics databases. Based on prior evidence of involvement with intermittent hypoxia and CAD, we tested pathway-specific polygenic risk scores for the HIF1 (hypoxia-inducible factor 1), VEGF (vascular endothelial growth factor), NFκB (nuclear factor kappa-light-chain-enhancer of activated B cells) and TNF (tumor necrosis factor) signaling pathways. Results: In a multivariable-adjusted logistic generalized additive model, elevated pathway-specific polygenic risk scores for the Kyoto Encyclopedia of Genes and Genomes VEGF pathway (39 genes) associated with protection for CAD in OSA (interaction odds ratio 0.86, P= 6×10 -4 ). By contrast, the genome-wide CAD PRS did not show evidence of statistical interaction with OSA. Conclusions: We find evidence that pathway-specific genetic risk of CAD differs between individuals with and without OSA in a qualitatively pathway-dependent manner. These results provide evidence that gene-by-environment interaction influences CAD risk in certain pathways among people with OSA, an effect that is not well-captured by the genome-wide PRS. This invites further study of how OSA interacts with genetic risk at the molecular level and suggests eventual personalization of OSA treatment to reduce CAD risk according to individual pathway-specific genetic risk profiles.