Comprehensive Bioinformatics Analysis to Reveal Key RNA Targets and Hub Competitive Endogenous RNA Network of Keratoconus

Abstract

Keratoconus (KC) is the most common corneal ectatic disease, with its pathological mechanisms unclear. We mainly performed bioinformatics approaches to reveal core RNA targets and hub competitive endogenous RNA (ceRNA) network and explored the potential regulatory mechanisms of ceRNA in KC. The high-throughput sequencing datasets GSE77938 and GSE151631 were downloaded from the Gene Expression Omnibus (GEO) database. The differential expression of mRNAs and lncRNAs was identified using the DESeq2 package. Functional enrichment analyses and protein–protein interaction (PPI) were executed. Then, the hub genes were filtered and molecular docking analysis was performed. Moreover, we predicted miRNAs through a website database and validated them using quantitative PCR (qPCR). Eventually, the lncRNA–miRNA–mRNA regulatory network was constructed by Cytoscape. We revealed that 428 intersected differentially expressed mRNA (DEGs) and 68 intersected differentially expressed lncRNA (DELs) were shared between the two datasets. Functional enrichment results innovatively showed that the ubiquitin-dependent protein catabolic process was upregulated in KC. The pathway enrichment showed that DEGs were mainly involved in NF-kB signaling and neurodegenerative diseases. In addition, we uncovered the top 20 hub genes in which FBXW11, FBXO9, RCHY1, and CD36 were validated by qPCR. Particularly, a small-molecule drug triptolide was predicted by molecular docking to be a candidate drug for treating KC. Moreover, we innovatively predicted and validated four core miRNAs (miR-4257, miR-4494, miR-4263, and miR-4298) and constructed a ceRNA network that contained 165 mRNA, eight lncRNAs, and four core miRNAs. Finally, we proposed a potential regulatory mechanism for KC. Overall, we uncovered a hub ceRNA network that might underlie a critical posttranslational regulatory mechanism in KC, in which miR-4257, miR-4494, miR-4263, and miR-4298 could be valuable biomarkers and provided core RNAs therapeutic targets for KC.