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.

Update on the genetics of keratoconus.

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

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.

PPIP5K2 and PCSK1 are Candidate Genetic Contributors to Familial Keratoconus.

Keratoconus (KC) is the most common corneal ectatic disorder affecting >300,000 people in the US. KC normally has its onset in adolescence, progressively worsening through the third to fourth decades of life. KC patients report significant impaired vision-related quality of life. Genetic factors play an important role in KC pathogenesis. To identify novel genes in familial KC patients, we performed whole exome and genome sequencing in a four-generation family. We identified potential variants in the PPIP5K2 and PCSK1 genes. Using in vitro cellular model and in vivo gene-trap mouse model, we found critical evidence to support the role of PPIP5K2 in normal corneal function and KC pathogenesis. The gene-trap mouse showed irregular corneal surfaces and pathological corneal thinning resembling KC. For the first time, we have integrated corneal tomography and pachymetry mapping into characterization of mouse corneal phenotypes which could be widely implemented in basic and translational research for KC diagnosis and therapy in the future.

Corneal Perforation After Corneal Cross-Linking in Keratoconus Associated With Potentially Pathogenic ZNF469 Mutations.

PURPOSE: To report a case of bilateral and repetitive corneal perforations after corneal cross-linking (CXL) for keratoconus in a woman harboring potentially pathogenic variants in the ZNF469 gene and to characterize the keratoconus phenotype in this woman and her daughter who shared the same ZNF469 mutations. METHODS: Clinical characterization of the proband and her daughter followed by sequencing of the genes associated with brittle cornea syndrome, ZNF469 and PRDM5, in both individuals. RESULTS: An Ashkenazi Jewish woman in her sixth decade presented with diffuse corneal thinning and progressive steepening consistent with keratoconus. After CXL, epithelium-off in the first eye and epithelium-on in the second, she developed spontaneous corneal perforations in each eye. Her daughter in her fourth decade demonstrated a similar pattern of diffuse corneal thinning and progressive corneal steepening but did not undergo CXL and did not develop corneal perforation. Screening of the ZNF469 and PRDM5 genes revealed 3 missense ZNF469 variants (c.2035G>A, c.10244G>C, and c.11119A>G) in cis arrangement on 1 allele of ZNF469 in both proband and her daughter. Although the 3 variants share low (<0.01) global minor allele frequencies, each has significantly higher minor allele frequencies (0.01-0.03) in the Ashkenazi Jewish population, leading to uncertainty regarding a pathogenic role for the identified variants. CONCLUSIONS: CXL may be associated with the development of corneal perforation in particular at-risk individuals with keratoconus. Identifying clinical and genetic risk factors, including screening of ZNF469 and PRDM5, may be useful in the prevention of significant complications after CXL.

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.

Differential Expression of Coding and Long Noncoding RNAs in Keratoconus-Affected Corneas.

Purpose: Keratoconus (KC) is the most common corneal ectasia. We aimed to determine the differential expression of coding and long noncoding RNAs (lncRNAs) in human corneas affected with KC. Methods: From the corneas of 10 KC patients and 8 non-KC healthy controls, 200 ng total RNA was used to prepare sequencing libraries with the SMARTer Stranded RNA-Seq kit after ribosomal RNA depletion, followed by paired-end 50-bp sequencing with Illumina Sequencer. Differential analysis was done using TopHat/Cufflinks with a gene file from Ensembl and a lncRNA file from NONCODE. Pathway analysis was performed using WebGestalt. Using the expression level of differentially expressed coding and noncoding RNAs in each sample, we correlated their expression levels in KC and controls separately and identified significantly different correlations in KC against controls followed by visualization using Cytoscape. Results: Using |fold change| ≥ 2 and a false discovery rate ≤ 0.05, we identified 436 coding RNAs and 584 lncRNAs with differential expression in the KC-affected corneas. Pathway analysis indicated the enrichment of genes involved in extracellular matrix, protein binding, glycosaminoglycan binding, and cell migration. Our correlation analysis identified 296 pairs of significant KC-specific correlations containing 117 coding genes enriched in functions related to cell migration/motility, extracellular space, cytokine response, and cell adhesion. Our study highlighted the potential roles of several genes (CTGF, SFRP1, AQP5, lnc-WNT4-2:1, and lnc-ALDH3A2-2:1) and pathways (TGF-ß, WNT signaling, and PI3K/AKT pathways) in KC pathogenesis. Conclusions: Our RNA-Seq-based differential expression and correlation analyses have identified many potential KC contributing coding and noncoding RNAs.

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.

TSC1 Mutations in Keratoconus Patients With or Without Tuberous Sclerosis.

Purpose: To test candidate genes TSC1 and TSC2 in a family affected by tuberous sclerosis complex (TSC) where proband was also diagnosed with bilateral keratoconus (KC) and to test the hypothesis that defects in the same gene may lead to a nonsyndromic KC. Methods: Next-generation sequencing of TSC1 and TSC2 genes was performed in a proband affected by TSC and KC. Identified mutation was confirmed by Sanger DNA sequencing. Whole exome sequencing (WES) was performed in patients with nonsyndromic KC. Sanger DNA sequencing was used to confirm WES results and to screen additional patients. RT-PCR was used to investigate TSC1 expression in seven normal human corneas and eight corneas from patients with KC. Various in silico tools were employed to model functional consequences of identified mutations. Results: A heterozygous nonsense TSC1 mutation g.132902703C>T (c.2293C>T, p.Gln765Ter) was identified in a patient with TSC and KC. Two heterozygous missense TSC1 variants g.132896322A>T (c.3408A>T, p.Asp1136Glu) and g.132896452G>A (c.3278G>A, p.Arg1093Gln) were identified in three patients with nonsyndromic KC. Two mutations were not present in The Genome Aggregation (GnomAD), The Exome Aggregation (ExAC), and 1000 Genomes (1000G) databases, while the third one was present in GnomAD and 1000G with minor allele frequencies (MAF) of 0.00001 and 0.0002, respectively. We found TSC1 expressed in normal corneas and KC corneas, albeit with various levels. Conclusions: Here for the first time we found TSC1 gene to be involved in bilateral KC and TSC as well as with nonsyndromic KC, supporting the hypothesis that diverse germline mutations of the same gene can cause genetic disorders with overlapping clinical features.

Genetics in Keratoconus: where are we?

Keratoconus (KC) is a non-inflammatory thinning and protrusion of the cornea in which the cornea assumes a conical shape. Complex etiology of this condition at present remains an enigma. Although environmental factors have been involved in KC pathogenesis, strong underlining genetic susceptibility has been proven. The lack of consistent findings among early genetic studies suggested a heterogeneity and complex nature of the genetic contribution to the development of KC. Recently, genome-wide linkage studies (GWLS) and genome-wide association studies (GWAS) were undertaken. Next-generation sequencing (NGS)-based genomic screens are also currently being carried out. Application of these recently developed comprehensive genetic tools led to a much greater success and increased reproducibility of genetic findings in KC. Involvement of the LOX gene identified through GWLS has been confirmed in multiple cohorts of KC patients around the world. KC susceptibility region located at the 2q21.3 chromosomal region near the RAB3GAP1 gene identified through GWAS was independently replicated. Rare variants in the ZNF469 gene (mutated in corneal dystrophy Brittle Cornea Syndrome) and in the TGFBI gene (mutated in multiple corneal epithelial-stromal TGFBI dystrophies) have been repeatedly identified in familial and sporadic KC patients of different ethnicities. Additional comprehensive strategies using quantitative endophenotypes have been successfully employed to bring further understanding to the genetics of KC. Additional genetic determinants including the COL5A1 gene have been identified in the GWAS of KC-related trait central corneal thickness. These recent discoveries confirmed the importance of the endophenotype approach for studying complex genetic diseases such as KC and showed that different connective tissue disorders may have the same genetic determinants.

Pseudouridine synthase 1 deficient mice, a model for Mitochondrial Myopathy with Sideroblastic Anemia, exhibit muscle morphology and physiology alterations.

Mitochondrial myopathy with lactic acidosis and sideroblastic anemia (MLASA) is an oxidative phosphorylation disorder, with primary clinical manifestations of myopathic exercise intolerance and a macrocytic sideroblastic anemia. One cause of MLASA is recessive mutations in PUS1, which encodes pseudouridine (Ψ) synthase 1 (Pus1p). Here we describe a mouse model of MLASA due to mutations in PUS1. As expected, certain Ψ modifications were missing in cytoplasmic and mitochondrial tRNAs from Pus1(-/-) animals. Pus1(-/-) mice were born at the expected Mendelian frequency and were non-dysmorphic. At 14 weeks the mutants displayed reduced exercise capacity. Examination of tibialis anterior (TA) muscle morphology and histochemistry demonstrated an increase in the cross sectional area and proportion of myosin heavy chain (MHC) IIB and low succinate dehydrogenase (SDH) expressing myofibers, without a change in the size of MHC IIA positive or high SDH myofibers. Cytochrome c oxidase activity was significantly reduced in extracts from red gastrocnemius muscle from Pus1(-/-) mice. Transmission electron microscopy on red gastrocnemius muscle demonstrated that Pus1(-/-) mice also had lower intermyofibrillar mitochondrial density and smaller mitochondria. Collectively, these results suggest that alterations in muscle metabolism related to mitochondrial content and oxidative capacity may account for the reduced exercise capacity in Pus1(-/-) mice.

Abnormal regulation of extracellular matrix and adhesion molecules in corneas of patients with keratoconus.

AIM: To identify changes in the expression of genes coding for extracellular matrix (ECM) proteins in patients with non-inflammatory corneal disorder keratoconus (KC), patients with corneal scarring, and normal controls. MATERIALS AND METHODS: Total RNA extracted from corneal tissue of 13 KC patients, 2 patients with corneal scaring and 4 normal controls was analyzed using Human Extracellular Matrix & Adhesion Molecules Profiler PCR Array. Statistically significant changes in gene expression were identified using the Data Analysis software. RESULTS: Comparison of KC and control corneas with thresholds of 1.5 or greater fold change and a p-value of 0.05 or lower, revealed 21 differentially expressed genes, 16 genes were downregulated and 5 were upregulated. Among transcripts downregulated in KC patients we identified THBS1, ADAMTS1, SPP1, several collagens and integrins. We found TGFBI (BIGH3) gene was the most significantly upregulated transcript. CONCLUSION: Development of keratoconus results in deregulation of gene expression of extracellular matrix and adhesion molecules. CLINICAL SIGNIFICANCE: Downregulation of collagens and upregulation of TGFBI repeatedly identified in KC patients may be used as clinical markers of the disease.

C.57 C > T Mutation in MIR 184 is Responsible for Congenital Cataracts and Corneal Abnormalities in a Five-generation Family from Galicia, Spain.

A c.57 C > T mutation in the seed region of MIR184 located at the 15q25.1 chromosomal region has been independently associated with autosomal dominant keratoconus with early-onset anterior polar cataract in the Northern Irish family and with autosomal dominant EDICT (Endothelial Dystrophy, Iris hypoplasia, Congenital cataracts, and stromal Thinning) syndrome. In this study we report a five-generation family originating in Galicia, Spain with early onset cataracts and variable corneal abnormalities which include non-ectatic corneal thinning and severe early-onset keratoconus. We identified a heterozygous c.57 C > T mutation in miR-184 in the proband and two additional affected relatives on the maternal side. This finding represents a third independent occurrence of this mutation in familiar ocular disease thus strengthening the link between miR-184 abnormalities and inherited eye defects.

Optical coherence tomography combined with videokeratography to differentiate mild keratoconus subtypes.

PURPOSE: To develop parameters using a combination of optical coherence tomography (OCT) and videokeratography to detect early keratoconus. METHODS: Videokeratography, wavefront analysis, and measured OCT indices were performed on 180 normal eyes, 46 eyes with moderate keratoconus, 54 eyes with early keratoconus, 7 eyes with forme fruste keratoconus, and 16 eyes with keratoconus "suspect" to determine the most sensitive parameters for separating these groups. RESULTS: A combination of videokeratography and OCT indices (inferior-superior [I-S] value and minimum pachymetry) was statistically the most significant in separating the keratoconus groups from normal eyes (P < .001). Using a newly derived index, the minimum pachymetry divided by the I-S value (pachymetry/asymmetry [PA]/I-S index) with a cut-off of 100, 100% of early and forme fruste keratoconus could be identified as being abnormal with 7 normals misclassified (misclassification rate 2.7%). By adding keratoconus "suspect" to the analysis and an I-S value of 1.2 as a cut-off point, 5 "suspects" were classified as normal and 11 normals as abnormal (misclassification rate 7.8%). The PA/I-S index, with a cut-off point of 100, reduced this misclassification rate to 4.4%. CONCLUSIONS: These results suggest that OCT combined with videokeratography may be more useful for differentiating mild forms of keratoconus than videokeratography alone.

Genetic association of COL5A1 variants in keratoconus patients suggests a complex connection between corneal thinning and keratoconus.

PURPOSE: Single nucleotide polymorphisms (SNPs) located near or within the COL5A1 gene, at 9q34.2-q34.3 chromosomal region have been reported in association with central corneal thickness (CCT). Using family linkage analysis, we identified a keratoconus susceptibility locus at 9q34. These findings led us to perform an association study between COL5A1 variation and keratoconus susceptibility. METHODS: A Caucasian case-control cohort of 222 keratoconus patients and 3324 controls was selected as the discovery panel. An independent case-control panel of 304 cases and 518 controls and a family panel of 186 subjects were replicated for genotyping and association. Forty-four SNPs (21 for discovery and 23 for fine-mapping) spanning 300 kilobases in and around COL5A1 were genotyped and tested for genetic association. Logistic regression models implemented in PLINK were used to test for association in case controls. Generalized estimating equation models accounting for familial correlations implemented in genome-wide interaction analyses with family data were used for association testing in families. RESULTS: Two CCT associated SNPs (rs1536482 and rs7044529 near and within COL5A1) were identified in the keratoconus discovery cohort (P values of 6.5 × 10(-3) and 7.4 × 10(-3)). SNP rs1536482 was replicated in the second case-control sample (P = 0.02), and SNP rs7044529 was replicated in a keratoconus family panel (P = 0.03). Meta P values of rs1536482 and rs7044529 in the keratoconus cohorts were 1.5 × 10(-4) (odds ratio [OR] = 1.30) and 2.9 × 10(-3) (OR = 1.39). After Bonferroni correction, the association of SNP rs1536482 remained significant (P = 6.5 × 10(-3)). CONCLUSIONS: SNPs in the COL5A1 region, which regulate normal variation in CCT, may play a role in the thinning associated with keratoconus.

An association between the calpastatin (CAST) gene and keratoconus.

PURPOSE: Keratoconus (KC) is a genetically heterogeneous corneal dystrophy. Previously, we performed 2 genome-wide linkage scans in a 4-generation autosomal dominant pedigree and repeatedly mapped a KC locus to a genomic region located on chromosome 5q overlapping the gene encoding the inhibitor of calpains, calpastatin (CAST). To test whether variants in CAST gene are involved in genetic susceptibility to KC, we performed genetic testing of polymorphic markers in CAST gene in family and case-control panels of patients with KC. METHODS: We genotyped single-nucleotide polymorphisms (SNPs) located in CAST gene in 262 patients in 40 white KC families and in a white case-control panel with 304 cases and 518 controls. Generalized estimating equation models accounting for familial correlations implemented in GWAF program were used for association testing in families. Logistic regression models implemented in PLINK were performed to test the associations in case-control samples. RESULTS: Genetic testing of the first set of 7 SNPs in familial samples revealed 2 tentative nominally significant markers (rs4869307, P = 0.03; rs27654, P = 0.07). Additional genotyping of 12 tightly spaced SNPs identified CAST SNP rs4434401 to be associated with KC in both familial and case-control panels with P values of 0.005 and 0.05, respectively, and with combined meta P value of familial and case-control cohorts of 0.002 or after Bonferroni correction of 0.04. CONCLUSIONS: Linkage analysis and genetic association support involvement of CAST gene in the genetic susceptibility to KC. In silico analysis of CAST expression suggests differential regulation of calpain/calpastatin system in cornea as a potential mechanism of functional defect.

Genome-wide association analyses identify multiple loci associated with central corneal thickness and keratoconus.

Central corneal thickness (CCT) is associated with eye conditions including keratoconus and glaucoma. We performed a meta-analysis on >20,000 individuals in European and Asian populations that identified 16 new loci associated with CCT at genome-wide significance (P < 5 × 10(-8)). We further showed that 2 CCT-associated loci, FOXO1 and FNDC3B, conferred relatively large risks for keratoconus in 2 cohorts with 874 cases and 6,085 controls (rs2721051 near FOXO1 had odds ratio (OR) = 1.62, 95% confidence interval (CI) = 1.4-1.88, P = 2.7 × 10(-10), and rs4894535 in FNDC3B had OR = 1.47, 95% CI = 1.29-1.68, P = 4.9 × 10(-9)). FNDC3B was also associated with primary open-angle glaucoma (P = 5.6 × 10(-4); tested in 3 cohorts with 2,979 cases and 7,399 controls). Further analyses implicate the collagen and extracellular matrix pathways in the regulation of CCT.

Variation in the lysyl oxidase (LOX) gene is associated with keratoconus in family-based and case-control studies.

PURPOSE: Keratoconus is a bilateral noninflammatory progressive corneal disorder with complex genetic inheritance and a common cause for cornea transplantation in young adults. A genomewide linkage scan in keratoconus families identified a locus at 5q23.2, overlapping the gene coding for the lysyl oxidase (LOX). LOX encodes an enzyme responsible for collagen cross-linking in a variety of tissues including the cornea. Corneal collagen cross-linking with long-wave ultraviolet light and riboflavin is a promising new treatment for keratoconus. To determine whether LOX is a genetic determinant of the pathogenesis of keratoconus, we analyzed association results of LOX polymorphisms in two independent case-control samples and in keratoconus families. METHODS: Association results were analyzed of single-nucleotide polymorphisms (SNPs) in the LOX gene from a Genome-Wide Association Study (GWAS) investigation in two independent panels of patients with keratoconus and controls and in keratoconus families. RESULTS: Evidence of association was found at SNPs rs10519694 and rs2956540 located in intron 4 of LOX in the GWAS discovery case-control panel with P values of 2.3×10(-3) and 7×10(-3), respectively. The same two SNPs were found to be associated with keratoconus by family-based association testing with P values of 2.7×10(-3) and 7.7×10(-4), respectively. Meta P values of 4.0×10(-5) and 4.0×10(-7) were calculated for SNPs rs10519694 and rs2956540 by analyzing case-control and family samples simultaneously. Sequencing of LOX exons in a subset of keratoconus patients identified two polymorphisms, rs1800449 and rs2288393, located in LOX transcripts I and II, associated with keratoconus in case-control and family samples with a meta P value of 0.02. CONCLUSIONS: Results provided strong genetic evidence that LOX variants lead to increased susceptibility to developing of keratoconus.