Secreted frizzled-related protein 1 (SFRP1) is highly upregulated in keratoconus epithelium: a novel finding highlighting a new potential focus for keratoconus research and treatment

Wnt/β-catenin signaling in corneal epithelium development, homeostasis, and pathobiology

The corneal epithelium is located on the most anterior surface of the eyeball and protects against external stimuli. The development of the corneal epithelium and the maintenance of corneal homeostasis are essential for the maintenance of visual acuity. It has been discovered recently via the in-depth investigation of ocular surface illnesses that the Wnt/β-catenin signaling pathway is necessary for the growth and stratification of corneal epithelial cells as well as the control of endothelial cell stability. In addition, the Wnt/β-catenin signaling pathway is directly linked to the development of common corneal illnesses such as keratoconus, fungal keratitis, and corneal neovascularization. This review mainly summarizes the role of the Wnt/β-catenin signaling pathway in the development, homeostasis, and pathobiology of cornea, hoping to provide new insights into the study of corneal epithelium and the treatment of related diseases.

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

OBJECTIVE

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

METHODS

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

RESULTS

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

CONCLUSION

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

Characteristics of Autophagy-Related Genes, Diagnostic Models, and Their Correlation with Immune Infiltration in Keratoconus

Purpose

Keratoconus (KTCN) is one of the most common degenerative keratopathies, significantly affecting vision and even leading to blindness. This study identifies potential biomarkers of KTCN based on the characterization of autophagy-related genes (ARGs) and the construction of a diagnostic model; and explores their relevance to immune infiltrating cells in KTCN.

Methods

Gene Expression Omnibus (GEO) data were downloaded and ARGs were acquired from GeneCards and Molecular Signatures Database (MSigDB). Autophagy-related differential expression genes (ARDEGs) were discovered through the integration of differentially expressed genes (DEGs) with ARGs, while hub genes of KTCN were discovered by protein-protein interaction (PPI) network analysis. The probable biological roles of these hub ARDEGs were examined using functional enrichment analysis, and a KTCN diagnostic model was generated using the least absolute shrinkage and selection operator (LASSO) regression analysis. We also employed the CIBERSORTx and ssGSEA algorithms to identify potential regulatory pathways to compare the abundance of immune cell infiltrates and their association with hub genes. Finally, the hub gene expression levels were confirmed using validation datasets as well as blood samples from KTCN and healthy individuals.

Results

In this study, we identified 12 hub ARDEGs, of which 9 genes were substantially distinct between KTCN patients and normal groups. The LASSO risk score was used to generate the nomogram, and the calibration curve evaluated the model's effective diagnostic performance (C index of 0.961). Patients with KTCN had greater percentages of M2 Macrophages and Gamma delta T cells, according to CIBERSORTx and ssGSEA. The outcomes of the bioinformatics analysis were supported by the DDIT3 and BINP3 expression levels in KTCN patients and healthy controls, according to the qRT-PCR data.

Conclusion

Five biomarkers (CFTR, PLIN2, DDIT3, BAG3, and BNIP3) and diagnostic models offer fresh perspectives on identifying and managing KTCN.

Comprehensive Transcriptome Analysis of Patients With Keratoconus Highlights the Regulation of Immune Responses and Inflammatory Processes

Keratoconus (KTCN), characterized by the steeper curvature of the cornea and the thinner central corneal thickness, was a degenerative corneal ectasia with ambiguous etiology and mechanism. We aim to reveal underlying pathological mechanisms of KTCN by multi-level transcriptomic, integrative omics analyses. We performed RNA-sequencing on corneal epithelial samples from seven patients and seven control donors, as well as peripheral matched blood samples from three of the patients and three control donors. After RNA extraction, library construction, and sequencing, differentially expressed genes and splicing events were identified, followed by over-representation analysis. In total, 547 differential expressed genes were identified in KTCN corneal epithelial samples, which were mainly enriched in immune responses and inflammatory processes. WGCNA module analysis, the transcriptomic analysis of peripheral blood samples, multiple omics data, and the meta-analysis of GEO samples confirmed the involvement of immune and inflammatory factors. Besides, 190 and 1,163 aberrant splicing events were identified by rMATS combined with CASH methods in corneal epithelial and blood samples with KCTN. In conclusion, this comprehensive transcriptome analysis of KTCN patients based on patients’ tissue and blood samples uncovered a significant association between immune-inflammatory genes and pathways with KCTN, highlighting the contribution of the perturbed immune signaling to the pathogenesis of KCTN. Our study suggested the importance of measures to control inflammation in the treatment of KCTN.

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.

Transcriptomic and Immunohistochemical Analysis of Progressive Keratoconus Reveal Altered WNT10A in Epithelium and Bowman's Layer

Purpose

To identify global gene expression changes in the corneal epithelium of keratoconus (KC) patients compared to non-KC myopic controls.

Methods

RNA-sequencing was performed on corneal epithelium samples of five progressive KC and five myopic control patients. Selected results were validated using TaqMan quantitative PCR (qPCR) on 31 additional independent samples, and protein level validation was conducted using western blot analysis on a subset. Immunohistochemistry was performed on tissue microarrays containing cores from over 100 KC and control cases. WNT10A transcript levels in corneal epithelium were correlated with tomographic indicators of KC disease severity in 15 eyes. Additionally, WNT10A was overexpressed in vitro in immortalized corneal epithelial cells.

Results

WNT10A was found to be underexpressed in KC epithelium at the transcript (ratio KC/control = 0.59, P = 0.02 per RNA-sequencing study; ratio = 0.66, P = 0.03 per qPCR) and protein (ratio = 0.07, P = 0.06) levels. Immunohistochemical analysis also indicated WNT10A protein was decreased in Bowman's layer of KC patients. In contrast, WNT10A transcript level positively correlated with increased keratometry (Kmax ρ = 0.57, P = 0.02). Finally, WNT10A positively regulated COL1A1 expression in corneal epithelial cells.

Conclusions

A specific Wnt ligand, WNT10A, is reduced at the mRNA and protein level in KC epithelium and Bowman's layer. This ligand positively regulates collagen type I expression in corneal epithelial cells. The results suggest that WNT10A expression in the corneal epithelium may play a role in progressive KC.

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.

Mapping Keratoconus Molecular Substrates by Multiplexed High-Resolution Proteomics of Unpooled Corneas

Keratoconus (KCN) is a leading cause for cornea grafting worldwide. Keratoconus is a multifactorial disease that causes progressive thinning of the cornea and whose etiology is poorly understood. Several studies have used proteomics on patient tear fluids to identify potential biomarkers. However, proteome of the cornea itself has not been investigated fully. We report here new findings from a case-control study using multiplexed mass spectrometry (MS) on individual (unpooled) corneas to gain deeper insights into proteins and biomarkers relevant to keratoconus. We employed a high-pressure approach to extract total protein from individual corneas from five cases and five controls, followed by trypsin digestion and tandem mass tag (TMT) labeling. The MS-derived data were searched using the Human NCBI RefSeq protein database v92, with peptides and proteins filtered at 1% false discovery rate. A total of 3132 proteins were detected, of which 627 were altered significantly (p ≤ 0.05) in keratoconus corneas. The increases were overwhelmingly in the mTOR/PI3/AKT signal-mediated regulations of cell survival and proliferation, nonsense-mediated decay of transcripts, and proteasomal pathways. The decreases were in several extracellular matrix proteins and in many members of the complement system. Importantly, this multiplexed proteomic study of keratoconus corneas identified, to our knowledge, the largest number of corneal proteins. The novel findings include changes in pathways that regulate transcript stability, proteasomal degradation, and the complement system in corneas with keratoconus. These observations offer new prospects toward future discovery of novel molecular targets for diagnostic and therapeutic innovations for patients with keratoconus.

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.

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.

Characterization of Corneal Epithelial Cells in Keratoconus

Purpose

We studied the cellular characteristics of epithelial cells in the cone and extraconal periphery of corneas in keratoconus eyes.

Methods

This prospective observational study was conducted at Narayana Nethralaya Eye Institute. A total of 83 and 42 eyes with keratoconus and normal topography, respectively, were included in the study. Corneal epithelial cells were collected and analyzed for apoptosis, proliferation, epithelial-mesenchymal transition, and differentiation status using molecular and biochemical tools. Statistical analysis was performed using the Student's t-test.

Results

Corneal epithelial cells from the cone showed significantly higher expression of proapoptotic marker BAX (P < 0.005) compared to controls. Significantly elevated expression of cell cycle markers CYCLIN D1 (P < 0.005) and Ki67 (P < 0.005) were noted in the extraconal region compared to controls. Cells of the cone showed significantly higher ZO-1 (P < 0.005) and lower vimentin (P < 0.005) compared to controls. Significantly lower expression of the differentiation marker CK3/12 (P < 0.05) was observed in cones compared to controls.

Conclusions

Cones of keratoconic corneas show enhanced cell death, poor differentiation, proliferation and epithelial-mesenchymal transition. The cellular changes of the corneal epithelial cells in the cone and extraconal region differ significantly in a keratoconus corneas.

Translational Relevance

Characterization of patient-specific corneal epithelial cellular status in keratoconus has the potential to determine the optimal treatment and therapeutic outcomes paving the way towards personalized treatment in the future.

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.

Modeling Keratoconus Using Induced Pluripotent Stem Cells

PURPOSE

To model keratoconus (KC) using induced pluripotent stem cells (iPSC) generated from fibroblasts of both KC and normal human corneal stroma by a viral method.

METHODS

Both normal and KC corneal fibroblasts from four human donors were reprogramed directly by delivering reprogramming factors in a single virus using 2A "self-cleaving" peptides, using a single polycistronic lentiviral vector coexpressing four transcription factors (Oct 4, Sox2, Klf4, and Myc) to yield iPSC. These iPS cells were characterized by immunofluorescence detection using of stem cell markers (SSEA4, Oct4, and Sox2). The mRNA sequencing was performed and the datasets were analyzed using ingenuity pathways analysis (IPA) software.

RESULTS

The generated stem cell-like clones expressed the pluripotency markers, SSEA4, Oct4, Sox2, Tra-1-60, and also expressed pax6. Our transcriptome analysis showed 4300 genes, which had 2-fold change and 870 genes with a q-value of <0.05 in keratoconus iPSC compared to normal iPSC. One of the genes that showed difference in KC iPSC was FGFR2 (down-regulated by 2.4 fold), an upstream target of Pi3-Kinase pathway, was further validated in keratoconus corneal sections and also KC iPSC-derived keratocytes (down regulated by 2.0-fold). Both normal and KC-derived keratocytes expressed keratocan, signature marker for keratocytes. KC iPSC-derived keratocytes showed adverse growth and proliferation and was further confirmed by using Ly2924002, a PI3k inhibitor, which severely affected the growth and differentiation in normal iPSC.

CONCLUSIONS

Based on our result, we propose a model for KC in which inhibition FGFR2-Pi3-Kinase pathway affects the AKT phosphorylation, and thus affecting the keratocytes survival signals. This inhibition of the survival signals could be a potential mechanism for the KC-specific decreased cell survival and apoptosis of keratocytes.

Proteomic and gene expression patterns of keratoconus

Keratoconus is a progressive corneal thinning disease associated with significant tissue remodeling activities and activation of a variety of signaling networks. However, it is not understood how differential gene and protein expression direct function in keratoconus corneas to drive the underlying pathology, ectasia. Research in the field has focused on discovering differentially expressed genes and proteins and quantifying their levels and activities in keratoconus patient samples. In this study, both microarray analysis of total ribonucleic acid (RNA) and whole proteome analyses are carried out using corneal epithelium and tears from keratoconus patients and compared to healthy controls. A number of structural proteins, signaling molecules, cytokines, proteases, and enzymes have been found to be deregulated in keratoconus corneas. Together, the data provide clues to the complex process of corneal degradation which suggest novel ways to clinically diagnose and manage the disease. This review will focus on discussing these recent advances in the knowledge of keratoconus biology from a gene expression and function point-of-view.

Increased expression of secreted frizzled-related protein-1 and microtubule-associated protein light chain 3 in keratoconus

PURPOSE

To study the expression of secreted frizzled-related protein-1 (SFRP-1) and microtubule-associated protein light chain 3 (LC3), an autophagy marker, in keratoconus.

METHODS

Under an institutional review board-approved protocol, de-identified and/or surgically discarded normal donor (n = 10) and keratoconus corneas (n = 10) were obtained. The corneal samples were fixed in formalin and embedded in paraffin. Immunohistochemical staining using SFRP-1 and LC3 antibodies was performed.

RESULTS

The majority of expression of SFRP-1 was seen in the epithelium; however, in 3 tissues that showed high expression, staining was also present in the stroma and endothelium. Like SFRP-1, the LC3 expression in keratoconus tissues occurred at 3 different levels: low, medium, and high. Collectively these data suggest that there are differences in the expression levels of SFRP-1 and LC3 in keratoconus tissue compared with the normal tissue. Low expressivity of SFRP-1 seemed to correspond to low expressivity of LC3, whereas medium to high expressivity of SFRP-1 corresponded to medium to high expressivity of LC3.

CONCLUSIONS

Increased expression of SFRP-1 and LC3 was observed in keratoconus corneas. Keratocyte autophagy seen with keratoconus may play a role in the pathogenesis of keratoconus.

Expression of SFRP Family Proteins in Human Keratoconus Corneas

We investigated the expression of the secreted frizzled-related proteins (SFRPs) in keratoconus (KC) and control corneas. KC buttons (∼8 mm diameter) (n = 15) and whole control corneas (n = 7) were fixed in 10% formalin or 2% paraformaldehyde and subsequently paraffin embedded and sectioned. Sections for histopathology were stained with hematoxylin and eosin, or Periodic Acid Schiff's reagent. A series of sections was also immunolabelled with SFRP 1 to 5 antibodies, visualised using immunofluorescence, and examined with a Zeiss LSM700 scanning laser confocal microscope. Semi-quantitative grading was used to compare SFRP immunostaining in KC and control corneas. Overall, KC corneas showed increased immunostaining for SFRP1 to 5, compared to controls. Corneal epithelium in all KC corneas displayed heterogeneous moderate to strong immunoreactivity for SFRP1 to 4, particularly in the basal epithelium adjacent to cone area. SFRP3 and 5 were localised to epithelial cell membranes in KC and control corneas, with increased SFRP3 cytoplasmic expression observed in KC. Strong stromal expression of SFRP5, including extracellular matrix, was seen in both KC and control corneas. In control corneas we observed differential expression of SFRP family proteins in the limbus compared to more central cornea. Taken together, our results support a role for SFRPs in maintaining a healthy cornea and in the pathogenesis of epithelial and anterior stromal disruption observed in KC.

Study of VSX1 mutations in patients with keratoconus in southwest Iran using PCR-single-strand conformation polymorphism/heteroduplex analysis and sequencing method

OBJECTIVE

Keratoconus (KC) is an eye disorder in which the cornea is swollen, thinned and deformed. Despite extensive studies, the pathophysiological processes and genetic etiology of KC are unknown. The disease incidence is approximately 1 in 2,000, and it is the most common cause of corneal transplantation in the USA. Many genes are involved in the disease, but evidence suggests a major role for VSX1 in the etiology of KC. This study aimed to determine the frequency of mutations in exons 2, 3 and 4 of the VSX1 gene in Chaharmahal va Bakhtiari province in the southwest of Iran.

STUDY DESIGN

In this experimental study, mutations in 3 exons, namely exons 2, 3 and 4, of VSX1 were investigated in 50 patients with KC and 50 healthy control subjects. DNA was extracted using a standard phenol-chloroform method. PCR-single-strand conformational polymorphism/heteroduplex analysis was performed, followed by DNA sequencing to confirm the identified motility shifts.

RESULTS

H244R mutations were found in 1 patient and also in 1 healthy control subject. Furthermore, 12 polymorphisms were identified in patients with KC and 7 in healthy control subjects [rs6138482 and c.546A>G (rs12480307)].

CONCLUSION

Our investigation showed that KC-related VSX1 mutations were found in a very small proportion of the studied patients from Iran. Further investigations on other genes are needed to clarify their roles in KC pathogenesis.