Multi-trait genome-wide association study identifies new loci associated with optic disc parameters

Diagnostic and therapeutic value of human serpin family proteins

SERPIN (serine proteinase inhibitors) is an acronym for the superfamily of structurally similar proteins found in animals, plants, bacteria, viruses, and archaea. Over 1500 SERPINs are known in nature, while only 37 SERPINs are found in humans, which participate in inflammation, coagulation, angiogenesis, cell viability, and other pathophysiological processes. Both qualitative or quantitative deficiencies or overexpression and/or abnormal accumulation of SERPIN can lead to diseases commonly referred to as "serpinopathies". Hence, strategies involving SERPIN supplementation, elimination, or correction are utilized and/or under consideration. In this review, we discuss relationships between certain SERPINs and diseases as well as putative strategies for the clinical explorations of SERPINs.

Human cytomegalovirus infection coopts chromatin organization to diminish TEAD1 transcription factor activity

Human cytomegalovirus (HCMV) infects up to 80% of the world's population. Here, we show that HCMV infection leads to widespread changes in human chromatin accessibility and chromatin looping, with hundreds of thousands of genomic regions affected 48 hours after infection. Integrative analyses reveal HCMV-induced perturbation of Hippo signaling through drastic reduction of TEAD1 transcription factor activity. We confirm extensive concordant loss of TEAD1 binding, active H3K27ac histone marks, and chromatin looping interactions upon infection. Our data position TEAD1 at the top of a hierarchy involving multiple altered important developmental pathways. HCMV infection reduces TEAD1 activity through four distinct mechanisms: closing of TEAD1-bound chromatin, reduction of YAP1 and phosphorylated YAP1 levels, reduction of TEAD1 transcript and protein levels, and alteration of TEAD1 exon-6 usage. Altered TEAD1-based mechanisms are highly enriched at genetic risk loci associated with eye and ear development, providing mechanistic insight into HCMV's established roles in these processes.

A multi-ancestry GWAS of Fuchs corneal dystrophy highlights the contributions of laminins, collagen, and endothelial cell regulation

Fuchs endothelial corneal dystrophy (FECD) is a leading indication for corneal transplantation, but its molecular etiology remains poorly understood. We performed genome-wide association studies (GWAS) of FECD in the Million Veteran Program followed by multi-ancestry meta-analysis with the previous largest FECD GWAS, for a total of 3970 cases and 333,794 controls. We confirm the previous four loci, and identify eight novel loci: SSBP3, THSD7A, LAMB1, PIDD1, RORA, HS3ST3B1, LAMA5, and COL18A1. We further confirm the TCF4 locus in GWAS for admixed African and Hispanic/Latino ancestries and show an enrichment of European-ancestry haplotypes at TCF4 in FECD cases. Among the novel associations are low frequency missense variants in laminin genes LAMA5 and LAMB1 which, together with previously reported LAMC1, form laminin-511 (LM511). AlphaFold 2 protein modeling, validated through homology, suggests that mutations at LAMA5 and LAMB1 may destabilize LM511 by altering inter-domain interactions or extracellular matrix binding. Finally, phenome-wide association scans and colocalization analyses suggest that the TCF4 CTG18.1 trinucleotide repeat expansion leads to dysregulation of ion transport in the corneal endothelium and has pleiotropic effects on renal function.

Assessment of Causality Between Diet-Derived Antioxidants and Primary Open-Angle Glaucoma: A Mendelian Randomization Study

Purpose

This study aimed to investigate the genetic causal relationships among diet-derived circulating antioxidants, primary open-angle glaucoma (POAG), and glaucoma-related traits using two-sample Mendelian randomization (MR).

Methods

Genetic variants associated with diet-derived circulating antioxidants (retinol, ascorbate, β-carotene, lycopene, α-tocopherol, and γ-tocopherol) were assessed as absolute and metabolic instrumental variables. POAG and glaucoma-related traits data were derived from a large, recently published genome-wide association study database; these traits included intraocular pressure (IOP), macular retinal nerve fiber layer (mRNFL) thickness, macular ganglion cell-inner plexiform layer (mGCIPL) thickness, and vertical cup-to-disc ratio (vCDR). MR analyses were performed per outcome for each exposure.

Results

We found no causal association between six diet-derived antioxidants and POAG using the International Glaucoma Genetics Consortium data. For absolute antioxidants, the odds ratios (ORs) ranged from 1.011 (95% confidence interval [CI], 0.854-1.199; P = 0.895) per natural log-transformed β-carotene to 1.052 (95% CI, 0.911-1.215; P = 0.490) for 1 µmol/L of ascorbate. For antioxidant metabolites, the OR ranged from 0.998 (95% CI, 0.801-1.244; P = 0.989) for ascorbate to 1.210 (95% CI, 0.870-1.682; P = 0.257) for γ-tocopherol, using log-transformed levels. A similar result was obtained with the FinnGen Biobank. Furthermore, our results showed no significant genetic association between six diet-derived antioxidants and glaucoma-related traits.

Conclusions

Our study did not support a causal association among six diet-derived circulating antioxidants, POAG, and glaucoma-related traits. This suggests that the intake of antioxidants may not have a preventive effect on POAG and offers no protection to retinal nerve cells.

Translational Relevance

This study provides valid evidence regarding the use of diet-derived antioxidants for glaucoma patients.

WITHDRAWN: Genome-wide risk prediction of primary open-angle glaucoma across multiple ancestries

This manuscript has been withdrawn by medRxiv following a formal request by the QIMR Berghofer Medical Research Institute Research Integrity Office owing to lack of author consent.

Using Drosophila to identify naturally occurring genetic modifiers of amyloid beta 42- and tau-induced toxicity

Alzheimer's disease is characterized by 2 pathological proteins, amyloid beta 42 and tau. The majority of Alzheimer's disease cases in the population are sporadic and late-onset Alzheimer's disease, which exhibits high levels of heritability. While several genetic risk factors for late-onset Alzheimer's disease have been identified and replicated in independent studies, including the ApoE ε4 allele, the great majority of the heritability of late-onset Alzheimer's disease remains unexplained, likely due to the aggregate effects of a very large number of genes with small effect size, as well as to biases in sample collection and statistical approaches. Here, we present an unbiased forward genetic screen in Drosophila looking for naturally occurring modifiers of amyloid beta 42- and tau-induced ommatidial degeneration. Our results identify 14 significant SNPs, which map to 12 potential genes in 8 unique genomic regions. Our hits that are significant after genome-wide correction identify genes involved in neuronal development, signal transduction, and organismal development. Looking more broadly at suggestive hits (P < 10-5), we see significant enrichment in genes associated with neurogenesis, development, and growth as well as significant enrichment in genes whose orthologs have been identified as significantly or suggestively associated with Alzheimer's disease in human GWAS studies. These latter genes include ones whose orthologs are in close proximity to regions in the human genome that are associated with Alzheimer's disease, but where a causal gene has not been identified. Together, our results illustrate the potential for complementary and convergent evidence provided through multitrait GWAS in Drosophila to supplement and inform human studies, helping to identify the remaining heritability and novel modifiers of complex diseases.

Repeat polymorphisms underlie top genetic risk loci for glaucoma and colorectal cancer

Many regions in the human genome vary in length among individuals due to variable numbers of tandem repeats (VNTRs). To assess the phenotypic impact of VNTRs genome-wide, we applied a statistical imputation approach to estimate the lengths of 9,561 autosomal VNTR loci in 418,136 unrelated UK Biobank participants and 838 GTEx participants. Association and statistical fine-mapping analyses identified 58 VNTRs that appeared to influence a complex trait in UK Biobank, 18 of which also appeared to modulate expression or splicing of a nearby gene. Non-coding VNTRs at TMCO1 and EIF3H appeared to generate the largest known contributions of common human genetic variation to risk of glaucoma and colorectal cancer, respectively. Each of these two VNTRs associated with a >2-fold range of risk across individuals. These results reveal a substantial and previously unappreciated role of non-coding VNTRs in human health and gene regulation.

Large-scale multitrait genome-wide association analyses identify hundreds of glaucoma risk loci

Glaucoma, a leading cause of irreversible blindness, is a highly heritable human disease. Previous genome-wide association studies have identified over 100 loci for the most common form, primary open-angle glaucoma. Two key glaucoma-associated traits also show high heritability: intraocular pressure and optic nerve head excavation damage quantified as the vertical cup-to-disc ratio. Here, since much of glaucoma heritability remains unexplained, we conducted a large-scale multitrait genome-wide association study in participants of European ancestry combining primary open-angle glaucoma and its two associated traits (total sample size over 600,000) to substantially improve genetic discovery power (263 loci). We further increased our power by then employing a multiancestry approach, which increased the number of independent risk loci to 312, with the vast majority replicating in a large independent cohort from 23andMe, Inc. (total sample size over 2.8 million; 296 loci replicated at P < 0.05, 240 after Bonferroni correction). Leveraging multiomics datasets, we identified many potential druggable genes, including neuro-protection targets likely to act via the optic nerve, a key advance for glaucoma because all existing drugs only target intraocular pressure. We further used Mendelian randomization and genetic correlation-based approaches to identify novel links to other complex traits, including immune-related diseases such as multiple sclerosis and systemic lupus erythematosus.

Multi-ancestry GWAS of Fuchs corneal dystrophy highlights roles of laminins, collagen, and endothelial cell regulation

Fuchs endothelial corneal dystrophy (FECD) is a leading indication for corneal transplantation, but its molecular pathophysiology remains poorly understood. We performed genome-wide association studies (GWAS) of FECD in the Million Veteran Program (MVP) and meta-analyzed with the previous largest FECD GWAS, finding twelve significant loci (eight novel). We further confirmed the TCF4 locus in admixed African and Hispanic/Latino ancestries, and found an enrichment of European-ancestry haplotypes at TCF4 in FECD cases. Among the novel associations are low frequency missense variants in laminin genes LAMA5 and LAMB1 which, together with previously reported LAMC1, form laminin-511 (LM511). AlphaFold 2 protein modeling suggests that mutations at LAMA5 and LAMB1 may destabilize LM511 by altering inter-domain interactions or extracellular matrix binding. Finally, phenome-wide association scans and co-localization analyses suggest that the TCF4 CTG18.1 trinucleotide repeat expansion leads to dysregulation of ion transport in the corneal endothelium and has pleiotropic effects on renal function.

Mendelian Randomization and Transcriptome-Wide Association Analysis Identified Genes That Were Pleiotropically Associated with Intraocular Pressure

BACKGROUND

Intraocular pressure (IOP) is a major modifiable risk factor for glaucoma. However, the mechanisms underlying the controlling of IOP remain to be elucidated.

OBJECTIVE

To prioritize genes that are pleiotropically associated with IOP.

METHODS

We adopted a two-sample Mendelian randomization method, named summary-based Mendelian randomization (SMR), to examine the pleiotropic effect of gene expression on IOP. The SMR analyses were based on summarized data from a genome-wide association study (GWAS) on IOP. We conducted separate SMR analyses using Genotype-Tissue Expression (GTEx) and Consortium for the Architecture of Gene Expression (CAGE) expression quantitative trait loci (eQTL) data. Additionally, we performed a transcriptome-wide association study (TWAS) to identify genes whose cis-regulated expression levels were associated with IOP.

RESULTS

We identified 19 and 25 genes showing pleiotropic association with IOP using the GTEx and CAGE eQTL data, respectively. RP11-259G18.3 (P_SMR = 2.66 × 10^-6), KANSL1-AS1 (P_SMR = 2.78 × 10^-6), and RP11-259G18.2 (P_SMR = 2.91 × 10^-6) were the top three genes using the GTEx eQTL data. LRRC37A4 (P_SMR = 1.19 × 10^-5), MGC57346 (P_SMR = 1.19 × 10^-5), and RNF167 (P_SMR = 1.53 × 10^-5) were the top three genes using the CAGE eQTL data. Most of the identified genes were found in or near the 17q21.31 genomic region. Additionally, our TWAS analysis identified 18 significant genes whose expression was associated with IOP. Of these, 12 and 4 were also identified by the SMR analysis using the GTEx and CAGE eQTL data, respectively.

CONCLUSIONS

Our findings suggest that the 17q21.31 genomic region may play a critical role in the regulation of IOP.

Comparison of two multi-trait association testing methods and sequence-based fine mapping of six additive QTL in Swiss Large White pigs

BACKGROUND

Genetic correlations between complex traits suggest that pleiotropic variants contribute to trait variation. Genome-wide association studies (GWAS) aim to uncover the genetic underpinnings of traits. Multivariate association testing and the meta-analysis of summary statistics from single-trait GWAS enable detecting variants associated with multiple phenotypes. In this study, we used array-derived genotypes and phenotypes for 24 reproduction, production, and conformation traits to explore differences between the two methods and used imputed sequence variant genotypes to fine-map six quantitative trait loci (QTL).

RESULTS

We considered genotypes at 44,733 SNPs for 5,753 pigs from the Swiss Large White breed that had deregressed breeding values for 24 traits. Single-trait association analyses revealed eleven QTL that affected 15 traits. Multi-trait association testing and the meta-analysis of the single-trait GWAS revealed between 3 and 6 QTL, respectively, in three groups of traits. The multi-trait methods revealed three loci that were not detected in the single-trait GWAS. Four QTL that were identified in the single-trait GWAS, remained undetected in the multi-trait analyses. To pinpoint candidate causal variants for the QTL, we imputed the array-derived genotypes to the sequence level using a sequenced reference panel consisting of 421 pigs. This approach provided genotypes at 16 million imputed sequence variants with a mean accuracy of imputation of 0.94. The fine-mapping of six QTL with imputed sequence variant genotypes revealed four previously proposed causal mutations among the top variants.

CONCLUSIONS

Our findings in a medium-size cohort of pigs suggest that multivariate association testing and the meta-analysis of summary statistics from single-trait GWAS provide very similar results. Although multi-trait association methods provide a useful overview of pleiotropic loci segregating in mapping populations, the investigation of single-trait association studies is still advised, as multi-trait methods may miss QTL that are uncovered in single-trait GWAS.

Mendelian Randomization Implicates Bidirectional Association between Myopia and Primary Open-Angle Glaucoma or Intraocular Pressure

PURPOSE

Observational studies suggest that myopic eyes carry a greater risk of primary open-angle glaucoma (POAG); however, the evidence for this association is inconsistent. This may be the result of confounding factors that arise from myopia that complicate clinical tests for glaucoma. This study used Mendelian randomization (MR) analysis to determine genetic causal associations among myopia, glaucoma, and glaucoma-related traits that overcome the effects of external confounders.

DESIGN

Bidirectional genetic associations between myopia and refractive spherical equivalent (RSE), POAG, and POAG endophenotypes were investigated.

PARTICIPANTS

Data from the largest publicly available genetic banks (n = 216,257-542,934) were analyzed.

METHODS

Multiple MR models and multivariate genomic structural modeling to identify significant mediators for the relationship between myopia and POAG.

MAIN OUTCOME MEASURES

Genetic causal associations between myopia and POAG and POAG endophenotypes.

RESULTS

We found consistent bidirectional genetic associations between myopia and POAG and between myopia and intraocular pressure (IOP) using multiple MR models at Bonferroni-corrected levels of significance. Intraocular pressure showed the most significant mediation effect on RSE and POAG (Sobel test, 0.13; 95% confidence interval, 0.09-0.17; P = 1.37 × 10^-8).

CONCLUSIONS

A strong bidirectional genetic causal link exists between myopia and POAG that is mediated mainly by IOP. Our findings suggest that IOP-lowering treatment for glaucoma may be beneficial in myopic eyes, despite the challenges of establishing a clear clinical diagnosis.

FINANCIAL DISCLOSURE(S)

Proprietary or commercial disclosure may be found after the references.

Development and Evaluation of a Rules-based Algorithm for Primary Open-Angle Glaucoma in the VA Million Veteran Program

The availability of electronic health record (EHR)-linked biobank data for research presents opportunities to better understand complex ocular diseases. Developing accurate computable phenotypes for ocular diseases for which gold standard diagnosis includes imaging remains inaccessible in most biobank-linked EHRs. The objective of this study was to develop and validate a computable phenotype to identify primary open-angle glaucoma (POAG) through accessing the Department of Veterans Affairs (VA) Computerized Patient Record System (CPRS) and Million Veteran Program (MVP) biobank. Accessing CPRS clinical ophthalmology data from VA Medical Center Eye Clinic (VAMCEC) patients, we developed and iteratively refined POAG case and control algorithms based on clinical, prescription, and structured diagnosis data (ICD-CM codes). Refinement was performed via detailed chart review, initially at a single VAMCEC (n = 200) and validated at two additional VAMCECs (n = 100 each). Positive and negative predictive values (PPV, NPV) were computed as the proportion of CPRS patients correctly classified with POAG or without POAG, respectively, by the algorithms, validated by ophthalmologists and optometrists with access to gold-standard clinical diagnosis data. The final algorithms performed better than previously reported approaches in assuring the accuracy and reproducibility of POAG classification (PPV >83% and NPV >97%) with consistent performance in Black or African American and in White Veterans. Applied to the MVP to identify cases and controls, genetic analysis of a known POAG-associated locus further validated the algorithms. We conclude that ours is a viable approach to use combined EHR-genetic data to study patients with complex diseases that require imaging confirmation.

Repurposing Antihypertensive Drugs for the Prevention of Glaucoma: A Mendelian Randomization Study

Purpose

Several antihypertensive drugs have been used for the treatment of glaucoma. However, the effect of hypertension and antihypertensive drugs on glaucoma is still unclear.

Methods

Leveraging large-scale genome-wide association study summary statistics for glaucoma (Ncase = 4737, Ncontrol = 458,196), blood pressure (BP) (N = 422,771), and intraocular pressure (IOP) (N = 31,269), the genetic correlation and causal relationship of genetically assessed IOP, systolic blood pressure (SBP), diastolic blood pressure (DBP), and 12 types of antihypertensive drugs with glaucoma were evaluated using linkage disequilibrium score (LDSC) regression, univariate mendelian randomization (MR), and multivariable MR.

Results

LDSC results showed a suggestive association of glaucoma with SBP (Rg = 0.12, P = 0.0076) and DBP (Rg = 0.17, P = 0.02). In univariate MR, genetically elevated BP in participants was not identified to lead to an increased glaucoma risk (SBP: odds ratio [OR], 1.05 [95% confidence interval {CI}, 0.91-1.21]; P = 0.52; DBP: OR, 1.07 [95% CI, 0.93-1.23]; P = 0.34). The results of univariate MR were replicated in multivariable MR (SBP: OR, 0.95 [95% CI, 0.71-1.29]; P  = 0.75; DBP: OR, 1.13 [95% CI, 0.85-1.51]; P  = 0.41). Furthermore, there was insufficient evidence to suggest that antihypertensive drugs were associated with glaucoma.

Conclusions

Together, controlling BP may not help prevent and treat glaucoma, and antihypertensive drugs may neither treat nor worsen glaucoma.

Translational Relevance

Treating with antihypertensive drugs should not be used as an intervention for patients with glaucoma.

Vision-related convergent gene losses reveal SERPINE3's unknown role in the eye

Despite decades of research, knowledge about the genes that are important for development and function of the mammalian eye and are involved in human eye disorders remains incomplete. During mammalian evolution, mammals that naturally exhibit poor vision or regressive eye phenotypes have independently lost many eye-related genes. This provides an opportunity to predict novel eye-related genes based on specific evolutionary gene loss signatures. Building on these observations, we performed a genome-wide screen across 49 mammals for functionally uncharacterized genes that are preferentially lost in species exhibiting lower visual acuity values. The screen uncovered several genes, including SERPINE3, a putative serine proteinase inhibitor. A detailed investigation of 381 additional mammals revealed that SERPINE3 is independently lost in 18 lineages that typically do not primarily rely on vision, predicting a vision-related function for this gene. To test this, we show that SERPINE3 has the highest expression in eyes of zebrafish and mouse. In the zebrafish retina, serpine3 is expressed in Müller glia cells, a cell type essential for survival and maintenance of the retina. A CRISPR-mediated knockout of serpine3 in zebrafish resulted in alterations in eye shape and defects in retinal layering. Furthermore, two human polymorphisms that are in linkage with SERPINE3 are associated with eye-related traits. Together, these results suggest that SERPINE3 has a role in vertebrate eyes. More generally, by integrating comparative genomics with experiments in model organisms, we show that screens for specific phenotype-associated gene signatures can predict functions of uncharacterized genes.

Glaucoma Heritability: Molecular Mechanisms of Disease

Glaucoma is one of the world’s leading causes of irreversible blindness. A complex, multifactorial disease, the underlying pathogenesis and reasons for disease progression are not fully understood. The most common form of glaucoma, primary open-angle glaucoma (POAG), was traditionally understood to be the result of elevated intraocular pressure (IOP), leading to optic nerve damage and functional vision loss. Recently, researchers have suggested that POAG may have an underlying genetic component. In fact, studies of genetic association and heritability have yielded encouraging results showing that glaucoma may be influenced by genetic factors, and estimates for the heritability of POAG and disease-related endophenotypes show encouraging results. However, the vast majority of the underlying genetic variants and their molecular mechanisms have not been elucidated. Several genes have been suggested to have molecular mechanisms contributing to alterations in key endophenotypes such as IOP (LMX1B, MADD, NR1H3, and SEPT9), and VCDR (ABCA1, ELN, ASAP1, and ATOH7). Still, genetic studies about glaucoma and its molecular mechanisms are limited by the multifactorial nature of the disease and the large number of genes that have been identified to have an association with glaucoma. Therefore, further study into the molecular mechanisms of the disease itself are required for the future development of therapies targeted at genes leading to POAG endophenotypes and, therefore, increased risk of disease.

Automated AI labeling of optic nerve head enables insights into cross-ancestry glaucoma risk and genetic discovery in >280,000 images from UKB and CLSA

Cupping of the optic nerve head, a highly heritable trait, is a hallmark of glaucomatous optic neuropathy. Two key parameters are vertical cup-to-disc ratio (VCDR) and vertical disc diameter (VDD). However, manual assessment often suffers from poor accuracy and is time intensive. Here, we show convolutional neural network models can accurately estimate VCDR and VDD for 282,100 images from both UK Biobank and an independent study (Canadian Longitudinal Study on Aging), enabling cross-ancestry epidemiological studies and new genetic discovery for these optic nerve head parameters. Using the AI approach, we perform a systematic comparison of the distribution of VCDR and VDD and compare these with intraocular pressure and glaucoma diagnoses across various genetically determined ancestries, which provides an explanation for the high rates of normal tension glaucoma in East Asia. We then used the large number of AI gradings to conduct a more powerful genome-wide association study (GWAS) of optic nerve head parameters. Using the AI-based gradings increased estimates of heritability by ∼50% for VCDR and VDD. Our GWAS identified more than 200 loci associated with both VCDR and VDD (double the number of loci from previous studies) and uncovered dozens of biological pathways; many of the loci we discovered also confer risk for glaucoma.

Elevated dimethylarginine, ATP, cytokines, metabolic remodeling involving tryptophan metabolism and potential microglial inflammation characterize primary open angle glaucoma

Glaucoma of which primary open angle glaucoma (POAG) constitutes 75%, is the second leading cause of blindness. Elevated intra ocular pressure and Nitric oxide synthase (NOS) dysfunction are hallmarks of POAG. We analyzed clinical data, cytokine profile, ATP level, metabolomics and GEO datasets to identify features unique to POAG. N9 microglial cells are used to gain mechanistic insights. Our POAG cohort showed elevated ATP in aqueous humor and cytokines in plasma. Metabolomic analysis showed changes in 21 metabolites including Dimethylarginine (DMAG) and activation of tryptophan metabolism in POAG. Analysis of GEO data sets and previously published proteomic data sets bins genes into signaling and metabolic pathways. Pathways from reanalyzed metabolomic data from literature significantly overlapped with those from our POAG data. DMAG modulated purinergic signaling, ATP secretion and cytokine expression were inhibited by N-Ethylmaleimide, NO donors, BAPTA and purinergic receptor inhibitors. ATP induced elevated intracellular calcium level and cytokines expression were inhibited by BAPTA. Metabolomics of cell culture supernatant from ATP treated sets showed metabolic deregulation and activation of tryptophan metabolism. DMAG and ATP induced IDO1/2 and TDO2 were inhibited by N-Ethylmaleimide, sodium nitroprusside and BAPTA. Our data obtained from clinical samples and cell culture studies reveal a strong association of elevated DMAG, ATP, cytokines and activation of tryptophan metabolism with POAG. DMAG mediated ATP signaling, inflammation and metabolic remodeling in microglia might have implications in management of POAG.

Genetic variation affects morphological retinal phenotypes extracted from UK Biobank optical coherence tomography images

Optical Coherence Tomography (OCT) enables non-invasive imaging of the retina and is used to diagnose and manage ophthalmic diseases including glaucoma. We present the first large-scale genome-wide association study of inner retinal morphology using phenotypes derived from OCT images of 31,434 UK Biobank participants. We identify 46 loci associated with thickness of the retinal nerve fibre layer or ganglion cell inner plexiform layer. Only one of these loci has been associated with glaucoma, and despite its clear role as a biomarker for the disease, Mendelian randomisation does not support inner retinal thickness being on the same genetic causal pathway as glaucoma. We extracted overall retinal thickness at the fovea, representative of foveal hypoplasia, with which three of the 46 SNPs were associated. We additionally associate these three loci with visual acuity. In contrast to the Mendelian causes of severe foveal hypoplasia, our results suggest a spectrum of foveal hypoplasia, in part genetically determined, with consequences on visual function.

Genetics of primary open-angle glaucoma and its endophenotypes

Glaucoma is a neurodegenerative disorder characterized by the loss of retinal ganglion cells and optic nerve fibers, resulting in the loss of visual field. Primary open-angle glaucoma (POAG) is the most prevalent subtype of glaucoma. Recent genome-wide association studies (GWASs) identified more than 100 variants associated with POAG and multiple loci associated with endophenotypes including the disc area, vertical cup-to-disc ratio (VCDR), and intraocular pressure (IOP). Especially, several GWASs reported the association between VCDR and variants near CDKN2B/CDKN2B-AS1, ATOH7, and CHEK2, and between IOP and variants near TMCO1, CAV1/CAV2, GAS7, and ARHGEF12. However, the effect of each variant on endophenotypes is modest; therefore, it is useful to construct a genetic risk score (GRS) based on the effect on endophenotypes by combining susceptible genetic variants. Several studies demonstrated that higher GRS was closely associated with endophenotypes including the VCDR, IOP, and age of diagnosis. Henceforth, by quantifying GRS, identification of high risk group before the disease onset, prediction of visual prognosis and early intervention may be possible.