Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Endocrinol., 12 February 2026

Sec. Clinical Diabetes

Volume 17 - 2026 | https://doi.org/10.3389/fendo.2026.1758234

This article is part of the Research TopicHarnessing Machine Learning for Enhanced Biomedical Diagnosis and Early Disease Detection: Bridging Data Science and HealthcareView all 14 articles

Development and validation of a nomogram model integrating noninvasive detection of radial pulse wave for predicting diabetic foot risk in type 2 diabetes mellitus

Yuting Huo&#x;Yuting Huo1†Wen Li&#x;Wen Li1†Ai XuAi Xu1Siman WangSiman Wang1Yuqing WangYuqing Wang2Jin XuJin Xu1Rui GuoRui Guo1Cheng Zhao*Cheng Zhao2*Haixia Yan*Haixia Yan1*
  • 1School of Traditional Chinese Medicine, Shanghai University of Traditional Chinese Medicine (TCM), Shanghai, China
  • 2Department of Vascular Disease, Shanghai Traditional Chinese Medicine (TCM)- Integrated Hospital, Shanghai University of Traditional Chinese Medicine (TCM), Shanghai, China

Background: This study aimed to develop and evaluate a nomogram model based on Traditional Chinese Medicine (TCM) pulse mapping parameters to predict the risk of diabetic foot (DF) in patients with type 2 diabetes mellitus (T2DM).

Methods: A retrospective analysis was performed on clinical data and pulse mapping parameters from 485 T2DM patients hospitalized in the Departments of Vascular Diseases and Endocrinology at Shanghai TCM-Integrated Hospital between July 26, 2023, and April 9, 2025. Patients were divided into a DF group (n = 242) and a T2DM-only group (n = 243) according to the presence of DF. Pulse mapping parameters and clinical indicators were compared between the two groups. Least Absolute Shrinkage and Selection Operator (LASSO) regression and multivariate logistic regression analyses were conducted to identify risk factors associated with DF development in T2DM. The risk prediction nomogram was constructed using R Studio. Model performance was evaluated using the Receiver Operating Characteristic (ROC) curve and its Area Under the Curve (AUC). Accuracy, Precision, Recall, and F1-score were calculated. Calibration curves were plotted to assess the agreement between predicted probabilities and actual outcomes, and decision curve analysis (DCA) was used to evaluate the model’s clinical utility.

Results: Gender (OR = 0.501, p = 0.002), BMI (OR = 0.919, p = 0.003), hypertension (OR = 2.482, p < 0.001), systolic blood pressure (OR = 0.972, p < 0.001), blood urea nitrogen (BUN) (OR = 1.099, p = 0.011), and the pulse mapping parameter H5/H1 (OR = 0.002, p < 0.001) were identified as independent risk factors for DF in T2DM. The nomogram model based on these factors achieved an AUC of 0.786 (95% CI: 0.694–0.877) in the test set, with an Accuracy of 0.7113, Precision of 0.6481, Recall of 0.7955, and F1-score of 0.7143. The Hosmer–Lemeshow test demonstrated good calibration (χ² = 6.1672, p = 0.6285), showing close agreement between predicted and observed outcomes. The DCA confirmed favorable clinical applicability.

Conclusion: This nomogram model provides an individualized and visualized tool for predicting DF risk in T2DM patients, offering TCM-based reference indicators for early DF risk assessment.

1 Introduction

Diabetes mellitus (DM) is one of the most prevalent chronic metabolic diseases worldwide and represents a major public health challenge (1). Type 2 diabetes mellitus (T2DM) accounts for over 90% of all DM cases (2). The number of people with DM has steadily increased over the past three decades, especially in China, which has one of the largest diabetic populations. Diabetic foot (DF) is one of the most serious complications of DM, characterized by foot ulcers, infections, and deep tissue destruction that can lead to amputation, severely impairing patients’ quality of life. Therefore, early diagnosis, prevention, and treatment of DF are vital for reducing the medical burden and mortality rate.

The main pathological changes in DF involve microcirculatory disorders and neurovascular damage caused by metabolic disturbances (3). Although the International Working Group on the DF (IWGDF) guidelines emphasize early risk stratification and management, most existing prediction systems rely on static indicators (e.g., ankle–brachial index, vibration perception threshold). These structure-based assessments fail to capture the critical stage when microcirculatory dysfunction progresses to tissue necrosis—the optimal window for clinical intervention.

Pulse taking is an important skill for diagnosing diseases in TCM. Objective pulse diagnosis research results enable its clinical application. Studies conducted by Yiming Hao et al. reported the noninvasive and convenient pulse-taking diagnosis has the potential to become a low-cost and accurate method to monitor the development of diabetes (4). Jiekee Lim et al. found that the new prediction model that uses pulse wave parameters and TCM clinical indices offers a non-invasive and cost-effective way to diagnose PCOS (5).

At present, the relationship between pulse wave and DF has not been reported. Current efforts to objectify and standardize TCM pulse diagnosis have achieved progress (6, 7), providing a basis for DF risk prediction model development. The nomogram model is widely applied in medical and biological research. It serves as a visual tool for multivariate predictive modeling, enabling clinicians, researchers, and patients to efficiently assess risks, prognoses, or diagnostic outcomes (8, 9). Hence, this study explores DF risk factors in T2DM patients by integrating noninvasive detection of radial pulse wave to construct a predictive nomogram model. The goal is to offer a novel, simple, non-invasive, and cost-effective early warning tool for DF, promote the application of TCM pulse diagnosis in chronic disease management, reduce DF-related disability, and improve patient outcomes.

2 Data and methods

2.1 Study participants

In this exploratory phase of the study, to mitigate the potential bias introduced by class imbalance in the training data and to enhance the model’s ability to identify DF features while improving training stability and efficiency, we applied dataset balancing techniques during model development. This approach facilitates a more stable preliminary evaluation of model performance in the development stage. To avoid the risk of introducing artifactual noise or undermining the biological fidelity of the data, we refrained from using synthetic oversampling methods (e.g., SMOTE) that generate data points not grounded in actual clinical observations.

Final model validation was conducted using an independent test set that retained the original observed prevalence of the condition. This strategy ensures that performance metrics—such as sensitivity, specificity, and positive predictive value—reflect the model’s generalizability to real-world populations, thereby providing clinically relevant evidence with greater external validity for decision-making.

Due to the retrospective design of this study, DF diagnoses documented in the electronic medical records often lacked standardized staging information, such as Wagner classification. Consequently, stratification of all DF patients according to disease severity was not feasible. The distribution of patients for whom staging information was available is presented in Supplementary Material 1.

This retrospective study included 485 patients with T2DM who were admitted to the Departments of Vascular Medicine and Endocrinology at Shanghai Integrated Traditional Chinese and Western Medicine Hospital between July 2023, and April 2025. Participants were included in the study if they were aged 25 years or older, had a confirmed diagnosis of T2DM, and had complete clinical data. Based on the presence of DF, patients were divided into a DF group (n = 242) and a T2DM group (n = 243).

Diagnostic criteria for T2DM and DF were aligned with the China Guidelines for the Prevention and Treatment of Diabetes (2024 Edition) (10). The detailed criteria with patient inclusion and exclusion criteria are described in Supplementary Materials 2.The diagnosis for T2DM was validated through clinical evaluation and laboratory findings, including HbA1c and fasting blood glucose levels, with a documented disease duration of at least one year. The study flowchart is shown in Figure 1.

Figure 1
Flowchart detailing a study of 618 patients with Type 2 Diabetes Mellitus from July 2023 to April 2025. Data was collected via pulse wave analysis and clinical indicators, with missing values addressed through simple imputation. After excluding patients due to tumors, cognitive impairments, and incomplete data, 485 participants remained: 242 in the diabetic foot group and 243 in the control group. Using LASSO and logistic regression, six predictors were identified. A predictive nomogram model was constructed and validated using bootstrap resampling with one thousand iterations, assessing ROC, calibration curve, and DCA.

Figure 1. Flowchart of the study.

The study was reviewed and approved by the Ethics Committee of Shanghai Integrated Traditional Chinese and Western Medicine Hospital, affiliated with Shanghai University of Traditional Chinese Medicine (Approval No. 2022-004-1). Informed consent is obtained from all participants before data collection.

2.2 Measurement indicators

2.2.1 Pulse waveform parameter acquisition

Pulse waveform parameters were collected using the digital pulse analysis system for TCM (11, 12). The standardized procedure was as follows: participants were placed in a sitting or supine position with the forearm naturally extended and maintained at heart level. The wrist was kept straight and relaxed, with the palm facing upward. A pulse pillow was placed under the wrist joint, and a pressure sensor was positioned over the guan region of the cunkou (13) (The cunkou is located at the radial side of the wrist, in the depression between the styloid process of the radius and the radial artery. The guan region is a fixed landmark within the cunkou pulse). The optimal pulse-taking pressure was determined dynamically through a stepwise pressure increase method. Once the optimal pressure was identified, participants were instructed to maintain quiet and steady breathing. Pulse signals were continuously recorded for 60 seconds under this optimal pressure. All pulse wave measurements were conducted by a single operator who had received standardized training, and all assessments were performed between 10:00 AM and 4:00 PM. Patients were examined in the supine position following at least 10 minutes of rest to minimize the influence of diurnal variations, body posture, and recent physical activity.

The system automatically extracted and analyzed time-domain pulse waveform parameters, including T1, T4, T5, T, H1, H3, H4, H5, W1, W2, As, and Ad. The following parameter ratios were also calculated: H3/H1, H4/H1, H5/H1, T1/T, W1/T, W2/T, and As/Ad (14, 15), as shown in Figure 2.

Figure 2
Graph depicting height (h) in millimeters versus time (t) in seconds. Peaks and troughs are labeled as H1, H3, H4, and H5 indicating height variations. Time intervals are marked as T1, T4, T5, and total duration T. Widths W1 and W2 are indicated near the first peak.

Figure 2. Time-domain parameters of the pulse wave.

2.2.2 Clinical data

Clinical data were obtained from the hospital’s electronic medical record system. Collected variables included: (1)demographic and clinical characteristics: gender, age, body mass index(BMI), history of smoking and drinking; (2) blood pressure value: systolic blood pressure (SBP), and diastolic blood pressure (DBP), hypertension (HTN); (3)serological markers (preoperative laboratory testing): triglycerides (TG), total cholesterol (TC), low-density lipoprotein (LDL), high-density lipoprotein (HDL), serum creatinine (Cr), blood urea nitrogen (BUN), uric acid (UA), FBG, homocysteine (Hcy), and D-dimer levels.

2.3 Statistical methods

To ensure sufficient statistical power and generalization ability, sample size was estimated based on R² shrinkage and the C-statistic for a prediction model with 6 predictors(Variables selected according to the results of the multiple regression analysis: Gender, BMI, HTN, SBP, BUN, H5/H1). The target shrinkage factor was ≥0.90 (16), with an acceptable difference in R² of 0.05. Based on Chinese guidelines and epidemiological data, the prevalence of DF in patients with T2DM is approximately 20%. The target C-statistic and R² were set at 0.63 and 0.11.

All statistical analyses and data visualizations were performed using R (version 4.4.2) and JD_DCPM (V6.11, Jingding Medical Technology Co., Ltd.).

Calculation using the “pmsampsize” package in R indicated a minimum required sample size of 461 participants, including 93 events. This corresponds to 15.37 Events Per Predictor (EPV), which exceeds the recommended threshold of EPV ≥10 (17). The study ultimately enrolled 485 patients with T2DM, of whom 242 were in the DF group, ensuring an adequate number of outcome events. Therefore, the sample size is sufficient for both model development and internal validation.

Data entry was conducted using Microsoft Excel 2021 with double entry and verification. Parameters with more than 10% missing data were excluded. For parameters with less than 10% missing data, simple imputation was performed using R Studio. For continuous variables exhibiting a non-normal distribution, missing values will be imputed using the median of the observed (non-missing) values. For continuous variables with a normal distribution, the mean of the non-missing values will be used for imputation (the number and percentage of missing values for each predictor variable are presented in Supplementary Table 3). This imputation procedure is implemented using the ‘transcan’ function from the Hmisc package in RStudio.

Statistical analyses were performed using IBM SPSS 27.0. Normally distributed continuous variables are presented as mean ± standard deviation (x¯ ± S) and were compared using independent samples t-tests. Non-normally distributed data are reported as median (first quartile, third quartile) [M (Q1, Q3)] and were analyzed using the Mann–Whitney U test. Categorical variables are expressed as frequencies and percentages were compared using the χ² test.

A significance threshold of α = 0.05 was set, and p < 0.05 was considered statistically significant. Univariate logistic regression was first performed for preliminary variable screening. Variables showing statistical significance were further analyzed using LASSO regression combined with multivariate logistic regression (backward stepwise method) to identify independent risk factors for DF in T2DM patients (p < 0.05 considered significant).

2.4 Model construction

The R software with the rms package was used to construct the nomogram prediction model based on factors independently associated with DF risk in T2DM patients. The dependent variable was the presence or absence of DF. Data were randomly divided into a training set (70%) for model construction and a test set (30%) for model evaluation. Internal validation of the nomogram was conducted using the Bootstrap method with 1000 repetitions.

2.5 Model evaluation

Predictive performance was assessed using Accuracy, Precision, Recall, F1-score, and the Area Under the ROC Curve (AUC). Model calibration was evaluated using the Hosmer–Lemeshow test and calibration curves to assess the consistency between predicted probabilities and observed outcomes. Clinical applicability was analyzed using Decision Curve Analysis (DCA).

3 Results

3.1 Comparison of clinical data between groups

Comparisons between the T2DM and DF groups showed no significant differences in TC, LDL, UA, or D-dimer levels (p > 0.05). The DF group had significantly higher proportions of males, smokers, drinkers, and patients with HTN (coded as 0 = No, 1 = Yes in the table) (p < 0.05). Age, TG, Cr, BUN, and Hcy levels were significantly higher in the DF group (p < 0.05). BMI, DBP, SBP, and HDL were significantly lower in the DF group (p < 0.05) (Table 1).

Table 1
www.frontiersin.org

Table 1. Comparison of clinical data between T2DM and DF groups.

3.2 Comparison of pulse waveform parameters between T2DM groups with and without concurrent DF

When comparing the T2DM and DF groups, no significant differences were observed in H3/H1, T1/T4, T5/T4, W1/T, and W2/T ratios (p > 0.05). However, the As/Ad and T1/T ratios were significantly higher in the DF group (p < 0.05), while the H4/H1 and H5/H1 ratios were significantly lower (p < 0.05), as shown in Table 2.

Table 2
www.frontiersin.org

Table 2. Comparison of pulse mapping parameters between T2DM and DF groups.

3.3 Independent risk factors for DF in T2DM

Univariate logistic regression analysis was conducted on the 17 parameters that showed significant differences between groups. The results indicated that gender, smoking, alcohol consumption, HTN, age, BMI, SBP, DBP, TG, HDL, Cr, BUN, Hcy, and the pulse parameters As/Ad, H4/H1, H5/H1, and T1/T might be influencing factors for DF risk prediction in T2DM (p < 0.05) (Table 3).

Table 3
www.frontiersin.org

Table 3. Results of univariate logistic regression analysis.

Variables with statistical significance in the univariate analysis were included in the LASSO regression model for further screening. The analysis results are shown in Figures 3 and 4. The optimal lambda (λ) value was determined through 10-fold cross-validation. The two vertical dashed lines represent lambda.min (11 features) and lambda.1se (6 features). To achieve a parsimonious model with stable performance, lambda.1se was selected, resulting in six retained variables: gender, BMI, HTN, SBP, BUN, and H5/H1.

Figure 3
Plot showing binomial deviance versus log-transformed lambda values. Red dots represent data points with error bars. Two dashed vertical lines at approximately -5.3 and -2.9 mark specific values, with deviance increasing beyond these points.

Figure 3. Determination of the optimal value by 10-fold cross-validation in the LASSO regression.

Figure 4
Plot showing coefficients in a LASSO regression against the logarithm of Lambda. Each colored line represents different coefficients, converging near zero as Lambda increases. Dashed vertical lines likely indicate specific values of interest for Lambda.

Figure 4. LASSO coefficient path plot.

The variables selected by LASSO regression were incorporated into multivariate logistic regression analysis (backward stepwise method). The results identified gender, BMI, HTN, SBP, BUN, and the pulse waveform parameter H5/H1 as independent factors associated with the risk of T2DM complicated by DF (p < 0.05) (Table 4).

Table 4
www.frontiersin.org

Table 4. Results of multivariate logistic regression analysis for DF in T2DM.

3.4 Establishment and evaluation of the nomogram model

A nomogram model for predicting DF risk in T2DM was developed based on the optimal variables identified in the multivariate logistic regression analysis (Figure 5). Points were assigned for gender, BMI, HTN, SBP, BUN, and the pulse parameter H5/H1. The total score was calculated by summing the points, and the corresponding DF risk probability for each T2DM patient was determined based on the total score. The total score ranged from 0 to 200 points: <60 points indicated low risk, 60–120 points indicated medium risk, and >120 points indicated high risk.

Figure 5
A medical scoring chart with multiple scales indicating points based on different variables. Variables include Gender, BMI (Body Mass Index), HTN (Hypertension), SBP (Systolic Blood Pressure), BUN (Blood Urea Nitrogen), and H5H1. Each variable has a corresponding point scale. Total points and diagnostic possibilities are shown on separate axes.

Figure 5. Nomogram model for predicting DF risk in patients with T2DM.

The model underwent internal validation using the bootstrap method. Model discrimination was assessed with the ROC curve, calibration was evaluated using the calibration curve, and clinical applicability was examined through DCA. The AUC values for predicting the risk of T2DM complicated with DF were 0.750 [95% CI (0.702, 0.799)] in the training set and 0.786 [95% CI (0.694, 0.877)] in the validation set, demonstrating good discriminative ability of the model. The higher AUC in the validation set compared with the training set indicated that the model was not overfitted and exhibited stable predictive performance for new samples, showing good generalizability. The recall rate in the validation set was 0.7955, indicating that the model identified 79.55% of actual DF cases, effectively reducing the risk of missed diagnoses. These results suggest a satisfactory diagnostic capability, supporting the model’s use for early DF screening (Figure 6). The accuracy, precision, recall, and F1-score were 0.701 and 0.7113, 0.7278 and 0.6481, 0.6616 and 0.7955, and 0.6931 and 0.7143 for the training and validation sets, respectively (Table 5).

Figure 6
ROC curve graph comparing test and train data. The x-axis shows 1 minus specificity, and the y-axis shows sensitivity. The red line represents test data with an AUC of 0.786 (confidence interval: 0.694-0.877). The blue line represents train data with an AUC of 0.750 (confidence interval: 0.702-0.799). A diagonal dashed line indicates random performance.

Figure 6. ROC curves of the nomogram model.

Table 5
www.frontiersin.org

Table 5. Evaluation Metrics of the nomogram model.

The Hosmer–Lemeshow test result was χ² = 6.1672, p = 0.6285. The calibration curves showed close alignment with the ideal curve, indicating good agreement between predicted and observed values and confirming the model’s predictive accuracy (Figure 7). The DCA results showed that using this model for DF risk assessment in T2DM patients provided net clinical benefit within a threshold probability range of approximately 0% to 90%, suggesting favorable clinical applicability (Figure 8).

Figure 7
Two calibration plots compare actual probability versus predicted probability. Plot (A) shows near-perfect calibration with lines close to the ideal diagonal. Plot (B) deviates more, illustrating discrepancies. Both display statistical metrics like Dxy and Brier scores for model evaluation.

Figure 7. Calibration curves of the nomogram model. (A) Train set. (B) Test set.

Figure 8
Two decision curve graphs labeled A and B. Both graphs plot net benefit against threshold probability, ranging from 0 to 0.8, with cost-benefit ratios displayed on the x-axis. Graph A shows a decreasing red line labeled “Nomo” and a gray line labeled “All,” while graph B shows a similar trend but ends with an upward spike in the red line. Both graphs also include a baseline labeled “None.

Figure 8. DCA of the nomogram model. (A) Train set. (B) Test set.

4 Discussion

Data from the Global Burden of Disease 2021 database show that the disease burden of T2DM in China has continued to increase from 1990 to 2021 and is expected to keep rising in the next five years, presenting major challenges for prevention and control (18). DF is one of the complications that severely threaten the health and survival of DM patients. Early identification and timely intervention are essential for improving the quality of life of DF patients.

Clinical prediction models use mathematical modeling to assign weights to risk factors associated with specific clinical outcomes and visualize the results through tools such as nomograms (19). These models provide clinicians with evidence-based decision support, facilitating the development of individualized treatment strategies through precise risk assessment (20). This study incorporated TCM pulse mapping parameters into a DF risk prediction model for T2DM, aiming to provide a non-invasive, economical, and convenient method for DF risk prediction.

The results of this study indicate that the pulse parameter H5/H1 (OR = 0.002, p < 0.001) is a protective factor against DF risk in T2DM, showing a significant negative effect on risk prediction. This parameter offers a complementary perspective to established DF risk assessment indicators such as the Ankle-Brachial Index (ABI) and Toe-Brachial Index (TBI). While ABI and TBI primarily evaluate structural arterial stenosis or occlusion through direct limb pressure measurements (21), representing static anatomical anatomical assessments, the pulse wave parameter H5/H1 is extracted from the time-domain analysis of the radial artery waveform, reflecting dynamic functional information. Existing research has demonstrated that the radial artery pulse wave plays an important role in predicting adverse cardiovascular and microvascular events in diabetic patients (22). According to hemodynamic principles, H5/H1 is a quantitative indicator of the relative height of the dicrotic wave in the pulse waveform, reflecting large artery compliance and aortic function. H5/H1 decreases when arterial elasticity is reduced or aortic valve insufficiency occurs (23).Therefore, it offers dynamic functional insights into central arterial elasticity and hemodynamics. This represents a shift in assessment perspective—from “peripheral structure” to “central function”—which distinguishes the present study from a TCM viewpoint (24, 25). Existing prediction systems rely predominantly on static indicators, which often fail to capture the critical transition from microcirculatory dysfunction to tissue necrosis (26). In contrast, pulse wave parameters may provide enhanced sensitivity for detecting early functional deterioration within the vascular system. Previous studies have shown that lower extremity atherosclerosis is a common macrovascular complication of T2DM and plays a key role in the pathogenesis of DF (27). These findings are consistent with the pathophysiological mechanism of DM (28), Future studies involving the simultaneous measurement of pulse wave parameters and ABI/TBI will help elucidate the relationship between central arterial function and peripheral arterial structure, thereby establishing a more integrated risk assessment pathway from “functional warning” to “structural diagnosis”.

The results also show that patients with HTN (OR = 2.482) have a 2.482-fold higher risk of developing DF than those without HTN, consistent with previous studies (29). The coexistence of HTN and T2DM is closely associated with DF occurrence, aligning with findings from a multicenter study in Bangladesh (30). However, SBP was negatively correlated with DF risk, which differs from some existing studies. This may be explained by several factors: (1) treatment-related bias—patients in this study likely received antihypertensive therapy, and varying degrees of blood pressure control may have influenced the results; (2) the J-curve phenomenon—when SBP falls below 120 mmHg, toe microcirculation perfusion pressure may decrease, exceeding the lower limit of autoregulation and resulting in lower limb hypoperfusion, thereby increasing DF risk. The optimal blood pressure management strategy for T2DM remains controversial and inconclusive (31). Although elevated blood pressure is a major risk factor for cardiovascular disease in T2DM and blood pressure control can reduce cardiovascular risk, excessive lowering may also lead to adverse outcomes (32). Therefore, individualized blood pressure management is essential for patients with high DF risk, and clinicians should remain vigilant regarding the J-curve risk associated with overly aggressive blood pressure reduction or low-pressure states.

The findings indicate that a higher BMI is associated with a reduced risk of DF, which contrasts with certain prior studies reporting an increased risk of DF with elevated BMI (30). This observation partially supports the so-called “obesity paradox in diabetes” phenomenon (33, 34). However, this association should be interpreted with caution, as BMI does not differentiate between fat mass and lean muscle mass, thereby limiting its objectivity and comprehensiveness as an indicator of obesity. In individuals with long-standing diabetes, sarcopenic obesity—defined by the coexistence of excess adiposity and muscle atrophy—is frequently observed and may mask the true metabolic risk that BMI purports to reflect (35). Indeed, evidence suggests that predicted fat mass demonstrates greater discriminative accuracy than conventional metrics such as BMI and can independently predict the risk of T2MD (36). Furthermore, a linear relationship between predicted body fat percentage and T2DM incidence has been documented (37), indicating that fat accumulation may represent a distinct risk factor independent of BMI. The apparent protective effect of higher BMI observed in our cohort may be attributable to confounding factors or compensatory mechanisms, such as improved nutritional status, enhanced insulin sensitivity, or more intensive medical management among individuals with higher BMI (38).

Furthermore, the heterogeneous effects of weight loss interventions underscore the complexity of the relationship between obesity and clinical outcomes. For instance, bariatric surgery demonstrates more benefits in individuals with early-stage T2DM patients, whereas its efficacy is comparatively limited in patients with advanced complications such as DF (39, 40). Therefore, the associations between BMI and other obesity-related indicators warrant further investigation (41), with particular attention on the potential role of sarcopenic obesity in the pathogenesis of DF. Future multi-center, prospective cohort studies should integrate assessments of muscle mass, inflammatory biomarkers, and metabolomic profiles to comprehensively elucidate the intricate interplay between distinct obesity phenotypes and diabetic complications.

This study also found that male T2DM patients were more likely to develop DF than female patients, which agrees with previous research. This difference may be attributed to the protective effects of estrogen in women, which supports cardiovascular health through multiple mechanisms such as regulating vascular function, protecting the vascular endothelium, and preventing atherosclerosis (42). Combined with gender differences in lifestyle and self-care behaviors, understanding the protective mechanisms of estrogen may provide new insights for developing sex-specific DF prevention strategies.

The model constructed in this study includes six predictive variables (gender, BMI, HTN, SBP, BUN, and pulse parameter H5/H1). The total score is obtained by summing the points for each variable, and the predicted probability of DF occurrence is calculated based on the total score. This study quantified the contribution of each variable to DF risk, with H5/H1 serving as a key predictive variable. Clinicians can use this model to rapidly assess individual DF risk, addressing the issue of equal weighting of indicators in traditional scoring systems and providing a visual and individualized basis for clinical decision-making.

The model demonstrates good diagnostic performance and clinical applicability. It can support early identification of high-risk DF populations in DM, which is important for early intervention and for reducing DF-related disability and mortality rates. The recall of 0.7955 on the test set indicates that approximately 20.45% of DF cases may remain undetected. Therefore, this model should not be used as a substitute for comprehensive clinical evaluation, but rather as an auxiliary screening tool. In practice, the risk threshold could be adjusted to enhance sensitivity and minimize missed diagnoses, albeit at the expense of an increased rate of false positives.

This study has several limitations. First, the single-center, retrospective design may limit the model’s generalizability, and external validation using multicenter data was not conducted. Second, there is currently no standardized protocol for TCM pulse acquisition devices, and intraclass correlation coefficient analysis was not performed to evaluate the reproducibility of pulse wave measurements. Third, the absence of uniform and complete disease severity staging in medical records prevented patient stratification by severity level, which may impair the model’s ability to distinguish risk across clinical subgroups. Fourth, key clinical indicators recommended by the IWGDF guidelines—such as duration of diabetes, ankle–brachial index, and monofilament testing—were not systematically collected. This lack of comprehensive data may reduce the comparability of our model with established clinical assessment pathways, notably, diabetes duration was excluded due to incomplete retrospective data. Future studies should employ multicenter prospective designs, adopt standardized measurement protocols, and incorporate a broader range of clinical and body composition parameters to further validate and refine the predictive performance of the model.

In summary, the nomogram model constructed in this study for predicting DF risk in T2DM demonstrates strong predictive performance. It simplifies the risk assessment process and provides a practical tool for screening high-risk DF populations in primary healthcare settings. Future research could further integrate TCM diagnostic methods with electronic health records, artificial intelligence algorithms, and mobile health applications to establish an integrated “screening–warning–intervention” system for DF prevention and management, combining traditional Chinese and Western medicine approaches.

5 Conclusion

This study identified gender, BMI, HTN, SBP, BUN, and the pulse waveform parameter H5/H1 as independent predictors of T2DM patients complicated by DF. Based on these variables, this study developed a nomogram model that demonstrated high discriminative accuracy and clinical utility. By enabling personalized risk stratification, this tool may support more precise clinical decision – making for T2DM and DF. Prospective multicenter validation is warranted to confirm its generalizability and facilitate its integration into clinical practice.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

The study was reviewed and approved by the Ethics Committee of Shanghai Integrated Traditional Chinese and Western Medicine Hospital, affiliated with Shanghai University of Traditional Chinese Medicine. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

YH: Conceptualization, Data curation, Methodology, Writing – original draft, Investigation, Formal Analysis. WL: Formal Analysis, Methodology, Writing – original draft, Data curation, Conceptualization, Investigation. AX: Writing – review & editing, Resources, Data curation, Validation. SW: Software, Visualization, Writing – review & editing, Data curation. YW: Data curation, Writing – review & editing, Investigation. JX: Project administration, Supervision, Writing – review & editing. RG: Writing – review & editing, Software, Validation. CZ: Conceptualization, Funding acquisition, Supervision, Project administration, Writing – review & editing. HY: Supervision, Conceptualization, Project administration, Funding acquisition, Writing – review & editing.

Funding

The author(s) declared that financial support was received for this work and/or its publication. This work was supported by Shanghai Science and Technology Commission Science and Technology Plan Project (NO.21DZ2271000), National Administration of TCM High-Level Key Discipline Construction Project (zyyzdxk-2023069), Shanghai Science and Technology Commission Natural Science Foundation (NO.23ZR1460300).

Acknowledgments

This study was generously supported by Jingding Medical Tech, to whom we extend our sincere gratitude. We especially thank them for providing authorization and technical support for the JD_GBDR software. The team at Jingding Medical Tech offered invaluable assistance in data processing.

Conflict of interest

The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declared that generative AI was not used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

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/fendo.2026.1758234/full#supplementary-material

References

1. Saeedi P, Petersohn I, Salpea P, Malanda B, Karuranga S, Unwin N, et al. Global and regional diabetes prevalence estimates for 2019 and projections for 2030 and 2045: Results from the International Diabetes Federation Diabetes Atlas, 9th edition. Diabetes Res Clin Pract. (2019) 157:107843. doi: 10.1016/j.diabres.2019.107843

PubMed Abstract | Crossref Full Text | Google Scholar

2. Lu X, Xie Q, Pan X, Zhang R, Zhang X, Peng G, et al. Type 2 diabetes mellitus in adults: pathogenesis, prevention and therapy. Signal Transduct Target Ther. (2024) 9:262. doi: 10.1038/s41392-024-01951-9

PubMed Abstract | Crossref Full Text | Google Scholar

3. Raja JM, Maturana MA, Kayali S, Khouzam A, and Efeovbokhan N. Diabetic foot ulcer: A comprehensive review of pathophysiology and management modalities. World J Clin Cases. (2023) 11:1684–93. doi: 10.12998/wjcc.v11.i8.1684

PubMed Abstract | Crossref Full Text | Google Scholar

4. Hao Y, Cheng F, Pham M, Rein H, Patel D, Fang Y, et al. A noninvasive, economical, and instant-result method to diagnose and monitor type 2 diabetes using pulse wave: case-control study. JMIR Mhealth Uhealth. (2019) 7:e11959. doi: 10.2196/11959

PubMed Abstract | Crossref Full Text | Google Scholar

5. Lim J, Li J, Feng X, Feng L, Xiao X, Xia Y, et al. Machine learning-based evaluation of application value of traditional Chinese medicine clinical index and pulse wave parameters in the diagnosis of polycystic ovary syndrome. Eur J Integr Med. (2023) 64):102311. doi: 10.1016/j.eujim.2023.102311

Crossref Full Text | Google Scholar

6. Lyu Y, Huang W, Wu H, Hong J, Wang Y, Yan H, et al. A heart failure classification model from radial artery pulse wave using LSTM neural networks. BMC Med Inform Decis Mak. (2025) 25:318. doi: 10.1186/s12911-025-03167-5

PubMed Abstract | Crossref Full Text | Google Scholar

7. Yang J, Lü J, Qiu Z, Zhang M, and Yan H. Risk prediction of pulse wave for hypertensive target organ damage based on frequency-domain feature map. Med Eng Physics. (2024) 126):104161. doi: 10.1016/j.medengphy.2024.104161

PubMed Abstract | Crossref Full Text | Google Scholar

8. Chen J, Lin Y, Hu R, and Hu C. Prediction model for depression risk in middle-aged and elderly patients with metabolic syndrome: a nomogram and interpretable machine learning approach based on CHARLS. BMC Psychiatry. (2025) 25:987. doi: 10.1186/s12888-025-07434-7

PubMed Abstract | Crossref Full Text | Google Scholar

9. Ding J, Liu J, Cui X, Zhai Y, Zhang S, He J, et al. Development of a nomogram risk prediction model for subtherapeutic valproic acid plasma concentration: A retrospective cohort study. J Affect Disord. (2025) 393:120473. doi: 10.1016/j.jad.2025.120473

PubMed Abstract | Crossref Full Text | Google Scholar

10. Chinese Diabetes Society. Guidelines for the prevention and management of diabetes in China (2024 edition). Chin J Diabetes. (2025) 17:16–139.

Google Scholar

11. Zhang C, Liu L, Wu W, Wang Y, Yan H, Guo R, et al. Identifying coronary artery lesions by feature analysis of radial pulse wave: A case-control study. BioMed Res Int. (2021) 2021):5047501. doi: 10.1155/2021/5047501

PubMed Abstract | Crossref Full Text | Google Scholar

12. Yan J, Cai X, Chen S, Guo R, Yan H, and Wang Y. Ensemble learning-based pulse signal recognition: classification model development study. JMIR Med informatics. (2021) 9:e28039. doi: 10.2196/28039

PubMed Abstract | Crossref Full Text | Google Scholar

13. Tang Q, Chen Y, Liu S, Wang J, Zhang L, Qu H, et al. Pulse signal may be a key bridge connecting cardiology and pulse diagnosis of TCM. Heliyon. (2024) 10:e36785. doi: 10.1016/j.heliyon.2024.e36785

PubMed Abstract | Crossref Full Text | Google Scholar

14. Lim J, Li J, Feng X, Feng L, Xia Y, Xiao X, et al. Machine learning classification of polycystic ovary syndrome based on radial pulse wave analysis. BMC Complementary Med Therapies. (2023) 23:409. doi: 10.1186/s12906-023-04249-5

PubMed Abstract | Crossref Full Text | Google Scholar

15. Wu W, Zhang C, Ma X, Guo R, Yan J, Wang Y, et al. An early screening model of pulse detection technology for hepatic steatosis. Intelligent Med. (2023) 3:280–6. doi: 10.1016/j.imed.2023.03.002

Crossref Full Text | Google Scholar

16. Riley RD, Snell KI, Ensor J, Burke DL, Harrell Jr FE, Moons KG, et al. Minimum sample size for developing a multivariable prediction model: PART II - binary and time-to-event outcomes. Stat Med. (2019) 38:1276–96. doi: 10.1002/sim.7992

PubMed Abstract | Crossref Full Text | Google Scholar

17. Moons KG, de Groot JA, Bouwmeester W, Vergouwe Y, Mallett S, Altman DG, et al. Critical appraisal and data extraction for systematic reviews of prediction modelling studies: the CHARMS checklist. PloS Med. (2014) 11:e1001744. doi: 10.1371/journal.pmed.1001744

PubMed Abstract | Crossref Full Text | Google Scholar

18. Ye J, Wu Y, Yang S, Zhu D, Chen F, Chen J, et al. The global, regional and national burden of type 2 diabetes mellitus in the past, present and future: a systematic analysis of the Global Burden of Disease Study 2019. Front Endocrinol (Lausanne). (2023) 14:1192629. doi: 10.3389/fendo.2023.1192629

PubMed Abstract | Crossref Full Text | Google Scholar

19. Efthimiou O, Seo M, Chalkou K, Debray T, Egger M, and Salanti G. Developing clinical prediction models: a step-by-step guide. BMJ. (2024) 386:e078276. doi: 10.1136/bmj-2023-078276

PubMed Abstract | Crossref Full Text | Google Scholar

20. Slieker RC, van der Heijden AAWA, Siddiqui MK, Langendoen-Gort M, Nijpels G, Herings R, et al. Performance of prediction models for nephropathy in people with type 2 diabetes: systematic review and external validation study. BMJ. (2021) 374:n2134. doi: 10.1136/bmj.n2134

PubMed Abstract | Crossref Full Text | Google Scholar

21. Bandyk DF. The diabetic foot: Pathophysiology, evaluation, and treatment. Semin Vasc Surg. (2018) 31:43–8. doi: 10.1053/j.semvascsurg.2019.02.001

PubMed Abstract | Crossref Full Text | Google Scholar

22. Liao KM, Chang CW, Wang SH, Chang YT, Chen YC, and Wang GC. The first harmonic of radial pulse wave predicts major adverse cardiovascular and microvascular events in patients with type 2 diabetes. J Diabetes Complications. (2019) 33:107420. doi: 10.1016/j.jdiacomp.2019.107420

PubMed Abstract | Crossref Full Text | Google Scholar

23. Fei ZF. Modern pulse diagnosis in traditional Chinese medicine. 2nd ed. Beijing: People’s Medical Publishing House (2015) p. 152–5.

Google Scholar

24. Reference Values for Arterial Stiffness’ Collaboration. Determinants of pulse wave velocity in healthy people and in the presence of cardiovascular risk factors: ‘establishing normal and reference values’. Eur Heart J. (2010) 31:2338–50. doi: 10.1093/eurheartj/ehq165

PubMed Abstract | Crossref Full Text | Google Scholar

25. Liang Y, Yang Y, Chen P, Zhou T, and Wang Z. Correlation of arterial total compliance and carotid-femoral pulse wave velocity: an echocardiographic study. Chin J Med Ultrasound (Electronic Edition). (2017) 14:671–5.

Google Scholar

26. Su Y, Lin H, Chen S, and Xu Y. Clinical characteristics of lower extremity arterial disease in newly diagnosed diabetic patients. J Vasc Endovascular Surg. (2024) 10:1356–60. doi: 10.19418/j.cnki.issn2096-0646

Crossref Full Text | Google Scholar

27. Rümenapf G, Abilmona N, Morbach S, and Sigl M. Peripheral arterial disease and the diabetic foot syndrome: neuropathy makes the difference! A narrative review. J Clin Med. (2024) 13:2141. doi: 10.3390/jcm13072141

PubMed Abstract | Crossref Full Text | Google Scholar

28. Jia M and Gong W. Exploring pulse characteristics and their correlation with symptoms in type 2 diabetes mellitus with qi deficiency and phlegm turbidity pattern based on pressure pulse wave analysis. Modern J Integrated Traditional Chin Western Med. (2025) 34:24–31.

Google Scholar

29. ElSayed NA, Aleppo G, Aroda VR, Bannuru RR, Brown FM, Bruemmer D, et al. 10. Cardiovascular disease and risk management: standards of care in diabetes-2023. Diabetes Care. (2023) 46:S158–90. doi: 10.2337/dc23-S010

PubMed Abstract | Crossref Full Text | Google Scholar

30. Basu S, Sussman JB, Rigdon J, Steimle L, Denton BT, and Hayward RA. Benefit and harm of intensive blood pressure treatment: Derivation and validation of risk models using data from the SPRINT and ACCORD trials. PloS Med. (2017) 14:e1002410. doi: 10.1371/journal.pmed.1002410

PubMed Abstract | Crossref Full Text | Google Scholar

31. Buckley LF, Dixon DL, Wohlford GF 4th, Wijesinghe DS, Baker WL, and Van Tassell BW. Intensive versus standard blood pressure control in SPRINT-eligible participants of ACCORD-BP. Diabetes Care. (2017) 40:1733–8. doi: 10.2337/dc17-1366

PubMed Abstract | Crossref Full Text | Google Scholar

32. Pinzur M, Freeland R, and Juknelis D. The association between body mass index and foot disorders in diabetic patients. Foot Ankle Int. (2005) 26:375–7. doi: 10.1177/107110070502600506

PubMed Abstract | Crossref Full Text | Google Scholar

33. Carnethon MR, Rasmussen-Torvik LJ, and Palaniappan L. The obesity paradox in diabetes. Curr Cardiol Rep. (2014) 16:446. doi: 10.1007/s11886-013-0446-3

PubMed Abstract | Crossref Full Text | Google Scholar

34. Lin CC, Li CI, Liu CS, Lin WY, Lin CH, Chiang JI, et al. Obesity paradox in associations between body mass index and diabetes-related hospitalization and mortality in patients with type 2 diabetes: Retrospective cohort studies. Diabetes Metab. (2019) 45:564–72. doi: 10.1016/j.diabet.2019.02.007

PubMed Abstract | Crossref Full Text | Google Scholar

35. Zhou Y, Wang J, Yao Q, Jian Q, and Luo Z. Prevalence of sarcopenic obesity in patients with diabetes and adverse outcomes: A systematic review and meta-analysis. Clin Nutr ESPEN. (2023) 58:128–35. doi: 10.1016/j.clnesp.2023.09.920

PubMed Abstract | Crossref Full Text | Google Scholar

36. Liu L, Ban C, Jia S, Chen X, and He S. Association of predicted fat mass, predicted lean mass and predicted percent fat with diabetes mellitus in Chinese population: a 15-year prospective cohort. BMJ Open. (2022) 12:e058162. doi: 10.1136/bmjopen-2021-058162

PubMed Abstract | Crossref Full Text | Google Scholar

37. Tang J, Cai X, Liu A, Yu N, and Wang S. Association between predicted fat mass, predicted lean mass, predicted percent fat and type 2 diabetes mellitus in Japanese adults: a retrospective study. BMC Endocr Disord. (2024) 24:48. doi: 10.1186/s12902-024-01579-4

PubMed Abstract | Crossref Full Text | Google Scholar

38. Pan XF, Wang L, and Pan A. Epidemiology and determinants of obesity in China. Lancet Diabetes Endocrinol. (2021) 9:373–92. doi: 10.1016/S2213-8587(21)00045-0

PubMed Abstract | Crossref Full Text | Google Scholar

39. Renard E. Bariatric surgery in patients with late-stage type 2 diabetes: expected beneficial effects on risk ratio and outcomes. Diabetes Metab. (2009) 35:564–8. doi: 10.1016/S1262-3636(09)73467-6

PubMed Abstract | Crossref Full Text | Google Scholar

40. andforth L, Raverdy V, Sandforth A, Bauvin P, Chatelain E, Verkindt H, et al. Subphenotype-dependent benefits of bariatric surgery for individuals at risk for type 2 diabetes. Diabetes Care. (2025) 48:996–1006. doi: 10.2337/dc25-0160

PubMed Abstract | Crossref Full Text | Google Scholar

41. Lee DH, Keum N, Hu FB, Orav EJ, Rimm EB, Willett WC, et al. Comparison of the association of predicted fat mass, body mass index, and other obesity indicators with type 2 diabetes risk: two large prospective studies in US men and women. Eur J Epidemiol. (2018) 33:1113–23. doi: 10.1007/s10654-018-0433-5

PubMed Abstract | Crossref Full Text | Google Scholar

42. Ciarambino T, Crispino P, Leto G, Mastrolorenzo E, Para O, and Giordano M. Influence of gender in diabetes mellitus and its complication. Int J Mol Sci. (2022) 23:8850. doi: 10.3390/ijms23168850

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: diabetic foot, nomogram, predictive model, pulse mapping parameters, type 2 diabetes mellitus

Citation: Huo Y, Li W, Xu A, Wang S, Wang Y, Xu J, Guo R, Zhao C and Yan H (2026) Development and validation of a nomogram model integrating noninvasive detection of radial pulse wave for predicting diabetic foot risk in type 2 diabetes mellitus. Front. Endocrinol. 17:1758234. doi: 10.3389/fendo.2026.1758234

Received: 01 December 2025; Accepted: 09 January 2026; Revised: 02 January 2026;
Published: 12 February 2026.

Edited by:

Arvind Mukundan, National Chung Cheng University, Taiwan

Reviewed by:

Shihua Shi, Friedrich Miescher Institute for Biomedical Research (FMI), Switzerland
Andrii Sverstiuk, Ternopil State Medical University, Ukraine

Copyright © 2026 Huo, Li, Xu, Wang, Wang, Xu, Guo, Zhao and Yan. 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: Cheng Zhao, Y2hlbmd6aGFvXzc5QDE2My5jb20=; Haixia Yan, aGp5MjAxMm9rQDE2My5jb20=

These authors have contributed equally to this work

Disclaimer: 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.