Skip to main content

Two ferroptosis-specific expressed genes NOX4 and PARP14 are considered as potential biomarkers for the diagnosis and treatment of diabetic retinopathy and atherosclerosis

Abstract

Objectives

Both Diabetic retinopathy (DR) and Atherosclerosis (AS) are common complications in patients with diabetes, and they share major pathophysiological similarities and have a common pathogenesis. Studies performed to date have demonstrated that ferroptosis plays a vital part in the occurrence and development of DR and AS, but its mechanism in the two diseases remains poorly understood.

Methods

DR Chip data (GSE60436 and GSE102485) and AS chip data (GSE100927 and GSE57691) were obtained from the Gene Expression Omnibus (GEO) database. The screening of the differential expression genes (DEGs) was analyzed using the limma package, and the genes related to ferroptosis were obtained from the FerrDb V2 database. Two key genes (NOX4 and PARP14) were identified through external datasets validation and receiver operating characteristic (ROC) curve analysis. Gene Ontology (GO) and Gene Set Enrichment Analysis (GSEA) were used to conduct a functional enrichment analysis, and miRNA-mRNA networks were established. The CIBERSORT algorithm was applied to identify the immune cell infiltration between the disease group and control group. Next, the correlations between key genes and infiltrating immune cells were investigated by the Spearman method. Finally, the correlation between 2 key genes and ferroptosis markers was confirmed.

Results

Nine ferroptosis differentially expressed genes (DE-FRGs) between DR and AS were identified in this study. NOX4 and PARP14 were selected as key genes for further analysis by external datasets and ROC curve analysis. The key genes NOX4, PARP14 and their correlated genes (such as CYBA, NOX1, NOX3, CYBB, PARP9, PARP10, and PARP15) are mainly enriched in oxidoreductase activity, protein ADP-ribosylation, superoxide metabolic process, reactive oxygen species metabolic process, PID pathway, and VEGFA-VEGFR2 pathway. A miRNA-mRNA network was constructed, and we got 12 miRNAs correlated with the target gene NOX4, 38 miRNAs correlated with the target gene PARP14. Three common miRNAs (hsa-miR-1-3p, hsa-miR-129-2-3p, and hsa-miR-155-5p) were observed in the network. Immune infiltration analysis displayed that activated B cell, MDSC, and Type 17 T helper cell are the common immune cells involved in the immune infiltration process of DR and AS. The results revealed that there are significant correlations between two key genes and most ferroptosis marker genes no matter in DR or AS.

Conclusion

Ferroptosis-related genes NOX4 and PARP14 may be common biomarkers of DR and AS. Both were associated with immune infiltration in patients with DR and AS. Our data provide a theoretical basis for the early diagnosis and immunotherapy of the two diseases.

Introduction

Diabetic retinopathy (DR) is the most common and devastating ocular complication in diabetes [1]. Atherosclerosis (AS) is also one of the common complications in diabetic patients. The association between DR and AS has attracted a large portion of research attention in recent years [2,3,4]. The American Diabetes Association (ADA) points out that DR is an independent predictor of type 2 diabetes (T2D) subclinical cardiovascular disease [5]. In T2D, DR alone or with other microvascular complications is associated with carotid atherosclerosis in patients [6, 7]. Furthermore, DR is associated with the occurrence of carotid atherosclerosis in type 1 diabetes mellitus (T1D) patients without history of cardiovascular disease. The presence of DR can determine the risk of carotid atherosclerotic in T1D patients [8]. The occurrence and development of DR and AS are very insidious and difficult to find. When vision loss or vascular damage has occurred, the best time for diagnosis and treatment is missed. Thus, this has major implications to further investigate the pathogenesis of DR and AS and identify new biomarkers for diagnosis.

Ferroptosis is an unconventional pattern of programmed cell death dependent on ferric ions. The morphology of ferroptosis is mitochondria shrank, condensed mitochondrial membrane density, mitochondrial crista reduction or disappearance [9]. The occurrence of ferroptosis and the inactivation of GPX4 lead to the increase of intracellular lipid reactive oxygen species (ROS), which attack lipids, induce lipid peroxidation, and leads to cell death [3]. Endothelial dysfunction caused by ROS is the initial link of AS occurrence [10]. Recent studies revealed that ferroptosis is highly linked with the pathogenesis of AS [11,12,13]. In addition, studies have also shown the function of ferroptosis in DR. The key characteristics of the occurrence and development of DR are the damage of retinal pigment epithelial (RPE) cells, which causes blood-retinal barrier disrupted, and the increase of the permeability of human retinal capillary endothelial cells (HRCECs). In DR, ferroptosis is involved in RPE cells and HRCECs death process [14, 15]. However, the pathogenesis, signaling pathways, and biomarkers of DR and AS related to ferroptosis remain unknown.

Bioinformatics methods enable us to understand the potential molecular mechanisms of diseases at a genetic level more deeply by screening molecules and displaying differences between patients and healthy individuals. This study aimed to identify hub genes associated with the pathogenesis of DR and AS. We obtained ferroptosis-related DEGs from the GEO database and FerrDb V2s database. Using integrated bioinformatics, enrichment analysis, and immune infiltration analysis, we identified common ferroptosis-related biomarkers for DR and AS and the associated immune cell infiltration status. These findings provide new insights into the molecular mechanisms involved in the pathophysiological processes of DR and AS.

Materials and methods

Microarray data retrieval

Four datasets (GSE60436, GSE100927, GSE102485, and GSE57691) related to DR and AS were downloaded from NCBI GEO (http://www.ncbi.nlm.nih.gov/geo) [16]. The GSE60436 (Illuminai HumanWG-6 v3.0 expression beadchip) was generated by the GPL6884 platform that included 6 fibrovascular membranes (FVMs) of retina samples from proliferative diabetic retinopathy (PDR) and 3 control samples [17]. The GSE100927 (Agilent-039494 SurePrint G3 Human GE v2 8 × 60 K Microarray 039381 (Probe Name version)) is acquired from the GPL17077 platform and contained 69 atherosclerotic lesions arteries samples and 35 control samples [18]. To better identify the key genes, two external datasets (GSE102485 and GSE57691) were used. The GSE102485 (Illumina NextSeq 500 (Homo sapiens)) is acquired from the GPL18573 platform and included 22 neovascular proliferative membrane specimens from patients with PDR and 3 control retina samples [19]. The GSE57691 (Illumina HumanHT-12 V4.0 expression beadchip) is generated by the GPL10558 platform comprising 9 atherosclerosis specimens obtained from aortic occlusive disease (AOD) patients and 10 control aortic specimens of organ donors [20].

Acquisition of microarray data and identification of DE-FRGs

Microarray data were accessed from GEO using the “GEO query” R package. After data preprocessing, the normalized expression profile data were analyzed and DEGs were screened using the “Limma” package in R language, p < 0.05 and, | log2FC (Fold-change) |≥ 0.585 were used to define DEGs. The volcano plot and heatmap were created with R package “ggplot2” and “ComplexHeatmap” [21,22,23]. We used the FerrDb V2 database (http://www.zhounan.org/ferrdb) to acquire FRGs. Venn diagrams were created using the R package “VennDiagram” [24].

Identification of key genes in DR and AS

Wilcoxon test was used for comparing differences between the disease group and the control group, and p < 0.05 suggests the difference was statistically significant. Receiver operating characteristic (ROC) curve analysis was generated with the “pROC” package in R [25], and the area under the curve (AUC) was calculated to assess the diagnostic value.

miRNA-mRNA network construction

The miRNA target key genes prediction was performed by the TarBase V8.0 [26]. The mRNA-miRNA network was established and visualized through the Cytoscape (Version 3.7.1).

Functional enrichment analysis

The functional association of target genes was analyzed using GeneMANIA (http://genemania.org). Gene Ontology (GO) and pathway enrichment analysis are achieved through the Metascape website (http://metascape.org/). Gene Set Enrichment Analysis (GSEA) associated with two key genes was performed using the “clusterProfiler” package, with gene set c2 (cp.v7.2.symbols.gmt) [27, 28].

Analysis of immune infiltrating cells

The CIBERSORT method was adopted to investigate the immune cell infiltration between the disease group and the control group [29]. T-tests were used for immune infiltrating cell comparisons between the two groups. The Boxplot of immune cell proportions was visualized with the ggplot2 package (http://ggplot2.tidyverse.org). Correlation analysis between key genes and immune cells was identified using the Spearman correlation coefficient, and we visualized the results with lollipop plots.

Statistical analysis

R software (version 3.6.3) and GraphPad Prism 5 were used for statistical analysis in this work. Data were expressed as mean ± standard deviation and unpaired Student’s t-test was employed for group comparisons. p < 0.05 suggests the difference was statistically significant.

Results

Identification of DE-FRGs in DR and AS

Figure 1 illustrates the flowchart of dataset selection and the study design. 484 FRGs were acquired from the FerrDb V2 database, which was a database that commits to ferroptosis regulators and ferroptosis-related diseases association. The expressions of the DEGs in GSE60436 and GSE100927 were displayed in volcano plots and heatmaps (Fig. 2A–D). Through differential expression analysis of GSE60436 and GSE100927, 307 common genes were significantly differentially expressed in both DR and AS (Fig. 2E). We obtained 9 DE-FRGs including TLR4, PARP12, PARP14, CYBB, NOX4, ABCC5, BID, LIFR and SNCA by taking the intersection of DEGs and FRGs (Fig. 2F).

Fig. 1
figure 1

Flowchart of the bioinformatics analysis

Fig. 2
figure 2

Identification of the DE-FRGs. A, B Volcano plot of DEGs in GSE60436 and in GSE100927, Red dots represent the upregulated, green dots denote the downregulated genes. C, D Heatmap of DEGs in GSE60436 and in GSE100927, Red bricks represent the more highly expressed DEGs, blue bricks denote lower expression DEGs. E Venn diagram displaying the crosstalk genes between GSE60436 and GSE100927. F Venn diagram displaying the DE-FRGs by taking the intersection of crosstalk genes and FRGs

Identification of the key genes in DR and AS

To further identify the key genes and evaluate their diagnostic value, we observed the expressions of 9 DE-FRGs in two external validation datasets GSE102485 and GSE57691 to get robust diagnostic key genes. Among these genes, we found a significant difference in the expression of ABCC5, NOX4, PARP12, PARP14, and TLR4 in the DR group in comparison with the control in GSE102485 (Fig. 3A–E). In the dataset GES57691, the expression level of CYBB, LIFR, NOX4, and PARP14 showed a significant difference in the AS group when compared to controls (Fig. 3F–I). Of these, NOX4 and PARP14 showed a significant difference in the two disease groups compared with controls. The diagnostic accuracy of NOX4 and PARP14 in GSE57691 was assessed with ROC curves and AUC. It was found that both NOX4 and PARP14 had high diagnostic values with AUC = 0.80 and AUC = 1.00 (Fig. 3J–M). Collectively, NOX4 and PARP14 were considered key genes and were chosen for subsequent analysis.

Fig. 3
figure 3

Verification of key genes for DR and AS. A–E The expressions of partial DE-FRGs in the GSE102485 dataset. F–I The expressions of partial DE-FRGs in the GSE57691 dataset. J–M ROC curve and AUC evaluated the diagnostic values of key genes in GSE57691 dataset. *p < 0.05, **p < 0.01, ***p < 0.001

Construction of the miRNA-mRNA network

To accurately identify target genes for miRNAs, we built a miRNA-mRNA network using data from the TarBase V8.0 database and visualized the network through Cytoscape. As shown in Fig. 4, from the results, we got 12 miRNAs linked to the target gene NOX4, and 38 miRNAs linked to the target gene PARP14. Three common miRNAs (hsa-miR-1-3p, hsa-miR-129-2-3p, and hsa-miR-155-5p) were observed in the network.

Fig. 4
figure 4

The mRNA-miRNA network. Orange color denotes target genes and yellow color denotes miRNA. Yellow color shows the common miRNA

Functional enrichment analysis of key genes in DR and AS

The gene–gene interaction network and significantly enriched biological processes of NOX4 and PARP14 were identified through GeneMANIA. The significantly enriched biological processes included oxidoreductase activity, protein ADP-ribosylation, superoxide metabolic process, reactive oxygen species metabolic process, etc. (Fig. 5A). The Metascape analysis displayed the top clusters of enriched sets (Fig. 5B). The key genes showed enrichment in metabolic process, biological process involving in interspecies interaction between organisms, immune system process, etc. The pathway data were enriched in the PID pathway and VEGFA-VEGFR2 signaling pathway, etc. GSEA associated with two key genes was subsequently conducted to compare the different pathways between the disease group and controls. NOX4 and PARP14 were significantly enriched in histidine metabolism, phenylalanine metabolism, and phototransduction in DR (Fig. 5C, D), while significantly enriched in circadian rhythm and dilated cardiomyopathy in AS (Fig. 5E, F). In conclusion, NOX4, PARP14, and their related genes are mainly enriched in oxidative stress, metabolic processes, immune system, and VEGFA-VEGFR2 signaling pathway.

Fig. 5
figure 5

Functional Enrichment Analysis of key genes. A Gene interaction network was analyzed with GeneMANIA. The inner circle represents key genes, the outer circle represents genes that interacts with both key genes. B Significantly enriched GO terms and KEGG pathways of key genes. Heatmap colored by P-values. CF GSEA of NOX4 and PARP14 in DR and AS

Immunocyte infiltration and correlation in DR and AS

The CIBERSORT method was applied to evaluate the abundance degree of various immune cells. We further calculate the significant differences in the proportion of immune cells between disease and controls. Differences in activated CD4 T cell, natural killer T cell, T cells regulatory, Type 2 helper T cells (p < 0.001), activated CD8 T cell, effector memory CD8 T cells, MDSC, natural killer cell (p < 0.01), activated dendritic cell, eosinophil, immature B cell, mast cell, neutrophil, and Type 1 helper T cells (p < 0.05) between the DR and control samples were significant (Fig. 6A). Almost all the immune cells expression were significantly different between the AS and control group (Fig. 6B). Correlation analysis between key genes and immune cells was identified using Spearman correlation coefficient, and we visualized the results with lollipop plots. As illustrated in Fig. 6C, the expression of NOX4 demonstrated a strong positive correlation with infiltration abundance of natural killer T cell, natural killer cell, and T cells regulatory, whereas NOX4 expression and activated B cell, Type 17 T helper cell showed a negative correlation in DR. Figure 6E showed that positive correlations were observed between PARP14 expression and MDSC, activated dendritic cell, and activated CD4 T cell, while negatively correlated with Type 17 T helper cell in DR. Figure 6D revealed that there was a negative relationship between NOX4 expression and eosinophil, activated B cell and monocyte in AS. As Fig. 6F shows, the expression of PARP14 was positively correlated with immature B cells, activated B cells, and MDSC in AS. Collectively, activated B cells, MDSC, and Type 17 T helper cells are common immune cells involved in the immune infiltration process of DR and AS.

Fig. 6
figure 6

Analysis of immune infiltration in DR and AS. A Differential expression of infiltrating immune cells between the disease group and control group in GSE60436. B Differential expression of infiltrating immune cells between the disease group and control group in GSE100927. Significant differences are shown in red font. CF Correlation analysis between key genes and various immune cells. The more significant the correlation, the larger and the lighter color of the dot. *p < 0.05, **p < 0.01, ***p < 0.001

Correlation between key genes and ferroptosis marker genes

To provide further evidence on the relationship between key genes and ferroptosis, we then downloaded 9 ferroptosis marker genes from the FerrDb database and observed the connection between key genes and ferroptosis marker genes. The results revealed that there are significant correlations between key genes and most ferroptosis marker genes no matter in DR or AS. In all, among these ferroptosis marker genes, PTGS2, TF, CHAC1, SLC40A1, FTH1, HSPB1 are highly correlated with NOX4 and PARP14 (Fig. 7).

Fig. 7
figure 7

Correlation between key genes and ferroptosis marker genes A Correlation between key genes and ferroptosis marker genes in DR. B Correlation between key genes and ferroptosis marker genes in AS. The more significant the correlation, the larger the dot size. Red denotes a positive correlation and blue denotes a negative correlation. *p < 0.05

Discussion

This study employed various bioinformatics methods. Firstly, microarray datasets of DR and AS were obtained from GEO, and 307 common DEGs were identified and screened. These were intersected with 484 FRGs extracted from the FERDB database to identify 9 DE-FRGs, including TLR4, PARP12, PARP14, CYBB, NOX4, ABCC5, BID, LIFR, and SNCA. External dataset validation revealed significant differences in NOX4 and PARP14 between disease and control groups of DR and AS, and these two genes showed high diagnostic value. By establishing a miRNA-mRNA network related to NOX4 and PARP14, three miRNA (hsa-miR-1-3p, hsa-miR-129–2-3p, and hsa-miR-155-5p) associated with both diseases were observed. Enrichment analysis showed that key genes were enriched in metabolic processes, biological processes of biointeractions, immune system processes, etc. Pathway data were enriched in the PID pathway, VEGFA-VEGFR2 signaling pathway, etc. Immune infiltration analysis showed that activated B cells, MDSC, and Th17 cells were common immune cells involved in the immune infiltration process of DR and AS. The results of this study contribute to a better understanding of ferroptosis and its role in DR and AS, providing a theoretical basis for immune therapy strategies for DR and AS.

Studies have confirmed the correlation between DR and AS. Retinal microangiogenesis is a major sign of PDR, and angiogenesis is also a common manifestation of progressive AS. The mechanisms of microvascular dysfunction share commonalities with those of large vessel injury. Some biological pathways and mechanisms affect both the development of DR and AS, such as oxidative stress injury mechanisms, chronic inflammation, and growth factors [30]. Ferroptosis is a non-conventional mode of programmed cell death dependent on iron ions. Studies have shown that ferroptosis plays a crucial role in the regulation of diabetes, ocular diseases, and cardiovascular diseases [31]. However, the exact role of ferroptosis in the pathogenesis of DR and AS is not clear, and there is currently a lack of effective treatment strategies. This study identified two key genes related to iron drooping in DR and AS, NOX4 and PARP14.

NOX4 is a regulator of various metabolic processes. In the eye, NOX4 is the main source of ROS in retinal endothelial cells. In the retina of T1D and T2D mice, the expression level of NOX4 is significantly increased. It has also been found to be associated with plaque stability, cell apoptosis, and plaque bleeding in AS [32]. PARP14 (also known as ARTD8 or BAL2) is a member of the PARP family, composed of macro1, macro2, macro3, WWE, and catalytic structural domains. PARP14 uses nicotinamide adenine dinucleotide as a metabolic substrate to mono ADP-ribosylate target proteins, participating in cell responses and immune system signaling pathways. Therefore, PARP14 is considered an attractive target for the treatment of tumors and allergic inflammation [33]. It has been reported that PARP is involved in lipid metabolism and regulates lipid metabolism and homeostasis through transcription factors. The activation of PARP family members is usually a destructive signal of lipid metabolism, and changes in lipid metabolism related to PARP play a core role in the occurrence of cardiovascular aging, atherosclerosis, obesity, fatty liver, hyperlipidemia, skin lesions, type 2 diabetes and its complications [34]. In addition, through enrichment analysis, NOX4, PARP14, and their related genes are mainly enriched in oxidative stress, metabolic processes, the immune system, and the VEGFA-VEGFR2 signaling pathway, etc. Studies have confirmed that destructive oxidative stress is one of the important mechanisms for the occurrence and development of DR and AS [35, 36].

This study observed three miRNAs (hsa-miR-1-3p, hsa-miR-129–2-3p, and hsa-miR-155-5p) related to both NOX4 and PARP14 through the establishment of a miRNA-mRNA network. It has been reported that miRNA plays a core role in cell proliferation, cell apoptosis, and angiogenesis in DR and AS. Researchers have found that hsa-miR-1-3p is a key inhibitor of the occurrence and deterioration of various human malignant tumors [37, 38]. Recent studies have found that hsa-miR-129–2-3p can be used as a biomarker for predicting the severity of DR and early diagnosis of DR [39]. Overexpression of miR-129-2-3p can regulate inflammation and cell apoptosis, promoting wound healing in diabetic patients [40]. However, there are no reports in AS. It has been reported that hsa-miR-155-5p is a medium connecting cellular aging and vascular aging, and its expression level in hVSMCs of AS patients is significantly down-regulated. These findings suggest that hsa-miR-155-5p plays a crucial role in cellular aging [41]. In addition, hsa-miR-155-5p is involved in the pathogenesis of diabetic nephropathy, diabetic neuropathy, and other diabetic complications, but there are few reports in DR. The enrichment analysis results of this study also show that NOX4, PARP14, and their related genes are enriched in the immune system. Our study found that activated B cells, MDSC, and Th17 cells are common immune cells involved in the immune infiltration process of DR and AS. Previous studies have shown that various activities of immune cells such as CD4 T cells, B cells, macrophages, etc., are involved in the pathogenesis of DR and AS [42, 43].

Identifying disease biomarkers through bioinformatics is a common and effective method. Sun et al. [44] used bioinformatics analysis to mine PDR-related genes and pathways. They identified a pathway and eight core genes associated with PDR. The main pathway involved is the cytokines-cytokines receptors interaction, and the core genes are tumor necrosis factor (TNF), tumor necrosis factor receptor superfamily member 12A (TNFRSF12A), C–C chemokine 20 (CCL20), chemokine (C–X–C motif) ligand 2 (CXCL2), oncostatin M (OSM), interleukin 10 (IL10), interleukin 15 (IL15), and interleukin 1B (IL1B). Zhao et al. [45] used bioinformatics methods to explore the gene expression changes in the retina of streptozotocin-induced diabetic rats, obtaining six differential genes. Our study focused on genes related to ferroptosis in DR and AS, so the core genes and pathways were not consistent with the above studies. A common approach is to analyze the correlation between two diseases through comprehensive bioinformatics analysis. Mo et al. [46] discussed the molecular mechanism between osteoporosis and AS. This study is similar to our research method, including identifying DEGs through the GEO database, performing enrichment analysis through GO and KEGG, and constructing a miRNA-gene network. The same method was also applied to the relationship between vascular dementia (VaD) and AS and the evaluation of potential biomarkers. The study identified seven hub genes and screened out two miRNAs significantly associated with AS and VaD [47]. Our study used the above methods to identify ferroptosis-related biomarkers in DR and AS, aiming to analyze the role of ferroptosis in the occurrence and development of DR and AS and to explore related targets. So far, there have been few reports on this topic. This study has certain innovations. The results of this study provide a new perspective on the common molecular mechanism between DR and AS. These common pathways and hub genes may provide promising clues for further experimental research.

Of course, this study inevitably has some limitations. First, only two datasets were included in this study for DR and AS, ignoring the impact of population heterogeneity in different countries on the results. Second, this study is a secondary analysis of existing data, lacking experimental and clinical data validation. Further research includes testing the role of two key genes through animal experiments and collecting clinical samples.

Conclusion

This study provided an in-depth understanding of the underlying pathogenesis related to ferroptosis and identified new biomarkers NOX4 and PARP14 for DR and AS, the presence of validated biomarkers may provide novel avenues for the early diagnosis of the two diseases in the future.

Data availability statement

Publicly available datasets were utilized in this study. The data used to support the results of this study are available from the online website mentioned above, and further inquiries can be directed to the corresponding author.

Data availability

Data will be made available on request.

References

  1. Kropp M, et al. Diabetic retinopathy as the leading cause of blindness and early predictor of cascading complications-risks and mitigation. EPMA J. 2023;14(1):21–42.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Castelblanco E, et al. Diabetic retinopathy predicts cardiovascular disease independently of subclinical atherosclerosis in individuals with type 2 diabetes: a prospective cohort study. Front Cardiovasc Med. 2022;9: 945421.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Fang Q, et al. Subclinical atherosclerosis associates with diabetic retinopathy incidence: a prospective study. Acta Diabetol. 2022;59(8):1041–52.

    Article  PubMed  Google Scholar 

  4. Eid M, et al. Diabetic retinopathy as a predictor of angiographic coronary atherosclerosis severity in patients with type 2 diabetes mellitus. Diabetes Metab Syndr Obes. 2022;15:1485–94.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Simó R, et al. Diabetic retinopathy as an independent predictor of subclinical cardiovascular disease: baseline results of the PRECISED study. BMJ Open Diabetes Res Care. 2019;7(1): e000845.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Gao L, et al. Proliferative diabetic retinopathy in patients with type 2 diabetes correlates with the presence of atherosclerosis cardiovascular disease. Diabetol Metab Syndr. 2021;13(1):48.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Zhang C, et al. Association between atherosclerosis and diabetic retinopathy in chinese patients with type 2 diabetes mellitus, diabetes. Metab Syndr Obes. 2020;13:1911–20.

    Article  CAS  Google Scholar 

  8. Carbonell M, et al. Diabetic retinopathy is associated with the presence and burden of subclinical carotid atherosclerosis in type 1 diabetes. Cardiovasc Diabetol. 2018;17(1):66.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Mou Y, et al. Ferroptosis, a new form of cell death: opportunities and challenges in cancer. J Hematol Oncol. 2019;12(1):34.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Maiocchi SL, et al. Myeloperoxidase: a versatile mediator of endothelial dysfunction and therapeutic target during cardiovascular disease. Pharmacol Ther. 2021;221: 107711.

    Article  CAS  PubMed  Google Scholar 

  11. Lin L, et al. Autophagy, pyroptosis, and ferroptosis: new regulatory mechanisms for atherosclerosis. Front Cell Dev Biol. 2022;9: 809955.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Bai T, et al. Inhibition of ferroptosis alleviates atherosclerosis through attenuating lipid peroxidation and endothelial dysfunction in mouse aortic endothelial cell. Free Radic Biol Med. 2020;160:92–102.

    Article  CAS  PubMed  Google Scholar 

  13. Ouyang S, et al. Ferroptosis: the potential value target in atherosclerosis. Cell Death Dis. 2021;12(8):782.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Zhang J, et al. TRIM46 contributes to high glucose-induced ferroptosis and cell growth inhibition in human retinal capillary endothelial cells by facilitating GPX4 ubiquitination. Exp Cell Res. 2021;407(2): 112800.

    Article  ADS  CAS  PubMed  Google Scholar 

  15. Tang X, et al. Astragaloside-IV alleviates high glucose-induced ferroptosis in retinal pigment epithelial cells by disrupting the expression of miR-138-5p/Sirt1/Nrf2. Bioengineered. 2022;13(4):8240–54.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Barrett T, et al. NCBI GEO: archive for functional genomics data sets—update. Nucleic Acids Res. 2013;41(D1):D991-995.

    Article  CAS  PubMed  Google Scholar 

  17. Ishikawa K, et 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.

    Article  CAS  PubMed  Google Scholar 

  18. Steenman M, et al. Identification of genomic differences among peripheral arterial beds in atherosclerotic and healthy arteries. Sci Rep. 2018;8(1):3940.

    Article  ADS  PubMed  PubMed Central  Google Scholar 

  19. Li Y, et al. Induced expression of VEGFC, ANGPT, and EFNB2 and their receptors characterizes neovascularization in proliferative diabetic retinopathy. Invest Ophthalmol Vis Sci. 2019;60(13):4084–96.

    Article  CAS  PubMed  Google Scholar 

  20. Biros E, et al. Differential gene expression in human abdominal aortic aneurysm and aortic occlusive disease. Oncotarget. 2015;6(15):12984–96.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Ritchie ME, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7): e47.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Gustavsson EK, et al. ggtranscript: an R package for the visualization and interpretation of transcript isoforms using ggplot2. Bioinformatics. 2022;38(15):3844–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Gu Z, et al. Complex heatmaps reveal patterns and correlations in multidimensional genomic data. Bioinformatics. 2016;32(18):2847–9.

    Article  CAS  PubMed  Google Scholar 

  24. Chen H, et al. VennDiagram: a package for the generation of highly-customizable Venn and Euler diagrams in R. BMC Bioinformatics. 2011;12:35.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Robin X, et al. pROC: an open-source package for R and S+ to analyze and compare ROC curves. BMC Bioinformatics. 2011;12:77.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Karagkouni D, et al. DIANA-TarBase v8: a decade-long collection of experimentally supported miRNA-gene interactions. Nucleic Acids Res. 2018;46(D1):D239–45.

    Article  CAS  PubMed  Google Scholar 

  27. Subramanian A, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005;102(43):15545–50.

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  28. Yu G, et al. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16(5):284–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Newman AM, et al. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods. 2015;12(5):453–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Ellis TP, et al. Diabetic retinopathy and atherosclerosis: is there a link? Curr Diabetes Rev. 2013;9(2):146–60.

    Article  PubMed  Google Scholar 

  31. Stockwell BR, et al. The chemistry and biology of ferroptosis. Cell Chem Biol. 2020;27(4):365–75.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Zhou T, et al. microRNA-363-3p reduces endothelial cell inflammatory responses in coronary heart disease via inactivation of the NOX4-dependent p38 MAPK axis. Aging (Albany NY). 2021;13(8):11061–82.

    Article  CAS  PubMed  Google Scholar 

  33. Qin W, et al. Research progress on PARP14 as a drug target. Front Pharmacol. 2019;10:172.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Szántó M, et al. PARPs in lipid metabolism and related diseases. Prog Lipid Res. 2021;84: 101117.

    Article  PubMed  Google Scholar 

  35. Bianchetti G, et al. Investigation of DHA-induced regulation of redox homeostasis in retinal pigment epithelium cells through the combination of metabolic imaging and molecular biology. Antioxidants (Basel). 2022;11(6):1072.

    Article  CAS  PubMed  Google Scholar 

  36. Kotozaki Y, et al. Human plasma xanthine oxidoreductase activity in cardiovascular disease: evidence from a population-based study. Biomedicines. 2023;11(3):754.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Mao Z, et al. Post-transcriptional dysregulation of microRNA and alternative polyadenylation in colorectal cancer. Front Genet. 2020;11:64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Song Y, et al. Dysregulation of pseudogenes/lncRNA-Hsa-miR-1-3p-PAICS pathway promotes the development of NSCLC. J Oncol. 2022;2022:4714931.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Wu J, et al. A 3-miRNA risk scoring signature in early diabetic retinopathy. J Clin Med. 2023;12(5):1777.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Umehara T, et al. Identification of specific miRNAs in neutrophils of type 2 diabetic mice: overexpression of miRNA-129-2-3p accelerates diabetic wound healing. Diabetes. 2019;68(3):617–30.

    Article  CAS  PubMed  Google Scholar 

  41. Nguyen DDN, et al. Intracellular and exosomal microRNAome profiling of human vascular smooth muscle cells during replicative senescence. Am J Physiol Heart Circ Physiol. 2021;321(4):H770–83.

    Article  CAS  PubMed  Google Scholar 

  42. Wu MY, et al. New insights into the role of inflammation in the pathogenesis of atherosclerosis. Int J Mol Sci. 2017;18(10):2034.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Tang X, et al. Sustained upregulation of endothelial nox4 mediates retinal vascular pathology in type 1 diabetes. Diabetes. 2023;2(1):112–25.

    Article  Google Scholar 

  44. Sun H, et al. Mining the proliferative diabetic retinopathy-associated genes and pathways by integrated bioinformatic analysis. Int Ophthalmol. 2020;40(2):269–79.

    Article  PubMed  Google Scholar 

  45. Zhao W, et al. Bioinformatic analysis of retinal gene function and expression in diabetic rats. Exp Ther Med. 2017;14(3):2485–92.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Mo L, et al. Integrated bioinformatic analysis of the shared molecular mechanisms between osteoporosis and atherosclerosis. Front Endocrinol (Lausanne). 2022;22: 950030.

    Article  Google Scholar 

  47. Li D, et al. Bioinformatic identification of potential biomarkers and therapeutic targets in carotid atherosclerosis and vascular dementia. Front Neurol. 2023;13:1091453.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Funding

None.

Author information

Authors and Affiliations

Authors

Contributions

CL: Conceptualization, methodology, data curation, writing—original draft, writing. QC: conceptualization, supervision, project administration.

Corresponding author

Correspondence to QinHua Cai.

Ethics declarations

Ethical approval and consent to participate

GEO database is a publicly accessible database and includes patient data that have obtained ethical approval. There are no other required ethical statements.

Competing interests

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Li, C., Cai, Q. Two ferroptosis-specific expressed genes NOX4 and PARP14 are considered as potential biomarkers for the diagnosis and treatment of diabetic retinopathy and atherosclerosis. Diabetol Metab Syndr 16, 61 (2024). https://doi.org/10.1186/s13098-024-01301-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13098-024-01301-3

Keywords