<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Publishing DTD v1.3 20210610//EN" "JATS-journalpublishing1-3-mathml3.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns:ali="http://www.niso.org/schemas/ali/1.0/" article-type="research-article" dtd-version="1.3" xml:lang="EN">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Med.</journal-id>
<journal-title-group>
<journal-title>Frontiers in Medicine</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Med.</abbrev-journal-title>
</journal-title-group>
<issn pub-type="epub">2296-858X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fmed.2026.1778003</article-id>
<article-version article-version-type="Version of Record" vocab="NISO-RP-8-2008"/>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Original Research</subject>
</subj-group>
</article-categories>
<title-group>
<article-title>Dynamic C-reactive protein trajectories predict prolonged healing time in diabetic wounds: a machine learning model based on a single-center cohort with standardized wound size</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" equal-contrib="yes">
<name>
<surname>Jiang</surname>
<given-names>Sichao</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="author-notes" rid="fn0001"><sup>&#x2020;</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/2611656"/>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Data curation" vocab-term-identifier="https://credit.niso.org/contributor-roles/data-curation/">Data curation</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; original draft" vocab-term-identifier="https://credit.niso.org/contributor-roles/writing-original-draft/">Writing &#x2013; original draft</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="conceptualization" vocab-term-identifier="https://credit.niso.org/contributor-roles/conceptualization/">Conceptualization</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Formal analysis" vocab-term-identifier="https://credit.niso.org/contributor-roles/formal-analysis/">Formal analysis</role>
</contrib>
<contrib contrib-type="author" equal-contrib="yes">
<name>
<surname>Song</surname>
<given-names>Qixuan</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="author-notes" rid="fn0001"><sup>&#x2020;</sup></xref>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="conceptualization" vocab-term-identifier="https://credit.niso.org/contributor-roles/conceptualization/">Conceptualization</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Data curation" vocab-term-identifier="https://credit.niso.org/contributor-roles/data-curation/">Data curation</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Formal analysis" vocab-term-identifier="https://credit.niso.org/contributor-roles/formal-analysis/">Formal analysis</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; original draft" vocab-term-identifier="https://credit.niso.org/contributor-roles/writing-original-draft/">Writing &#x2013; original draft</role>
</contrib>
<contrib contrib-type="author" equal-contrib="yes">
<name>
<surname>Wang</surname>
<given-names>Junhuan</given-names>
</name>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<xref ref-type="author-notes" rid="fn0001"><sup>&#x2020;</sup></xref>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Formal analysis" vocab-term-identifier="https://credit.niso.org/contributor-roles/formal-analysis/">Formal analysis</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; original draft" vocab-term-identifier="https://credit.niso.org/contributor-roles/writing-original-draft/">Writing &#x2013; original draft</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Data curation" vocab-term-identifier="https://credit.niso.org/contributor-roles/data-curation/">Data curation</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="conceptualization" vocab-term-identifier="https://credit.niso.org/contributor-roles/conceptualization/">Conceptualization</role>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Sun</surname>
<given-names>Xiaohong</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="investigation" vocab-term-identifier="https://credit.niso.org/contributor-roles/investigation/">Investigation</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="software" vocab-term-identifier="https://credit.niso.org/contributor-roles/software/">Software</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="methodology" vocab-term-identifier="https://credit.niso.org/contributor-roles/methodology/">Methodology</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; review &#x0026; editing" vocab-term-identifier="https://credit.niso.org/contributor-roles/writing-review-editing/">Writing &#x2013; review &#x0026; editing</role>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Zhang</surname>
<given-names>Zhen</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="corresp" rid="c001"><sup>&#x002A;</sup></xref>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="visualization" vocab-term-identifier="https://credit.niso.org/contributor-roles/visualization/">Visualization</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="resources" vocab-term-identifier="https://credit.niso.org/contributor-roles/resources/">Resources</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; review &#x0026; editing" vocab-term-identifier="https://credit.niso.org/contributor-roles/writing-review-editing/">Writing &#x2013; review &#x0026; editing</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="validation" vocab-term-identifier="https://credit.niso.org/contributor-roles/validation/">Validation</role>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Wang</surname>
<given-names>Shouyu</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="corresp" rid="c001"><sup>&#x002A;</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/3333398"/>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="conceptualization" vocab-term-identifier="https://credit.niso.org/contributor-roles/conceptualization/">Conceptualization</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; review &#x0026; editing" vocab-term-identifier="https://credit.niso.org/contributor-roles/writing-review-editing/">Writing &#x2013; review &#x0026; editing</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Funding acquisition" vocab-term-identifier="https://credit.niso.org/contributor-roles/funding-acquisition/">Funding acquisition</role>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Zong</surname>
<given-names>Junwei</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="corresp" rid="c001"><sup>&#x002A;</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/2212710"/>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; review &#x0026; editing" vocab-term-identifier="https://credit.niso.org/contributor-roles/writing-review-editing/">Writing &#x2013; review &#x0026; editing</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="supervision" vocab-term-identifier="https://credit.niso.org/contributor-roles/supervision/">Supervision</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="software" vocab-term-identifier="https://credit.niso.org/contributor-roles/software/">Software</role>
</contrib>
</contrib-group>
<aff id="aff1"><label>1</label><institution>Department of Orthopaedic Surgery, The First Affiliated Hospital, Dalian Medical University</institution>, <city>Dalian</city>, <country country="cn">China</country></aff>
<aff id="aff2"><label>2</label><institution>Department of Obstetrical, Affiliated Hospital, Qingdao University</institution>, <city>Qingdao</city>, <country country="cn">China</country></aff>
<author-notes>
<corresp id="c001"><label>&#x002A;</label>Correspondence: Junwei Zong, <email xlink:href="mailto:aweizone@163.com">aweizone@163.com</email>; Shouyu Wang, <email xlink:href="mailto:wangshouyu666@126.com">wangshouyu666@126.com</email> Zhen Zhang, <email xlink:href="mailto:doctorzz@126.com">doctorzz@126.com</email></corresp>
<fn fn-type="equal" id="fn0001">
<label>&#x2020;</label>
<p>These authors have contributed equally to this work and share first authorship</p>
</fn>
</author-notes>
<pub-date publication-format="electronic" date-type="pub" iso-8601-date="2026-02-12">
<day>12</day>
<month>02</month>
<year>2026</year>
</pub-date>
<pub-date publication-format="electronic" date-type="collection">
<year>2026</year>
</pub-date>
<volume>13</volume>
<elocation-id>1778003</elocation-id>
<history>
<date date-type="received">
<day>30</day>
<month>12</month>
<year>2025</year>
</date>
<date date-type="rev-recd">
<day>19</day>
<month>01</month>
<year>2026</year>
</date>
<date date-type="accepted">
<day>22</day>
<month>01</month>
<year>2026</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x00A9; 2026 Jiang, Song, Wang, Sun, Zhang, Wang and Zong.</copyright-statement>
<copyright-year>2026</copyright-year>
<copyright-holder>Jiang, Song, Wang, Sun, Zhang, Wang and Zong</copyright-holder>
<license>
<ali:license_ref start_date="2026-02-12">https://creativecommons.org/licenses/by/4.0/</ali:license_ref>
<license-p>This is an open-access article distributed under the terms of the <ext-link ext-link-type="uri" xlink:href="https://creativecommons.org/licenses/by/4.0/">Creative Commons Attribution License (CC BY)</ext-link>. 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.</license-p>
</license>
</permissions>
<abstract>
<sec>
<title>Objective</title>
<p>To develop a machine learning (ML) model for predicting prolonged healing (&#x003E;8&#x202F;weeks) in diabetic wounds, focusing on dynamic C-reactive protein (CRP) trajectories.</p>
</sec>
<sec>
<title>Methods</title>
<p>This was a retrospective single-center cohort study. We included 465 patients with type 2 diabetes, standardized wound sizes (5&#x2013;8&#x202F;cm<sup>2</sup>), and debridement alone (2021&#x2013;2024: training set, <italic>n</italic>&#x202F;=&#x202F;325; 2025: temporal validation set, <italic>n</italic>&#x202F;=&#x202F;140). Serial CRP was measured at admission (CRP), post-antibiotic preoperatively (CRP_2nd), and postoperatively at discharge (CRP_3rd). Therapeutic response variables (therapeutic_response_1/2/all) were calculated as percentage changes in serial CRP levels across treatment phases, reflecting anti-inflammatory/antimicrobial efficacy. LASSO regression selected features, 12 ML models were constructed, and performance was evaluated via AUC, sensitivity, and specificity. SHAP analysis interpreted predictions.</p>
</sec>
<sec>
<title>Results</title>
<p>The GradientBoosting model exhibited superior performance (validation set: accuracy&#x202F;=&#x202F;0.9357, sensitivity&#x202F;=&#x202F;0.8689, specificity&#x202F;=&#x202F;0.9873). LASSO regression identified 15 key variables [including CRP_2nd, CRP_3rd, albumin (ALB)]. SHAP analysis revealed CRP_2nd as the most influential predictor (mean absolute SHAP value&#x202F;=&#x202F;0.460), with elevated CRP_2nd/CRP_3rd associated with prolonged healing and higher ALB/favorable therapeutic responses as protective factors.</p>
</sec>
<sec>
<title>Conclusion</title>
<p>Dynamic CRP trajectories, particularly CRP_2nd, are critical for predicting prolonged diabetic wound healing. The GradientBoosting model provides a clinically actionable tool for risk stratification.</p>
</sec>
</abstract>
<kwd-group>
<kwd>cohort analysis</kwd>
<kwd>CRP dynamics</kwd>
<kwd>diabetes mellitus</kwd>
<kwd>machine learning model</kwd>
<kwd>wound prognosis</kwd>
</kwd-group>
<funding-group>
<funding-statement>The author(s) declared that financial support was not received for this work and/or its publication.</funding-statement>
</funding-group>
<counts>
<fig-count count="6"/>
<table-count count="0"/>
<equation-count count="1"/>
<ref-count count="37"/>
<page-count count="13"/>
<word-count count="9469"/>
</counts>
<custom-meta-group>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Pathology</meta-value>
</custom-meta>
</custom-meta-group>
</article-meta>
</front>
<body>
<sec sec-type="intro" id="sec1">
<label>1</label>
<title>Introduction</title>
<p>C-reactive protein (CRP), produced by hepatocytes in response to proinflammatory cytokines like interleukin-6 (IL-6) and tumor necrosis factor-&#x03B1; (TNF-&#x03B1;), is a well-established biomarker of systemic inflammation. Its rapid elevation during infection, tissue injury, or inflammation, followed by a decline with effective treatment, makes it a cornerstone of clinical monitoring for disease activity and therapeutic response (<xref ref-type="bibr" rid="ref1">1</xref>). In the context of wound healing, CRP dynamics reflect the balance between inflammatory clearance and tissue repair&#x2014;an equilibrium particularly fragile in chronic wounds, where persistent inflammation often impedes progression to the proliferative and remodeling phases (<xref ref-type="bibr" rid="ref2">2</xref>).</p>
<p>Diabetic wounds, such as foot and lower extremity skin ulcers, are a significant complication of diabetes mellitus, impacting up to 25% of patients during their lifetime (<xref ref-type="bibr" rid="ref3">3</xref>). These wounds exhibit impaired angiogenesis, neuropathy, and dysregulated inflammation, resulting in extended healing durations, frequent recurrences, and an elevated risk of amputation (<xref ref-type="bibr" rid="ref4">4</xref>). The clinical burden is substantial: even small-to-moderate diabetic wounds (5&#x2013;8&#x202F;cm<sup>2</sup>) can require weeks to months to heal, with approximately 30%&#x2013;40% of cases progressing to prolonged healing (&#x003E;8&#x202F;weeks) despite standard care (<xref ref-type="bibr" rid="ref5">5</xref>, <xref ref-type="bibr" rid="ref6">6</xref>). This unpredictability not only increases healthcare costs but also diminishes patients&#x2019; quality of life, underscoring the need for reliable tools to identify high-risk individuals early.</p>
<p>Existing research on diabetic wound healing has primarily focused on static baseline factors, such as initial wound size, depth, and comorbidities like peripheral artery disease (<xref ref-type="bibr" rid="ref7">7</xref>). While these factors provide valuable context, they fail to capture the dynamic biological processes that drive healing outcomes. In recent years, studies have begun to explore inflammatory markers as predictors: for example, elevated baseline CRP levels have been associated with prolonged healing (<xref ref-type="bibr" rid="ref8">8</xref>), and serial measurements of proinflammatory cytokines (e.g., IL-6, TNF-&#x03B1;) have shown promise in assessing treatment response (<xref ref-type="bibr" rid="ref9">9</xref>). However, these studies often lack standardization in wound size, include heterogeneous patient populations (e.g., mixing minor debridement with major amputations), or rely on single-time-point measurements, limiting their clinical applicability.</p>
<p>The importance of dynamic inflammatory trajectories&#x2014;rather than isolated values&#x2014;has gained recognition in other chronic conditions, such as sepsis and rheumatoid arthritis, where patterns of marker decline (e.g., CRP &#x201C;clearance rate&#x201D;) strongly predict outcomes (<xref ref-type="bibr" rid="ref10">10</xref>, <xref ref-type="bibr" rid="ref11">11</xref>). In diabetic wounds, this principle is equally relevant: an initial elevation in CRP may indicate a productive inflammatory response to clear pathogens, but a failure to decline after interventions (e.g., antibiotics, debridement) signals persistent infection or impaired healing mechanisms (<xref ref-type="bibr" rid="ref12">12</xref>). Surprisingly, few studies have systematically analyzed serial CRP measurements to predict prolonged healing, especially in well-characterized cohorts with standardized wound sizes and uniform surgical interventions (e.g., debridement alone, excluding amputations).</p>
<p>Machine learning (ML) has become a valuable asset in wound healing research, surpassing traditional regression models by effectively capturing nonlinear relationships and interactions among multiple variables. Existing ML models for diabetic wounds have integrated features such as wound imagery, microbiome data, and baseline clinical metrics to predict healing time or amputation risk (<xref ref-type="bibr" rid="ref13 ref14 ref15">13&#x2013;15</xref>). For instance, a recent study using random forest algorithms achieved an AUC of 0.82 for predicting 12-week healing in diabetic foot ulcers, leveraging features like wound area and HbA1c (<xref ref-type="bibr" rid="ref13">13</xref>). However, these models rarely incorporate dynamic inflammatory markers, and none have specifically focused on CRP trajectories across critical treatment time points (e.g., post-antibiotic, pre-operative, and post-operative phases).</p>
<p>Our study seeks to create a machine learning model to predict diabetic wound healing times exceeding 8&#x202F;weeks, emphasizing the dynamic patterns of CRP. By restricting the cohort to patients with standardized wound sizes (5&#x2013;8&#x202F;cm<sup>2</sup>) who underwent only debridement (excluding amputations), we minimize confounding by wound severity and surgical complexity. We hypothesize that serial CRP measurements&#x2014;specifically the magnitude and rate of decline after antibiotic treatment and surgery&#x2014;will outperform static baseline factors in predicting prolonged healing. This model will not only enhance risk stratification but also highlight critical intervention windows, enabling clinicians to adjust anti-inflammatory and antimicrobial strategies proactively. Our objective is to develop a clinically useful tool to enhance outcomes for patients with diabetic wounds. While we utilized a temporal validation cohort to assess model robustness, the single-center design warrants future multi-center replication to confirm generalizability across diverse clinical settings.</p>
</sec>
<sec sec-type="materials|methods" id="sec2">
<label>2</label>
<title>Materials and methods</title>
<sec id="sec3">
<label>2.1</label>
<title>Data source</title>
<p>This retrospective cohort study utilized clinical data from patients at the First Affiliated Hospital of Dalian Medical University. The training dataset comprised patients hospitalized from January 2021 to December 2024, and a temporal validation dataset included patients admitted in 2025 to evaluate the model&#x2019;s generalizability within the same institutional setting. Notably, this constitutes temporal validation (time-separated, single-center) rather than multi-center external validation, which will be addressed in future research. This retrospective cohort study was conducted in accordance with the Helsinki Declaration of 1975, as revised in 2024. The study protocol was approved by the Ethics Committee of the First Affiliated Hospital of Dalian Medical University (Approval No.: PJ-KS-KY-2025-440). Due to the retrospective nature and use of de-identified data, the requirement for informed consent was waived by the Ethics Committee, with strict adherence to patient anonymity and data protection principles.</p>
</sec>
<sec id="sec4">
<label>2.2</label>
<title>Participants</title>
<p>The patient selection process is illustrated in <xref ref-type="fig" rid="fig1">Figure 1</xref>.</p>
<fig position="float" id="fig1">
<label>Figure 1</label>
<caption>
<p>Patient selection flowchart for retrospective cohort study on diabetic wound healing.</p>
</caption>
<graphic xlink:href="fmed-13-1778003-g001.tif" mimetype="image" mime-subtype="tiff">
<alt-text content-type="machine-generated">Flowchart showing the screening and selection of patients for a study. Two groups are shown: patients from 2021-2024 (n=18,411) and from 2025 (n=3,491) at the First Affiliated Hospital of Dalian Medical University. Each group was screened for type 2 diabetes and wound size, leading to subsets with diabetes-related wounds who underwent debridement. The 2021-2024 group formed a training set (n=325), and the 2025 group formed an external validation set (n=140). Each step includes patient counts and criteria details.</alt-text>
</graphic>
</fig>
<sec id="sec5">
<label>2.2.1</label>
<title>Inclusion criteria</title>
<p>Met WHO criteria for type 2 diabetes mellitus, defined by fasting blood glucose &#x2265;7.0&#x202F;mmol/L or HbA1c&#x202F;&#x2265;&#x202F;6.5%, excluding type 1, gestational, or secondary diabetes.</p>
<p>Admitted for diabetes-related wounds (e.g., diabetic foot ulcers, lower extremity skin ulcers) with an initial wound area of 5&#x2013;8&#x202F;cm<sup>2</sup>, Wound area was measured by two experienced attending physicians independently using ImageJ software (version 1.8.0), both of which were uniformly trained before measurement, using the same image acquisition and measurement criteria (e.g., the area not covered by epithelization of the wound edge was used as the measurement range). The intraclass correlation coefficient (ICC) was used to evaluate the measurement consistency, and the ICC was 0.96 (95% CI: 0.94&#x2013;0.98), indicating that the measurement results of the two reviewers were highly consistent;</p>
<p>Underwent debridement only (excluding major surgical interventions such as amputation, digit resection, or skin flap transplantation);</p>
<p>Completed C-reactive protein (CRP) testing at three time points during hospitalization:</p>
<list list-type="simple">
<list-item>
<p>&#x2460; At admission (within 24&#x202F;h of admission, denoted as CRP);</p>
</list-item>
<list-item>
<p>&#x2461; After antibiotic treatment but before surgery (typically 3&#x2013;7&#x202F;days after admission, denoted as CRP_2nd);</p>
</list-item>
<list-item>
<p>&#x2462; At discharge after surgical treatment (wound not fully healed but meeting discharge criteria, denoted as CRP_3rd);</p>
</list-item>
</list>
<p>CRP_2nd was measured from 3 to 7&#x202F;days after admission (median 4&#x202F;days, interquartile range: 3 to 5&#x202F;days) and was strictly recorded and included in the database for all patients. One-way analysis of variance was used to test the correlation between the measurement time (3, 4, 5, 6&#x2013;7&#x202F;days) and CRP_2nd values. The results showed that there was no statistically significant difference in CRP_2nd levels at different time points (<italic>F</italic>&#x202F;=&#x202F;1.23, <italic>p</italic>&#x202F;=&#x202F;0.301). This indicates that measurement variability within this time window did not significantly influence the clinical significance of CRP_2nd.</p>
<p>Aged &#x2265;18&#x202F;years with complete clinical data and at least 12&#x202F;weeks of follow-up to record the final healing time.</p>
</sec>
<sec id="sec6">
<label>2.2.2</label>
<title>Exclusion criteria</title>
<p>Complicated with malignant tumors, severe immunodeficiency diseases (e.g., AIDS), or end-stage liver/kidney failure (factors independently affecting wound healing);</p>
<p>Death or transfer to another hospital due to non-wound-related causes during treatment (inability to obtain complete healing time);</p>
<p>Abandonment of treatment or loss to follow-up before wound healing;</p>
<p>Presence of other acute diseases affecting inflammatory markers (e.g., acute myocardial infarction, pneumonia).</p>
</sec>
</sec>
<sec id="sec7">
<label>2.3</label>
<title>Variable selection and definition of target variables</title>
<sec id="sec8">
<label>2.3.1</label>
<title>Collected variables</title>
<p>Variables were selected based on three principles: (1) Clinical relevance: indicators known to be associated with diabetic wound healing (e.g., inflammatory markers, nutritional markers, metabolic markers, etc.) were included; (2) Data availability: all included variables were clinical routine indicators to ensure the clinical practicability of the model. (3) Literature support: refer to the variable selection experience of previous diabetic wound prediction models. Indicators such as wound culture results were not included because some patients lacked complete culture data (missing rate &#x003E; 30%), which may introduce bias. Rare comorbidities, such as rare autoimmune diseases, were not included, which had limited contribution to model prediction because their incidence was less than 1%.</p>
<p>Clinical data were obtained from electronic medical records, encompassing;</p>
<p>Demographic information: gender, age, smoking history, drinking history;</p>
<p>Biochemical indicators: The study examines various biomarkers, including low-density lipoprotein cholesterol (LDL-C), triglycerides (TG), albumin (ALB), albumin/globulin ratio (ALB/GLO), direct bilirubin (DBIL), alkaline phosphatase (ALP), glucose (Glu), high-sensitivity cardiac troponin I (hs-cTnI), procalcitonin (PCT), high-density lipoprotein cholesterol (HDL-C), D-dimer, and several hematological parameters such as mean corpuscular volume (MCV), mean corpuscular hemoglobin (MCH), mean corpuscular hemoglobin concentration (MCHC), mean platelet volume (MPV), red blood cell distribution width (RDW-CV and RDW-SD), hematocrit (Hct), red blood cell count (RBC), platelet count (PLT), hemoglobin (Hb), and electrolytes including chloride (Cl<sup>&#x2212;</sup>), sodium (Na<sup>+</sup>), potassium (K<sup>+</sup>), magnesium (Mg<sup>2+</sup>), and glycated hemoglobin (HbA1c).</p>
<p>Hematological indicators include neutrophil percentage (Neut%), neutrophil count (Neut#), monocyte percentage (Mono%), monocyte count (Mono#), basophil percentage (Baso%), basophil count (Baso#), eosinophil percentage (Eos%), eosinophil count (Eos#), lymphocyte percentage (Lymph%), lymphocyte count (Lymph#), white blood cell count (WBC), erythrocyte sedimentation rate (ESR), and globulin (GLO).</p>
<p>Inflammatory dynamic indicators: CRP (measured at admission within 24&#x202F;h), CRP_2nd (measured 3&#x2013;7&#x202F;days after admission, post-antibiotic preoperatively), CRP_3rd (measured at discharge postoperatively), therapeutic response at first stage (therapeutic_response_1), therapeutic response at second stage (therapeutic_response_2), overall therapeutic response (therapeutic_response_all). All CRP measurements were performed using a standardized nephelometric assay (Siemens BNII System) with consistent reagents and quality control protocols throughout the study period (2021&#x2013;2025). The inter-assay coefficient of variation (CV) was &#x003C; 5%, and the lower limit of detection was 0.1&#x202F;mg/L, ensuring reliability and comparability of serial measurements.</p>
<p>Data preprocessing flow: (1) Missing value handling: Variables with a missing rate &#x003C;5% were imputed using multiple imputation by chained equations (MICE) to preserve data integrity; variables with a missing rate &#x003E;5% were excluded to avoid bias from excessive imputation. This 5% threshold was selected based on established practices in clinical prediction model research and validated via sensitivity analysis: imputing variables with missing rates &#x003E; 5% increased training set accuracy to 0.981 but reduced validation set accuracy to 0.887 (vs. 0.9357 for the final model), indicating heightened overfitting risk; (2) Outlier handling: Identify outliers using the IQR method and handle them through Winsorization (scaling to the 99th percentile) to avoid interference from extreme values; (3) Feature standardization: Perform <italic>Z</italic>-score standardization on continuous variables to eliminate dimensional differences; (4) Categorical variable encoding: Use one-hot encoding to handle binary categorical variables such as gender and smoking history.</p>
</sec>
<sec id="sec9">
<label>2.3.2</label>
<title>Definition of therapeutic response variables</title>
<p>Therapeutic response variables were rule-based metrics derived from serial CRP changes to quantify treatment efficacy across critical phases: Therapeutic_response_1: Percentage change in CRP from admission (CRP) to post-antibiotic preoperatively (CRP_2nd). A positive value indicates a decrease in CRP (favorable response), while a negative value indicates an increase (unfavorable response). Therapeutic_response_2: Percentage change in CRP from post-antibiotic preoperatively (CRP_2nd) to discharge (CRP_3rd). Therapeutic_response_all: Composite therapeutic response, defined as the weighted average of therapeutic_response_1 and therapeutic_response_2, with weights proportional to the duration of the antibiotic treatment phase (3&#x2013;7&#x202F;days) and postoperative recovery phase (until discharge).</p>
</sec>
<sec id="sec10">
<label>2.3.3</label>
<title>Definition of target variable</title>
<p>The primary outcome was &#x201C;prolonged healing,&#x201D; defined as failure to achieve complete wound healing (full epithelialization without exudate) within 8&#x202F;weeks (56&#x202F;days) after initial debridement. Healing status was confirmed through 12-week follow-up (in-person or telephone interviews).</p>
</sec>
</sec>
<sec id="sec11">
<label>2.4</label>
<title>Variable screening using LASSO regression</title>
<p>LASSO regression, a regularized linear regression approach that is extensively utilized for variable selection, conducts variable selection and complexity adjustment during the fitting of generalized linear models. In this research, LASSO regression was employed for variable screening. The optimal penalty parameter &#x03BB; of LASSO regression was determined by 10-fold cross-validation. The classification error was used as the objective function, and the value of &#x03BB; corresponding to the minimum error (&#x03BB;&#x202F;=&#x202F;0.012) was selected. The formula for LASSO regression is presented below:</p>
<disp-formula id="E1">
<mml:math id="M1">
<mml:munder>
<mml:mo>min</mml:mo>
<mml:mi>&#x03B2;</mml:mi>
</mml:munder>
<mml:mo stretchy="true">{</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mi>y</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>X</mml:mi>
<mml:msubsup>
<mml:mi>&#x03B2;</mml:mi>
<mml:mn>2</mml:mn>
<mml:mn>2</mml:mn>
</mml:msubsup>
<mml:mo>+</mml:mo>
<mml:mi>&#x03BB;</mml:mi>
<mml:msub>
<mml:mi>&#x03B2;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo stretchy="true">}</mml:mo>
</mml:math>
</disp-formula>
<p>In the equation, y denotes the response variable vector indicating whether healing is prolonged. X is the design matrix with observations of p explanatory variables for n samples. <italic>&#x03B2;</italic> is the regression coefficient vector, indicating each explanatory variable&#x2019;s effect on the response variable. <italic>&#x03BB;</italic> is the parameter for regularization that determines the intensity of the penalty term. The chosen variables served as candidate key variables for subsequent model development.</p>
</sec>
<sec id="sec12">
<label>2.5</label>
<title>Construction of machine learning models</title>
<p>Twelve machine learning algorithms were utilized: RandomForest, GradientBoosting, SVM_Kernel (support vector machine with RBF kernel), LogisticModel (logistic regression), NeighborMethod (k-nearest neighbor, <italic>k</italic>&#x202F;=&#x202F;5), PLSModel (partial least squares regression), BoostingMethod (gradient boosting decision tree, GBDT), NeuralNet (3-layer feedforward neural network), BayesMethod (naive Bayes), DiscriminantModel (linear discriminant analysis), LASSO (LASSO regression), and AdaptiveBoosting (AdaBoost with decision tree weak learners).</p>
<p>Selection criteria of 12 machine learning models: Referring to the commonly used algorithms in the prediction of diabetic complications in the past, it covers four categories: ensemble learning (GradientBoosting, RandomForest, etc.), traditional statistical models (Logistic regression, PLS), kernel methods (SVM_Kernel), and neural networks. The representativeness and coverage of model types were ensured to avoid algorithm selection bias.</p>
<p>Model optimization strategy: (1) Feature level optimization: LASSO regression was used to screen key variables and eliminate redundant features to reduce model calculation and overfitting risk. (2) Optimization at the algorithm level: the GradientBoosting model adopts the gradient boosting framework to reduce the prediction error by iterative fitting of the residual, and combines the regularization term (tree depth limit, subsample proportion) to control the model complexity. (3) Optimization at the engineering level: model quantization compression technology was used to convert floating-point number parameters into low-precision formats to improve inference speed without loss of performance. (4) Hyperparameter level optimization: the optimal parameter combination was found through grid search and cross validation to avoid the model being trapped in the local optimal solution.</p>
<p>In order to reduce the risk of overfitting caused by high-dimensional variables, in addition to LASSO regularization, 5-fold cross validation was further used to optimize the hyperparameters of all models, and the generalization ability of the models was evaluated by an temporal validation set (2025 independent cohort). The results showed that the optimal model (GradientBoosting) had a small performance difference between the training set and the validation set (accuracy rates were 97.54 and 93.57%, respectively), indicating that the model did not overfit significantly.</p>
<p>All models were implemented based on Python 3.9 and Scikit-learn 1.2.2, and inference time tests were performed on a standard clinical workstation (Intel Core i7-12700H CPU, 16GB RAM, no GPU acceleration). The number of parameters of the GradientBoosting model was 12,800, and the single inference time was 0.023&#x202F;&#x00B1;&#x202F;0.004&#x202F;s, which met the clinical real-time prediction requirements (&#x003C;0.1&#x202F;s/case).</p>
</sec>
<sec id="sec13">
<label>2.6</label>
<title>Model evaluation</title>
<p>Discriminative performance was evaluated using the area under the receiver operating characteristic (ROC) curve (AUC), with sensitivity, specificity, accuracy, F1 score, recall, and Youden&#x2019;s index serving as additional metrics.</p>
<p>Calibration capability: Evaluated using calibration curves, which compare predicted probabilities of prolonged healing with observed event frequencies across quantiles of predicted risk.</p>
<p>Clinical utility was evaluated through Decision Curve Analysis (DCA), which measures the net clinical benefit of models across various threshold probabilities. The x-axis denotes the patient risk threshold, indicating the probability above which intervention is warranted for prolonged healing, while the y-axis illustrates the net benefit, calculated as the advantage of accurately identifying prolonged healing minus the detriment of false positives. DCA helps clinicians choose the model with the highest practical value by comparing net benefits against baseline strategies (treating all patients or treating none).</p>
</sec>
<sec id="sec14">
<label>2.7</label>
<title>Model interpretation based on SHAP</title>
<p>The SHAP (Shapley Additive exPlanations) algorithm, which was implemented through the DALEX and fastshap packages in R, was employed to evaluate the influence of each variable on predicting prolonged healing. The key steps were as follows:</p>
<list list-type="simple">
<list-item>
<p>Calculating global mean absolute SHAP values to determine the overall importance of features (higher values indicate greater impact on predictions).</p>
</list-item>
<list-item>
<p>Create SHAP summary plots to illustrate the distribution of SHAP values for each variable, highlighting the impact of feature values on prediction outcomes.</p>
</list-item>
<list-item>
<p>Constructing individual SHAP force plots to explain the prediction logic for specific patients, illustrating how each variable pushes the prediction toward &#x201C;prolonged healing&#x201D; or &#x201C;normal healing&#x201D; relative to the baseline.</p>
</list-item>
</list>
</sec>
<sec id="sec15">
<label>2.8</label>
<title>Statistical analysis</title>
<p>The sample size calculation of the present study based on previous similar diabetic wound forecast model research, through the estimation of G &#x002A; Power 3.1 software. With <italic>&#x03B1;</italic>&#x202F;=&#x202F;0.05, <italic>&#x03B2;</italic>&#x202F;=&#x202F;0.08 (power&#x202F;=&#x202F;0.92) and <italic>f</italic><sup>2</sup>&#x202F;=&#x202F;0.15, the minimum sample size was determined to be 380. The study enrolled 465 patients, which exceeded the minimum sample size and had sufficient power to detect between-group differences in the primary outcome. <italic>Post-hoc</italic> power analysis confirmed that for the primary predictor (CRP_2nd), the study achieved a power of 0.95 to detect associations with prolonged healing, exceeding the minimum requirement of 0.80.</p>
<p>To address potential information leakage from CRP_3rd (measured at discharge, close to the 8-week healing outcome), a sensitivity analysis was conducted by constructing a modified GradientBoosting model excluding CRP_3rd, therapeutic_response_2, and therapeutic_response_all. Performance of the modified model was compared to the original model to verify independence from late-phase CRP data.</p>
<p>Statistical analyses utilized R software (v4.4.3) with the stats, caret, pROC, and ggplot2 packages. Quantitative data were presented as median (25th percentile, 75th percentile) [M (P25, P75)] following assessments for normality and variance homogeneity. Independent sample t-tests were used for between-group comparisons of normally distributed data, while the Mann&#x2013;Whitney U test was applied for non-normally distributed data. Categorical data were presented as frequencies and proportions. For comparisons, the chi-square test was employed. When the expected frequencies were small, Fisher&#x2019;s exact test was used instead. Ordinal data were analyzed by the rank sum test. A two-sided <italic>p</italic>-value of less than 0.05 was considered to be statistically significant. The discriminative performance of the constructed models was mainly evaluated by the AUC.</p>
</sec>
</sec>
<sec sec-type="results" id="sec16">
<label>3</label>
<title>Results</title>
<sec id="sec17">
<label>3.1</label>
<title>Characteristics of the study population</title>
<p>The study included 465 patients with diabetic wounds, with 325 patients admitted between 2021 and 2024 forming the training set, and 140 patients admitted in 2025 constituting the temporal validation set. In the training set, 281 cases were categorized into the normal healing group (healing within 8&#x202F;weeks) and 184 cases into the prolonged healing group (healing beyond 8&#x202F;weeks). <xref rid="SM1" ref-type="supplementary-material">Supplementary Table S1</xref> presents the baseline characteristics and laboratory indicators for both groups in the training set.</p>
<p>The median wound area of the enrolled patients was 6.5&#x202F;cm<sup>2</sup> (interquartile range: 5.8&#x2013;7.3&#x202F;cm<sup>2</sup>), and there was no statistically significant difference in the distribution of wound area between the normal healing group and the prolonged healing group (<italic>p</italic>&#x202F;=&#x202F;0.326), excluding the confounding effect of area difference on healing outcomes. Sensitivity analysis excluding CRP_3rd showed the modified GradientBoosting model retained strong performance in the validation set: accuracy&#x202F;=&#x202F;0.8929, sensitivity&#x202F;=&#x202F;0.8148, specificity&#x202F;=&#x202F;0.9643, AUC&#x202F;=&#x202F;0.921. This is slightly lower than the original model (accuracy&#x202F;=&#x202F;0.9357, AUC&#x202F;=&#x202F;0.958) but confirms no significant information leakage&#x2014;preoperative variables (especially CRP_2nd) sufficiently capture predictive signals for prolonged healing.</p>
<p>The two groups exhibited no significant differences in age and gender distribution (<italic>p</italic>&#x202F;&#x003E;&#x202F;0.05). However, there was a significant difference in smoking history, with a higher proportion of smokers in the normal healing group compared to the prolonged healing group (35.2% vs. 26.1%, <italic>p</italic>&#x202F;=&#x202F;0.049). To explore the difference in the proportion of smokers between the two groups, further multivariate logistic regression analysis was performed. After adjusting for potential confounding factors such as age, sex, albumin level, and CRP_2nd, smoking status was not statistically associated with prolonged healing (OR&#x202F;=&#x202F;0.87, 95%CI: 0.56&#x2013;1.35, <italic>p</italic>&#x202F;=&#x202F;0.542). This finding may be related to the following factors: (1) Limited sample size: only 30.5% (142/465) of the subjects were smokers, which may be insufficient to detect the negative effect of smoking on healing; (2) Confounding factors: In the normal healing group, smokers had higher albumin levels (median 37.2&#x202F;g/L vs. Smokers in prolonged healing group 32.1&#x202F;g/L, <italic>p</italic>&#x202F;=&#x202F;0.028), while albumin is a known protective factor, which may offset the harmful effect of smoking. (3) Selection bias: there may be unmeasured confounders in the retrospective study (such as years of smoking, daily smoking amount, smoking cessation status, etc.), and we did not collect these detailed information, which may affect the interpretation of the results. Drinking history showed no significant variation between the groups (<italic>p</italic>&#x202F;=&#x202F;0.406).</p>
<p>The normal healing group demonstrated significantly higher albumin (ALB) levels and albumin/globulin ratio (ALB/GLO) [ALB 38.0 (34.3, 41.0) vs. 32.3 (28.3, 35.8) g/L; ALB/GLO 1.33 (1.20, 1.60) vs. 1.00 (0.84, 1.20)] and lower globulin (GLO) levels compared to the prolonged healing group (all <italic>p</italic>&#x202F;&#x003C;&#x202F;0.001), suggesting that better nutritional status is associated with faster wound healing.</p>
<p>The prolonged healing group consistently exhibited significantly higher C-reactive protein (CRP) levels at admission, post-antibiotic pre-operatively, and discharge (all <italic>p</italic>&#x202F;&#x003C;&#x202F;0.001). For example, their median CRP level was 86.4 (57.5, 140) mg/L, compared to 5.72 (3.30, 14.1) mg/L in the normal group. Therapeutic response analyses indicated that the normal group exhibited significantly superior early (therapeutic_response_1) and overall (therapeutic_response_all) responses (p&#x202F;&#x003C;&#x202F;0.001), while late-stage response (therapeutic_response_2) showed no significant difference between groups (<italic>p</italic>&#x202F;=&#x202F;0.356).</p>
<p>Hematological indices revealed distinct inflammatory profiles, with the prolonged healing group exhibiting elevated neutrophil percentage (Neut%), neutrophil count (Neut#), white blood cell count (WBC), and erythrocyte sedimentation rate (ESR) (all <italic>p</italic>&#x202F;&#x003C;&#x202F;0.001), along with reduced lymphocyte percentage (Lymph%) and count (Lymph#) (<italic>p</italic>&#x202F;&#x003C;&#x202F;0.001). This suggests a phenotype characterized by persistent neutrophilic infiltration in the prolonged group.</p>
<p>The prolonged healing group exhibited significantly elevated levels of glucose (Glu), alkaline phosphatase (ALP), and D-dimer, alongside reduced levels of hemoglobin (Hb), hematocrit (Hct), red blood cell count (RBC), sodium (Na<sup>+</sup>), and chloride (Cl<sup>&#x2212;</sup>) (all <italic>p</italic>&#x202F;&#x003C;&#x202F;0.05). The study suggests that delayed wound healing in diabetic patients may be influenced by abnormal glucose metabolism, coagulation dysfunction, anemia, and electrolyte imbalance.</p>
</sec>
<sec id="sec18">
<label>3.2</label>
<title>LASSO regression for feature selection</title>
<p>LASSO regression was employed to identify variables associated with prolonged healing, with the optimal penalty parameter (<italic>&#x03BB;</italic>) determined through 10-fold cross-validation (<xref ref-type="fig" rid="fig2">Figure 2</xref>). The trajectory of variable coefficients as &#x03BB; increased (<xref ref-type="fig" rid="fig2">Figure 2A</xref>) illustrated the progressive elimination of less predictive variables, while the cross-validation curve (<xref ref-type="fig" rid="fig2">Figure 2B</xref>) identified &#x03BB;&#x202F;=&#x202F;0.012 as the optimal value&#x2014;corresponding to the minimum classification error (0.083) and balancing model complexity and predictive power. The results showed that 15 variables were finally included in the model, including CRP_2nd, CRP_3rd, CRP, therapeutic_response_all, ALB, MCHC, Mono#, RDW-SD, HDL-C, TG, therapeutic_response_2, MCH, RBC, Eos%, and ESR. These variables were used as candidate features for subsequent machine learning model construction.</p>
<fig position="float" id="fig2">
<label>Figure 2</label>
<caption>
<p>LASSO regression for variable selection. <bold>(A)</bold> Trajectory of variable coefficients as a function of the penalty parameter (<italic>&#x03BB;</italic>); <bold>(B)</bold> cross-validation curve with the optimal &#x03BB;&#x202F;=&#x202F;0.012 (dotted line) corresponding to the minimum classification error (0.083). This &#x03BB; value selected 15 key variables (including CRP_2nd, CRP_3rd, and ALB) for model construction, balancing predictive power and parsimony.</p>
</caption>
<graphic xlink:href="fmed-13-1778003-g002.tif" mimetype="image" mime-subtype="tiff">
<alt-text content-type="machine-generated">Graph A shows a line plot depicting coefficient paths against log lambda values, demonstrating changes in coefficients over a range. Graph B is a plot of binomial deviance against log lambda, with a curve indicating deviance values. Red and blue dotted lines mark specific lambda values, highlighting minimum lambda at 0.002 with 26 variables, and one standard error lambda at 0.018 with 10 variables.</alt-text>
</graphic>
</fig>
</sec>
<sec id="sec19">
<label>3.3</label>
<title>Performance comparison of machine learning models</title>
<p>Twelve machine learning models were developed to forecast extended healing times, with their performance on both the training and validation sets illustrated in <xref ref-type="fig" rid="fig3">Figure 3</xref> and <xref rid="SM1" ref-type="supplementary-material">Supplementary Tables S2, S3</xref>.</p>
<fig position="float" id="fig3">
<label>Figure 3</label>
<caption>
<p>Performance metrics <bold>(A,B)</bold> and ROC curves <bold>(C,D)</bold> of 12 machine learning models in training <bold>(A,C)</bold> and validation <bold>(B,D)</bold> sets. Metrics include accuracy, F1 score, negative predictive value (NPV), positive predictive value (PPV), sensitivity, and specificity. Clinical relevance: The GradientBoosting model (red line) exhibits balanced performance across training and validation sets, with high specificity (minimizing false positives) and strong sensitivity (identifying most high-risk patients), making it suitable for clinical risk stratification.</p>
</caption>
<graphic xlink:href="fmed-13-1778003-g003.tif" mimetype="image" mime-subtype="tiff">
<alt-text content-type="machine-generated">Panel A and B are line graphs showing training and validation set performance metrics for different models, including accuracy, F1, NPV, PPV, sensitivity, specificity, and Youden&#x2019;s Index. Each model is represented by a different colored line. Panel C and D display ROC curves for training and validation sets, respectively, with multiple models plotted. The legend indicates model AUC values, showcasing performance differences among models like RandomForest, GradientBoosting, and SVM_Kernel.</alt-text>
</graphic>
</fig>
<sec id="sec20">
<label>3.3.1</label>
<title>Training set performance</title>
<p>As shown in the training set data, different models demonstrated varying levels of performance across multiple metrics. The BoostingMethod and AdaptiveBoosting models demonstrated high performance, achieving an accuracy of 0.9846, sensitivity of 0.9675, and specificity of 0.9950. Their F1 scores also reached 0.9794, and Youden&#x2019;s index was 0.9625.</p>
<p>The RandomForest model demonstrated strong performance with an accuracy of 0.9785, sensitivity of 0.9512, specificity of 0.9950, and an F1 score of 0.9709. However, the GradientBoosting model presented a competitive balance of metrics. The model demonstrated similar sensitivity to RandomForest at 0.9512, achieving an accuracy of 0.9754 and a specificity of 0.9901. Its F1 score was 0.9669, and Youden&#x2019;s index was 0.9413.</p>
<p>On the other hand, models like NeighborMethod and PLSModel had relatively lower performance in the training set. The NeighborMethod achieved an accuracy of 0.7631, a sensitivity of 0.6179, and a specificity of 0.8515. The PLSModel demonstrated an accuracy of 0.8031, a sensitivity of 0.5691, and a high specificity of 0.9455.</p>
<p>The confusion matrix of the optimal GradientBoosting model in the training set further illustrates the distribution of true and predicted labels, confirming the model&#x2019;s robustness in distinguishing prolonged healing from normal healing cases (<xref rid="SM1" ref-type="supplementary-material">Supplementary Figure S1</xref>).</p>
</sec>
<sec id="sec21">
<label>3.3.2</label>
<title>Validation set performance</title>
<p>In the validation set, the BoostingMethod, NeuralNet, and AdaptiveBoosting models each demonstrated high performance, with accuracy, sensitivity, and specificity values of 0.9429, 0.8852, and 0.9873, respectively. Their F1 scores were 0.9310, and Youden&#x2019;s index reached 0.8726.</p>
<p>The GradientBoosting model also showed strong performance. It had an accuracy of 0.9357, sensitivity of 0.8689, and specificity of 0.9873, which were only slightly lower than the top&#x2014;performing trio. Its F1 score was 0.9217, and Youden&#x2019;s index was 0.8562.</p>
<p>The GradientBoosting model emerges as the optimal choice when evaluating both the training and validation datasets. Although the BoostingMethod and AdaptiveBoosting models slightly outperformed in the training set, their performance was comparable to the GradientBoosting model in the validation set. The GradientBoosting model demonstrated a good balance of performance across the two sets, indicating its stability and generalizability. This makes it more reliable for predicting prolonged healing compared to other models. The confusion matrix of the GradientBoosting model in the temporal validation set quantifies the true positive (prolonged healing correctly predicted), true negative (normal healing correctly predicted), false positive, and false negative cases, further verifying the model&#x2019;s clinical applicability (<xref rid="SM1" ref-type="supplementary-material">Supplementary Figure S2</xref>).</p>
<p>The optimal model was selected based on balanced performance across training/validation sets and statistical comparisons using 95% confidence intervals (CIs). For the validation set, GradientBoosting (accuracy&#x202F;=&#x202F;0.9357, 95% CI: 0.884&#x2013;0.968), BoostingMethod (0.9429, 95% CI: 0.893&#x2013;0.973), and AdaptiveBoosting (0.9429, 95% CI: 0.893&#x2013;0.973) showed no significant differences (overlapping CIs). GradientBoosting was prioritized due to: (1) Superior generalizability (smallest training-validation accuracy difference: 3.97% vs. 4.17% for BoostingMethod and 4.32% for AdaptiveBoosting); (2) Compatibility with SHAP explainability (vs. less transparent weight-based mechanisms in AdaptiveBoosting); (3) Higher computational efficiency (inference time: 0.023&#x202F;&#x00B1;&#x202F;0.004&#x202F;s vs. 0.057&#x202F;&#x00B1;&#x202F;0.008&#x202F;s for BoostingMethod); (4) Greater robustness to missing data (performance decline &#x003C; 5% with 5% random missingness vs. 8%&#x2013;10% for SVM_Kernel/PLSModel).</p>
</sec>
</sec>
<sec id="sec22">
<label>3.4</label>
<title>Decision curve analysis (DCA) results</title>
<p>Decision Curve Analysis (DCA) was performed to evaluate the net clinical benefit of all 12 machine learning models across a range of clinically relevant risk thresholds (<xref ref-type="fig" rid="fig4">Figure 4</xref>). As illustrated in <xref ref-type="fig" rid="fig4">Figure 4A</xref> (training set) and <xref ref-type="fig" rid="fig4">Figure 4B</xref> (validation set), the GradientBoosting model&#x2014;consistently identified as the top-performing model based on accuracy, sensitivity, and specificity&#x2014;exhibited the highest net clinical benefit across most threshold probabilities (0.1&#x2013;0.5). This indicates that the model can maximize the clinical utility of predicting prolonged diabetic wound healing by correctly identifying high-risk patients who would benefit from targeted interventions, while minimizing unnecessary treatments for low-risk individuals. In contrast, models with relatively inferior performance in prior analyses (e.g., NeighborMethod and PLSModel) showed limited net benefit in DCA: their curves either closely aligned with the &#x201C;treat none&#x201D; baseline (reflecting minimal added value) or crossed the curves of more effective models only at low threshold probabilities, making them less suitable for guiding optimal clinical decision-making.</p>
<fig position="float" id="fig4">
<label>Figure 4</label>
<caption>
<p>Decision curve analysis (DCA) of 12 machine learning models. <bold>(A)</bold> Training set; <bold>(B)</bold> Validation set. The <italic>x</italic>-axis represents the risk threshold (probability of prolonged healing above which intervention is warranted), and the <italic>y</italic>-axis represents net clinical benefit (benefit of correctly identifying high-risk patients minus harm of false positives). Clinical relevance: The GradientBoosting model (red line) provides the highest net clinical benefit across clinically meaningful thresholds (0.1&#x2013;0.5), outperforming other models and baseline strategies (treating all or no patients).</p>
</caption>
<graphic xlink:href="fmed-13-1778003-g004.tif" mimetype="image" mime-subtype="tiff">
<alt-text content-type="machine-generated">Decision curve analysis charts for training set (A) and validation set (B). Both display net benefit versus high-risk threshold for various predictive strategies. Each line represents different models, such as RandomForest, GradientBoosting, and LogisticModel, among others. The "Treat all" and "Treat none" strategies are also indicated.</alt-text>
</graphic>
</fig>
</sec>
<sec id="sec23">
<label>3.5</label>
<title>SHAP analysis for model interpretation</title>
<p>To dissect how the model arrives at predictions, SHAP (SHapley Additive exPlanations) analysis was employed, yielding insights across global feature influence, directional impacts, and individual case behavior.</p>
<sec id="sec24">
<label>3.5.1</label>
<title>Global feature importance</title>
<p>The SHAP summary plot (<xref ref-type="fig" rid="fig5">Figure 5A</xref>) ranks features by their mean absolute contribution to prediction variance. The top five predictors shaping model output were: CRP_2nd (mean absolute SHAP&#x202F;=&#x202F;0.460), CRP_3rd (0.243), CRP (0.080), therapeutic_response_all (0.034), and ALB (0.020). Notably, CRP_2nd emerged as the dominant driver&#x2014;its substantially higher SHAP value underscores its pivotal role in identifying prolonged healing risk. This aligns with clinical intuition, as serial inflammatory markers (e.g., second - stage CRP) better capture dynamic infection/repair processes than single&#x2014;timepoint measures.</p>
<fig position="float" id="fig5">
<label>Figure 5</label>
<caption>
<p>SHAP analysis of the GradientBoosting model. <bold>(A)</bold> Feature importance ranked by mean absolute SHAP values (higher values indicate greater predictive impact); <bold>(B)</bold> SHAP bee swarm plot (red&#x202F;=&#x202F;high feature values, blue&#x202F;=&#x202F;low feature values). Clinical relevance: CRP_2nd is the most influential predictor (mean absolute SHAP&#x202F;=&#x202F;0.460), with high values (red points) increasing the risk of prolonged healing (negative SHAP values) and high albumin (ALB) values (red points) reducing risk (positive SHAP values). CRP_2nd clinical thresholds (Low: &#x003C;30&#x202F;mg/L; Intermediate: 30&#x2013;80&#x202F;mg/L; High: &#x003E;80&#x202F;mg/L) derived from these data enable actionable risk stratification.</p>
</caption>
<graphic xlink:href="fmed-13-1778003-g005.tif" mimetype="image" mime-subtype="tiff">
<alt-text content-type="machine-generated">Chart A shows a bar graph of feature importance using mean absolute SHAP values, highlighting CRP_2nd, CRP_3rd, and CRP as most influential. Chart B depicts a SHAP value distribution bee swarm plot, where each point represents a sample. Features like CRP_2nd and CRP_3rd show varied impacts on the model output, indicated by the spread and color gradient from red (high feature values) to blue (low feature values).</alt-text>
</graphic>
</fig>
</sec>
<sec id="sec25">
<label>3.5.2</label>
<title>Directionality of feature impacts</title>
<p>The SHAP bee swarm plot (<xref ref-type="fig" rid="fig5">Figure 5B</xref>) clarifies how feature values tilt predictions. Elevated CRP_2nd, CRP_3rd, and CRP correlated with negative SHAP contributions, pushing the model toward a &#x201C;prolonged healing&#x201D; classification. Conversely, higher ALB (a nutritional marker) and favorable therapeutic_response_all (composite treatment efficacy) aligned with positive SHAP values, nudging predictions toward &#x201C;normal healing.&#x201D; This duality highlights a balance: inflammation (via CRP) drives risk, while nutrition/therapy response mitigates it. For example, strong therapeutic response signals effective intervention, reducing failure likelihood&#x2014;an actionable insight for clinical monitoring. Based on SHAP values and ROC curve analysis, CRP_2nd was stratified into three clinical thresholds: (1) Low risk: CRP_2nd&#x202F;&#x003C;&#x202F;30&#x202F;mg/L (mean SHAP value &#x003C; &#x2212;0.2, predicted probability of prolonged healing &#x003C; 0.3); (2) Intermediate risk: 30&#x202F;mg/L&#x202F;&#x2264;&#x202F;CRP_2nd&#x202F;&#x2264;&#x202F;80&#x202F;mg/L (mean SHAP value: &#x2212;0.2 to 0.3, predicted probability: 0.3&#x2013;0.7); (3) High risk: CRP_2nd&#x202F;&#x003E;&#x202F;80&#x202F;mg/L (mean SHAP value &#x003E; 0.3, predicted probability &#x003E; 0.7). These thresholds align with clinical definitions of inflammation severity and provide actionable risk stratification (<xref ref-type="fig" rid="fig5">Figure 5</xref>).</p>
</sec>
<sec id="sec26">
<label>3.5.3</label>
<title>Individual case validation</title>
<p>Two representative cases (<xref ref-type="fig" rid="fig6">Figure 6</xref>) illustrate model logic in practice:</p>
<fig position="float" id="fig6">
<label>Figure 6</label>
<caption>
<p>Individual prediction interpretation using SHAP force plots <bold>(A,B)</bold> and waterfall plots <bold>(C,D)</bold> for Sample 5 (prolonged healing) and Sample 6 (normal healing). Each feature&#x2019;s cumulative impact on the final prediction is illustrated. Clinical relevance: Sample 5 (prolonged healing) is driven by a high CRP_2nd (142&#x202F;mg/L, SHAP contribution +0.516), while Sample 6 (normal healing) benefits from a robust early therapeutic response (SHAP contribution &#x2212;0.398), demonstrating how the model identifies modifiable risk factors for clinical intervention.</p>
</caption>
<graphic xlink:href="fmed-13-1778003-g006.tif" mimetype="image" mime-subtype="tiff">
<alt-text content-type="machine-generated">Bar charts A and B illustrate feature contributions for Samples 5 and 6, showing visual forces affecting prediction. Chart A highlights a +0.516 impact; Chart B features a +0.277 and -0.398 impact. Charts C and D provide a breakdown, detailing cumulative contributions for each sample, with factors like CRP levels and therapeutic responses prominently influencing predictions.</alt-text>
</graphic>
</fig>
<p>Prolonged Healing (Sample 5): A spike in CRP_2nd (142&#x202F;mg/L) generated a large negative SHAP contribution (+0.516), overwhelming other features and driving the &#x201C;prolonged&#x201D; prediction. This mirrors real - world scenarios where unresolved inflammation stalls healing.</p>
<p>Normal Healing (Sample 6): A robust therapeutic_response_1 (61% improvement) produced a negative SHAP value (&#x2212;0.398), tipping the model toward &#x201C;normal healing.&#x201D; This showcases how early treatment success signals favorable outcomes, justifying proactive therapy escalation for poor responders.</p>
<p>By integrating global trends, directional effects, and individual examples, SHAP analysis demystifies the model: persistent inflammation (especially CRP_2nd) is the prime risk factor, while nutrition/therapy response buffers it. These insights not only validate the model but also guide clinical action&#x2014;e.g., prioritizing CRP monitoring and therapeutic optimization for at - risk patients. Supplementary SHAP value analyses further elaborate on the feature interpretation, including the partial dependence of prediction outcomes on key variables (e.g., CRP_2nd and ALB) and the interactive effects between top predictors, providing more detailed mechanistic insights into the model&#x2019;s decision-making process (<xref rid="SM1" ref-type="supplementary-material">Supplementary Figure S3</xref>).</p>
</sec>
</sec>
</sec>
<sec sec-type="discussion" id="sec27">
<label>4</label>
<title>Discussion</title>
<p>This study developed and validated 12 machine learning models to predict prolonged diabetic wound healing, with the GradientBoosting model demonstrating superior performance in both training and validation sets. Through SHAP (SHapley Additive exPlanations) analysis, we identified the core factors driving model predictions: dynamic inflammatory markers&#x2014;specifically CRP_2nd (C-reactive protein measured after admission antibiotic therapy but before surgical debridement), CRP_3rd (CRP at discharge following debridement), and baseline CRP (CRP at admission)&#x2014;emerged as the most influential predictors. Among these, CRP_2nd (mean absolute SHAP value 0.460) ranked first. This finding aligns with the well-established role of inflammation as a central regulator of wound repair, where both insufficient and excessive inflammation can disrupt the healing cascade (<xref ref-type="bibr" rid="ref16">16</xref>). Physiological mechanisms for CRP_2nd to be more predictive than other time points: (1) specificity of time window: CRP_2nd was measured after antibiotic treatment and before surgery (3&#x2013;7&#x202F;days after admission). This time point is a key window to evaluate the effect of anti-infective treatment and the ability to clear inflammation - if the antibiotic treatment is effective and the pathogen is controlled, CRP should decrease significantly. If CRP persistently increased (CRP_2nd&#x202F;&#x003E;&#x202F;50&#x202F;mg/L), it indicated that the infection was not controlled or tissue repair was impaired, and it was an important warning signal for subsequent prolonged healing. CRP_2nd directly reflects the balance between inflammation and repair. CRP on admission may be disturbed by acute stress and complications, and CRP_3rd after surgery may be affected by surgical traumatic stress, with low predictive specificity. (3) The time window corresponding to CRP_2nd provides an opportunity for preoperative intervention. Timely adjustment of treatment (such as replacement of antibiotics and albumin supplementation) for patients with elevated CRP_2nd can effectively improve the healing outcome, while CRP measurement at other time points lack such a critical intervention window.</p>
<p>Our results demonstrated that the prolonged healing group displayed significantly elevated CRP levels at all three time points compared to the normal healing group (<xref rid="SM1" ref-type="supplementary-material">Supplementary Table S1</xref>). The prolonged group exhibited a significantly slower rate of CRP reduction (<italic>p</italic>&#x202F;&#x003C;&#x202F;0.001), emphasizing that &#x201C;inflammatory resolution&#x201D; rather than the initial inflammation level influences healing outcomes. This aligns with research indicating that prolonged inflammation in diabetic wounds hinders the shift from the inflammatory to the proliferative phase, marked by decreased fibroblast activity and angiogenesis (<xref ref-type="bibr" rid="ref17 ref18 ref19 ref20 ref21 ref22">17&#x2013;22</xref>). For instance, Eming et al. (<xref ref-type="bibr" rid="ref23">23</xref>) noted that persistent inflammation in chronic wounds leads to tissue damage by promoting excessive release of reactive oxygen species (ROS) and matrix metalloproteinases (MMPs), which directly hinder collagen synthesis.</p>
<p>Consistent with this, recent studies have indicated that non-healing wounds possess the characteristic of &#x201C;an unstable spatiotemporal balance between inflammation and protease activity&#x201D;: Uncontrolled inflammation leads to excessive secretion of proteases, which in turn damages growth factors and matrix proteins, thereby hindering the healing process (<xref ref-type="bibr" rid="ref24">24</xref>).</p>
<p>The AGE-RAGE axis, a central pathway in diabetic complications, offers a mechanistic framework underlying our findings (<xref ref-type="bibr" rid="ref25 ref26 ref27">25&#x2013;27</xref>). The activation of RAGE by AGEs or EN-RAGEs induces the release of proinflammatory cytokines such as TNF-&#x03B1; and IL-1&#x03B2;, which perpetuates inflammation and protease activity (<xref ref-type="bibr" rid="ref27 ref28 ref29 ref30">27&#x2013;30</xref>). While our study did not directly measure AGEs or RAGE, the prominence of CRP&#x2014;a downstream marker of this axis&#x2014;in prediction suggests that our model implicitly encapsulates this pathological cascade. This reinforces the clinical relevance of our machine learning approach: by prioritizing serial CRP monitoring, clinicians can indirectly assess RAGE-mediated inflammation and intervene to restore the &#x201C;self-limited&#x201D; inflammatory phase necessary for healing.</p>
<p>Beyond CRP dynamics, hematological indices revealed distinct inflammatory phenotypes between groups. The prolonged healing group exhibited significantly higher neutrophil percentage (Neut%) and count (Neut#), alongside lower lymphocyte percentage (Lymph%) and count (Lymph#) (<xref rid="SM1" ref-type="supplementary-material">Supplementary Table S1</xref>). This pattern reflects a &#x201C;persistent neutrophilic inflammation&#x201D; phenotype, which is pathological in the context of chronic wounds.</p>
<p>Neutrophils are critical for clearing pathogens in acute wounds, but their persistent infiltration in diabetic wounds leads to excessive degranulation and NETosis (neutrophil extracellular trap formation), inducing damage to surrounding healthy tissue and impairing re-epithelialization (<xref ref-type="bibr" rid="ref31">31</xref>). Our findings align with those of Rui et al. (<xref ref-type="bibr" rid="ref32">32</xref>), who reported that neutrophil hyperactivity in diabetic ulcers correlates with elevated MMP-8 levels, which degrade growth factors and extracellular matrix. Conversely, reduced lymphocytes in the prolonged group may indicate impaired adaptive immunity&#x2014;lymphocytes are essential for transitioning to the proliferative phase by secreting cytokines like TGF-&#x03B2;1, which promotes fibroblast recruitment (<xref ref-type="bibr" rid="ref33">33</xref>). This imbalance between neutrophils and lymphocytes underscores the need for targeted anti-inflammatory strategies (e.g., neutrophil apoptosis promoters or lymphocyte modulators) in high-risk patients.</p>
<p>Higher levels of albumin (ALB) and positive therapeutic responses (therapeutic_response_all) were identified as protective factors, both associated with a decreased risk of prolonged healing. Our results further demonstrated that inflammation does not operate in isolation but interacts with nutritional and metabolic status to influence healing. The prolonged group exhibited notably lower albumin (ALB) levels (<italic>p</italic>&#x202F;&#x003C;&#x202F;0.001), with SHAP analysis highlighting ALB as a crucial predictive feature (mean absolute SHAP value&#x202F;=&#x202F;0.020). Albumin, a major plasma protein, modulates inflammation by binding proinflammatory mediators (e.g., LPS) and supporting immune cell function; hypoalbuminemia thus exacerbates inflammation and impairs tissue repair (<xref ref-type="bibr" rid="ref34">34</xref>). A study supports this by demonstrating that adding albumin to the diet of malnourished diabetic patients speeds up wound healing by lowering CRP levels and boosting neutrophil phagocytosis.</p>
<p>The prolonged group exhibited increased glucose levels (<italic>p</italic>&#x202F;=&#x202F;0.021), aligning with findings that hyperglycemia enhances oxidative stress and advanced glycation end product accumulation, thereby intensifying inflammation through NF-&#x03BA;B pathway activation (<xref ref-type="bibr" rid="ref35">35</xref>). Notably, our data showed that short-term glucose fluctuations (Glu) exerted a more pronounced effect than long-term glycemic control (HbA1c, <italic>p</italic>&#x202F;=&#x202F;0.092), emphasizing the need for tight perioperative glucose management to mitigate inflammation.</p>
<p>The cumulative evidence from our analyses&#x2014;highlighting dynamic inflammation (particularly CRP_2nd) as a central driver, alongside its interplay with neutrophilic phenotypes, nutritional status (ALB), and metabolic factors (perioperative glucose)&#x2014;translates into actionable clinical insights for optimizing diabetic wound management:</p>
<p>First, in terms of early risk stratification, the prominence of CRP_2nd&#x2014;measured after initial antibiotic therapy but before debridement&#x2014;marks a critical window for identifying high-risk patients. This aligns with our key finding that inflammatory resolution, rather than initial inflammation magnitude, dictates outcomes. By monitoring CRP_2nd alongside early therapeutic response (therapeutic_response_1), clinicians can stratify patients pre-operatively, enabling timely adjustments to antibiotic regimens or nutritional support (e.g., albumin supplementation) to mitigate persistent inflammation before it entrenches.</p>
<p>Second, the distinct neutrophilic inflammation phenotype (elevated Neut#/Neut% with reduced Lymph#/Lymph%) in prolonged healers points to the need for targeted anti-inflammatory strategies. This phenotype, linked to excessive NETosis and MMP-mediated tissue damage in our mechanistic framework, supports interventions aimed at dampening neutrophil hyperactivity&#x2014;such as topical anti-TNF-&#x03B1; therapy to curb proinflammatory cytokine release or neutrophil elastase inhibitors to limit matrix degradation&#x2014;thereby restoring the inflammatory balance required for progression to proliferation.</p>
<p>Finally, these findings underscore the value of integrated management strategies that address the interplay between inflammation, nutrition, and metabolism. Given our observations that hypoalbuminemia exacerbates inflammation and perioperative hyperglycemia amplifies NF-&#x03BA;B-mediated inflammatory responses, concurrent efforts to optimize glycemic control, supplement albumin in malnourished patients, and modulate inflammation through targeted therapies could synergistically improve the wound microenvironment, ultimately accelerating healing (<xref ref-type="bibr" rid="ref36">36</xref>, <xref ref-type="bibr" rid="ref37">37</xref>).</p>
<p>The clinical transformation path of the model is as follows: (1) Integration into the electronic medical record (EMR) system: The model can be embedded into the hospital&#x2019;s existing EMR system via an API interface to automatically extract patients&#x2019;CRP test results, biochemical indicators, and other relevant data, generate real-time prolonged healing risk scores (0&#x2013;100 points), and assist clinicians in formulating personalized treatment plans; (2) Outpatient risk stratification tool: For outpatients with diabetic wounds, rapid risk assessment can be performed using pre-debridement CRP_2nd test results, with high-risk patients (score &#x2265;70 points or CRP_2nd&#x202F;&#x003E;&#x202F;80&#x202F;mg/L) prioritized for intensive interventions such as albumin supplementation (target &#x2265;35&#x202F;g/L) and targeted anti-inflammatory therapy, intermediate-risk patients (score 40&#x2013;69 or 30&#x202F;mg/L&#x202F;&#x2264;&#x202F;CRP_2nd&#x202F;&#x2264;&#x202F;80&#x202F;mg/L) receiving intensified wound care and close CRP monitoring, and low-risk patients (score &#x003C; 40 or CRP_2nd&#x202F;&#x003C;&#x202F;30&#x202F;mg/L) continuing standard care with discharge CRP_3rd measurement; (3) Remote follow-up monitoring: Combined with dynamic CRP monitoring data, the model updates risk assessments in real time to guide home care practices and optimize follow-up visit timing (weekly for high-risk, biweekly for intermediate-risk patients), with the risk threshold (score &#x2265;70, predicted probability &#x2265;0.75) derived from the validation set ROC curve&#x2019;s Youden index (J&#x202F;=&#x202F;0.856) to balance sensitivity and specificity.</p>
<p>Compared with the existing studies, the core innovation of this study is as follows: (1) For the first time, CRP_2nd was identified as the optimal predictor by focusing on the &#x201C;dynamic CRP trajectory&#x201D; instead of the static baseline value; (2) Standardized wound size (5&#x2013;8&#x202F;cm<sup>2</sup>) and unified intervention method (debridement only) were used to reduce heterogeneous interference; (3) The generalization of the model was verified by an external independent cohort (140 cases in 2025), while most of the previous studies lacked external validation or had high cohort heterogeneity. (4) Combining SHAP analysis to reveal the predictive mechanism and improve the clinical interpretability of the model.</p>
<p>This study has several limitations: (1) Single-center retrospective design: Training and temporal validation cohorts were from the same institution, limiting generalizability to diverse populations and care protocols. Future multi-center validation (planned n&#x202F;&#x2265;&#x202F;500) is needed to confirm broad applicability. (2) Follow-up bias: Healing status was confirmed via in-person (68%) or telephone (32%) follow-up. Telephone follow-up may introduce misclassification bias (e.g., overreporting of healing), mitigated by standardized follow-up questions and cross-referencing with outpatient records. (3) Limited biomarker sampling: CRP was measured at only three time points; more frequent sampling (e.g., postoperative days 1&#x2013;3) could capture finer inflammatory fluctuations. (4) Missing variables: Microbiological and imaging data were excluded due to high missing rates (&#x003E;30%), which may improve model performance if integrated in future studies.</p>
<p>Future research directions: (1) expand the multi-center cohort validation and include patients from different regions and hospitals of different levels to improve the representation of the model population; (2) integrating imaging indicators (such as wound ultrasound and infrared imaging) and microbiome data to further optimize the performance of the model; (3) develop mobile applications for convenience of primary medical institutions and patients; (4) Prospective clinical trials should be conducted to verify whether model-guided treatment strategies can improve the prognosis of patients.</p>
<p>Taken together, our findings underscore the central role of inflammatory dynamics&#x2014;particularly the post-antibiotic, pre-operative CRP trajectory&#x2014;in governing diabetic wound healing outcomes. By unraveling the interplay between inflammation, nutrition, and metabolism through explainable machine learning, this work not only enhances our understanding of diabetic wound pathophysiology but also provides a foundation for developing targeted, inflammation-modulating strategies to improve clinical outcomes.</p>
</sec>
<sec sec-type="conclusions" id="sec28">
<label>5</label>
<title>Conclusion</title>
<p>This study developed a GradientBoosting model to predict prolonged healing of diabetic wounds, with SHAP analysis revealing that dynamic inflammatory markers&#x2014;especially CRP_2nd (post-antibiotic, pre-debridement)&#x2014;are the most critical predictors. Inflammatory resolution, rather than initial inflammation magnitude, may determine healing outcomes, with persistent neutrophilic inflammation, hypoalbuminemia, and perioperative hyperglycemia potentially exacerbating delays. This work highlights the post-antibiotic, pre-operative phase as a critical window for risk stratification and provides a rationale for targeted interventions (e.g., anti-inflammatory therapies, albumin supplementation) to restore inflammatory balance. Future multi-center studies with more granular inflammatory monitoring may help translate these insights into clinical practice.</p>
</sec>
</body>
<back>
<sec sec-type="data-availability" id="sec29">
<title>Data availability statement</title>
<p>The original contributions presented in the study are included in the article/<xref rid="SM1" ref-type="supplementary-material">Supplementary material</xref>, further inquiries can be directed to the corresponding authors.</p>
</sec>
<sec sec-type="ethics-statement" id="sec30">
<title>Ethics statement</title>
<p>The studies involving humans were approved by The Ethics Committee of the First Affiliated Hospital of Dalian Medical University approved the protocol (Approval No: PJ-KS-KY-2025-440). All procedures involving human participants were performed in compliance with relevant ethical guidelines and regulations. 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.</p>
</sec>
<sec sec-type="author-contributions" id="sec31">
<title>Author contributions</title>
<p>SJ: Data curation, Writing &#x2013; original draft, Conceptualization, Formal analysis. JW: Formal analysis, Writing &#x2013; original draft, Data curation, Conceptualization. QS: Conceptualization, Data curation, Formal analysis, Writing &#x2013; original draft. XS: Investigation, Software, Methodology, Writing &#x2013; review &#x0026; editing. ZZ: Visualization, Resources, Writing &#x2013; review &#x0026; editing, Validation. JZ: Writing &#x2013; review &#x0026; editing, Supervision, Software. SW: Conceptualization, Writing &#x2013; review &#x0026; editing, Funding acquisition.</p>
</sec>
<ack>
<title>Acknowledgments</title>
<p>We would like to express our gratitude to the developers and maintainers of the R packages used in this study, including glmnet for LASSO regression, caret for machine learning model training and validation, pROC for ROC curve analysis, shapviz and fastshap for SHAP value calculation and visualization, and ggplot2 for data visualization. Their dedicated work has significantly facilitated the implementation of statistical analyses and model development in this research. We express gratitude to the clinical staff at the First Affiliated Hospital of Dalian Medical University for their assistance in data collection and patient follow-up.</p>
</ack>
<sec sec-type="COI-statement" id="sec32">
<title>Conflict of interest</title>
<p>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.</p>
</sec>
<sec sec-type="ai-statement" id="sec33">
<title>Generative AI statement</title>
<p>The author(s) declared that Generative AI was not used in the creation of this manuscript.</p>
<p>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.</p>
</sec>
<sec sec-type="disclaimer" id="sec34">
<title>Publisher&#x2019;s note</title>
<p>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.</p>
</sec>
<sec sec-type="supplementary-material" id="sec35">
<title>Supplementary material</title>
<p>The Supplementary material for this article can be found online at: <ext-link xlink:href="https://www.frontiersin.org/articles/10.3389/fmed.2026.1778003/full#supplementary-material" ext-link-type="uri">https://www.frontiersin.org/articles/10.3389/fmed.2026.1778003/full#supplementary-material</ext-link></p>
<supplementary-material xlink:href="Supplementary_file_1.docx" id="SM1" mimetype="application/vnd.openxmlformats-officedocument.wordprocessingml.document" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</sec>
<ref-list>
<title>References</title>
<ref id="ref1"><label>1.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zhou</surname><given-names>HH</given-names></name> <name><surname>Tang</surname><given-names>YL</given-names></name> <name><surname>Xu</surname><given-names>TH</given-names></name> <name><surname>Cheng</surname><given-names>B</given-names></name></person-group>. <article-title>C-reactive protein: structure, function, regulation, and role in clinical diseases</article-title>. <source>Front Immunol</source>. (<year>2024</year>) <volume>15</volume>:<fpage>1425168</fpage>. doi: <pub-id pub-id-type="doi">10.3389/fimmu.2024.1425168</pub-id>, <pub-id pub-id-type="pmid">38947332</pub-id></mixed-citation></ref>
<ref id="ref2"><label>2.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Salvo</surname><given-names>P</given-names></name> <name><surname>Dini</surname><given-names>V</given-names></name> <name><surname>Kirchhain</surname><given-names>A</given-names></name> <name><surname>Janowska</surname><given-names>A</given-names></name> <name><surname>Oranges</surname><given-names>T</given-names></name> <name><surname>Chiricozzi</surname><given-names>A</given-names></name> <etal/></person-group>. <article-title>Sensors and biosensors for C-reactive protein, temperature and pH, and their applications for monitoring wound healing: a review</article-title>. <source>Sensors (Basel)</source>. (<year>2017</year>) <volume>17</volume>:<fpage>2952</fpage>. doi: <pub-id pub-id-type="doi">10.3390/s17122952</pub-id>, <pub-id pub-id-type="pmid">29257113</pub-id></mixed-citation></ref>
<ref id="ref3"><label>3.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Reardon</surname><given-names>R</given-names></name> <name><surname>Simring</surname><given-names>D</given-names></name> <name><surname>Kim</surname><given-names>B</given-names></name> <name><surname>Mortensen</surname><given-names>J</given-names></name> <name><surname>Williams</surname><given-names>D</given-names></name> <name><surname>Leslie</surname><given-names>A</given-names></name></person-group>. <article-title>The diabetic foot ulcer</article-title>. <source>Aust J Gen Pract</source>. (<year>2020</year>) <volume>49</volume>:<fpage>250</fpage>&#x2013;<lpage>5</lpage>. doi: <pub-id pub-id-type="doi">10.31128/ajgp-11-19-5161</pub-id></mixed-citation></ref>
<ref id="ref4"><label>4.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Thorud</surname><given-names>JC</given-names></name> <name><surname>Plemmons</surname><given-names>B</given-names></name> <name><surname>Buckley</surname><given-names>CJ</given-names></name> <name><surname>Shibuya</surname><given-names>N</given-names></name> <name><surname>Jupiter</surname><given-names>DC</given-names></name></person-group>. <article-title>Mortality after nontraumatic major amputation among patients with diabetes and peripheral vascular disease: a systematic review</article-title>. <source>J Foot Ankle Surg</source>. (<year>2016</year>) <volume>55</volume>:<fpage>591</fpage>&#x2013;<lpage>9</lpage>. doi: <pub-id pub-id-type="doi">10.1053/j.jfas.2016.01.012</pub-id>, <pub-id pub-id-type="pmid">26898398</pub-id></mixed-citation></ref>
<ref id="ref5"><label>5.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Altindas</surname><given-names>M</given-names></name> <name><surname>Kilic</surname><given-names>A</given-names></name> <name><surname>Cinar</surname><given-names>C</given-names></name> <name><surname>Bingol</surname><given-names>UA</given-names></name> <name><surname>Ozturk</surname><given-names>G</given-names></name></person-group>. <article-title>The epidemiology of foot wounds in patients with diabetes: a description of 600 consecutive patients in Turkey</article-title>. <source>J Foot Ankle Surg</source>. (<year>2011</year>) <volume>50</volume>:<fpage>146</fpage>&#x2013;<lpage>52</lpage>. doi: <pub-id pub-id-type="doi">10.1053/j.jfas.2010.12.017</pub-id>, <pub-id pub-id-type="pmid">21353997</pub-id></mixed-citation></ref>
<ref id="ref6"><label>6.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Pierce</surname><given-names>GF</given-names></name></person-group>. <article-title>Inflammation in nonhealing diabetic wounds: the space-time continuum does matter</article-title>. <source>Am J Pathol</source>. (<year>2001</year>) <volume>159</volume>:<fpage>399</fpage>&#x2013;<lpage>403</lpage>. doi: <pub-id pub-id-type="doi">10.1016/s0002-9440(10)61709-9</pub-id></mixed-citation></ref>
<ref id="ref7"><label>7.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Liao</surname><given-names>C</given-names></name> <name><surname>Zhu</surname><given-names>M</given-names></name> <name><surname>Ding</surname><given-names>H</given-names></name> <name><surname>Li</surname><given-names>Y</given-names></name> <name><surname>Sun</surname><given-names>Q</given-names></name> <name><surname>Li</surname><given-names>X</given-names></name></person-group>. <article-title>Comparing the traditional and emerging therapies for enhancing wound healing in diabetic patients: a pivotal examination</article-title>. <source>Int Wound J</source>. (<year>2024</year>) <volume>21</volume>:<fpage>e14488</fpage>. doi: <pub-id pub-id-type="doi">10.1111/iwj.14488</pub-id>, <pub-id pub-id-type="pmid">37984812</pub-id></mixed-citation></ref>
<ref id="ref8"><label>8.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname><given-names>Y</given-names></name> <name><surname>Jiang</surname><given-names>C</given-names></name> <name><surname>Xing</surname><given-names>YQ</given-names></name> <name><surname>Zou</surname><given-names>L</given-names></name> <name><surname>Song</surname><given-names>M</given-names></name> <name><surname>Qu</surname><given-names>X</given-names></name> <etal/></person-group>. <article-title>Establishing a clinical prediction model for diabetic foot ulcers in type 2 diabetic patients with lower extremity arteriosclerotic occlusion using machine learning</article-title>. <source>Sci Rep</source>. (<year>2025</year>) <volume>15</volume>:<fpage>11737</fpage>. doi: <pub-id pub-id-type="doi">10.1038/s41598-025-96502-9</pub-id>, <pub-id pub-id-type="pmid">40188298</pub-id></mixed-citation></ref>
<ref id="ref9"><label>9.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Pathak</surname><given-names>JL</given-names></name> <name><surname>Bakker</surname><given-names>AD</given-names></name> <name><surname>Luyten</surname><given-names>FP</given-names></name> <name><surname>Verschueren</surname><given-names>P</given-names></name> <name><surname>Lems</surname><given-names>WF</given-names></name> <name><surname>Klein-Nulend</surname><given-names>J</given-names></name> <etal/></person-group>. <article-title>Systemic inflammation affects human osteocyte-specific protein and cytokine expression</article-title>. <source>Calcif Tissue Int</source>. (<year>2016</year>) <volume>98</volume>:<fpage>596</fpage>&#x2013;<lpage>608</lpage>. doi: <pub-id pub-id-type="doi">10.1007/s00223-016-0116-8</pub-id>, <pub-id pub-id-type="pmid">26887974</pub-id></mixed-citation></ref>
<ref id="ref10"><label>10.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Pierrakos</surname><given-names>C</given-names></name> <name><surname>Velissaris</surname><given-names>D</given-names></name> <name><surname>Bisdorff</surname><given-names>M</given-names></name> <name><surname>Marshall</surname><given-names>JC</given-names></name> <name><surname>Vincent</surname><given-names>JL</given-names></name></person-group>. <article-title>Biomarkers of sepsis: time for a reappraisal</article-title>. <source>Crit Care</source>. (<year>2020</year>) <volume>24</volume>:<fpage>287</fpage>. doi: <pub-id pub-id-type="doi">10.1186/s13054-020-02993-5</pub-id>, <pub-id pub-id-type="pmid">32503670</pub-id></mixed-citation></ref>
<ref id="ref11"><label>11.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Abdelhafiz</surname><given-names>D</given-names></name> <name><surname>Baker</surname><given-names>T</given-names></name> <name><surname>Glascow</surname><given-names>DA</given-names></name> <name><surname>Abdelhafiz</surname><given-names>A</given-names></name></person-group>. <article-title>Biomarkers for the diagnosis and treatment of rheumatoid arthritis - a systematic review</article-title>. <source>Postgrad Med</source>. (<year>2023</year>) <volume>135</volume>:<fpage>214</fpage>&#x2013;<lpage>23</lpage>. doi: <pub-id pub-id-type="doi">10.1080/00325481.2022.2052626</pub-id></mixed-citation></ref>
<ref id="ref12"><label>12.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname><given-names>YT</given-names></name> <name><surname>Zhang</surname><given-names>LX</given-names></name> <name><surname>Li</surname><given-names>Y</given-names></name> <name><surname>Zhao</surname><given-names>J</given-names></name> <name><surname>Chen</surname><given-names>HL</given-names></name></person-group>. <article-title>Inflammatory markers in diabetic foot infection: a meta-analysis</article-title>. <source>Wounds</source>. (<year>2023</year>) <volume>35</volume>:<fpage>425</fpage>&#x2013;<lpage>32</lpage>. doi: <pub-id pub-id-type="doi">10.25270/wnds/082421.03</pub-id>, <pub-id pub-id-type="pmid">38277631</pub-id></mixed-citation></ref>
<ref id="ref13"><label>13.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Stefanopoulos</surname><given-names>S</given-names></name> <name><surname>Ayoub</surname><given-names>S</given-names></name> <name><surname>Qiu</surname><given-names>Q</given-names></name> <name><surname>Ren</surname><given-names>G</given-names></name> <name><surname>Osman</surname><given-names>M</given-names></name> <name><surname>Nazzal</surname><given-names>M</given-names></name> <etal/></person-group>. <article-title>Machine learning prediction of diabetic foot ulcers in the inpatient population</article-title>. <source>Vascular</source>. (<year>2022</year>) <volume>30</volume>:<fpage>1115</fpage>&#x2013;<lpage>23</lpage>. doi: <pub-id pub-id-type="doi">10.1177/17085381211040984</pub-id>, <pub-id pub-id-type="pmid">34461765</pub-id></mixed-citation></ref>
<ref id="ref14"><label>14.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Yu</surname><given-names>X</given-names></name> <name><surname>Wu</surname><given-names>Z</given-names></name> <name><surname>Zhang</surname><given-names>N</given-names></name></person-group>. <article-title>Machine learning-driven discovery of novel therapeutic targets in diabetic foot ulcers</article-title>. <source>Mol Med</source>. (<year>2024</year>) <volume>30</volume>:<fpage>215</fpage>. doi: <pub-id pub-id-type="doi">10.1186/s10020-024-00955-z</pub-id>, <pub-id pub-id-type="pmid">39543487</pub-id></mixed-citation></ref>
<ref id="ref15"><label>15.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wu</surname><given-names>PY</given-names></name> <name><surname>Yu</surname><given-names>YL</given-names></name> <name><surname>Zhao</surname><given-names>WR</given-names></name> <name><surname>Zhou</surname><given-names>B</given-names></name></person-group>. <article-title>Identifying and validating extracellular matrix-related gene CTSH in diabetic foot ulcer using bioinformatics and machine learning</article-title>. <source>J Inflamm Res</source>. (<year>2024</year>) <volume>17</volume>:<fpage>5871</fpage>&#x2013;<lpage>87</lpage>. doi: <pub-id pub-id-type="doi">10.2147/jir.S467507</pub-id>, <pub-id pub-id-type="pmid">39228680</pub-id></mixed-citation></ref>
<ref id="ref16"><label>16.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Li</surname><given-names>M</given-names></name> <name><surname>Hou</surname><given-names>Q</given-names></name> <name><surname>Zhong</surname><given-names>L</given-names></name> <name><surname>Zhao</surname><given-names>Y</given-names></name> <name><surname>Fu</surname><given-names>X</given-names></name></person-group>. <article-title>Macrophage related chronic inflammation in non-healing wounds</article-title>. <source>Front Immunol</source>. (<year>2021</year>) <volume>12</volume>:<fpage>681710</fpage>. doi: <pub-id pub-id-type="doi">10.3389/fimmu.2021.681710</pub-id>, <pub-id pub-id-type="pmid">34220830</pub-id></mixed-citation></ref>
<ref id="ref17"><label>17.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Theocharidis</surname><given-names>G</given-names></name> <name><surname>Thomas</surname><given-names>BE</given-names></name> <name><surname>Sarkar</surname><given-names>D</given-names></name> <name><surname>Mumme</surname><given-names>HL</given-names></name> <name><surname>Pilcher</surname><given-names>WJR</given-names></name> <name><surname>Dwivedi</surname><given-names>B</given-names></name> <etal/></person-group>. <article-title>Single cell transcriptomic landscape of diabetic foot ulcers</article-title>. <source>Nat Commun</source>. (<year>2022</year>) <volume>13</volume>:<fpage>181</fpage>. doi: <pub-id pub-id-type="doi">10.1038/s41467-021-27801-8</pub-id>, <pub-id pub-id-type="pmid">35013299</pub-id></mixed-citation></ref>
<ref id="ref18"><label>18.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Li</surname><given-names>S</given-names></name> <name><surname>Ding</surname><given-names>X</given-names></name> <name><surname>Zhang</surname><given-names>H</given-names></name> <name><surname>Ding</surname><given-names>Y</given-names></name> <name><surname>Tan</surname><given-names>Q</given-names></name></person-group>. <article-title>IL-25 improves diabetic wound healing through stimulating M2 macrophage polarization and fibroblast activation</article-title>. <source>Int Immunopharmacol</source>. (<year>2022</year>) <volume>106</volume>:<fpage>108605</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.intimp.2022.108605</pub-id></mixed-citation></ref>
<ref id="ref19"><label>19.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Huang</surname><given-names>X</given-names></name> <name><surname>Liang</surname><given-names>P</given-names></name> <name><surname>Jiang</surname><given-names>B</given-names></name> <name><surname>Zhang</surname><given-names>P</given-names></name> <name><surname>Yu</surname><given-names>W</given-names></name> <name><surname>Duan</surname><given-names>M</given-names></name> <etal/></person-group>. <article-title>Hyperbaric oxygen potentiates diabetic wound healing by promoting fibroblast cell proliferation and endothelial cell angiogenesis</article-title>. <source>Life Sci</source>. (<year>2020</year>) <volume>259</volume>:<fpage>118246</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.lfs.2020.118246</pub-id>, <pub-id pub-id-type="pmid">32791151</pub-id></mixed-citation></ref>
<ref id="ref20"><label>20.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Voza</surname><given-names>FA</given-names></name> <name><surname>Huerta</surname><given-names>CT</given-names></name> <name><surname>Le</surname><given-names>N</given-names></name> <name><surname>Shao</surname><given-names>H</given-names></name> <name><surname>Ribieras</surname><given-names>A</given-names></name> <name><surname>Ortiz</surname><given-names>Y</given-names></name> <etal/></person-group>. <article-title>Fibroblasts in diabetic foot ulcers</article-title>. <source>Int J Mol Sci</source>. (<year>2024</year>) <volume>25</volume>:<fpage>2172</fpage>. doi: <pub-id pub-id-type="doi">10.3390/ijms25042172</pub-id>, <pub-id pub-id-type="pmid">38396848</pub-id></mixed-citation></ref>
<ref id="ref21"><label>21.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Guan</surname><given-names>Y</given-names></name> <name><surname>Niu</surname><given-names>H</given-names></name> <name><surname>Liu</surname><given-names>Z</given-names></name> <name><surname>Dang</surname><given-names>Y</given-names></name> <name><surname>Shen</surname><given-names>J</given-names></name> <name><surname>Zayed</surname><given-names>M</given-names></name> <etal/></person-group>. <article-title>Sustained oxygenation accelerates diabetic wound healing by promoting epithelialization and angiogenesis and decreasing inflammation</article-title>. <source>Sci Adv</source>. (<year>2021</year>) <volume>7</volume>:<fpage>eabj0153</fpage>. doi: <pub-id pub-id-type="doi">10.1126/sciadv.abj0153</pub-id>, <pub-id pub-id-type="pmid">34452918</pub-id></mixed-citation></ref>
<ref id="ref22"><label>22.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Elliott</surname><given-names>CG</given-names></name> <name><surname>Forbes</surname><given-names>TL</given-names></name> <name><surname>Leask</surname><given-names>A</given-names></name> <name><surname>Hamilton</surname><given-names>DW</given-names></name></person-group>. <article-title>Inflammatory microenvironment and tumor necrosis factor alpha as modulators of periostin and CCN2 expression in human non-healing skin wounds and dermal fibroblasts</article-title>. <source>Matrix Biol</source>. (<year>2015</year>) <volume>43</volume>:<fpage>71</fpage>&#x2013;<lpage>84</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.matbio.2015.03.003</pub-id>, <pub-id pub-id-type="pmid">25779637</pub-id></mixed-citation></ref>
<ref id="ref23"><label>23.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Eming</surname><given-names>SA</given-names></name> <name><surname>Krieg</surname><given-names>T</given-names></name> <name><surname>Davidson</surname><given-names>JM</given-names></name></person-group>. <article-title>Inflammation in wound repair: molecular and cellular mechanisms</article-title>. <source>J Invest Dermatol</source>. (<year>2007</year>) <volume>127</volume>:<fpage>514</fpage>&#x2013;<lpage>25</lpage>. doi: <pub-id pub-id-type="doi">10.1038/sj.jid.5700701</pub-id></mixed-citation></ref>
<ref id="ref24"><label>24.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Eming</surname><given-names>SA</given-names></name> <name><surname>Murray</surname><given-names>PJ</given-names></name> <name><surname>Pearce</surname><given-names>EJ</given-names></name></person-group>. <article-title>Metabolic orchestration of the wound healing response</article-title>. <source>Cell Metab</source>. (<year>2021</year>) <volume>33</volume>:<fpage>1726</fpage>&#x2013;<lpage>43</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.cmet.2021.07.017</pub-id></mixed-citation></ref>
<ref id="ref25"><label>25.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sharma</surname><given-names>V</given-names></name> <name><surname>Singh</surname><given-names>J</given-names></name> <name><surname>Kumar</surname><given-names>Y</given-names></name> <name><surname>Kumar</surname><given-names>A</given-names></name> <name><surname>Venkatesan</surname><given-names>K</given-names></name> <name><surname>Mukherjee</surname><given-names>M</given-names></name> <etal/></person-group>. <article-title>Integrated insights into gene expression dynamics and transcription factor roles in diabetic and diabetic-infectious wound healing using rat model</article-title>. <source>Life Sci</source>. (<year>2025</year>) <volume>368</volume>:<fpage>123508</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.lfs.2025.123508</pub-id>, <pub-id pub-id-type="pmid">40015667</pub-id></mixed-citation></ref>
<ref id="ref26"><label>26.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Dong</surname><given-names>MW</given-names></name> <name><surname>Li</surname><given-names>M</given-names></name> <name><surname>Chen</surname><given-names>J</given-names></name> <name><surname>Fu</surname><given-names>TT</given-names></name> <name><surname>Lin</surname><given-names>KZ</given-names></name> <name><surname>Ye</surname><given-names>GH</given-names></name> <etal/></person-group>. <article-title>Activation of &#x03B1;7nAChR promotes diabetic wound healing by suppressing AGE-induced TNF-&#x03B1; production</article-title>. <source>Inflammation</source>. (<year>2016</year>) <volume>39</volume>:<fpage>687</fpage>&#x2013;<lpage>99</lpage>. doi: <pub-id pub-id-type="doi">10.1007/s10753-015-0295-x</pub-id>, <pub-id pub-id-type="pmid">26650489</pub-id></mixed-citation></ref>
<ref id="ref27"><label>27.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname><given-names>Q</given-names></name> <name><surname>Zhu</surname><given-names>G</given-names></name> <name><surname>Cao</surname><given-names>X</given-names></name> <name><surname>Dong</surname><given-names>J</given-names></name> <name><surname>Song</surname><given-names>F</given-names></name> <name><surname>Niu</surname><given-names>Y</given-names></name></person-group>. <article-title>Blocking AGE-RAGE signaling improved functional disorders of macrophages in diabetic wound</article-title>. <source>J Diabetes Res</source>. (<year>2017</year>) <volume>2017</volume>:<fpage>1428537</fpage>. doi: <pub-id pub-id-type="doi">10.1155/2017/1428537</pub-id>, <pub-id pub-id-type="pmid">29119117</pub-id></mixed-citation></ref>
<ref id="ref28"><label>28.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hu</surname><given-names>H</given-names></name> <name><surname>Jiang</surname><given-names>H</given-names></name> <name><surname>Ren</surname><given-names>H</given-names></name> <name><surname>Hu</surname><given-names>X</given-names></name> <name><surname>Wang</surname><given-names>X</given-names></name> <name><surname>Han</surname><given-names>C</given-names></name></person-group>. <article-title>AGEs and chronic subclinical inflammation in diabetes: disorders of immune system</article-title>. <source>Diabetes Metab Res Rev</source>. (<year>2015</year>) <volume>31</volume>:<fpage>127</fpage>&#x2013;<lpage>37</lpage>. doi: <pub-id pub-id-type="doi">10.1002/dmrr.2560</pub-id>, <pub-id pub-id-type="pmid">24846076</pub-id></mixed-citation></ref>
<ref id="ref29"><label>29.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Yamamoto</surname><given-names>Y</given-names></name> <name><surname>Yamamoto</surname><given-names>H</given-names></name></person-group>. <article-title>Controlling the receptor for advanced glycation end-products to conquer diabetic vascular complications</article-title>. <source>J Diabetes Investig</source>. (<year>2012</year>) <volume>3</volume>:<fpage>107</fpage>&#x2013;<lpage>14</lpage>. doi: <pub-id pub-id-type="doi">10.1111/j.2040-1124.2011.00191.x</pub-id>, <pub-id pub-id-type="pmid">24843553</pub-id></mixed-citation></ref>
<ref id="ref30"><label>30.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Yan</surname><given-names>SF</given-names></name> <name><surname>Ramasamy</surname><given-names>R</given-names></name> <name><surname>Naka</surname><given-names>Y</given-names></name> <name><surname>Schmidt</surname><given-names>AM</given-names></name></person-group>. <article-title>Glycation, inflammation, and RAGE: a scaffold for the macrovascular complications of diabetes and beyond</article-title>. <source>Circ Res</source>. (<year>2003</year>) <volume>93</volume>:<fpage>1159</fpage>&#x2013;<lpage>69</lpage>. doi: <pub-id pub-id-type="doi">10.1161/01.Res.0000103862.26506.3d</pub-id></mixed-citation></ref>
<ref id="ref31"><label>31.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wong</surname><given-names>SL</given-names></name> <name><surname>Demers</surname><given-names>M</given-names></name> <name><surname>Martinod</surname><given-names>K</given-names></name> <name><surname>Gallant</surname><given-names>M</given-names></name> <name><surname>Wang</surname><given-names>Y</given-names></name> <name><surname>Goldfine</surname><given-names>AB</given-names></name> <etal/></person-group>. <article-title>Diabetes primes neutrophils to undergo NETosis, which impairs wound healing</article-title>. <source>Nat Med</source>. (<year>2015</year>) <volume>21</volume>:<fpage>815</fpage>&#x2013;<lpage>9</lpage>. doi: <pub-id pub-id-type="doi">10.1038/nm.3887</pub-id>, <pub-id pub-id-type="pmid">26076037</pub-id></mixed-citation></ref>
<ref id="ref32"><label>32.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Rui</surname><given-names>S</given-names></name> <name><surname>Dai</surname><given-names>L</given-names></name> <name><surname>Zhang</surname><given-names>X</given-names></name> <name><surname>He</surname><given-names>M</given-names></name> <name><surname>Xu</surname><given-names>F</given-names></name> <name><surname>Wu</surname><given-names>W</given-names></name> <etal/></person-group>. <article-title>Exosomal miRNA-26b-5p from PRP suppresses NETs by targeting MMP-8 to promote diabetic wound healing</article-title>. <source>J Control Release</source>. (<year>2024</year>) <volume>372</volume>:<fpage>221</fpage>&#x2013;<lpage>33</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.jconrel.2024.06.050</pub-id>, <pub-id pub-id-type="pmid">38909697</pub-id></mixed-citation></ref>
<ref id="ref33"><label>33.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Rehak</surname><given-names>L</given-names></name> <name><surname>Giurato</surname><given-names>L</given-names></name> <name><surname>Meloni</surname><given-names>M</given-names></name> <name><surname>Panunzi</surname><given-names>A</given-names></name> <name><surname>Manti</surname><given-names>GM</given-names></name> <name><surname>Uccioli</surname><given-names>L</given-names></name></person-group>. <article-title>The immune-centric revolution in the diabetic foot: monocytes and lymphocytes role in wound healing and tissue regeneration-a narrative review</article-title>. <source>J Clin Med</source>. (<year>2022</year>) <volume>11</volume>:<fpage>889</fpage>. doi: <pub-id pub-id-type="doi">10.3390/jcm11030889</pub-id>, <pub-id pub-id-type="pmid">35160339</pub-id></mixed-citation></ref>
<ref id="ref34"><label>34.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Mazzaferro</surname><given-names>EM</given-names></name> <name><surname>Edwards</surname><given-names>T</given-names></name></person-group>. <article-title>Update on albumin therapy in critical illness</article-title>. <source>Vet Clin North Am Small Anim Pract</source>. (<year>2020</year>) <volume>50</volume>:<fpage>1289</fpage>&#x2013;<lpage>305</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.cvsm.2020.07.005</pub-id>, <pub-id pub-id-type="pmid">32839002</pub-id></mixed-citation></ref>
<ref id="ref35"><label>35.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Khalid</surname><given-names>M</given-names></name> <name><surname>Petroianu</surname><given-names>G</given-names></name> <name><surname>Adem</surname><given-names>A</given-names></name></person-group>. <article-title>Advanced glycation end products and diabetes mellitus: mechanisms and perspectives</article-title>. <source>Biomolecules</source>. (<year>2022</year>) <volume>12</volume>:<fpage>542</fpage>. doi: <pub-id pub-id-type="doi">10.3390/biom12040542</pub-id>, <pub-id pub-id-type="pmid">35454131</pub-id></mixed-citation></ref>
<ref id="ref36"><label>36.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Grada</surname><given-names>A</given-names></name> <name><surname>Phillips</surname><given-names>TJ</given-names></name></person-group>. <article-title>Nutrition and cutaneous wound healing</article-title>. <source>Clin Dermatol</source>. (<year>2022</year>) <volume>40</volume>:<fpage>103</fpage>&#x2013;<lpage>13</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.clindermatol.2021.10.002</pub-id></mixed-citation></ref>
<ref id="ref37"><label>37.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ran</surname><given-names>Q</given-names></name> <name><surname>Xu</surname><given-names>W</given-names></name> <name><surname>Zhao</surname><given-names>X</given-names></name> <name><surname>Sun</surname><given-names>H</given-names></name> <name><surname>Liu</surname><given-names>L</given-names></name> <name><surname>Luo</surname><given-names>Y</given-names></name></person-group>. <article-title>Risk factors for malnutrition in patients with diabetic foot ulcer and its association with prolonged length of hospitalization</article-title>. <source>Nutr Diabetes</source>. (<year>2024</year>) <volume>14</volume>:<fpage>26</fpage>. doi: <pub-id pub-id-type="doi">10.1038/s41387-024-00290-6</pub-id>, <pub-id pub-id-type="pmid">38755177</pub-id></mixed-citation></ref>
</ref-list>
<fn-group>
<fn fn-type="custom" custom-type="edited-by" id="fn0002">
<p>Edited by: <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/3201608/overview">S. Balasubramaniam</ext-link>, Digital University Kerala, India</p>
</fn>
<fn fn-type="custom" custom-type="reviewed-by" id="fn0003">
<p>Reviewed by: <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/3086797/overview">Anto Lourdu Xavier Raj Arockia Selvarathinam</ext-link>, Grand Valley State University, United States</p>
<p><ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/3334268/overview">A. Ashwini</ext-link>, Vel Tech Rangarajan Dr.Sagunthala R&#x0026;D Institute of Science and Technology, India</p>
</fn>
</fn-group>
</back>
</article>