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.

Genetic and Environmental Risk Factors for Keratoconus.

Keratoconus, a progressive corneal ectasia, is a complex disease with both genetic and environmental risk factors. The exact etiology is not known and is likely variable between individuals. Conditions such as hay fever and allergy are associated with increased risk, while diabetes may be protective. Behaviors such as eye rubbing are also implicated, but direct causality has not been proven. Genetics plays a major role in risk for some individuals, with many large pedigrees showing autosomal inheritance patterns. Several genes have been implicated using linkage and follow-up sequencing in these families. Genome-wide association studies for keratoconus and for quantitative traits such as central corneal thickness have identified several genetic loci that contribute to a cumulative risk for keratoconus, even in people without a family history of the disease. Identification of risk genes for keratoconus is improving our understanding of the biology of this complex disease.

Association of Genetic Variation With Keratoconus.

Importance: Keratoconus is a condition in which the cornea progressively thins and protrudes in a conical shape, severely affecting refraction and vision. It is a major indication for corneal transplant. To discover new genetic loci associated with keratoconus and better understand the causative mechanism of this disease, we performed a genome-wide association study on patients with keratoconus. Objective: To identify genetic susceptibility regions for keratoconus in the human genome. Design, Setting, and Participants: This study was conducted with data from eye clinics in Australia, the United States, and Northern Ireland. The discovery cohort of individuals with keratoconus and control participants from Australia was genotyped using the Illumina HumanCoreExome single-nucleotide polymorphism array. After quality control and data cleaning, genotypes were imputed against the 1000 Genomes Project reference panel (phase III; version 5), and association analyses were completed using PLINK version 1.90. Single-nucleotide polymorphisms with P < 1.00 × 10-6 were assessed for replication in 3 additional cohorts. Control participants were drawn from the cohorts of the Blue Mountains Eye Study and a previous study of glaucoma. Replication cohorts were from a previous keratoconus genome-wide association study data set from the United States, a cohort of affected and control participants from Australia and Northern Ireland, and a case-control cohort from Victoria, Australia. Data were collected from January 2006 to March 2019. Main Outcomes and Measures: Associations between keratoconus and 6 252 612 genetic variants were estimated using logistic regression after adjusting for ancestry using the first 3 principal components. Results: The discovery cohort included 522 affected individuals and 655 control participants, while the replication cohorts included 818 affected individuals (222 from the United States, 331 from Australia and Northern Ireland, and 265 from Victoria, Australia) and 3858 control participants (2927 from the United States, 229 from Australia and Northern Ireland, and 702 from Victoria, Australia). Two novel loci reached genome-wide significance (defined as P < 5.00 × 10-8), with a P value of 7.46 × 10-9 at rs61876744 in patatin-like phospholipase domain-containing 2 gene (PNPLA2) on chromosome 11 and a P value of 6.35 × 10-12 at rs138380, 2.2 kb upstream of casein kinase I isoform epsilon gene (CSNK1E) on chromosome 22. One additional locus was identified with a P value less than 1.00 × 10-6 in mastermind-like transcriptional coactivator 2 (MAML2) on chromosome 11 (P = 3.91 × 10-7). The novel locus in PNPLA2 reached genome-wide significance in an analysis of all 4 cohorts (P = 2.45 × 10-8). Conclusions and Relevance: In this relatively large keratoconus genome-wide association study, we identified a genome-wide significant locus for keratoconus in the region of PNPLA2 on chromosome 11.

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

Emmanuelle Souzeau, who contributed to analysis of data, was inadvertently omitted from the author list in the originally published version of this Article. This has now been corrected in both the PDF and HTML versions of the Article.

Rare, potentially pathogenic variants in 21 keratoconus candidate genes are not enriched in cases in a large Australian cohort of European descent.

Many genes have been suggested as candidate genes for keratoconus based on their function, their proximity to associated polymorphisms or due to the identification of putative causative variants within the gene. However, very few of these genes have been assessed for rare variation in keratoconus more broadly. In contrast, VSX1 and SOD1 have been widely assessed, however, the vast majority of studies have been small and the findings conflicting. In a cohort of Australians of European descent, consisting of 385 keratoconus cases and 396 controls, we screened 21 keratoconus candidate genes: BANP, CAST, COL4A3, COL4A4, COL5A1, FOXO1, FNDC3B, HGF, IL1A, IL1B, ILRN, IMMP2L, MPDZ, NFIB, RAB3GAP1, RAD51, RXRA, SLC4A11, SOD1, TF and VSX1. The candidate genes were sequenced in these individuals by either whole exome sequencing or targeted gene sequencing. Variants were filtered to identify rare (minor allele frequency <1%), potentially pathogenic variants. A total of 164 such variants were identified across the two groups with no variants fulfilling these criteria in cases in IL1RN, BANP, IL1B, RAD51 or SOD1. The frequency of variants was compared between cases and controls using chi-square or Fishers' Exact tests for each gene with at least one rare potentially pathogenic variant identified in the case cohort. The number of rare potentially pathogenic variants per gene ranged from three (RXRA) to 102 (MPDZ), however for all genes, there was no difference in the frequency between the cases and controls. We conclude that rare potentially pathogenic variation in the 21 candidate genes assessed do not play a major role in keratoconus susceptibility and pathogenesis.

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

Central corneal thickness (CCT) is a highly heritable trait associated with complex eye diseases such as keratoconus and glaucoma. We perform a genome-wide association meta-analysis of CCT and identify 19 novel regions. In addition to adding support for known connective tissue-related pathways, pathway analyses uncover previously unreported gene sets. Remarkably, >20% of the CCT-loci are near or within Mendelian disorder genes. These included FBN1, ADAMTS2 and TGFB2 which associate with connective tissue disorders (Marfan, Ehlers-Danlos and Loeys-Dietz syndromes), and the LUM-DCN-KERA gene complex involved in myopia, corneal dystrophies and cornea plana. Using index CCT-increasing variants, we find a significant inverse correlation in effect sizes between CCT and keratoconus (r = -0.62, P = 5.30 × 10-5) but not between CCT and primary open-angle glaucoma (r = -0.17, P = 0.2). Our findings provide evidence for shared genetic influences between CCT and keratoconus, and implicate candidate genes acting in collagen and extracellular matrix regulation.

Identification of novel mutations causing pediatric cataract in Bhutan, Cambodia, and Sri Lanka.

BACKGROUND: Pediatric cataract is an important cause of blindness and visual impairment in children. A large proportion of pediatric cataracts are inherited, and many genes have been described for this heterogeneous Mendelian disease. Surveys of schools for the blind in Bhutan, Cambodia, and Sri Lanka have identified many children with this condition and we aimed to identify the genetic causes of inherited cataract in these populations. METHODS: We screened, in parallel, 51 causative genes for inherited cataracts in 33 probands by Ampliseq enrichment and sequencing on an Ion Torrent PGM. Rare novel protein coding variants were assessed for segregation in family members, where possible, by Sanger sequencing. RESULTS: We identified 24 rare (frequency <1% in public databases) or novel protein coding variants in 12 probands and confirmed segregation of variants with disease in the extended family where possible. Of these, six are predicted to be the cause of disease in the patient, with four other variants also highly likely to be pathogenic. CONCLUSION: This study found that 20%-30% of patients in these countries have a mutation in a known cataract causing gene, which is considerably lower than the 60%-70% reported in Caucasian cohorts. This suggests that additional cataract genes remain to be discovered in this cohort of Asian pediatric cataract patients.

Rare, Potentially Pathogenic Variants in ZNF469 Are Not Enriched in Keratoconus in a Large Australian Cohort of European Descent.

Purpose: The Zinc Finger Protein 469 (ZNF469) gene has been proposed as a candidate gene for keratoconus due to the association of an upstream polymorphism (rs9938149) with the disease in two independent studies, and the role of the gene in the autosomal recessive disease Brittle Cornea Syndrome. Coding variants in ZNF469 have been assessed for association with keratoconus in several small studies, with conflicting results. We assessed rare, potentially pathogenic variants in ZNF469 for enrichment in keratoconus patients in a cohort larger than all previous studies combined. Methods: ZNF469 was sequenced in 385 Australian keratoconus patients of European descent, 346 population controls, and 230 ethnically matched screened controls by either whole exome sequencing or targeted gene sequencing. The frequency of rare and very rare potentially pathogenic variants was compared between cases and controls using χ2 or Fisher's exact tests and further explored using a gene based test (Sequence Kernel Association Test [SKAT]), weighting on the rarity of variants. Results: A total of 49 rare, including 33 very rare, potentially pathogenic variants were identified across all groups. No enrichment of rare or very rare potentially pathogenic variants in ZNF469 was observed in our cases compared to the control groups following analysis using χ2 or Fisher's exact tests. This finding was further supported by the SKAT results, which found no significant difference in the frequency of variants predicted to be damaging between cases and either control group (P = 0.06). Conclusions: Rare variants in ZNF469 do not contribute to keratoconus susceptibility and do not account for the association at rs9938149.

WNT10A exonic variant increases the risk of keratoconus by decreasing corneal thickness.

Keratoconus is a degenerative eye condition which results from thinning of the cornea and causes vision distortion. Treatments such as ultraviolet (UV) cross-linking have proved effective for management of keratoconus when performed in early stages of the disease. The central corneal thickness (CCT) is a highly heritable endophenotype of keratoconus, and it is estimated that up to 95% of its phenotypic variance is due to genetics. Genome-wide association efforts of CCT have identified common variants (i.e. minor allele frequency (MAF) >5%). However, these studies typically ignore the large set of exonic variants whose MAF is usually low. In this study, we performed a CCT exome-wide association analysis in a sample of 1029 individuals from a population-based study in Western Australia. We identified a genome-wide significant exonic variant rs121908120 (P = 6.63 × 10(-10)) in WNT10A. This gene is 437 kb from a gene previously associated with CCT (USP37). We showed in a conditional analysis that the WNT10A variant completely accounts for the signal previously seen at USP37. We replicated our finding in independent samples from the Brisbane Adolescent Twin Study, Twin Eye Study in Tasmania and the Rotterdam Study. Further, we genotyped rs121908120 in 621 keratoconus cases and compared the frequency to a sample of 1680 unscreened controls from the Queensland Twin Registry. We found that rs121908120 increases the risk of keratoconus two times (odds ratio 2.03, P = 5.41 × 10(-5)).