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

ORIGINAL RESEARCH article

Front. Pharmacol., 07 November 2025

Sec. Ethnopharmacology

Volume 16 - 2025 | https://doi.org/10.3389/fphar.2025.1666278

A randomized controlled trial of Shengji ointment combined with bromelain in promoting healing of tendon-exposed diabetic foot wounds: integrated 16S rDNA sequencing and metabolomics analysis

Xu Sun,,&#x;Xu Sun1,2,3Jielin Song,,&#x;Jielin Song2,3,4Junchao Sun,,
Junchao Sun1,2,3*Zhaohui Zhang,,
Zhaohui Zhang1,2,3*
  • 1Department of Traditional Chinese Medicine Surgery, The Second Affiliated Hospital of Tianjin University of Traditional Chinese Medicine, Tianjin, China
  • 2Tianjin Institute of Traditional Chinese Medicine Surgery, Tianjin, China
  • 3TCM Institute of Sore and Ulcer, Tianjin University of Traditional Chinese Medicine, Tianjin, China
  • 4Graduate School, Tianjin University of Traditional Chinese Medicine, Tianjin, China

Background: Diabetic foot tendon exposures are prone to infection, necrosis, and prolonged treatment cycle, not only hinder wound healing, but may also lead to amputation and even life-threatening. Utilizing Shengji ointment (a Traditional Chinese Medicine, TCM) combined with bromelain for treating diabetic foot tendon exposure wounds has demonstrated preliminary efficacy. However, the underlying mechanisms, as well as the changes in wound microbiota and metabolites before and after treatment, still warrant further investigation.

Methods: This study used a randomized controlled trial design, with 60 patients randomly assigned to either the experimental group or the control group in a 1:1 ratio. The treatment cycle was 4 weeks. The experimental group was treated with Shengji ointment combined with bromelain, while the control group was treated with hydrocolloid dressing. The clinical efficacy of the two groups was evaluated through a controlled study, while wound exudates were collected only from the experimental group before and after treatment for 16S rDNA sequencing and metabolomics analysis to evaluate changes in wound microbiota and metabolites.

Results: Following the treatment, the granulation tissue coverage, wound healing rate, and Maryland Foot Function Score in the experimental group were markedly superior to those of the control group. Post-treatment analysis revealed significant changes in the wound microbiota composition of the trial group, with a reduction in potential pathogenic bacteria, including Erysipelatoclostridium, Lachnoclostridium, Pontimonas, Hydrogenovibrio, and Aquabacterium, alongside an increase in beneficial bacteria, such as Cetobacterium and Allisonella. Furthermore, 4,034 differential metabolites were identified, with 1907 upregulated and 2,127 downregulated, involving key metabolic pathways such as phenylpropanoid biosynthesis, tyrosine metabolism, and amino acid biosynthesis. Correlation analysis indicated a strong negative association between Pontimonas, Hydrogenovibrio, Aquabacterium, and the majority of the differential metabolites.

Conclusion: The Shengji ointment combined with bromelain has shown notable effectiveness in treating diabetic foot wounds with tendon exposure by modulating microbial composition (reducing pathogenic bacteria and increasing beneficial bacteria) and optimizing the metabolic environment (targeting key metabolic pathways), providing valuable insights for further exploration of its mechanisms.

1 Introduction

Diabetic foot (DF), one of the most prevalent and serious complications of diabetes, is frequently associated with lower limb neuropathy and peripheral vascular disease (Ogbera et al., 2008; Baroni et al., 2014; Bowling et al., 2015; Mariam et al., 2017; Huang et al., 2024). It typically manifests as skin infections, ulcers, and deep tissue damage below the ankle, potentially leading to foot gangrene or even amputation in severe cases (Lipsky, 2024; Dalla Paola and Faglia, 2006; Stancu et al., 2022; Armstrong et al., 2023; Matijević et al., 2023). Diabetic patients face a lifetime risk of 15%–25% for developing diabetic foot ulcers (DFU). For those with a history of previous foot injuries or infections, the incidence of DFU can ascend to 19%–34%, with approximately 17% of such patients requiring amputation (Armstrong et al., 2017; Lipsky et al., 2020). Diabetic foot tendon exposure typically occurs in the mid to late stage of the ulcerative period, with a Wagner classification ranging from 2 to 5 and a TEXAS classification of grade 3 (Wagner, 1981; Lavery et al., 1996). It is manifested as the gradual degeneration and necrosis of the tendon, along with the spread of infection along the tendon towards the proximal end, ultimately resulting in amputation or death, imposing substantial burdens on individuals, healthcare systems, and economies worldwide (Boulton et al., 2005; Raghav et al., 2018; Teobaldi et al., 2018; Behzadifar et al., 2023; McDermott et al., 2023; Zhao et al., 2023a). The current prevailing treatment modality involves surgical excision of the exposed necrotic or degenerated tendon. However, this approach may lead to a partial loss of foot function and inflict damage on the newly formed collateral microcirculation in the surrounding tissues (Sun et al., 2024). Consequently, there is an imperative need to explore a mild and efficacious tendon removal approach for optimizing treatment strategies.

In Traditional Chinese Medicine clinical practice, the Huafu Zaisheng method involves the application of bromelain combined with Shengji ointment to treat diabetic foot wounds with exposed tendons (Zhang et al., 2012a; 2012b). This approach selectively removes degenerated and necrotic tendon tissue without harming healthy structures, thereby promoting the formation of granulation tissue to cover the exposed tendons. Ultimately, it accelerates wound healing while contributing to the preservation of foot functionality. Our previous study has demonstrated that the method of Huafu Zaisheng facilitates tissue regeneration and in situ repair (Zhang et al., 2012c). The exudate from the wound is the outcome of material alterations resulting from the action of the active ingredients of the drug on the wound, representing variations in the microenvironment of the wound (Löffler et al., 2013; Widgerow et al., 2015; Junka et al., 2017; Dong et al., 2022; Gao et al., 2022; Gould and Mahmoudi, 2024). Consequently, investigating the patterns and trends of exudate changes holds significant importance for uncovering the treatment mechanisms. The study suggests that this approach may promote wound healing by modifying wound pH levels and increasing the concentrations of amino acids, proteins, and lysozyme in the wound exudate (Zhang et al., 2011). However, the specific alterations in the microbial community and the underlying metabolic mechanisms in wound exudate before and after intervention with the Huafu Zaisheng method remain inadequately understood.

This study employed 16S rDNA sequencing in conjunction with non-targeted metabolomics to examine the alterations in wound microbiota and metabolites in patients before and after treatment. The objective was to ascertain the correlation between these changes and to elucidate the potential mechanism through which this method facilitates wound healing in diabetic foot tendon exposure. This study was registered at the Chinese Clinical Trial Registry (ChiCTR2000039327) on 23 October 2020.

2 Materials and methods

2.1 Research design and source

This study was a randomized controlled trial conducted from 1 December 2020, to 31 December 2021, designed to evaluate the clinical efficacy and potential mechanism of the Huafu Zaisheng method in treating patients with DFU. All participants were recruited from the Second Affiliated Hospital of Tianjin University of Traditional Chinese Medicine. Prior to enrollment, each patient voluntarily provided written informed consent. The research protocol received approval from the hospital’s Ethics Committee (Approval number: 2020-006-01).

2.2 Inclusion and exclusion criteria

2.2.1 Inclusion criteria

Patients who met the following criteria were included. (1) Patients who met the diagnostic criteria of diabetic foot and Wagner grade 3–4 with tendons exposed; (2) Age between 18 and 85 years; (3) Fasting blood glucose ≤10 mmol/L; (4) Targeted ulcer debridement area between 1 and 20 cm2 (for patients with multiple lesions, the largest ulcer will be the target lesion); (5) An ankle‒brachial index ≥0.5 on the side of the limb with the ulcer; (6) The ulcer has blood, pus, or sticky secretions; (7) Voluntary participation and signing an informed consent form.

2.2.2 Exclusion criteria

Patients who met any of the following criteria were excluded. (1) DFU caused by electrical, chemical, radioactive, neoplastic, or varicose veins, among other reasons, or malignant lesions within the ulcer; (2) There were clinical signs of a systemic infection, such as cellulitis, fever, increased white blood cells, or a positive bacterial culture; (3) Severe uncontrolled hypertension with a systolic blood pressure of ≥160 mmHg or a diastolic blood pressure of ≥110 mmHg; (4) Serum albumin <28 g/L; (5) Hemoglobin <90 g/L; (6) Platelet count <50 × 109/L; (7) Severe heart, liver, or kidney injury, in which case medical treatment may seriously affect patient safety; (8) Women who are pregnant, lactating, recently pregnant, or planning to get pregnant in the near future; (9) Patients with cognitive impairment who could not fully understand the research content or give informed consent; (10) Patients who were allergic to some components of the study drug; (11) Participation in other types of clinical trials within the last month; (12) Poor compliance, inability to complete the study or failure to comply with the study regulations.

2.3 Grouping and intervention

2.3.1 Grouping and sample size calculation

Based on preliminary clinical trial results demonstrating efficacy rates of 83% (treatment group) and 50% (control group), we performed sample size estimation using PASS software (version 21.0). With parameters set at a two-sided α = 0.05, 80% statistical power, and 1:1 allocation ratio, the initial calculation yielded a requirement of 58 participants. Accounting for a potential 10% dropout rate, the final estimated sample size was determined as 64 cases (32 per group).

This study employed a parallel-group randomized controlled trial design. Using the minimization method, a 1:1 dynamic randomization was implemented via the MagMinDA system, with the allocation bias probability set at 1.0 to rigorously balance baseline characteristics between groups. The random allocation sequence was dynamically generated in real-time by the system, with automated encrypted group assignment to guarantee allocation concealment. Neither participants nor researchers were blinded to treatment assignments after randomization. Eligible participants were enrolled by independent study coordinators who had no access to the allocation sequence. Wound exudate samples were collected from the treatment group before (PRE_T) and after (POST_T) therapy for subsequent 16S rDNA sequencing and untargeted metabolomics analysis.

2.3.2 Intervention measures

In the experimental group (Huafu Zaisheng method), bromelain powder was evenly applied to the exposed tendon surface of the wound, followed by external application of Shengji ointment-coated defatted cosmetic cotton dressing. The control group received standard wound care using Comfeel® Hydrocolloid Dressing. The treatment duration for both groups was 4 weeks, with daily dressing changes. Bromelain was obtained from Olive Branch Pharmaceutical Co., Ltd. (Shantou, Guangdong Province, China; specification: 10,000 units; National Drug Approval No.: H44024825; Batch No.: 200,913). Shengji ointment was purchased from Tianjin Darentang Jingwanhong Pharmaceutical Co., Ltd. (Tianjin, China; specification: 30 g/box; National Drug Approval No.: Z12020345; Batch No.: 206,600). Shengji ointment contains Gypsum Fibrosum, Crinis Carbonisatus, Rehmannia glutinosa, Angelica sinensis, calamine, and beeswax as its key components. Comfeel® Hydrocolloid Dressing (debridement gel, 20 g/tube; Coloplast A/S, Denmark; China Medical Device Import Registration Certificate No.: SFDA(I)2008-3640,515; Batch No.: 103,579).

2.4 Outcome indicators

2.4.1 Primary outcome indicator

The primary outcome indicator, granulation tissue coverage rate, is calculated as the percentage of the area covered by granulation tissue relative to the total wound area (Wang et al., 2023; 2024).

2.4.2 Secondary outcome indicator

The secondary outcome indicators include the wound healing rate and the Maryland Foot Function Score. The wound healing rate was defined as the percentage of the difference between the original wound area and the unhealed area relative to the original wound area (Li et al., 2024; Shahat et al., 2024; Thng et al., 2024; Tomoi et al., 2024; Xu et al., 2024). The Maryland Foot Function Score was used to evaluate foot function, including the presence or absence of pain. The maximum score is 100 points, with scores >89 points indicating excellent function, 75–89 points indicating good function, 50–74 points indicating fair function, and <50 points indicating poor function.

2.5 Collection of exudate from wound

A non-invasive technique was utilized. A sterile cotton swab was rotated for 5 s at the granulation tissue’s center to ensure sufficient exudate collection. The swab was subsequently placed in sterile packaging and stored at −80 °C until further analysis.

2.6 16S rDNA analysis

Total microbial DNA was isolated from wound exudate samples using the CTAB method. DNA quality was assessed by 0.8% agarose gel electrophoresis, and concentration was measured using NanoDrop spectrophotometry. The V3-V4 hypervariable regions of bacterial 16S rDNA genes were amplified using primers 341F/806R under standard PCR conditions. Amplification products were verified by 2% agarose gel electrophoresis alongside negative controls. After quantification and purification, equimolar pools were prepared for Illumina NovaSeq 6,000 sequencing (2 × 250 bp paired-end). Raw sequences were processed using the DADA2 pipeline for quality filtering (Q-score ≥20), chimera removal, and ASV generation. Taxonomic classification was performed against the SILVA database (v138) using a Naïve Bayes classifier with 99% confidence threshold. Microbial diversity was assessed using Shannon and Chao1 indices (alpha diversity) and Bray-Curtis dissimilarity (beta diversity). Differential abundance analysis was conducted using Mann-Whitney U tests with FDR correction (p < 0.05). Biomarker identification employed LEfSe analysis. Functional potential was predicted using PICRUSt2 based on MetaCyc pathways.

2.7 Untargeted metabolomics analysis

Metabolites were extracted from wound exudate samples using ice-cold methanol (4:1 v/v methanol:sample) with overnight protein precipitation at −20 °C. After centrifugation (12,000 × g, 10 min, 4 °C), supernatants were concentrated under nitrogen gas and stored at −80 °C until analysis. Chromatographic separation was performed on a UHPLC system (ACQUITY BEH C18 column; 100 × 2.1 mm, 1.8 µm) with a 20-min gradient (5%–95% acetonitrile/0.1% formic acid) at 0.4 mL/min. Metabolites were analyzed via Q Exactive HF-X Orbitrap MS in both ionization modes (70,000 resolution, m/z 70–1,050), with data-dependent MS/MS fragmentation (top 3 precursors, HCD at 28 eV). System stability was monitored using pooled QC samples injected every 10 runs. Raw LC-MS/MS data were processed in XCMS (signal-to-noise >5, RT deviation <0.3 min), with CAMERA for adduct annotation. Metabolites were identified by mass matching (±5 ppm) against HMDB/KEGG (Level 1) and MS/MS spectral matching against an in-house library of 800 standards (Level 2; mass error <0.02 Da, dot product >70%). PCA validated QC stability, while OPLS-DA (VIP≥1.0, P < 0.05) identified differential metabolites (q < 0.05, fold change >2). Enriched KEGG pathways (P < 0.01) were inferred via clusterProfiler.

2.8 Statistical analysis

Results are expressed as mean ± standard deviation (SD) or median with interquartile range (IQR). Normality (Shapiro-Wilk test) and variance homogeneity (Levene’s test) were assessed. Normally distributed data were compared using independent t-tests, while non-normal data were analyzed with Mann-Whitney U tests. Spearman’s correlation was used to examine relationships between metabolites and microbial genera. For handling missing data due to patient dropouts, complete case analysis was performed as primary approach, where participants with missing follow-up data were excluded from analysis. All analyses were conducted in GraphPad Prism 7.0, with p < 0.05 considered statistically significant.

3 Results

During the study, 79 individuals were assessed for eligibility, of whom 13 were excluded. Thirty-three participants were enrolled in each group, and 3 from each group discontinued the intervention. Ultimately, 30 patients per group completed the treatment (Figure 1).

Figure 1
CONSORT 2010 flow diagram for a clinical trial. Enrollment: 79 assessed for eligibility; 13 excluded (5 severe infection, 8 consent withdrawal). Randomization: 66 participants. Allocation: 33 assigned to Shengji Ointment with Bromelain, 33 to Comfeel Plus wound dressing. Follow-up: 3 failed in each group, reasons include withdrawal and lack of communication. Analysis: 30 completed in each group.

Figure 1. CONSORT flowchart of participant enrollment.

3.1 Clinical characteristics and efficacy analysis

An analysis of baseline characteristics, including gender, age, ankle-brachial index, Wagner grading, and disease location, revealed no statistically significant differences between the two groups (P > 0.05) (Table 1), confirming the comparability of treatment outcomes. After 4 weeks of treatment, the experimental group demonstrated significantly higher granulation tissue coverage rate, wound healing rate, and Maryland Foot Function Score compared to the control group (P < 0.05). As the data were non-normally distributed, results are presented as median (IQR) (see Table 2; Figure 2).

Table 1
www.frontiersin.org

Table 1. Baseline characteristics analysis of the two groups.

Table 2
www.frontiersin.org

Table 2. Clinical efficacy analysis of the two groups.

Figure 2
Three box plots labeled A, B, and C compare two groups: HFZS and KHE. Plot A shows granulation coverage percentage, Plot B shows wound healing rate percentage, and Plot C shows the Maryland Foot Function Score percentage. Each plot includes individual data points and indicates a statistically significant difference with asterisk-marked dashed lines between the groups.

Figure 2. Comparison of therapeutic efficacy between two groups. (A) The granulation tissue coverage rates after treatment. (B) The wound healing rates after treatment. (C) Maryland foot function score after treatment. HFZS: Huafu Zaisheng method. KHE: Comfeel® Plus wound dressing. “*” indicates P < 0.05.

3.2 16S rDNA analysis of wound microbiota

3.2.1 Analysis of diversity and composition of wound microbiota before and after treatment with Huafu Zaisheng method

A Venn diagram was generated based on the ASV abundance table (Figure 3A). The analysis identified 1,013 shared ASVs between the two groups, with 1,588 ASVs unique to the pre-treatment group and 1,232 ASVs unique to the post-treatment group. These results indicate that the composition of the wound microbiota was altered following treatment with the Huafu Zaisheng method. To evaluate the α-diversity of wound microbiota before and after treatment, indices including Chao1, Observed species, Goods coverage, Shannon, Simpson, and Pielou e were analyzed. The results showed no significant differences in these indices between the pre- and post-treatment groups (p > 0.05). For simplicity, the Chao1 index, as a representative metric, is presented in Figure 3D, while the results for the other indices are provided in the Supplementary Material. The Simpson index rarefaction curve indicates that as the sampling depth increases, the curve approaches a plateau, suggesting that the sequencing depth was sufficient to capture the majority of microbial diversity (Figure 3C). The Rank Abundance curve spans a wide range on the horizontal axis, indicating high species richness. Meanwhile, the curve exhibits a certain degree of steepness on the vertical axis, suggesting the presence of dominant microbial species in the community (Figure 3E). The beta diversity of wound microbiota before and after treatment was analyzed using PCA. Differences in the shape and distribution observed in the PCA suggest potential variations in microbial community composition between the pre- and post-treatment groups (Figure 3B).

Figure 3
Panel A depicts a Venn diagram with overlapping circles for pre-treatment and post-treatment groups, indicating shared and unique elements. Panel B shows a principal component analysis scatter plot with two groups, pre-treatment and post-treatment, displayed in different colors. Panel C is a line graph showing chao1 diversity index values versus sequence number for pre-treatment and post-treatment. Panel D presents violin plots comparing dispersion in post-treatment and pre-treatment groups. Panel E is a line plot of relative abundance against species rank for various samples, with different colors representing pre- and post-treatment groups.

Figure 3. Diversity and composition analysis of wound microbiota before and after Treatment. (A) ASV Distribution Venn Diagram. (B) PCA plot illustrating beta diversity differences. (C) Rarefaction curves of Simpson index showing sequencing depth adequacy. (D) Representative Chao1 index for alpha diversity. (E) Rank Abundance curve depicting species richness and dominance.

3.2.2 Hierarchical changes in bacterial abundance after Huafu Zaisheng method treatment

Differential analysis was performed for species at each taxonomic level, and significant differential species (P < 0.05) were identified. At the mesophyll level, Thiotrichales and Polyangiales showed a significant decrease after treatment compared to before treatment (Figure 4A). At the family level, Defluviitaleaceae and Piscirickettsiaceae exhibited a significant decrease after treatment compared to before treatment (Figure 4B). At the genus level, after excluding unclassified genera, a significant decrease was observed in Erysipelatoclostridium, Lachnoclostridium, Defluviitaleaceae_UCG-011, Pontimonas, Aquabacterium, Hydrogenovibrio, and Prevotellaceae_UCG-001 after treatment compared to before treatment, while Cetobacterium and Allisonella significantly increased after treatment (Figure 4C). At the species level, after removing unclassified taxa, a significant decrease was observed in Hydrogenovibrio crunogenus and Lactobacillus amylovorus after treatment, while Serratia rubidaea showed an increase (Figure 4D).

Figure 4
Bar charts showing relative abundance of various microbial taxa before (PRE_T) and after treatment (POST_T). Chart A depicts Thiotrichales and Polyangiales showing minimal change. Chart B shows Defluviitaleaceae with an increase and Piscirickettsiaceae with no significant change. Chart C includes several genera like Erysipelatoclostridium and Lactobacillus with varying abundances across treatments. Chart D illustrates species like Hydrogenovibrio crunogenus and Lactobacillus amylovorus, also indicating changes in abundance pre and post-treatment. Error bars indicate variability.

Figure 4. An analysis of the disparities in the wound microbiome before and after treatment at various classification levels. (A) Relative abundance of species at the order level. (B) Relative abundance of species at the family level. (C) Relative abundance of species at the genus level. (D) Relative abundance of species at the species level.

LEfSe was applied to identify differential bacterial taxa with a Linear Discriminant Analysis (LDA) score >1 and p < 0.05 between pre-treatment and post-treatment groups. The cladogram (Figure 5A) illustrates the phylogenetic relationships between distinct microbiota at the phylum to genus level. The letters “p”, “c”, “o”, “f”, and “g” represent phyla, class, order, family, and genus, respectively. The LEfSe analysis results (Figure 5B) indicate that the predominant bacteria in the pre-treatment group were Thiotrichales, Polyangiale, Defluviitaleaceae, Piscirickettsiaceae, Lachnoclostridium, Defluviitaleaceae_UCG-011, Pontimonas, Aquabacterium, Hydrogenovibrio, and Prevotellaceae_UCG-001. These bacteria may be associated with pre-treatment dysbiosis. Following treatment, the group exhibited a notable increase in Cetobacterium and Allisonella, which indicates that the treatment may facilitate wound healing by reducing the prevalence of potential pathogens and promoting the enrichment of bacteria associated with tissue repair.

Figure 5
Cladogram and bar chart comparing microbial taxa from pre-treatment (green) and post-treatment (red) groups. The cladogram shows a circular arrangement of taxa with highlighted differences. The bar chart displays taxa and their LDA scores, indicating relative abundances between the two groups. Key taxa include g_Lachnoclostridium, f_Defluviitaleaceae, and g_Cetobacterium.

Figure 5. Differential species analysis and LDA effect size visualization. (A) Differential species annotation cladogram. (B) Histogram of LDA Effect Size.

3.2.3 Functional prediction of wound microbiome

Based on PICRUSt2 functional predictions, we compared the metabolic functions of the wound microbiome before and after treatment with Huafu Zaisheng method. Analysis based on the COG database indicated that post-treatment, there was a significant downregulation in functions including uridine phosphorylase, asparagine synthetase A, and cobalamin biosynthesis protein (P < 0.05) (Figure 6A). Furthermore, KEGG pathway analysis demonstrated that post-treatment, there was a significant upregulation in functions such as the tricarboxylic acid cycle VII and pyruvate fermentation to butanoate, while the function of mannose biosynthesis was significantly downregulated (P < 0.05) (Figure 6B).

Figure 6
Two bar charts labeled A and B compare mean proportions of various biochemical components. Chart A lists components like

Figure 6. Differential microbial function analysis based on COG (A) and pathway (B).

3.3 Untargeted metabolomics analysis

3.3.1 Data quality control

The metabolite data were acquired using a Q-Exactive high-resolution mass spectrometer with a resolution greater than 30,000. For multi-sample analysis, peak alignment was performed based on small variations in mass-to-charge ratio (m/z) (±0.01) and retention time (rt) (±0.5 min) to ensure accurate quantification of metabolites. Data processing was carried out using XCMS software, and metabolite abundance was visualized by color intensity, with darker colors indicating higher metabolite concentrations. In positive ion mode, the m/z range of 200–600 and retention time between 200 and 400 s exhibited the highest number of detected peaks, suggesting a higher abundance of metabolites in this range and confirming the stability and reliability of the data quality (Figure 7A). In negative ion mode, the m/z range of 0–400, with peak intensities observed in several retention time windows (e.g., 0–100s, 200–300s, and 550–650s), reflected the variation in metabolite abundances across different metabolites (Figure 7B). Overall, the peak alignment ensures precise metabolite quantification, while the distinct peak patterns across different time windows and m/z ranges confirm the stability, consistency, and high reproducibility of the data. These results provide a solid foundation for further metabolomics analysis.

Figure 7
Six-panel image showing data visualizations. Panels A and B depict heatmaps of retention time versus m/z in positive and negative modes, with color gradients indicating density. Panels C and D show PCA plots for positive and negative modes, respectively, with three classes: POST_T, PRE_T, and QC, differentiated by color and enclosed in ellipses. Panels E and F display PCA plots highlighting the POST_T and PRE_T classes. Each plot indicates separation and clustering patterns among the classes.

Figure 7. Data quality analysis. (A,B) The m/z-rt distribution plot of positive/negative ion metabolites, with rt on the x-axis and m/z on the y-axis. (C,D) PCA plots of positive/negative ions including quality control (QC) samples. (E,F) PCA plots of positive/negative ions for pre- and post-treatment groups.

PCA was performed to analyze the metabolite ions identified before and after Huafu Zaisheng method. The quality control samples clustered tightly, indicating stable detection quality (Figures 7C,D). In the PCA plot, some overlap of samples was observed, which may be attributed to biological variability and technical noise. Nevertheless, the metabolite trend plots show a clear separation between the pre- and post-treatment samples, suggesting the presence of metabolic differences between the two groups (Figures 7E,F). These findings provide a reliable foundation for further analysis, indicating that the treatment induced changes in the metabolite profile.

3.3.2 Analysis of differential metabolites

A PLS-DA was performed to model and predict the metabolic expression profiles of samples before and after treatment with the Huafu Zaisheng method. The Variable Importance in Projection (VIP) scores were used to assess the contribution of individual metabolites to the classification. Metabolites with VIP scores ≥1.0 were considered to significantly contribute to the classification and were identified as potential metabolic biomarkers. Separate PLS-DA models were established for each group, and the robustness of the models was validated by 200 permutation tests, evaluating the statistical significance of the explanatory power (R2) and predictive ability (Q2) of the models. The PLS-DA score plots (Figures 8A,B) showed partial overlap in the distribution of samples before and after treatment but revealed a clear separation trend overall, indicating significant metabolic differences between the two groups. The permutation test results (Figures 8C,D) demonstrated that the R2 values were consistently higher than the Q2 values, and the regression line of Q2 intersected the y-axis below zero. These findings indicate that the model exhibits good predictive performance without overfitting. The robustness of the model further supports the reliability and significance of the identified differential metabolites in distinguishing between the two groups. Compared to the pre-treatment group, a total of 4,034 differential metabolites were identified after treatment, including 1,907 upregulated and 2,127 downregulated metabolites (Figure 8F). A heatmap of the top 30 differential metabolites is presented in Figure 8E.

Figure 8
Six-panel image showing various analytical plots. Panel A: PCA plot in positive mode with red and teal clusters, PC1 at 48.61%, and PC2 at 14.14%. Panel B: PCA plot in negative mode with similar layout, PC1 at 47.90%, and PC2 at 13.38%. Panel C: Scatter plot for positive mode with values and correlation shown by R2 and Q2. Panel D: Similar scatter plot for negative mode. Panel E: Heatmap displaying metabolites across POST_T and PRE_T groups. Panel F: Volcano plot indicating significant changes, with 2,127 downregulated and 1,907 upregulated points.

Figure 8. Analysis of differential metabolites. (A,B) PLS-DA score plots of positive and negative ions before and after treatment. (C,D) Permutation test results of the PLS-DA score plots for positive and negative ions. (E,F) Heatmap of the top 30 most significant differential metabolites and volcano plot representing the complete set of differential metabolites.

3.3.3 Analysis of differential pathways

After the KEGG enrichment analysis of differential metabolites, the first 20 enrichment pathways were presented (Figure 9A). It mainly involves the biosynthesis of phenylpropanoids, Tyrosine metabolism, and the biosynthesis of amino acids. The biosynthesis of phenylpropanoids is one of the key differential metabolic pathways. Metabolites related to the biosynthesis of phenylpropanoids, such as p-Coumaroylquinic acid and Esculin, decreased significantly after treatment (P < 0.05) (Figure 9B).

Figure 9
Panel A shows a dot plot of the top 20 KEGG enrichment pathways with varying sizes and colors representing count and -log10 p-values respectively.

Figure 9. Analysis of differential pathways. (A) Bubble map of KEGG enrichment analysis for differential metabolites. (B) Histogram of quantitative outcomes of p-Coumaroylquinic acid and Esculin before and after treatment.

3.4 Correlation analysis of sore flora and metabolites

In this study, correlation analysis takes P < 0.05 and ∣r∣> 0.6 as significance criteria. The analysis of differential genus-level microbial taxa before and after treatment with the Huafu Zaisheng method revealed that Defluviitaleaceae_UCG-011 and Lachnoclostridium exhibited the strongest positive correlation (Figure 10A). Further analysis of differential metabolites before and after treatment showed that scopolin and lactobionic acid exhibited the strongest positive correlation (Figure 10B). In the correlation analysis between differential genera and metabolites, Pontimonas, Hydrogenovibrio, and Aquabacterium showed a significant negative correlation with most differential metabolites, suggesting a pivotal role in metabolite degradation or the regulation of metabolic balance (Figure 10C).

Figure 10
Three correlation heatmaps titled A, B, and C, display relationships between various bacterial genera and compounds. Each map uses a color gradient from red (positive correlation) to blue (negative correlation). Asterisks indicate significance levels: * for p < 0.05, ** for p < 0.01, *** for p < 0.001. Panel A shows correlations among bacterial genera. Panel B illustrates correlations among compounds. Panel C presents correlations between bacterial genera and compounds. The data emphasizes key correlations and significance.

Figure 10. Correlation analysis of sore flora and metabolites. (A) Correlation analysis of differential genus-level microbial taxa. (B) Correlation analysis of differential metabolites. (C) Correlation analysis between differential microbial taxa and differential metabolites. “*” indicates P < 0.05.

4 Discussion

This study, based on clinical research, investigates the efficacy of Huafu Zaisheng method in treating diabetic foot tendon-exposed wounds and explores the mechanisms of microbial communities and metabolites through 16S rDNA sequencing and untargeted metabolomics. The findings indicated that the treatment markedly enhanced the formation of granulation tissue and expedited the healing process, thereby substantiating its clinical efficacy to a considerable extent. Bromelain, a complex mixture of proteolytic enzymes including various proteases such as stem bromelain, fruit bromelain, and ananain, has been demonstrated to promote wound healing through multiple mechanisms, including anti-inflammatory, antioxidant, fibrinolytic, and immunomodulatory effects (Hale et al., 2005; Azarkan et al., 2020). Based on preliminary UPLC fingerprinting analysis, Shengji Ointment contains bioactive components such as azelaic acid, 3-N-butyl phthalate, n-butylphthalide, palmitoleic acid, linolenic acid, ligustilide, linoleic acid, palmitic acid, and oleic acid, which may synergistically contribute to its wound-healing properties (Zhao et al., 2023b). The synergistic wound-healing mechanism of bromelain and Shengji ointment may involve multifaceted cooperative effects: On one hand, the proteolytic activity of bromelain effectively degrades necrotic tissue in wounds, reducing its physical barrier effect, thereby significantly enhancing the penetration of active components in Shengji ointment. On the other hand, Shengji ointment stimulates nutrient secretion and creates a moist wound environment, which not only helps maintain bromelain’s bioactivity but also promotes liquefaction and drainage of necrotic tissue. Notably, existing research has demonstrated that nitric oxide (NO) plays pivotal regulatory roles in diabetic wound healing, including vascular homeostasis modulation, anti-inflammatory effects, and antibacterial activity (Malone-Povolny et al., 2019). Our preliminary experimental results indicate that the combination of bromelain and Shengji ointment significantly increases NO production at wound sites compared to mono-therapy groups, suggesting this phenomenon may serve as a key molecular mechanism underlying their synergistic wound-healing effects (Zhang et al., 2012a). The treatment altered the composition of the local microbial community in treating diabetic foot ulcers with tendon exposure, resulting in a significant reduction of potential pathogenic bacteria and an increase in potentially beneficial bacteria. Anaerobic bacteria, commonly present in deeper lesions, are often linked to fever, foul-smelling wounds, and ulcer depth and duration (Sasikumar et al., 2018; Villa et al., 2024). The observed reduction in anaerobic bacteria, such as Thiotrichales and Defluviitaleaceae, following treatment may suggest a decrease in bacterial virulence and potentially a reduction in wound depth (Coluccio et al., 2024). Erysipelatoclostridium is closely associated with tissue inflammation and severe infections caused by compromised immunity (Zakham et al., 2019; Milosavljevic et al., 2021). Its reduction in wound tissue after treatment may help alleviate local inflammatory responses, thereby enhancing the environment for tissue repair. Previous research has indicated that Lachnoclostridium is positively correlated with the risk of type 2 diabetes. (Li and Li, 2023). Cetobacterium exhibits strong carbohydrate metabolism capabilities (Garibay-Valdez et al., 2024). The reduction in Lachnoclostridium and the increase in Cetobacterium may indicate a shift in the microbial community from a pathological to a healthier state. However, the specific functions of certain bacterial groups, including Defluviitaleaceae, Piscirickettsiaceae, Pontimonas, Aquabacterium, Hydrogenovibrio, and Prevotellaceae_UCG-001, remain insufficiently evidenced. The PICRUSt2 functional prediction analysis indicated that the Huafu Zaisheng method significantly altered the functional potential of wound microbiota. Notably, pathways like the tricarboxylic acid cycle VII and pyruvate fermentation to butanoate were upregulated, whereas mannose biosynthesis was downregulated. These changes suggest that these pathways may be critical to how the Huafu Zaisheng method promotes healing of diabetic tendon-exposed wounds, potentially offering valuable targets for future mechanistic investigations.

Non-targeted metabolomics were used to evaluate the changes of wound metabolites before and after the treatment of decomposition and regeneration. We focus on the first 30 different metabolites. KEGG enrichment analysis shows that different metabolites mainly involve the biosynthesis of phenylpropanoid, Tyrosine metabolism, biosynthesis of amino acids and other pathways. A study has demonstrated that phenylpyruvate, an essential intermediate metabolite in the phenylpropanoid biosynthesis pathway, is significantly accumulated in diabetic foot ulcers. This accumulation contributes to the perpetuation of chronic inflammation and impaired diabetic wound healing by enhancing inflammasome activity and the release of proinflammatory factors (Lv et al., 2023). Based on enrichment results and literature reports, the Huafu Zaisheng method may regulate phenylpyruvate metabolism to alleviate chronic inflammation, while also influencing the overall pathway of phenylpropanoid biosynthesis, thereby optimizing the metabolic microenvironment of tendon-exposed wounds in diabetic foot ulcers and promoting wound healing and tissue regeneration. Amino acids serve as the fundamental building blocks for protein synthesis. They contribute to the production of structural proteins like collagen, which supports the formation of new tissues (Barbul, 2008). Additionally, amino acids play a critical role in regulating cellular metabolic pathways, thereby facilitating the proliferation and migration of fibroblasts and keratinocytes. These processes collectively promote wound healing and tissue regeneration (Kim et al., 2015; Albaugh et al., 2017). The enrichment results of this study suggest that enhancing amino acid biosynthesis may represent a key mechanism by which the Huafu Zaisheng method promotes wound healing. A study has indicated that key genes of tyrosine metabolism, namely, tyrosinase (TYR), tyrosinase-related protein 1 (TYRP1) and dopachrome tautomerase (DCT), are involved in melanogenesis and play an important role in the acute phase of wound healing (Zhu et al., 2021). This study also revealed notable alterations in the tyrosine metabolic pathway, indicating that the Huafu Zaisheng method may potentially target this pathway to facilitate the healing of diabetic foot tendon-exposed wounds.

The study identified a significant reduction in the abundance of Pontimonas, Hydrogenovibrio, and Aquabacterium on the wound following treatment, along with a notable negative correlation between these bacteria and the majority of the differential metabolites. The results of this study indicate that the Huafu Zaisheng method may facilitate wound healing and tissue regeneration by influencing the structure of the microbial community and mitigating the detrimental effects of specific bacteria on metabolites, thereby enhancing the metabolic microenvironment. Further investigation is required to elucidate the precise function of these flora in the metabolic regulation of tendon-exposed wounds in diabetic feet, as well as to identify potential therapeutic targets.

The limitations of this study should be acknowledged. First, the lack of follow-up data precludes evaluation of the intervention’s long-term efficacy and safety, necessitating further investigation. Moreover, while wound exudate analysis provides insights, it may not fully reflect the wound tissue’s metabolic profile or complete microbial composition. An important limitation is the absence of microbial/metabolomic data from hydrogel-treated controls. Previous research has confirmed that certain bioactive hydrogels can increase beneficial microbiota (e.g., Corynebacterium_1 and Lactobacillus) and restore skin microbial balance in diabetic wounds (Hao et al., 2022). Consequently, our inability to compare microbial profiles between treatment groups may obscure potential differences in wound microenvironment modulation and hinder definitive conclusions about our intervention’s specific effects. Additionally, the study population was restricted to a single center with specific inclusion criteria, which may limit the generalisability of findings to broader clinical settings or diverse patient groups. Another limitation is that some differential microbiota and metabolites lack robust literature support, hindering interpretation of their biological relevance. Critically, correlation analyses cannot establish causality between microbial/metabolite changes and therapeutic effects. Future studies should employ targeted metabolomics and animal models to address these mechanistic gaps.

5 Conclusion

Generally, the treatment of tendon-exposed wounds of diabetic foot by Huafu Zaisheng method may facilitate wound healing by regulating the composition of the bacterial community and the metabolic environment. This study concentrated on the alterations in microflora at the generic level. The potential pathogenic bacteria (Erysipelatoclostridium, Lachnoclostridium, Defluviitaleaceae_UCG-011, Pontimonas, Aquabacterium, Hydrogenovibrio, Prevotellaceae_UCG-001) exhibited a notable decline, whereas the possible beneficial bacteria (Cetobacterium and Allisonella) demonstrated a marked increase. The regulation of pivotal metabolic pathways, including phenylpropanoid biosynthesis, tyrosine metabolism, and amino acid biosynthesis, may serve as a mechanism for facilitating healing. In particular, the reduction in the abundance of Pontimonas, Hydrogenovibrio, and Aquabacterium upregulates the expression of beneficial metabolites and thus promotes wound healing. These findings offer new insights into the mechanism of treatment for tendon-exposed wounds of the diabetic foot by Huafu Zaisheng method.

Data availability statement

The metabolomics data presented in this study have been deposited in the Metabolights repository under accession number MTBLS12713, and the 16S r RNA gene sequencing data have been deposited in the NCBI Sequence Read Archive (SRA) under accession number PRJNA1290628.

Ethics statement

The studies involving humans were approved by The Ethics Committee of the Second Affiliated Hospital of Tianjin University of Traditional Chinese Medicine. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

XS: Conceptualization, Data curation, Investigation, Visualization, Writing – original draft, Writing – review and editing. JSo: Conceptualization, Data curation, Investigation, Visualization, Writing – original draft, Writing – review and editing. JSu: Supervision, Writing – original draft, Writing – review and editing. ZZ: Funding acquisition, Supervision, Writing – original draft, Writing – review and editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This research was funded by the National Key Research and Development Program of China (Grant No. 2019YFC1709303), the Construction Project of the TCM Institute of Sore and Ulcer, Tianjin University of Traditional Chinese Medicine (Grant No. 202206), the Construction Project of the Institute of Traditional Chinese Medicine and Integrated Chinese and Western Medicine under the Tianjin Municipal Health Commission (Grant No. 202455), and the Research Planning Projects of the Tianjin Municipal Education Commission (Grant Nos 2024ZD014 and 2024KJ048).

Acknowledgements

We thank all of the subjects who participated in this study.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

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

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2025.1666278/full#supplementary-material

References

Albaugh, V. L., Mukherjee, K., and Barbul, A. (2017). Proline precursors and collagen synthesis: biochemical challenges of nutrient supplementation and wound healing. J. Nutr. 147, 2011–2017. doi:10.3945/jn.117.256404

PubMed Abstract | CrossRef Full Text | Google Scholar

Armstrong, D. G., Boulton, A. J. M., and Bus, S. A. (2017). Diabetic foot ulcers and their recurrence. N. Engl. J. Med. 376, 2367–2375. doi:10.1056/NEJMra1615439

PubMed Abstract | CrossRef Full Text | Google Scholar

Armstrong, D. G., Tan, T.-W., Boulton, A. J. M., and Bus, S. A. (2023). Diabetic foot ulcers: a review. JAMA 330, 62–75. doi:10.1001/jama.2023.10578

PubMed Abstract | CrossRef Full Text | Google Scholar

Azarkan, M., Maquoi, E., Delbrassine, F., Herman, R., M’Rabet, N., Calvo Esposito, R., et al. (2020). Structures of the free and inhibitors-bound forms of bromelain and ananain from Ananas comosus stem and in vitro study of their cytotoxicity. Sci. Rep. 10, 19570. doi:10.1038/s41598-020-76172-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Barbul, A. (2008). Proline precursors to sustain mammalian collagen synthesis. J. Nutr. 138, 2021S–2024S. doi:10.1093/jn/138.10.2021S

PubMed Abstract | CrossRef Full Text | Google Scholar

Baroni, A., Russo, T., and Piccolo, V. (2014). Relationship between local neuroimmune impairment and diabetic foot: the immunocompromised district theory. Int. J. Dermatol 53, 263–266. doi:10.1111/ijd.12328

PubMed Abstract | CrossRef Full Text | Google Scholar

Behzadifar, S., Barras, A., Plaisance, V., Pawlowski, V., Szunerits, S., Abderrahmani, A., et al. (2023). Polymer-based nanostructures for pancreatic beta-cell imaging and non-invasive treatment of diabetes. Pharmaceutics 15, 1215. doi:10.3390/pharmaceutics15041215

PubMed Abstract | CrossRef Full Text | Google Scholar

Boulton, A. J. M., Vileikyte, L., Ragnarson-Tennvall, G., and Apelqvist, J. (2005). The global burden of diabetic foot disease. Lancet 366, 1719–1724. doi:10.1016/S0140-6736(05)67698-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Bowling, F. L., Rashid, S. T., and Boulton, A. J. M. (2015). Preventing and treating foot complications associated with diabetes mellitus. Nat. Rev. Endocrinol. 11, 606–616. doi:10.1038/nrendo.2015.130

PubMed Abstract | CrossRef Full Text | Google Scholar

Coluccio, A., Lopez Palomera, F., and Spero, M. A. (2024). Anaerobic bacteria in chronic wounds: roles in disease, infection and treatment failure. Wound Repair Regen. 32, 840–857. doi:10.1111/wrr.13208

PubMed Abstract | CrossRef Full Text | Google Scholar

Dalla Paola, L., and Faglia, E. (2006). Treatment of diabetic foot ulcer: an overview strategies for clinical approach. Curr. Diabetes Rev. 2, 431–447. doi:10.2174/1573399810602040431

PubMed Abstract | CrossRef Full Text | Google Scholar

Dong, H. F., Huang, X., You, S., and Li, X. H. (2022). Research advances on burn blister fluid. Zhonghua Shao Shang Yu Chuang Mian Xiu Fu Za Zhi 38, 994–998. doi:10.3760/cma.j.cn501120-20211109-00380

PubMed Abstract | CrossRef Full Text | Google Scholar

Gao, L., Liu, X., Zhao, W., Li, C., Wang, F., Gao, J., et al. (2022). Extracellular-matrix-mimicked 3D nanofiber and hydrogel interpenetrated wound dressing with a dynamic autoimmune-derived healing regulation ability based on wound exudate. Biofabrication 15, 015021. doi:10.1088/1758-5090/acaa01

PubMed Abstract | CrossRef Full Text | Google Scholar

Garibay-Valdez, E., Olivas-Bernal, C. A., Vargas-Albores, F., Martínez-Porchas, M., García-Godínez, D. M., Medina-Félix, D., et al. (2024). Deciphering the gut microbiota of zebrafish, the Most used fish as a biological model: a meta-analytic approach. Comp. Biochem. Physiology Part A Mol. & Integr. Physiology 297, 111713. doi:10.1016/j.cbpa.2024.111713

PubMed Abstract | CrossRef Full Text | Google Scholar

Gould, L., and Mahmoudi, M. (2024). Analysis of biogenic amines and small molecule metabolites in human diabetic wound ulcer exudate. ACS Pharmacol. Transl. Sci. 7, 2894–2899. doi:10.1021/acsptsci.4c00418

PubMed Abstract | CrossRef Full Text | Google Scholar

Hale, L. P., Greer, P. K., Trinh, C. T., and James, C. L. (2005). Proteinase activity and stability of natural bromelain preparations. Int. Immunopharmacol. 5, 783–793. doi:10.1016/j.intimp.2004.12.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Hao, M., Ding, C., Sun, S., Peng, X., and Liu, W. (2022). Chitosan/Sodium alginate/velvet antler blood peptides hydrogel promotes diabetic wound healing via regulating angiogenesis, inflammatory response and skin flora. J. Inflamm. Res. 15, 4921–4938. doi:10.2147/JIR.S376692

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, J., Yang, J., Qi, H., Xu, M., Xu, X., and Zhu, Y. (2024). Prediction models for amputation after diabetic foot: systematic review and critical appraisal. Diabetol. Metab. Syndr. 16, 126. doi:10.1186/s13098-024-01360-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Junka, A., Wojtowicz, W., Ząbek, A., Krasowski, G., Smutnicka, D., Bakalorz, B., et al. (2017). Metabolic profiles of exudates from chronic leg ulcerations. J. Pharm. Biomed. Anal. 137, 13–22. doi:10.1016/j.jpba.2017.01.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, D.-G., Kim, E.-Y., Kim, Y.-R., and Kong, I.-S. (2015). Construction of chimeric human epidermal growth factor containing short collagen-binding domain moieties for use as a wound tissue healing agent. J. Microbiol. Biotechnol. 25, 119–126. doi:10.4014/jmb.1405.05073

PubMed Abstract | CrossRef Full Text | Google Scholar

Lavery, L. A., Armstrong, D. G., and Harkless, L. B. (1996). Classification of diabetic foot wounds. J. Foot Ankle Surg.: Off. Publ. Am. Coll. Foot Ankle Surg. 35, 528–531. doi:10.1016/s1067-2516(96)80125-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, H., and Li, C. (2023). Causal relationship between gut microbiota and type 2 diabetes: a two-sample Mendelian randomization study. Front. Microbiol. 14, 1184734. doi:10.3389/fmicb.2023.1184734

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, X., Xiao, X., Zhang, Y., Long, R., Kankala, R. K., Wang, S., et al. (2024). Microneedles based on hyaluronic acid-polyvinyl alcohol with antibacterial, anti-inflammatory, and antioxidant effects promote diabetic wound healing. Int. J. Biol. Macromol. 282, 137185. doi:10.1016/j.ijbiomac.2024.137185

PubMed Abstract | CrossRef Full Text | Google Scholar

Lipsky, B. A.International consensus group on diagnosing and treating the infected diabetic foot (2024). A report from the international consensus on diagnosing and treating the infected diabetic foot. Diabetes Metab. Res. Rev. 20 (Suppl. 1), S68–S77. doi:10.1002/dmrr.453

PubMed Abstract | CrossRef Full Text | Google Scholar

Lipsky, B. A., Senneville, É., Abbas, Z. G., Aragón-Sánchez, J., Diggle, M., Embil, J. M., et al. (2020). Guidelines on the diagnosis and treatment of foot infection in persons with diabetes (IWGDF 2019 update). Diabetes Metab. Res. Rev. 36 (Suppl. 1), e3280. doi:10.1002/dmrr.3280

PubMed Abstract | CrossRef Full Text | Google Scholar

Löffler, M. W., Schuster, H., Bühler, S., and Beckert, S. (2013). Wound fluid in diabetic foot ulceration: more than just an undefined soup? Int. J. Low. Extrem Wounds 12, 113–129. doi:10.1177/1534734613489989

PubMed Abstract | CrossRef Full Text | Google Scholar

Lv, D., Cao, X., Zhong, L., Dong, Y., Xu, Z., Rong, Y., et al. (2023). Targeting phenylpyruvate restrains excessive NLRP3 inflammasome activation and pathological inflammation in diabetic wound healing. Cell Rep. Med. 4, 101129. doi:10.1016/j.xcrm.2023.101129

PubMed Abstract | CrossRef Full Text | Google Scholar

Malone-Povolny, M. J., Maloney, S. E., and Schoenfisch, M. H. (2019). Nitric oxide therapy for diabetic wound healing. Adv. Healthc. Mater. 8, e1801210. doi:10.1002/adhm.201801210

PubMed Abstract | CrossRef Full Text | Google Scholar

Mariam, T. G., Alemayehu, A., Tesfaye, E., Mequannt, W., Temesgen, K., Yetwale, F., et al. (2017). Prevalence of diabetic foot ulcer and associated factors among adult diabetic patients who attend the diabetic Follow-Up clinic at the University of Gondar Referral Hospital, North West Ethiopia, 2016: Institutional-Based cross-sectional Study. J. Diabetes Res. 2017, 2879249. doi:10.1155/2017/2879249

PubMed Abstract | CrossRef Full Text | Google Scholar

Matijević, T., Talapko, J., Meštrović, T., Matijević, M., Erić, S., Erić, I., et al. (2023). Understanding the multifaceted etiopathogenesis of foot complications in individuals with diabetes. World J. Clin. Cases 11, 1669–1683. doi:10.12998/wjcc.v11.i8.1669

PubMed Abstract | CrossRef Full Text | Google Scholar

McDermott, K., Fang, M., Boulton, A. J. M., Selvin, E., and Hicks, C. W. (2023). Etiology, epidemiology, and disparities in the burden of diabetic foot ulcers. Diabetes Care 46, 209–221. doi:10.2337/dci22-0043

PubMed Abstract | CrossRef Full Text | Google Scholar

Milosavljevic, M. N., Kostic, M., Milovanovic, J., Zaric, R. Z., Stojadinovic, M., Jankovic, S. M., et al. (2021). Antimicrobial treatment of Erysipelatoclostridium ramosum invasive infections: a systematic review. Rev. Inst. Med. Trop. Sao Paulo 63, e30. doi:10.1590/S1678-9946202163030

PubMed Abstract | CrossRef Full Text | Google Scholar

Ogbera, O. A., Osa, E., Edo, A., and Chukwum, E. (2008). Common clinical features of diabetic foot ulcers: perspectives from a developing nation. Int. J. Low. Extrem Wounds 7, 93–98. doi:10.1177/1534734608318236

PubMed Abstract | CrossRef Full Text | Google Scholar

Raghav, A., Khan, Z. A., Labala, R. K., Ahmad, J., Noor, S., and Mishra, B. K. (2018). Financial burden of diabetic foot ulcers to world: a progressive topic to discuss always. Ther. Adv. Endocrinol. Metab. 9, 29–31. doi:10.1177/2042018817744513

PubMed Abstract | CrossRef Full Text | Google Scholar

Sasikumar, K., Vijayakumar, C., Jagdish, S., Kadambari, D., Raj Kumar, N., Biswas, R., et al. (2018). Clinico-microbiological profile of septic diabetic foot with special reference to anaerobic infection. Cureus 10, e2252. doi:10.7759/cureus.2252

PubMed Abstract | CrossRef Full Text | Google Scholar

Shahat, M., Hassan, A., Khalil, M. S., Taha, A. G., and Elbadawy, A. (2024). Exploring the other side of the River: early and midterm outcomes of endovascular pedal arch revascularization in patients with chronic limb-threatening ischemia. J. Endovasc. Ther., 15266028241289034. doi:10.1177/15266028241289034

PubMed Abstract | CrossRef Full Text | Google Scholar

Stancu, B., Ilyés, T., Farcas, M., Coman, H. F., Chiş, B. A., and Andercou, O. A. (2022). Diabetic foot complications: a retrospective cohort Study. Int. J. Environ. Res. Public Health 20, 187. doi:10.3390/ijerph20010187

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, X., Jing, J., Dai, R., Zhu, C., Sun, Y., Sun, J., et al. (2024). Shengji ointment combined with bromelain promotes granulation of exposed tendons in diabetic foot ulcers: a multicenter, randomized, positive-controlled clinical trial. Heliyon 10, e39716. doi:10.1016/j.heliyon.2024.e39716

PubMed Abstract | CrossRef Full Text | Google Scholar

Teobaldi, I., Stoico, V., Perrone, F., Bruti, M., Bonora, E., and Mantovani, A. (2018). Honey dressing on a leg ulcer with tendon exposure in a patient with type 2 diabetes. Endocrinol. Diabetes Metab. Case Rep. 2018, 18-0117–0117. doi:10.1530/EDM-18-0117

PubMed Abstract | CrossRef Full Text | Google Scholar

Thng, C. B., Wong, K. L., Wong, A. W.-J., Chew, K. Y., Leow, K., Alvis, L. T. C., et al. (2024). Early results of an orthoplastic multidisciplinary approach to diabetic foot ulcer soft tissue reconstruction. J. Wound Care 33, 926–932. doi:10.12968/jowc.2023.0008

PubMed Abstract | CrossRef Full Text | Google Scholar

Tomoi, Y., Takahara, M., Soga, Y., Hirano, T., Tsunoda, K., Katsuki, T., et al. (2024). Impact of high Lipoprotein(a) levels on delayed wound healing in patients with chronic limb-threatening ischemia after peripheral endovascular therapy. Circ. J. 89, 550–556. doi:10.1253/circj.CJ-24-0383

PubMed Abstract | CrossRef Full Text | Google Scholar

Villa, F., Marchandin, H., Lavigne, J.-P., Schuldiner, S., Cellier, N., Sotto, A., et al. (2024). Anaerobes in diabetic foot infections: pathophysiology, epidemiology, virulence, and management. Clin. Microbiol. Rev. 37, e0014323. doi:10.1128/cmr.00143-23

PubMed Abstract | CrossRef Full Text | Google Scholar

Wagner, F. W. (1981). The dysvascular foot: a system for diagnosis and treatment. Foot Ankle 2, 64–122. doi:10.1177/107110078100200202

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, H.-P., Hu, X.-L., Zhou, H.-Q., Chen, X.-S., and Wang, M.-X. (2023). The safety and effectiveness of a self-made negative pressure suction device in the treatment of chronic wounds. Skin. Res. Technol. 29, e13415. doi:10.1111/srt.13415

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, R., Na, H., Cheng, S., and Zeng, Y. (2024). Effect of Compound Xuejie Powder combined with negative pressure wound therapy on diabetic foot ulcers and blood flow in the dorsal foot artery. Am. J. Transl. Res. 16, 4996–5003. doi:10.62347/XSXY3749

PubMed Abstract | CrossRef Full Text | Google Scholar

Widgerow, A. D., King, K., Tocco-Tussardi, I., Banyard, D. A., Chiang, R., Awad, A., et al. (2015). The burn wound exudate-an under-utilized resource. Burns 41, 11–17. doi:10.1016/j.burns.2014.06.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, J., Tong, Y., Lin, M., Zhang, Z., Li, T., and Zhang, F. (2024). Inhibition of cytokine-like protein 1 transcription hinders wound-healing process in diabetic rats. Diabet. Med. 42, e15459. doi:10.1111/dme.15459

PubMed Abstract | CrossRef Full Text | Google Scholar

Zakham, F., Pillonel, T., Brunel, A.-S., Zambelli, P.-Y., Greub, G., Croxatto, A., et al. (2019). Molecular diagnosis and enrichment culture identified a septic pseudoarthrosis due to an infection with Erysipelatoclostridium ramosum. Int. J. Infect. Dis. 81, 167–169. doi:10.1016/j.ijid.2019.02.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Z., Xu, Q., Ma, J., Liu, X., and Zhu, C. (2011). Experimental Study on the treatment of Huafu Zaisheng in diabetic model rats with Tendon necrosis. Shi Zhen Guo Yi Guo Yao 22, 2787–2790. doi:10.3969/j.issn.1008-0805.2011.11.099

CrossRef Full Text | Google Scholar

Zhang, Z., Xu, Q., Ma, J., Liu, X., and Zhu, C. (2012a). The preliminary mechanism of Huafu Zaisheng method in promoting diabetic wound healing. Chin. J. Basic Med. Tradit. Chin. Med. 18, 378–379+391. doi:10.19945/j.cnki.issn.1006-3250.2012.04.014

CrossRef Full Text | Google Scholar

Zhang, Z., Xu, Q., Ma, J., Liu, X., Zhu, C., Lu, X., et al. (2012b). Clinical research on method of Huafu Zaisheng in treating diabetic foot with tendons worsened. Chin. J. Exp. Tradit. Med. Formulae 18, 276–281. doi:10.13422/j.cnki.syfjx.2012.12.005

CrossRef Full Text | Google Scholar

Zhang, Z., Xu, Q., Ma, J., Liu, X., Zhu, C., Lu, X., et al. (2012c). Preliminary mechanism of the method of Huafu Zaisheng in treating diabetic foot with tendons necrosis. Liaoning J. Tradit. Chin. Med. 39, 1686–1690. doi:10.13192/j.ljtcm.2012.09.29.zhangchh.082

CrossRef Full Text | Google Scholar

Zhao, Y., Dai, X., Sun, X., Zhang, Z., Gao, H., and Gao, R. (2023a). Combination of Shengji ointment and bromelain in the treatment of exposed tendons in diabetic foot ulcers: study protocol for a non-blind, randomized, positive control clinical trial. BMC Complement. Med. Ther. 23, 359. doi:10.1186/s12906-023-04128-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, Y., Li, Z.-H., Sheng, S., Dai, X.-Y., Li, Q.-N., Cao, W.-Y., et al. (2023b). Predictive factors and clinical efficacy of chinese medicine shengji ointment in the treatment of diabetic foot ulcers in the elderly: a prospective study. Front. Pharmacol. 14, 1236229. doi:10.3389/fphar.2023.1236229

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, H.-J., Fan, M., and Gao, W. (2021). Identification of potential hub genes associated with skin wound healing based on time course bioinformatic analyses. BMC Surg. 21, 303. doi:10.1186/s12893-021-01298-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: tendon-exposed wounds in diabetic foot, 16s rDNA sequencing, metabolomics, Shengji ointment, bromelain, Huafu Zaisheng method, diabetic foot ulcer

Citation: Sun X, Song J, Sun J and Zhang Z (2025) A randomized controlled trial of Shengji ointment combined with bromelain in promoting healing of tendon-exposed diabetic foot wounds: integrated 16S rDNA sequencing and metabolomics analysis. Front. Pharmacol. 16:1666278. doi: 10.3389/fphar.2025.1666278

Received: 16 July 2025; Accepted: 21 October 2025;
Published: 07 November 2025.

Edited by:

Javier Echeverria, University of Santiago, Chile

Reviewed by:

Pei Yuen Ng, National University of Malaysia, Malaysia
Lutfi Chabib, Islamic University of Indonesia, Indonesia

Copyright © 2025 Sun, Song, Sun and Zhang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Junchao Sun, NTQyMDk0ODM2QHFxLmNvbQ==; Zhaohui Zhang, enpoNDVAYWxpeXVuLmNvbQ==

These authors have contributed equally to this work

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.