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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Inflammatory profile of keratoconic corneal epithelium.

BACKGROUND: Recent studies have presented inflammatory features on keratoconus (KC) and many inflammatory markers are described in the tears of patients with this disease. The KC pathogenesis is still unknown just like the correlation with inflammatory patterns. However, environmental and genetic issues may be part of the progress of KC. In addition, some systemic features, such as allergy and obesity, seem to be related to the progression of KC. Our purpose was to evaluate the neuropeptides vasoactive intestinal peptide (VIP), neuropeptide Y (NPY), chemokines ligand 2 (CCL-2) and 5 (CCL-5), and interleukins 6 (IL-6) and 8 (IL-8) on corneal epithelial cells and blood of patients with KC and in healthy controls. In addition, the neutrophil-to-lymphocyte ratio (NLR) was evaluated to predict inflammation. METHODS: This including prospective observational study included 32 KC patients who underwent corneal crosslinking (CXL) and 32 control patients who underwent photorefractive keratectomy (PRK). Patients' corneal epithelial cells were removed surgically, and blood (buffy coat) was analyzed. Samples in triplicate were evaluated on rt-PCR for neuropeptides (VIP e NPY), interleukins (IL-6 e IL-8), and chemokines (CCL-2 and CCL-5). RESULTS: Our study showed statistically higher CCL-5 and IL-8 on corneal epithelial cells in patients with KC. Blood cells were statistically higher in VIP and NPY in the KC group. Interleukin-8 on blood cells was statistically significant in KC'S group; for CCL-2 and CCL-5 they were statistically lower in patients with KC compared with controls. NLR showed no difference between the groups. CONCLUSIONS: Our data support the findings of other studies that suggested altering KC status, such as inflammatory corneal disease. The presence of IL-8 in the cornea and blood samples of KC's group suggested systemic disease with a possible local or repercussion action. Further studies are warranted to elucidate KC pathogenesis and its correlation to systemic disease.

Meta-Analysis of Keratoconus Transcriptomic Data Revealed Altered RNA Editing Levels Impacting Keratin Genomic Clusters.

Introduction: Keratoconus (KC) is an ocular disorder with a multifactorial origin. Transcriptomic analyses (RNA-seq) revealed deregulations of coding (mRNA) and non-coding RNAs (ncRNAs) in KC, suggesting that mRNA-ncRNA co-regulations can promote the onset of KC. The present study investigates the modulation of RNA editing mediated by the adenosine deaminase acting on dsRNA (ADAR) enzyme in KC. Materials: The level of ADAR-mediated RNA editing in KC and healthy corneas were determined by two indexes in two different sequencing datasets. REDIportal was used to localize known editing sites, whereas new putative sites were de novo identified in the most extended dataset only and their possible impact was evaluated. Western Blot analysis was used to measure the level of ADAR1 in the cornea from independent samples. Results: KC was characterized by a statistically significant lower RNA-editing level compared to controls, resulting in a lower editing frequency, and less edited bases. The distribution of the editing sites along the human genome showed considerable differences between groups, particularly relevant in the chromosome 12 regions encoding for Keratin type II cluster. A total of 32 recoding sites were characterized, 17 representing novel sites. JUP, KRT17, KRT76, and KRT79 were edited with higher frequencies in KC than in controls, whereas BLCAP, COG3, KRT1, KRT75, and RRNAD1 were less edited. Both gene expression and protein levels of ADAR1 appeared not regulated between diseased and controls. Conclusions: Our findings demonstrated an altered RNA-editing in KC possibly linked to the peculiar cellular conditions. The functional implications should be further investigated.

The association between VSX1 exon3 gene variants and keratoconus in Malaysian patients.

Purpose: This case-control study aims to examine possible associations of VSX1 exon3 gene variants with the development of keratoconus (KC) in Malaysian patients. Methods: A case-control study was done on 42 keratoconus cases, 127 family member controls, and 96 normal controls. Results: Three gene variants, p.A182A, p.P237P, and p.R217H showed significant associations with keratoconus (P < 0.05). While p.A182A and p.P227P were more prevalent than in the family and normal controls (OR 3.14-4.05), the reverse was observed with p.R217H (OR 0.086-1.59). With Haploview analysis, p.A182A and p.P237P were shown to be in linkage disequilibrium (LD) (LOD (logarithm of the odds score) score of 2.0, r2 of 0.957, and 95% confidence interval (CI) of 0.96-1.00). Conclusion: The study results suggest that the p.A182A and p.P237P variants could have contributed to the development of keratoconus in some Malaysians and that these two variants are likely to be co-inherited. In contrast, the p.R217H variant appeared to confer some protection against the development of keratoconus.

Investigation of Potential Crucial Genes and Key Pathways in Keratoconus: An Analysis of Gene Expression Omnibus Data.

Keratoconus is one of the most common causes leading to visual impairment in young adult population. The pathogenesis of keratoconus remains poorly understood. The aim of this study was to identify the potential key genes and pathways associated with keratoconus and to further analyze its molecular mechanism. Two RNA-sequencing datasets of keratoconus and paired normal corneal tissues from the Gene Expression Omnibus database were obtained. Differentially expressed genes (DEGs) were identified, and the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were conducted. The protein-protein interaction (PPI) network of the DEGs was established, and the hub genes and significant gene modules of PPI were further constructed. Lastly, the GO and KEGG analyses of the hub gene were performed. In total, 548 common DEGs were identified. GO enrichment analysis showed that the DEGs were primarily associated with regulation of cell adhesion, the response to molecule of bacterial origin, lipopolysaccharide and biotic stimulus, collagen-containing extracellular matrix, extracellular matrix, and structure organization. KEGG pathway analysis showed that these DEGs were mainly involved in the TNF signaling pathway, IL-17 signaling pathway, Rheumatoid arthritis, Cytokine-cytokine receptor interaction. The PPI network was constructed with 146 nodes and 276 edges, and 3 significant modules are selected. Finally, top 10 hub genes were identified from the PPI network. The results revealed that extracellular matrix remodeling and immune inflammatory response could be the key links of keratoconus, TNF, IL6, IL1A, IL1B, CCL3, MMP3, MMP9, MMP1, and TGFB1 may be potential crucial genes, and TNF signaling pathway and IL-17 signaling pathway were the potential pathways accounting for pathogenesis and development of keratoconus.