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

Omics in Keratoconus: From Molecular to Clinical Practice

Keratoconus (KC) is a progressive ocular disorder marked by structural and functional alterations of the cornea, leading to significant visual impairment. Recent studies indicate that these corneal changes are linked to molecular and cellular mechanisms that disrupt and degrade the extracellular matrix. This degradation is influenced by proteinases that contribute to a loss of homeostasis and an imbalance in the antioxidant/oxidative state within the cornea, fostering oxidative stress, inflammation, and apoptosis. Although these biological processes have been identified primarily through molecular biology research, omics technologies have significantly advanced our understanding of the physiological and pathological phenomena associated with KC. Omics studies encompassing genomics, transcriptomics, proteomics, epigenomics, and metabolomics, have emerged as critical tools in elucidating the complex biological landscape of various diseases, including ocular conditions. The integrative application of these studies has demonstrated their potential in personalizing medicine across diverse fields such as oncology, neurology, and ophthalmology. This review aims to describe findings from omics research applied to keratoconus, highlighting the genomic, transcriptomic, proteomic, epigenomic, and metabolomic aspects derived from ocular and other biological samples. Notably, the molecular insights gained from these studies hold promise for identifying biomarkers of keratoconus, which could enhance diagnostic accuracy and therapeutic strategies. The exploration of these biomarkers may facilitate improved management and treatment options for patients, contributing to personalized care in keratoconus management.

Genotype-Phenotype Correlations in Corneal Dystrophies: Advances in Molecular Genetics and Therapeutic Insights

Corneal dystrophies are a group of predominantly rare inherited disorders. They are by definition bilateral, relatively symmetrical, and without systemic involvement, affecting corneal transparency and/or refraction. Traditional classification of corneal dystrophies is based on slit-lamp appearance, affected corneal layer and histological features. Molecular genetics has provided ultimate proof for the existence of distinct corneal dystrophies and discarded duplicates in their terminology. Currently, there are at least 16 genes with identified pathogenic variants implicated in corneal dystrophies. Herein, we summarise contemporary knowledge on genotype-phenotype correlations of corneal dystrophies, including a critical review of some reported variants, along with the understanding of the underlying pathogenic dystrophic process; essential knowledge for the development of targeted therapies.

Genome-wide association studies, Polygenic Risk Scores and Mendelian randomisation: an overview of common genetic epidemiology methods for ophthalmic clinicians

Genetic information will be increasingly integrated into clinical eye care within the current generation of ophthalmologists. For monogenic diseases such as retinoblastoma, genetic studies have been relatively straightforward as these conditions result from pathogenic variants in a single gene resulting in large physiological effects. However, most eye diseases result from the cumulative effects of multiple genetic variants and environmental factors. In such diseases, because each variant usually has an individually small effect, genetic studies for complex diseases are comparatively more challenging. This article aims to provide an overview of three genetic epidemiology methods for polygenic (or complex) diseases: genome-wide association studies (GWAS), Polygenic Risk Scores (PRS) and Mendelian randomisation (MR). A GWAS systematically conducts association analyses of a trait of interest against millions of genetic variants, usually in the form of single nucleotide polymorphisms, across the genome. GWAS findings can then be used for PRS construction and MR analyses. To construct a PRS, the cumulative effect of many genetic variants associated with a trait from a prior GWAS is calculated and taken as a quantitative representation of an individual's genetic risk of a complex disease. MR studies analyse an outcome measure against the genetic variants of an exposure, and are particularly useful in investigating causal relations between two traits where randomised controlled trials are not possible or ethical. In addition to explaining the principles of these three genetic epidemiology concepts, this article provides a minimally technical description of their basic methodology that is accessible to the non-expert reader.

Pilot Programme for Keratoconus Screening and Management in Adolescents with Down Syndrome

Background/Objectives: To screen a group of adolescents with Down Syndrome (Trisomy 21) for keratoconus and assess the feasibility of setting up a national screening service. Methods: Twenty-seven patients with Down Syndrome between 9 and 18 years of age attended our pilot keratoconus screening clinic. We recorded demographics, medical history, risk factors, best-corrected distance visual acuity, clinical examination results and corneal tomography results. The presence of keratoconus was confirmed by one of three corneal specialists based on clinical and tomographic findings. Tomographic analysis included zonal Kmax, thinnest point, inferior-superior asymmetry (IS Values), Belin/Ambrosio deviation value (BAD-D) and anterior and posterior elevation maps. Results: Early keratoconus was detected on tomography in 8 out of 54 eyes (15%) at the first review. These eyes were listed for crosslinking. The mean age of diagnosis was 14.6. Corneas in the Down Syndrome screening group were thinner and steeper (mean central corneal thickness (CCT) 479 µm vs. 536 µm and mean Kmax 49.2D vs. 45.8D, respectively) than healthy, age-matched controls from the literature. Conclusions: Fifteen percent of eyes (5 out of 27 patients) screened had tomographic evidence of keratoconus requiring treatment at their first review. We found an increased incidence of keratoconus in European individuals with Down Syndrome. Screening this vulnerable, high-risk population with corneal tomography can diagnose early keratoconus and enable corneal crosslinking to safely and effectively stabilise the disease. We advocate tomographic keratoconus screening for individuals with Down Syndrome in their mid-teens.

Genetic Estimates of Correlation and Causality Between Keratoconus and Osteoarthritis

PURPOSE

The relationship between keratoconus (KC) and osteoarthritis (OA) has been hypothesized, yet the causal link remained obscure. We aimed to evaluate the genetic correlation and causality between KC and OA through a genetic perspective.

METHODS

Linkage Disequilibrium Score Regression and bidirectional 2-sample Mendelian randomization (MR) analysis were conducted. Data were used from genome-wide association study on KC (4669 cases and 116,547 controls), OA (encompassing 24,955 patients with knee OA, 15,704 patients with hip OA, 39,427 patients with knee or hip OA, and 378,169 control participants), and KC data in the FinnGen consortium for replication and meta-analysis. Simple model, weighted model, inverse-variance weighted, weighted median, and MR-Egger regression were employed to assess the causal effects. MR pleiotropy residual sum and outlier method, MR-Egger method, and leave-one-out analysis were also applied to detect pleiotropy.

RESULTS

No statistically significant genetic correlation between KC and OA (all P > 0.05) were observed. MR estimates obtained from all 5 methods after meta-analysis revealed nonsignificant odds ratio values of KC-related traits to knee OA, hip OA, and OA (all adjusted P > 0.05). No evidence of causal relationships of knee and hip OA with KC-related traits was detected after meta-analysis of all 5 MR methods (all adjusted P > 0.05).

CONCLUSIONS

The large MR analysis indicated that KC may not be causally associated with a risk of OA.

Looking to the Future of Viral Vectors in Ocular Gene Therapy: Clinical Review

Eye diseases can significantly affect the quality of life of patients due to decreased visual acuity. Although modern ophthalmological diagnostic methods exist, some diseases of the visual system are asymptomatic in the early stages. Most patients seek advice from an ophthalmologist as a result of rapidly progressive manifestation of symptoms. A number of inherited and acquired eye diseases have only supportive treatment without eliminating the etiologic factor. A promising solution to this problem may be gene therapy, which has proven efficacy and safety shown in a number of clinical studies. By directly altering or replacing defective genes, this therapeutic approach will stop as well as reverse the progression of eye diseases. This review examines the concept of gene therapy and its application in the field of ocular pathologies, emphasizing the most recent scientific advances and their potential impacts on visual function status.

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.

Evaluation of Scleral Thickness in Patients With Keratoconus

OBJECTIVES

To evaluate scleral thickness of patients with keratoconus by anterior segment optical coherence tomography (AS-OCT).

MATERIALS AND METHODS

Fifty-two eyes of 42 patients with keratoconus (group 1) and 42 right eyes of 42 healthy individuals (group 2) were included. Scleral thickness measurements were taken with AS-OCT 6, 4, and 2 mm behind the scleral spur in four gaze positions: superior, inferior, temporal, and nasal. The data including central corneal thickness (CCT) by Pentacam from patients with keratoconus was also analyzed and its relationship with scleral thickness was evaluated.

RESULTS

In group 1, mean scleral thickness from 6 mm posterior to scleral spur in superior, inferior, temporal, and nasal quadrants were 397±35.3 (320-482), 408±29.5 (341-467), 393±40.8 (297-469), and 403±35.8 (310-482) μm, respectively. In group 2, mean scleral thickness from 6 mm posterior to scleral spur in superior, inferior, temporal, and nasal quadrants were 440±42.2 (346-552), 461±33.9 (390-563), 451±46.6 (340-549) and 448±34.3 (350-521) μm, respectively. The mean scleral thickness in group 1 was significantly lower in all quadrants compared with group 2 ( P <0.001). The mean CCT in group 1 was 478.22±87.55 (306-632) μm. Positive correlation was found between scleral thickness measurements taken at 4 mm (ϱ=0.312, P =0.037) and 2 mm (ϱ=0.308, P =0.039) posterior from the scleral spur in the inferior quadrant and CCT.

CONCLUSION

Scleral thickness was significantly lower in patients with keratoconus. Scleral thickness and CCT were positively correlated in inferior quadrant. Even though it is not yet known whether scleral changes are progressive or not, scleral thickness measurement may also find a place in the diagnosis and follow-up of patients with keratoconus in the future.

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.

Allergic disease and keratoconus: A two-sample univariable and multivariable Mendelian randomization study

Background

There is accumulating evidence that allergy is a risk factor for keratoconus. Nonetheless the association between allergic disease and keratoconus remains controversial. We performed a two-sample Mendelian randomization (MR) study to determine the putative causal association of 4 allergic diseases (allergic conjunctivitis, allergic asthma, allergic rhinitis and atopic dermatitis) with keratoconus.

Methods

Summary statistics were obtained from genome-wide association studies (GWAS) of allergic conjunctivitis (AC) (20,958 cases and 356,319 controls), allergic asthma (AA) (9631 cases and 210,122 controls), allergic rhinitis (AR) (11,009 cases and 359,149 controls), atopic dermatitis (AD) (13,473 cases and 336,589 controls), keratoconus (KC) (2116 cases and 24,626 controls) and 91 circulating inflammatory cytokines (n = 14,824). Two-sample univariable and multivariable MR analyses were performed. A two-step MR was then applied to determine whether systemic inflammatory cytokines mediated the effect of allergic disease on keratoconus.

Results

The causal odds ratio (OR) estimate of genetically determined KC was 1.66 (95% CI: 1.32-2.08; P < 0.001) for AC, 1.29 (95% CI: 1.10-1.51, P = 0.0014) for AA, 1.39 (95% CI: 1.15-1.68; P < 0.001) for AR and 1.30 (95% CI: 1.17-1.45, P < 0.001) for AD. Multivariable MR indicated a suggestive association between AC and KC after conditioning on other allergic diseases (OR 1.61; 95% CI: 1.10-2.34; P adjusted = 0.054). Two-step MR revealed that the effect was not mediated by systemic inflammatory cytokines.

Conclusions

Our findings provide evidence of a potential causal relationship between AC and KC. The effect of AC on KC may be mediated via other systemic inflammatory cytokines not included in the present study, or by alternative mechanisms. These findings may offer insight for prevention and intervention strategies to lower the risk of KC in patients with AC.

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.

Causal Association Between Atopic Dermatitis and Keratoconus: A Mendelian Randomization Study

Purpose

Although many studies have indicated that atopic dermatitis (AD) could contribute to the risk of keratoconus (KC), the causality between AD and KC remains controversial. This study aimed to explore the potential causal associations between AD and KC.

Methods

Instrumental variables for both exposures and outcomes were obtained from large-scale genome-wide association study summary statistics from previous meta-analyses. Mendelian randomization (MR) was applied to infer causal associations between AD and KC. Our main analyses were conducted by inverse-variance weighted (IVW) method multiplicative random effect model, complemented with additional five models and sensitivity analyses. Reverse MR analysis was applied to determine the direction of the causal association between AD and KC.

Results

Both IVW and weighted median methods revealed a causal effect of AD on KC (IVW odds ratio [OR], 1.475; P = 4.16 × 10-4; weighted median OR, 1.351; P = 7.65 × 10-3). The weighted mode, simple mode, and MR Egger methods demonstrated consistent direction of causality. Evidence from all sensitivity analyses further supported these associations. Reverse MR analyses did not suggest causal effects of KC on AD.

Conclusions

This study supported a significant causal effect of AD on KC, and reverse MR analysis proved that the causal association was unilateral.

Translational Relevance

This study provides valid evidence that regular ophthalmic examinations are recommended for patients with AD to detect and prevent KC at an early stage.

Polygenic Risk Scores and Genetically Complex Eye Disease

The success of genome-wide association studies (GWASs) in uncovering genetic variants associated with complex eye diseases has paved the way for the development of risk prediction approaches based on disease genetics. Derived from GWAS data, polygenic risk scores (PRSs) have been emerging as a promising indicator of an individual's genetic liability to disease. In this review, we recap the current progress of PRS development and utility across a range of common eye diseases. While illustrating the prediction accuracy of PRSs and their valuable role in risk stratification for certain eye diseases, we also address PRSs' uncertain implementation in clinical settings at this stage, particularly in circumstances where limited treatment options are available. Finally, we discuss obstacles in translating PRSs into practice, including barriers to clinical impact, issues when working with different ancestry groups, and communicating risk scores, as well as projections for future improvements.

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 Links Between Systemic Disorders and Keratoconus in European Population

PURPOSE

To establish the presence of a causal linkage between prevalent systemic diseases and keratoconus (KC).

DESIGN

Mendelian randomization (MR) analysis.

METHODS

After an exhaustive screening process, genetic variants linked to various systemic diseases were identified as instrumental variables at the genome-wide significance level. Subsequently, MR analyses were conducted to elucidate their potential causal connection with KC (N = 26,742). The encompassed systemic ailments comprise diabetes, hay fever/allergic rhinitis/eczema, obstructive sleep apnea, thyroid dysfunction, aortic aneurysm, major depressive disorder, inflammatory bowel disease (including Crohn's disease and ulcerative colitis), and mitral valve prolapse. Our study adheres to the principles of Strengthening the Reporting of Observational Studies in Epidemiology Using MR guidelines.

RESULTS

Using inverse variance weighting as the primary MR analysis method, our findings revealed that hay fever/allergic rhinitis/eczema (odds ratio, 10.144; 95% CI, 2.441-42.149; P = .001) and ulcerative colitis (odds ratio, 1.147; 95% CI, 1.054-1.248; P = .002) were associated with an increased risk of KC within the largest population under scrutiny. Conversely, the prolonged hyperglycemic state did not exhibit a potentially protective effect in delaying the pathogenesis of KC, and no correlation was observed between the two (odds ratio, 0.320; 95% CI, 0.029-3.549; P = .353). Also, obstructive sleep apnea, thyroid function, aortic aneurysm, major depressive disorder, Crohn's disease, and mitral valve prolapse did not exhibit a causal association with KC (P > .05 for all comparisons).

CONCLUSIONS

This study indicates an increased risk of KC related to hay fever/allergic rhinitis/eczema and ulcerative colitis, with diabetes not providing a protective effect. These findings may potentially contribute some insights to inform clinical interventions.

Targeted next-generation sequencing analysis in Italian patients with keratoconus

OBJECTIVE

To report variants in 26 candidate genes and describe the clinical features of Italian patients with keratoconus (KC).

SUBJECTS/METHODS

Sixty-four patients with a confirmed diagnosis of KC were enrolled in this genetic association study. Patients were classified into two study groups according to whether they had a confirmed diagnosis of progressive or stable KC. A purpose-developed Next Generation Sequencing (NGS) panel was used to identify and analyse the coding exons and flanking exon/intron boundaries of 26 genes known to be associated with KC and corneal dystrophies. Interpretation of the pathogenic significance of variants was performed using in silico predictive algorithms.

RESULT

The targeted NGS research identified a total of 167 allelic variants of 22 genes in the study population; twenty-four patients had stable keratoconus (n. 54 variants) and forty patients had progressive disease (n. 113 variants). We identified genetic variants of certain pathogenic significance in five patients with progressive KC; in addition, eight novel genetic variants were found in eight patients with progressive KC. Mutations of FLG, LOXHD1, ZNF469, and DOCK9 genes were twice more frequently identified in patients with progressive than stable disease. Filaggrin gene variants were found in 49 patients (76% of total), of whom 32 patients (80% of progressive KC group) had progressive disease.

CONCLUSIONS

Targeted NGS research provided new insights into the causative effect of candidate genes in the clinical phenotype of keratoconus. Filaggrin mutations were found to represent a genetic risk factor for development of progressive disease in Italy.

Mendelian randomization analysis does not reveal a causal influence between keratoconus and three major mental disorders

Background

Observational studies have suggested at a possible link between keratoconus (KC) and various mental disorders, but the exact direction of causation in these associations remains unclear. This study aims to investigate the potential causal link between KC and three prominent mental conditions: Anxiety, Depression, and Schizophrenia.

Methods

Using instrumental variables identified from Genome-wide association study (GWAS) data of European individuals, we conducted bidirectional two-sample Mendelian Randomization (MR) analyses to explore potential causal relationships between KC and the three major mental disorders. We primarily employed the Inverse-Variance Weighted (IVW) method to evaluate causality. In addition, we performed four supplementary MR methods (MR-Egger, Weighted Median, Simple Mode, and Weighted Mode). Furthermore, we conducted various sensitivity analyses to assess heterogeneity, horizontal pleiotropy, and result stability.

Results

Our findings did not reveal any concrete evidence of a causal link between KC and the three major mental disorders, namely anxiety, depression, and schizophrenia [anxiety: odds ratio (OR)=0.997, 95% confidence interval (CI)=0.988-1.008, p = 0.621; depression: OR=1.008, 95% CI=0.999-1.017, p = 0.084; schizophrenia: OR=1.002, 95% CI= 0.984-1.020, p = 0.840]. Similarly, the three major mental disorders were not caustically associated with KC [anxiety: OR=1.014, 95% CI=0.635-1.620, p = 0.953; depression: OR=1.109, 95% CI= 0.749-1.643, p = 0.604; schizophrenia: OR= 0.969, 95% CI= 0.884-1.062, p = 0.497]. The sensitivity analyses indicated that the results remained robust, with no signs of pleiotropy or heterogeneity.

Conclusions

Our study does not support a genetically determined significant causal connection between KC and the three major mental disorders. The increased occurrence of mental disorders observed in KC patients in observational reports likely arises from factors that can be modified. Further research is warranted to unveil the underlying mechanisms behind the associations observed in observational studies.

Gene‒Environment Interaction Between CAST Gene and Eye-Rubbing in the Chinese Keratoconus Cohort Study: A Case-Only Study

Purpose

Keratoconus (KC), characterized by progressive corneal protrusion and thinning, is a complex disease influenced by the combination of genetic and environmental factors. The purpose of this study was to explore potential gene‒environment interaction between the calpastatin (CAST) gene and eye-rubbing in KC.

Methods

A case-only study including 930 patients (676 patients with eye-rubbing and 254 patients without eye-rubbing) from the Chinese Keratoconus (CKC) cohort study was performed in the present study. Genotyping of single nucleotide polymorphism (SNP) was conducted using the Illumina Infinium Human Asian Screening Array (ASA) Beadchip. The gene‒environment interactions between CAST gene and eye-rubbing were analyzed using PLINK version 1.90. The interactions between CAST genotypes and eye-rubbing were analyzed by logistic regression models. The SNP-SNP-environment interactions were analyzed using generalized multifactor dimensionality reduction (GMDR).

Results

Three SNPs in CAST gene, namely, rs26515, rs27991, and rs9314177, reached the significance threshold for interactions (defined as P < 2.272 × 10-3). Notably, the minor alleles of these three SNPs exhibited negative interactions with eye-rubbing in KC. The results of logistic regression models revealed that the minor allele homozygotes and heterozygotes of rs26515, rs27991, and rs9314177 also exhibited negative interactions with eye-rubbing. Furthermore, GMDR analysis revealed the significant SNP-SNP-environment interactions among rs26515, rs27991, rs9314177, and eye-rubbing in KC.

Conclusions

This study identified rs26515, rs27991, and rs9314177 in CAST gene existed gene-environment interactions with eye-rubbing in KC, which is highly important for understanding the underlying biological mechanisms of KC and guiding precision prevention and proper management.

Whole-exome sequencing screening for candidate genes and variants associated with primary sporadic keratoconus in Chinese patients

The pathogenesis of keratoconus (KC) is complex, and genetic factors play an important role. The purpose of this study was to screen and analyse candidate genes and variants in Chinese patients with primary sporadic KC. Whole-exome sequencing (WES) was performed to identify candidate genes and variants in 105 unrelated Chinese patients with primary sporadic KC. Through a series of screening processes, 54 candidate variants in 26 KC candidate genes were identified in 53 KC patients (53/105, 50.5%). These 54 candidate variants included 10 previously identified variants in 9 KC candidate genes and 44 novel variants in 20 KC candidate genes. The previously identified variants occurred in 25.7% (27/105) of patients. Of these, 4 variants (COL6A5, c.5014T > G; CAST, c.1814G > A; ZNF469, c.946G > A; and MPDZ, c.3836A > G) were identified for the first time in Chinese KC patients. The novel variants occurred in 33.3% (35/105) of patients. Of the 26 screened KC candidate genes, 11 KC candidate genes (CAT, COL12A1, FLG, HKDC1, HSPG2, PLOD1, ITGA2, TFAP2B, USH2A, WNT10A, and COL6A5) were found to be potentially pathogenic in Chinese KC patients for the first time. Gene Ontology (GO) biological process (BP) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were performed on the 26 KC candidate genes using the Database for Annotation, Visualization, and Integrated Discovery (DAVID). The results showed that the KC candidate genes were significantly enriched in biological processes such as collagen fibril organization and extracellular matrix (ECM) organization and in ECM-receptor interaction and protein digestion and absorption pathways. The results further expand the spectrum of KC candidate variants and provide a basis for further KC gene studies.

Inflammation and keratoconus: A comprehensive bidirectional Mendelian randomization analysis

An increasing body of evidence supports the involvement of inflammation and immune responses in the occurrence and development of keratoconus (KC). However, the causal relationship between inflammatory factors and KC remains unclear. We employed a 2-way Mendelian randomization (MR) approach to investigate the interaction between KC and inflammatory factors. Instrumental variables for 41 circulating inflammatory regulators and 12 matrix metalloproteinases (MMPs) were selected from genome-wide association studies of European ancestry. Summary statistics for KC were obtained from a genome-wide association study comprising 2116 cases and 24,626 controls of European ancestry. The primary analytical method for assessing causality was the inverse-variance weighted method. Two additional MR methods (MR-Egger and weighted median) were employed to complement the inverse-variance weighted results. In addition, several sensitivity analyses were conducted to evaluate heterogeneity, horizontal pleiotropy, and stability. Our findings indicated that genetically predicted higher levels of macrophage inflammatory protein-1β (odds ratio = 1.126, 95% confidence interval: 1.029-1.232, P = .01) and MMP-13 (odds ratio = 1.211, 95% confidence interval: 1.070-1.371, P = .003) were positively associated with an elevated risk of KC. Conversely, genetically predicted KC was associated with increased levels of interferon-gamma, interleukin-4, and MMP-1. Our current study provided suggestive evidence supporting causal associations of macrophage inflammatory protein-1β and MMP-13 with the risk of KC. In addition, KC appeared to affect the expression of interferon-gamma, interleukin-4, and MMP-1.

Systemic and Ocular Associations of Keratoconus

Introduction

Keratoconus (KC) is the most prevalent corneal ectasia in the world and its pathogenesis is influenced by both ocular and systemic factors. This review explores the multifaceted associations between keratoconus and systemic health conditions, ocular characteristics, and various other environmental/exogenous factors, aiming to illuminate how these relationships influence the pathophysiology of the disease.

Areas Covered

This review will summarize the fundamental attributes of KC, review and discuss the systemic and ocular association of KC including molecular biomarkers, and provide an organized overview of the parallel alterations occurring within various biological pathways in KC.

Expert Opinion

Despite the substantial volume of research on keratoconus, the precise etiology of the disease remains elusive. Further studies are necessary to deepen our understanding of this intricate disorder and improve its management.

Identification of hub genes and molecular pathways in keratoconus by integrating bioinformatics and literature mining at the RNA level

OBJECTIVE

Keratoconus (KC) is a condition characterized by progressive corneal steepening and thinning. However, its pathophysiological mechanism remains vague. We mainly performed literature mining to extract bioinformatic and related data on KC at the RNA level. The objective of this study was to explore the potential pathological mechanisms of KC by identifying hub genes and key molecular pathways at the RNA level.

METHODS

We performed an exhaustive search of the PubMed database and identified studies that pertained to gene transcripts derived from diverse corneal layers in patients with KC. The identified differentially expressed genes were intersected, and overlapping genes were extracted for further analyses. Significantly enriched genes were screened using "Gene Ontology" (GO) and "Kyoto Encyclopedia of Genes and Genomes" (KEGG) analysis with the "Database for Annotation, Visualization, and Integrated Discovery" (DAVID) database. A protein-protein interaction (PPI) network was constructed for the significantly enriched genes using the STRING database. The PPI network was visualized using the Cytoscape software, and hub genes were screened via betweenness centrality values. Pathways that play a critical role in the pathophysiology of KC were discovered using the GO and KEGG analyses of the hub genes.

RESULTS

68 overlapping genes were obtained. Fifty genes were significantly enriched in 67 biological processes, and 16 genes were identified in 7 KEGG pathways. Moreover, 14 nodes and 32 edges were identified via the PPI network constructed using the STRING database. Multiple analyses identified 4 hub genes, 12 enriched biological processes, and 6 KEGG pathways. GO enrichment analysis showed that the hub genes are mainly involved in the positive regulation of apoptotic process, and KEGG analysis showed that the hub genes are primarily associated with the interleukin-17 (IL-17) and tumor necrosis factor (TNF) pathways. Overall, the matrix metalloproteinase 9, IL-6, estrogen receptor 1, and prostaglandin-endoperoxide synthase 2 were the potential important genes associated with KC.

CONCLUSION

Four genes, matrix metalloproteinase 9, IL-6, estrogen receptor 1, and prostaglandin endoperoxide synthase 2, as well as IL-17 and TNF pathways, are critical in the development of KC. Inflammation and apoptosis may contribute to the pathogenesis of KC.

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.

Association of Novel Loci With Keratoconus Susceptibility in a Chinese Genome-Wide Association Study

Purpose

Keratoconus (KC) is a progressive corneal disease that can lead to corneal blindness if not properly managed. The purpose of this study was to identify genetic associations with KC in China and to investigate whether these genetic variants are associated with corneal thickness and corneal curvature in KC cases.

Methods

A genome-wide association study was conducted on 853 patients with KC and 6248 controls. The KC cases were genotyped with the Illumina Infinium Human Asian Screening Array BeadChip, and the controls were genotyped with the Illumina Infinium Human Global Screening Array BeadChip. Genetic associations with KC, as well as correlations between the positive variants and corneal parameters including central corneal thickness (CCT) and mean keratometry (Km), were compared using PLINK version 1.90.

Results

Our present study identified four single-nucleotide polymorphisms (SNPs) within four risk loci (PTGER3: rs2300163, EYA1: rs1077435, ASS1: rs141365191, and CHTF8: rs3743680) associated with KC in Chinese patients that reached genome-wide significance. Among the identified SNPs with P < 1.00 × 10-4, seven SNPs (FOSL2-PLB1: rs12622211, RXRA-COL5A1: rs3118515, rs3132306, rs1536482, rs3118520, KAT6B: rs192187772, RAP2A-IPO5: rs41361245) were observed to be associated with CCT, and one SNP (USP13: rs6767552) was found to be associated with Km.

Conclusions

In the first genome-wide association study of KC with a relatively large study population in China, we identified four SNPs in four risk loci associated with the disease. The findings enriched the understanding of genetic susceptibility to KC and provided new insights into the genetic etiology of the disease.

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.

Collagen type XII is undetectable in keratoconus Bowman's layer

PURPOSE

Corneal biomechanical failure is the hallmark of keratoconus (KC); however, the cause of this failure remains elusive. Collagen type XII (COL12A1), which localises to Bowman's layer (BL), is thought to function in stress-bearing areas, such as BL. Given the putative protective role of COL12A1 in biomechanical stability, this study aims to characterise COL12A1 expression in all corneal layers involved in KC.

METHODS

TaqMan quantitative PCR was performed on 31 corneal epithelium samples of progressive KC and myopic control eyes. Tissue microarrays were constructed using full-thickness corneas from 61 KC cases during keratoplasty and 18 non-KC autopsy eyes and stained with an antibody specific to COL12A1. Additionally, COL12A1 was knocked out in vitro in immortalised HEK293 cells.

RESULTS

COL12A1 expression was reduced at transcript levels in KC epithelium compared with controls (ratio: 0.58, p<0.03). Immunohistochemical studies demonstrated that COL12A1 protein expression in BL was undetectable, with reduced expression in KC epithelium, basement membrane and stroma.

CONCLUSIONS

The apparent absence of COL12A1 in KC BL, together with the functional importance that COL12A1 is thought to have in stress bearing areas, suggests that COL12A1 may play a role in the pathogenesis of KC. Further studies are necessary to investigate the mechanisms that lead to COL12A1 dysregulation in KC.

Genetic variants in the FOXO1 and ZNF469 genes are associated with keratoconus in Sweden: a case-control study

BACKGROUND

Keratoconus (KC) is characterized by pathological thinning and bulging of the cornea that may lead to visual impairment. The etiology of sporadic KC remains enigmatic despite intensive research in recent decades. The purpose of this study was to examine the relationship between previously highlighted genetic variants associated with KC and sporadic KC in a Swedish cohort.

METHODS

A total of 176 patients (age 16-70 years) with sporadic KC diagnosed by Scheimpflug-topography (Pentacam) were included. The control group (n = 418; age 70 years) was a subsample originating from the Gothenburg H70 Birth Cohort Studies of ageing. Extraction of DNA from blood samples was performed according to standard procedures, and genotyping was performed using competitive allele specific PCR (KASP) technology. A total of 11 single nucleotide polymorphisms (SNPs) were selected for analysis.

RESULTS

Statistically significant associations (p = 0.005) were found between the SNPs rs2721051 and rs9938149 and sporadic KC. These results replicate earlier research that found associations between genetic variants in the FOXO1 and BANP-ZNF469 genes and sporadic KC in other populations.

CONCLUSION

Genetic variations in the FOXO1 and BANP-ZNF469 genes may be involved in the pathogenesis of sporadic KC.

Japanese 17q12 Deletion Syndrome with Complex Clinical Manifestations

17q12 deletion syndrome is a rare chromosomal anomaly with variable phenotypes, caused by the heterozygous deletion of chromosome 17q12. We herein report a 35-year-old Japanese patient with chromosomal 17q12 deletion syndrome identified by de novo deletion of the 1.46 Mb segment at the 17q12 band by genetic analyses. He exhibited a wide range of phenotypes, such as maturity-onset diabetes of the young (MODY) type 5, structural or functional abnormalities of the kidney, liver, and pancreas; facial dysmorphic features, electrolyte disorders; keratoconus, and acquired perforating dermatosis. This case report provides valuable resources concerning the clinical spectrum of rare 17q12 deletion syndrome.

Association of torpedo maculopathy and keratoconus in a young patient: A multimodal imaging study

PURPOSE

To report and document a case of torpedo maculopathy found in a patient affected by keratoconus.Case report: An healthy 16-year-old male patient, affected by keratoconus in both eyes, was referred to the cornea service of our hospital for a follow-up visit.During the dilated fundus examination of the left eye, an oval, well-demarcated, hypopigmented lesion was observed in the juxtafoveal temporal region, pointing towards the center of the macula. Multimodal imaging of the lesion was performed, and the diagnosis of Torpedo Maculopathy was established based on the clinical picture.

CONCLUSION

This is the first case of torpedo maculopathy described in a patient affected by keratoconus. This association may be merely fortuitous or the result of developmental abnormalities affecting both corneal and retinal structures.

Causal relationships between type 2 diabetes, glycemic traits and keratoconus

Purpose

The relationship between diabetes mellitus and keratoconus remains controversial. This study aimed to assess the potential causal relationships among type 2 diabetes, glycemic traits, and the risk of keratoconus.

Methods

We used a two-sample Mendelian randomization (MR) design based on genome-wide association summary statistics. Fasting glucose, proinsulin levels, adiponectin, hemoglobin A1c (HbA1c) and type 2 diabetes with and without body mass index (BMI) adjustment were used as exposures and keratoconus was used as the outcome. MR analysis was performed using the inverse-variance weighted method, MR-Egger regression method, weighted-mode method, weighted median method and the MR-pleiotropy residual sum and outlier test (PRESSO).

Results

Results showed that genetically predicted lower fasting glucose were significantly associated with a higher risk of keratoconus [IVW: odds ratio (OR) = 0.382; 95% confidence interval (CI) = 0.261-0.560; p = 8.162 × 10-7]. Genetically predicted lower proinsulin levels were potentially linked to a higher risk of keratoconus (IVW: OR = 0.739; 95% CI = 0.568-0.963; p = 0.025). In addition, genetically predicted type 2 diabetes negatively correlated with keratoconus (IVW: BMI-unadjusted: OR = 0.869; 95% CI = 0.775-0.974, p = 0.016; BMI-adjusted: OR = 0.880, 95% CI = 0.789-0.982, p = 0.022). These associations were further corroborated by the evidence from all sensitivity analyses.

Conclusion

These findings provide genetic evidence that higher fasting glucose levels are associated with a lower risk of keratoconus. However, further studies are required to confirmed this hypothesis and to understand the mechanisms underlying this putative causative relationship.

Identification of the regulatory circuit governing corneal epithelial fate determination and disease

The transparent corneal epithelium in the eye is maintained through the homeostasis regulated by limbal stem cells (LSCs), while the nontransparent epidermis relies on epidermal keratinocytes for renewal. Despite their cellular similarities, the precise cell fates of these two types of epithelial stem cells, which give rise to functionally distinct epithelia, remain unknown. We performed a multi-omics analysis of human LSCs from the cornea and keratinocytes from the epidermis and characterized their molecular signatures, highlighting their similarities and differences. Through gene regulatory network analyses, we identified shared and cell type-specific transcription factors (TFs) that define specific cell fates and established their regulatory hierarchy. Single-cell RNA-seq (scRNA-seq) analyses of the cornea and the epidermis confirmed these shared and cell type-specific TFs. Notably, the shared and LSC-specific TFs can cooperatively target genes associated with corneal opacity. Importantly, we discovered that FOSL2, a direct PAX6 target gene, is a novel candidate associated with corneal opacity, and it regulates genes implicated in corneal diseases. By characterizing molecular signatures, our study unveils the regulatory circuitry governing the LSC fate and its association with corneal opacity.

Contribution of common and rare variants to Asian neovascular age-related macular degeneration subtypes

Neovascular age-related macular degeneration (nAMD), along with its clinical subtype known as polypoidal choroidal vasculopathy (PCV), are among the leading causes of vision loss in elderly Asians. In a genome-wide association study (GWAS) comprising 3,128 nAMD (1,555 PCV and 1,573 typical nAMD), and 5,493 controls of East Asian ancestry, we identify twelve loci, of which four are novel ([Formula: see text]). Substantial genetic sharing between PCV and typical nAMD is noted (rg = 0.666), whereas collagen extracellular matrix and fibrosis-related pathways are more pronounced for PCV. Whole-exome sequencing in 259 PCV patients revealed functional rare variants burden in collagen type I alpha 1 chain gene (COL1A1; [Formula: see text]) and potential enrichment of functional rare mutations at AMD-associated loci. At the GATA binding protein 5 (GATA5) locus, the most significant GWAS novel loci, the expressions of genes including laminin subunit alpha 5 (Lama5), mitochondrial ribosome associated GTPase 2 (Mtg2), and collagen type IX alpha 3 chain (Col9A3), are significantly induced during retinal angiogenesis and subretinal fibrosis in murine models. Furthermore, retinoic acid increased the expression of LAMA5 and MTG2 in vitro. Taken together, our data provide insights into the genetic basis of AMD pathogenesis in the Asian population.

Phenotype and genotype of concurrent keratoconus and Fuchs endothelial corneal dystrophy

PURPOSE

To characterise the phenotype and genotype of concurrent keratoconus and Fuchs endothelial corneal dystrophy (KC + FECD).

METHODS

We recruited 20 patients with concurrent KC + FECD for a retrospective observational case series from the United Kingdom and the Czech Republic. We compared eight parameters of corneal shape (Pentacam, Oculus) with two groups of age-matched controls who had either isolated keratoconus (KC) or isolated FECD. We genotyped probands for an intronic triplet TCF4 repeat expansion (CTG18.1) and the ZEB1 variant c.1920G >T p.(Gln640His).

RESULTS

The median age at diagnosis of patients with KC + FECD was 54 (interquartile range 46 to 66) years, with no evidence of KC progression (median follow-up 84 months, range 12 to 120 months). The mean (standard deviation (SD)) of the minimum corneal thickness, 493 (62.7) μm, was greater than eyes with KC, 458 (51.1) μm, but less than eyes with FECD, 590 (55.6) μm. Seven other parameters of corneal shape were more like KC than FECD. Seven (35%) probands with KC + FECD had a TCF4 repeat expansion of ≥50 compared to five controls with isolated FECD. The average of the largest TCF4 expansion in cases with KC + FECD (46 repeats, SD 36 repeats) was similar to the age-matched controls with isolated FECD (36 repeats, SD 28 repeats; p = 0.299). No patient with KC + FECD harboured the ZEB1 variant.

CONCLUSIONS

The KC + FECD phenotype is consistent with KC but with superimposed stromal swelling from endothelial disease. The proportion of cases with a TCF4 expansion is similar in concurrent KC + FECD and age-matched controls with isolated FECD.

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.

Genetic Correlations Among Corneal Biophysical Parameters and Anthropometric Traits

Purpose

The genetic architecture of corneal dysfunction remains poorly understood. Epidemiological and clinical evidence suggests a relationship between corneal structural features and anthropometric measures. We used global and local genetic similarity analysis to identify genomic features that may underlie structural corneal dysfunction.

Methods

We assembled genome-wide association study summary statistics for corneal features (central corneal thickness, corneal hysteresis [CH], corneal resistance factor [CRF], and the 3 mm index of keratometry) and anthropometric traits (body mass index, weight, and height) in Europeans. We calculated global genetic correlations (rg) between traits using linkage disequilibrium (LD) score regression and local genetic covariance using ρ-HESS, which partitions the genome and performs regression with LD regions. Finally, we identified genes located within regions of significant genetic covariance and analyzed patterns of tissue expression and pathway enrichment.

Results

Global LD score regression revealed significant negative correlations between height and both CH (rg = -0.12; P = 2.0 × 10-7) and CRF (rg = -0.11; P = 6.9 × 10-7). Local analysis revealed 68 genomic regions exhibiting significant local genetic covariance between CRF and height, containing 2874 unique genes. Pathway analysis of genes in regions with significant local rg revealed enrichment among signaling pathways with known keratoconus associations, including cadherin and Wnt signaling, as well as enrichment of genes modulated by copper and zinc ions.

Conclusions

Corneal biophysical parameters and height share a common genomic architecture, which may facilitate identification of disease-associated genes and therapies for corneal ectasias.

Translational Relevance

Local genetic covariance analysis enables the identification of associated genes and therapeutic targets for corneal ectatic disease.

Keratoconus International Consortium (KIC)- advancing keratoconus research

CLINICAL RELEVANCE

The Keratoconus International Consortium (KIC) will allow better understanding of keratoconus.

BACKGROUND

Keratoconus is a disorder characterised by corneal elevation and thinning, leading to reduced vision. The current gaps in understanding of this disease will be discussed and the need for a multi-pronged and multi-centre engagement to enhance our understanding of keratoconus will be highlighted.

DESIGN

KIC has been established to address the gaps in our understanding of keratoconus with the aim of collecting baseline as well as longitudinal data on several fields.

PARTICIPANTS

Keratoconus and control (no corneal condition) subjects from different sites globally will be recruited in the study.

METHODS

KIC collects data using an online, secure database, which enables standardised data collection at member sites. Data fields collected include medical history, clinical features, quality of life and economic burden questionnaires and possible genetic sample collection from patients of different ethnicities across different geographical locations.

RESULTS

There are currently 40 Australian and international clinics or hospital departments who have joined the KIC. Baseline data has so far been collected on 1130 keratoconus patients and indicates a median age of 29.70 years with 61% being male. A total of 15.3% report a positive family history of keratoconus and 57.7% self-report a history of frequent eye rubbing.

CONCLUSION

The strength of this consortium is its international, collaborative design and use of a common data collection tool. Inclusion and analyses of cross-sectional and longitudinal data will help answer many questions that remain in keratoconus, including factors affecting progression and treatment outcomes.

The association between VSX1 exon3 gene variants and keratoconus in Malaysian patients

Purpose

This case-control study aims to examine possible associations of VSX1 exon3 gene variants with the development of keratoconus (KC) in Malaysian patients.

Methods

A case-control study was done on 42 keratoconus cases, 127 family member controls, and 96 normal controls.

Results

Three gene variants, p.A182A, p.P237P, and p.R217H showed significant associations with keratoconus (P < 0.05). While p.A182A and p.P227P were more prevalent than in the family and normal controls (OR 3.14-4.05), the reverse was observed with p.R217H (OR 0.086-1.59). With Haploview analysis, p.A182A and p.P237P were shown to be in linkage disequilibrium (LD) (LOD (logarithm of the odds score) score of 2.0, r2 of 0.957, and 95% confidence interval (CI) of 0.96-1.00).

Conclusion

The study results suggest that the p.A182A and p.P237P variants could have contributed to the development of keratoconus in some Malaysians and that these two variants are likely to be co-inherited. In contrast, the p.R217H variant appeared to confer some protection against the development of keratoconus.

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.

Bioinformatics analysis of key candidate genes and pathways in Chinese patients with keratoconus

Keratoconus (KC) is a multifactorial disease in which genetic factors played important roles in its pathogenesis. The purpose of the current study was to identify the key candidate genes and pathways in Chinese patients with KC through bioinformatics analysis. Totally, we identified 71 candidate genes by analyzing the results of whole exome sequencing on 51 Chinese patients with KC, combining with previous reports on differential expression at transcription and protein levels in KC. Gene enrichment analysis with GeneCodis demonstrated that two significantly enriched terms including 21 genes in biological process (BP) were detected, and six significantly enriched terms containing 14 genes in Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway were discovered. The STRING was utilized to construct the protein-protein interaction (PPI) network of identified genes. The result showed that a PPI network consisted of 14 nodes with 14 edges was constructed, and two gene modules were obtained. Eight hub genes (LAMB3, LAMA3, LAMA1, ITGA6, ITGA3, COL6A3, COL6A2, and COL6A1) were identified as key candidate genes for KC by cytoHubba in Cytoscape. Functional enrichment analysis with ClueGO and CluePedia indicated that the ECM-receptor interaction was the key pathway accounted for KC. The findings might provide novel insights on the genetic basis of KC.

Intergenic variants, rs1200114 and rs1200108 are genetically associated along with a decreased ATP1B1 expression in Fuchs Endothelial Corneal Dystrophy

Fuchs endothelial corneal dystrophy (FECD) is an age-related, bilateral corneal condition, characterized by apoptosis of the terminally differentiated endothelial cells. A genome-wide association study (GWAS) conducted in the European population in 2017, identified a new single nucleotide polymorphism (SNP), rs1200114 in the intergenic region between long intergenic non-protein coding RNA 970 (LINC00970) and ATPase Na+/K+ transporting subunit beta 1 (ATP1B1). The major focus of the current study is to understand the genetic association of this intergenic variant, rs1200114 with FECD in the Indian population. Sanger sequencing followed by statistical analysis indicated a significant difference in the allelic frequency between controls and cases (P = 0.01) with the minor allele 'G' of rs1200114 imparting a 1.64 fold increased risk for the disease. Luciferase reporter assay revealed no significant difference in the luciferase activity between allele 'A' and 'G' of rs1200114. However, quantitative RT-PCR assay revealed lower expression of ATP1B1 in FECD subjects compared with controls (P = 0.007). Therefore, to find whether another nearby SNP imparts regulatory effect, tag SNP association analysis was carried out; which revealed a significant association of another SNP, rs1200108, present in the intergenic region between LINC00970 and ATP1B1 with FECD (P = 0.009). The protective allele 'A' of rs1200108 displayed reduced reporter activity as opposed to the risk allele 'G' (P = 0.014). Furthermore, haplotype 'A-A' of rs1200108 - rs1200114 was present at a higher frequency in control subjects, suggesting it as a protective haplotype. Altogether, this study inferred the genetic association of rs1200114 and rs1200108 along with the decreased expression of ATP1B1 related to FECD pathogenesis in the Indian population.

Cytoprotective Effects of Human Platelet Lysate during the Xeno-Free Culture of Human Donor Corneas

We evaluated the suitability of 2% human platelet lysate medium (2%HPL) as a replacement for 2% fetal bovine serum medium (2%FBS) for the xeno-free organ culture of human donor corneas. A total of 32 corneas from 16 human donors were cultured in 2%FBS for 3 days (TP1), then evaluated using phase contrast microscopy (endothelial cell density (ECD) and cell morphology). Following an additional 25-day culture period (TP2) in either 2%FBS or 2%HPL, the pairs were again compared using microscopy; then stroma and Descemet membrane/endothelium (DmE) were processed for next generation sequencing (NGS). At TP2 the ECD was higher in the 2%HPL group (2179 ± 288 cells/mm²) compared to 2%FBS (2113 ± 331 cells/mm²; p = 0.03), and endothelial cell loss was lower (ECL HPL = -0.7% vs. FBS = -3.8%; p = 0.01). There were no significant differences in cell morphology between TP1 and 2, or between 2%HPL and 2%FBS. NGS showed the differential expression of 1644 genes in endothelial cells and 217 genes in stromal cells. It was found that 2%HPL led to the upregulation of cytoprotective, anti-inflammatory and anti-fibrotic genes (HMOX1, SERPINE1, ANGPTL4, LEFTY2, GADD45B, PLIN2, PTX3, GFRA1/2), and the downregulation of pro-inflammatory/apoptotic genes (e.g., CXCL14, SIK1B, PLK5, PPP2R3B, FABP5, MAL, GATA3). 2%HPL is a suitable xeno-free substitution for 2%FBS in human cornea organ culture, inducing less ECL and producing potentially beneficial alterations in gene expression.

Big data in corneal diseases and cataract: Current applications and future directions

The accelerated growth in electronic health records (EHR), Internet-of-Things, mHealth, telemedicine, and artificial intelligence (AI) in the recent years have significantly fuelled the interest and development in big data research. Big data refer to complex datasets that are characterized by the attributes of "5 Vs"-variety, volume, velocity, veracity, and value. Big data analytics research has so far benefitted many fields of medicine, including ophthalmology. The availability of these big data not only allow for comprehensive and timely examinations of the epidemiology, trends, characteristics, outcomes, and prognostic factors of many diseases, but also enable the development of highly accurate AI algorithms in diagnosing a wide range of medical diseases as well as discovering new patterns or associations of diseases that are previously unknown to clinicians and researchers. Within the field of ophthalmology, there is a rapidly expanding pool of large clinical registries, epidemiological studies, omics studies, and biobanks through which big data can be accessed. National corneal transplant registries, genome-wide association studies, national cataract databases, and large ophthalmology-related EHR-based registries (e.g., AAO IRIS Registry) are some of the key resources. In this review, we aim to provide a succinct overview of the availability and clinical applicability of big data in ophthalmology, particularly from the perspective of corneal diseases and cataract, the synergistic potential of big data, AI technologies, internet of things, mHealth, and wearable smart devices, and the potential barriers for realizing the clinical and research potential of big data in this field.

Clinical-epidemiological characteristics of keratoconus in Asturias

INTRODUCTION

Keratoconus is the most frequent corneal ectasia and second most common cause of keratoplasty worldwide. The clinical and epidemiological characteristics of keratoconus have rarely been studied in Spain, and to date, no study has been done in Asturias.

MATERIAL AND METHODS

A cross-sectional study was conducted to analyze all cases of keratoconus diagnosed in the first consultations of the Cornea and Surface section of the Hospital Universitario Central de Asturias (Oviedo) and Hospital Universitario de Cabueñes (Gijón) between 2017 and 2020. Demographic data, visual acuity and corneal topography were collected.

RESULTS

A total of 112 patients (42.0% women and 58.0% men) were included. The mean age was 36.84 ± 15.59, with 39.3% being older than 40 years. 31.8%, 28.0%, 15.0%, and 25.2% were in stages I, II, III and IV, Amsler-Krumeich scale, respectively. Corneal segments were implanted in 7.5% of patients, cross-linking was performed in 4.7%, and keratoplasty was performed in 13.2%.

CONCLUSION

Approximately 40% of keratoconus cases are detected at advanced stages and another 40% are detected at an age over 40 years. 13% of these patients require keratoplasty, which is why it is frequently performed in our setting. Therefore, early diagnosis before irreversible visual loss is necessary.

Trio-based exome sequencing broaden the genetic spectrum in keratoconus

Keratoconus (KC) is a complex corneal disorder with genetic factors involving in its pathogenesis. The genetic etiology of KC has not been fully elucidated. In this study, we aimed to expand the genetic spectrum in KC by trio-based exome sequencing. Trio-based exome sequencing was conducted in 20 patients with KC and their unaffected parents to broaden the genetic spectrum of the disease. With a series of filtering criteria, de novo, recessive homozygous, and compound heterozygous variants in candidate genes were identified, and the candidate genes were classified for further analysis. Finally, we identified 60 variants in 32 candidate genes through trio-based exome sequencing. Among the candidate genes, 10 genes (ARHGEF10, ARHGEF17, ASPM, FLNA, NDRG1, NEB, PLS3, STARD8, SYNE1, TTN) were classified as cytoskeleton-related genes, 4 genes (COL28A1, SDK1, STAB1, TENM2) were classified as cell adhesion-related genes, and 18 genes (APLP2, BCORL1, CCNB3, FOXN1, FUT8, GALNT10, HEPH, HHIP, HMGB3, HS6ST2, JADE3, KIAA0040, MCF2L, MYOF, QRICH2, RPS6KA6, SMARCA1, TNRC6A) were classified into other genes group. Additionally, the candidate rare deleterious variants in TTN were highly repeated in 25% trios. In conclusion, the study provided new insights into the genetic spectrum of KC which might underlie the genetic etiology for the disease. The findings would improve our understanding of pathogenesis in KC and provide critical clues to future functional validation.

Cataract surgery in patients with underlying keratoconus: focused review

An underlying diagnosis of keratoconus (KC) can complicate cataract surgery. In this study, the results of a focused review of the literature pertaining to cataract surgery in patients with KC are detailed. Topics essential for the appropriate management of this patient population are discussed. First, the individual and shared epidemiology and pathophysiology of cataract and KC are reviewed. Then, the theory and approach to intraocular lens power calculation are discussed, highlighting particularities and pitfalls of this exercise when performed in patients with KC. Finally, several special-although not uncommon-management scenarios and questions are addressed, such as surgical planning in cases where corneal stabilization or tissue replacement interventions are also necessitated.

Association of 2 Lysyl Oxidase Gene Single Nucleotide Polymorphisms with Keratoconus: A Nationwide Registration Study

Purpose

Keratoconus (KC) is the most common primary ectatic corneal disease, characterized by progressive thinning of the cornea, affecting its shape and structure and leading to visual loss. Lysyl oxidase is an important component of the extracellular matrix and contributes to the homeostasis of corneal stromal extracellular matrix via enzymatic reaction. This nationwide registration study aims to examine the association of KC with 2 known single nucleotide polymorphisms, rs2956540 and rs10519694, in a population of Iranian descent.

Design

Case-control.

Participants

One hundred seventy-eight subjects with KC and 180 clinically healthy subjects participated in the study.

Methods

Genomic DNA was extracted from peripheral blood samples, and their genotypes were determined using tetra-primer amplification refractory mutation system-polymerase chain reaction.

Main Outcome Measures

Allele frequency for rs2956540 and rs10519694.

Results

Genotype frequency was significantly different between cases and controls for rs2956540 (P value = 0.019). The rs2956540 C allele carriers were significantly more frequent among KC cases than healthy controls (P value_chi-square = 0.015, P value_{Fisher exact} = 0.017). There was a significant difference in genotype frequency between groups for rs10519694 (P value = 0.001). T allele carriers were significantly more frequent among KC patients (P value_chi-square = 0.002, P value_{Fisher exact} = 0.001). Sex stratification revealed no significant differences in genotype frequency between males and females in cases and controls. Fitting the general linear model showed that rs10519694 could be considered a predictor for the development of KC (P value = 0.001); however, this was not observed for rs2956540 (P value = 0.323).

Conclusions

rs2956540 and rs10519694 are associated with KC in a population of Iranian descent. rs10519694 could potentially be used for KC risk prediction.

The Underlying Relationship between Keratoconus and Down Syndrome

Keratoconus (KC) is one of the most significant corneal disorders worldwide, characterized by the progressive thinning and cone-shaped protrusion of the cornea, which can lead to severe visual impairment. The prevalence of KC varies greatly by ethnic groups and geographic regions and has been observed to be higher in recent years. Although studies reveal a possible link between KC and genetics, hormonal disturbances, environmental factors, and specific comorbidities such as Down Syndrome (DS), the exact cause of KC remains unknown. The incidence of KC ranges from 0% to 71% in DS patients, implying that as the worldwide population of DS patients grows, the number of KC patients may continue to rise significantly. As a result, this review aims to shed more light on the underlying relationship between KC and DS by examining the genetics relating to the cornea, central corneal thickness (CCT), and mechanical forces on the cornea, such as vigorous eye rubbing. Furthermore, this review discusses KC diagnostic and treatment strategies that may help detect KC in DS patients, as well as the available DS mouse models that could be used in modeling KC in DS patients. In summary, this review will provide improved clinical knowledge of KC in DS patients and promote additional KC-related research in these patients to enhance their eyesight and provide suitable treatment targets.

Family-based exome sequencing identifies candidate genes related to keratoconus in Chinese families

Background: Keratoconus (KC) is a complex corneal disorder with a strong genetic component. The present study aimed to identify candidate genes related to KC in Chinese families.

Methods: Family-based exome sequencing was performed in ten patients suffering from KC who belong to five families with two affected members in each. The candidate rare variants were identified with multi-step bioinformatics analysis. The STRING website was used to perform the protein interaction of the identified genes.

Results: Our analyses identified 32 candidate rare variants in 13 genes by family-based exome sequencing. The molecular analyses of identified genes showed that EPCAM directly interacted with CTNNB1 of the Hippo signaling pathway and focal adhesion pathway, and directly interacted with CTNNB1, CDH1 of the WNT signaling pathway. SHROOM3 directly interacted with ROCK2, ROCK1 of the focal adhesion pathway. SYNE1 directly interacted with MUSK of the extracellular matrix organization pathway. TEK directly interacted with VEGFA, SHC1, PIK3R1, GRB2 of the focal adhesion pathway. TTN directly interacted with CAPN3 of the extracellular matrix organization pathway.

Conclusion: The EPCAM, SHROOM3, SYNE1, TEK, and TTN genes were potential high-risk candidate pathogenic genes of familial KC. The findings might significantly improve our understanding of the genetic etiology of the disease, providing novel insights on KC pathogenesis.

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.

Genome-Wide Association Study Identifies Two Common Loci Associated with Pigment Dispersion Syndrome/Pigmentary Glaucoma and Implicates Myopia in its Development

PURPOSE

To identify genetic variants associated with pigment dispersion syndrome (PDS) and pigmentary glaucoma (PG) in unrelated patients and to further understand the genetic and potentially causal relationships between PDS and associated risk factors.

DESIGN

A 2-stage genome-wide association meta-analysis with replication and subsequent in silico analyses including Mendelian randomization.

PARTICIPANTS

A total of 574 cases with PG or PDS and 52 627 controls of European descent.

METHODS

Genome-wide association analyses were performed in 4 cohorts and meta-analyzed in 3 stages: (1) a discovery meta-analysis was performed in 3 cohorts, (2) replication was performed in the fourth cohort, and (3) all 4 cohorts were meta-analyzed to increase statistical power. Two-sample Mendelian randomization was used to determine whether refractive error and intraocular pressure exert causal effects over PDS.

MAIN OUTCOME MEASURES

The association of genetic variants with PDS and whether myopia exerts causal effects over PDS.

RESULTS

Significant association was present at 2 novel loci for PDS/PG. These loci and follow-up analyses implicate the genes gamma secretase activator protein (GSAP) (lead single nucleotide polymorphism [SNP]: rs9641220, P = 6.0×10-10) and glutamate metabotropic receptor 5 (GRM5)/TYR (lead SNP: rs661177, P = 3.9×10-9) as important factors in disease risk. Mendelian randomization showed significant evidence that negative refractive error (myopia) exerts a direct causal effect over PDS (P = 8.86×10-7).

CONCLUSIONS

Common SNPs relating to the GSAP and GRM5/TYR genes are associated risk factors for the development of PDS and PG. Although myopia is a known risk factor, this study uses genetic data to demonstrate that myopia is, in part, a cause of PDS and PG.