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.

The Association of a Single Nucleotide Variant in COL5A1 to Early Onset Keratoconus and Pectus Excavatum-Convergence of Extracellular Matrix Pathologies.

Keratoconus is a bilateral ocular condition characterized by irregularities and the thinning of the cornea. Decreased central corneal thickness is a hallmark of the condition, and numerous genes have played a role in altering corneal thickness and the subsequent development of keratoconus. Variants in the structural and regulatory genes of the extracellular matrix have been highly associated with keratoconus, as well as with pectus excavatum, a chest wall deformity commonly seen in connective tissue disorders. This report describes a patient with a c.1720-11T>A intronic variant in the collagen-encoding gene, COL5A1, who was diagnosed with early-onset keratoconus and demonstrated a significant pectus excavatum. This report associates a COL5A1 variant with these seemingly unrelated phenotypic associations, further advancing the literature on the topic.

Identification of Keratoconus-Related Phenotypes in Three Ppip5k2 Mouse Models.

Purpose: It is necessary to establish a mouse model of keratoconus (KC) for research and therapy. We aimed to determine corneal phenotypes in 3 Ppip5k2 mouse models. Methods: Central corneal thickness (CCT) was determined using spectral domain optical coherence tomography (SD-OCT) in Ppip5k2+/K^ (n = 41 eyes), Ppip5k2K^/K^ (n = 17 eyes) and 2 knock-in mice, Ppip5k2S419A/+ (n = 54 eyes) and Ppip5k2S419A/S419A (n = 18 eyes), and Ppip5k2D843S/+ (n = 42 eyes) and Ppip5k2D843S/D843S (n = 44 eyes) at 3 and 6 months. Pachymetry maps were generated using the Mouse Corneal Analysis Program (MCAP) to process OCT images. Slit lamp biomicroscopy was used to determine any corneal abnormalities, and, last, hematoxylin and eosin (H&E) staining using corneal sections from these animals was used to examine morphological changes. Results: CCT significantly decreased from 3 to 6 months in the Ppip5k2+/K^ and Ppip5k2K^/K^ mice compared to their littermate controls. OCT-based pachymetry maps revealed abnormally localized thinning in all three models compared to their wild-type (WT) controls. Slit lamp examinations revealed corneal abnormalities in the form of bullous keratopathy, stromal edema, stromal scarring, deep corneal neovascularization, and opacities in the heterozygous/homozygous mice of the three models in comparison with their controls. Corneal histological abnormalities, such as epithelial thickening and stromal layer damage, were observed in the heterozygous/homozygous mice of the three models in comparison with the WT controls. Conclusions: We have identified phenotypic and histological changes in the corneas of three mouse lines that could be relevant in the development of animal models of KC.

Bioinformatics analysis of signature genes related to cell death in keratoconus.

Keratoconus is corneal disease in which the progression of conical dilation of cornea leads to reduced visual acuity and even corneal perforation. However, the etiology mechanism of keratoconus is still unclear. This study aims to identify the signature genes related to cell death in keratoconus and examine the function of these genes. A dataset of keratoconus from the GEO database was analysed to identify the differentially expressed genes (DEGs). A total of 3558 DEGs were screened from GSE151631. The results of Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that they mainly involved in response to hypoxia, cell-cell adhesion, and IL-17 signaling pathway. Then, the cell death-related genes datasets were intersected with the above 3558 DEGs to obtain 70 ferroptosis-related DEGs (FDEGs), 32 autophagy-related DEGs (ADEGs), six pyroptosis-related DEGs (PDEGs), four disulfidptosis-related DEGs (DDEGs), and one cuproptosis-related DEGs (CDEGs). After using Least absolute shrinkage and selection operator (LASSO), Random Forest analysis, and receiver operating characteristic (ROC) curve analysis, one ferroptosis-related gene (TNFAIP3) and five autophagy-related genes (CDKN1A, HSPA5, MAPK8IP1, PPP1R15A, and VEGFA) were screened out. The expressions of the above six genes were significantly decreased in keratoconus and the area under the curve (AUC) values of these genes was 0.944, 0.893, 0.797, 0.726, 0.882 and 0.779 respectively. GSEA analysis showed that the above six genes mainly play an important role in allograft rejection, asthma, and circadian rhythm etc. In conclusion, the results of this study suggested that focusing on these genes and autoimmune diseases will be a beneficial perspective for the keratoconus etiology research.

The Chinese keratoconus (CKC) cohort study.

The Chinese keratoconus (CKC) cohort study is a population-based longitudinal prospective cohort study in the Chinese population involving a clinical database and biobanks. This ongoing study focuses on the prevention of KC progression and is the first to involve the effect of gene‒environment interactions on KC progression. The CKC cohort is hospital-based and dynamic and was established in Zhengzhou, China; KC patients (n = 1114) from a large geographical area were enrolled from January 2019 to June 2023, with a mean age of 22.23 years (6‒57 years). Demographic details, socioeconomic characteristics, lifestyle, disease history, surgical history, family history, and visual and social function data are being collected using questionnaires. General physical examination, eye examination, biological specimen collection, and first-degree relative data were collected and analyzed in the present study. The primary focus of the present study was placed on gene, environment and the effect of gene‒environment interactions on KC progression. The follow-up of the CKC cohort study is expected to include data collection at 3 months, 6 months, and 1 year after the initial examination and then at the annual follow-up examinations. The first follow-up of the CKC cohort study was recorded. A total of 918 patients completed the follow-up by June 1, 2023, with a response rate of 82.40%. Aside from the younger age of patients who were followed up, no significant differences were found between patients who were followed up and patients who were not.

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.

The candidate proteins associated with keratoconus: A meta-analysis and bioinformatic analysis.

PURPOSE: Keratoconus (KC) is a multifactorial disorder. This study aimed to conduct a systematic meta-analysis to exclusively explore the candidate proteins associated with KC pathogenesis. METHODS: Relevant literature published in the last ten years in Pubmed, Web of Science, Cochrane, and Embase databases were searched. Protein expression data were presented as the standard mean difference (SMD) and 95% confidence intervals (CI). The meta-analysis is registered on PROSPERO, registration number CRD42022332442 and was conducted in accordance with the Preferred Reporting Items for Systematic reviews and Meta-Analyses statement (PRISMA). GO and KEGG enrichment analysis were performed, as well as the miRNAs and chemicals targeting the candidate proteins were predicted. PPI was analyzed to screen the hub proteins, and their expression was verified by RT-qPCR. RESULTS: A total of 21 studies were included in the meta-analysis, involving 346 normal eyes and 493 KC eyes. 18 deregulated proteins with significant SMD values were subjected to further analysis. In which, 7 proteins were up-regulated in KC compared with normal controls, including IL6 (SMD 1.54, 95%CI [0.85, 2.24]), IL1B (SMD 2.07, 95%CI [0.98, 3.16]), TNF (SMD 2.1, 95%CI [0.24, 3.96]), and MMP9 (SMD 1.96, 95%CI [0.68, 3.24]). While 11 proteins were down-regulated in KC including LOX (SMD 2.54, 95%CI [-4.51, -0.57]). GO and KEGG analysis showed that the deregulated proteins were involved in inflammation, extracellular matrix (ECM) remodeling, and apoptosis. MMP9, IL6, LOX, TNF, and IL1B were regarded as hub proteins according to the PPI analysis, and their transcription changes in stromal fibroblasts of KC were consistent with the results of the meta-analysis. Moreover, 10 miRNAs and two natural polyphenols interacting with hub proteins were identified. CONCLUSION: This study obtained 18 candidate proteins and demonstrated altered cytokine profiles, ECM remodeling, and apoptosis in KC patients through meta-analysis and bioinformatic analysis. It will provide biomarkers for further understanding of KC pathogenesis, and potential therapeutic targets for the drug treatment of KC.

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.

Impact of a 50bp insertion/deletion polymorphism of the superoxide dismutase-1 on oxidative stress status and risk of keratoconus.

Keratoconus (KC) is characterized by the predominant primary ectatic disease, affecting the cornea, necessitating corneal transplants in some cases. While some loci associated with KC risk have been identified, the understanding of the disease remains limited. Superoxide dismutase (SOD) enzymes play a crucial role in countering the reactive oxygen species and providing protection against oxidative stress (OS). Accordingly, the objective of this study was to investigate a potential association of a 50 nucleotide base pairs (bp) insertion/deletion (I/D) within the SOD1 promoter, and the located 1684 bp upstream of the SOD1 ATG, with KC in the Iranian population. Additionally, an assessment was conducted on SOD activity and the total antioxidant capacity (TAC), as determined by the ferric reducing-antioxidant power assay, along with malondialdehyde (MDA) levels. In this case-control study, genomic DNA was extracted from the blood cells of KC (n = 402) and healthy (n = 331) individuals. The genotype of this gene was determined using the PCR technique. Furthermore, the amount of SOD enzyme activity and the MDA and TAC levels were measured in the serum of the study groups. The (I/I) genotype was present in 84.23%, the (I/D) genotype in 15.06%, and the (D/D) genotype in 0.69% of both groups. A statistically significant relationship was seen between different genotypes and TAC, MDA, and SOD1 activity indices (P < 0.05). Individuals with the D/D genotype exhibited a decrease in total antioxidant capacity, an increase in the amount of MDA, and a decrease in SOD1 enzyme activity (P < 0.05). Moreover, the logistic regression analysis of KC development indicated that elevated levels of MDA increased the risk of KC incidence in the patient group compared to the healthy group, while a higher activity of SOD1 and greater values of TAC decreased the KC risk. The removal of the 50 bp fragment reduced SOD1 activity and elevated OS levels, thereby impacting the oxidant-antioxidant balance. This could potentially play a significant role in individuals afflicted by KC.

Progressive Keratoconus in a Patient With Severe Pectus Excavatum and a Cartilage Oligomeric Matrix Protein Gene Mutation: A Case Report.

INTRODUCTION: Keratoconus is a progressive ocular disorder associated with numerous systemic diseases, many of which affect the musculoskeletal system. Although the etiology and pathophysiology of the disorder remain elusive, recent studies suggest a significant role of genetic predisposition in the pathogenesis of keratoconus. This case report aims to elucidate a potential genetic association in a patient presenting with keratoconus, severe pectus excavatum, generalized muscular weakness, and skeletal deformities. CASE DESCRIPTION: A 31-year-old Iranian man presented with progressively diminishing vision in both eyes over the years, eventually diagnosed with keratoconus. The patient's history and further examination indicated generalized muscular weakness, skeletal deformities, and severe pectus excavatum with cardiac and large vessel displacement. Whole-exome sequencing identified two heterozygous gene variants: one in the Cartilage Oligomeric Matrix Protein (COMP) gene and another in the Regulating Synaptic Membrane Exocytosis 1 gene. The patient's systemic and ocular symptoms, combined with the gene variants identified, suggested a connective tissue systemic disorder, potentially within the clinical spectrum of COMPopathies. CONCLUSION: This is the first documented case of bilateral progressive keratoconus associated with severe pectus excavatum, generalized musculoskeletal dystrophy, and a COMP gene mutation. It highlights the necessity of continued search into the pathogenic genes of keratoconus, particularly in cases with coexisting systemic manifestations, to further our understanding of the etiology and pathogenesis of this complex disease.

Circular RNA expression profile identifies potential circulating biomarkers for keratoconus.

Early diagnosis is important for improving the outcomes of keratoconus (KC). Stable expression and a closed-loop structure of circular RNAs (circRNAs) make them ideal for the diagnosis and treatment of diseases. However, the expression pattern and potential function of circRNAs in KC is not studied yet. Hence, this study explored the circRNA expression profile of KC corneas through transcriptome sequencing and circRNA expression profile analysis. The diagnostic potential of blood circRNAs for KC was explored by analysing the circRNAs' expression levels of fifty paired blood samples from patients with KC and normal controls. The results showed that 107 significantly upregulated and 145 significantly downregulated circRNAs (|fold change| ≥ 2.0, p-value <0.05) were identified in KC tissues. Eight top differently expressed circRNAs were further validated in more cornea samples. Among them, five circRNAs expressed in peripheral blood, and four circRNAs (circ_0006156, circ_0006117, circ_0000284 and circ_0001801) showed significant downregulation in KC patients' peripheral blood too. The blood circ_0000284 expression levels of early, moderate, and advanced KC patients both were significantly lower than the controls. The blood circ_0006117 expression levels present a positive correlation with corrected distance visual acuity values, and a negative correlation with back elevation values of KC eyes. Notably, the expression levels of these circRNAs distinguished KC patients from their healthy counterparts, with the area under the curve (AUC) of circ_0000284, circ_0001801, and circ_0006117 being 0.7306, 0.6871 and 0.6701, respectively. Further, the AUC value for five circRNAs under the logistic regression model was 0.8203, indicating that they can function as effective biomarkers for the KC diagnostics. In conclusion, the expression of circRNAs showed a relationship with KC, with four significantly differentially expressed circRNAs demonstrating potential as biomarkers for the disease.

Animal Models for the Study of Keratoconus.

Keratoconus (KC) is characterized by localized, central thinning and cone-like protrusion of the cornea. Its precise etiology remains undetermined, although both genetic and environmental factors are known to contribute to disease susceptibility. Due to KC's complex nature, there is currently no ideal animal model to represent both the corneal phenotype and underlying pathophysiology. Attempts to establish a KC model have involved mice, rats, and rabbits, with some additional novel animals suggested. Genetic animal models have only been attempted in mice. Similarly, spontaneously occurring animal models for KC have only been discovered in mice. Models generated using chemical or environmental treatments have been attempted in mice, rats, and rabbits. Among several methods used to induce KC in animals, ultraviolet radiation exposure and treatment with collagenase are some of the most prevalent. There is a clear need for an experimental model animal to elucidate the underlying mechanisms behind the development and progression of keratoconus. An appropriate animal model could also aid in the development of treatments to slow or arrest the disorder.

Association between Polymorphism rs61876744 in PNPLA2 Gene and Keratoconus in a Saudi Cohort.

The genetic etiology of Keratoconus (KC) in Middle Eastern Arabs of Saudi origin is still unclear. A recent genome-wide study identified two significant loci in the region of PNPLA2 (rs61876744) and CSNK1E (rs138380) for KC that may be associated with KC in the Saudi population. In addition, polymorphisms in the apolipoprotein E (APOE) gene, namely, rs429358 and rs7412, responsible for APOE allelic variants ε2, ε3, and ε4, may influence KC via oxidative stress mechanism(s). Thus, we investigated the possible association of polymorphisms rs61876744, rs138380, rs429358, rs7412, and APOE genotypes in KC patients of the Saudi population. This study included 98 KC cases and 167 controls. Polymorphisms rs6187644 and rs138380 were genotyped using TaqMan assays, and rs429358 and rs7412 were genotyped via Sanger sequencing. Although the allele frequency of rs61876744(T) in PNPLA2 was a protective effect against KC (odds ratio (OR) = 0.64, 95% confidence interval (CI) = 0.44-0.93), the p-value (p = 0.020) was not significant for multiple testing correction (p = 0.05/4 = 0.015). However, rs6187644 genotype showed a modestly significant protective effect in the dominant model (OR = 0.53, 95% CI = 0.32-0.88, p = 0.013). Polymorphisms rs138380, rs429358, and rs7412 showed no significant allelic or genotype association with KC. However, the ε2-carriers (ε2/ε2 and ε2/ε3 genotypes) exhibited a greater than 5-fold increased risk of KC, albeit non-significantly (p = 0.055). Regression analysis showed no significant effect of age, gender, and the four polymorphisms on KC. Our results suggest that polymorphism rs6187644 in PNPLA2 might be associated with KC in the Middle Eastern Arabs of Saudi origin but warrant a large-scale association analysis at this locus.

Identification of the immune-associated characteristics and predictive biomarkers of keratoconus based on single-cell RNA-sequencing and bulk RNA-sequencing.

Background: Whether keratoconus (KC) is an inflammatory disease is currently debated. Hence, we aimed to investigate the immune-related features of KC based on single-cell RNA sequencing (scRNA-seq) and bulk RNA sequencing (bulk RNA-seq) data. Methods: scRNA-seq data were obtained from the Genome Sequence Archive (GSA), bulk RNA-seq data were obtained from the Gene Expression Omnibus (GEO), and immune-associated genes(IAGs) were obtained from the ImmPort database. Cell clusters of KC were annotated, and different cell clusters were then selected. The IAG score of each cell was calculated using the AUCell package. Three bulk RNA-seq datasets were merged and used to identify the differentially expressed genes (DEGs), biological functions, and immune characteristics. Weighted gene coexpression network analysis (WGCNA) was used to select the IAG score-related hub genes. Based on scRNA-seq and bulk RNA-seq analyses, three machine learning algorithms, including random forest (RF), support vector machine (SVM), and least absolute shrinkage and selection operator (LASSO) regression analysis, were used to identify potential prognostic markers for KC. A predictive nomogram was developed based on prognostic markers. Results: Six cell clusters were identified in KC, and decreased corneal stromal cell-5 (CSC-5) and increased CSC-6 were found in KC. CSC and immune cell clusters had the highest IAG scores. The bulk RNA-seq analysis identified 1362 DEGs (553 upregulated and 809 downregulated) in KC. We found different immune cell populations and differentially expressed cytokines in KC. More than three key IAG score-related modules and 367 genes were identified. By integrating the scRNA-seq and bulk RNA-seq analyses, 250 IAGs were selected and then incorporated into three machine learning models, and 10 IAGs (CEP112, FYN, IFITM1, IGFBP5, LPIN2, MAP1B, RNASE1, RUNX3, SMIM10, and SRGN) were identified as potential prognostic genes that were significantly associated with cytokine and matrix metalloproteinase(MMP)1-14 expression. Finally, a predictive nomogram was constructed and validated. Conclusion: Taken together, our results identified CSCs and immune cell clusters that may play a key role during KC progression by regulating immunological features and maintaining cell stability.

Integrated analysis of murine cornea identifies JAK/STAT signaling pathway upregulated specifically in female Vitamin A Deficient mice.

The Keratoconus (KC) is a corneal ectatic disease with unclear etiology. There are increasing studies that reported its association with a variety of inflammatory mechanisms. Vitamin A(VA) is an important nutrient related to inflammation regulation, and its deficiency may cause abnormalities of the ocular surface. However, the proportion of Vitamin A deficiency(VAD) was found surprisingly high among KC patients in our clinic practice. The aim of this study is to explore the effects of VAD on the transcriptome of corneas with the help of the VAD murine model and transcriptomics techniques. Blood samples of KC patients and non-KC controls (NC) were collected and the serum VA concentrations were measured and analyzed. A total of 52 NC and 39 KC were enrolled and the comparison of serum VA showed that the proportion of VAD in KC patients was 48.7% versus 1.9% in NC group. The further analysis of gender differences showed the proportion of VAD in female KC was 88.9% versus 36.7% in KC male patients. To explore the influence of VAD on cornea, the VAD mice fed with VAD diets were used. The RNA sequencing was employed to compare the corneal transcriptomic characteristics between the VAD female mice, NC female mice, VAD male mice and NC male mice. The transcriptome analysis revealed that the upregulated differential genes were mainly enriched in the immune response related pathways in VAD female mice versus NC female mice, especially the genes of JAK-STAT signaling pathway. The downstream molecules of JAK-STAT pathway were also significant after corneal mechanical scratching in female VAD mice. While, the differential genes between VAD male mice and NC male mice were estrogen signaling pathway instead of JAK-STAT pathway. This study indicates that VAD affects the transcriptomics of murine cornea with gender differences, which specifically affects the inflammatory status of the female murine cornea.

Identification of genetic variants in two families with Keratoconus.

BACKGROUND: This research investigated the genetic characteristic of two Chinese families with keratoconus (KC). METHODS: For all people in the two families with KC, their history, clinical data, and peripheral blood were collected. One hundred healthy participants without KC and 112 sporadic KC patients were recruited as the controls. Whole exome sequencing of the genomic DNA and polymerase chain reaction were conducted for all the controls and family members to verify the variants. Functional analyses of the variants was performed using the software programs. RESULTS: A missense tuberous sclerosis 1 (TSC1) variant g.135797247A > G (c.622A > G, p.Ser208Gly) was detected in family 1. A single nucleotide polymorphism (SNP) rs761232139 (p.Gly235Arg) in aldehyde dehydrogenase 3 family member A1 (ALDH3A1) gene was detected in family 2. The variant c.622A > G in TSC1 and the SNP rs761232139 in ALDH3A1 were predicted as being probably damaging. CONCLUSIONS: Novel variant c.622A > G in TSC1 and SNP rs761232139 in ALDH3A1 have been detected in families with KC. These two findings may play a role in the pathogenesis of KC.

Whole-Exome Sequencing of 24 Spanish Families: Candidate Genes for Non-Syndromic Pediatric Keratoconus.

Keratoconus is a corneal dystrophy that is one of the main causes of corneal transplantation and for which there is currently no effective treatment for all patients. The presentation of this disease in pediatric age is associated with rapid progression, a worse prognosis and, in 15-20% of cases, the need for corneal transplantation. It is a multifactorial disease with genetic variability, which makes its genetic study difficult. Discovering new therapeutic targets is necessary to improve the quality of life of patients. In this manuscript, we present the results of whole-exome sequencing (WES) of 24 pediatric families diagnosed at the University Hospital La Paz (HULP) in Madrid. The results show an oligogenic inheritance of the disease. Genes involved in the structure, function, cell adhesion, development and repair pathways of the cornea are proposed as candidate genes for the disease. Further studies are needed to confirm the involvement of the candidate genes described in this article in the development of pediatric keratoconus.

Exosomes and their miRNA/protein profile in keratoconus-derived corneal stromal cells.

Keratoconus (KC) is a corneal thinning disorder and a leading cause of corneal transplantation worldwide. Exosomes are small, secreted extracellular vesicles (30-150 nm) that mediate cellular communication via their protein, lipid, and nucleic acid content. We aimed to characterize the exosomes secreted by primary corneal fibroblasts from subjects with or without KC. Using human keratoconus stromal fibroblast cells (HKC, n = 4) and healthy stromal fibroblasts (HCF, n = 4), we collected and isolated exosomes using serial ultracentrifugation. Using nanoparticle tracking analysis (NTA) with ZetaView®, we compared the size and concentration of isolated exosomes. Different exosomal markers were identified and quantified using a transmission electron microscope (TEM) (CD81) and Western blot (CD9 and CD63). Exosomal miRNA profiles were determined by qRT-PCR using Exiqon Human panel I miRNA assays of 368 pre-selected miRNAs. Proteomic profiles were determined using a label-free spectral counting method with mass spectrometry. Differential expression analysis for miRNAs and proteins was done using student's t-test with a significance cutoff of p-value ≤0.05. We successfully characterized exosomes isolated from HCFs using several complementary techniques. We found no significant differences in the size, quantity, or morphology between exosomes secreted by HCFs with or without KC. Expression of CD81 was confirmed by immuno-EM, and expression of CD63 and CD9 with western blots in all exosome samples. We detected the expression of 72-144 miRNAs (threshold cycle Ct < 36) in all exosome samples. In HKC-derived exosome samples, miR-328-3p, miR-532-5p, miR-345-5p, and miR-424-5p showed unique expression, while let-7c-5p and miR-665 have increased expression. Protein profiling identified 157 proteins in at least half of the exosome samples, with 38 known exosomal proteins. We identified 12 up- and 2 down-regulated proteins in HKC-derived exosomes. The proteins are involved in membrane-bounded vesicles, cytoskeletal, calcium binding, and nucleotide binding. These proteins are predicted to be regulated by NRF2, miR-205, and TGF-ß1, which are involved in KC pathogenesis. We successfully characterized the HKC-derived exosomes and profiled their miRNA and protein contents, suggesting their potential role in KC development. Further studies are necessary to determine if and how these exosomes with differential protein/miRNA profiles contribute to the pathogenesis of KC.

Role of Fibroblast Growth Factor Receptor 2 (FGFR2) in Corneal Stromal Thinning.

Purpose: To determine the role of fibroblast growth factor receptor 2 (FGFR2)-mediated signaling in keratocytes during corneal development, a keratocyte-specific FGFR2-knockout (named FGFR2cKO) mouse model was generated, and its phenotypic characteristics were determined. Methods: A FGFR2cKO mouse model was generated by the following method: FGFR2 flox mice were crossed with the inducible keratocyte specific-Cre mice (Kera-rtTA/tet-O-Cre). Both male and female FGFR2cKO- and control mice (1 to 3-months-old) were analyzed for changes in corneal topography and pachymetry maps using the optical coherence tomography (OCT) method. The comparative TUNEL assay and immunohistochemical analyses were performed using corneas of FGFR2cKO and control mice to determine apoptotic cells, and expression of collagen-1 and fibronectin. Transmission electron microscopic analysis was conducted to determine collagen structures and their diameters in corneas of FGFR2cKO and control mice. Results: OCT-analyses of corneas of FGFR2cKO mice (n = 24) showed localized central thinning and an increased corneal steepness compared to control mice (n = 23). FGFR2cKO mice further showed a decreased expression in collagen-1, decreased collagen diameters, acute corneal hydrops, an increased fibronectin expression, and an increased number of TUNEL-positive cells suggesting altered collagen structures and keratocytes' apoptosis in the corneas of FGFR2cKO mice compared to control mice. Conclusions: The FGFR2cKO mice showed several corneal phenotypes (as described above in the results) that are also exhibited by the human keratoconus corneas. The results suggested that the FGFR2cKO mouse model serves to elucidate not only the yet unknown role of FGFR2-mediated signaling in corneal physiology but also serves as a model to determine molecular mechanism of human keratoconus development.

Screening and identification of genes related to ferroptosis in keratoconus.

Corneal keratoconus (KC) is a dilated (ectatic) corneal disease characterized by a central thinning of the cornea, which causes protrusion into a conical shape that seriously affects vision. However, due to the complex etiology of keratoconus, its entire mechanism remains unclear and there is no mechanism-directed treatment method. Ferroptosis is a novel programmed cell death mechanism related to lipid peroxidation, stress, and amino acid metabolism, which plays a crucial role in various diseases. This study aimed to explore the relationship between keratoconus and ferroptosis, to provide new insights into the mechanism of keratoconus development, and potential treatment options based on further elucidation of this mechanism. The corresponding mRNA microarray expression matrix data of KC patients were obtained from GEO database (GSE204791). Weighted co-expression network analysis (WGCNA) and support vector machine recursive feature elimination (SVM-RFE) were selected to screen hub genes, which were overlapped with ferroptosis genes (FRGs) from FerrDb. GO and GSEA were performed to analyze differential pathways, ssGSEA was used to determine immune status, and then, feasible drugs were predicted by gene-drug network. Additionally, we predicted the miRNA and IncRNA of hub genes to identify the underlying mechanism of disease so as to predict treatment for the disease. The epithelial transcriptome from keratoconus tissue mRNA microarray data (GSE204791) was extracted for the main analysis, including eight epithelial cells and eight epithelial control cells. The differential genes that were overlapped by WGCAN, SVM-RFE and FRGs were mainly related to oxidative stress, immune regulation, cellular inflammation, and metal ion transport. Through further analysis, aldo-keto reductase family 1 member C3 (AKR1C3) was selected, and negatively correlated with mature CD56 natural killer (NK) cells and macrophages. Then, gene-drug interaction network analysis and miRNA prediction were performed through the website. It was concluded that four immune-related drugs (INDOMETHACIN, DAUNORUBICIN, DOXORUBICIN, DOCETAXEL) and a miRNA (has-miR-184) were screened to predict potential drugs and targets for disease treatment. To our knowledge, this was the first report of KC being associated with ferroptosis and prompted search for differential genes to predict drug targets of gene immunotherapy. Our findings provided insight and a solid basis for the analysis and treatment of KC.