Cross-ancestry genome-wide association analysis of corneal thickness strengthens link between complex and Mendelian eye diseases

Pilot work of the 10K Chinese People Genomic Diversity Project along the Silk Road suggests a complex east-west admixture landscape and biological adaptations

Genomic sources from China are underrepresented in the population-specific reference database. We performed whole-genome sequencing or genome-wide genotyping on 1,207 individuals from four linguistically diverse groups (1,081 Sinitic, 56 Mongolic, 40 Turkic, and 30 Tibeto-Burman people) living in North China included in the 10K Chinese People Genomic Diversity Project (10K_CPGDP) to characterize the genetic architecture and adaptative history of ethnic groups in the Silk Road Region of China. We observed a population split between Northwest Chinese minorities (NWCMs) and Han Chinese since the Upper Paleolithic and later Neolithic genetic differentiation within NWCMs. The observed population substructures among ethnically/linguistically diverse NWCMs suggested that differentiated admixture events contributed to the differences in their genomic and phenotypic diversity. We estimated that the Dongxiang, Tibetan, and Yugur people inherited more than 10% of the Western Eurasian ancestry, which is much greater than that of the Salar and Tu people (<7%), while Han neighbors showed less West Eurasian ancestry (∼1%-3%). Male-biased admixture introduced Western Eurasian ancestry in the Dongxiang, Tibetan, and Yugur populations. We found that the eastern-western admixture in NWCMs occurred ∼800-1,100 years ago, coinciding with intensive economic and cultural exchanges during the Tang and Song dynasties. Additionally, we identified the signatures of natural selection associated with cardiovascular system diseases or lipid metabolism and developmental/neurogenetic disorders. Moreover, the EPAS1 gene showed relatively high population branch statistic values in NWCMs. The well-fitted demographical models presented the vast landscape of complex admixture processes of the Silk Road people, and the newly reported functionally important variations suggested the importance of including ethnolinguistically diverse populations in Chinese genetic studies for uncovering the genetic basis of complex traits/diseases.

Mussel-Inspired In Situ Photodynamic Antibacterial Coating for Postoperative Management of Artificial Corneas

Artificial corneal (AC) implants offer hope to millions with corneal blindness, including 5 million in China. Titanium is one of the materials commonly used in the fabrication of artificial corneal scaffolders because of its light texture, high mechanical properties, and high biosafety. However, postoperative bacterial infections, especially from Pseudomonas aeruginosa and Staphylococcus aureus (S. aureus), remain a significant challenge due to the bioinert nature of titanium materials, leading to high infection rates. In this study, we introduce an innovative in situ photodynamic coating technology designed to manage postoperative infections in artificial corneas. Inspired by mussel adhesive proteins, this coating employs a composite of APTES-TA formed by Schiff base and Michael addition reactions between 3-aminopropyltriethoxysilane (APTES) and tannic acid (TA), integrated with the bacterial targeting capabilities of 4-carboxyphenylboronic acid (CPBA) and the photo-Fenton activity of FeOOH (iron(III) hydroxide). The design of the AC@APTES-TA-CPBA-FeOOH coating leverages the dynamic boronate ester bonds, which interact specifically with bacteria in tears, effectively capturing them on the surface of the artificial cornea. The coating exhibits a photocatalytic Fenton-like effect, which confers it with an exceptional bactericidal efficiency of over 99% in vitro. Furthermore, it demonstrates excellent protective functionality for mouse corneas in vivo experiments.

Developing and validating a comprehensive polygenic risk score to enhance keratoconus risk prediction

PURPOSE

This study aimed to develop and validate a comprehensive polygenic risk score (PRS) for keratoconus, enhancing the predictive accuracy for identifying individuals at increased risk, which is crucial for preventing keratoconus-associated visual impairment such as post-Laser-assisted in situ keratomileusis (LASIK) ectasia.

METHODS

We applied a multi-trait analysis approach (MTAG) to genome-wide association study data on keratoconus and quantitative keratoconus-related traits and used this to construct PRS models for keratoconus risk using several PRS methodologies. We evaluated the predictive performance of the PRSs in two biobanks: Estonian Biobank (EstBB; 375 keratoconus cases and 17 902 controls) and UK Biobank (UKB: 34 keratoconus cases and 1000 controls). Scores were compared using the area under the curve (AUC) and odds ratios (ORs) for various PRS models.

RESULTS

The PRS models demonstrated significant predictive capabilities in EstBB, with the SBayesRC model achieving the highest OR of 2.28 per standard deviation increase in PRS, with a model containing age, sex and PRS showing good predictive accuracy (AUC = 0.72). In UKB, we found that adding the best-performing PRS to a model containing corneal measurements increased the AUC from 0.84 to 0.88 (P = 0.012 for difference), with an OR of 4.26 per standard deviation increase in the PRS. These models showed improved predictive capability compared to previous keratoconus PRS.

CONCLUSION

The PRS models enhanced prediction of keratoconus risk, even with corneal measurements, showing potential for clinical use to identify individuals at high risk of keratoconus, and potentially help reduce the risk of post-LASIK ectasia.

Whole-genome sequencing identifies novel loci for keratoconus and facilitates risk stratification in a Han Chinese population

BACKGROUND

Keratoconus (KC) is a prevalent corneal condition with a modest genetic basis. Recent studies have reported significant genetic associations in multi-ethnic cohorts. However, the situation in the Chinese population remains unknown. This study was conducted to identify novel genetic variants linked to KC and to evaluate the potential applicability of a polygenic risk model in the Han Chinese population.

METHODS

A total of 830 individuals diagnosed with KC and 779 controls from a Chinese cohort were enrolled and genotyped by whole-genome sequencing (WGS). Common and rare variants were respectively subjected to single variant association analysis and gene-based burden analysis. Polygenic risk score (PRS) models were developed using top single-nucleotide polymorphisms (SNPs) identified from a multi-ethnic meta-analysis and then evaluated in the Chinese cohort.

RESULTS

The characterization of germline variants entailed correction for population stratification and validation of the East Asian ancestry of the included samples via principal component analysis. For rare protein-truncating variants (PTVs) with minor allele frequency (MAF) < 5%, ZC3H11B emerged as the top prioritized gene, albeit failing to reach the significance threshold. We detected three common variants reaching genome-wide significance (P ≤ 5 × 10-8), all of which are novel to KC. Our study validated three well known predisposition loci, COL5A1, EIF3A and FNDC3B. Additionally, a significant correlation of allelic effects was observed for suggestive SNPs between the largest multi-ethnic meta-genome-wide association study (GWAS) and our study. The PRS model, generated using top SNPs from the meta-GWAS, stratified individuals in the upper quartile, revealing up to a 2.16-fold increased risk for KC.

CONCLUSIONS

Our comprehensive WGS-based GWAS in a large Chinese cohort enhances the efficiency of array-based genetic studies, revealing novel genetic associations for KC and highlighting the potential for refining clinical decision-making and early prevention strategies.

Mechanical Strain of Corneal Epithelium Influences the Expression of Genes Implicated in Keratoconus

Purpose

Although mechanical injury to the cornea (e.g. chronic eye rubbing) is a known risk factor for keratoconus progression, how it contributes to loss of corneal integrity is not known. Here, we set out to determine how eye rubbing can influence keratoconus progression by exploring the expression of known disease markers in mechanically stressed corneal epithelial cells.

Methods

To explore the effects of mechanical stress on the expression of genes implicated in keratoconus (e.g. WNT10A, COL12A1, and TGFB1), we measured their expression using an in vitro model that simulates eye rubbing by cyclic stretching of an immortalized human corneal epithelial cell line (hTCEpi) for 16 hours. We further examined the influence of WNT10A expression in hTCEpi cells using loss-of-function approaches.

Results

Mechanical strain led to a marked reduction in WNT10A mRNA and protein expression, as well as decreased collagen XII mRNA and protein expression, in hTCEpi cells. Reduced expression of WNT10A protein in WNT10A knockdown cells resulted in reduced protein expression of collagens I and XII, and reduced mRNA expression of MMP9 and TGFB1. Conversely, primary keratocytes treated with recombinant WNT10A protein increased TGFB1 mRNA expression.

Conclusions

We provide a molecular explanation for how mechanical strain results in reduced expression of WNT10A in the corneal epithelium, which, in turn, leads to depletion of collagen type I and XII, and TGFβ1 expression. These results provide a molecular link among mechanical strain, WNT10A expression, and the biomechanical failure of the keratoconus cornea.

Genetic Variants Associated With a History of Glaucoma Surgery in Japanese Patients With Primary Open Angle Glaucoma

PRCIS

The genetic risk score (GRS) of genetic variants associated with intraocular pressure (IOP) elevation, but not those associated with optic nerve vulnerability, was associated with a history of glaucoma surgery in Japanese patients with primary open angle glaucoma (POAG).

OBJECTIVE

To investigate genetic variants associated with a history of glaucoma surgery in Japanese patients with POAG.

METHODS

Japanese patients with POAG (n = 468), including normal tension glaucoma (n = 246) and high tension glaucoma (n = 222), and control subjects (n = 246) were genotyped for 22 genetic variants predisposing to POAG, which can be classified into those associated with IOP elevation (IOP-related genetic variants) and optic nerve vulnerability independent of IOP (optic nerve-related genetic variants). The unweighted and weighted GRSs of 17 IOP-related, 5 optic nerve-related, and all 22 genetic variants were calculated, and the association between the GRS and a history of glaucoma surgery was evaluated.

RESULTS

There was a significant association (odds ratio 1.13 per unweighted GRS, 95% CI: 1.03 to 1.24, P = 0.0093) between IOP-related unweighted GRS and a history of glaucoma surgery. A significant association (odds ratio 1.09 per 0.1 weighted GRS, 95% CI: 1.04 to 1.14, P = 0.00022) was also found between IOP-related weighted GRS and a history of glaucoma surgery. The IOP-related GRS was positively correlated with the need for glaucoma surgery. The mean of IOP-related unweighted and weighted GRS in patients with POAG with a history of glaucoma surgery were significantly higher ( P = 0.013 and P = 0.00031, respectively) than those in patients with POAG without a history of glaucoma surgery.

CONCLUSIONS

IOP-related, but not optic nerve-related, genetic variants were associated with a history of glaucoma surgery in Japanese patients with POAG. These results indicate that IOP elevation induced by IOP-related genetic variants rather than optic nerve vulnerability induced by optic nerve-related genetic variants may play an important role in requiring glaucoma surgery.

Central corneal thickness and the risk of primary open-angle glaucoma: a Mendelian randomisation mediation analysis

BACKGROUND

The association of central corneal thickness (CCT) with primary open-angle glaucoma (POAG) remains uncertain. Although several observational studies assessing this relationship have reported an inverse association between CCT and POAG, this could be the result of collider bias. In this study, we leveraged human genetic data to assess through Mendelian randomisation (MR) the effect of CCT on POAG risk and whether this effect is mediated by intraocular pressure (IOP) changes.

METHODS

We used 24 single-nucleotide polymorphisms (SNPs) associated with CCT (p value<5×10-8) from a genome-wide association study (GWAS) (N=17 803) provided by the International Glaucoma Genetics Consortium and 53 SNPs associated with IOP (p value<5×10-8) from a GWAS of the UK Biobank (UKBB) (N=97 653). We related these instruments to POAG using a GWAS meta-analysis of 8283 POAG cases and 753 827 controls from UKBB and FinnGen.

RESULTS

MR analysis suggested a positive association between CCT and POAG (OR of POAG per 50 µm increase in CCT: 1.38; 95% CI: 1.18 to 1.61; p value<0.01). MR mediation analysis showed that 28.4% of the total effect of CCT on POAG risk was mediated through changes in IOP. The primary results were consistent with estimates of pleiotropy-robust MR methods.

CONCLUSION

Contrary to most observational studies, our results showed that a higher CCT is associated with an increased risk of POAG.

In Search of Mouse Models for Exfoliation Syndrome

PURPOSE

Exfoliation syndrome (XFS) is a systemic connective tissue disorder with elusive pathophysiology. We hypothesize that a mouse model with elastic fiber defects caused by lack of lysyl oxidase like 1 (LOXL1 encoded by Loxl1), combined with microfibril deficiency due to Fbn1 mutation (encoding fibrillin-1, Fbn1C1041G/+) will display ocular and systemic phenotypes of XFS.

METHODS

Loxl1-/- was crossed with Fbn1C1041G/+ to create double mutant (dbm) mice. Intraocular pressure (IOP), visual acuity (VA), electroretinogram (ERG), and biometry were characterized in 4 genotypes (wt, Fbn1C1041G/+, Loxl1-/-, dbm) at 16 weeks of age. Optic nerve (ON) area was measured by ImageJ, and axon counting was achieved by AxonJ. Deep whole-body phenotyping was performed in wt and dbm mice. Two-tailed Student t test was used for statistical analysis.

RESULTS

There was no difference in IOP between the 4 genotypes. VA was significantly reduced only in dbm mice. The majority of biometric parameters showed significant differences in all 3 mutant genotypes compared with wt, and dbm had exacerbated anomalies compared with single mutants. Dbm mice showed reduced retinal function and significantly enlarged ON area compared with wt. Dbm mice exhibited severe systemic phenotypes related to abnormal elastic fibers, such as pelvic organ prolapse and cardiovascular and pulmonary abnormalities.

CONCLUSIONS

Ocular and systemic findings in dbm mice support functional overlap between fibrillin-1 and LOXL1, 2 prominent components of exfoliation material. Although no elevated IOP or reduction of axon numbers was detected in dbm mice at 16 weeks of age, their reduced retinal function and enlarged ON area indicate early retinal ganglion cell dysfunction. Dbm mice also provide insight on the link between XFS and systemic diseases in humans. NOTE: Publication of this article is sponsored by the American Ophthalmological Society.

Comprehensive Evaluation of the Genetic Basis of Keratoconus: New Perspectives for Clinical Translation

Purpose

Keratoconus (KC) is a corneal disorder with complex etiology, apparently involving both genetic and environmental factors, characterized by progressive thinning and protrusion of the cornea. We aimed to identify novel genetic regions associated with KC susceptibility, elucidate relevant genes for disease development, and explore the translational implications for therapeutic intervention and risk assessment.

Methods

We conducted a genome-wide association study (GWAS) that integrated previously published data with newly generated genotyping data from an independent European cohort. To evaluate the clinical translation of our results, we performed functional annotation, gene prioritization, polygenic risk score (PRS), and drug repositioning analyses.

Results

We identified two novel genetic loci associated with KC, with rs2806689 and rs807037 emerging as lead variants (P = 1.71E-08, odds ratio [OR] = 0.88; P = 1.93E-08, OR = 1.16, respectively). Most importantly, we identified 315 candidate genes influenced by confirmed KC-associated variants. Among these, MINK1 was found to play a pivotal role in KC pathogenesis through the WNT signaling pathway. Moreover, we developed a PRS model that successfully differentiated KC patients from controls (P = 7.61E-16; area under the curve = 0.713). This model has the potential to identify individuals at high risk for developing KC, which could be instrumental in early diagnosis and management. Additionally, our drug repositioning analysis identified acetylcysteine as a potential treatment option for KC, opening up new avenues for therapeutic intervention.

Conclusions

Our study provides valuable insights into the genetic and molecular basis of KC, offering new targets for therapy and highlighting the clinical utility of PRS models in predicting disease risk.

Polygenic Prediction of Keratoconus and its Measures: Cross-Sectional and Longitudinal Analyses in Community-Based Young Adults

PURPOSE

This study evaluates the performance of a multitrait polygenic risk score (PRS) in an independent cohort to predict incident or progression of keratoconus.

DESIGN

Prospective cross-sectional and cohort study METHODS: Setting: Single-center; Study population: 1478 community-based young adults (18-30 years; 51% female), including 609 (52% female) who returned for an 8-year follow-up; Observation procedures: Scheimpflug imaging (Pentacam, Oculus), genotyping and development of a multitrait PRS previously validated to predict keratoconus in older adults.; Main outcome measure: Belin/Ambrόsio enhanced ectasia display (BAD-D) score and keratoconus, defined as BAD-D ≥2.6, were each analyzed against the PRS using linear and logistic regression, respectively.

RESULTS

Prevalence of keratoconus was 2.5% (95% confidence interval [CI] = 1.9-3.6) in the cross-sectional cohort. Each z-score increase in PRS was associated with worse BAD-D z-score by 0.13 (95%CI = 0.08-0.18) and 1.6 increased odds of keratoconus. The 8-year keratoconus incidence was 2.6% (95%CI = 1.3-4.0). Participants in the highest PRS decile were more likely to have incident keratoconus compared to the rest of the cohort (odds ratio = 3.85, 95%CI = 1.21-12.22). For each z-score increase in PRS, 8-year change in BAD-D z-score worsened by 0.11 (95%CI = 0.04-0.17).

CONCLUSIONS

A PRS for keratoconus could be useful in predicting incident keratoconus and progression, demonstrating its potential utility in clinical settings to identify patients at high risk of postsurgery ectasia or those who may benefit most from keratoconus intervention.

Causal relationship between central corneal thickness and open-angle glaucoma: Evidence from Mendelian randomization

Glaucoma is the leading cause of irreversible blindness worldwide. Previous observational studies have suggested a relationship between central corneal thickness (CCT) and glaucoma; however, the results are inconsistent. This study aimed to investigate whether CCT is associated with a risk for developing open-angle glaucoma (OAG). We employed two-sample Mendelian randomization to assess the relationship between CCT and OAG, namely, primary open-angle glaucoma (POAG) and suspected glaucoma. Genetic instruments composed of variants associated with CCT at genome-wide significance (P < 5 × 10-8) were obtained from published genome-wide association studies from Iglesias et al. for discovery and Bonnemaijer et al. for replication. Summary-level statistics for these instruments for the OAG were obtained from the FinnGen Project (Release 10). Inverse-variance-weighted regression of genetic susceptibility predicted that increased CCT was positively associated with an increased risk for POAG (odds ratio [OR], 1.005; 95% confidence interval [CI], 1.002-1.008; P = 0.001) and suspected glaucoma (OR, 1.006; 95% CI, 1.003-1.009; P < 0.001). In the replication sample of CCT, increased CCT was also positively associated with an increased risk for POAG (OR, 1.004; 95% CI, 1.000-1.008; P = 0.029) and suspected glaucoma (OR, 1.005; 95% CI, 1.001-1.008; P = 0.013). We found genetic evidence supporting a potential causal association between increased CCT and the risk of POAG and suspected glaucoma in the European population. This findings indicates the clinical significance of CCT in the diagnosis and treatment of glaucoma. Further studies are needed to elucidate the underlying mechanisms of this causal relationship.

Genome-wide association study of the common retinal disorder epiretinal membrane: Significant risk loci in each of three American populations

Epiretinal membrane (ERM) is a common retinal condition characterized by the presence of fibrocellular tissue on the retinal surface, often with visual distortion and loss of visual acuity. We studied European American (EUR), African American (AFR), and Latino (admixed American, AMR) ERM participants in the Million Veteran Program (MVP) for genome-wide association analysis-a total of 38,232 case individuals and 557,988 control individuals. We completed a genome-wide association study (GWAS) in each population separately, and then results were meta-analyzed. Genome-wide significant (GWS) associations were observed in all three populations studied: 31 risk loci in EUR subjects, 3 in AFR, and 2 in AMR, with 48 in trans-ancestry meta-analysis. Many results replicated in the FinnGen sample. Several GWS variants associate to alterations in gene expression in the macula. ERM showed significant genetic correlation to multiple traits. Pathway enrichment analyses implicated collagen and collagen-adjacent mechanisms, among others. This well-powered ERM GWAS identified novel genetic associations that point to biological mechanisms for ERM.

Identification of Keratoconus-Related Phenotypes in Three Ppip5k2 Mouse Models

Purpose

It is necessary to establish a mouse model of keratoconus (KC) for research and therapy. We aimed to determine corneal phenotypes in 3 Ppip5k2 mouse models.

Methods

Central corneal thickness (CCT) was determined using spectral domain optical coherence tomography (SD-OCT) in Ppip5k2+/K^ (n = 41 eyes), Ppip5k2K^/K^ (n = 17 eyes) and 2 knock-in mice, Ppip5k2S419A/+ (n = 54 eyes) and Ppip5k2S419A/S419A (n = 18 eyes), and Ppip5k2D843S/+ (n = 42 eyes) and Ppip5k2D843S/D843S (n = 44 eyes) at 3 and 6 months. Pachymetry maps were generated using the Mouse Corneal Analysis Program (MCAP) to process OCT images. Slit lamp biomicroscopy was used to determine any corneal abnormalities, and, last, hematoxylin and eosin (H&E) staining using corneal sections from these animals was used to examine morphological changes.

Results

CCT significantly decreased from 3 to 6 months in the Ppip5k2+/K^ and Ppip5k2K^/K^ mice compared to their littermate controls. OCT-based pachymetry maps revealed abnormally localized thinning in all three models compared to their wild-type (WT) controls. Slit lamp examinations revealed corneal abnormalities in the form of bullous keratopathy, stromal edema, stromal scarring, deep corneal neovascularization, and opacities in the heterozygous/homozygous mice of the three models in comparison with their controls. Corneal histological abnormalities, such as epithelial thickening and stromal layer damage, were observed in the heterozygous/homozygous mice of the three models in comparison with the WT controls.

Conclusions

We have identified phenotypic and histological changes in the corneas of three mouse lines that could be relevant in the development of animal models of KC.

High phosphate and calcium induce osteoblastic phenotype switching and calcification of corneal epithelial cells in a Runx2-dependent and synergistic manner; a possible mechanism of chronic kidney disease-associated corneal calcification

Patients with advanced chronic kidney disease (CKD) have elevated circulating calcium × phosphate product levels and exhibit soft tissue calcification. Besides the cardiovascular system, calcification is commonly observed in the cornea in CKD patients on hemodialysis. Cardiovascular calcification is a cell-mediated, highly regulated process, and we hypothesized that a similar regulatory mechanism is implicated in corneal calcification with the involvement of corneal epithelial cells (CECs). We established a mouse model of CKD-associated corneal calcification by inducing CKD in DBA/2J mice with an adenine and high phosphate diet. CKD was associated with aorta and corneal calcification as detected by OsteoSense staining and corneal Ca measurement (1.67-fold elevation, p < 0.001). In vitro, excess phosphate and Ca induced human CEC calcification in a dose-dependent and synergistic manner, without any influence on cell viability. High phosphate and Ca-containing osteogenic medium (OM; 2.5 mmol/L excess phosphate and 0.6 mmol/L excess Ca over control) increased the protein expression of Runx2 and induced its nuclear translocation. OM increased the expression of the bone-specific Ca-binding protein osteocalcin (130-fold increase, p < 0.001). Silencing of Runx2 attenuated OM-induced CEC calcification. Immunohistology revealed upregulation of Runx2 and overlapping between the Runx2 and the Alizarin red positive areas of calcification in the cornea of CKD mice. This work sheds light on the mechanism of CKD-induced corneal calcification and provides tools to test calcification inhibitors for the prevention of this detrimental process.

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.

SMAD4-Dependent Signaling Pathway Involves in the Pathogenesis of TGFBR2-Related CE-like Phenotype

(1) Background: Our previous data indicated that disturbance of the Transforming Growth Factor beta (TGFB) signaling pathway via its Type-2 Receptor (TGFBR2) can cause a Corneal Ectasia (CE)-like phenotype. The purpose of this study is to elucidate whether the SMAD4-dependent signaling pathway is involved in the TGFBR2-related CE-like pathogenesis. (2) Methods: Smad4 was designed to be conditionally knocked out from keratocytes. Novel triple transgenic mice, KerartTA; Tet-O-Cre; Smad4flox/flox (Smad4kera-cko), were administered with doxycycline (Dox). Optical Coherence Tomography (OCT) was performed to examine Central Corneal Thickness (CCT), Corneal Radius, Anterior Chamber and CE-like phenotype and compared to the littermate Control group (Smad4Ctrl). (3) Results: The OCT revealed normal cornea in the Smad4Ctrl and a CE-like phenotype in the Smad4kera-cko cornea, in which the overall CCT in Smad4kera-cko was thinner than that of Smad4Ctrl at P42 (n = 6, p < 0.0001) and showed no significant difference when compared to that in Tgfbr2kera-cko. Furthermore, the measurements of the Anterior Chamber and Corneal Radius indicated a substantial ectatic cornea in the Smad4kera-cko compared to Smad4Ctrl. The H&E staining of Smad4kera-cko mimics the finding in the Tgfbr2kera-cko. The positive immunostaining of cornea-specific marker K12 indicating the cell fate of cornea epithelium remained unchanged in Smad4kera-cko and the Proliferating Cell Nuclear Antigen (PCNA) immunostaining further indicated an enhanced proliferation in the Smad4kera-cko. Both immunostainings recapitulated the finding in Tgfbr2kera-cko. The Masson's Trichrome staining revealed decreased collagen formation in the corneal stroma from both Smad4kera-cko and Tgfbr2kera-cko. The collagen type 1 (Col1a1) immunostaining further confirmed the reduction in collagen type 1 formation in Smad4kera-cko. (4) Conclusions: The aforementioned phenotypes in the Smad4kera-cko strain indicated that the SMAD4-dependent signaling pathway is involved in the pathogenesis of the CE-like phenotype observed in Tgfbr2kera-cko.

Genetic investigation in a four-generation Chinese family with congenital fibrosis of extraocular muscles and keratoconus

Here, we have reported the genetic and clinical characteristics of four generations of a family patient from China with congenital fibrosis of extraocular muscles 1 (CFEOM1) and keratoconus (KC). The history of diseases, clinical observations, and blood samples of all family members were collected. A total of 100 healthy participants were recruited as normal controls. The whole exome sequencing of the genomic DNA and polymerase chain reaction were performed on samples obtained from the controls and their family members to verify the gene variants. The functional analyses of the variants were performed by using different software. Two single nucleotide polymorphisms were detected in the proband and other patients in his families, including a heterozygous missense variation, g.39726207C > T (c.2860C > T, p.R954W, rs121912585), in the third highly conserved coiled-coil domain of KIF21A, and a heterozygous missense variant, g.30664732A > C (c.136A > C, p.S46R, rs200111443) in TGFBR2. The variant p.R954W in KIF21A was predicted to be pathogenic using software, whereas p.S46R in TGFBR2 was predicted to be of uncertain significance (VUS). Thus, KC might have occurred in the proband and his daughter because of a combination of genetic mutations and involuntary eye rubbing induced by CFEOM1. This is the first case of concomitant KC in a family having CFEOM1. Thus, the study provides new information about patients with KC having CFEOM1. Furthermore, the study suggests that attention should be paid to the early detection and diagnosis of KC in patients with CFEOM1.

Causal relationships between Gut microbiota and primary open-angle Glaucoma: A Mendelian randomization and mediation analysis of Glaucoma endophenotypes

Primary open-angle glaucoma (POAG) is a widespread condition responsible for irreversible blindness, and its prevalence is expected to increase substantially in the coming decades. Despite its significance, the exact cause of POAG remains elusive, necessitating a comprehensive exploration of its pathogenesis. Emerging research suggests a potential link between alterations in gut microbiota composition and POAG. However, establishing causality in these associations remains a challenge. In this study, we employed Mendelian randomization (MR) analysis to investigate the potential causal relationships between gut microbiota (GM) and POAG. Significant bacteria taxa were further analyzed with POAG endophenotypes. We utilized data from genome-wide association studies (GWAS) for GM and POAG, as well as for glaucoma endophenotypes, including intraocular pressure (IOP), retinal nerve fiber layer (RNFL) thickness, vertical cup-to-disc ratio (VCDR), and central corneal thickness (CCT). Univariable, multivariable MR and mediation effect analysis were conducted. Our analysis revealed that certain taxa, including phylum Euryarchaeota, genus Odoribacter, Rumnicoccaceae UCG009, Ruminiclostridium9, unknown genus id.2071, and Eubacterium rectale group, were associated with an increased risk of POAG. On the other hand, family Victivallaceae, Lacchnospiraceae, genus Lachnoclostridium, Oscillospira, Ruminococcaceae UCG011, Alloprevotella, and Faecalibacterium were found to be associated with a decreased risk of POAG. Furthermore, some of these taxa were found to be connected to glaucoma endophenotypes. Through further multivariable MR analysis, it was determined that IOP, VCDR, and CCT might played mediating roles between GM and POAG. In conclusion, this study utilizes MR analysis to elucidate potential causal associations between GM and POAG, providing insights into specific GM taxa that influence POAG risk and related endophenotypes. These findings emphasize the potential role of the gut microbiota in the pathogenesis of POAG and pave the way for future research and therapeutic interventions.

Integrating genetic regulation and single-cell expression with GWAS prioritizes causal genes and cell types for glaucoma

Primary open-angle glaucoma (POAG), characterized by retinal ganglion cell death, is a leading cause of irreversible blindness worldwide. However, its molecular and cellular causes are not well understood. Elevated intraocular pressure (IOP) is a major risk factor, but many patients have normal IOP. Colocalization and Mendelian randomization analysis of >240 POAG and IOP genome-wide association study (GWAS) loci and overlapping expression and splicing quantitative trait loci (e/sQTLs) in 49 GTEx tissues and retina prioritizes causal genes for 60% of loci. These genes are enriched in pathways implicated in extracellular matrix organization, cell adhesion, and vascular development. Analysis of single-nucleus RNA-seq of glaucoma-relevant eye tissues reveals that the POAG and IOP colocalizing genes and genome-wide associations are enriched in specific cell types in the aqueous outflow pathways, retina, optic nerve head, peripapillary sclera, and choroid. This study nominates IOP-dependent and independent regulatory mechanisms, genes, and cell types that may contribute to POAG pathogenesis.

Multi-ancestry genome-wide meta-analysis identifies novel basal cell carcinoma loci and shared genetic effects with squamous cell carcinoma

Basal cell carcinoma (BCC) is one of the most common malignancies worldwide, yet its genetic determinants are incompletely defined. We perform a European ancestry genome-wide association (GWA) meta-analysis and a Hispanic/Latino ancestry GWA meta-analysis and meta-analyze both in a multi-ancestry GWAS meta-analysis of BCC, totaling 50,531 BCC cases and 762,234 controls from four cohorts (GERA, Mass-General Brigham Biobank, UK Biobank, and 23andMe research cohort). Here we identify 122 BCC-associated loci, of which 36 were novel, and subsequently fine-mapped these associations. We also identify an association of the well-known pigment gene SLC45A2 as well as associations at RCC2 and CLPTM1L with BCC in Hispanic/Latinos. We examine these BCC loci for association with cutaneous squamous cell carcinoma (cSCC) in 16,407 SCC cases and 762,486 controls of European ancestry, and 33 SNPs show evidence of association. Our study findings provide important insights into the genetic basis of BCC and cSCC susceptibility.

Progressive Keratoconus in a Patient With Severe Pectus Excavatum and a Cartilage Oligomeric Matrix Protein Gene Mutation: A Case Report

INTRODUCTION

Keratoconus is a progressive ocular disorder associated with numerous systemic diseases, many of which affect the musculoskeletal system. Although the etiology and pathophysiology of the disorder remain elusive, recent studies suggest a significant role of genetic predisposition in the pathogenesis of keratoconus. This case report aims to elucidate a potential genetic association in a patient presenting with keratoconus, severe pectus excavatum, generalized muscular weakness, and skeletal deformities.

CASE DESCRIPTION

A 31-year-old Iranian man presented with progressively diminishing vision in both eyes over the years, eventually diagnosed with keratoconus. The patient's history and further examination indicated generalized muscular weakness, skeletal deformities, and severe pectus excavatum with cardiac and large vessel displacement. Whole-exome sequencing identified two heterozygous gene variants: one in the Cartilage Oligomeric Matrix Protein (COMP) gene and another in the Regulating Synaptic Membrane Exocytosis 1 gene. The patient's systemic and ocular symptoms, combined with the gene variants identified, suggested a connective tissue systemic disorder, potentially within the clinical spectrum of COMPopathies.

CONCLUSION

This is the first documented case of bilateral progressive keratoconus associated with severe pectus excavatum, generalized musculoskeletal dystrophy, and a COMP gene mutation. It highlights the necessity of continued search into the pathogenic genes of keratoconus, particularly in cases with coexisting systemic manifestations, to further our understanding of the etiology and pathogenesis of this complex disease.

Understanding the complex genetics and molecular mechanisms underlying glaucoma

Glaucoma is the leading cause of irreversible blindness worldwide. Currently the only effective treatment for glaucoma is to reduce the intraocular pressure, which can halt the progression of the disease. Highlighting the importance of identifying individuals at risk of developing glaucoma and those with early-stage glaucoma will help patients receive treatment before sight loss. However, some cases of glaucoma do not have raised intraocular pressure. In fact, glaucoma is caused by a variety of different mechanisms and has a wide range of different subtypes. Understanding other risk factors, the underlying mechanisms, and the pathology of glaucoma might lead to novel treatments and treatment of underlying diseases. In this review we present the latest research into glaucoma including the genetics and molecular basis of the disease.

Colocalization of corneal resistance factor GWAS loci with GTEx e/sQTLs highlights plausible candidate causal genes for keratoconus postnatal corneal stroma weakening

Background: Genome-wide association studies (GWAS) for corneal resistance factor (CRF) have identified 100s of loci and proved useful to uncover genetic determinants for keratoconus, a corneal ectasia of early-adulthood onset and common indication of corneal transplantation. In the current absence of studies to probe the impact of candidate causal variants in the cornea, we aimed to fill some of this knowledge gap by leveraging tissue-shared genetic effects. Methods: 181 CRF signals were examined for evidence of colocalization with genetic signals affecting steady-state gene transcription and splicing in adult, non-eye, tissues of the Genotype-Tissue Expression (GTEx) project. Expression of candidate causal genes thus nominated was evaluated in single cell transcriptomes from adult cornea, limbus and conjunctiva. Fine-mapping and colocalization of CRF and keratoconus GWAS signals was also deployed to support their sharing causal variants. Results and discussion: 26.5% of CRF causal signals colocalized with GTEx v8 signals and nominated genes enriched in genes with high and specific expression in corneal stromal cells amongst tissues examined. Enrichment analyses carried out with nearest genes to all 181 CRF GWAS signals indicated that stromal cells of the limbus could be susceptible to signals that did not colocalize with GTEx's. These cells might not be well represented in GTEx and/or the genetic associations might have context specific effects. The causal signals shared with GTEx provide new insights into mediation of CRF genetic effects, including modulation of splicing events. Functionally relevant roles for several implicated genes' products in providing tensile strength, mechano-sensing and signaling make the corresponding genes and regulatory variants prime candidates to be validated and their roles and effects across tissues elucidated. Colocalization of CRF and keratoconus GWAS signals strengthened support for shared causal variants but also highlighted many ways into which likely true shared signals could be missed when using readily available GWAS summary statistics.

Transforming growth factor beta receptor 2 (Tgfbr2) deficiency in keratocytes results in corneal ectasia

PURPOSE

We hypothesized that Transforming growth factor beta receptor 2 (Tgfbr2) deletion in keratocyte (Tgfbr2kera-cko), the corneal stroma cell, can result in corneal thinning and generate a potential model for Cornea Ectasia (CE).

METHODS

Corneal thickness of Tgfbr2kera-cko and Tgfbr2Ctrl was examined with Optical Coherence Tomography (OCT) at post-natal (P) days 42 and 70, respectively. Histological H&E staining, transmission electron micrograph (TEM), and immunofluorescence staining (IFS) were harnessed to examine corneal cell morphology, proliferation, differentiation, and collagen fibrils.

RESULTS

Slit-Lamp revealed that corneas were transparent in both Tgfbr2kera-cko and Tgfbr2Ctrl, however, Tgfbr2kera-cko cornea was 33.5% and 42.9% thinner as compared with those of Tgfbr2Ctrl at P42 and P70, respectively. H&E and semithin section staining with toluidine blue-O confirmed that Tgfbr2kera-cko cornea has a thinner stroma. In contrast, the epithelium in Tgfbr2kera-cko was substantially thicker. The cell proliferation marker Ki67 expression level increased ∼9% in Tgfbr2kera-cko corneal epithelium as compared with that in Tgfbr2Ctrl, however, the Krt14 and Krt12 expression pattern was not obviously changed in Tgfbr2kera-cko corneal epithelium. It was noticed that Col1a1 expression was substantially reduced in Tgfbr2kera-cko as compared with that in Tgfbr2Ctrl. TEM showed that keratocytes were unhealthy and stromal collagen fibril density was significantly reduced in Tgfbr2kera-cko as compared with that in Tgfbr2Ctrl cornea. Moreover, mechanical eye-rubbing on Tgfbr2kera-cko resulted in corneal hydrops and edema.

CONCLUSION

Tgfbr2 in keratocytes is indispensable for the corneal stroma at postnatal homeostasis. The cornea phenotype manifested in these Tgfbr2kera-cko mice resembles corneal ectasia disease in humans.

A multiethnic genome-wide analysis of 19,420 individuals identifies novel loci associated with axial length and shared genetic influences with refractive error and myopia

Introduction: Long axial length (AL) is a risk factor for myopia. Although family studies indicate that AL has an important genetic component with heritability estimates up to 0.94, there have been few reports of AL-associated loci. Methods: Here, we conducted a multiethnic genome-wide association study (GWAS) of AL in 19,420 adults of European, Latino, Asian, and African ancestry from the Genetic Epidemiology Research on Adult Health and Aging (GERA) cohort, with replication in a subset of the Consortium for Refractive Error and Myopia (CREAM) cohorts of European or Asian ancestry. We further examined the effect of the identified loci on the mean spherical equivalent (MSE) within the GERA cohort. We also performed genome-wide genetic correlation analyses to quantify the genetic overlap between AL and MSE or myopia risk in the GERA European ancestry sample. Results: Our multiethnic GWA analysis of AL identified a total of 16 genomic loci, of which 5 are novel. We found that all AL-associated loci were significantly associated with MSE after Bonferroni correction. We also found that AL was genetically correlated with MSE (rg = -0.83; SE, 0.04; p = 1.95 × 10-89) and myopia (rg = 0.80; SE, 0.05; p = 2.84 × 10-55). Finally, we estimated the array heritability for AL in the GERA European ancestry sample using LD score regression, and found an overall heritability estimate of 0.37 (s.e. = 0.04). Discussion: In this large and multiethnic study, we identified novel loci, associated with AL at a genome-wide significance level, increasing substantially our understanding of the etiology of AL variation. Our results also demonstrate an association between AL-associated loci and MSE and a shared genetic basis between AL and myopia risk.

Cornea plana in a family from Pakistan: Case series and literature review on the principles of management

Cornea plana (CP) is a rare ocular condition existing in two distinct clinical and hereditary forms: a milder, autosomal dominant type I and a more severe, autosomal recessive type II. The condition is more commonly found in Finnish, Saudi, and Czech families. We report three brothers from a consanguineous marriage that presented with complaints of decreased vision of varying degrees. All three of them have blue, thick, and hazy corneas with shallow anterior chamber depths. The additional features of CP type II were seen in the older two brothers including arcus lipoids, ill-demarcated limbus, and an accommodative squint. They were managed by the correction of refractive errors through spectacles and detailed counseling with follow-up visits to look for progressive complications. The management is mainly centered around optically or surgically correcting the developmental anomalies. This is complimented with proper genetic counseling and regular follow-up visits to look for and manage complications. There are, however, novel therapies that can be considered in these patients including corneal transplants or corneal stromal stem cellular therapies.

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 (PSMR = 2.66 × 10-6), KANSL1-AS1 (PSMR = 2.78 × 10-6), and RP11-259G18.2 (PSMR = 2.91 × 10-6) were the top three genes using the GTEx eQTL data. LRRC37A4 (PSMR = 1.19 × 10-5), MGC57346 (PSMR = 1.19 × 10-5), and RNF167 (PSMR = 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.

Pierre CL, Bimschleger H, Gao J, Cheng R, Mohammadi P, Chen H, Palmer AA, Polesskaya O, Jablonski MM. Genome-wide association study finds multiple loci associated with intraocular pressure in HS rats

Elevated intraocular pressure (IOP) is influenced by environmental and genetic factors. Increased IOP is a major risk factor for most types of glaucoma, including primary open angle glaucoma (POAG). Investigating the genetic basis of IOP may lead to a better understanding of the molecular mechanisms of POAG. The goal of this study was to identify genetic loci involved in regulating IOP using outbred heterogeneous stock (HS) rats. HS rats are a multigenerational outbred population derived from eight inbred strains that have been fully sequenced. This population is ideal for a genome-wide association study (GWAS) owing to the accumulated recombinations among well-defined haplotypes, the relatively high allele frequencies, the accessibility to a large collection of tissue samples, and the large allelic effect size compared to human studies. Both male and female HS rats (N = 1,812) were used in the study. Genotyping-by-sequencing was used to obtain ∼3.5 million single nucleotide polymorphisms (SNP) from each individual. SNP heritability for IOP in HS rats was 0.32, which agrees with other studies. We performed a GWAS for the IOP phenotype using a linear mixed model and used permutation to determine a genome-wide significance threshold. We identified three genome-wide significant loci for IOP on chromosomes 1, 5, and 16. Next, we sequenced the mRNA of 51 whole eye samples to find cis-eQTLs to aid in identification of candidate genes. We report 5 candidate genes within those loci: Tyr, Ctsc, Plekhf2, Ndufaf6 and Angpt2. Tyr, Ndufaf6 and Angpt2 genes have been previously implicated by human GWAS of IOP-related conditions. Ctsc and Plekhf2 genes represent novel findings that may provide new insight into the molecular basis of IOP. This study highlights the efficacy of HS rats for investigating the genetics of elevated IOP and identifying potential candidate genes for future functional testing.

Genetic variants associated with glaucomatous visual field loss in primary open-angle glaucoma

Primary open-angle glaucoma (POAG) is characterized by a progressive optic neuropathy with visual field loss. To investigate the genetic variants associated with visual field loss in POAG, Japanese POAG patients (n = 426) and control subjects (n = 246) were genotyped for 22 genetic variants predisposing to POAG that can be classified into those associated with intraocular pressure (IOP) elevation (IOP-related genetic variants) and optic nerve vulnerability independent of IOP (optic nerve-related genetic variants). The genetic risk score (GRS) of the 17 IOP-related and five optic nerve-related genetic variants was calculated, and the associations between the GRS and the mean deviation (MD) of automated static perimetry as an indicator of the severity of visual field loss and pattern standard deviation (PSD) as an indicator of the focal disturbance were evaluated. There was a significant association (Beta = - 0.51, P = 0.0012) between the IOP-related GRS and MD. The severity of visual field loss may depend on the magnitude of IOP elevation induced by additive effects of IOP-related genetic variants. A significant association (n = 135, Beta = 0.65, P = 0.0097) was found between the optic nerve-related, but not IOP-related, GRS and PSD. The optic nerve-related (optic nerve vulnerability) and IOP-related (IOP elevation) genetic variants may play an important role in the focal and diffuse visual field loss respectively. To our knowledge, this is the first report to show an association between additive effects of genetic variants predisposing to POAG and glaucomatous visual field loss, including severity and focal/diffuse disturbance of visual field loss, in POAG.

Optic neuropathy associated with TGFβ dysregulation in mice with a glaucoma-causative mutation of ADAMTS10

Glaucoma is a neurodegenerative disease that causes irreversible blindness due to loss of retinal ganglion cells (RGCs) and their axons. We previously identified a G661R mutation of ADAMTS10 (A Disintegrin And Metalloproteinase with ThromboSpondin type 1 motif 10) as the disease-causing mutation in a beagle model of glaucoma. ADAMTS10 is a secreted matrix metalloproteinase that belongs to the ADAMTS family which is involved in extracellular matrix (ECM) turnover. Previous studies have shown that ADAMTS10 binds fibrillin microfibrils, promotes their formation, and influences their fibrillin isoform composition. Here, we established a mouse model carrying the G661R mutation of ADAMTS10 (ADAMTS10G661R/G661R) to investigate its ocular phenotypes related to glaucoma and to explore possible functions of ADAMTS10. We found that ADAMTS10 was expressed in the inner retina and along RGC axons in the optic nerve. However, ADAMTS10 was not colocalized with fibrillin microfibrils in these tissues, suggesting fibrillin-independent function for ADAMTS10. In electroretinogram experiments, we found that ADAMTS10G661R/G661R mice had reduced amplitude of retinal responses to dim light stimulus, indicating RGC dysfunction. The reduced RGC function coincided with RGC axon structural changes manifested as smaller optic nerves and fewer optic nerve axons, which may contribute to glaucoma. The reduced number of optic nerve axons found for ADAMTS10G661R/G661R mice occurred early, suggesting developmental deficits. Subsequent experiments found increased apoptosis in the retina of ADAMTS10G661R/G661R mice during postnatal development, which could result in fewer RGCs produced, accounting for fewer optic nerve axons in adulthood. Consistent with a protective effect of transforming growth factor β (TGFβ) signaling against apoptosis during retinal development as shown previously by others, we found increased apoptosis accompanied by decreased TGFβ signaling in the developing retina of ADAMTS10G661R/G661R mice, suggesting a novel role for ADAMTS10 in regulating TGFβ signaling which could involve direct interaction between ADAMTS10 and latent TGFβ.

Heritability of Corneal Parameters in Nuclear Families With Keratoconus

Purpose

This study aimed to investigate the heritability of corneal parameters obtained by Pentacam in nuclear families with keratoconus (KC).

Methods

A total of 82 patients with KC and their biological parents (n = 164) were recruited in the current study. All subjects underwent corneal tomography with Pentacam. Family units were analyzed to calculate the heritability of corneal parameters by linear mixed effects model using the R statistical software.

Results

The pachymetry at apex, pupil, and thinnest point were all significantly heritable at 43.26%, 42.63%, and 43.09%, respectively. The heritability of flat meridian keratometry, steep meridian keratometry, and mean keratometry in the anterior surface were 10.36%, 9.05%, and 10.21%, respectively, and that of flat meridian keratometry, steep meridian keratometry, and mean keratometry in the posterior surface were 8.44%, 9.67%, and 9.06%, respectively. The posterior radius of curvature had higher heritability in comparison with anterior radius of curvature (19.16% vs. 14.37%). Moreover, among combined topometric indices, the heritability of index of vertical asymmetry was the highest (19.49%), and that of central keratoconus index was the lowest (6.64%).

Conclusions

The present study demonstrated a substantial heritability of corneal parameters in nuclear families with KC. The pachymetric indices are heritable and may be suitable as KC endophenotypes, suggesting a necessity to discover the genes associated with corneal thickness in KC.

Translational Relevance

The pachymetric indices are heritable and may be suitable as KC endophenotypes, indicating that the pachymetric indices might be a corneal characteristic to predict the occurrence of KC.

The genetic basis for adult onset glaucoma: Recent advances and future directions

Glaucoma, a diverse group of eye disorders that results in the degeneration of retinal ganglion cells, is the world's leading cause of irreversible blindness. Apart from age and ancestry, the major risk factor for glaucoma is increased intraocular pressure (IOP). In primary open-angle glaucoma (POAG), the anterior chamber angle is open but there is resistance to aqueous outflow. In primary angle-closure glaucoma (PACG), crowding of the anterior chamber angle due to anatomical alterations impede aqueous drainage through the angle. In exfoliation syndrome and exfoliation glaucoma, deposition of white flaky material throughout the anterior chamber directly interfere with aqueous outflow. Observational studies have established that there is a strong hereditable component for glaucoma onset and progression. Indeed, a succession of genome wide association studies (GWAS) that were centered upon single nucleotide polymorphisms (SNP) have yielded more than a hundred genetic markers associated with glaucoma risk. However, a shortcoming of GWAS studies is the difficulty in identifying the actual effector genes responsible for disease pathogenesis. Building on the foundation laid by GWAS studies, research groups have recently begun to perform whole exome-sequencing to evaluate the contribution of protein-changing, coding sequence genetic variants to glaucoma risk. The adoption of this technology in both large population-based studies as well as family studies are revealing the presence of novel, protein-changing genetic variants that could enrich our understanding of the pathogenesis of glaucoma. This review will cover recent advances in the genetics of primary open-angle glaucoma, primary angle-closure glaucoma and exfoliation glaucoma, which collectively make up the vast majority of all glaucoma cases in the world today. We will discuss how recent advances in research methodology have uncovered new risk genes, and how follow up biological investigations could be undertaken in order to define how the risk encoded by a genetic sequence variant comes into play in patients. We will also hypothesise how data arising from characterising these genetic variants could be utilized to predict glaucoma risk and the manner in which new therapeutic strategies might be informed.

Replication of gene polymorphisms associated with periodontitis-related traits in an elderly cohort: the Washington Heights/Inwood Community Aging Project Ancillary Study of Oral Health

AIM

We sought to replicate findings from published genome-wide association studies (GWAS), linking specific candidate gene loci with periodontitis-related clinical/microbial traits.

MATERIALS AND METHODS

In the published GWAS, a total of 2196 single nucleotide polymorphisms associated with periodontitis-related traits at a p ≤ 5 × 10^-6 and mapped to 136 gene loci. The replication cohort included 1124 individuals, 65-98 years old (67% female, 45% Hispanic, 30% Black, 23% White) with available genome-wide genotypes and full-mouth periodontal status. Microbial profiles using checkerboard DNA-DNA hybridization and 16SrRNA sequencing were available from 912 and 739 participants, respectively.

RESULTS

Using gene-specific p-values after linkage disequilibrium pruning, the following gene/phenotype associations replicated successfully: CLEC19A with edentulism and %teeth with pocket depth (PD) ≥4 mm; IL37, HPVC1, TRPS1, ABHD12B, LDLRAD4 (C180rF1), TGM3, and GRK5 with %teeth with PD ≥4 mm; DAB2IP with presence of PD ≥6 mm; KIAA1715(LNPK), ROBO2, RAB28, LINC01017, NELL1, LDLRAD4(C18orF1), and CRYBB2P1 with %teeth with clinical attachment level (CAL) ≥3 mm; RUNX2 and LAMA2 with %teeth with CAL ≥5 mm; and KIAA1715(LNPK) with high colonization by Aggregatibacter actinomycetemcomitans. In addition, CLEC19A, IQSEC1, and EMR1 associated with microbial abundance based on checkerboard data, LBP and NCR2 with abundance based on sequencing data, and NCR2 with microbial diversity based on sequencing data.

CONCLUSIONS

Several gene loci identified in published GWAS as associated with periodontitis-related phenotypes replicated successfully in an elderly cohort.

Association of Novel Loci With Keratoconus Susceptibility in a Multitrait Genome-Wide Association Study of the UK Biobank Database and Canadian Longitudinal Study on Aging

Importance

Keratoconus can be a debilitating corneal ectasia in which the cornea thins, bulges, and steepens into a conical shape. Early features of keratoconus include myopia and irregular astigmatism, which affect vision and can be treated with contact lenses, collagen cross-linking, or, in advanced cases, corneal transplant. Recent estimates of the prevalence of keratoconus based on results of Scheimpflug imaging in young adults are as high as 1.2%. However, obtaining very large keratoconus data sets for a genome-wide association study (GWAS) is problematic because few population studies include Scheimpflug imaging and because severe keratoconus is relatively rare.

Objective

To identify novel keratoconus loci using corneal resistance factor (CRF) and central corneal thickness (CCT).

Design, Setting, and Participants

This multitrait GWAS used European ancestry CRF data from UK Biobank (UKB) (n = 105 427) and the Canadian Longitudinal Study on Aging (CLSA) (n = 18 307) and European ancestry CCT data from the International Glaucoma Genetics Consortium (IGGC) (n = 17 803). The CRF and CCT variants in published keratoconus data sets (4669 cases and 116 547 controls) were compared. The data set from UKB was compiled March 24, 2020; data were released from the CLSA in July 2020; and IGGC data were available from May 1, 2018.

Main Outcomes and Measures

Association of CRF and CCT variants with keratoconus risk.

Results

The GWAS included 4 cohorts: 105 427 UKB European ancestry (56 134 women [53.2%] and 49 293 men [46.7%]; mean [SD] age, 57 [8] years), 5029 UKB South Asian ancestry (2368 women [47.1%] and 2661 men [52.9%]; mean [SD] age, 54 [8] years), 902 UKB East Asian ancestry (622 women [68.9%] and 280 men [31.0%]; mean [SD] age, 53 [8] years), and 18 307 CLSA European ancestry (9260 women [50.6%] and 9047 men [49.4%]; mean [SD] age, 63 [10] years) participants. A total of 369 CRF and 233 CCT loci were identified, including 36 novel CRF loci and 114 novel CCT loci. Twenty-nine CRF loci and 24 CCT loci were associated with keratoconus. Polygenic risk scores (PRS) were constructed using CRF- and CCT-associated variants and published keratoconus variants. The PRS result showed that adding a CRF- or CCT-based PRS to the keratoconus PRS from previously published variants improved the prediction area under the receiver operating characteristic curve (from 0.705 to 0.756 for CRF and from 0.715 to 0.755 for CCT).

Conclusions and Relevance

These findings support the use of multitrait modeling of corneal parameters in a relatively large data set to identify new keratoconus risk loci and enhance polygenic risk score models.

Human and rat skeletal muscle single-nuclei multi-omic integrative analyses nominate causal cell types, regulatory elements, and SNPs for complex traits

Skeletal muscle accounts for the largest proportion of human body mass, on average, and is a key tissue in complex diseases and mobility. It is composed of several different cell and muscle fiber types. Here, we optimize single-nucleus ATAC-seq (snATAC-seq) to map skeletal muscle cell-specific chromatin accessibility landscapes in frozen human and rat samples, and single-nucleus RNA-seq (snRNA-seq) to map cell-specific transcriptomes in human. We additionally perform multi-omics profiling (gene expression and chromatin accessibility) on human and rat muscle samples. We capture type I and type II muscle fiber signatures, which are generally missed by existing single-cell RNA-seq methods. We perform cross-modality and cross-species integrative analyses on 33,862 nuclei and identify seven cell types ranging in abundance from 59.6% to 1.0% of all nuclei. We introduce a regression-based approach to infer cell types by comparing transcription start site-distal ATAC-seq peaks to reference enhancer maps and show consistency with RNA-based marker gene cell type assignments. We find heterogeneity in enrichment of genetic variants linked to complex phenotypes from the UK Biobank and diabetes genome-wide association studies in cell-specific ATAC-seq peaks, with the most striking enrichment patterns in muscle mesenchymal stem cells (∼3.5% of nuclei). Finally, we overlay these chromatin accessibility maps on GWAS data to nominate causal cell types, SNPs, transcription factor motifs, and target genes for type 2 diabetes signals. These chromatin accessibility profiles for human and rat skeletal muscle cell types are a useful resource for nominating causal GWAS SNPs and cell types.

Rare single nucleotide variants in COL5A1 promoter do not play a major role in keratoconus susceptibility associated with rs1536482

BACKGROUND

Keratoconus is a chronic degenerative disorder of the cornea characterized by thinning and cone-shaped protrusions. Although genetic factors play a key role in keratoconus development, the etiology is still under investigation. The occurrence of single-nucleotide polymorphisms (SNPs) associated with keratoconus in Russian patients is poorly studied. The purpose of this study was to validate whether three reported keratoconus-associated SNPs (rs1536482 near the COL5A1 gene, rs2721051 near the FOXO1 gene, rs1324183 near the MPDZ gene) are also actual for a Russian cohort of patients. Additionally, we investigated the COL5A1 promoter sequence for single-nucleotide variants (SNVs) in a subgroup of keratoconus patients with at least one rs1536482 minor allele (rs1536482+) to assess the role of these SNVs in keratoconus susceptibility associated with rs1536482.

METHODS

This case-control study included 150 keratoconus patients and two control groups (main and additional, 205 and 474 participants, respectively). We performed PCR targeting regions flanking SNVs and the COL5A1 promoter, followed by Sanger sequencing of amplicons. The additional control group was genotyped using an SNP array.

RESULTS

The minor allele frequency was significantly different between the keratoconus and control cohorts (main and combined) for rs1536482, rs2721051, and rs1324183 (p-value < 0.05). The rare variants rs1043208782 and rs569248712 were found in the COL5A1 promoter in two out of 94 rs1536482+ keratoconus patients.

CONCLUSION

rs1536482, rs2721051, and rs1324183 were associated with keratoconus in a Russian cohort. SNVs in the COL5A1 promoter do not play a major role in keratoconus susceptibility associated with rs1536482.

A Mouse Model of Brittle Cornea Syndrome caused by mutation in Zfp469

Brittle cornea syndrome (BCS) is a rare recessive condition characterised by extreme thinning of the cornea and sclera. BCS results from loss-of-function mutations in the poorly understood genes ZNF469 or PRDM5. In order to determine the function of ZNF469 and to elucidate pathogenic mechanisms, we used genome editing to recapitulate a human ZNF469 BCS mutation in the orthologous mouse gene Zfp469. Ophthalmic phenotyping showed that homozygous Zfp469 mutation causes significant central and peripheral corneal thinning arising from reduced stromal thickness. Expression of key components of the corneal stroma in primary keratocytes from Zfp469^BCS/BCS mice is affected, including decreased Col1a1 and Col1a2 expression. This alters the collagen type I/collagen type V ratio and results in collagen fibrils with smaller diameter and increased fibril density in homozygous mutant corneas, correlating with decreased biomechanical strength in the cornea. Cell-derived matrices generated by primary keratocytes show reduced deposition of collagen type I, offering an in vitro model for stromal dysfunction. Work remains to determine whether modulating ZNF469 activity will have therapeutic benefit in BCS or in conditions such as keratoconus in which the cornea thins progressively.

This article has an associated First Person interview with the first author of the paper.

Summary: A mouse model of brittle cornea syndrome was created to elucidate molecular mechanisms underlying the pathology of this rare connective tissue disorder in which extremely thin corneas rupture, causing irreversible blindness.

Mendelian randomization analysis identified genes pleiotropically associated with central corneal thickness

OBJECTIVE

To prioritize genes that were pleiotropically or potentially causally associated with central corneal thickness (CCT).

METHODS

We applied the summary data-based Mendelian randomization (SMR) method integrating summarized data of genome-wide association study (GWAS) on CCT and expression quantitative trait loci (eQTL) data to identify genes that were pleiotropically associated with CCT. We performed separate SMR analysis using CAGE eQTL data and GTEx eQTL data. SMR analyses were done for participants of European and East Asian ancestries, separately.

RESULTS

We identified multiple genes showing pleiotropic association with CCT in the participants of European ancestry. CLIC3 (ILMN_1796423; PSMR = 4.15 × 10- 12), PTGDS (ILMN_1664464; PSMR = 6.88 × 10- 9) and C9orf142 (ILMN_1761138; PSMR = 8.09 × 10- 9) were the top three genes using the CAGE eQTL data, and RP11-458F8.4 (ENSG00000273142.1; PSMR = 5.89 × 10- 9), LCNL1 (ENSG00000214402.6; PSMR = 5.67 × 10- 8), and PTGDS (ENSG00000107317.7; PSMR = 1.92 × 10- 7) were the top three genes using the GTEx eQTL data. No genes showed significantly pleiotropic association with CCT in the participants of East Asian ancestry after correction for multiple testing.

CONCLUSION

We identified several genes pleiotropically associated with CCT, some of which represented novel genes influencing CCT. Our findings provided important leads to a better understanding of the genetic factors influencing CCT, and revealed potential therapeutic targets for the treatment of primary open-angle glaucoma and keratoconus.

Clinical, Demographic, and Tomographic Aspects Related to Iris Mammillations Among Patients With Keratoconus: A Cross-Sectional Study

PURPOSE

Iris mammillations (IM) were previously described in patients with keratoconus, but the clinical implications of this finding have never been studied. Our aim was to investigate demographic, tomographic, and clinical characteristics potentially associated with the presence of IM among patients with keratoconus.

METHODS

This was a cross-sectional study performed among patients with keratoconus in a public-affiliated university hospital. All patients under follow-up were considered eligible to participate in the study. Participants were evaluated by 2 trained ophthalmologists and submitted to corneal tomography (Pentacam). Selected demographic, clinical, and tomographic characteristics were assessed and compared among participants with IM (IM group) and without IM (No-IM group) using the Wilcoxon test or 2-tailed Fisher exact test, as appropriate.

RESULTS

The study population consisted of 106 subjects and 19 (17.9%) presented with IM. The median age and interquartile range were 18 years old (14-24) in the IM group and 20 years old (17-24) in the No-IM group (P = 0.135). The female proportion was 47.3% in the IM group and 52.8% in the No-IM group (P = 0.801). Median (interquartile ranges) pachymetric values of the right eyes were 498 (466-525) for the IM group and 459 (421-482) for the No-IM group (P = 0.005). For the left eyes, the values were 490 (456-523) in the IM group and 450 (418-485) in the No-IM group (P = 0.024).

CONCLUSIONS

Subjects with keratoconus presenting with IM have thicker corneas than those without IM. Follow-up studies should be performed to evaluate the clinical implications of this finding.

A large multiethnic GWAS meta-analysis of cataract identifies new risk loci and sex-specific effects

Cataract is the leading cause of blindness among the elderly worldwide and cataract surgery is one of the most common operations performed in the United States. As the genetic etiology of cataract formation remains unclear, we conducted a multiethnic genome-wide association meta-analysis, combining results from the GERA and UK Biobank cohorts, and tested for replication in the 23andMe research cohort. We report 54 genome-wide significant loci, 37 of which were novel. Sex-stratified analyses identified CASP7 as an additional novel locus specific to women. We show that genes within or near 80% of the cataract-associated loci are significantly expressed and/or enriched-expressed in the mouse lens across various spatiotemporal stages as per iSyTE analysis. Furthermore, iSyTE shows 32 candidate genes in the associated loci have altered gene expression in 9 different gene perturbation mouse models of lens defects/cataract, suggesting their relevance to lens biology. Our work provides further insight into the complex genetic architecture of cataract susceptibility.

Genetic Determinants of Intraocular Pressure

Intraocular pressure (IOP) is the cardinal and only modifiable risk factor for glaucoma, the leading cause of irreparable blindness worldwide. Twin and family studies estimate the heritability of IOP to be 40-70%, and linkage studies for IOP have identified numerous loci. Mutations in MYOC can cause markedly elevated IOP and aggressive glaucoma often requiring surgical intervention. However, the majority of the genetic basis for raised IOP and glaucoma in populations is complex, and recent large genome-wide association studies (GWASs) have identified over 100 common variants that contribute to IOP variation. In combination, these loci are predictive for primary open-angle glaucoma in independent populations, achieving an area under the receiver operating characteristic curve of 76% for high-pressure primary open-angle glaucoma; this suggests the possibility of targeted screening in the future. Additionally, GWAS findings have identified important biological pathways underlying IOP regulation, including lymphangiogenesis and lipid metabolism, providing novel targets for new therapies. Expected final online publication date for the Annual Review of Vision Science, Volume 7 is September 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

IL1A and IL1B gene polymorphisms and keratoconus susceptibility: evidence from an updated meta-analysis

Background: Several single-nucleotide polymorphisms (SNPs) in IL1B genes have been associated with KTCN. However, the results of these studies were not conclusive. This meta-analysis association study is aimed to quantitatively estimate the association of IL1B rs16944 (g.4490T>C) and rs1143627 (g.4970C>T), and IL1A rs2071376 (c.615 + 169C>A) polymorphisms with KTCN susceptibility.Materials and Methods: Systematic literature search was performed in Web of Science, MEDLINE, PubMed, Scopus, and Google Scholar databases. The odds ratios (ORs) and 95% confidence intervals (CI) were calculated assuming different contrasted genetic models.Results: The reference T allele of IL1B (g.4490T>C) polymorphism was significantly associated with decreased KTCN risk under all assessed genetic models. Regarding the reference C allele of IL1B (g.4970C>T) polymorphism, decreased risk of KTCN was found. The reference C allele of IL1A (c.615 + 169C>A) polymorphism conferred a decreased risk of KTCN under heterozygous codominant (AC vs. AA), homozygous codominant (CC vs. AA), and dominant (AC+CC vs. AA) genetic models. The pooling estimates showed that the T C haplotype was associated with a significant increase in KTCN risk. In contrast, the T T haplotype was correlated with a decreased risk of KTCN. With the assumption of a prior probability of 0.25, the false-positive report probability (FPRP) values were less than 0.2, indicating the observed significant associations were notable.Conclusion: These findings propose that the studied IL1B polymorphisms and the IL1A variation have opposite effects on KTCN susceptibility. More large-scale replication studies are warranted to illuminate the precise role of these SNPs on the etiology of eye disorders.

Evaluating the association between single nucleotide polymorphisms in the stonin 2 (STON2) gene and keratoconus in a Han Chinese population

Background

A recent genome-wide association study (GWAS) identified a significant association between the single nucleotide polymorphism (SNP) rs2371597 in the stonin 2 gene (STON2) and keratoconus (KCTN) susceptibility. The current study further explored the association between STON2 and KCTN susceptibility in an independent Han Chinese population.

Methods

Three SNPs (rs2371597, rs8004137, and rs8008602) located in the STON2 gene were examined in 164 Han Chinese patients with KCTN and 239 age- and gender-matched healthy subjects. The TaqMan SNP genotyping assays were performed, and the LDlink, RegulomeDB, and PLINK package were applied for data analyses. The gene expression levels of STON2 were investigated in various murine organ tissues using quantitative real-time polymerase chain reaction (qRT-PCR).

Results

The SNP rs2371597 was significantly associated with KCTN risk in this Han Chinese population. The frequency of the C allele in KCTN patients was significantly higher than that in healthy subjects [34.8% vs. 26.6%; odds ratio (OR) =1.47; 95% confidence interval (CI): 1.08 to 2.02; P=0.01409]. The genotype distribution of the SNP rs2371597 was also significantly different between KCTN patients and controls. The other two genotyped SNPs allele and genotypic frequencies were not remarkably different between the KCTN group and the control group. However, the haplotype CAT formed by the three SNPs was substantially associated with the risk of KCTN (P=0.04101). Also, gene expression pattern analysis showed a relatively higher expression of STON2 in the cornea in comparison to other tissues.

Conclusions

The current study demonstrated that SNPs in the STON2 gene were associated with an increased risk of developing KCTN in this Han Chinese population, suggesting that the STON2 gene may play an important role in the etiology of KCTN.

Resequencing of candidate genes for Keratoconus reveals a role for Ehlers-Danlos Syndrome genes

The involvement of genetic factors in the pathogenesis of KC has long been recognized but the identification of variants affecting the underlying protein functions has been challenging. In this study, we selected 34 candidate genes for KC based on previous whole-exome sequencing (WES) and the literature, and resequenced them in 745 KC patients and 810 ethnically matched controls from Belgium, France and Italy. Data analysis was performed using the single variant association test as well as gene-based mutation burden and variance components tests. In our study, we detected enrichment of genetic variation across multiple gene-based tests for the genes COL2A1, COL5A1, TNXB, and ZNF469. The top hit in the single variant association test was obtained for a common variant in the COL12A1 gene. These associations were consistently found across independent subpopulations. Interestingly, COL5A1, TNXB, ZNF469 and COL12A1 are all known Ehlers-Danlos Syndrome (EDS) genes. Though the co-occurrence of KC and EDS has been reported previously, this study is the first to demonstrate a consistent role of genetic variants in EDS genes in the etiology of KC. In conclusion, our data show a shared genetic etiology between KC and EDS, and clearly confirm the currently disputed role of ZNF469 in disease susceptibility for KC.

A multi-ethnic genome-wide association study implicates collagen matrix integrity and cell differentiation pathways in keratoconus

Keratoconus is characterised by reduced rigidity of the cornea with distortion and focal thinning that causes blurred vision, however, the pathogenetic mechanisms are unknown. It can lead to severe visual morbidity in children and young adults and is a common indication for corneal transplantation worldwide. Here we report the first large scale genome-wide association study of keratoconus including 4,669 cases and 116,547 controls. We have identified significant association with 36 genomic loci that, for the first time, implicate both dysregulation of corneal collagen matrix integrity and cell differentiation pathways as primary disease-causing mechanisms. The results also suggest pleiotropy, with some disease mechanisms shared with other corneal diseases, such as Fuchs endothelial corneal dystrophy. The common variants associated with keratoconus explain 12.5% of the genetic variance, which shows potential for the future development of a diagnostic test to detect susceptibility to disease.

Smarce1 and Tensin 4 Are Putative Modulators of Corneoscleral Stiffness

The biomechanical properties of the cornea and sclera are important in the onset and progression of multiple ocular pathologies and vary substantially between individuals, yet the source of this variation remains unknown. Here we identify genes putatively regulating corneoscleral biomechanical tissue properties by conducting high-fidelity ocular compliance measurements across the BXD recombinant inbred mouse set and performing quantitative trait analysis. We find seven cis-eQTLs and non-synonymous SNPs associating with ocular compliance, and show by RT-qPCR and immunolabeling that only two of the candidate genes, Smarce1 and Tns4, showed significant expression in corneal and scleral tissues. Both have mechanistic potential to influence the development and/or regulation of tissue material properties. This work motivates further study of Smarce1 and Tns4 for their role(s) in ocular pathology involving the corneoscleral envelope as well as the development of novel mouse models of ocular pathophysiology, such as myopia and glaucoma.

Molecular Genetics of Glaucoma: Subtype and Ethnicity Considerations

Glaucoma, the world’s leading cause of irreversible blindness, is a complex disease, with differential presentation as well as ethnic and geographic disparities. The multifactorial nature of glaucoma complicates the study of genetics and genetic involvement in the disease process. This review synthesizes the current literature on glaucoma and genetics, as stratified by glaucoma subtype and ethnicity. Primary open-angle glaucoma (POAG) is the most common cause of glaucoma worldwide, with the only treatable risk factor (RF) being the reduction of intraocular pressure (IOP). Genes associated with elevated IOP or POAG risk include: ABCA1, AFAP1, ARHGEF12, ATXN2, CAV1, CDKN2B-AS1, FOXC1, GAS7, GMDS, SIX1/SIX6, TMCO1, and TXNRD2. However, there are variations in RF and genetic factors based on ethnic and geographic differences; it is clear that unified molecular pathways accounting for POAG pathogenesis remain uncertain, although inflammation and senescence likely play an important role. There are similar ethnic and geographic complexities in primary angle closure glaucoma (PACG), but several genes have been associated with this disorder, including MMP9, HGF, HSP70, MFRP, and eNOS. In exfoliation glaucoma (XFG), genes implicated include LOXL1, CACNA1A, POMP, TMEM136, AGPAT1, RBMS3, and SEMA6A. Despite tremendous progress, major gaps remain in resolving the genetic architecture for the various glaucoma subtypes across ancestries. Large scale carefully designed studies are required to advance understanding of genetic loci as RF in glaucoma pathophysiology and to improve diagnosis and treatment options.

Update on the genetics of keratoconus

In the past few years we have seen a great acceleration of discoveries in the field of keratoconus including new treatments, diagnostic tools, genomic and molecular determinants of disease risk. Recent genome-wide association studies (GWAS) of keratoconus cases and population wide studies of variation in central corneal thickness and in corneal biomechanical properties confirmed already identified genes and found many new susceptibility variants and biological pathways. Recent findings in genetic determinants of familial keratoconus revealed functionally important variants and established first mouse model of keratoconus. Latest transcriptomic and expression studies started assessing novel non-coding RNA targets in addition to identifying tissue specific effects of coding genes. First genomic insights into better prediction of treatment outcomes are bringing the advent of genomic medicine into keratoconus clinical practice.

Fine-mapping and cell-specific enrichment at corneal resistance factor loci prioritize candidate causal regulatory variants

Corneal resistance factor (CRF) is altered during corneal diseases progression. Genome-wide-association studies (GWAS) indicated potential CRF and disease genetics overlap. Here, we characterise 135 CRF loci following GWAS in 76029 UK Biobank participants. Enrichment of extra-cellular matrix gene-sets, genetic correlation with corneal thickness (70% (SE = 5%)), reported keratoconus risk variants at 13 loci, all support relevance to corneal stroma biology. Fine-mapping identifies a subset of 55 highly likely causal variants, 91% of which are non-coding. Genomic features enrichments, using all associated variants, also indicate prominent regulatory causal role. We newly established open chromatin landscapes in two widely-used human cornea immortalised cell lines using ATAC-seq. Variants associated with CRF were significantly enriched in regulatory regions from the corneal stroma-derived cell line and enrichment increases to over 5 fold for variants prioritised by fine-mapping-including at GAS7, SMAD3 and COL6A1 loci. Our analysis generates many hypotheses for future functional validation of aetiological mechanisms.