Decoding Primary Open-Angle Glaucoma: A Multi-Omics Approach to Identify Druggable Effector Genes

Genomewide association studies (GWAS) have identified numerous primary open angle glaucoma (POAG) risk loci, yet most reside in non-coding regions with unclear function. Mapping these loci to effector genes can elucidate disease mechanisms, identify functionally conserved variants, improve cross-ancestry risk prediction by reducing population-specific noise, and uncover shared therapeutic targets.

Here, we integrate European POAG GWAS with six types of multi-omics molecular Quantitative Trait Locis (xQTLs) using multi-trait colocalization to identify candidate effector variants and evaluate their cross-population relevance using genetic risk score (GRS) analysis, and their therapeutic potential through drug target prioritization.

We identified 25 POAG effector variants colocalized with at least one xQTLs. In non-European populations, effector variants showed stronger effect size correlations with Europeans than non-colocalized variants (Pearson r2 = African 0.85 vs. 0.71; East Asian 0.81 vs. 0.69; and Latin American 0.91 vs. 0.75). Effector variants also had smaller allele frequency variations across populations (average interquartile range [IQR] = 0.15 vs. 0.20). The genetic risk score based on effector variants performed comparably to the genome-wide significant single-nucleotide polymorphism (SNP)-based GRS in non-European populations. Drug prioritization identified zinc, copper, sunitinib, probucol, and astemizole as potential common therapeutic agents for POAG and its subtypes.

Our findings offer deeper insight into the molecular mechanisms underlying glaucoma and effector variants for developing more robust GRS models and broadly effective therapeutic strategies for POAG.