Skip to main content
Download PDF

Interpretable machine learning models for detecting peripheral neuropathy and lower extremity arterial disease in diabetics: an analysis of critical shared and unique risk factors

BMC Medical Informatics and Decision Making volume 24, Article number: 200 (2024) Cite this article

Abstract

Background

Diabetic peripheral neuropathy (DPN) and lower extremity arterial disease (LEAD) are significant contributors to diabetic foot ulcers (DFUs), which severely affect patients’ quality of life. This study aimed to develop machine learning (ML) predictive models for DPN and LEAD and to identify both shared and distinct risk factors.

Methods

This retrospective study included 479 diabetic inpatients, of whom 215 were diagnosed with DPN and 69 with LEAD. Clinical data and laboratory results were collected for each patient. Feature selection was performed using three methods: mutual information (MI), random forest recursive feature elimination (RF-RFE), and the Boruta algorithm to identify the most important features. Predictive models were developed using logistic regression (LR), random forest (RF), and eXtreme Gradient Boosting (XGBoost), with particle swarm optimization (PSO) used to optimize their hyperparameters. The SHapley Additive exPlanation (SHAP) method was applied to determine the importance of risk factors in the top-performing models.

Results

For diagnosing DPN, the XGBoost model was most effective, achieving a recall of 83.7%, specificity of 86.8%, accuracy of 85.4%, and an F1 score of 83.7%. On the other hand, the RF model excelled in diagnosing LEAD, with a recall of 85.7%, specificity of 92.9%, accuracy of 91.9%, and an F1 score of 82.8%. SHAP analysis revealed top five critical risk factors shared by DPN and LEAD, including increased urinary albumin-to-creatinine ratio (UACR), glycosylated hemoglobin (HbA1c), serum creatinine (Scr), older age, and carotid stenosis. Additionally, distinct risk factors were pinpointed: decreased serum albumin and lower lymphocyte count were linked to DPN, while elevated neutrophil-to-lymphocyte ratio (NLR) and higher D-dimer levels were associated with LEAD.

Conclusions

This study demonstrated the effectiveness of ML models in predicting DPN and LEAD in diabetic patients and identified significant risk factors. Focusing on shared risk factors may greatly reduce the prevalence of both conditions, thereby mitigating the risk of developing DFUs.

Peer Review reports

Introduction

Diabetes has become a growing global health concern, affecting around 451 million people in 2017, with projections to rise to 693 million by 2045 [1]. In China, the prevalence of diabetes has surged from less than 1% in the 1980s to approximately 10.9% in 2013 and 12.4% in 2018, making it the country with the world’s largest diabetic population [2]. Diabetic peripheral neuropathy (DPN) and lower extremity arterial disease (LEAD) are prevalent complications of diabetes, with occurrence rates of about 50% and 3-20%, respectively [3,4,5,6]. Both complications serve as extrinsic risk factors for diabetic foot ulcers (DFUs) [7], leading to higher rates of amputation, increased mortality, and substantial economic burdens for patients with diabetes [8]. Unfortunately, patients with DPN or LEAD may be asymptomatic in their early stages, and many patients already have these complications at the time of initial diagnosis [8, 9]. Therefore, early identification and management of DPN and LEAD are crucial in preventing DFUs among diabetic patients.

At present, the diagnosis of DPN and LEAD mainly relies on physical examination of the peripheral nervous system, electromyography (EMG), ankle–brachial index, lower limb vascular ultrasound, etc. [5, 10]. However, these methods require well-trained endocrinologists and specialized diagnostic equipment, which are often scarce in underdeveloped regions. To address this challenge, researchers are actively exploring the development of practical, accessible, and cost-effective clinical diagnosis models for DPN and LEAD based on clinical features and routinely measured lab parameters. Recent studies have demonstrated that machine learning (ML) models, by utilizing medical history, physical examinations, and basic lab tests, could effectively predict DPN and LEAD [4, 11]. Moreover, ML algorithms achieved high accuracy in identifying DPN through the analysis of immune biomarkers or microcirculatory parameters [12, 13]. Additionally, a model based on support vector machine (SVM) has been shown to accurately predict DPN severity in about 76% of cases, utilizing general patient information and responses from a neuropathy disability score questionnaire [14]. Despite these advances, most studies have concentrated on developing one model for either DPN or LEAD, without considering common risk factors for both conditions in a single study. Given the high prevalence of DPN and LEAD in developing countries and their potential to lead to DFUs, which can significantly increase mortality rates [15, 16], it is crucial to identify and target shared risk factors for both conditions. This strategy could facilitate early and concurrent interventions, potentially diminishing the prevalence and severity of these diseases.

In this study, we employed logistic regression (LR), random forest (RF), and eXtreme Gradient Boosting (XGBoost) to develop diagnostic models for both DPN and LEAD among diabetic individuals, utilizing demographic, clinical, and laboratory information. This research spanned the interdisciplinary fields of medicine, biostatistics, and ML. To identify shared and unique risk factors for DPN and LEAD, we used the SHapley Additive exPlanation (SHAP) method to prioritize risk factors within the most effective models for each condition.

Contributions of this work

The major contributions of this study are outlined as follows:

  1. (1)

    We constructed ML models for DPN and LEAD detection based on accessible demographic, clinical, and laboratory data, minimizing the need for specialized tests and advanced medical facilities. This approach is especially beneficial for areas with limited healthcare resources.

  2. (2)

    We utilized three feature selection methods—mutual information (MI), random forest recursive feature elimination (RF-RFE), and the Boruta algorithm—to identify the most significant features. This strategy effectively reduces overfitting and enhances the robustness of our models.

  3. (3)

    To optimize the performance of each ML model, we applied particle swarm optimization (PSO) to fine-tune hyperparameters.

  4. (4)

    The SHAP method was applied to elucidate the contribution of each feature to the risk of developing DPN and LEAD in the best-performing models. This analysis identified both shared and distinct risk factors for DPN and LEAD, deepening our insight into their pathophysiological foundations. Concentrating on shared risk factors may significantly reduce the prevalence of these conditions and subsequently the risk of DFUs.

Methods

Study design and participants

This is a cross-sectional study conducted at Tongji Hospital from January 2022 to March 2023. We collected clinical characteristics and laboratory data of 712 diabetic inpatients who underwent EMG and lower limb vascular ultrasound examinations. Patients were excluded from the study if they had diabetic ketoacidosis (DKA), hyperosmotic hyperglycemia syndrome (HHS), autoimmune diseases, infectious diseases, malignant tumors, or if more than 30% of their data was missing. Ultimately, 479 diabetic patients were enrolled in this study. This research was performed in compliance with the Code of Ethics of the World Medical Association (Declaration of Helsinki) and received approval from the Institutional Ethics Committee (K-2023-022).

Diagnostic criteria

All diabetic patients enrolled in our study underwent screening for DPN and LEAD during hospitalization. The diagnosis of diabetic complications was made by two qualified endocrinologists based on the locally recognized criteria [17]. A Dantec® Keypoint® G4 EMG/NCS/EP Workstation with 8-Channel Amplifier (Dantec Medical A/S, Denmark) was used to test the motor conduction velocity of bilateral ulnar, median, and common peroneal nerves, as well as the sensory conduction velocity of bilateral radial, median, and superficial fibular nerves. The diagnosis of DPN relied on typical symptoms, neurological examination, EMG, and exclusion of other causes of peripheral neuropathy. In addition, Philips IU22 Doppler ultrasonic color imaging system (Philips, USA) equipped with a 3-D array probe (5-12 MHz) was applied to examine the bilateral common femoral arteries, superficial femoral artery, popliteal artery, and dorsal artery. The diagnosis of LEAD was based on arterial lumen stenosis, severe blood flow filling deficiency, or arterial occlusion.

Data collection

Demographic information was collected on admission, including age, gender, body mass index (BMI), diabetes duration, diabetes type, smoking history, systolic blood pressure (SBP) and diastolic blood pressure (DBP). BMI was calculated by dividing an individual’s weight in kilograms by the square of their height in meters. SBP and DBP were measured three times while the person was at rest, and then the average of the readings was recorded. Clinical routine examinations were performed, and laboratory parameters were obtained from fasting blood and spot urine samples collected the next morning after admission. The routine laboratory tests included hematology (neutrophil, lymphocyte, and platelet counts, neutrophil-to-lymphocyte ratio (NLR)), liver and kidney function tests (alanine aminotransferase (ALT), aspartate aminotransferase (AST), serum albumin, total bilirubin (TBiL), direct bilirubin (DBiL), serum urea nitrogen (SUN), serum uric acid (SUA), serum creatinine (Scr)), glucose metabolism (fasting blood glucose (FBG), glycosylated hemoglobin (HbA1c)), lipid profiles (total cholesterol (TC), high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), and triglycerides (TG)), islet function (fasting insulin and fasting C-peptide), C-reactive protein (CRP), D-dimer, 25-hydroxy vitamin D3 (25-OH VitD), ferritin, neuron specific enolase (NSE), urinary microalbumin and creatinine, and the urinary albumin-to-creatinine ratio (UACR). We calculated the estimated glomerular filtration rate (eGFR) using a formula provided in a previous study [18]. In addition, we collected the results of carotid artery ultrasound. Bilateral common carotid artery, internal carotid artery, and external carotid artery were examined using a Philips IU22 Doppler ultrasonic color imaging system (Philips, USA) equipped with a 3-D array probe (7-12 MHz). The dataset, encompassing all features and their respective values, was detailed in Supplementary Table S1.

Missing data

Comprehensive demographic information was available for all participants, as each patient underwent the hospitalization process. Missing data for laboratory parameters were below 30%. We addressed these missing values using the most recent available measurements. Any remaining missing values were imputed using the median.

Data balancing

In constructing models for LEAD, we encountered a class imbalance issue due to the low proportion of patients with LEAD in the overall population. To address this problem, we used the imbalanced-learn package in Python to employ a random undersampling technique, achieving a 1:3 ratio between the LEAD and non-LEAD groups [19]. This approach helped us to maintain a more balanced ratio and decrease the number of non-LEAD cases to three times the number of LEAD cases, thus mitigating the imbalance and enhancing the reliability of our models.

Feature selection strategy

The feature selection process was conducted using three distinct methods: MI, RF-RFE, and the Boruta algorithm. MI quantifies the dependency between variables by capturing all types of relationships, both linear and nonlinear. For feature selection, MI assesses the dependency of each feature on the target label to identify the most informative features for prediction [20]. RF-RFE utilizes a RF to iteratively build models, systematically removing the least important features in each round. This method emphasizes features that significantly affect model performance [21]. The Boruta algorithm employs a RF classifier to evaluate features against their randomized “shadow” versions, ensuring only essential features are retained for accurate model predictions [22].

For both MI and RF-RFE, the top 15 features were identified independently. The Boruta algorithm categorized features as confirmed, tentative, or rejected, selecting features that were either confirmed or tentative. Only features chosen by at least two of these three methods were used to develop ML models. This approach reduces redundancy and enhances the predictive accuracy of the ML models.

ML model construction and interpretation

The model was constructed and interpreted using Python (version 3.9.6, Python Software Foundation, USA). The workflow for constructing and interpreting ML models is illustrated in Fig. 1. First, the dataset was randomly divided into two subsets: 80% designated for training the model and the remaining 20% reserved for testing. Then, three distinct ML models—LR, RF, and XGBoost—were developed to predict DPN and LEAD based on selected features. To optimize these models and select the most suitable hyperparameters, we employed PSO. PSO is a computational method that mimics a swarm of particles navigating through the parameter space to find optimal solutions. In ML, PSO enhances model parameterization by representing each particle as a potential solution that is continuously refined through both individual and collective experiences within the swarm. This strategy efficiently identifies the best parameter combinations, significantly improving model performance [23]. The set and optimal value of hyperparameters were displayed in Supplementary Table S2.

Fig. 1
figure 1

ML model development and evaluation process

Full size image

The effectiveness of each model was evaluated using various metrics, including recall, specificity, precision, accuracy, and the F1 score. We also calculated the area under the receiver operating curve (AUC) for the test sets to assess the performance of each model. Furthermore, the SHAP method was employed to interpret the contribution of each predictor within the optimal models. Through SHAP analysis, we gained a detailed understanding of how each feature influences the model’s output, providing a comprehensive insight into the model’s decision-making process [24].

Statistical analysis

The statistical analyses were performed using SPSS (version 27.0, IBM, USA). For data adhering to a normal distribution, values were depicted as mean ± standard deviation. Differences among these values were examined using the independent Student’s t-test. Conversely, for data not following a normal distribution, variables were presented as medians (interquartile range, IQR), and the Mann-Whitney U test was employed to evaluate disparities in their distributions. Categorical data were represented as n (%) and analyzed for distribution differences via the Chi-square (χ2) test or Fisher’s exact test, as appropriate. A p-value < 0.05 was considered statistically significant.

Results

Clinical features of patients

In this study, we initially enrolled 712 diabetic inpatients. After applying exclusion criteria, 479 patients qualified for inclusion. Among them, 215 were diagnosed with DPN, and 69 with LEAD. The median age of participants was 50 years (IQR: 48-56), and the male-to-female ratio was 0.58. All these cases were utilized to develop models for diagnosing DPN. To correct for the imbalance in sample sizes, a one-to-three random undersampling strategy was employed for LEAD cases versus non-LEAD controls, resulting in 69 LEAD cases and 207 non-LEAD cases being selected to construct LEAD prediction models (Fig. 2). According to univariate analysis, out of the 38 features, 24 exhibited significant discrepancies between patients with and without DPN, whereas 17 features displayed differences between those with and without LEAD (Tables 1 and 2). Patients with DPN or LEAD were found to be older compared to those without these complications. Men were more likely to develop DPN or LEAD than women. Moreover, patients with DPN and LEAD exhibited increased levels of HbA1c, SUN, Scr, FBG, D-dimer, NLR, urinary microalbumin, and UACR compared to those without these complications. In contrast, the levels of serum albumin, eGFR, TC, and LDL were found to be lower in patients with DPN and LEAD. Additionally, a positive association was observed between the presence of carotid stenosis and the occurrence of DPN and LEAD.

Fig. 2
figure 2

Flowchart of patient enrollment. DPN, diabetic peripheral neuropathy. LEAD, lower extremity arterial disease. DKA, diabetic ketoacidosis. HHS, hyperosmotic hyperglycemia syndrome. ML, machine learning

Full size image
Table 1 Clinical features in diabetic patients with or without DPN
Full size table
Table 2 Clinical features in diabetic patients with or without LEAD
Full size table

Selected features

For DPN, a consensus was reached on eight features selected by all three feature selection methods. An additional four features were agreed upon by two of the methods, resulting in a total of 12 distinct features that were incorporated into the models, as outlined in Supplementary Table S3. Similarly, for LEAD, unanimous selection was achieved for five features across all methods, with another four features chosen by two of the methods. Thus, a total of nine features were integrated into the ML models, as detailed in Supplementary Table S4.

Diagnostic performance of LR, RF, and XGBoost in detecting DPN

The diagnostic performances of three models for detecting DPN were shown in Figs. 3A and 4A. Among these models, XGBoost demonstrated the highest diagnostic efficacy, achieving an AUC of 0.903, a recall of 83.7%, a specificity of 86.8%, an accuracy of 85.4%, a precision of 83.7%, and an F1 score of 83.7%. Additionally, RF showed the highest specificity, at 90.6%.

Fig. 3
figure 3

ROC curves of LR, RF, and XGBoost models for detecting DPN and LEAD. A ROC curves for DPN; B ROC curves for LEAD. ROC, receiver operating characteristic. LR, logistic regression. RF, random forest. DPN, diabetic peripheral neuropathy. LEAD, lower extremity arterial disease

Full size image
Fig. 4
figure 4

Performance of LR, RF, and XGBoost for detecting DPN and LEAD. A Performance for DPN models; B Performance for LEAD models. LR, logistic regression. RF, random forest. DPN, diabetic peripheral neuropathy. LEAD, lower extremity arterial disease

Full size image

Diagnostic performance of LR, RF, and XGBoost in detecting LEAD

The performances of the LEAD models were presented in Figs. 3B and 4B. The RF model outperformed the others with the highest AUC of 0.923, recall of 85.7%, specificity of 92.9%, accuracy of 91.9%, precision of 80.0%, and an F1 score of 82.8%, followed by XGBoost and LR.

Critical shared and unique risk factors for DPN and LEAD through SHAP analysis

SHAP was applied to evaluate the importance of features within the optimal ML models for DPN and LEAD, with a prioritized list vividly illustrating their respective impacts. Figure 5A and B presented the rankings of critical features in the XGBoost model for DPN and the RF model for LEAD, respectively. This analysis highlighted the importance of both shared and unique risk factors. Common risk factors identified for both conditions include increased UACR, elevated HbA1c, elevated Scr, advanced age, carotid stenosis, high FBG, and reduced eGFR. Unique to DPN were decreased serum albumin and lower lymphocyte count, whereas LEAD was specifically associated with increased NLR and higher D-dimer levels.

Fig. 5
figure 5

Feature importance of SHAP values for XGBoost model in detecting DPN and for RF model in detecting LEAD. A SHAP values of XGBoost model in detecting DPN; B SHAP values of RF model in detecting LEAD. SHAP, SHapley Additive exPlanation. RF, random forest. DPN, diabetic peripheral neuropathy. LEAD, lower extremity arterial disease

Full size image

Discussion

This study constructed three different ML models for predicting DPN and LEAD among diabetic patients, utilizing basic clinical and laboratory data. We discovered that the XGBoost model demonstrated superior diagnostic performance in detecting DPN, whereas the RF model excelled in identifying LEAD. Furthermore, SHAP analysis identified the top five important risk factors common to both conditions: elevated UACR, HbA1c, Scr, advanced age, and carotid stenosis. Additionally, it pinpointed unique risk factors for each condition: a decrease in serum albumin and lymphocyte count were significant for DPN, while increased NLR and D-dimer were key indicators for LEAD.

ML models are significantly advancing the field of medical diagnostics. Recent advancements in predicting DPN and LEAD were summarized Table 3. For DPN detection, Metsker et al. [25] developed ML models using age, gender, and 27 laboratory tests. Among these models, the artificial neural network (ANN) achieved the highest recall at 0.809, the LR had the highest precision at 0.683, while Linear Regression displayed both the highest F1 score at 0.730 and the highest accuracy at 0.747. Another study demonstrated that, using demographic, clinical, and laboratory data, both RF and SVM models significantly distinguished DPN in individuals with T2DM. The accuracy, sensitivity, and specificity were 67.8%, 68.09%, and 67.44% for RF, and 67.8%, 68.89%, and 66.67% for SVM, respectively [26]. By contrast, our study showed that the XGBoost model had the highest diagnostic performance, with an accuracy of 85.4%, a sensitivity of 83.7%, and a specificity of 86.8%, which were much higher than those reported in previous research. For LEAD, our RF model showed superior performance, aligning with previous findings that highlighted the RF model’s enhanced predictive capabilities over the LR model [4]. Of note, the improved performance in previous studies was attributed to the inclusion of the ankle-brachial pressure index (ABI), a common indicator for diagnosing LEAD. Our study, however, relied solely on clinical data and routine laboratory tests to construct ML models. The remarkable performance of our models can be attributed to our methods of feature selection and hyperparameter optimization. We combined three different methods—MI, RF-RFE, and the Boruta algorithm—to identify the most significant features. This approach significantly reduces overfitting and enhances the robustness of our models. Besides, PSO was applied to optimize hyperparameters. Unlike traditional methods such as grid search, PSO does not rely on fixed parameter value range and step size, making it particularly effective for complex optimization challenges with large parameter spaces.

Table 3 Comparative analysis of the proposed work with previous studies for DPN and LEAD prediction models
Full size table

SHAP, a game theory-based method, was used in this study to identify key risk factors for DPN and LEAD. The analysis revealed that the primary risk factors common to both conditions were increased UACR, HbA1c, Scr, advanced age, and carotid stenosis. Notably, UACR was ranked as the most crucial predictor for DPN and the third most significant for LEAD. This finding was consistent with a large retrospective cohort study that identified UACR as a crucial predictor for DPN [27]. Additionally, a 30% or greater increase in UACR was reported to be a risk factor for the onset of DPN [28]. Previous studies also discovered that UACR served as a biomarker for the early detection of LEAD [29] and a risk factor for mortality in LEAD patients [30]. This underscored the critical importance of regular UACR monitoring to prevent DPN and LEAD, thereby potentially reducing the risk of DFUs.

As expected, HbA1c and older age were critical shared risk factors for both conditions, aligning with previous studies [27, 31,32,33]. Unlike FBG, which can fluctuate significantly due to various factors, HbA1c provides a more stable measure of blood glucose levels over the preceding three months. Chronic hyperglycemia in diabetes contributed to the development of DPN through mechanisms such as increased oxidative stress and inflammation [34]. These processes disrupt blood flow to peripheral nerves and impair nerve function. Chronic hyperglycemia can also cause damage to endothelial cells and thicken the intima-media layer in blood vessels, particularly in the lower extremities [35]. Additionally, as individuals age, the key components of the extracellular matrix, particularly elastic fibers, are subjected to degradation and fragmentation. Age-related increases in cross-linking between collagen fibers could further contribute to the development of arterial stiffness [36], which may diminish blood flow to nerves and affect their repair capabilities, potentially increasing the prevalence of DPN [37]. These findings underscored the importance of maintaining good glycemic control, especially in older patients.

Scr is a key marker for kidney function, with elevated levels often indicating renal damage. Impaired kidney function can affect the microcirculation in distant organs [38], potentially compromising blood flow to peripheral nerves and arteries, which increases the risk of DPN and LEAD. Moreover, carotid stenosis was also recognized as a significant risk factor for both conditions. While the direct link between carotid stenosis and DPN is less studied, recent findings suggested that carotid atherosclerosis, the primary cause of carotid stenosis, could independently predict small fiber nerve dysfunction in individuals with T2DM [39]. Furthermore, a cross-sectional study of 653 patients with LEAD found that 415 (63.5%) had carotid stenosis [40], implying that carotid stenosis may be a contributing risk factor for LEAD. Therefore, diabetic patients should also pay more attention on kidney function and neck vascular health to reduce the prevalence of DPN and LEAD.

Furthermore, unique risk factors were also identified. For DPN, decreased serum albumin was a critical predictor. Among patients with T2DM, a serum albumin level below 36.75g/L was independently associated with impaired peripheral nerve function, with a sensitivity of 65.6% and a specificity of 78.0% for detecting abnormal function in those with albuminuria [41]. Recent studies further supported the inverse relationship between serum albumin levels and the prevalence of DPN among T2DM patients [42, 43]. These findings suggest that serum albumin may play a protective role against the development of DPN, potentially due to its antioxidant, anti-inflammatory, and anti-atherosclerotic properties [42]. Another unique but often overlooked risk factor for DPN was lymphocyte count. Both serum albumin and lymphocyte count are indicators for nutritional status [44], highlighting the importance for patients with DPN to closely monitor and manage their nutrition.

For LEAD, elevated NLR was identified as a unique key risk factor, consistent with previous studies [45, 46], which discovered that NLR was positively related with the prevalence of LEAD. D-dimer was identified as another crucial predictor for LEAD. In a prospective cohort study, patients with LEAD had significant higher levels of D-dimer than those without LEAD [47]. In addition, the levels of D-dimer were observed to increase with the severity of LEAD [48]. Elevated D-dimer levels may reflect the extent of atherosclerosis, as they indicate ongoing fibrin formation and degradation [49].

Conclusion

Our study underscored the potential of ML models in predicting DPN and LEAD among diabetic patients. We found that XGBoost showed superior performance in identifying DPN, whereas RF model was more effective for diagnosing LEAD. SHAP analysis revealed the top five most critical risk factors common to both conditions, including elevated UACR, HbA1c, Scr, advanced age, and carotid stenosis. Additionally, unique predictors were identified for each condition: decreased serum albumin and lymphocyte count were associated with DPN, whereas increased NLR and D-dimer levels were linked to LEAD. These insights underscored the complexity of managing DPN and LEAD, emphasizing the need for personalized and comprehensive treatment strategies. Implementing these insights could enhance early detection and management of these diabetic complications, particularly beneficial in regions with limited medical resources. Prioritizing the management of shared risk factors, like glycemic control, renal function, and macrovascular health, may reduce the frequency of DPN and LEAD, thereby decreasing the risk of DFUs. Patients with DPN should also focus on maintaining good nutritional health. For future progress, research should be expanded to include a broader and more diverse population, and investigate the feasibility of developing a unified ML model capable of predicting both DPN and LEAD in individuals with diabetes.

Availability of data and materials

Some or all datasets generated during and analyzed during the current study are not publicly available but are available from the corresponding author on reasonable request.

Abbreviations

DPN:

Diabetic peripheral neuropathy

LEAD:

Lower extremity arterial disease

DFUs:

Diabetic foot ulcers

T2DM:

Type 2 diabetes mellitus

MI:

Mutual information

RF-RFE:

Random forest recursive feature elimination

ML:

Machine learning

LR:

Logistic regression

RF:

Random forest

XGBoost:

EXtreme Gradient Boosting

SVM:

Support vector machine

PSO:

Particle swarm optimization

SHAP:

SHapley Additive exPlanation

ROC:

Receiver operating characteristic

AUC:

Area under the receiver operating curve

UACR:

Urinary albumin-to-creatinine ratio

EMG:

Electromyography

DKA:

Diabetic ketoacidosis

HHS:

Hyperosmotic hyperglycemia syndrome

BMI:

Body mass index

SBP:

Systolic blood pressure

DBP:

Diastolic blood pressure

NLR:

Neutrophil-to-lymphocyte ratio

ALT:

Alanine aminotransferase

AST:

Aspartate aminotransferase

TBiL:

Total bilirubin

DBiL:

Direct bilirubin

SUN:

Serum urea nitrogen

SUA:

Serum uric acid

Scr:

Serum creatinine

FBG:

Fasting blood glucose

HbA1c:

Glycosylated hemoglobin

TC:

Total cholesterol

HDL:

High-density lipoprotein

LDL:

Low-density lipoprotein

TG:

Triglycerides

CRP:

C-reactive protein

25-OH VitD:

25-Hydroxy vitamin D3

NSE:

Neuron specific enolase

eGFR:

Estimated glomerular filtration rate

References

  1. Cho NH, Shaw JE, Karuranga S, Huang Y, da Rocha Fernandes JD, Ohlrogge AW, et al. IDF Diabetes Atlas: Global estimates of diabetes prevalence for 2017 and projections for 2045. Diabetes Res Clin Pract. 2018;138:271–81. https://doi.org/10.1016/j.diabres.2018.02.023.

    Article  CAS  PubMed  Google Scholar 

  2. Wang L, Peng W, Zhao Z, Zhang M, Shi Z, Song Z, et al. Prevalence and treatment of diabetes in China, 2013–2018. JAMA. 2021;326:2498–506. https://doi.org/10.1001/jama.2021.22208.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Feldman EL, Callaghan BC, Pop-Busui R, Zochodne DW, Wright DE, Bennett DL, et al. Diabetic neuropathy. Nat Rev Dis Primers. 2019;5:41. https://doi.org/10.1038/s41572-019-0092-1.

    Article  PubMed  Google Scholar 

  4. Gao JM, Ren ZH, Pan X, Chen YX, Zhu W, Li W, et al. Identifying peripheral arterial disease in the elderly patients using machine-learning algorithms. Aging Clin Exp Res. 2022;34:679–85. https://doi.org/10.1007/s40520-021-01985-x.

    Article  PubMed  Google Scholar 

  5. Firnhaber JM, Powell CS. Lower extremity peripheral artery disease: diagnosis and treatment. Am Fam Physician. 2019;99:362–9.

    PubMed  Google Scholar 

  6. Buso G, Aboyans V, Mazzolai L. Lower extremity artery disease in patients with type 2 diabetes. Eur J Prev Cardiol. 2019;26:114–24. https://doi.org/10.1177/2047487319880044.

    Article  PubMed  Google Scholar 

  7. Abbas ZG, Boulton AJM. Diabetic foot ulcer disease in African continent: ‘From clinical care to implementation’ - Review of diabetic foot in last 60 years - 1960 to 2020. Diabetes Res Clin Pract. 2022;183:109155. https://doi.org/10.1016/j.diabres.2021.109155.

    Article  PubMed  Google Scholar 

  8. Hicks CW, Selvin E. Epidemiology of peripheral neuropathy and lower extremity disease in diabetes. Curr Diab Rep. 2019;19:86. https://doi.org/10.1007/s11892-019-1212-8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Abián MF, Vanesa BB, Diego BG, Manuel GS, Maria VC, Raquel VS, et al. Frequency of lower extremity artery disease in type 2 diabetic patients using pulse oximetry and the ankle-brachial index. Int J Med Sci. 2021;18:2776–82. https://doi.org/10.7150/ijms.58907.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Yang K, Wang Y, Li YW, Chen YG, Xing N, Lin HB, et al. Progress in the treatment of diabetic peripheral neuropathy. Biomed Pharmacother. 2022;148:112717. https://doi.org/10.1016/j.biopha.2022.112717.

    Article  CAS  PubMed  Google Scholar 

  11. Lian X, Qi J, Yuan M, Li X, Wang M, Li G, et al. Study on risk factors of diabetic peripheral neuropathy and establishment of a prediction model by machine learning. BMC Med Inform Decis Mak. 2023;23:146. https://doi.org/10.1186/s12911-023-02232-1.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Allwright M, Karrasch JF, O’Brien JA, Guennewig B, Austin PJ. Machine learning analysis of the UK Biobank reveals prognostic and diagnostic immune biomarkers for polyneuropathy and neuropathic pain in diabetes. Diabetes Res Clin Pract. 2023;201: 110725. https://doi.org/10.1016/j.diabres.2023.110725.

    Article  CAS  PubMed  Google Scholar 

  13. Zhang X, Sun Y, Ma Z, Lu L, Li M, Ma X. Machine learning models for diabetic neuropathy diagnosis using microcirculatory parameters in type 2 diabetes patients. Int Angiol. 2023;42:191–200. https://doi.org/10.23736/s0392-9590.23.05008-3.

    Article  PubMed  Google Scholar 

  14. Kazemi M, Moghimbeigi A, Kiani J, Mahjub H, Faradmal J. Diabetic peripheral neuropathy class prediction by multicategory support vector machine model: a cross-sectional study. Epidemiol Health. 2016;38:e2016011. https://doi.org/10.4178/epih.e2016011.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Armstrong DG, Tan TW, Boulton AJM, Bus SA. Diabetic foot ulcers: a review. JAMA. 2023;330:62–75. https://doi.org/10.1001/jama.2023.10578.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Yang L, Rong GC, Wu QN. Diabetic foot ulcer: challenges and future. World J Diabetes. 2022;13:1014–34. https://doi.org/10.4239/wjd.v13.i12.1014.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Chinese Diabetes Society. Guideline for the prevention and treatment of type 2 diabetes mellitus in China (2020 edition). Int J Endocrinol Metab. 2021;41:482–548. https://doi.org/10.3760/cma.j.cn121383-20210825-08063.

  18. Levey AS, Stevens LA, Schmid CH, Zhang YL, Castro AF 3rd, Feldman HI, et al. A new equation to estimate glomerular filtration rate. Ann Intern Med. 2009;150:604–12. https://doi.org/10.7326/0003-4819-150-9-200905050-00006.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Tong T, Ledig C, Guerrero R, Schuh A, Koikkalainen J, Tolonen A, et al. Five-class differential diagnostics of neurodegenerative diseases using random undersampling boosting. Neuroimage Clin. 2017;15:613–24. https://doi.org/10.1016/j.nicl.2017.06.012.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Gonzalez-Lopez J, Ventura S, Cano A. Distributed Selection of Continuous Features in Multilabel Classification Using Mutual Information. IEEE Trans Neural Netw Learn Syst. 2020;31:2280–93. https://doi.org/10.1109/tnnls.2019.2944298.

    Article  PubMed  Google Scholar 

  21. Simic V, Ebadi Torkayesh A, Ijadi Maghsoodi A. Locating a disinfection facility for hazardous healthcare waste in the COVID-19 era: a novel approach based on Fermatean fuzzy ITARA-MARCOS and random forest recursive feature elimination algorithm. Ann Oper Res. 2022:1-46. https://doi.org/10.1007/s10479-022-04822-0.

  22. Zhou H, Xin Y, Li S. A diabetes prediction model based on Boruta feature selection and ensemble learning. BMC Bioinformatics. 2023;24:224. https://doi.org/10.1186/s12859-023-05300-5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Albadr MAA, Tiun S, Ayob M, Al-Dhief FT. Particle swarm optimization-based extreme learning machine for COVID-19 detection. Cognit Comput. 2022:1-16. https://doi.org/10.1007/s12559-022-10063-x.

  24. Nohara Y, Matsumoto K, Soejima H, Nakashima N. Explanation of machine learning models using shapley additive explanation and application for real data in hospital. Comput Methods Programs Biomed. 2022;214:106584. https://doi.org/10.1016/j.cmpb.2021.106584.

    Article  PubMed  Google Scholar 

  25. Metsker O, Magoev K, Yakovlev A, Yanishevskiy S, Kopanitsa G, Kovalchuk S, et al. Identification of risk factors for patients with diabetes: diabetic polyneuropathy case study. BMC Med Inform Decis Mak. 2020;20:201. https://doi.org/10.1186/s12911-020-01215-w.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Rashid M, Alkhodari M, Mukit A, Ahmed KIU, Mostafa R, Parveen S, et al. Machine learning for screening microvascular complications in type 2 diabetic patients using demographic, clinical, and laboratory profiles. J Clin Med. 2022;11. https://doi.org/10.3390/jcm11040903.

  27. Schallmoser S, Zueger T, Kraus M, Saar-Tsechansky M, Stettler C, Feuerriegel S. Machine learning for predicting micro- and macrovascular complications in individuals with prediabetes or diabetes: retrospective cohort study. J Med Internet Res. 2023;25:e42181. https://doi.org/10.2196/42181.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Zhong M, Yang YR, Zhang YZ, Yan SJ. Change in urine albumin-to-creatinine ratio and risk of diabetic peripheral neuropathy in type 2 diabetes: a retrospective cohort study. Diabetes Metab Syndr Obes. 2021;14:1763–72. https://doi.org/10.2147/dmso.S303096.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Lee YH, Kweon SS, Choi JS, Rhee JA, Nam HS, Jeong SK, et al. Determining the optimal cut-off value of the urinary albumin-to-creatinine ratio to detect atherosclerotic vascular diseases. Kidney Blood Press Res. 2012;36:290–300. https://doi.org/10.1159/000343418.

    Article  CAS  PubMed  Google Scholar 

  30. Amrock SM, Weitzman M. Multiple biomarkers for mortality prediction in peripheral arterial disease. Vasc Med. 2016;21:105–12. https://doi.org/10.1177/1358863x15621797.

    Article  CAS  PubMed  Google Scholar 

  31. Dagliati A, Marini S, Sacchi L, Cogni G, Teliti M, Tibollo V, et al. Machine learning methods to predict diabetes complications. J Diabetes Sci Technol. 2018;12:295–302. https://doi.org/10.1177/1932296817706375.

    Article  PubMed  Google Scholar 

  32. Liu J, Yuan X, Liu J, Yuan G, Sun Y, Zhang D, et al. Risk factors for diabetic peripheral neuropathy, peripheral artery disease, and foot deformity among the population with diabetes in Beijing, China: a multicenter, cross-sectional study. Front Endocrinol (Lausanne). 2022;13:824215. https://doi.org/10.3389/fendo.2022.824215.

    Article  PubMed  Google Scholar 

  33. Wang W, Ji Q, Ran X, Li C, Kuang H, Yu X, et al. Prevalence and risk factors of diabetic peripheral neuropathy: a population-based cross-sectional study in China. Diabetes Metab Res Rev. 2023;39:e3702. https://doi.org/10.1002/dmrr.3702.

    Article  CAS  PubMed  Google Scholar 

  34. Zhu J, Hu Z, Luo Y, Liu Y, Luo W, Du X, et al. Diabetic peripheral neuropathy: pathogenetic mechanisms and treatment. Front Endocrinol (Lausanne). 2023;14:1265372. https://doi.org/10.3389/fendo.2023.1265372.

    Article  PubMed  Google Scholar 

  35. Dokun AO, Chen L, Lanjewar SS, Lye RJ, Annex BH. Glycaemic control improves perfusion recovery and VEGFR2 protein expression in diabetic mice following experimental PAD. Cardiovasc Res. 2014;101:364–72. https://doi.org/10.1093/cvr/cvt342.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Arribas SM, Hinek A, González MC. Elastic fibres and vascular structure in hypertension. Pharmacol Ther. 2006;111:771–91. https://doi.org/10.1016/j.pharmthera.2005.12.003.

    Article  CAS  PubMed  Google Scholar 

  37. Kilo S, Berghoff M, Hilz M, Freeman R. Neural and endothelial control of the microcirculation in diabetic peripheral neuropathy. Neurology. 2000;54:1246–52. https://doi.org/10.1212/wnl.54.6.1246.

    Article  CAS  PubMed  Google Scholar 

  38. Krishnan S, Suarez-Martinez AD, Bagher P, Gonzalez A, Liu R, Murfee WL, et al. Microvascular dysfunction and kidney disease: challenges and opportunities? Microcirculation. 2021;28:e12661. https://doi.org/10.1111/micc.12661.

    Article  PubMed  Google Scholar 

  39. Guo S, Jing Y, Li C, Zhu D, Wang W. Carotid atherosclerosis: an independent risk factor for small fiber nerve dysfunction in patients with type 2 diabetes mellitus. J Diabetes Investig. 2023;14:289–96. https://doi.org/10.1111/jdi.13936.

    Article  CAS  PubMed  Google Scholar 

  40. Li Z, Yang H, Zhang W, Wang J, Zhao Y, Cheng J. Prevalence of asymptomatic carotid artery stenosis in Chinese patients with lower extremity peripheral arterial disease: a cross-sectional study on 653 patients. BMJ Open. 2021;11:e042926. https://doi.org/10.1136/bmjopen-2020-042926.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Li L, Liu B, Lu J, Jiang L, Zhang Y, Shen Y, et al. Serum albumin is associated with peripheral nerve function in patients with type 2 diabetes. Endocrine. 2015;50:397–404. https://doi.org/10.1007/s12020-015-0588-8.

    Article  CAS  PubMed  Google Scholar 

  42. Yan P, Tang Q, Wu Y, Wan Q, Zhang Z, Xu Y, et al. Serum albumin was negatively associated with diabetic peripheral neuropathy in Chinese population: a cross-sectional study. Diabetol Metab Syndr. 2021;13:100. https://doi.org/10.1186/s13098-021-00718-4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Hu Y, Wang J, Zeng S, Chen M, Zou G, Li Y, et al. Association between serum albumin levels and diabetic peripheral neuropathy among patients with type 2 diabetes: effect modification of body mass index. Diabetes Metab Syndr Obes. 2022;15:527–34. https://doi.org/10.2147/dmso.S347349.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Borda MG, Salazar-Londoño S, Lafuente-Sanchis P, Patricio Baldera J, Venegas LC, Tarazona-Santabalbina FJ, et al. Neutrophil-to-lymphocyte ratio and lymphocyte count as an alternative to body mass index for screening malnutrition in older adults living in the community. Eur J Nutr. 2024 https://doi.org/10.1007/s00394-024-03392-0.

  45. Liu N, Sheng J, Pan T, Wang Y. Neutrophil to lymphocyte ratio and platelet to lymphocyte ratio are associated with lower extremity vascular lesions in chinese patients with type 2 diabetes. Clin Lab. 2019;65. https://doi.org/10.7754/Clin.Lab.2018.180804.

  46. Santoro L, Ferraro PM, Nesci A, D’Alessandro A, Macerola N, Forni F, et al. Neutrophil-to-lymphocyte ratio but not monocyte-to-HDL cholesterol ratio nor platelet-to-lymphocyte ratio correlates with early stages of lower extremity arterial disease: an ultrasonographic study. Eur Rev Med Pharmacol Sci. 2021;25:3453–9. https://doi.org/10.26355/eurrev_202105_25826.

    Article  CAS  PubMed  Google Scholar 

  47. Unlü Y, Karapolat S, Karaca Y, Kiziltunç A. Comparison of levels of inflammatory markers and hemostatic factors in the patients with and without peripheral arterial disease. Thromb Res. 2006;117:357–64. https://doi.org/10.1016/j.thromres.2005.03.019.

    Article  CAS  PubMed  Google Scholar 

  48. Lee AJ, Fowkes FG, Lowe GD, Rumley A. Fibrin D-dimer, haemostatic factors and peripheral arterial disease. Thromb Haemost. 1995;74:828–32.

    Article  CAS  PubMed  Google Scholar 

  49. Moresco RN, Silla L. D-dimer and inflammatory markers in the peripheral arterial disease. Thromb Res. 2007;119:797–8. https://doi.org/10.1016/j.thromres.2006.08.002.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work was supported by grants from the National Natural Science Foundation of China (No. 82072024).

Author information

Authors and Affiliations

  1. Department of Endocrinology and Metabolism, Shanghai Tenth People’s Hospital, Tongji University School of Medicine, Shanghai, China

    Ya Wu, Lijie Zhu & Xiaoyun Xie

  2. School of Medicine, Anhui University of Science and Technology, Huainan, China

    Danmeng Dong

  3. School of Advanced Technology, Xi’an Jiaotong-Liverpool University, Suzhou, China

    Zihong Luo

  4. Department of Geriatrics, Tongji Hospital, Tongji University School of Medicine, Shanghai, China

    Yang Liu

Authors
  1. Ya Wu
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Danmeng Dong
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Lijie Zhu
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Zihong Luo
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Yang Liu
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Xiaoyun Xie
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

Ya Wu: Data acquisition and analysis, Methodology, Writing – original draft. Danmeng Dong and Lijie Zhu: Data acquisition. Zihong Luo: Methodology. Yang Liu: Data curation. Xiaoyun Xie: Conceptualization, Writing – review & editing.

Corresponding author

Correspondence to Xiaoyun Xie.

Ethics declarations

Ethics approval and consent to participate

This study was performed in compliance with the Code of Ethics of the World Medical Association (Declaration of Helsinki), and received approval from the Institutional Ethics Committee of Tongji hospital, Tongji University School of Medicine (K-2023-022). Due to the retrospective nature of the study, the Ethics Committee of Tongji hospital has waived informed consent for the study. All patient data were anonymized, with no identifiable personal information. The study was registered at www.chictr.org.cn, identifier ChiCTR2400087019.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

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

Supplementary Information

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

Wu, Y., Dong, D., Zhu, L. et al. Interpretable machine learning models for detecting peripheral neuropathy and lower extremity arterial disease in diabetics: an analysis of critical shared and unique risk factors. BMC Med Inform Decis Mak 24, 200 (2024). https://doi.org/10.1186/s12911-024-02595-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12911-024-02595-z

Keywords

BMC Medical Informatics and Decision Making

ISSN: 1472-6947

Contact us