Identification of miRNA-mRNA Regulatory Networks Associated with Diabetic Retinopathy using Bioinformatics Analysis

Introduction: Diabetic retinopathy (DR) is a major complication of diabetes and a leading cause of visual loss. This study aimed to explore biomarkers for DR that may provide additional reference to DR pathogenesis and development. Methods: The differentially expressed genes (DEGs) between the DR and control samples in the GSE53257 dataset were identified. Logistics analyses were performed to identify DR-associated miRNAs and genes, and correlation analysis was performed to determine the correlation between them in GSE160306. Results: A total of 114 DEGs in DR were identified in GSE53257. Three genes, including ATP5A1 (down), DAUFV2 (down), and OXA1L (down), were differentially expressed between DR and control samples in GSE160306. Univariate logistics analysis identified that ATP5A1 (OR=0.007, p = 1.40E-02), NDUFV2 (OR = 0.003, p = 6.40E-03), and OXA1L (OR = 0.093, p = 3.08E-02) were DR-associated genes. ATP5A1 and OXA1L were regulated by multiple miRNAs, of which hsa-let-7b-5p (OR = 26.071, p = 4.40E-03) and hsa-miR-31-5p (OR = 4.188, p = 5.09E-02) were related to DR. ATP5A1 and OXA1L were closely correlated with each other in DR. Conclusion: The hsa-miR-31-5p-ATP5A1 and hsa-let-7b-5p-OXA1L axes might play novel and important roles in the pathogenesis and development of DR.