Diversity is key for cross-ancestry transferability of glaucoma genetic risk scores in Hispanic Veterans in the Million Veteran Program.

A major goal of precision medicine is to stratify patients based on their genetic risk for a disease to inform future screening and intervention strategies. For conditions like primary open-angle glaucoma (POAG), the genetic risk architecture is complicated with multiple variants contributing small effects on risk. Following the tepid success of genome-wide association studies for high-effect disease risk variant discovery, genetic risk scores (GRS), which collate effects from multiple genetic variants into a single measure, have shown promise for disease risk stratification. We assessed the application of GRS for POAG risk stratification in Hispanic-descent (HIS) and European-descent (EUR) Veterans in the Million Veteran Program. Unweighted and cross-ancestry meta-weighted GRS were calculated based on 127 genomic variants identified in the most recent report of cross-ancestry POAG meta-analyses. We found that both GRS types were associated with POAG case-control status and performed similarly in HIS and EUR Veterans. This trend was also seen in our subset analysis of HIS Veterans with less than 50% EUR global genetic ancestry. Our findings highlight the importance of evaluating GRS based on known POAG risk variants in different ancestry groups and emphasize the need for more multi-ancestry POAG genetic studies.

Consequences of a Rare Complement Factor H Variant for Age-Related Macular Degeneration in the Amish.

Purpose: Genetic variants in the complement factor H gene (CFH) have been consistently implicated in age-related macular degeneration (AMD) risk. However, their functional effects are not fully characterized. We previously identified a rare, AMD-associated variant in CFH (P503A, rs570523689) in 19 Amish individuals, but its functional consequences were not investigated. Methods: We performed genotyping for CFH P503A in 1326 Amish individuals to identify additional risk allele carriers. We examined differences for age at AMD diagnosis between carriers and noncarriers. In blood samples from risk allele carriers and noncarriers, we quantified (i) CFH RNA expression, (ii) CFH protein expression, and (iii) C-reactive protein (CRP) expression. Potential changes to the CFH protein structure were interrogated computationally with Phyre2 and Chimera software programs. Results: We identified 39 additional carriers from Amish communities in Ohio and Indiana. On average, carriers were younger than noncarriers at AMD diagnosis, but this difference was not significant. CFH transcript and protein levels in blood samples from Amish carriers and noncarriers were also not significantly different. CRP levels were also comparable in plasma samples from carriers and noncarriers. Computational protein modeling showed slight changes in the CFH protein conformation that were predicted to alter interactions between the CFH 503 residue and other neighboring residues. Conclusions: In total, we have identified 58 risk allele carriers for CFH P503A in the Ohio and Indiana Amish. Although we did not detect significant differences in age at AMD diagnosis or expression levels of CFH in blood samples from carriers and noncarriers, we observed modest structural changes to the CFH protein through in silico modeling. Based on our functional and computational observations, we hypothesize that CFH P503A may affect CFH binding or function rather than expression, which would require additional research to confirm.

Glaucoma Genetic Risk Scores in the Million Veteran Program.

PURPOSE: Primary open-angle glaucoma (POAG) is a degenerative eye disease for which early treatment is critical to mitigate visual impairment and irreversible blindness. POAG-associated loci individually confer incremental risk. Genetic risk score(s) (GRS) could enable POAG risk stratification. Despite significantly higher POAG burden among individuals of African ancestry (AFR), GRS are limited in this population. A recent large-scale, multi-ancestry meta-analysis identified 127 POAG-associated loci and calculated cross-ancestry and ancestry-specific effect estimates, including in European ancestry (EUR) and AFR individuals. We assessed the utility of the 127-variant GRS for POAG risk stratification in EUR and AFR Veterans in the Million Veteran Program (MVP). We also explored the association between GRS and documented invasive glaucoma surgery (IGS). DESIGN: Cross-sectional study. PARTICIPANTS: MVP Veterans with imputed genetic data, including 5830 POAG cases (445 with IGS documented in the electronic health record) and 64 476 controls. METHODS: We tested unweighted and weighted GRS of 127 published risk variants in EUR (3382 cases and 58 811 controls) and AFR (2448 cases and 5665 controls) Veterans in the MVP. Weighted GRS were calculated using effect estimates from the most recently published report of cross-ancestry and ancestry-specific meta-analyses. We also evaluated GRS in POAG cases with documented IGS. MAIN OUTCOME MEASURES: Performance of 127-variant GRS in EUR and AFR Veterans for POAG risk stratification and association with documented IGS. RESULTS: GRS were significantly associated with POAG (P < 5 × 10-5) in both groups; a higher proportion of EUR compared with AFR were consistently categorized in the top GRS decile (21.9%-23.6% and 12.9%-14.5%, respectively). Only GRS weighted by ancestry-specific effect estimates were associated with IGS documentation in AFR cases; all GRS types were associated with IGS in EUR cases. CONCLUSIONS: Varied performance of the GRS for POAG risk stratification and documented IGS association in EUR and AFR Veterans highlights (1) the complex risk architecture of POAG, (2) the importance of diverse representation in genomics studies that inform GRS construction and evaluation, and (3) the necessity of expanding diverse POAG-related genomic data so that GRS can equitably aid in screening individuals at high risk of POAG and who may require more aggressive treatment.

Diversity in Polygenic Risk of Primary Open-Angle Glaucoma.

Glaucoma is the leading cause of irreversible blindness worldwide. Primary open-angle glaucoma (POAG), the most common glaucoma subtype, is more prevalent and severe in individuals of African ancestry. Unfortunately, this ancestral group has been historically under-represented among genetic studies of POAG. Moreover, both genetic and polygenic risk scores (GRS, PRS) that are typically based on genetic data from European-descent populations are not transferable to individuals without a majority of European ancestry. Given the aspirations of leveraging genetic information for precision medicine, GRS and PRS demonstrate clinical potential but fall short, in part due to the lack of diversity in these studies. Prioritizing diversity in the discovery of risk variants will improve the performance and utility of GRS and PRS-derived risk estimation for disease stratification, which could bring about earlier POAG intervention and treatment for a disease that often goes undetected until significant damage has occurred.

The GGLEAM Study: Understanding Glaucoma in the Ohio Amish.

Glaucoma leads to millions of cases of visual impairment and blindness around the world. Its susceptibility is shaped by both environmental and genetic risk factors. Although over 120 risk loci have been identified for glaucoma, a large portion of its heritability is still unexplained. Here we describe the foundation of the Genetics of GLaucoma Evaluation in the AMish (GGLEAM) study to investigate the genetic architecture of glaucoma in the Ohio Amish, which exhibits lower genetic and environmental heterogeneity compared to the general population. To date, we have enrolled 81 Amish individuals in our study from Holmes County, Ohio. As a part of our enrollment process, 62 GGLEAM study participants (42 glaucoma-affected and 20 unaffected individuals) received comprehensive eye examinations and glaucoma evaluations. Using the data from the Anabaptist Genealogy Database, we found that 80 of the GGLEAM study participants were related to one another through a large, multigenerational pedigree containing 1586 people. We plan to integrate the health and kinship data obtained for the GGLEAM study to interrogate glaucoma genetics and pathophysiology in this unique population.

Statistical driver genes as a means to uncover missing heritability for age-related macular degeneration.

BACKGROUND: Age-related macular degeneration (AMD) is a progressive retinal disease contributing to blindness worldwide. Multiple estimates for AMD heritability (h2) exist; however, a substantial proportion of h2 is not attributable to known genomic loci. The International AMD Genomics Consortium (IAMDGC) gathered the largest dataset of advanced AMD (ADV) cases and controls available and identified 34 loci containing 52 independent risk variants defining known AMD h2. To better define AMD heterogeneity, we used Pathway Analysis by Randomization Incorporating Structure (PARIS) on the IAMDGC data and identified 8 statistical driver genes (SDGs), including 2 novel SDGs not discovered by the IAMDGC. We chose to further investigate these pathway-based risk genes and determine their contribution to ADV h2, as well as the differential ADV subtype h2. METHODS: We performed genomic-relatedness-based restricted maximum-likelihood (GREML) analyses on ADV, geographic atrophy (GA), and choroidal neovascularization (CNV) subtypes to investigate the h2 of genotyped variants on the full DNA array chip, 34 risk loci (n = 2758 common variants), 52 variants from the IAMDGC 2016 GWAS, and the 8 SDGs, specifically the novel 2 SDGs, PPARA and PLCG2. RESULTS: Via GREML, full chip h2 was 44.05% for ADV, 46.37% for GA, and 62.03% for CNV. The lead 52 variants' h2 (ADV: 14.52%, GA: 8.02%, CNV: 13.62%) and 34 loci h2 (ADV: 13.73%, GA: 8.81%, CNV: 12.89%) indicate that known variants contribute ~ 14% to ADV h2. SDG variants account for a small percentage of ADV, GA, and CNV heritability, but estimates based on the combination of SDGs and the 34 known loci are similar to those calculated for known loci alone. We identified modest epistatic interactions among variants in the 2 SDGs and the 52 IAMDGC variants, including modest interactions between variants in PPARA and PLCG2. CONCLUSIONS: Pathway analyses, which leverage biological relationships among genes in a pathway, may be useful in identifying additional loci that contribute to the heritability of complex disorders in a non-additive manner. Heritability analyses of these loci, especially amongst disease subtypes, may provide clues to the importance of specific genes to the genetic architecture of AMD.

Pathway Analysis Integrating Genome-Wide and Functional Data Identifies PLCG2 as a Candidate Gene for Age-Related Macular Degeneration.

Purpose: Age-related macular degeneration (AMD) is the worldwide leading cause of blindness among the elderly. Although genome-wide association studies (GWAS) have identified AMD risk variants, their roles in disease etiology are not well-characterized, and they only explain a portion of AMD heritability. Methods: We performed pathway analyses using summary statistics from the International AMD Genomics Consortium's 2016 GWAS and multiple pathway databases to identify biological pathways wherein genetic association signals for AMD may be aggregating. We determined which genes contributed most to significant pathway signals across the databases. We characterized these genes by constructing protein-protein interaction networks and performing motif analysis. Results: We determined that eight genes (C2, C3, LIPC, MICA, NOTCH4, PLCG2, PPARA, and RAD51B) "drive" the statistical signals observed across pathways curated in the Kyoto Encyclopedia of Genes and Genomes (KEGG), Reactome, and Gene Ontology (GO) databases. We further refined our definition of statistical driver gene to identify PLCG2 as a candidate gene for AMD due to its significant gene-level signals (P < 0.0001) across KEGG, Reactome, GO, and NetPath pathways. Conclusions: We performed pathway analyses on the largest available collection of advanced AMD cases and controls in the world. Eight genes strongly contributed to significant pathways from the three larger databases, and one gene (PLCG2) was central to significant pathways from all four databases. This is, to our knowledge, the first study to identify PLCG2 as a candidate gene for AMD based solely on genetic burden. Our findings reinforce the utility of integrating in silico genetic and biological pathway data to investigate the genetic architecture of AMD.

Rare variants and loci for age-related macular degeneration in the Ohio and Indiana Amish.

Age-related macular degeneration (AMD) is a leading cause of blindness in the world. While dozens of independent genomic variants are associated with AMD, about one-third of AMD heritability is still unexplained. To identify novel variants and loci for AMD, we analyzed Illumina HumanExome chip data from 87 Amish individuals with early or late AMD, 79 unaffected Amish individuals, and 15 related Amish individuals with unknown AMD affection status. We retained 37,428 polymorphic autosomal variants across 175 samples for association and linkage analyses. After correcting for multiple testing (n = 37,428), we identified four variants significantly associated with AMD: rs200437673 (LCN9, p = 1.50 × 10-11), rs151214675 (RTEL1, p = 3.18 × 10-8), rs140250387 (DLGAP1, p = 4.49 × 10-7), and rs115333865 (CGRRF1, p = 1.05 × 10-6). These variants have not been previously associated with AMD and are not in linkage disequilibrium with the 52 known AMD-associated variants reported by the International AMD Genomics Consortium based on physical distance. Genome-wide significant linkage peaks were observed on chromosomes 8q21.11-q21.13 (maximum recessive HLOD = 4.03) and 18q21.2-21.32 (maximum dominant HLOD = 3.87; maximum recessive HLOD = 4.27). These loci do not overlap with loci previously linked to AMD. Through gene ontology enrichment analysis with ClueGO in Cytoscape, we determined that several genes in the 1-HLOD support interval of the chromosome 8 locus are involved in fatty acid binding and triglyceride catabolic processes, and the 1-HLOD support interval of the linkage region on chromosome 18 is enriched in genes that participate in serine-type endopeptidase inhibitor activity and the positive regulation of epithelial to mesenchymal transition. These results nominate novel variants and loci for AMD that require further investigation.