Autophagy and senescence of rat retinal precursor cells under high glucose

Abstract

Backgrounds:

Diabetic retinopathy (DR) is a common diabetic ocular disease characterized by retinal ganglion cell (RGC) changes. An abnormal environment, hyperglycemia, may progressively alter the structure and function of RGCs, which is a primary pathological feature of retinal neurodegeneration in DR. Accumulated studies confirmed autophagy and senescence play a vital role in DR; however, the underlying mechanisms need to be clarified.

Methods:

This study included the microarray expression profiling dataset GSE60436 from Gene Expression Omnibus (GEO) to conduct the bioinformatics analysis. The R software was used to identify autophagy-related genes (ARGs) that were differentially expressed in fibrovascular membranes (FVMs) and normal retinas. Co-expression and tissue-specific expression were elicited for the filtered genes. The genes were then analyzed by ontology (GO) enrichment analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis and Gene Set Enrichment Analysis (GSEA). R28 cells were cultured with high glucose, detected by reverse transcription-quantitative (RT-qPCR) and stained by apoptosis kit.

Results:

In the retina, 31 differentially expressed ARGs (24 up-regulated genes) were discovered and enriched. The enrichment results revealed that differentially expressed ARGs were significantly enriched in autophagy, apoptosis, aging, and neural function. Four hub genes (i.e., TP53, CASP1, CCL2, and CASP1) were significantly up-regulated. Upregulation of cellular autophagy and apoptosis level was detected in the hyperglycemia model in vitro.

Conclusions:

Our results provide evidence for the autophagy and cellular senescence mechanisms involved in retinal hyperglycemia injury, and the protective function of autophagy is limited. Further study may favour understanding the disease progression and neuroprotection of DR.

Introduction

Diabetic retinopathy (DR) is a common microvascular complication of diabetes mellitus (DM), one of the world’s fastest-growing diseases that has become a leading cause of blindness in the working-age population (1). By 2040, 642 million people are expected to suffer from DM (2). Patients will face physical and mental suffering at an advanced stage. Within five years, almost half of patients with untreated proliferative diabetic retinopathy (PDR), the most resistant form of DR, will go blind (3, 4). A series of pathological changes are proved in PDR, including oxidative stress, endoplasmic reticulum stress, chronic retinal inflammation (5), retinal angiogenesis, retinal neuron injury (6), and other retinal disorders.

Retinal ganglion cells (RGCs) are bridges between photoreceptors and the brain to transmit visual information with the help of bipolar cells. The morphological and functional integrity of RGCs plays a vital role in maintaining normal retinal function. RGCs have their bodies in different retinal layers, making RGCs vulnerable to the retinal microenvironment and other cells (e.g., Müller cells, glial cells, astrocytes, and microglia). Excessive glucose in the retinal microenvironment (i.e., hyperglycemia) and reactive oxygen species (ROS) accumulation can alter the structure and function of RGCs, causing damage to them that progresses and leads to DR (7). RGCs injury is the onset of retinal neurodegeneration in DR. In advanced DR with proliferative fibrovascular membranes, retinal vascular occlusion, non-perfusion, retinal ischemia, and hypoxia can be observed clinically. These pathological changes are highly suggestive of elevated autophagy in retinal tissues (8, 9).

Autophagy refers to any intracellular process involving the lysosome’s degradation of cytosolic components. There are three distinct autophagic pathways: macroautophagy (autophagy), microautophagy and Chaperone-mediated autophagy. Autophagy is essential for survival, differentiation, development, and homeostasis and is widely involved in ocular diseases, such as DR and age-related macular degeneration (10–16).

Previous studies mainly focused on pericytes and RPE cells’ autophagy in DR. Wang et al. (17) found that the high expression of the autophagy-related gene (ARG) MAPK3 in retinal pigment epithelial cells might be involved in outer blood-retina barrier destruction in DR. Lin et al. (18) described in detail the effect of oxidative stress-activated autophagy on RGCs; however, it remains controversial whether the autophagy in RGCs is beneficial. It is worth noting that there is increasing evidence that autophagy is related to senescence, another common pathological change in endocrine and neurological diseases (19–21). Conclusively, the role of the double-edged sword of autophagy in DR is complicated; particularly significant is that the role and mechanism of RGCs autophagy in hyperglycemia and DR remain unclear.

This study investigated the potential mechanism for autophagy in DR via in silico and in vitro approaches. Our results showed that R28 cells under high glucose conditions exhibited significant upregulation of some ARGs, accompanied by significant apoptosis and senescence, consistent with the results of bioinformatics analysis.

Materials and methods

Microarray data and autophagy-related genes collection

The gene expression dataset GSE60436 was acquired from the Gene Expression Omnibus (GEO) database (http://www.ncbi.nlm.nih.gov/geo/) (22). Six fibrovascular membranes (FVMs) samples from PDR patients and three normal retinal tissue samples from donated eyes were included in this dataset. The original research summarized the raw gene expression data using the Illumina Bead Studio software. The thresholding of the signal values was normalized to the 75th percentile, and the baseline was transformed to the median using Genespring GX 11.0 software. The datasets were subjected to Welch’s t-tests and multidimensional false-discovery control (the cutoff value between FVMs and retina was 0.5%).

In this study, the expression data were normalized using the ‘limma’ package in R software (4.2.0 GUI). The 222 ARGs were obtained from The Human Autophagy Database (http://www.autophagy.lu/index.html), which provides a comprehensive and up-to-date list of human genes and proteins involved directly or indirectly in autophagy as described in the literature.

Differential expression analysis of autophagy-related genes

Using the ‘limma’ package to screen the DEGs, P value 0.05 and |log2FoldChange| > 1 filters were used to narrow the results. R software was used to visualize the DEGs and autophagy-related DEGs. The expression in the retina of differentially expressed ARGs was evaluated and visualized using BioGPS website (http://biogps.org) (23–25) and GraphpadPrism software (version 9.0.1 for macOS). Expression values from Affymetrix chips relate to fluorescence intensity, which was mostly summarized using various data processing algorithms, gcrma, in BioGPS.

Construction of co-expression network

The Search Tool for the Retrieval of Interacting Genes/Proteins (STRING v.11.5, https://string-db.org/) database and Cytoscape (v.3.8.0) were used to establish and perform the visualization of the protein-protein interaction (PPI) network. Ten hub genes were identified by the cytoHubba plugin in Cytoscape with the Maximal Clique Centrality (MCC) method. The interactions between the differentially expressed ARGs were analyzed by GeneMANIA (http://genemania.org/) (26), a gene function prediction tool for the genetic interaction (GI) network. Additionally, the bulk tissue expression of the top ten hub genes was analyzed on GTEx Portal (https://www.gtexportal.org/) (27, 28).

Gene ontology and kyoto encyclopedia of genes and genomes pathway enrichment

Based on the differentially expressed ARGs, gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment were performed using the Database for Annotation, Visualization and Integrated Discovery (DAVID) v2022q2 (https://david.ncifcrf.gov/) (29, 30). The GO analysis contained biological processes (BPs), cellular components (CCs), and molecular functions (MFs).

Gene set enrichment analysis for pathways

Gene Set Enrichment Analysis (GSEA) was used to investigate the probable biological functions and pathways involved in DR in FVMs and normal retinas. GSEA was carried out using WebGestalt (WEB-based Gene SeT AnaLysis Toolkit, http://www.webgestalt.org/) (31–34). The top ten significant terms were shown.

Cell line and culture

Rat retinal precursor R28 cells were cultured in Dulbeccos Modified Eagle Medium (DMEM, BasalMedia, China) containing 10% fetal bovine serum (FBS, Gibco, USA) and 1% Antibiotic-Antimycotic (Gibco, USA) at 37°C with 5% CO₂. Cells were rinsed with PBS (Gibco, USA) and treated with 0.05% trypsin (Gibco, USA) for passage at 1: 3.

In subsequent experiments, the logarithmic growth phase cells were fed with a medium containing either 30 mM D-glucose (HG group) or a medium containing 5.6 mM D-glucose (LG group) for 48 hours, as appropriate.

RNA extraction and reverse transcription-quantitative PCR

Cells were placed in TRIzol (TaKaRa Biotechnology, China) reagent for RNA isolation. Chloroform (Macklin Reagent, Beijing, China) was applied to the cells for 15 seconds before being incubated at room temperature for 10 minutes. The top layer was transferred to a fresh microcentrifuge tube and mixed with isopropanol after being centrifuged at 10,000 × g for 10 minutes at 4°C. The material was centrifuged, rinsed, and then RNase-free water was used to elute the total RNA. Using a NanoVue Plus spectrophotometer (GE Healthcare, USA), the size and purity of the RNA were evaluated for quality. The extracted total RNA was treated to reverse transcription using a RevertAid Master Mix kit to convert isolated RNA into cDNA (Thermo Scientific, USA).

On a StepOne Plus Real-Time PCR System (Applied Biosystems, USA), cDNA from the cells was utilized as the template for reverse transcription-quantitative PCR (RT-qPCR) amplification using TB Green Premix Ex TaqII (TaKaRa Biotechnology, China). In separate experiments, RT-qPCR was performed to assess the transcript levels of TP53, CASP3, CCL2, HIF1A, CASP1, CDKN1A, CDKN2A, and β-actin in cell samples using the gene-specific primers shown in Supplementary Table S1. Target gene expression levels were calculated by 2^-△△ cycle threshold (Ct) method and normalized to that of β-actin.

Flow cytometry

As reported previously, apoptosis was evaluated by flow cytometry (35). Briefly, R28 cells were harvested with EDTA-free trypsin (Gibco, USA) and centrifuged at 300 × g for 5 minutes at 4 °C. The cells were washed with pre-cooled PBS twice and then stained with Annexin V-FITC Apoptosis Detection Kit (Vazyme, China). Each sample was resuspended with 100 μL Binding Buffer and stained with 5 μL Annexin V-FITC and 5 μL PI Staining Solution at dark at room temperature for 10 minutes. After adding 400 μL Binding Buffer, apoptotic cells were detected using a FACS Dxp AthenaTM (Cytek, Fremont, USA) and analyzed using FlowJo v10.0 (FlowJo Software, USA).

Statistical analysis

Statistical analyses were carried out by GraphpadPrism software. The statistical significance for RT-qPCR and flow cytometry were compared using Two-way and One-way ANOVA, respectively. The difference was statistically significant when P < 0.05.

Results

Data collection and differentially expressed autophagy-related genes identification

As shown in Figure 1, we firstly downloaded the expression profiling by array dataset GSE60436 (GPL6884, Illumina HumanWG-6 v3.0 expression BeadChip) from the GEO database. This dataset included six FVMs (three active FVMs: age 57.3 ± 2.1 years, one man and two women, duration of diabetes 15.7 ± 2.1 years, glycosylated hemoglobin 8.3%; three inactive FVMs: age 48.0 ± 13.0 years, one man and two women, duration of diabetes 12.3 ± 2.5 years, glycosylated hemoglobin 6.1%) and three normal retinas as control. After normalizing the expression data with R (Figure 2A), we found 1,442 up-regulated genes and 1,618 down-regulated genes in FVMs compared to control samples (Figure 2B). Among these DEGs, 31 ARGs were screened (Figure 2C), seven were down-regulated, and 24 were up-regulated (Table 1). The down-regulated autophagy-related DEGs were GABARAPL1, PRKAR1A, HSPB8, CTSD, MAPK8IP1, ITPR1, and CX3CL1. The top ten DEGs with the most significant upregulation were CDKN2A, CDKN1A, ITGB1, DRAM1, CASP1, CCL2, CASP4, BIRC5, CXCR4, and CASP3. To gain further insight into the tissue distribution of genes of interest, we searched the BioGPS database. The results indicated that the autophagy-related DEGs were enriched in human retina tissue, particularly CCL2, CDKN1A, CTSD, HIF1A, HSPB8 and other genes (Figure 3).

Figure 1

Figure 2

Figure 3

Co-expression network construction and hub genes identification

We used data from STRING and GeneMANIA to establish PPI and GI networks for investigating the interactions at the protein and gene levels. We used the MCC method to obtain the hub genes based on all the autophagy-related DEGs. These genes have extensive interactions at the protein (Figures 4A, B) and gene levels (Figure 4C). The top ten hub genes with descending rank were TP53, CASP3, CCL2, HIF1A, CASP1, FAS, CXCR4, TNFSF10, CDKN2A, and ITGB1 (Figure 4D). All these ten genes were significantly up-regulated in FVMs, which indicated autophagy might be overactivated in DR (Figure 5). To deeper understand the bulk tissue expression of the top ten hub genes, the GTEx Portal was accessed (Supplementary Figure S1).

Figure 4

Figure 5

Functional and pathway enrichment

The GO, KEGG and GSEA enrichment was carried out to study the functions, subcellular localizations and pathways of the autophagy-related DEGs. For BPs, the DEGs were mainly enriched at autophagy, apoptosis, mitophagy, immune response, and other terms (Figure 6A). For CCs, in addition to inflammasomes and autophagosomes, DEGs were enriched in the neuronal cell body, axonal growth cone, nuclear body, and cellular membranes (Figure 6B). MFs revealed the DEGs involved in the death receptor binding and binding of p53 and other proteins and enzymes (Figure 6C). Furthermore, autophagy, apoptosis, cellular senescence, necroptosis, and neurodegeneration pathway were popular KEGG terms for DEGs (Figure 6D). The autophagy, apoptosis, cellular senescence, and p53 signaling were also confirmed by GSEA pathway enrichment (Supplementary Figure S2). To further validate the apoptosis level, we performed flow cytometry. The results showed that the apoptosis level of R28 cells was significantly increased after high glucose challenge (Supplementary Figure S3). Compared with the control group, the levels of late apoptosis and total apoptosis were significantly increased in the high glucose stimulation group (P < 0.01). Although not statistically different, the HG also had higher levels of early apoptosis than the LG.

Figure 6

Moreover, the down-regulated genes were primarily enriched in neuron protection, cell adhesion, and metabolism terms (Supplementary Table S2). In this study, most autophagy-related DEGs were up-regulated and associated with angiogenesis, inflammation, fibrosis, senescence, and programmed cell death (Supplementary Table S3). Cellular senescence and cell death (i.e., autophagy, apoptosis, and necroptosis) were speculated to get involved in DR occurrence and development. It is worth noting that some neuronal terms were enriched, which suggested that retinal nerve damage is also one of the main pathological changes of DR.

Differentially expressed autophagy-related genes expression validation

It was found that RGCs autophagy is highly involved in DR development. Our results mentioned above suggested autophagy was up-regulated in DR samples. In this present study, we incubated R28 cells with high glucose (30 mM) for 48 hours to mimic the retinal hyperglycemic environment and measured hub gene expression by RT-qPCR. The results were generally consistent with expectations that TP53, CASP3, CCL2, and CASP1 mRNA were significantly higher in the HG group than in the LG group (Figure 7). However, HIF1A remains unchanged statistically, although a slight uptrend was observed. In addition to the top five hub genes, we detected two senescence-associated genes, CDKN1A and CDKN2A, in expression. Interestingly, CDKN2A, rather than CDKN1A, was significantly up-regulated in the HG group (Figure 7).

Figure 7

Discussion

In this study, we investigated the RGCs autophagy in DR. DR is a complex and gradual chronic complication of diabetic patients. Previous studies have shown several metabolic pathways involved in DR pathological changes. It is worth noting that these pathways’ abnormal signaling may lead to microvascular complications (e.g., blood-retinal barrier destruction, microvascular leakage, and edema) and neuron disorders (e.g., neurodegeneration and neuroinflammation). The neurovascular unit (NVU), which is made up of neurons, vascular tissue, and supporting glial cells, is found in the retina (36). The ocular NVU may be affected if the blood-retinal barrier (BRB) breaks down. DR has traditionally been associated with retinal angiopathy or vasculopathy, however there is evidence of neurological impairments (37). As a result, diabetic retinal neurodegeneration (DRN) and diabetic retinal vasculopathy (DRV) are currently thought to be two connected parts of DR. The BRB breakdown and retinal vascular disorders in DR are widely studied. Notably, Wang et al. (17) have primarily demonstrated the role of RPE autophagy in BRB damage in DR. DRN is now considered an initial event contributing to DRV and is detectable before the clinical appearance of DRV (38). In-depth exploration of DRN is more helpful to directly understand the mechanism and prevention methods of blindness caused by DR.

In this study, we focused on RGCs loss, one of the most significant features of DRN. We screened a GEO dataset with an autophagy gene list from The Human Autophagy Database and found that 31 ARGs were significantly differentially expressed. The top ten hub genes were screened out via the MCC method (TP53, CASP3, CCL2, HIF1A, CASP1, FAS, CXCR4, TNFSF10, CDKN2A, and ITGB1). Additionally, these genes were highly enriched in retinal tissue, which suggested that retinal autophagy levels might be altered in DR. Further analyses showed that the up-regulated ARGs were associated with autophagy, apoptosis, cellular senescence, necroptosis, and neurodegeneration. Our results were consistent with clinical findings that long-term and chronic hyperglycemia may cause a significant and extensive reduction of RGCs and thinning of the RGCs layer in DM (39). These preliminary results confirmed the theory of nerve damage in DR.

To further verify the genes of interest expression in RGCs, we construct an in vitro high glucose model using R28 cells to mimic the retina environment of DR. The RT-qPCR results of the top five hub genes and two senescence-related genes (i.e., CDKN1A and CDKN2A) showed that TP53, CASP3, CCL2, CASP1, and CDKN2A were significantly up-regulated in HG than LG (Figure 7). In response to various cellular stressors, the tumor protein p53 (TP53) controls the expression of target genes, causing cell cycle arrest, apoptosis, autophagy, senescence, DNA repair, and adjustments in cellular metabolism (40). Previous studies found higher TP53 expression in diabetic retinal tissue than in control samples and considered TP53 polymorphism a genetic marker for diabetic complications (41–43). The caspase family has been proven to play a vital role in DR, particularly CASP3 and CASP1. CASP3 and CASP9 showed increased expression in RGCs, suggesting early apoptosis of RGCs (44). In addition to RGCs, animal studies showed that increased CASP3 could reduce the rod outer segments length, increase early apoptosis in photoreceptors in DM, and increase apoptosis of dopaminergic amacrine cells in the retina with hyperglycemia (45–47). The direct activation of CASP1 participated in signaling transformation in the canonical pathway of pyroptosis. Wang et al. (48) have confirmed that pyroptosis and inflammation improve DR’s occurrence and development. CCL2 gets involved in neuroinflammation. Treatment with CCL2 could promote macrophage-like cells recruited on the retina surface, a characteristic of vision-threatening retinal vascular diseases like DR and retinal vein occlusion (49). The regulation between HIF1A and TP53 is bidirectional and complex, and one of the possible reasons there was no statistical difference in HIF1A expression is that it was suppressed by highly expressed TP53 (50). The CDKN2A was highly expressed in HG rather than LG; however, no significant difference was observed for CDKN1A. This result indicated that the R28 cells in HG suffered cellular senescence, and cells have entered the senescence maintenance stage under a long-term high-glucose environment.

Apoptosis is another meaningful programmed cell death and shares some regulatory genes with autophagy. Our results showed a significant increase in apoptotic cells in the HG group. Appropriate apoptosis and autophagy are protective mechanisms of cells to scavenge unhealthy cells and metabolic wastes to avoid carcinogenesis and aging. Interestingly, as we showed, excessive activation of apoptosis and autophagy under long-term chronic hyperglycemia did not alleviate aging, but further aggravated RGCs loss. Additionally, the period of autophagy’s protective effect might be shortened with acute injury (51). Although inhibition of autophagy decreased the apoptosis of RGCs in diabetic retinas (9), in high glucose, the self-rescue effect of RGCs via autophagy is limited. Therefore, harnessing the protective effects of apoptosis and autophagy in short- and long-term hyperglycemic environments in DR is a complex topic that requires more exploration.

We acknowledge the limitations that this study used RGCs to verify the differences between FVMs and normal retinal tissue. Due to the difficulty of obtaining human retinal tissue, there are few sequencing results on the retina and RGCs. By mining the differentially expressed genes in FVMs and normal retinas, we obtained genes potentially altered in retinal under a high glucose environment and inferred the damage targets of high glucose to RGCs. Like the studies by Wang et al. (17, 48), we extrapolated the conclusions of the initial results to provide a new perspective on the pathogenesis of DR. This study preliminarily confirmed that autophagy and aging of RGCs are involved in high glucose injury combined with sequencing results and gene expression verification. Although R28 cells display both glial and neuronal cell markers, they are widely used in studies. They can respond to various stimuli with advantages including unlimited lifespan, not being tumorigenic in vivo, and closely simulating a retinal explant (52).

Conclusion

In summary, we investigated the potential mechanism of autophagy and senescence in DR via bioinformatics approaches and validation in vivo. We found that under high glucose conditions, the levels of some autophagy- and senescence-related genes in R28 cells were significantly increased. These results initially revealed that autophagy and senescence-related genes (e.g., TP53 and CDKN2A) regulated retinal cells involved in DR pathogenesis and challenged the protective effect of enhanced autophagy in retinal cells under high glucose environment. Further studies could provide novel mechanistic insight into DR.

References

• 1

  TeoZLThamYCYuMCheeMLRimTHCheungNet al. Global prevalence of diabetic retinopathy and projection of burden through 2045: Systematic review and meta-analysis. Ophthalmology (2021) 128(11):1580–91. doi: 10.1016/j.ophtha.2021.04.027

  • CrossRef
  • Google Scholar

• 2

  OgurtsovaKda Rocha FernandesJDHuangYLinnenkampUGuariguataLChoNHet al. IDF diabetes atlas: Global estimates for the prevalence of diabetes for 2015 and 2040. Diabetes Res Clin Pract (2017) 128:40–50. doi: 10.1016/j.diabres.2017.03.024

  • CrossRef
  • Google Scholar

• 3

  FerrisFL3rd. Results of 20 years of research on the treatment of diabetic retinopathy. Prev Med (1994) 23(5):740–2. doi: 10.1006/pmed.1994.1127

  • CrossRef
  • Google Scholar

• 4

  AntonettiDAKleinRGardnerTW. Diabetic retinopathy. N Engl J Med (2012) 366(13):1227–39. doi: 10.1056/NEJMra1005073

  • CrossRef
  • Google Scholar

• 5

  TangJKernTS. Inflammation in diabetic retinopathy. Prog Retin Eye Res (2011) 30(5):343–58. doi: 10.1016/j.preteyeres.2011.05.002

  • CrossRef
  • Google Scholar

• 6

  DuhEJSunJKStittAW. Diabetic retinopathy: current understanding, mechanisms, and treatment strategies. JCI Insight (2017) 2(14). doi: 10.1172/jci.insight.93751

  • CrossRef
  • Google Scholar

• 7

  NgDSChiangPPTanGCheungCGChengCYCheungCYet al. Retinal ganglion cell neuronal damage in diabetes and diabetic retinopathy. Clin Exp Ophthalmol (2016) 44(4):243–50. doi: 10.1111/ceo.12724

  • CrossRef
  • Google Scholar

• 8

  TangLXuGTZhangJF. Inflammation in diabetic retinopathy: possible roles in pathogenesis and potential implications for therapy. Neural Regener Res (2023) 18(5):976–82. doi: 10.4103/1673-5374.355743

  • CrossRef
  • Google Scholar

• 9

  ParkHLKimJHParkCK. Different contributions of autophagy to retinal ganglion cell death in the diabetic and glaucomatous retinas. Sci Rep (2018) 8(1):13321. doi: 10.1038/s41598-018-30165-7

  • CrossRef
  • Google Scholar

• 10

  KangQYangC. Oxidative stress and diabetic retinopathy: Molecular mechanisms, pathogenetic role and therapeutic implications. Redox Biol (2020) 37:101799. doi: 10.1016/j.redox.2020.101799

  • CrossRef
  • Google Scholar

• 11

  ChaiPNiHZhangHFanX. The evolving functions of autophagy in ocular health: A double-edged sword. Int J Biol Sci (2016) 12(11):1332–40. doi: 10.7150/ijbs.16245

  • CrossRef
  • Google Scholar

• 12

  ChenCLChenYHLiangCMTaiMCLuDWChenJT. Glucosamine-induced autophagy through AMPK(-)mTOR pathway attenuates lipofuscin-like autofluorescence in human retinal pigment epithelial cells In vitro. Int J Mol Sci (2018) 19(5). doi: 10.3390/ijms19051416

  • CrossRef
  • Google Scholar

• 13

  KasettiRBMaddineniPKiehlbauchCPatilSSearbyCCLevineBet al. Autophagy stimulation reduces ocular hypertension in a murine glaucoma model via autophagic degradation of mutant myocilin. JCI Insight (2021) 6(5). doi: 10.1172/jci.insight.143359

  • CrossRef
  • Google Scholar

• 14

  IshikawaMTakasekiSYoshitomiTCoveyDFZorumskiCFIzumiY. The neurosteroid allopregnanolone protects retinal neurons by effects on autophagy and GABRs/GABAA receptors in rat glaucoma models. Autophagy (2021) 17(3):743–60. doi: 10.1080/15548627.2020.1731270

  • CrossRef
  • Google Scholar

• 15

  RoySKernTSSongBStuebeC. Mechanistic insights into pathological changes in the diabetic retina: Implications for targeting diabetic retinopathy. Am J Pathol (2017) 187(1):9–19. doi: 10.1016/j.ajpath.2016.08.022

  • CrossRef
  • Google Scholar

• 16

  RosaMDDistefanoGGaglianoCRuscianoDMalaguarneraL. Autophagy in diabetic retinopathy. Curr Neuropharmacol (2016) 14(8):810–25. doi: 10.2174/1570159x14666160321122900

  • CrossRef
  • Google Scholar

• 17

  WangNWeiLLiuDZhangQXiaXDingLet al. Identification and validation of autophagy-related genes in diabetic retinopathy. Front Endocrinol (Lausanne) (2022) 13:867600. doi: 10.3389/fendo.2022.867600

  • CrossRef
  • Google Scholar

• 18

  LinWJKuangHY. Oxidative stress induces autophagy in response to multiple noxious stimuli in retinal ganglion cells. Autophagy (2014) 10(10):1692–701. doi: 10.4161/auto.36076

  • CrossRef
  • Google Scholar

• 19

  ZhangCYTanXHYangHHJinLHongJRZhouYet al. COX-2/sEH dual inhibitor alleviates hepatocyte senescence in NAFLD mice by restoring autophagy through Sirt1/PI3K/AKT/mTOR. Int J Mol Sci (2022) 23(15). doi: 10.3390/ijms23158267

  • CrossRef
  • Google Scholar

• 20

  ZagkouSMaraisVZeghoudiNGuillouELEskelinenELPanasyukGet al. Design and evaluation of autophagy-inducing particles for the treatment of abnormal lipid accumulation. Pharmaceutics (2022) 14(7). doi: 10.3390/pharmaceutics14071379

  • CrossRef
  • Google Scholar

• 21

  GonzalezCDResnikRVaccaroMI. Secretory autophagy and its relevance in metabolic and degenerative disease. Front Endocrinol (Lausanne) (2020) 11:266. doi: 10.3389/fendo.2020.00266

  • CrossRef
  • Google Scholar

• 22

  IshikawaKYoshidaSKobayashiYZhouYNakamaTNakaoSet al. Microarray analysis of gene expression in fibrovascular membranes excised from patients with proliferative diabetic retinopathy. Invest Ophthalmol Vis Sci (2015) 56(2):932–46. doi: 10.1167/iovs.14-15589

  • CrossRef
  • Google Scholar

• 23

  WuCJinXTsuengGAfrasiabiCSuAI. BioGPS: Building your own mash-up of gene annotations and expression profiles. Nucleic Acids Res (2016) 44(D1):D313–316. doi: 10.1093/nar/gkv1104

  • CrossRef
  • Google Scholar

• 24

  WuCMacleodISuAI. BioGPS and MyGene.info: organizing online, gene-centric information. Nucleic Acids Res (2013) 41(Database issue):D561–565. doi: 10.1093/nar/gks1114

  • CrossRef
  • Google Scholar

• 25

  WuCOrozcoCBoyerJLegliseMGoodaleJBatalovSet al. BioGPS: An extensible and customizable portal for querying and organizing gene annotation resources. Genome Biol (2009) 10(11):R130. doi: 10.1186/gb-2009-10-11-r130

  • CrossRef
  • Google Scholar

• 26

  Warde-FarleyDDonaldsonSLComesOZuberiKBadrawiRChaoPet al. The GeneMANIA prediction server: biological network integration for gene prioritization and predicting gene function. Nucleic Acids Res (2010) 38(Web Server issue):W214–220. doi: 10.1093/nar/gkq537

  • CrossRef
  • Google Scholar

• 27

  ConsortiumGT. The genotype-tissue expression (GTEx) project. Nat Genet (2013) 45(6):580–5. doi: 10.1038/ng.2653

  • CrossRef
  • Google Scholar

• 28

  ConsortiumGT. The GTEx consortium atlas of genetic regulatory effects across human tissues. Science (2020) 369(6509):1318–30. doi: 10.1126/science.aaz1776

  • CrossRef
  • Google Scholar

• 29

  DennisGJr.ShermanBTHosackDAYangJGaoWLaneHCet al. DAVID: Database for annotation, visualization, and integrated discovery. Genome Biol (2003) 4(5):P3.

  • Google Scholar

• 30

  ShermanBTHaoMQiuJJiaoXBaselerMWLaneHCet al. DAVID: A web server for functional enrichment analysis and functional annotation of gene lists (2021 update). Nucleic Acids Res (2022). doi: 10.1093/nar/gkac194

  • CrossRef
  • Google Scholar

• 31

  ZhangBKirovSSnoddyJ. WebGestalt: An integrated system for exploring gene sets in various biological contexts. Nucleic Acids Res (2005) 33(Web Server issue):W741–748. doi: 10.1093/nar/gki475

  • CrossRef
  • Google Scholar

• 32

  WangJDuncanDShiZZhangB. WEB-based GEne SeT AnaLysis toolkit (WebGestalt): update 2013. Nucleic Acids Res (2013) 41(Web Server issue):W77–83. doi: 10.1093/nar/gkt439

  • CrossRef
  • Google Scholar

• 33

  WangJVasaikarSShiZGreerMZhangB. WebGestalt 2017: a more comprehensive, powerful, flexible and interactive gene set enrichment analysis toolkit. Nucleic Acids Res (2017) 45(W1):W130–7. doi: 10.1093/nar/gkx356

  • CrossRef
  • Google Scholar

• 34

  LiaoYWangJJaehnigEJShiZZhangB. WebGestalt 2019: Gene set analysis toolkit with revamped UIs and APIs. Nucleic Acids Res (2019) 47(W1):W199–205. doi: 10.1093/nar/gkz401

  • CrossRef
  • Google Scholar

• 35

  MaCZhuangZSuQHeJLiH. Curcumin has anti-proliferative and pro-apoptotic effects on tongue cancer in vitro: A study with bioinformatics analysis and in vitro experiments. Drug Des Devel Ther (2020) 14:509–18. doi: 10.2147/dddt.S237830

  • CrossRef
  • Google Scholar

• 36

  AntonettiDALiethEBarberAJGardnerTW. Molecular mechanisms of vascular permeability in diabetic retinopathy. Semin Ophthalmol (1999) 14(4):240–8. doi: 10.3109/08820539909069543

  • CrossRef
  • Google Scholar

• 37

  SoniDSagarPTakkarB. Diabetic retinal neurodegeneration as a form of diabetic retinopathy. Int Ophthalmol (2021) 41(9):3223–48. doi: 10.1007/s10792-021-01864-4

  • CrossRef
  • Google Scholar

• 38

  CarrascoEHernandezCMirallesAHuguetPFarresJSimoR. Lower somatostatin expression is an early event in diabetic retinopathy and is associated with retinal neurodegeneration. Diabetes Care (2007) 30(11):2902–8. doi: 10.2337/dc07-0332

  • CrossRef
  • Google Scholar

• 39

  TakahashiHGotoTShojiTTanitoMParkMChiharaE. Diabetes-associated retinal nerve fiber damage evaluated with scanning laser polarimetry. Am J Ophthalmol (2006) 142(1):88–94. doi: 10.1016/j.ajo.2006.02.016

  • CrossRef
  • Google Scholar

• 40

  YinYStephenCWLucianiMGFahraeusR. p53 stability and activity is regulated by Mdm2-mediated induction of alternative p53 translation products. Nat Cell Biol (2002) 4(6):462–7. doi: 10.1038/ncb801

  • CrossRef
  • Google Scholar

• 41

  MohammadHMFSamiMMMakarySToraihEAMohamedAOEl-GhaieshSH. Neuroprotective effect of levetiracetam in mouse diabetic retinopathy: Effect on glucose transporter-1 and GAP43 expression. Life Sci (2019) 232:116588. doi: 10.1016/j.lfs.2019.116588

  • CrossRef
  • Google Scholar

• 42

  SlominskiBSkrzypkowskaMRyba-StanislawowskaMMysliwiecMTrzonkowskiP. Associations of TP53 codon 72 polymorphism with complications and comorbidities in patients with type 1 diabetes. J Mol Med (Berl) (2021) 99(5):675–83. doi: 10.1007/s00109-020-02035-1

  • CrossRef
  • Google Scholar

• 43

  GuoDFangLYuXWangCWangYGuoW. Different roles of TP53 codon 72 polymorphism in type 2 diabetes and its complications: Evidence from a case-control study on a Chinese han population. Int J Gen Med (2021) 14:4259–68. doi: 10.2147/IJGM.S322840

  • CrossRef
  • Google Scholar

• 44

  OshitariTYamamotoSHataNRoyS. Mitochondria- and caspase-dependent cell death pathway involved in neuronal degeneration in diabetic retinopathy. Br J Ophthalmol (2008) 92(4):552–6. doi: 10.1136/bjo.2007.132308

  • CrossRef
  • Google Scholar

• 45

  ChangMLWuCHJiang-ShiehYFShiehJYWenCY. Reactive changes of retinal astrocytes and Muller glial cells in kainate-induced neuroexcitotoxicity. J Anat (2007) 210(1):54–65. doi: 10.1111/j.1469-7580.2006.00671.x

  • CrossRef
  • Google Scholar

• 46

  KumarBGuptaSKSrinivasanBPNagTCSrivastavaSSaxenaRet al. Hesperetin rescues retinal oxidative stress, neuroinflammation and apoptosis in diabetic rats. Microvasc Res (2013) 87:65–74. doi: 10.1016/j.mvr.2013.01.002

  • CrossRef
  • Google Scholar

• 47

  GastingerMJSinghRSBarberAJ. Loss of cholinergic and dopaminergic amacrine cells in streptozotocin-diabetic rat and Ins2Akita-diabetic mouse retinas. Invest Ophthalmol Vis Sci (2006) 47(7):3143–50. doi: 10.1167/iovs.05-1376

  • CrossRef
  • Google Scholar

• 48

  WangNDingLLiuDZhangQZhengGXiaXet al. Molecular investigation of candidate genes for pyroptosis-induced inflammation in diabetic retinopathy. Front Endocrinol (Lausanne) (2022) 13:918605. doi: 10.3389/fendo.2022.918605

  • CrossRef
  • Google Scholar

• 49

  RajeshADrohoSLavineJA. Macrophages in close proximity to the vitreoretinal interface are potential biomarkers of inflammation during retinal vascular disease. J Neuroinflamm (2022) 19(1):203. doi: 10.1186/s12974-022-02562-3

  • CrossRef
  • Google Scholar

• 50

  ZhangCLiuJWangJZhangTXuDHuWet al. The interplay between tumor suppressor p53 and hypoxia signaling pathways in cancer. Front Cell Dev Biol (2021) 9:648808. doi: 10.3389/fcell.2021.648808

  • CrossRef
  • Google Scholar

• 51

  YangSZhangJChenL. The cells involved in the pathological process of diabetic retinopathy. BioMed Pharmacother (2020) 132:110818. doi: 10.1016/j.biopha.2020.110818

  • CrossRef
  • Google Scholar

• 52

  SeigelGM. Review: R28 retinal precursor cells: The first 20 years. Mol Vis (2014) 20:301–6.

  • Google Scholar