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.

ImportanceKeratoconus is an important cause of visual loss in young adults, but little is known about its genetic causes. Understanding the genetic determinants of corneal biomechanical factors may in turn teach us about keratoconus etiology.ObjectivesTo identify genetic associations with corneal biomechanical properties and to examine whether these genetic variants are associated with keratoconus.Design, Setting, and ParticipantsA stage 1 discovery and replication genome-wide association study (GWAS) of corneal biomechanical properties was performed in 2 cross-sectional populations (6645 participants from the European Prospective Investigation into Cancer and Nutrition [EPIC]-Norfolk Eye Study and 2384 participants from the TwinsUK study). In stage 2, the association of genetic determinants identified in stage 1 with keratoconus was examined in a case-control study. A total of 752 patients with keratoconus were compared with 974 TwinsUK participants (undergoing direct sequencing) or 13 828 EPIC-Norfolk participants (undergoing genotyping and imputation) who were not part of the stage 1 analysis. Data were collected from March 1, 1993, through March 13, 2017, and analyzed from November 1, 2015, through February 1, 2018.ExposuresIn stage 1, allele dosage at genome-wide single-nucleotide polymorphisms (SNPs); in stage 2, allele dosage at SNPs with genome-wide significance (P < 5 × 10-8) in stage 1 and not previously reported as associated with corneal disease.Main Outcomes and MeasuresIn stage 1, corneal hysteresis (CH) and corneal resistance factor (CRF), measured with the Ocular Response Analyzer (ORA); in stage 2, association with keratoconus compared with controls.ResultsAmong 6645 participants in the discovery cohort (3635 women (54.7%); mean age, 69 years [range, 48-92 years]), 7 genome-wide significant loci associated with CH or CRF were identified that were independently replicated. Two further suggestive loci were identified after meta-analysis. To date, 5 of the identified loci, at ANAPC1, ADAMTS8, ADAMTS17, ABCA6, and COL6A1, have not previously been reported as associated with corneal disease. The ABCA6 locus (rs77542162) was associated with keratoconus using the TwinsUK (odds ratio [OR], 0.50; 95% CI, 0.27-0.92; P = .03) and EPIC-Norfolk controls (OR, 0.39; 95% CI, 0.22-0.70; P = .002). The other loci were associated with keratoconus using TwinsUK (OR per effect allele for ADAMTS8, 0.51 [95% CI, 0.37-0.71; P = 7.9 × 10-5]; for COL6A1, 1.65 [95% CI, 1.05-2.59; P = .03]) or EPIC-Norfolk (OR per effect allele for ANAPC1, 0.78 [95% CI, 0.68-0.89; P = 3.7 × 10-4]; for ADAMTS17, 0.82 [95% CI, 0.68-0.99; P = .04]) controls.Conclusions and RelevanceFive loci that are associated with corneal biomechanical properties and that have suggestive associations with keratoconus were reported. These findings suggest the role of type VI collagen, extracellular matrix, and connective-tissue development for corneal biomechanics and keratoconus and the role of CH and CRF as biomarkers for keratoconus.

The aim of this study was to identify the molecular genetic cause of disease in posterior polymorphous corneal dystrophy (PPCD) probands of diverse origin and to assess the utility of massively parallel sequencing in the detection of ZEB1 mutations. We investigated a total of 12 families (five British, four Czech, one Slovak and two Swiss). Ten novel and two recurrent disease-causing mutations in ZEB1, were identified in probands by Sanger (n = 5), exome (n = 4) and genome (n = 3) sequencing. Sanger sequencing was used to confirm the mutations detected by massively parallel sequencing, and to perform segregation analysis. Genome sequencing revealed that one proband harboured a novel similar to 0.34 Mb heterozygous de novo deletion spanning exons 1-7 and part of exon 8. Transcript analysis confirmed that the ZEB1 transcript is detectable in blood-derived RNA samples and that the disease-associated variant c.482-2A > G leads to aberrant pre-mRNA splicing. De novo mutations, which are a feature of PPCD3, were found in the current study with an incidence rate of at least 16.6%. In general, massively parallel sequencing is a time-efficient way to detect PPCD3-associated mutations and, importantly, genome sequencing enables the identification of full or partial heterozygous ZEB1 deletions that can evade detection by both Sanger and exome sequencing. These findings contribute to our understanding of PPCD3, for which currently, 49 pathogenic variants have been identified, all of which are predicted to be null alleles.

PurposeThe majority of anterior corneal dystrophies are caused by dominant mutations in TGFBI (transforming growth factor β-induced) collectively known as the epithelial-stromal TGFBI dystrophies. Most cases of epithelial basement membrane dystrophy (EBMD) are thought to result from a degenerative (nongenetic) process; however, a minority of cases are associated with specific TGFBI mutations. We evaluated the spectrum of TGFBI mutations and associated phenotypes in a United Kingdom cohort with typical epithelial-stromal TGFBI dystrophies and an EBMD cohort.MethodsWe recruited 68 probands with a clinical diagnosis of epithelial-stromal TGFBI dystrophy and 23 probands with bilateral EBMD. DNA was extracted from peripheral leukocytes, and TGFBI was bi-directly Sanger sequenced.ResultsNine TGFBI mutations were identified. The most common occurred at the mutation hot-spot residues R124 and R555 in 61 probands; these individuals had a genotype-phenotype correlation consistent with prior reports. Four probands with lattice corneal dystrophy carried a mutation in exon 14: p.(A620D), p.(V625D), and p.(H626R). We identified a p.(G623D) mutation in five probands, including two probands from the EBMD cohort. These subjects typically had an onset of severe recurrent corneal epithelial erosion in the fourth decade with mild diffuse or geographic subepithelial corneal opacities and only small anterior stromal lattice structures in older individuals. Symptoms of painful epithelial erosion improved markedly following phototherapeutic keratectomy.ConclusionsThere was a strong correlation between genotype and phenotype for the majority of TGFBI mutations. In this cohort, the p.(G623D) mutation caused a greater proportion of TGFBI-associated disease than anticipated, associated with variable phenotypes including individuals diagnosed with EBMD.