A genome-wide association study identifies a potential novel gene locus for keratoconus, one of the commonest causes for corneal transplantation in developed countries

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.

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.

Sequence variants contributing to dysregulated inflammatory responses across keratoconic cone surface in adolescent patients with keratoconus

Background

Keratoconus (KTCN) is the most common corneal ectasia resulting in a conical shape of the cornea. Here, genomic variation in the corneal epithelium (CE) across the keratoconic cone surface in patients with KTCN and its relevance in the functioning of the immune system were assessed.

Methods

Samples from four unrelated adolescent patients with KTCN and two control individuals were obtained during the CXL and PRK procedures, respectively. Three topographic regions, central, middle, and peripheral, were separated towards the whole-genome sequencing (WGS) study embracing a total of 18 experimental samples. The coding and non-coding sequence variation, including structural variation, was assessed and then evaluated together with the previously reported transcriptomic outcomes for the same CE samples and full-thickness corneas.

Results

First, pathway enrichment analysis of genes with identified coding variants pointed to "Antigen presentation" and "Interferon alpha/beta signaling" as the most overrepresented pathways, indicating the involvement of inflammatory responses in KTCN. Both coding and non-coding sequence variants were found in genes (or in their close proximity) linked to the previously revealed KTCN-specific cellular components, namely, "Actin cytoskeleton", "Extracellular matrix", "Collagen-containing extracellular matrix", "Focal adhesion", "Hippo signaling pathway", and "Wnt signaling" pathways. No genomic heterogeneity across the corneal surface was found comparing the assessed topographic regions. Thirty-five chromosomal regions enriched in both coding and non-coding KTCN-specific sequence variants were revealed, with a most representative 5q locus previously recognized as involved in KTCN.

Conclusion

The identified genomic features indicate the involvement of innate and adaptive immune system responses in KTCN pathogenesis.

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

Importance

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

Objective

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

Design, Setting, and Participants

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

Main Outcomes and Measures

Association of CRF and CCT variants with keratoconus risk.

Results

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

Conclusions and Relevance

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

REL and BHLHE40 Variants Are Associated with IL-12 and IL-10 Responses and Tuberculosis Risk

The major human genes regulating Mycobacterium tuberculosis-induced immune responses and tuberculosis (TB) susceptibility are poorly understood. Although IL-12 and IL-10 are critical for TB pathogenesis, the genetic factors that regulate their expression in humans are unknown. CNBP, REL, and BHLHE40 are master regulators of IL-12 and IL-10 signaling. We hypothesized that common variants in CNBP, REL, and BHLHE40 were associated with IL-12 and IL-10 production from dendritic cells, and that these variants also influence adaptive immune responses to bacillus Calmette-Guérin (BCG) vaccination and TB susceptibility. We characterized the association between common variants in CNBP, REL, and BHLHE40, innate immune responses in dendritic cells and monocyte-derived macrophages, BCG-specific T cell responses, and susceptibility to pediatric and adult TB in human populations. BHLHE40 single-nucleotide polymorphism (SNP) rs4496464 was associated with increased BHLHE40 expression in monocyte-derived macrophages and increased IL-10 from peripheral blood dendritic cells and monocyte-derived macrophages after LPS and TB whole-cell lysate stimulation. SNP BHLHE40 rs11130215, in linkage disequilibrium with rs4496464, was associated with increased BCG-specific IL-2+CD4+ T cell responses and decreased risk for pediatric TB in South Africa. SNPs REL rs842634 and rs842618 were associated with increased IL-12 production from dendritic cells, and SNP REL rs842618 was associated with increased risk for TB meningitis. In summary, we found that genetic variations in REL and BHLHE40 are associated with IL-12 and IL-10 cytokine responses and TB clinical outcomes. Common human genetic regulation of well-defined intermediate cellular traits provides insights into mechanisms of TB pathogenesis.

Systematically Displaying the Pathogenesis of Keratoconus via Multi-Level Related Gene Enrichment-Based Review

Keratoconus (KC) is an etiologically heterogeneous corneal ectatic disorder. To systematically display the pathogenesis of keratoconus (KC), this study reviewed all the reported genes involved in KC, and performed an enrichment analysis of genes identified at the genome, transcription, and protein levels respectively. Combined analysis of multi-level results revealed their shared genes, gene ontology (GO), and pathway terms, to explore the possible pathogenesis of KC. After an initial search, 80 candidate genes, 2,933 transcriptional differential genes, and 947 differential proteins were collected. The candidate genes were significantly enriched in extracellular matrix (ECM) related terms, Wnt signaling pathway and cytokine activities. The enriched GO/pathway terms of transcription and protein levels highlight the importance of ECM, cell adhesion, and inflammatory once again. Combined analysis of multi-levels identified 13 genes, 43 GOs, and 12 pathways. The pathogenic relationships among these overlapping factors maybe as follows. The gene mutations/variants caused insufficient protein dosage or abnormal function, together with environmental stimulation, leading to the related functions and pathways changes in the corneal cells. These included response to the glucocorticoid and reactive oxygen species; regulation of various signaling (P13K-AKT, MAPK and NF-kappaB), apoptosis and aging; upregulation of cytokines and collagen-related enzymes; and downregulation of collagen and other ECM-related proteins. These undoubtedly lead to a reduction of extracellular components and induction of cell apoptosis, resulting in the loosening and thinning of corneal tissue structure. This study, in addition to providing information about the genes involved, also provides an integrated insight into the gene-based etiology and pathogenesis of KC.

Anterior pituitary, sex hormones, and keratoconus: Beyond traditional targets

"The Diseases of the Horny-coat of The Eye", known today as keratoconus, is a progressive, multifactorial, non-inflammatory ectatic corneal disorder that is characterized by steepening (bulging) and thinning of the cornea, irregular astigmatism, myopia, and scarring that can cause devastating vision loss. The significant socioeconomic impact of the disease is immeasurable, as patients with keratoconus can have difficulties securing certain jobs or even joining the military. Despite the introduction of corneal crosslinking and improvements in scleral contact lens designs, corneal transplants remain the main surgical intervention for treating keratoconus refractory to medical therapy and visual rehabilitation. To-date, the etiology and pathogenesis of keratoconus remains unclear. Research studies have increased exponentially over the years, highlighting the clinical significance and international interest in this disease. Hormonal imbalances have been linked to keratoconus, both clinically and experimentally, with both sexes affected. However, it is unclear how (molecular/cellular signaling) or when (age/disease stage(s)) those hormones affect the keratoconic cornea. Previous studies have categorized the human cornea as an extragonadal tissue, showing modulation of the gonadotropins, specifically luteinizing hormone (LH) and follicle-stimulating hormone (FSH). Studies herein provide new data (both in vitro and in vivo) to further delineate the role of hormones/gonadotropins in the keratoconus pathobiology, and propose the existence of a new axis named the Hypothalamic-Pituitary-Adrenal-Corneal (HPAC) axis.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSION

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

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

Background

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

Methods

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

Results

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

Conclusions

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

Omics analyses in keratoconus: from transcriptomics to proteomics

Purpose of Review

To summarize the recent advances in transcriptomics and proteomics studies of keratoconus using advanced genome-wide gene and protein expression profiling techniques.

Recent Findings

Second-generation sequencing including RNA sequencing has been widely used to characterize the genome-wide gene expression in corneal tissues or cells affected by keratoconus. Due to different sample types, sequencing platforms, and analysis pipeline, different lists of genes have been identified to be differentially expressed in KC-affected samples. Gene ontology and pathway/network analyses have indicated the involvement of genes related with extracellular matrix, WNT-signaling, TGFβ pathway, and NRF2-regulated network. High throughput proteomics studies using mass spectrometry have uncovered many KC-related protein molecules in pathways related with cytoskeleton, cell matrix, TGFβ signaling, and extracellular matrix remodeling, consistent with gene expression profiling.

Summary

Both transcriptomics and proteomics studies using genome-wide gene/protein expression profiling techniques have identified significant genes/proteins that may contribute to the pathogenesis of keratoconus. These molecules may be involved in functional categories related with extracellular matrix and TGFβ signaling. It is necessary to perform comprehensive gene/protein expression studies using larger sample size, same type of samples, up-to-date platform and bioinformatics tools.

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.

Association of a Novel KIF26B Gene Polymorphism with Susceptibility to Schizophrenia and Breast Cancer: A Case-Control Study

Background

KIF26B gene is found to play essential roles in regulating different aspects of cell proliferation and development of the nervous system. We aimed to determine if rs12407427 T/C polymorphism could affect susceptibility to schizophrenia (SZN) and breast cancer (BC), the two genetically correlated diseases.

Methods

The current case-control study was performed from Aug 2018 to Dec 2018. Briefly, 159 female pathologically confirmed BC cases referring to Alzahra Hospital, Isfahan, Iran, and 102 psychologically confirmed SZN patients (60 males and 42 females) admitted to Baharan Hospital, Zahedan, Iran, were enrolled. Using the salting-out method, genomic DNA was extracted, and variants were genotyped using allele-specific amplification refractory mutation system polymerase chain reaction (ARMS-PCR) method.

Results

The results revealed a significant association between the KIF26B rs12407427 codominant CT (P=0.001), CC (P=0.0001), dominant CT+CC, and recessive CC (P=0.001) genotypes with the risk of developing SZN. Significant correlations were also found regarding rs12407427 and BC susceptibility in different inheritance models, including over-dominant CT (P=0.026), dominant CT+CC (P=0.001), recessive CC (P=0.009), and codominant CT and CC (P=0.001) genotypes. The over-presence of the C allele was also correlated with an increased risk for SZN (P=0.0001) and BC (P=0.0001). Finally, computational analysis predicted that T/C variation in this polymorphism could change the binding sites in proteins involved in splicing.

Conclusion

rs12407427 T/C as a de novo KIF26B variant might be a novel genetic biomarker for SZN and/or BC susceptibility in a sample of the Iranian population.

Cascade of interactions between candidate genes reveals convergent mechanisms in keratoconus disease pathogenesis

Keratoconus is a progressive thinning, steepening and distortion of the cornea which can lead to loss of vision if left untreated. Keratoconus has a complex multifactorial etiology, with genetic and environmental components contributing to the disease pathophysiology. Studies have observed high concordance between monozygotic twins, discordance between dizygotic twins, and high familial segregation indicating the presence of a very strong genetic component in the pathogenesis of keratoconus. The use of genome-wide linkage studies on families and twins, genome-wide association studies (GWAS) on case-controls, next-generation sequencing (NGS)-based genomic screens on both familial and non-familial cohorts have led to the identification of keratoconus candidate genes with much greater success and increased resproducibility of genetic findings. This review focuses on candidate genes identified till date and attempts to understand their role in biological processes underlying keratoconus pathogenesis. In addition, using these genes I propose molecular pathways that could contribute to keratoconus pathogenesis. The pathways identified the presence of direct cross-talk between known candidate genes of keratoconus and remarkably, 28 known candidate genes have a direct relationship among themselves that involves direct protein-protein binding, regulatory activities such as activation and inhibition, chaperone, transcriptional activation/co-activation, and enzyme catalysis. This review attempts to describe these relationships and cross-talks in the context of keratoconus pathogenesis.

Corneal epithelial thickness analysis of forme fruste keratoconus with optical coherence tomography

AIM

To explore the significance of corneal epithelial thickness analysis in diagnosing early keratoconus.

METHODS

There were 26 clinical keratoconus, 21 forme fruste keratoconus, 40 high corneal astigmatism (ΔK) and 40 low ΔK eyes involved in the study. Fourier-domain optical coherence tomography was used to measure the corneal epithelial thickness of four groups. The morphological features of topographic map and the thickness of corneal epithelial thinnest point were analyzed. The distribution curve of corneal epithelial thickness at 45°, 90°, and 135° axial directions that are through the pupil center was also analyzed. One-way ANOVA was performed to compare the data.

RESULTS

The topographic map of forme fruste keratoconus corneal epithelial thickness was uniformity shape; crater shape existed only in clinical keratoconus group; and central island shape mainly existed in high ΔK group. The thinnest point of corneal epithelial thickness of forme fruste keratoconus group was significantly lower than that of low ΔK group (P=0.022). The thickness of corneal epithelium in the forme fruste keratoconus at 90° was thinner than that in the low astigmatism group at -1, and -2 mm points (P _{-1 mm}=0.015, P _{-2 mm}=0.036).

CONCLUSION

The analysis of the thinnest point in forme fruste keratoconus corneal epithelium appears earlier than corneal epithelial remodeling. The topographic map of corneal epithelium in high ΔK eyes appears in central island shape, and can be used for the differential diagnosis of early keratoconus.

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.

The Enigma of Environmental Factors in Keratoconus

Keratoconus is an ectatic disorder with highly complex and varied causes including genetic variations and environmental factors. Its prevalence varies widely between regions and countries. Many environmental factors have been proposed to be associated with keratoconus, but the interpretation of their individual contributions is difficult due to the presence of many confounding variables. The current literature was reviewed to evaluate the strength of the associations and the causative effects of environmental factors on keratoconus. Ethnicity and consanguinity have been revealed as important determinants for geographical variations in keratoconus prevalence. Eye rubbing, atopy, floppy eyelid syndrome, contact lens wear, pregnancy, and thyroid hormone disturbances are likely associated with keratoconus. The first 4 factors can induce ocular surface inflammation, matrix metalloproteinase release, and keratocyte apoptosis, consistent with the postulated etiology of keratoconus. The associations of keratoconus with UV exposure, cigarette smoking, personality, and sex were less convincing once confounding factors were considered. Future studies powered for multivariate analysis of factors discussed will hopefully shed light on what is truly important in the development and progression of keratoconus.

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.

Exome sequencing identification of susceptibility genes in Chinese patients with keratoconus

PURPOSE

Keratoconus (KC) is a corneal ectasia disease with complex genetic heterogeneity. The present study aimed to identify susceptibility genes in Chinese patients with KC.

METHODS

Exome sequencing (ES) was performed in 28 Chinese KC patients to search for susceptibility genes of the disease. The candidate variants were filtered out by multi-step bioinformatics analysis and validated by Sanger sequencing. Another 100 individuals with KC were also recruited to verify those variants by Sanger sequencing.

RESULTS

By filtering out nonsynonymous variants located in exon, selecting variants which were presented in two or more samples and applying public databases to remove common variants, along with the inclusion of missense SNVs located in differential expressed genes and protein damaging variants (stop gain/stop loss SNVs and InDels), we have identified 6 SNVs (4 missense SNVs: c.1168 T > C in TRANK1, c.341A>T in ERMP1, c.4346 T > C in SDK2, c.1730A>C in COL6A1; 2 stop gain SNVs: c.1138 C > T in CNBD1, c.241 C > T in KRT82) and 2 InDels (c.193_195del in NSUN5, c.1690_1698del in COL9A3) as candidate variants for KC. The verifying results showed that c.341A>T in ERMP1 and c.193_195del in NSUN5 was found in one and two samples, respectively.

CONCLUSIONS

Our study suggested that a total of six SNVs in six genes and two InDels in two genes might be considered as candidate variants in Chinese patients with KC.

Keratoconus-susceptibility gene identification by corneal thickness genome-wide association study and artificial intelligence IBM Watson

Keratoconus is a common ocular disorder that causes progressive corneal thinning and is the leading indication for corneal transplantation. Central corneal thickness (CCT) is a highly heritable characteristic that is associated with keratoconus. In this two-stage genome-wide association study (GWAS) of CCT, we identified a locus for CCT, namely STON2 rs2371597 (P = 2.32 × 10⁻¹³), and confirmed a significant association between STON2 rs2371597 and keratoconus development (P = 0.041). Additionally, strong STON2 expression was observed in mouse corneal epithelial basal cells. We also identified SMAD3 rs12913547 as a susceptibility locus for keratoconus development using predictive analysis with IBM’s Watson question answering computer system (P = 0.001). Further GWAS analyses combined with Watson could effectively reveal detailed pathways underlying keratoconus development.

Yoshikatsu Hosoda et al. study the genetic basis for central corneal thickness (CCT) that is associated with keratoconus. They identify two susceptibility loci, STON2 rs2371597 and SMAD3 rs12913547, using two-step genome-wide association study (GWAS) and predictive analysis with IBM’s Watson question answering computer system, respectively.

Replication of the Association Between Keratoconus and Polymorphisms in PNPLA2 and MAML2 in a Han Chinese Population

Keratoconus (KC) is a complex ocular disease that is affected by both genetic and non-genetic triggers. A recent genome-wide association study (GWAS) identified a genome-wide significant locus for KC in the region of PNPLA2 (rs61876744), as well as a suggestive signal in the MAML2 (rs10831500) locus. In order to validate their findings, here we performed a replication study of the Han Chinese population, with 120 sporadic KC cases and 206 gender and age matched control subjects, utilizing the TaqMan SNP genotyping assays. SNP rs10831500, as well as two proxy SNPs for rs61876744, named rs7942159 and rs28633403, were subjected to genotyping. However, we did not find a significant difference (P > 0.05) in all the three genotyped SNPs between KC cases and the controls. A further meta-analysis on four previous cohorts of white patients and this Han Chinese cohort showed a significant genetic heterogeneity within the replicated loci. Thus, the current study suggests that SNP rs61876744 (or its proxy SNPs) and rs10831500 might not be associated with KC susceptibility in this Han Chinese cohort, and a large-scale association analysis focusing on the loci is therefore warranted in further investigations.

Absence of the c.169+50delTAAACAG mutation of SOD1 gene in a sample of keratoconus patients in Brazilian population

OBJECTIVE

To determine the presence of the 7-bp deletion c.169+50delTAAACAG in intron 2 of Superoxide Dismutase-1 gene in keratoconic patients from the State of São Paulo, Brazil, which promotes splicing variations, resulting in non-functional Superoxide Dismutase-1 antioxidant proteins, which may damage the corneal structure.

RESULTS

A group of 35 keratoconic patients, from whom 35 peripheral blood samples and 58 samples of corneal fragments were evaluated, and a control group of 89 individuals, from whom 41 blood samples and 149 samples of corneal fragments were collected. After the amplification of DNA fragments by polymerase chain reaction, mutational screening analysis was performed by enzymatic digestion, followed by direct sequencing. The absence of the 7-bp c.169+50delTAAACAG mutation in intron 2 of Superoxide Dismutase-1 gene was detected in the analyzed subjects of the 2 groups, both in the cornea and peripheral blood samples. Then, according to our results, there is no involvement of c.169+50delTAAACAG deletion in the pathogenesis of keratoconus in this population, once it was not detected. But we emphasize that studies involving this deletion must be continued in an attempt to elucidate this issue.

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.

Accumulation of sequence variants in genes of Wnt signaling and focal adhesion pathways in human corneas further explains their involvement in keratoconus

Background

Keratoconus (KTCN) is a protrusion and thinning of the cornea, resulting in loss of visual acuity. The etiology of KTCN remains unclear. The purpose of this study was to assess the potential involvement of new genetic variants in KTCN etiology based on both the genomic and transcriptomic findings recognized in the same corneal tissues.

Methods

Corneal tissues derived from five unrelated Polish individuals with KTCN were examined using exome sequencing (ES), followed by enrichment analyses. For comparison purposes, the datasets comprising ES data of five randomly selected Polish individuals without ocular abnormalities and five Polish patients with high myopia were used. Expression levels of selected genes from the overrepresented pathways were obtained from the previous RNA-Seq study.

Results

Exome capture discovered 117 potentially relevant variants that were further narrowed by gene overrepresentation analyses. In each of five patients, the assessment of functional interactions revealed rare (MAF ≤ 0.01) DNA variants in at least one gene from Wnt signaling (VANGL1, WNT1, PPP3CC, LRP6, FZD2) and focal adhesion (BIRC2, PAK6, COL4A4, PPP1R12A, PTK6) pathways. No genes involved in pathways enriched in KTCN corneas were overrepresented in our control sample sets.

Conclusions

The results of this first pilot ES profiling of human KTCN corneas emphasized that accumulation of sequence variants in several genes from Wnt signaling and/or focal adhesion pathways might cause the phenotypic effect and further points to a complex etiology of KTCN.

The keratoconus enigma: A review with emphasis on pathogenesis

PURPOSE

To review the literature on the etiopathogenesis of keratoconus (KC).

METHODS

A literature search was conducted using PUBMED and Google Scholar for keratoconus. The authors analyzed epidemiology studies, reviews, and case reports.

RESULTS

Atopy and ocular surface inflammation are a common features of KC and should lead to a reconsideration of the traditional definition of KC as a "non inflammatory" condition. Co-morbidities suggest that KC may be the ocular manifestation of a systemic disease. Finally, KC shows higher prevalence in certain ethnicities, which calls into question the status of KC as a rare disease, at least in these communities.

CONCLUSION

We believe that future studies should test whether selected, high prevalence populations exhibit specific genetic background and/or ethno-specific environmental risk factors.

Multiple Differentially Methylated Regions Specific to Keratoconus Explain Known Keratoconus Linkage Loci

Purpose

Keratoconus (KTCN) is a complex eye disorder resulting in loss of visual function. Its development is affected by genetic and environmental components. The aim of this study was to unravel the role of epigenetic factors in KTCN.

Methods

To verify if DNA methylation may play a role in KTCN development, reduced representation bisulfite sequencing of five KTCN and five non-KTCN human corneas was performed.

Results

Multiple KTCN-specific differentially methylated regions were detected and many of them overlap previously identified KTCN linkage loci (3p14.3, 5q35.2, 13q32.3, 15q24.1, and 20p13) and chromosome arms that have been linked to KTCN (2q, 4q, 5p, 9p, 14q, and 17q). Reanalysis of the previously described RNA sequencing dataset of 25 KTCN and 25 non-KTCN human corneas revealed that 12 genes downregulated in KTCN and 6 upregulated genes overlapped or were located in the near vicinity of the identified differentially methylated regions. Particularly interesting were the DNA methylation changes in WNT3 and WNT5A encoding Wnt ligands, as they provide a potential explanation for the Wnt signaling pathway dysregulation observed in KTCN.

Conclusions

We presented the results of data analysis from the first study of DNA methylation changes in human KTCN corneas compared to non-KTCN samples. We were able to identify genomic regions with distinct patterns of DNA hypo- and hypermethylation and link them to previously found KTCN susceptibility loci as well as transcriptomic disruption of Wnt signaling pathway observed in KTCN.

Genetic Variants Associated With Corneal Biomechanical Properties and Potentially Conferring Susceptibility to Keratoconus in a Genome-Wide Association Study

Importance

Keratoconus 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.

Objectives

To identify genetic associations with corneal biomechanical properties and to examine whether these genetic variants are associated with keratoconus.

Design, Setting, and Participants

A 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.

Exposures

In 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 Measures

In 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.

Results

Among 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 Relevance

Five 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.

[Search for genetic markers for precise diagnostics of keratoconus]

Keratoconus is a chronic disorder of the cornea, characterized by its progressive thinning, stretching, and conical protrusion. Diagnostics of subclinical keratoconus, as well as its early stages (forme fruste), is a complex problem. The presence of these forms of keratoconus in a patient is one of the reasons for the development of keratectasia after laser refractive surgery. Currently, the role of genetic factors in keratoconus development has been proven. This indicates the possibility of diagnostics of subclinical and forme fruste keratoconus using genetic markers. Knowledge about the patient's genetic susceptibility to keratoconus would allow correcting the tactics of treatment of refractive anomalies and avoiding serious side effects. The studies of causal mutations indicate the genetic heterogeneity of keratoconus, which complicates the development of a diagnostic panel. Selection of candidate variants from the currently known ones based on clear criteria may be one of the approaches for diagnostic markers search. In this review, we have analyzed articles on keratoconus markers in order to form a list of candidate variants for genotyping in the Russian population. The selection criteria took into account the complexes of symptoms in which a marker was found, populations in which a particular marker was investigated, the presence and results of replication studies. The analysis included markers in VSX1, SOD1, ZEB1, LOX, CAST, DOCK9, TGFBI, HGF, MAP3K19, KCND3, COL4A3, COL4A4, COL5A1, FNDC3B, FOXO1, BANP-ZNF469, MPDZ-NF1B, WNT10A genes. Based on the results of the analysis, the following candidate variants were selected for genotyping in the Russian population of patients with keratoconus: rs1536482 and rs7044529 in the COL5A1 gene, rs5745752 and rs2286194 in the HGF gene, rs4954218 in the MAP3K19 gene, rs4839200 near the KCND3 gene, rs2721051 near the FOXO1 gene, rs1324183 between the MPDZ and the NF1B genes, and rs121908120 in the WNT10A gene.

Genetic Variation in Human Gene Regulatory Factors Uncovers Regulatory Roles in Local Adaptation and Disease

Differences in gene regulation have been suggested to play essential roles in the evolution of phenotypic changes. Although DNA changes in cis-regulatory elements affect only the regulation of its corresponding gene, variations in gene regulatory factors (trans) can have a broader effect, because the expression of many target genes might be affected. Aiming to better understand how natural selection may have shaped the diversity of gene regulatory factors in human, we assembled a catalog of all proteins involved in controlling gene expression. We found that at least five DNA-binding transcription factor classes are enriched among genes located in candidate regions for selection, suggesting that they might be relevant for understanding regulatory mechanisms involved in human local adaptation. The class of KRAB-ZNFs, zinc-finger (ZNF) genes with a Krüppel-associated box, stands out by first, having the most genes located on candidate regions for positive selection. Second, displaying most nonsynonymous single nucleotide polymorphisms (SNPs) with high genetic differentiation between populations within these regions. Third, having 27 KRAB-ZNF gene clusters with high extended haplotype homozygosity. Our further characterization of nonsynonymous SNPs in ZNF genes located within candidate regions for selection, suggests regulatory modifications that might influence the expression of target genes at population level. Our detailed investigation of three candidate regions revealed possible explanations for how SNPs may influence the prevalence of schizophrenia, eye development, and fertility in humans, among other phenotypes. The genetic variation we characterized here may be responsible for subtle to rough regulatory changes that could be important for understanding human adaptation.

Association of KIF26B and COL4A4 gene polymorphisms with the risk of keratoconus in a sample of Iranian population

PURPOSE

Keratoconus (KTCN) is a congenital corneal eye disorder which correlates with abnormal distribution of the collagen fiber and causes loss of visual acuity. COLA4A gene has a substantive role in collagen synthesis, whereas KIF26B as a new candidate gene belonging to kinesin superfamily (KIFs) has been suggested to be associated with this disease. So, in this preliminary study, we simultaneously evaluated the effects of two single nucleotide polymorphisms, 222855rs7C/T and rs12407427C/T, on KTCN susceptibility in a sample of Iranian population.

METHODS

The present case-control study consists of 144 patients confirmed with KTCN and 153 healthy controls. The variants are genotyped by using amplification refractory mutation system-polymerase chain reaction method.

RESULTS

The findings disclosed that rs2228557C/T and rs12407427C/T polymorphisms significantly increased the risk of KTCN in measured (codominant1; p = 0.0001, codominant2; p = 0.0001, codominant3; p = 0.0006, dominant; p = 0.0001, over-dominant; p = 0.0005) and (codominant1; p = 0.0001, codominant3; p = 0.0005, recessive; p = 0.0001) inheritance patterns, respectively.

CONCLUSION

Our results did prove a statistical association of both rs2228557 and rs12407427 genotypes (TT and CT + CC) and allele (T) with KTCN susceptibility in Iranian population. Further studies in other ethnicities are required to verify our results.

Genetic Aspects of Keratoconus: A Literature Review Exploring Potential Genetic Contributions and Possible Genetic Relationships with Comorbidities

INTRODUCTION

Keratoconus (KC) is a complex, genetically heterogeneous, multifactorial degenerative disorder that is accompanied by corneal ectasia which usually progresses asymmetrically. With an incidence of approximately 1 per 2000 and 2 cases per 100,000 population presenting annually, KC follows an autosomal recessive or dominant pattern of inheritance and is, apparently, associated with genes that interact with environmental, genetic, and/or other factors. This is an important consideration in refractive surgery in the case of familial KC, given the association of KC with other genetic disorders and the imbalance between dizygotic twins. The present review attempts to identify the genetic loci contributing to the different KC clinical presentations and relate them to the common genetically determined comorbidities associated with KC.

METHODS

The PubMed, MEDLINE, Google Scholar, and GeneCards databases were screened for KC-related articles published in English between January 2006 and November 2017. Keyword combinations of "keratoconus," "risk factor(s)," "genetics," "genes," "genetic association(s)," and "cornea" were used. In total, 217 articles were retrieved and analyzed, with greater weight placed on the more recent literature. Further bibliographic research based on the 217 articles revealed another 124 relevant articles that were included in this review. Using the reviewed literature, an attempt was made to correlate genes and genetic risk factors with KC characteristics and genetically related comorbidities associated with KC based on genome-wide association studies, family-based linkage analysis, and candidate-gene approaches.

RESULTS

An association matrix between known KC-related genes and KC symptoms and/or clinical signs together with an association matrix between identified KC genes and genetically related KC comorbidities/syndromes were constructed.

CONCLUSION

Twenty-four genes were identified as potential contributors to KC and 49 KC-related comorbidities/syndromes were found. More than 85% of the known KC-related genes are involved in glaucoma, Down syndrome, connective tissue disorders, endothelial dystrophy, posterior polymorphous corneal dystrophy, and cataract.

An analytical enrichment-based review of structural genetic studies on keratoconus

Keratoconus is a progressive bilateral corneal protrusion that leads to irregular astigmatism and impairment of vision. Keratoconus is an etiologically heterogeneous corneal dystrophy and both environmental and genetic factors play a role in its etiopathogenesis. In this analytical review, we have studied all the genes that are structurally associated with keratoconus and have tried to explain the function of each gene and its association with other eye disorders in a concise way. In addition, using gene set enrichment analysis, it was attempted to find the most important impaired metabolic pathways in keratoconus. Several genetic studies have been carried out on keratoconus and several genes have been identified as risk factors involved in the etiology of the disease. In the current study, 16 studies, including nine association studies, five genome-wide association studies, one linkage study, and one meta-analysis, were reviewed and based on the 19 genes found, enrichment was performed and the most important metabolic pathways involved in the disease were identified. The enrichment results indicated that the two pathways, interleukin 1 processing and assembly of collagen fibrils, are significantly associated with the disease. Obviously, the results of this study, in addition to providing information about the genes involved in the disease, can provide an integrated insight into the gene-based etiology of keratoconus and therapeutic opportunities thereof.

Analysis of multiple genetic loci reveals MPDZ-NF1B rs1324183 as a putative genetic marker for keratoconus

OBJECTIVE

To investigate the associations between 16 single-nucleotide polymorphisms (SNPs) in 14 genetic loci and keratoconus in an independent Chinese cohort.

METHODS

This cross-sectional, case-control association study included a Chinese cohort of 133 patients with keratoconus and 371 control subjects. In a recent meta-analysis study, we identified association of 16 SNPs in 14 gene loci with keratoconus. In this study, we genotyped these 16 SNPs in all the patients and controls and analysed their association with keratoconus, its clinical severities and progression profiles. We also analysed the genotype-phenotype correlation between individual SNPs and steep keratometry, flat keratometry (Kf), average keratometry (Avg K) and best-fit sphere diameter (BFS) of the anterior and posterior corneal surface.

RESULTS

Among the 16 selected SNPs, rs1324183 in the MPDZ-NF1B locus showed a significant association with keratoconus (OR=2.22; 95% CI 1.42 to 3.45, p=4.30×10^-4), especially severe keratoconus (OR=5.10, 95% CI 1.63 to 15.93, p=0.005). The rs1324183 A allele was positively associated with anterior Kf (p=0.008), anterior Avg K (p=0.017), posterior Kf (p=0.01) and negatively associated with apex pachymetry (p=0.007) and anterior BFS (p=0.023) in keratoconus. The other 15 SNPs had no significant association with keratoconus or genotype-phenotype correlations.

CONCLUSIONS

This study confirmed the association of SNP rs1324183 in MPDZ-NF1B with keratoconus and revealed the association of this SNP with keratoconus severity and corneal parameters. It is thus a putative genetic marker for monitoring the progression of keratoconus to a severe form and facilitating early intervention.

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.

Corneal Cross-Linking for Pediatric Keratcoconus Review

PURPOSE

To comprehensively review the available published literature for cross-linking in the pediatric population.

METHODS

Review of the literature published in English in PubMed.

RESULTS

Two hundred ten publications were considered. One hundred fifteen were considered relevant to this review.

CONCLUSIONS

Studies of cross-linking in pediatric patients are sparse, with relatively short follow-up times, and mostly on small groups of patients. Treatment with cross-linking halts progression of keratoconus in the pediatric population, and early treatment seems to be cost-effective compared with later penetrating keratoplasty. Long-term effects and regression rates remain unclear, and further studies are needed in this population.

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.

RNA-Seq analysis and comparison of corneal epithelium in keratoconus and myopia patients

Keratoconus is a common degenerative corneal disease that can lead to significant visual morbidity, and both genetic and environmental factors have been implicated in its pathogenesis. We compared the transcriptome of keratoconus and control epithelium using RNA-Seq. Epithelial tissues were obtained prior to surgery from keratoconus and myopia control patients, undergoing collagen cross-linking and photorefractive keratectomy, respectively. We identified major differences in keratoconus linked to cell-cell communication, cell signalling and cellular metabolism. The genes associated with the Hedgehog, Wnt and Notch1 signaling pathways were down-regulated in keratoconus. We also identified plasmolipin and Notch1 as being significantly reduced in keratoconus for both gene and protein expression (p < 0.05). Plasmolipin is a novel protein identified in human corneal epithelium, and has been demonstrated to have a key role in epithelial cell differentiation in other tissues. This study shows altered gene and protein expression of these three proteins in keratoconus, and further studies are clearly warranted to confirm the functional role of these proteins in the pathogenesis of keratoconus.

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 - What is New?

Background:

Keratoconus is characterized as a bilateral, progressive, non-inflammatory thinning of the cornea resulting in blurred vision due to irregular astigmatism. Keratoconus has a multifactorial etiology, with multiple genetic and environmental components contributing to the disease pathophysiology. Several genomic loci and genes have been identified that highlight the complex molecular etiology of this disease.

Conclusion:

The review focuses on current knowledge of these genetic risk factors associated with keratoconus.

A review of keratoconus: Diagnosis, pathophysiology, and genetics

We discuss new approaches to the early detection of keratoconus and recent investigations regarding the nature of its pathophysiology. We review the current evidence for its complex genetics and evaluate the presently identified genes/loci and potential candidate gene/loci. In addition, we highlight current research methodologies that may be used to further elucidate the pathogenesis of keratoconus.

Genetic associations for keratoconus: a systematic review and meta-analysis

Genetic associations for keratoconus could be useful for understanding disease pathogenesis and discovering biomarkers for early detection of the disease. We conducted a systematic review and meta-analysis to summarize all reported genetic associations for the disease. We searched in the MEDLINE, Embase, Web of Science, and HuGENET databases for genetic studies of keratoconus published from 1950 to June 2016. The summary odds ratio and 95% confidence intervals of all polymorphisms were estimated using the random-effect model. Among 639 reports that were retrieved, 24 fulfilled required criteria as eligible studies for meta-analysis, involving a total of 53 polymorphisms in 28 genes/loci. Results of our meta-analysis lead to the prioritization of 8 single-nucleotide polymorphisms (SNPs) in 6 genes/loci for keratoconus in Whites. Of them 5 genes/loci were originally detected in genome-wide association studies, including FOXO1 (rs2721051, P = 5.6 × 10^-11), RXRA-COL5A1 (rs1536482, P = 2.5 × 10^-9), FNDC3B (rs4894535, P = 1.4 × 10^-8), IMMP2L (rs757219, P = 6.1 × 10^-7; rs214884, P = 2.3 × 10^-5), and BANP-ZNF469 (rs9938149, P = 1.3 × 10^-5). The gene COL4A4 (rs2229813, P = 1.3 × 10^-12; rs2228557, P = 4.5 × 10^-7) was identified in previous candidate gene studies. We also found SNPs in 10 genes/loci that had a summary P value < 0.05. Sensitivity analysis indicated that the results were robust. Replication studies and understanding the roles of these genes in keratoconus are warranted.

Novel Zinc Finger Protein Gene 469 (ZNF469) Variants in Advanced Keratoconus

PURPOSE

Common polymorphic variants upstream of Zinc finger protein gene 469 (ZNF469) have been associated with central corneal thickness. Rare ZNF469 variants have been shown in keratoconus patients. The aim of the current study was to investigate the frequency of ZNF 469 gene variants in rapidly progressive advance keratoconus patients who underwent corneal transplant surgery by the age of 30, compared to their frequency in the normal Turkish population.

METHODS

A search in a patient database was performed to identify patients with a rapidly progressive keratoconus requiring corneal transplant surgery by the age of 30 in at least one eye. Twenty-six advance keratoconus patients (study group) and 109 health subjects (control group) were included in the study. Blood samples were donated, and genomic DNA was extracted. The entire coding sequence of the ZNF469 gene including the 84 bp of the putative intron was amplified using PCR primers and analyzed using next generation sequencing (NGS).

RESULTS

Fifteen single nucleotide polymorphisms previously reported and registered to the dbSNP database were detected in the study group. The allele frequencies of these polymorphisms were higher in the keratoconus group compared to the control group and to the ExAC genome database. Three new missense heterozygote variants and one new synonym variant were detected in keratoconus group. According to prediction software, the P873T and Q2188H variants were shown to be non-tolerated, whereas G3424S could be tolerated. The synonymous variant R1060R is not predicted to lead to abnormal splicing by Human Splicing Finder in silico analysis.

CONCLUSION

New detected ZNF 469 P873T and Q2188H heterozygote coding variants in isolated advance keratoconus patients may be associated with the disease pathogenesis.

Non-invasive objective and contemporary methods for measuring ocular surface inflammation in soft contact lens wearers - A review

Contact lens wear is one of the primary risk factors for the development of ocular surface inflammatory events. The purpose of this review is to examine and summarize existing knowledge on the mechanisms of contact lens related ocular surface inflammation and the evidence for the effectiveness of current objective methods to measure ocular surface inflammation. Contact lens wear is postulated to trigger an inflammatory response on the ocular surface due to mechanical, chemical, hypoxic stress, or by the introduction of microbes and their toxins. Apart from the traditional signs of inflammation, such as swelling, oedema, redness and heat, on the ocular surface, other methods to measure ocular surface inflammation in sub-clinical levels include tear inflammatory mediator concentrations, conjunctival cell morphology, and corneal epithelial dendritic cell density and morphology. Tear inflammatory mediator concentrations are up- or down-regulated during contact lens wear, with or without the presence of associated inflammatory events. There is higher conjunctival cell metaplasia observed with contact lens wear, but changes in goblet cell density are inconclusive. Dendritic cell density is seen to increase soon after initiating soft contact lens wear. The long term effects of contact lens wear on dendritic cell migration in the cornea and conjunctiva, including the lid wiper area, require further investigation. Currently patient factors, such as age, smoking, systemic diseases and genetic profile are being studied. A better understanding of these mechanisms may facilitate the development of new management options and strategies to minimize ocular surface inflammation related to contact lens wear.

Replication of SNP associations with keratoconus in a Czech cohort

Introduction

Keratoconus is a relatively frequent disease leading to severe visual impairment. Existing therapies are imperfect and clinical management may benefit from improved understanding of mechanisms leading to this disease. We aim to investigate the replication of 11 single nucleotide polymorphisms (SNPs) with keratoconus.

Methods

SNPs from loci previously found in association with keratoconus were genotyped in 165 keratoconus cases of Caucasian Czech origin (108 males and 57 females) and 193 population and gender-matched controls. They included rs1536482 (COL5A1), rs4839200 (KCND3), rs757219 and rs214884 (IMMP2L), rs1328083 and rs1328089 (DAOA), rs2721051 (FOXO1), rs4894535 (FNDC3B), rs4954218 (MAP3K19, RAB3GAP1), rs9938149 (ZNF469) and rs1324183 (MPDZ). A case-control association analysis was assessed using Fisher’s exact tests.

Results

The strongest association was found for rs1324183 (allelic test OR = 1.58; 95% CI, 1.10–2.24, p = 0.01). Statistically significant values were also obtained for rs2721051 (allelic test OR = 1.72; 95% CI, 1.07–2.77, p = 0.025) and rs4954218 (allelic test OR = 1.53; 95% CI, 1.01–2.34; p = 0.047) which showed an opposite effect direction compared to previously reported one.

Conclusion

Independent replication of association between two SNPs and keratoconus supports the association of these loci with the risks for the disease development, while the effect of rs4954218 warrants further investigation. Understanding the role of the genetic factors involved in keratoconus etiopathogenesis may facilitate development of novel therapies and an early detection.

Genomic strategies to understand causes of keratoconus

Keratoconus (KTCN) is a degenerative disorder of the eye characterized by the conical shape and thinning of the cornea. The abnormal structure of KTCN-affected cornea results in loss of visual acuity. While many studies examine how environmental factors influence disease development, finding the genetic triggers has been a major emphasis of KTCN research. This paper focuses on genomic strategies that were implemented for finding candidate genes, including linkage and association studies, and presents different approaches of mutation screening. The advantages and limitations of particular tools are discussed based on literature and personal experience. Since etiology underlying KTCN is complex, numerous findings indicating heterogeneity of genetic factors involved KTCN etiology are presented.

Application of multiplex immunoassay technology to investigations of ocular disease

Eye-derived fluids, including tears, aqueous humour and vitreous humour often contain molecular signatures of ocular disease states. These signatures can be composed of cytokines, chemokines, growth factors, proteases and soluble receptors. However, the small quantities (<10 µl) of these fluids severely limit the detection of these proteins by traditional enzyme-linked immunosorbent assay or Western blot. To maximise the amount of information generated from the analysis of these specimens, many researchers have employed multiplex immunoassay technologies for profiling the expression or modification of multiple proteins from minute sample volumes.

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.