The Clinical Usefulness of a Glaucoma Polygenic Risk Score in 4 Population-Based European Ancestry Cohorts

Abstract

ObjectiveWe used a polygenic risk score (PRS) to identify high-risk groups for primary open-angle glaucoma (POAG) within population-based cohorts. DesignSecondary analysis of four prospective population-based studies. ParticipantsWe included four European-ancestry cohorts: the United States (US) based Nurses’ Health Study (NHS), NHS2, and the Health Professionals Follow-up Study, and the Rotterdam Study (RS) in the Netherlands. The US cohorts included female nurses and male health professionals aged 55+ years. The RS included residents aged 45 years or older living in Rotterdam, the Netherlands. MethodsPRS weights were estimated by applying the Lassosum method on imputed genotype and phenotype data from the UK-Biobank. This resulted in 144,020 variants, single nucleotide polymorphism (SNPs) and indels, with non-zero betas that we used to calculate a PRS in the target populations. Using multivariable Cox proportional hazard models, we estimated the relationship between the standardized PRS and relative risk for POAG. Additionally, POAG prediction was tested by calculating these models’ concordance (Harrell’s C-statistic). Finally, we assessed the association between PRS tertiles and glaucoma-related traits. Main Outcome MeasuresThe relative risk for POAG and Harrell’s C-statistic (the equivalent of an area-under-the-curve for longitudinal models). ResultsAmong 1,046 cases and 38,809‬ controls, the relative risk (95% confidence interval) for POAG for participants in the highest PRS quintile was 3.99 (3.08, 5.18) in the US cohorts, and 4.89 (2.93, 8.17) in the Rotterdam Study, compared with participants with median genetic risk (3rd quintile). In restricted cubic spline analyses, the relation between continuous PRS and POAG risk increased exponentially in the US and Rotterdam cohorts (Pspline<0.05). Combining age, sex, intraocular pressure >25 mmHg, and family history resulted in a meta-analyzed concordance of 0.75 (0.73, 0.75). Adding the PRS to this model improved the concordance to 0.82 (0.80, 0.84). In a meta-analysis of all cohorts, cases in the highest tertile had a larger cup-disc ratio at diagnosis, by 0.11 (0.07, 0.15), and a 2.07-fold increased risk of requiring glaucoma surgery (1.19, 3.60). ConclusionsIncorporating a PRS into a POAG predictive model improves identification concordance from 0.75 up to 0.82, supporting its potential for guiding more cost-effective screening strategies.