Hui Guo
Kui Xiao2†
Jianchun Zong- 1Department of Emergency, The Second Affiliated Hospital, Chongqing Medical University, Chongqing, China
- 2Department of Plastic Surgery, Guangzhou Red Cross Hospital, Jinan University, Guangzhou, China
- 3Department of Critical Medicine, Wusong Hospital, Zhongshan Hospital, Fudan University, Shanghai, China
Background: Diabetic foot ulcers are the most common and serious complication of diabetes mellitus, the high morbidity, mortality, and disability of which greatly diminish the quality of life of patients and impose a heavy socioeconomic burden. Thus, it is urgent to identify potential biomarkers and targeted drugs for diabetic foot ulcers.
Methods: In this study, we downloaded datasets related to diabetic foot ulcers from gene expression omnibus. Dysregulation of mitophagy-related genes was identified by differential analysis and weighted gene co-expression network analysis. Multiple machine algorithms were utilized to identify hub mitophagy-related genes, and a novel artificial neural network model for assisting in the diagnosis of diabetic foot ulcers was constructed based on their transcriptome expression patterns. Finally, potential drugs that can target hub mitophagy-related genes were identified using the Enrichr platform and molecular docking methods.
Results: In this study, we identified 702 differentially expressed genes related to diabetic foot ulcers, and enrichment analysis showed that these genes were associated with mitochondria and energy metabolism. Subsequently, we identified hexokinase-2, small ribosomal subunit protein us3, and l-lactate dehydrogenase A chain as hub mitophagy-related genes of diabetic foot ulcers using multiple machine learning algorithms and validated their diagnostic performance in a validation cohort independent of the present study (The areas under roc curve of hexokinase-2, small ribosomal subunit protein us3, and l-lactate dehydrogenase A chain are 0.671, 0.870, and 0.739, respectively). Next, we constructed a novel artificial neural network model for the molecular diagnosis of diabetic foot ulcers, and the diagnostic performance of the training cohort and validation cohort was good, with areas under roc curve of 0.924 and 0.840, respectively. Finally, we identified retinoic acid and estradiol as promising anti-diabetic foot ulcers by targeting hexokinase-2 (−6.6 and −7.2 kcal/mol), small ribosomal subunit protein us3 (−7.5 and −8.3 kcal/mol), and l-lactate dehydrogenase A chain (−7.6 and −8.5 kcal/mol).
Conclusion: The present study identified hexokinase-2, small ribosomal subunit protein us3 and l-lactate dehydrogenase A chain, and emphasized their critical roles in the diagnosis and treatment of diabetic foot ulcers through multiple dimensions, providing promising diagnostic biomarkers and targeted drugs for diabetic foot ulcers.
1 Introduction
Diabetes mellitus (DM) is a chronic metabolic disease characterized by hyperglycemia (Brody, 2012), and the number of patients with DM continues to increase as the global society ages (Hart et al., 2017). Among the many diabetes-induced complications, the problem of diabetic foot ulcers (DFU) appears to be particularly serious and is widely regarded as one of the most challenging (GBD, 2019 Diseases and Injuries Collaborators, 2020). Recent studies have shown that diabetic foot ulcers are the main contributor to non-traumatic lower limb amputations, with diabetic foot ulcers leading to 40% to 60% of all non-traumatic amputations (Jude et al., 2001). And diabetic foot accounts for 19%∼34% of diabetic patients, with a high mortality rate of 50%∼70% within 5 years (McDermott et al., 2023), and more than 1/3 of diabetic foot ulcers require amputation surgery (Armstrong et al., 2023). The high morbidity, mortality, and disability of the diabetic foot greatly diminishes patients’ quality of life and imposes a heavy socioeconomic burden (Armstrong et al., 2017). Additionally, DFU is also a complex mixture of neuropathy and infection (Villa et al., 2024). Previous studies have shown that foot infections also increase the risk of lower limb amputation in DM patients (Pitocco et al., 2019; Ruder, 2024; Wei et al., 2024), and its prevalence has risen from 4% to 7% (Pitocco et al., 2019). Diabetic foot infection (DFI) is a pathological condition that is primarily caused by the invasion and multiplication of microorganisms in host tissues resulting in an inflammatory response that usually results in tissue damage (Coye et al., 2024; Senneville et al., 2024; Villa et al., 2024; Wang B. et al., 2024). DFI usually occur after the breakage of the skin’s protective coating, and the most common breakage is the DFU, which involves at least the epidermis and part of the dermis (Senneville et al., 2024). Due to the fact that DFUs are often associated with vasculopathy and recurrent microbial infections, which contributes to high patient morbidity and mortality (Xiong et al., 2020). Thus, the main aim of DFU treatment is to control infection and promote wound healing (Primadhi et al., 2023).
However, the current clinical diagnosis and treatment of DFU also face challenges, with the clinical diagnosis of DFU being overly dependent on the clinical experience of clinicians and lacking effective treatments (Armstrong et al., 2017; Matos et al., 2018; Purwanti et al., 2024; Wang F. et al., 2024; Wang K. et al., 2024). At the present stage, the treatment of DFU patients mainly includes debridement, wound dressing, anti-infection treatment, peripheral vasculopathy treatment, strict glycemic control, and amputation (McIllhatton et al., 2021; Orlando et al., 2021; Ahmed et al., 2022; Slomski, 2022; Jeffcoate et al., 2024), but the therapeutic efficacy of these methods is still unsatisfactory and has not yet achieved a satisfactory clinical results, that is, they do not fully control the underlying metabolic lesions (Bardill et al., 2022; Chen et al., 2024; Yi et al., 2024; Zhang G. et al., 2024; Zhang N. et al., 2024). Thus, it is urgent to further identify novel biomarkers of DFU to improve the current dilemma of clinical diagnosis and treatment of DFU.
Mitochondria are the energy powerhouses of the cell, responsible not only for ATP production, but also for regulating calcium homeostasis, responding to oxidative stress, and modulating cell death processes (Li A. et al., 2022). Through its own specialized autophagic program, mitochondrial autophagy precisely removes destroyed mitochondria, and this mechanism is critical for maintaining mitochondrial functionality and ensuring intracellular homeostasis (Onishi et al., 2021). However, mitochondrial autophagy, an important mitochondrial quality control mechanism, leads to the gradual accumulation of mitochondrial DNA (mtDNA) mutations in response to reactive oxygen species (ROS) stress, as well as a decrease in intracellular mitochondrial membrane potential and depolarization damage, eventually leading to cell death (Lu et al., 2023). Previous studies have shown a strong association between DFU and mitophagy (Qi et al., 2023; Wang et al., 2023). The accumulation of ROS can damage the extracellular matrix and endothelial cells, impede the formulation of new blood vessels, and delay the wound healing process (Li et al., 2018). At the same time, the increase in ROS may also exacerbate the inflammatory response, making it difficult for diabetic foot ulcers to heal (Huang et al., 2020). Therefore, reducing the release of ROS appears to be crucial for improving the healing of DFU. Furthermore, in diabetic patients, a persistent hyperglycemic state may trigger impairment of mitochondrial function, which in turn exacerbates the cellular stress response and promotes cell death (Wronka et al., 2022). Enhanced oxidative stress and weakened protective mechanisms in diabetic patients can lead to the accumulation of ROS (Rendra et al., 2019). In this context, mitophagy activators present potential as novel therapeutic approaches for diabetes (Shan et al., 2022). By eliminating dysfunctional mitochondria, mitochondrial autophagy contributes to the reduction of ROS production and the alleviation of oxidative stress (Rovira-Llopis et al., 2017; Zhang et al., 2019), which is particularly important for the healing of DFU. Therefore, it is necessary to develop mitophagy-related genes for the diagnosis and treatment of DFU through whole transcriptome analysis.
In recent decades, machine learning has evolved from a peripheral technology to become a cornerstone of healthcare data analytics (Haug and Drazen, 2023; Shibue, 2023), which is highly effective in identifying disease-relevant biomarkers, making accurate clinical diagnoses, and determining optimal therapeutic regimens in the field of biomedicine (Le et al., 2021). Due to advances in high-throughput sequencing technology and the integration of machine learning in the medical field, many new methods for identifying molecular targets for various diseases are now available (Wu et al., 2023). Clinicians can now utilize biomarkers for early diagnosis of DFU to initiate timely interventions (Wang et al., 2021; Li Y. et al., 2022; Shi et al., 2023; Shi et al., 2024). It has been previously shown that the integration of bioinformatics and machine learning has become an important tool for identifying novel biomarkers associated with DFU (Khandakar et al., 2021; Chemello et al., 2022; Guan et al., 2024; Hettinger et al., 2024).
In this study, we systematically investigated the changes in the expression patterns of mitophagy-related genes in DFU based on DFU whole transcriptome expression data. Immediately after that, we found that hexokinase-2 (HK2), small ribosomal subunit protein us3 (RPS3), and l-lactate dehydrogenase A chain (LDHA) can be indicated using multiple bioinformatic algorithms, that can be used as hub mitophagy-related genes of DFU (DMGs) in DFU. Meanwhile, we validated the importance of DMGs in DFU through multiple dimensions. In conclusion, the present study takes mitophagy as the starting point, and through the combination of multiple bioinformatic analysis techniques and machine learning algorithms to identify novel targets that can be used for the diagnosis and treatment of DFU, which will provide promising targets for the clinical diagnosis and treatment of DFU.
2 Materials and methods
2.1 Data downloading and processing
In this study, the bulk RNA-seq datasets relevant to DFU were acquired from the gene expression omnibus (GEO) (https://www.ncbi.nlm.nih.gov/geo/), with Table 1 displaying the fundamental information of these datasets. The training cohort comprised 33 samples from GSE134431 and GSE80178, while the validation cohort included 32 samples from GSE7014 and GSE68183. “ComBat” from the R package “sva” was employed to remove batch effects within both the training and validation cohorts (Leek et al., 2012). The effectiveness of this de-batching effect was assessed through principal component analysis (PCA).
Table 1. The basic information regarding the dataset in the present study.
2.2 Gene expression difference analysis
In this study, we compared the whole transcriptome expression data of DFUs and Ctrls using the R package “limma” to identify differentially expressed genes (DEGs) in both clusters at a P-value of < 0.05. P-values were computed using the Wilcoxon rank-sum test.
2.3 Weighted gene co-expression network analysis
We used the “WGCNA” R package to identify DEGs associated with DFUs (Langfelder and Horvath, 2008). We first determined the appropriate soft threshold (β) and transformed the expression data matrix into a topological overlap matrix (TOM). Based on the TOM difference metric, we set the minimum module size to 200 and compute the module eigengenes (ME). Subsequently, gene significance (GS) and module significance (MS) were calculated by correlating the identified modules with the phenotypic data of the DFUs, which helped to analyze the correlation between each module and the DFUs.
2.4 Downloading of mitophagy-related genes
The mitophagy-related genes (MRGs) used in this study were obtained from the GeneCards database (https://www.genecards.org/). Specifically, we firstly entered the GeneCards database homepage, entered “mitophagy” as the keyword in the search panel, selected the gene type as “protein coding.” The score was set to greater than 2.5, resulting in a total of 417 MRGs (Supplementary Table S1).
2.5 Functional enrichment analysis
In order to identify the primary biological processes and functions the 22 DFU-related MRGs were enriched, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were conducted utilizing the R package “clusterProfler” (Yu et al., 2012).
2.6 Identification of DMGs using machine learning algorithms
We utilized three machine learning algorithms to identify DMGs: support vector machine recursive feature elimination (SVM-RFE), least absolute shrinkage and selection operator (LASSO), and random forest (RF). We applied the RF algorithm from the “randomForest” package, the LASSO algorithm from the “glmnet” package, and the SVM-RFE algorithm from the “e1071” package for screening the 22 DFU-related MRGs to identify potential candidate genes (Noble, 2006; Vasquez et al., 2016; Paul et al., 2018).
2.7 Evaluation of the diagnostic performance of DMGs for DFU
We utilized the R software package “pROC” to further evaluate the diagnostic potential of the identified DMGs for DFU. Using the “roc” function of the “pROC” package, we conducted receiver operating characteristic (ROC) analyses. The final area under the ROC curve (AUC) results were calculated using the “ci” function of “pROC.”
2.8 Construction of the ANN model
We used the R package “neuralnet” to develope an ANN model. Two data sets, GSE134431 and GSE80178, were used as the training set of the ANN, and two data sets, GSE7014 and GSE68183, were used as the validation set of the ANN. The resulting ANN model featured three main layers: an input layer, a hidden layer, and an output layer.
2.9 Gene set enrichment analysis
Based on the median value of the DMGs expression, we clustered the DFU cohort and conducted gene set enrichment analysis (GSEA) to ascertain the correlation between these genes and biological functions.
2.10 Immune infiltration analysis
We assessed immune cell infiltration in DFU patients using the R packages “GSVA” and “GSEABase” (Hänzelmann et al., 2013). First, genetic information of 28 immune cell genomes (Supplementary Table S2) was obtained from the TISIDB database (http://cis.hku.hk/TISIDB/). Subsequently, ssGSEA scores were calculated for each sample based on these 28 immune cell genomes. The infiltration levels of different immune cells were visualized using the “pheatmap” R package.
2.11 Molecular docking
We utilized the Enrichr platform (https://maayanlab.cloud/Enrichr/) to identify drugs targeting DMGs. The molecular structures of the drugs were sourced from the PubChem database (https://pubchem.ncbi.nlm.nih.gov/). Simultaneously, we retrieved the 3D structures of HK2 (PDB ID: 6Q6N; resolution: 1.63 Å), RPS3 (PDB ID: 4ZVV; resolution: 2.20 Å), and LDHA (PDB ID: 5W8J; resolution: 1.55 Å) from the PDB (https://www.rcsb.org/). Protein and drug files were converted to PDBQT format, excluding water molecules but including polar hydrogen atoms. Using AutoDock tools, we prepared the ligand and protein files, conducted docking simulations, and visualized the interactions with PyMOL.
2.12 Statistical analysis
In this study, we performed all statistical analyses of bioinformatics in R language. The analysis of variance (ANOVA) method was employed to statistically analyze multi-group data, while the wilcoxon rank sum test was used to compare two groups. Statistical significance was defined as P < 0.05 for all statistical analyses.
3 Results
3.1 Differential analysis of gene expression
In this study, we used the 33 samples contained in GSE134431 and GSE80178 as the training cohort, but due to the batch effect between different datasets, we need to remove the batch effect caused by the merging of GSE134431 and GSE80178 in the pre-processing of the data, to make the data in the same dimension and comparable. The data pattern of the training cohort before merging is shown in Figures 1A, B, which shows that the data are not in the same dimension in multiple aspects. After removing the batch effect, the data pattern is shown in Figures 1C, D. The results show that the samples are independent of each other and evenly dispersed, meaning that the data are in the same dimension and comparable, and can be included in subsequent analyses. Next, we used the corrected matrix for differential gene expression analysis (Figure 1E), resulting in a total of 1,665 DEGs, of which 854 were upregulated genes and 811 were down-regulated genes (Supplementary Table S3). The whole expression landscape of these DEGs was shown in Figure 1F, that is, their expression was significantly different between the control and DFU group. Finally, we visualized these 1,665 DEGs on chromosomes and showed that they were evenly distributed on chromosomes and underwent deregulated expression (Figure 1G).
Figure 1. Identification of DEGs based on bulk RNA-seq. (A,B) Distributional characteristics of the unprocessed data. (C,D) Distribution characteristics of the processed data. (E) Volcano plot showing the DEGs between DFUs and Controls. (F) Heatmap showing the overall landscape of the DEGs between DFUs and Controls. (G) Chromosome map showing the distribution of DEGs on chromosomes.
3.2 Identification of DFU-related DEGs
Dysregulated expression of genes may be caused by multiple factors, that is, to identify DFU-related DEGs, we further performed weighted gene co-expression network analysis. Based on the optimal soft threshold power β = 10 (unscaled R2 = 0.9), the genes were divided into four independent co-expression modules (Figures 2A, B). The clustering dendrogram depicted the sample grouping, with DFU samples being distinctly separated from the control samples (Figure 2C). The correlogram of module-trait relationships highlighted that the turquoise module showed the strongest correlation with DFU (Figure 2D). The turquoise module contained a total of 702 DEGs and had the highest gene percentage (Supplementary Table S4), that is, a total of 42.2% of the dysregulated expression of the DEGs identified in the results mentioned above were related to DFU (Figure 2E), and were thus included in the subsequent analysis as DFU-related DEGs.
Figure 2. Identification of DFU-related DEGs based on bulk RNA-seq. (A) Construction of the unscaled co-expression network. (B) Merging of similar modules. (C) Sample dendrogram and trait heatmap. (D) The correlation between modules and disease phenotypes. (E) Statistics on the proportion of the number of module genes.
3.3 The expression pattern of MRGs in DFU altered
To investigate whether the expression patterns of MRGs in DFUs were altered, we obtained 417 MRGs from the GeneCards database (https://www.genecards.org/) for subsequent studies (Figure 3A). We intersected these 417 MRGs with the 702 DFU-related DEGs identified in the analysis mentioned above, resulting in a total of 22 DFU-related MRGs (Figure 3B), implying that the deregulated expression of these genes affects the onset and progression of DFUs to a certain extent in the form of mitophagy. Immediately after that, we further validated the dysregulated expression of these genes, and the box plot showed that there was a significant difference in their expression between control and DFU groups (Figure 3C). However, genes are not independent of each other, and there would be a certain relationship between their own expression, so we performed correlation analysis on the expression of these 22 DFU-related MRGs, and the results showed that there was a strong expression relationship among them internally (Figures 3D, E). The dysregulated expression of DFU-related MRGs and the synergistic interactions between them suggest that the expression pattern of MRGs in DFU is altered resulting in their ability to contribute to the onset and progression of DFU.
Figure 3. Expression patterns of MRGs in DFU patients. (A) Heatmap showing the overall expression landscape of the 417 MRGs. (B) Venn plots showing the overlapped genes between 417 MRGs and 702 DFU-related DEGs. (C) Box plot showing the expression levels of DFU-related MRGs. (D) Gene relationship network diagram of DFU-related MRGs. (E) Correlation analysis of DFU-related MRGs. Red and green colors represent positive and negative correlations, respectively. The correlation coefficient was expressed as the area of the pie chart (*p < 0.05, **p < 0.01, ***p < 0.001).
3.4 Functional enrichment analysis
To ascertain the biological functions and signaling pathways in which MRGs are mainly enriched to contribute to the progression of DFU, we performed enrichment analysis on the 22 DFU-related MRGs mentioned above. Among them, GO results showed that they were mainly related to mitochondria/autophagy and energy generation and metabolism, such as “mitochondrial outer membrane,” “mitochondrial matrix,” “apoptotic mitochondrial changes,” “regulation of mitochondrial membrane permeability,” “macroautophagy,” “regulation of autophagy,” “aldehyde dehydrogenase [NAD(P)+] activity,” “aldehyde dehydrogenase (NAD+) activity,” “protein phosphatase binding,” “electron transfer activity,” “monosaccharide metabolic process,” “hexose metabolic process,” “reactive oxygen species metabolic process,” “glucose metabolic process,” “generation of precursor metabolites and energy,” “response to oxygen levels,” and “response to hypoxia” (Figure 4A). Additionally, KEGG results were similar to GO results, showing that they were associated with energy generation, amino acid metabolism, and cancer, that is, “Glycolysis/Gluconeogenesis,” “Pyruvate metabolism,” “Cysteine and methionine metabolism,” “Carbon metabolism,” “p53 signaling pathway,” “HIF-1 signaling pathway,” “NF-kappa B signaling pathway,” “Bladder cancer,” “Colorectal cancer,” “Small cell lung cancer,” “Prostate cancer,” and “Central carbon metabolism in cancer” (Figure 4B).
Figure 4. Enrichment analysis of DFU-related MRGs. (A) GO enrichment analysis results of DFU-related MRGs. (B) KEGG enrichment analysis results of DFU-related MRGs.
3.5 Identification of DMGs using multiple machine learning algorithms
In recent years, machine learning has become a major tool for healthcare data analysis and mining, and it is highly effective in identifying disease-related biomarkers and making accurate clinical diagnoses. Clinicians can now utilize biomarkers for early diagnosis of DFU and timely intervention. Thus, the present study analyzed the 22 DFU-associated MRGs mentioned above using multiple machine learning algorithms (Supplementary Table S5). First, we identified 7 candidate genes using the LASSO algorithm (Figure 5A); 7 candidate genes using the SVM-RFE algorithm (Figure 5B); and 11 candidate genes using the RF algorithm (Figure 5C). Subsequently, we crossed the identified candidate genes (Figure 5D), and found that HK2, RPS3, and LDHA could serve as DMGs (Table 2).
Figure 5. Identification of DMGs using three machine learning algorithms. (A) LASSO logistic regression was used to select candidate genes, accompanied by fine-tuning the penalty parameter through a process of tenfold cross-validation. (B) The SVM-RFE algorithm was utilized for candidate gene selection. (C) The RF algorithm was utilized for candidate gene selection. (D) Venn diagram showing the overlapped candidate genes by RF, SVM-RFE, and LASSO algorithms.
Table 2. Scanning of candidate machines by three machine learning algorithms.
3.6 Diagnostic performance of DMGs
The screened DMGs were significantly differentially expressed, suggesting that DMGs may play a potential role in contributing DFU (Figures 6A–C). Furthermore, AUC of DMGs was 0.955 of HK2, 0.909 of RPS3, and 0.822 of LDHA, respectively (Figure 6D). The results of the validation cohort showed that DMGs were also dysregulated, where RPS3 was down-regulated in DFU and HK2 and LDHA were both up-regulated in DFU (Figures 6E–G). DMGs also showed good performance in DFU in the validation set, with AUCs of 0.671 for HK2, 0.870 for RPS3 and 0.739 for LDHA (Figure 6H). Additionally, we also constructed an artificial neural network (ANN) model based on the DMGs. Two data sets, GSE134431 and GSE80178, were used as the training set of the ANN, and two data sets, GSE7014 and GSE68183, were used as the validation set of the ANN. The resulting ANN model featured three main layers: an input layer, a hidden layer, and an output layer (Figure 6I). Table 3 presents the results of the ANN prediction. The training set achieved a prediction accuracy of 90.9%, while the validation set yielded 87.5%. Furthermore, the ANN model has an AUC value of 0.924 in the training set (Figure 6J), and 0.840 in the validation set (Figure 6K). Overall, the ANN model constructed based on the transcriptome of DMGs can assist in the prediction of DFU.
Figure 6. Diagnostic performance of DMGs in DFU. (A–C) The expression of the DMGs in the training cohort. (D) The ROC of the DMGs in the training cohort. (E–G) The expression of the DMGs in the validation cohort. (H) The ROC of the DMGs in the validation cohort. (I) Construction of ANN based on the DMGs. (J) The AUC of the training set with a value of 0.924. (K) The AUC of the validation set with a value of 0.840.
Table 3. ANN diagnosis effect for the training and validation cohort.
3.7 GSEA analysis
We assessed signaling pathways associated with DMGs using GSEA analysis. The results showed that LDHA is enriched in type II diabetes mellitus, rheumatoid arthritis, ribosome, RNA degradation, valine, leucine and isoleucine degradation, fat digestion and absorption, and circadian rhythm (Figures 7A, B). RPS3 is enriched in ribosome, tyrosine metabolism, glycine serine and threonine metabolism, terpenoid backbone biosynthesis, and type II diabetes mellitus (Figures 7C, D). HK2 is enriched in IL-17 signaling pathway, type II diabetes mellitus, terpenoid backbone biosynthesis, systemic lupus erythematosus, and tyrosine metabolism (Figures 7E, F). As mentioned above, DMGs are enriched in diabetes-related biological pathways, such as “Type II diabetes mellitus,” which creates the foundation and necessary conditions for the development of DFU. Additionally, we found that DMGs are also associated with inflammatory signaling pathways and amino acid metabolism signaling pathways, such as “IL-17 signaling pathway,” “Valine, leucine and isoleucine degradation,” “Glycine, serine and threonine metabolism,” “Histidine metabolism,” and “Tyrosine metabolism.” There exists a strong relationship between amino acid metabolism and energy metabolism in organisms (e.g., glycolysis, TCA cycle, respiratory chain, etc.). Amino acids can be metabolized to produce ATP and also provide intermediates for other metabolic pathways to maintain normal life activities. However, dysregulation of amino acid metabolism can also lead to disturbed energy metabolism in organisms, which can activate mitochondrial abnormalities including mitochondrial autophagy, and thus induce inflammatory diseases.
Figure 7. The GSEA of DMGs in DFU. (A,B) The GSEA of LDHA in DFU. (C,D) The GSEA of RPS3 in DFU. (E,F) The GSEA of HK2 in DFU.
3.8 Immune cell infiltration analysis
We conducted an assessment of the immunological profile of DFU samples utilizing the ssGSEA algorithm (Supplementary Table S6). Figure 8A showed overall landscape of immune cell infiltration between the DFU and Control groups. Compared with Control samples, DFU samples have higher “Type 17 T helper cell,” “CD56dim natural killer cell,” “Activated dendritic cell” and “Eosinophil” “Neutrophil,” and have lower “Effector memory CD4 T cell,” “Immature B cell,” “Natural killer cell,” “CD56bright natural killer cell” (Figure 8B). Subsequently, we explored the correlation between the expression of DMGs and the aforementioned dysregulated immune cells through correlation analysis. The results showed that there was a significant correlation between the expression of DMGs and the abundance of these immune cells, for example, the expression of HK2 and LDHA was positively correlated with the abundance of Neutrophil, Monocyte, Eosinophil and Activated CD4 T cell, whereas RPS3 expression was negatively correlated with the abundance of Neutrophil and Central memory CD8 T cell (Figure 8C). The significant correlation suggests that the dysregulated expression of DMGs may contribute to the disruption of the immune microenvironment in DFU patients, thereby worsening the inflammatory features that characterize DFU patients.
Figure 8. Immune cell infiltration analysis. (A) Overall characterization of immune cell abundance. (B) The difference of immune cell infiltration abundance between the DFU samples and control samples. (C) Correlation between DMGs expression and immune cell infiltration abundance (*p < 0.05, **p < 0.01, ***p < 0.001).
3.9 Molecular docking
Previous studies has shown that drug discovery typically commences with identifying disease targets, with target-based drug discovery being the predominant approach for developing new drugs. Thus, targeting DMGs may provide a new effective therapeutic option for DFU treatment. Specifically, we identified two drugs capable of targeting DMGs based on the DSigDB database, namely, retinoic acid and estradiol. The results showed that both drugs were able to target DMGs and both had higher binding energies, sequentially LDHA-Retinoic acid (−7.6 kcal/mol) and LDHA-Estradiol (−8.5 kcal/mol) (Figures 9A, B), RPS3-Retinoic acid (−7.5 kcal/mol) and RPS3- Estradiol (−8.3 kcal/mol) (Figures 9C, D), HK2-Retinoic acid (−6.6 kcal/mol), and HK2-Estradiol (−7.2 kcal/mol) (Figures 9E, F). Notably, retinoic acid and estradiol have been shown in previous reports to be used in the treatment of skin-related diseases, for example, retinoic acid and its derivatives have potential for the treatment of severe skin diseases (Szymański et al., 2020). Estradiol promotes the production of extracellular matrix (ECM) (Soldano et al., 2010; Aida-Yasuoka et al., 2013; Baker Frost et al., 2021), which helps to enhance skin elasticity and reduce wrinkles (Rzepecki et al., 2019). These also provide the basis for the use of retinoic acid and estradiol in the treatment or alleviation of DFU symptoms.
Figure 9. Molecular docking simulations. (A,B) The structure of the complex formed by the docking of LDHA with Retinoic acid (A), Estradiol (B). (C,D) The structure of the complex formed by the docking of RPS3 with Retinoic acid (C), Estradiol (D). (E,F) The structure of the complex formed by the docking of HK2 with Retinoic acid (E), Estradiol (F).
4 Discussion
DFU is one of the primary causes of disability and death in DM patients. Previous studies have shown a strong relationship between DFU and mitophagy (Qi et al., 2023; Wang et al., 2023), and importantly mitophagy activators present potential as a novel therapeutic approach for diabetes (Shan et al., 2022), that is, through elimination of malfunctioning mitochondria, mitophagy helps to reduce the production of ROS and alleviate oxidative stress (Rovira-Llopis et al., 2017; Zhang et al., 2019), which is particularly important for the healing of DFUs. Therefore, the aim of the present study was to investigate the potential of MRGs as DFU biomarkers thereby contributing to the clinical diagnosis and treatment of DFU.
The current study was based on bulk RNA-seq data from DFU. Specifically, we firstly obtained 1,665 genes that were deregulated in DFU by differential analysis, of which 854 were up-regulated and 811 were down-regulated. Immediately after that, we used weighted gene co-expression network analysis to find that 42.2% of the mentioned above 1,665 DEGs, which are DFU-related DEGs, had their deregulated expression associated with the occurrence of DFU. Meanwhile, in order to investigate the role that mitophagy contributes to the occurrence of DFU, we isolated the MRGs in DFU-related DEGs using the MRGs annotated in the GeneCards database for MRGs with potential regulation in DFU, known as DFU-related MRGs. It is important to notice that the genes themselves do not operate and function independently, and their internal expression also has a strong correlation. The results of correlation analysis showed that almost all of the 22 DFU-related MRGs identified had significant correlation in their internal expression, implying that these genes, as the core genes of mitophagy factors, can largely contribute to the occurrence and development of DFU. Secondly, to verify the reliability of these genes, we enriched and analyzed them, and the results showed that they were significantly correlated with mitophagy, energy generation and metabolism, amino acid metabolism, and cancer, which means that our results have a certain reliability. In order to identify DMGs in DFU, we utilized multiple machine learning, and the results showed that HK2, RPS3 and LDHA could be indicated by all the machine learning algorithms.
It is worth mentioning that HK2, RPS3 and LDHA have been reported in previous studies. For example, Hexokinase, the first key enzyme in the glucose metabolic pathway responsible for the conversion of glucose to glucose 6-phosphate, is encoded by HK2. It was previously shown that HK2 expression is insulin-responsive, and in diabetes mellitus, HK2-associated unscheduled glycolysis may be a key initiator of insulin resistance and the development of vascular complications (Rabbani et al., 2022). It has also been demonstrated that inhibition of HK2 up-regulation can block excessive glycolytic responses, which may result in preservation of mitochondrial function and inhibition of diabetic renal fibrosis (Hasegawa et al., 2023). RPS3 is an essential component of the 40S subunit of eukaryotic ribosomes and is directly involved in ribosome maturation and translation initiation. RPS3 is expressed in human pancreatic islet cells in response to the transfer of cytokines from the cytoplasmic to the nuclear compartment. Furthermore, RPS3 is involved in NF-κB signaling, and a genetic increase in the activity of the NF-κB subunit c-Rel leads to protection against human islet cell death (Mokhtari et al., 2009). Meanwhile, ethanol may affect oxidative stress in β-cells by increasing the expression of the RPS3 gene, thereby decreasing their metabolic activity (Shin et al., 2004). Lactate dehydrogenase (LDH) is a metabolic enzyme that catalyzes the interconversion of pyruvate and lactate. LDH is a tetramer consisting of two subunits: muscular and cardiac, encoded by LDHA and LDHB, respectively (Zhang et al., 2022). LDHA catalyzes the reversible conversion of pyruvate to lactate while oxidizing NADH to NAD in anaerobic glycolysis. The inhibition of LDHA under normal physiological conditions and its inappropriate up-regulation in the diabetic milieu is well documented and is predominantly enriched in human islet α cells (Sanchez et al., 2021). In addition, the pathogenesis of diabetic kidney disease (DKD) is associated with LDHA-mediated lactic acidosis, which leads to mitochondrial abnormalities and renal fibrosis in DKD patients (Walbridge et al., 1987). In summary, the DMGs are associated with energy metabolism abnormalities.
Certainly, the performance of DMGs for DFU still needs both single-gene as well as multi-gene tests, including ROC and ANN (Mandair et al., 2023). ANN, as a form of artificial intelligence, has found extensive application in clinical medicine (Kufel et al., 2023; Mandair et al., 2023; Popovic et al., 2023). However, the current clinical methods for diagnosing DFU have limitations, particularly in terms of molecular-level diagnosis, which is often reliant on the clinical experience (Voelker, 2023). Thus, we constructed the ANN model based on DMGs with high prediction accuracy. The AUC values were all more than 0.8. Overall, the ANN constructed based on the transcriptome of DMGs can assist in the prediction of DFU.
Prior research has shown that drug discovery typically commences with identifying disease targets, with target-based drug discovery being the predominant approach for developing new drugs (Paananen and Fortino, 2020; Trajanoska et al., 2023). We identified two drugs capable of targeting DMGs based on the DSigDB database, namely, retinoic acid and estradiol. Retinoic acid and its derivatives have therapeutic potential for severe skin diseases (Szymański et al., 2020). Estradiol promotes extracellular matrix (ECM) production (Soldano et al., 2010; Aida-Yasuoka et al., 2013; Baker Frost et al., 2021), contributes to skin elasticity and reduces wrinkles (Rzepecki et al., 2019). In addition, the results of molecular docking showed a high binding affinity between retinoic acid and estradiol and DMGs, implying that they may be promising therapeutic agents against DFU.
There are also several limitations of this study that need to be clarified. Firstly, the number of DFU samples in the GEO database that could be incorporated into the present study was small, which also limited our possibility to perform additional cohort analyses and validation. Secondly, the clinical information of DFU samples was incomplete, which led to the lack of clinical phenotyping research on DFU in the present study. Finally, this study still needs more biological and clinical experiments for its subsequent analysis and validation.
5 Conclusion
In the present study, we performed a series of bioinformatic analyses based on DFU- related transcriptomics data. Specifically, we identified 702 DFU-related DEGs, including 22 MRGs, and the enrichment analysis results showed that these genes are related to mitochondria and energy metabolism. Immediately following this, we used multiple machine learning algorithms (RF, Lasso and SVM-RFE) to collectively identify HK2, RPS3 and LDHA to serve as DMGs. Meanwhile, we constructed a novel ANN model for DFU diagnosis, and the ROC curves of the model showed good performance. Additionally, the results of ssGSEA showed that DMGs could regulate the immune microenvironment of DFU patients. Finally, we found that retinoic acid and estradiol may be promising drugs against DFU. In conclusion, this study identified DMGs through multiple analytical techniques and verified their importance in DFU through multiple dimensions, providing promising targets for the clinical diagnosis and treatment of DFU.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author contributions
HG: Conceptualization, Data curation, Formal Analysis, Software, Writing–original draft. KX: Investigation, Methodology, Validation, Visualization, Writing–original draft. YZ: Data curation, Methodology, Writing–review and editing. JZ: Funding acquisition, Resources, Supervision, Writing–review and editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was supported by Chongqing Municipal Science and Health Joint Medical Research Project (2021MSXM267).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2024.1420136/full#supplementary-material
Abbreviations
DM, Diabetes mellitus; DFU, Diabetic foot ulcers; DFI, Diabetic foot infection; MRGs, Mitophagy-related genes; DMGs, Hub mitophagy-related genes of DFU; mtDNA, Mitochondrial DNA; ROS, Reactive oxygen species; PCA, Principal component analysis; DEGs, Differentially expressed genes; TOM, Topological overlap matrix; ME, Module eigengene; GS, Gene significance; MS, Module significance; GO, Gene ontology; KEGG, Kyoto encyclopedia of genes and genomes; SVM-RFE, Support vector machine recursive feature elimination; LASSO, Least absolute shrinkage and selection operator; RF, Random forest; ROC, Receiver operating characteristic; AUC, Area under roc curve; ANN, Artificial neural network; GSEA, Gene set enrichment analysis; ssGSEA, single sample gene set enrichment analysis; HK, Hexokinase; HK2, Hexokinase-2; RPS3, Small ribosomal subunit protein us3; LDHA, L-lactate dehydrogenase A chain; DKD, Diabetic kidney disease; ECM, Extracellular matrix.
References
Ahmed, R., Augustine, R., Chaudhry, M., Akhtar, U. A., Zahid, A. A., Tariq, M., et al. (2022). Nitric oxide-releasing biomaterials for promoting wound healing in impaired diabetic wounds: state of the art and recent trends. Biomed. Pharmacother. 149, 112707. doi:10.1016/j.biopha.2022.112707
Aida-Yasuoka, K., Peoples, C., Yasuoka, H., Hershberger, P., Thiel, K., Cauley, J. A., et al. (2013). Estradiol promotes the development of a fibrotic phenotype and is increased in the serum of patients with systemic sclerosis. Arthritis Res. Ther. 15, R10. doi:10.1186/ar4140
Armstrong, D. G., Boulton, A. J. M., and Bus, S. A. (2017). Diabetic foot ulcers and their recurrence. N. Engl. J. Med. 376, 2367–2375. doi:10.1056/NEJMra1615439
Armstrong, D. G., Tan, T. W., Boulton, A. J. M., and Bus, S. A. (2023). Diabetic foot ulcers: a review. Jama 330, 62–75. doi:10.1001/jama.2023.10578
Baker Frost, D., Savchenko, A., Ogunleye, A., Armstrong, M., and Feghali-Bostwick, C. (2021). Elucidating the cellular mechanism for E2-induced dermal fibrosis. Arthritis Res. Ther. 23, 68. doi:10.1186/s13075-021-02441-x
Bardill, J. R., Laughter, M. R., Stager, M., Liechty, K. W., Krebs, M. D., and Zgheib, C. (2022). Topical gel-based biomaterials for the treatment of diabetic foot ulcers. Acta Biomater. 138, 73–91. doi:10.1016/j.actbio.2021.10.045
Chemello, G., Salvatori, B., Morettini, M., and Tura, A. (2022). Artificial intelligence methodologies applied to technologies for screening, diagnosis and care of the diabetic foot: a narrative review. Biosens. (Basel) 12, 985. doi:10.3390/bios12110985
Chen, Y., Zhang, Y., Jiang, M., Ma, H., and Cai, Y. (2024). HMOX1 as a therapeutic target associated with diabetic foot ulcers based on single-cell analysis and machine learning. Int. Wound J. 21, e14815. doi:10.1111/iwj.14815
Coye, T. L., Crisologo, P. A., Suludere, M. A., Malone, M., Oz, O. K., and Lavery, L. A. (2024). The infected diabetic foot: modulation of traditional biomarkers for osteomyelitis diagnosis in the setting of diabetic foot infection and renal impairment. Int. Wound J. 21, e14770. doi:10.1111/iwj.14770
GBD 2019 Diseases and Injuries Collaborators (2020). Global burden of 369 diseases and injuries in 204 countries and territories, 1990-2019: a systematic analysis for the Global Burden of Disease Study 2019. Lancet 396, 1204–1222. doi:10.1016/S0140-6736(20)30925-9
Guan, H., Wang, Y., Niu, P., Zhang, Y., Zhang, Y., Miao, R., et al. (2024). The role of machine learning in advancing diabetic foot: a review. Front. Endocrinol. (Lausanne) 15, 1325434. doi:10.3389/fendo.2024.1325434
Hänzelmann, S., Castelo, R., and Guinney, J. (2013). GSVA: gene set variation analysis for microarray and RNA-seq data. BMC Bioinforma. 14, 7. doi:10.1186/1471-2105-14-7
Hart, T., Milner, R., and Cifu, A. (2017). Management of a diabetic foot. Jama 318, 1387–1388. doi:10.1001/jama.2017.11700
Hasegawa, K., Sakamaki, Y., Tamaki, M., and Wakino, S. (2023). PCK1 protects against mitoribosomal defects in diabetic nephropathy in mouse models. J. Am. Soc. Nephrol. 34, 1343–1365. doi:10.1681/ASN.0000000000000156
Haug, C. J., and Drazen, J. M. (2023). Artificial intelligence and machine learning in clinical medicine, 2023. N. Engl. J. Med. 388, 1201–1208. doi:10.1056/NEJMra2302038
Hettinger, G., Mitra, N., Thom, S. R., and Margolis, D. J. (2024). An improved clinical and genetics-based prediction model for diabetic foot ulcer healing. Adv. Wound Care (New Rochelle) 13, 281–290. doi:10.1089/wound.2023.0194
Huang, X., Liang, P., Jiang, B., Zhang, P., Yu, W., Duan, M., et al. (2020). Hyperbaric oxygen potentiates diabetic wound healing by promoting fibroblast cell proliferation and endothelial cell angiogenesis. Life Sci. 259, 118246. doi:10.1016/j.lfs.2020.118246
Jeffcoate, W., Boyko, E. J., Game, F., Cowled, P., Senneville, E., and Fitridge, R. (2024). Causes, prevention, and management of diabetes-related foot ulcers. Lancet Diabetes Endocrinol. 12, 472–482. doi:10.1016/S2213-8587(24)00110-4
Jude, E. B., Oyibo, S. O., Chalmers, N., and Boulton, A. J. (2001). Peripheral arterial disease in diabetic and nondiabetic patients: a comparison of severity and outcome. Diabetes Care 24, 1433–1437. doi:10.2337/diacare.24.8.1433
Khandakar, A., Chowdhury, M. E. H., Ibne Reaz, M. B., Md Ali, S. H., Hasan, M. A., Kiranyaz, S., et al. (2021). A machine learning model for early detection of diabetic foot using thermogram images. Comput. Biol. Med. 137, 104838. doi:10.1016/j.compbiomed.2021.104838
Kufel, J., Bargieł-Łączek, K., Kocot, S., Koźlik, M., Bartnikowska, W., Janik, M., et al. (2023). What is machine learning, artificial neural networks and deep learning? examples of practical applications in medicine. Diagn. (Basel) 13, 2582. doi:10.3390/diagnostics13152582
Langfelder, P., and Horvath, S. (2008). WGCNA: an R package for weighted correlation network analysis. BMC Bioinforma. 9, 559. doi:10.1186/1471-2105-9-559
Le, N. Q. K., Kha, Q. H., Nguyen, V. H., Chen, Y. C., Cheng, S. J., and Chen, C. Y. (2021). Machine learning-based radiomics signatures for EGFR and KRAS mutations prediction in non-small-cell lung cancer. Int. J. Mol. Sci. 22, 9254. doi:10.3390/ijms22179254
Leek, J. T., Johnson, W. E., Parker, H. S., Jaffe, A. E., and Storey, J. D. (2012). The sva package for removing batch effects and other unwanted variation in high-throughput experiments. Bioinformatics 28, 882–883. doi:10.1093/bioinformatics/bts034
Li A., A., Gao, M., Liu, B., Qin, Y., Chen, L., Liu, H., et al. (2022). Mitochondrial autophagy: molecular mechanisms and implications for cardiovascular disease. Cell Death Dis. 13, 444. doi:10.1038/s41419-022-04906-6
Li, X., Xie, X., Lian, W., Shi, R., Han, S., Zhang, H., et al. (2018). Exosomes from adipose-derived stem cells overexpressing Nrf2 accelerate cutaneous wound healing by promoting vascularization in a diabetic foot ulcer rat model. Exp. Mol. Med. 50, 29–14. doi:10.1038/s12276-018-0058-5
Li Y., Y., Ju, S., Li, X., Li, W., Zhou, S., Wang, G., et al. (2022). Characterization of the microenvironment of diabetic foot ulcers and potential drug identification based on scRNA-seq. Front. Endocrinol. (Lausanne) 13, 997880. doi:10.3389/fendo.2022.997880
Lu, Y., Li, Z., Zhang, S., Zhang, T., Liu, Y., and Zhang, L. (2023). Cellular mitophagy: mechanism, roles in diseases and small molecule pharmacological regulation. Theranostics 13, 736–766. doi:10.7150/thno.79876
Mandair, D., Reis-Filho, J. S., and Ashworth, A. (2023). Biological insights and novel biomarker discovery through deep learning approaches in breast cancer histopathology. NPJ Breast Cancer 9, 21. doi:10.1038/s41523-023-00518-1
Matos, M., Mendes, R., Silva, A. B., and Sousa, N. (2018). Physical activity and exercise on diabetic foot related outcomes: a systematic review. Diabetes Res. Clin. Pract. 139, 81–90. doi:10.1016/j.diabres.2018.02.020
McDermott, K., Fang, M., Boulton, A. J. M., Selvin, E., and Hicks, C. W. (2023). Etiology, epidemiology, and disparities in the burden of diabetic foot ulcers. Diabetes Care 46, 209–221. doi:10.2337/dci22-0043
McIllhatton, A., Lanting, S., Lambkin, D., Leigh, L., Casey, S., and Chuter, V. (2021). Reliability of recommended non-invasive chairside screening tests for diabetes-related peripheral neuropathy: a systematic review with meta-analyses. BMJ Open Diabetes Res. Care 9, e002528. doi:10.1136/bmjdrc-2021-002528
Mokhtari, D., Barbu, A., Mehmeti, I., Vercamer, C., and Welsh, N. (2009). Overexpression of the nuclear factor-κB subunit c-Rel protects against human islet cell death in vitro. Am. J. Physiol. Endocrinol. Metab. 297, E1067–E1077. doi:10.1152/ajpendo.00212.2009
Noble, W. S. (2006). What is a support vector machine? Nat. Biotechnol. 24, 1565–1567. doi:10.1038/nbt1206-1565
Onishi, M., Yamano, K., Sato, M., Matsuda, N., and Okamoto, K. (2021). Molecular mechanisms and physiological functions of mitophagy. Embo J. 40, e104705. doi:10.15252/embj.2020104705
Orlando, G., Prior, Y., Reeves, N. D., and Vileikyte, L. (2021). Patient and provider perspective of smart wearable technology in diabetic foot ulcer prevention: a systematic review. Med. Kaunas. 57, 1359. doi:10.3390/medicina57121359
Paananen, J., and Fortino, V. (2020). An omics perspective on drug target discovery platforms. Brief. Bioinform 21, 1937–1953. doi:10.1093/bib/bbz122
Paul, A., Mukherjee, D. P., Das, P., Gangopadhyay, A., Chintha, A. R., and Kundu, S. (2018). Improved random forest for classification. IEEE Trans. Image Process 27, 4012–4024. doi:10.1109/TIP.2018.2834830
Pitocco, D., Spanu, T., Di Leo, M., Vitiello, R., Rizzi, A., Tartaglione, L., et al. (2019). Diabetic foot infections: a comprehensive overview. Eur. Rev. Med. Pharmacol. Sci. 23, 26–37. doi:10.26355/eurrev_201904_17471
Popovic, D., Glisic, T., Milosavljevic, T., Panic, N., Marjanovic-Haljilji, M., Mijac, D., et al. (2023). The importance of artificial intelligence in upper gastrointestinal endoscopy. Diagn. (Basel) 13, 2862. doi:10.3390/diagnostics13182862
Primadhi, R. A., Septrina, R., Hapsari, P., and Kusumawati, M. (2023). Amputation in diabetic foot ulcer: a treatment dilemma. World J. Orthop. 14, 312–318. doi:10.5312/wjo.v14.i5.312
Purwanti, O. S., Nursalam, N., and Pandin, M. G. R. (2024). Early detection of diabetic neuropathy based on health belief model: a scoping review. Front. Endocrinol. (Lausanne) 15, 1369699. doi:10.3389/fendo.2024.1369699
Qi, X., Cai, E., Xiang, Y., Zhang, C., Ge, X., Wang, J., et al. (2023). An immunomodulatory hydrogel by hyperthermia-assisted self-cascade glucose depletion and ROS scavenging for diabetic foot ulcer wound therapeutics. Adv. Mater 35, e2306632. doi:10.1002/adma.202306632
Rabbani, N., Xue, M., and Thornalley, P. J. (2022). Hexokinase-2-Linked glycolytic overload and unscheduled glycolysis-driver of insulin resistance and development of vascular complications of diabetes. Int. J. Mol. Sci. 23, 2165. doi:10.3390/ijms23042165
Rendra, E., Riabov, V., Mossel, D. M., Sevastyanova, T., Harmsen, M. C., and Kzhyshkowska, J. (2019). Reactive oxygen species (ROS) in macrophage activation and function in diabetes. Immunobiology 224, 242–253. doi:10.1016/j.imbio.2018.11.010
Rovira-Llopis, S., Bañuls, C., Diaz-Morales, N., Hernandez-Mijares, A., Rocha, M., and Victor, V. M. (2017). Mitochondrial dynamics in type 2 diabetes: pathophysiological implications. Redox Biol. 11, 637–645. doi:10.1016/j.redox.2017.01.013
Ruder, K. (2024). Diabetic foot infections and amputations are all too common-here’s what could move the needle. Jama 331, 998–1000. doi:10.1001/jama.2024.1016
Rzepecki, A. K., Murase, J. E., Juran, R., Fabi, S. G., and McLellan, B. N. (2019). Estrogen-deficient skin: the role of topical therapy. Int. J. Womens Dermatol 5, 85–90. doi:10.1016/j.ijwd.2019.01.001
Sanchez, P. K. M., Khazaei, M., Gatineau, E., Geravandi, S., Lupse, B., Liu, H., et al. (2021). LDHA is enriched in human islet alpha cells and upregulated in type 2 diabetes. Biochem. Biophys. Res. Commun. 568, 158–166. doi:10.1016/j.bbrc.2021.06.065
Senneville, É., Albalawi, Z., van Asten, S. A., Abbas, Z. G., Allison, G., Aragón-Sánchez, J., et al. (2024). IWGDF/IDSA guidelines on the diagnosis and treatment of diabetes-related foot infections (IWGDF/IDSA 2023). Diabetes Metab. Res. Rev. 40, e3687. doi:10.1002/dmrr.3687
Shan, Z., Fa, W. H., Tian, C. R., Yuan, C. S., and Jie, N. (2022). Mitophagy and mitochondrial dynamics in type 2 diabetes mellitus treatment. Aging (Albany NY) 14, 2902–2919. doi:10.18632/aging.203969
Shi, H., Yuan, X., Liu, G., and Fan, W. (2023). Identifying and validating GSTM5 as an immunogenic gene in diabetic foot ulcer using bioinformatics and machine learning. J. Inflamm. Res. 16, 6241–6256. doi:10.2147/JIR.S442388
Shi, H., Yuan, X., Yang, X., Huang, R., Fan, W., and Liu, G. (2024). A novel diabetic foot ulcer diagnostic model: identification and analysis of genes related to glutamine metabolism and immune infiltration. BMC Genomics 25, 125. doi:10.1186/s12864-024-10038-2
Shibue, K. (2023). Artificial intelligence and machine learning in clinical medicine. N. Engl. J. Med. 388, 2398.
Shin, J. S., Kwon, Y. S., Lee, J. J., and Kim, C. W. (2004). Isolation of ethanol-induced genes in pancreatic beta-cells by representational difference analysis (RDA). Exp. Mol. Med. 36, 36–42. doi:10.1038/emm.2004.5
Slomski, A. (2022). Intensive glycemic treatment reduces diabetic foot ulcers. Jama 327, 1118. doi:10.1001/jama.2022.3754
Soldano, S., Montagna, P., Brizzolara, R., Sulli, A., Parodi, A., Seriolo, B., et al. (2010). Effects of estrogens on extracellular matrix synthesis in cultures of human normal and scleroderma skin fibroblasts. Ann. N. Y. Acad. Sci. 1193, 25–29. doi:10.1111/j.1749-6632.2009.05296.x
Szymański, Ł., Skopek, R., Palusińska, M., Schenk, T., Stengel, S., Lewicki, S., et al. (2020). Retinoic acid and its derivatives in skin. Cells 9, 2660. doi:10.3390/cells9122660
Trajanoska, K., Bhérer, C., Taliun, D., Zhou, S., Richards, J. B., and Mooser, V. (2023). From target discovery to clinical drug development with human genetics. Nature 620, 737–745. doi:10.1038/s41586-023-06388-8
Vasquez, M. M., Hu, C., Roe, D. J., Chen, Z., Halonen, M., and Guerra, S. (2016). Least absolute shrinkage and selection operator type methods for the identification of serum biomarkers of overweight and obesity: simulation and application. BMC Med. Res. Methodol. 16, 154. doi:10.1186/s12874-016-0254-8
Villa, F., Marchandin, H., Lavigne, J. P., Schuldiner, S., Cellier, N., Sotto, A., et al. (2024). Anaerobes in diabetic foot infections: pathophysiology, epidemiology, virulence, and management. Clin. Microbiol. Rev., e0014323. doi:10.1128/cmr.00143-23
Walbridge, D. G., Knutson, J. R., and Brand, L. (1987). Nanosecond time-resolved fluorescence measurements during protein denaturation. Anal. Biochem. 161, 467–478. doi:10.1016/0003-2697(87)90476-3
Wang, B., Du, L., Dong, B., Kou, E., Wang, L., and Zhu, Y. (2024). Current knowledge and perspectives of phage therapy for combating refractory wound infections. Int. J. Mol. Sci. 25, 5465. doi:10.3390/ijms25105465
Wang, F., Zhang, X., Zhang, J., Xu, Q., Yu, X., Xu, A., et al. (2024). Recent advances in the adjunctive management of diabetic foot ulcer: focus on noninvasive technologies. Med. Res. Rev. 44, 1501–1544. doi:10.1002/med.22020
Wang, K., Wang, Y., Shi, W., Shen, K., Tao, K., Ling, R., et al. (2024). Diagnosis and treatment of diabetic foot ulcer complicated with lower extremity vasculopathy: consensus recommendation from the Chinese Medical Association (CMA), Chinese Medical Doctor Association (CMDA). Diabetes Metab. Res. Rev. 40, e3776. doi:10.1002/dmrr.3776
Wang, S., Long, H., Hou, L., Feng, B., Ma, Z., Wu, Y., et al. (2023). The mitophagy pathway and its implications in human diseases. Signal Transduct. Target Ther. 8, 304. doi:10.1038/s41392-023-01503-7
Wang, Y., Shao, T., Wang, J., Huang, X., Deng, X., Cao, Y., et al. (2021). An update on potential biomarkers for diagnosing diabetic foot ulcer at early stage. Biomed. Pharmacother. 133, 110991. doi:10.1016/j.biopha.2020.110991
Wei, T., Pan, T., Peng, X., Zhang, M., Guo, R., Guo, Y., et al. (2024). Janus liposozyme for the modulation of redox and immune homeostasis in infected diabetic wounds. Nat. Nanotechnol. doi:10.1038/s41565-024-01660-y
Wronka, M., Krzemińska, J., Młynarska, E., Rysz, J., and Franczyk, B. (2022). The influence of lifestyle and treatment on oxidative stress and inflammation in diabetes. Int. J. Mol. Sci. 23, 15743. doi:10.3390/ijms232415743
Wu, Z., Liu, P., Huang, B., Deng, S., Song, Z., Huang, X., et al. (2023). A novel Alzheimer’s disease prognostic signature: identification and analysis of glutamine metabolism genes in immunogenicity and immunotherapy efficacy. Sci. Rep. 13, 6895. doi:10.1038/s41598-023-33277-x
Xiong, Y., Chen, L., Yu, T., Yan, C., Zhou, W., Cao, F., et al. (2020). Inhibition of circulating exosomal microRNA-15a-3p accelerates diabetic wound repair. Aging (Albany NY) 12, 8968–8986. doi:10.18632/aging.103143
Yi, W. J., Yuan, Y., Bao, Q., Zhao, Z., Ding, H. S., and Song, J. (2024). Analyzing immune cell infiltration and copper metabolism in diabetic foot ulcers. J. Inflamm. Res. 17, 3143–3157. doi:10.2147/JIR.S452609
Yu, G., Wang, L. G., Han, Y., and He, Q. Y. (2012). clusterProfiler: an R package for comparing biological themes among gene clusters. Omics 16, 284–287. doi:10.1089/omi.2011.0118
Zhang G., G., Samarawickrama, P. N., Gui, L., Ma, Y., Cao, M., Zhu, H., et al. (2024). Revolutionizing diabetic foot ulcer care: the senotherapeutic approach. Aging Dis., 0. doi:10.14336/AD.2024.0065
Zhang N., N., Liu, Y., Yan, W., and Liu, F. (2024). The effect of negative pressure wound therapy on the outcome of diabetic foot ulcers: a meta-analysis. Int. Wound J. 21, e14886. doi:10.1111/iwj.14886
Zhang, W., Wang, C., Hu, X., Lian, Y., Ding, C., and Ming, L. (2022). Inhibition of LDHA suppresses cell proliferation and increases mitochondrial apoptosis via the JNK signaling pathway in cervical cancer cells. Oncol. Rep. 47, 77. doi:10.3892/or.2022.8288
Keywords: diabetic foot ulcers, mitophagy-related genes, machine learning algorithms, artificial neural network, molecular docking method
Citation: Guo H, Xiao K, Zheng Y and Zong J (2024) Integrating bioinformatics and multiple machine learning to identify mitophagy-related targets for the diagnosis and treatment of diabetic foot ulcers: evidence from transcriptome analysis and drug docking. Front. Mol. Biosci. 11:1420136. doi: 10.3389/fmolb.2024.1420136
Received: 19 April 2024; Accepted: 20 June 2024;
Published: 09 July 2024.
Edited by:
Raffaele Capasso, University of Naples Federico II, ItalyReviewed by:
Marta Laranjo, University of Evora, PortugalAndaç Batur Çolak, Istanbul Commerce University, Türkiye
Copyright © 2024 Guo, Xiao, Zheng and Zong. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jianchun Zong, zongjianchuncq@163.com
†These authors have contributed equally to this work