Skip to main content

ORIGINAL RESEARCH article

Front. Endocrinol., 19 November 2024
Sec. Clinical Diabetes

Sequencing of messenger RNA in the healing process of diabetes foot ulcer

Guili WangGuili Wang1Ding WuDing Wu1Donglin LuDonglin Lu1Huifang WuHuifang Wu1Yunmin CaiYunmin Cai2Qingyi MengQingyi Meng1Zhaoxuan Liu*Zhaoxuan Liu1*
  • 1Department of Vascular Surgery, Affiliated Jinan Central Hospital of Shandong First Medical University, Jinan, China
  • 2Department of Vascular and Wound Center, Jinshan Hospital, Fudan University, Shanghai, China

Background: Transcriptome analysis of skin wound tissues from diabetic foot ulcer (DFU) patients to assess changes in the microenvironment during wound healing is performed by messenger RNA (mRNA) sequencing.

Methods: All 5 patients with initial DFU area ≥ 3 cm2 were selected for wound specimen collection at two time points of 0% and 50% wound healing. A total of 10 skin wound samples were obtained for mRNA sequencing. According to the sequencing results, quantitative polymerase chain reaction (qPCR) validation was performed on 12 relevant genes related to angiogenesis, fibroblast proliferation, and wound inflammation. All patients received electrospun poly (L-lactide-co-caprolactone) and formulated porcine fibrinogen (PLCL/Fg) dressing for DFU treatment.

Results: The mRNA sequencing results of DFU skin specimens showed that compared to the 0% and 50% wound healing time points, there were 4347 differentially expressed genes, including 2827 upregulated genes and 1520 downregulated genes. Enrichment analysis of the differentially expressed genes using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) revealed that the upregulated genes were mainly associated with biological processes such as cell adhesion, adhesion junctions, epidermal development, and skin barrier formation. The qPCR analysis results indicated that the increased expression of fibroblast growth factor, vascular endothelial growth factor, and CD200 gene was related to DFU healing.

Conclusion: The healing process of DFU wounds involves the interaction of multiple factors, especially in inflammation control, angiogenesis, and fibroblast proliferation.

1 Introduction

In 2020, approximately 9.3% of adults globally had diabetes; and it is estimated that by 2045, this number will rise to nearly 11% (1). Diabetes is one of the most common chronic diseases in the world today, with high mortality rates and complex complications. Diabetic foot ulcer (DFU) is a severe complication of diabetes and one of the most common complications in diabetic foot patients. Around 15% of diabetes patients will develop DFU, with 84% of these patients undergoing some degree of amputation. In addition to facing significant physical and mental challenges, diabetes patients also bear the financial burden of DFU treatment (2, 3). While clinical methods such as offloading, debridement, and infection control can help treat DFU, the disease remains prevalent and poses risks of amputation or death for patients (4). The pathophysiological relationship between diabetes and impaired skin wound healing is complex. Cells that play a key role in wound healing, when their activity is diminished, can lead to impaired tissue repair in diabetic wounds (5). Keratinocytes and fibroblasts isolated from diabetic foot ulcer wounds exhibit lower proliferative potential and weaker growth factor production capabilities (6). Wound healing is a complex natural process that involves the coordinated actions of various types of cells, enzymes, cytokines, proteins, and hormones. It involves a series of cellular interactions and biochemical processes to repair damaged skin wounds in the human body (7). The process of skin wound repair includes hemostasis, inflammation, cell migration/cell proliferation, and remodeling. The orderly progression of the healing process facilitates the rapid closure of wounds, resulting in well-formed skin morphology or acceptable scarring in the case of acute wounds (8). Chronic wound healing is relatively slow, which is common among diabetic patients. Factors such as persistent inflammation and peripheral nerve damage can exacerbate non-healing wound conditions, increasing the risk of complications. The mechanisms underlying delayed wound healing are multifactorial, including prolonged inflammatory phases and delayed proliferation and remodeling phases. Research by Declue et al. indicates that the delay in diabetic wound healing is associated with the excessive release of pro-inflammatory cytokines such as IL-1β, IL-6, and TNF-α (9). The study by Qiu et al. revealed that diabetic patients with high blood sugar levels exhibit impaired cell proliferation during the wound healing process, along with reduced production of collagen and growth factors (10). Reduced angiogenesis, accompanying decreased levels of growth factors like VEGF and TGF-β1, is also an important reason for impaired wound healing in diabetic patients (11). Rapid and effective wound healing can reduce the adverse complications of DFU. Considering the significant impact of DFU on medical costs and patient survival, there is an urgent need to find new or more effective treatment methods (12). All 5 patients received electrospun poly (L-lactide-co-caprolactone) and formulated porcine fibrinogen (PLCL/Fg) dressing treatment for DFU. Our research group aims to collect skin wound tissues from DFU patients, conduct mRNA sequencing, and perform qPCR analysis to explore the microenvironmental changes and mechanisms of DFU healing.

2 Materials and methods

2.1 Research protocol

The study was conducted at the Vascular & Wound Clinic of Jinshan Hospital, Fudan University, Shanghai, China, from July 1, 2023, to December 26, 2023. All enrolled patients had type 2 diabetes mellitus with concomitant diabetic foot ulcers (DFU). Ethical approval was obtained from the ethics committee of Jinshan Hospital affiliated with Fudan University (authorization no. JIEC 2021-S44). The experiment followed the principles of the Helsinki Declaration (13), and all participants were required to sign a written informed consent form upon enrollment.

2.2 Materials preparation

We collected 5 DFU patients’ wound specimens that met the inclusion criteria: 1. Initial wound area ≥ 3 cm2, 2. Age between 30-80 years old, 3. Duration of DFU wound > 1 month, 4. Stable blood sugar control during hospitalization, 5. Ankle-brachial index (ABI) of the affected limb with DFU > 0.6. Exclusion criteria included patients with wound infections, osteomyelitis, or other diseases that could affect wound healing, such as deep vein thrombosis, systemic lupus erythematosus, or any other systemic inflammatory diseases. Patients with severe limb ischemia, defined as unable to palpate the dorsalis pedis artery or ABI < 0.6, were also excluded. Standard wound care that ulcer wound cleaning and dressing changes were performed, the all wounds were covered with PLCL/Fg dressing of equal size to the wound. At the beginning of the experiment and during each dressing change follow-up, the principal researcher took photos and recorded the wounds. The trial lasted for 12 weeks after the patient was enrolled. According to the specific condition of the wound, participants received PLCL/Fg dressing 1 to 3 outpatient dressing changes per week. Samples were taken at 0% and 50% wound healing, with tissues from the bottom and edges of each DFU wound included at each sampling time. A total of 10 wound samples were obtained from the 5 patients for mRNA sequencing and qPCR analysis. The measurement of wound size was conducted using ImageJ software.

2.3 RNA extraction and library construction

Use TRIzol reagent to extract total RNA according to the manual. Use NanoDrop 2000 spectrophotometer (Thermo Scientific, USA) for RNA purity and quantification, and use Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA) to evaluate RNA integrity. Construct a transcriptome library using the VAHTS Universal V5 RNA-seq Library Prep kit according to the product manual. RNA Quality Check: Ensure the RNA meets the quality and quantity criteria for library preparation. mRNA Isolation: Enrich for mRNA using magnetic beads linked to oligo(dT) primers. First-Strand cDNA Synthesis: mRNA is reverse transcribed to cDNA using random hexamer primers. Second-Strand cDNA Synthesis: The cDNA is synthesized into double-stranded DNA. End Repair and A-Tailing: The double-stranded cDNA is repaired and tailed with an ‘A’ nucleotide to facilitate adapter ligation. Adapter Ligation: Dual indexing adapters are ligated to the cDNA fragments. PCR Amplification: The adapters-ligated cDNA is amplified to enrich the library for sequencing. Library Qualification: The library is quantified and its quality is evaluated using qPCR and/or bioanalyzer before sequencing. These detailed steps ensure comprehensive RNA extraction, assessment, and library preparation, facilitating high-quality RNA sequencing data.

2.4 mRNA sequencing

The Illumina Novaseq 6000 sequencing platform was used to sequence the libraries, generating 150 base pair paired-end reads. Each sample obtained 52.21 ± 3.27 million raw reads. The raw data in fastq format was processed using the fastp software to obtain clean reads by removing low-quality reads for subsequent data analysis (14). The HISAT2 software was used for reference genome alignment, and gene expression levels (FPKM) were calculated (15). The read count of each gene was obtained using HTSeq-count.112 The gene count data was analyzed using the R software (version 3.2.0) for Principal Component Analysis (PCA) and visualization to evaluate the intra-group similarity and inter-group differences of the samples. DESeq2 software was employed for differential gene expression analysis, where genes meeting the criteria of q < 0.05 and fold change > 2 or fold change < 0.5 were defined as differential expression genes (DEGs) (16). DEGs were subjected to hierarchical cluster analysis using R software to display the gene expression patterns in different sample groups. The top 30 genes were visualized using the ggradar package in R software to show the expression changes of upregulated or downregulated genes in a radar plot. Subsequently, gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of DEGs were performed based on hypergeometric distribution algorithm to identify significantly enriched functional terms (17, 18). Differential gene significant enrichment results were visualized using bar charts, chord diagrams, or enrichment analysis bubble charts in R software.

2.5 qPCR experimental

Validation of 12 related genes in vascular neogenesis, fibroblast proliferation, and wound inflammation based on mRNA sequencing differential gene results and previous research on hot genes was performed using qPCR. The specific genes validated were: ACTA2, VEGFC, PDGFA, FGFR2, COL4A2, COL7A1, TNF, IL6, IL10, CD68, CD200, and CD163.

2.5.1 RNA extraction

Total RNA was extracted and quantified using NanoDrop 2000 spectrophotometer (Thermo Scientific, USA) to determine concentration and OD260/OD280 ratio, followed by RNA integrity assessment via agarose gel electrophoresis.

2.5.2 Reverse transcription

The target RNA was reverse transcribed into cDNA using the TransScript All-in-One First-Strand cDNA Synthesis SuperMIX for qPCR kit. The reverse transcription system included: total RNA, 0.5 μg; 5×TransScript All-in-one SuperMix for qPCR, 2 μl; gDNA Remover, 0.5 μl; Nuclease-free H2O up to 10 μl. The reaction conditions were: 42°C for 15 minutes, followed by 85°C for 5 seconds. After reverse transcription, the cDNA was diluted with 90 μl Nuclease-free H2O and stored at -20°C for later use.

2.5.3 Quantitative PCR

The PCR reaction was conducted on a LightCycler® 480 II real-time PCR system (Roche, Swiss) using the PerfectStartTM Green qPCR SuperMix kit (Contains all components required for qPCR reactions, including Taq DNA polymerase, dNTPs, buffer, and fluorescent dye). The reaction components included: 2×PerfectStartTM Green qPCR SuperMix, 5 μl; 10 μM Forward primer, 0.2 μl; 10 μM Reverse primer, 0.2 μl; cDNA, 1 μl; Nuclease-free H2O, 3.6 μl. The PCR program consisted of: 94°C for 30 seconds; followed by 45 cycles of 94°C for 5 seconds, 60°C for 30 seconds. After cycling, a melting curve analysis was performed to verify product specificity by gradually increasing the temperature from 60°C to 97°C with fluorescence signal collection at each degree.

2.5.4 The expression levels of mRNAs were normalized to GAPDH and were calculated using the 2-ΔΔCt method

3 Results

3.1 DFU treatment

The baseline characteristics and DFU features of the 5 enrolled patients are shown in Table 1. After 6-12 weeks of treatment, all patients’ DFU wounds healed without any PLCL/Fg dressing-related complications. Two patients achieved complete healing of DFU wounds by the sixth week of treatment (Characteristic 3 and 5). One patient showed substantially healed DFU wounds by the end of the 12-week experimental period, with only local granulation tissue still exposed (Characteristic 2). The Figure 1 shows that wound conditions of DFU at two sampling time points and the healing status of the wound.

Table 1
www.frontiersin.org

Table 1. Demographic and ulcer characteristics of the participant.

Figure 1
www.frontiersin.org

Figure 1. The wound condition of DFU at wound healing 0% 50% and 100%.

3.2 Cluster analysis and screening of differentially expressed genes

For DFU wound samples at two time points, PCA shows significant differences. The samples of 0% wound healing (blue dots) cluster well, while those of 50% wound healing (red dots) are more dispersed but can still be classified effectively (Figure 2A). Cluster analysis is used to analyze the differences between samples and study their similarity. The analysis results accurately reflect the origin of the experimental samples, with samples from the same source or time point shown as darker blue and clustered together in Figure 2B. The 5 samples from 0% wound healing and 5 samples from 50% wound healing show significant differences between them, while demonstrating good similarity within each group. Figure 2C illustrates that there are 4347 DEGs between the samples of 0% wound healing and 50% wound healing, including 2827 upregulated genes and 1520 downregulated genes. The distribution of differentially expressed genes (DEGs) reveals significant differences between the two groups of samples, which are reflected in a volcano plot as shown in Figure 2D. The top 10 most significantly upregulated and downregulated genes are separately labeled. Figure 2E is a heatmap of differentially expressed gene grouping and clustering based on the expression levels of DEGs. It can reflect the expression patterns of DEGs. The results indicate that there are clear clustering differences in the expression of DEGs between the two groups of samples, confirming both the similarity within groups and the differences between groups.

Figure 2
www.frontiersin.org

Figure 2. Clustering analysis and differential gene screening of two groups. (A) Principal component analysis of 10 samples. (B) Cluster analysis of 10 samples. C1-C5 represent wound healing 0%, K1-K5 represent wound healing 50%. (C) Overview of differential genes. (D) Volcano plot of differential genes, gray represents non-significant differential genes, red and blue represent significant differential genes. (E) Heatmap of differential gene clustering. Labeling the top 10 upregulated and top 10 downregulated genes.

3.3 GO enrichment analysis

Perform GO enrichment analysis on the obtained differential genes and describe the associated functions. Calculate the results to obtain a significant enrichment p-value (calculate each item in biological process, cellular component, and molecular function using Fisher exact test), the lower the value, the more statistically significant the difference. Figures 3A, B respectively list the top 30 items of upregulated and downregulated genes in the GO analysis. The lengths of the bands in three colors represent the importance and difference of the corresponding enriched items. GO enrichment analysis and chord diagram showing the top 10 categories with the smallest q-value or p-value. Similarly, the analysis is based on both upregulated differential genes and downregulated differential genes. The chord diagram in Figure 3C displays the upregulated genes and their corresponding GO biological processes, such as cell adhesion, junction organization, epidermis development, and skin barrier establishment. Additionally, the chord diagram in Figure 3D shows the downregulated genes and their corresponding GO biological processes, such as inflammatory response, immune system processes, and immune response.

Figure 3
www.frontiersin.org

Figure 3. GO enrichment analysis of differential genes. (A) Top 30 GO enrichment terms for upregulated differential genes. The y-axis represents the GO term names, and the x-axis represents -log10p-value. (B) Top 30 GO enrichment terms for downregulated differential genes. (C) Chord diagram analysis of the top 10 GO enrichment terms for upregulated differential genes. The left side shows genes with larger |logFC| in each category, the right side reflects the composition of categories, and the lines in the middle indicate the correspondence between categories and genes. (D) Chord diagram analysis of the top 10 GO enrichment terms for downregulated differential genes.

3.4 KEGG enrichment analysis

KEGG enrichment analysis top 20 (filtering entries corresponding to PopHits ≥ 5 pathways, sorted in descending order by -log10p-value for each entry), as shown in Figure 4. The x-axis represents the enrichment score, with the size of the bubble indicating the number of differential protein coding genes contained in each entry. The bubble color changes from blue-white-yellow-red, with smaller enrichment p-values (red) indicating greater significance for that entry. The chord diagram of KEGG enrichment analysis displays the top 10 classification entries with the smallest q-value or p-value. Similarly, it focuses on analyzing both upregulated and downregulated differential genes. In Figure 4C, the chord diagram shows the upregulated genes and their corresponding KEGG entries, such as tight junctions, vascular smooth muscle contraction, Wnt signaling pathway, and pluripotent stem cell regulation pathway. In Figure 4D, the chord diagram shows the downregulated genes and their corresponding KEGG entries, such as lysosomes, phagosomes, and cell cytotoxicity mediated by natural killer cells.

Figure 4
www.frontiersin.org

Figure 4. KEGG enrichment analysis of differential genes. (A) Results of the top 20 entries for upregulated differential genes. (B) Results of the top 20 entries for downregulated differential genes. (C) Chord diagram analysis of the top 10 KEGG enrichment terms for upregulated differential genes. The left side shows genes with larger logFC in each category, the right side reflects the composition of categories, and the lines in the middle indicate the correspondence between categories and genes. (D) Chord diagram analysis of the top 10 KEGG enrichment terms for downregulated differential genes.

3.5 qPCR analysis

Based on mRNA sequencing results, we selected 12 representative differential genes for qPCR validation. They are: vascular neogenesis ACTA2, VEGFC, and PDGFA; fibroblast proliferation FGFR2, COL4A2, and COL7A1; wound inflammation TNF, IL6, and IL10; macrophage phenotype CD68, CD200, and CD163. Specific quantitative analysis results are shown in Figure 5. Among them, there is a statistical difference between the two groups for six representative genes in the directions of vascular neogenesis and fibroblast proliferation, with P-values less than 0.001. CD200 is a marker gene for M2 macrophages, and its expression is higher at the 50% wound healing time point. For the selected inflammatory representative genes, their expression decreases with wound healing, and there is a statistical difference between the two groups for IL6 (P < 0.05).

Figure 5
www.frontiersin.org

Figure 5. qPCR analysis results of 12 representative differential genes. (A) Genes related to angiogenesis. (B) Genes related to fibroblast proliferation. (C) Genes related to wound inflammation. (D) Genes related to macrophage phenotype. *p<0.05, ***p<0.001.

4 Discussion

We sampled wounds from 5 patients treated with PLCL/Fg dressings at two time points, 0% wound healing and 50% wound healing, to study changes in the wound microenvironment through mRNA sequencing and qPCR analysis. Principal component analysis and clustering analysis revealed significant mRNA expression changes between the two groups of wounds, with good similarity between samples within each group. As DFU patients’ wounds healed, there were significant differences in their wound microenvironment compared to the initial state. The study of differentially expressed genes showed significant differences between the two groups of samples, with 4347 statistically significant DEGs, including 2827 upregulated genes and 1520 downregulated genes. The heat map of clustered differentially expressed genes (Figure 2E), based on their expression levels, showed distinct clustering differences between the two groups of samples, confirming both intra-group sample similarity and inter-group sample differences.

GO enrichment analysis and KEGG pathway analysis revealed functional differences in differentially expressed genes between 0% and 50% wound healing, involving biological processes, cellular components, and molecular functions (Figure 3). GO enrichment analysis and chord diagrams were conducted based on upregulated and downregulated differentially expressed genes. The chord diagram in Figure 3C shows upregulated genes and their corresponding GO biological processes, such as cell adhesion, focal adhesion, epidermal development, and establishment of skin barrier. The enhancement of cellular junction manifests as enhanced intercellular communication or structural reinforcement. The upregulation of genes related to epidermal development may indicate the growth or regeneration process of the epidermis, which further affects the establishment of the skin barrier function and is associated with the enhancement of skin protection functions (19). While the chord diagram in Figure 3D displays downregulated genes and their corresponding GO biological processes, such as inflammatory response, immune system processes, and immune response. Chord diagrams for KEGG pathway analysis (Figures 4C, D) highlight pathway enrichment of upregulated and downregulated differentially expressed genes. Upregulated genes and their corresponding KEGG entries include tight junctions, vascular smooth muscle contraction, Wnt signaling pathway, and pluripotency of stem cells regulation. The enhancement of tight junctions may affect the barrier function between cells, and the increased contraction of vascular smooth muscle may be related to vascular function. The regulation of stem cell function may influence tissue regeneration and repair (20). While downregulated genes and their corresponding KEGG entries, such as lysosomes, phagosomes, and natural killer cell-mediated cytotoxicity. With the healing of DFU wounds, the ability of wound cells to clear foreign substances and metabolic waste decreases (21). It is evident that as DFU wounds heal, pathways such as inflammation and wound clearance are gradually downregulated, while pathways related to wound repair and vascular regeneration are enhanced. These findings are consistent with the study by Rastogi et al., where the use of PLCL/Fg dressings alleviated diabetes-induced wound inflammation by inhibiting inflammatory signaling pathways such as IL-17 and TNF, promoting the transition of macrophages from M1 to M2 phenotype, thereby accelerating wound healing (22).

Previous studies have shown that angiogenesis is an important process in wound healing (2325). The changes in the wound microenvironment during DFU healing are multifaceted, involving wound nutrition support, clearance of metabolic waste, and skin structure regeneration (26, 27). Faster and better angiogenesis is crucial for skin tissue reconstruction (28). We selected 12 differentially expressed genes from three key research directions for qPCR quantitative validation. They are: vascular neogenesis genes ACTA2, VEGFC, and PDGFA; fibroblast proliferation genes FGFR2, COL4A2, and COL7A1; wound inflammation genes TNF, IL6, and IL10; and macrophage phenotype genes CD68, CD200, and CD163. Specific quantitative analysis results are shown in Figure 5. Among them, there is a statistical difference between the six representative genes in the study directions of angiogenesis and fibroblast proliferation, with P<0.001. CD200, a marker gene for M2-type macrophages, shows higher expression at 50% wound healing time point. For the selected inflammatory representative genes, their expressions decrease as the wound heals, and IL6 shows statistical difference between the two groups (P<0.05). It can be observed that during the DFU wound healing process, there is downregulation of inflammation genes and M1 macrophage expression genes, and upregulation of genes related to angiogenesis, fibroblasts, collagen fibers, and M2 macrophages. The results of this study indicate that the angiogenesis-promoting genes α-SMA, VEGFC, and PDGFA are significantly upregulated at the 50% wound healing time point. Additionally, fibroblast proliferation and collagen deposition are crucial for promoting cell proliferation, differentiation, and wound healing (29, 30). Firstly, the changes in the microenvironment of DFU wounds during healing are multifaceted. Secondly, the PLCL/Fg skin regeneration rapidly initiating wound healing from multiple aspects, including angiogenesis, fibroblast proliferation, and alleviation of wound inflammation.

5 Conclusion

Through mRNA sequencing analysis, we found that the changes in the wound microenvironment during DFU healing are multifaceted. Pathways such as cell adhesion, tight junctions, epidermal development, and vascular smooth muscle contraction are upregulated, while pathways related to inflammation response and immune response are downregulated. PLCL/Fg dressings can effectively improve the microenvironment of DFU wounds, promoting rapid wound healing.

Data availability statement

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

Ethics statement

Ethical approval was obtained from the ethics committee of Jinshan Hospital affiliated with Fudan University (authorization no. JIEC 2021-S44). The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation in this study was provided by the participants’ legal guardians/next of kin.

Author contributions

GW: Writing – review & editing, Writing – original draft, Visualization, Validation, Supervision, Software, Resources, Project administration, Methodology, Investigation, Formal analysis, Data curation, Conceptualization. ZL: Writing – review & editing, Resources, Methodology, Investigation, Formal analysis. DW: Writing – original draft, Validation, Resources, Methodology, Conceptualization. DL: Writing – original draft, Resources, Methodology, Investigation, Formal analysis. HW: Writing – review & editing, Validation, Investigation, Data curation. YC: Writing – review & editing, Validation, Investigation, Conceptualization. QM: Writing – review & editing, Writing – original draft, Supervision, Project administration, Methodology, Data curation.

Funding

The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.

Conflict of interest

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

Publisher’s note

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

References

1. Pinchevsky Y, Butkow N, Raal FJ, Chirwa T, Rothberg A. Demographic and clinical factors associated with development of type 2 diabetes: A review of the literature. Int J Gen Med. (2020) 13:121–9. doi: 10.2147/ijgm.s226010

PubMed Abstract | PubMed Abstract | Crossref Full Text | Google Scholar

2. Liu Y, Zhou S, Gao Y, Zhai Y. Electrospun nanofibers as a wound dressing for treating diabetic foot ulcer. Asian J Pharm Sci. (2019) 14:130–43. doi: 10.1016/j.ajps.2018.04.004

PubMed Abstract | PubMed Abstract | Crossref Full Text | Google Scholar

3. Simeni Njonnou SR, Boombhi J, Etoa Etoga MC, Tiodoung Timnou A, Jingi AM, Nkem Efon K, et al. Prevalence of diabetes and associated risk factors among a group of prisoners in the yaoundé Central prison. J Diabetes Res. (2020) 2020:5016327. doi: 10.1155/2020/5016327

PubMed Abstract | PubMed Abstract | Crossref Full Text | Google Scholar

4. Amin N, Doupis J. Diabetic foot disease: From the evaluation of the “foot at risk” to the novel diabetic ulcer treatment modalities. World J Diabetes. (2016) 7:153–64. doi: 10.4239/wjd.v7.i7.153

PubMed Abstract | PubMed Abstract | Crossref Full Text | Google Scholar

5. Kim JB, Chun KW, Han SK, Kim WK. Effect of human bone marrow stromal cell allograft on proliferation and collagen synthesis of diabetic fibroblasts in vitro. J Plast Reconstr Aesthet Surg. (2010) 63:1030–5. doi: 10.1016/j.bjps.2009.04.006

PubMed Abstract | PubMed Abstract | Crossref Full Text | Google Scholar

6. Robson MC. Cytokine manipulation of the wound. Clinics Plast Surg. (2003) 30:57–65. doi: 10.1016/s0094-1298(02)00073-1

PubMed Abstract | PubMed Abstract | Crossref Full Text | Google Scholar

7. Eming SA, Martin P, Tomic-Canic M. Wound repair and regeneration: mechanisms, signaling, and translation. Sci Transl Med. (2014) 6:265sr6. doi: 10.1126/scitranslmed.3009337

PubMed Abstract | PubMed Abstract | Crossref Full Text | Google Scholar

8. Broughton G 2nd, Janis JE, Attinger CE. The basic science of wound healing. Plast Reconstr Surg. (2006) 117:12s–34s. doi: 10.1097/01.prs.0000225430.42531.c2

PubMed Abstract | PubMed Abstract | Crossref Full Text | Google Scholar

9. Declue C, Shornick LP. The cytokine milieu of diabetic wounds. Diabetes Manage. (2015) 5:525–37. doi: 10.2217/dmt.15.44

Crossref Full Text | Google Scholar

10. Qiu Z, Kwon AH, Kamiyama Y. Effects of plasma fibronectin on the healing of full-thickness skin wounds in streptozotocin-induced diabetic rats. J Surg Res. (2007) 138:64–70. doi: 10.1016/j.jss.2006.06.034

PubMed Abstract | PubMed Abstract | Crossref Full Text | Google Scholar

11. Costa PZ, Soares R. Neovascularization in diabetes and its complications. Unraveling the angiogenic paradox. Life Sci. (2013) 92:1037–45. doi: 10.1016/j.lfs.2013.04.001

PubMed Abstract | PubMed Abstract | Crossref Full Text | Google Scholar

12. Elraiyah T, Prutsky G, Domecq JP, Tsapas A, Nabhan M, Frykberg RG, et al. A systematic review and meta-analysis of off-loading methods for diabetic foot ulcers. J Vasc Surg. (2016) 63:59S–68S.e1-2. doi: 10.1016/j.jvs.2015.10.006

PubMed Abstract | PubMed Abstract | Crossref Full Text | Google Scholar

13. World Medical Association. World Medical Association Declaration of Helsinki: ethical principles for medical research involving human subjects. Jama. (2013) 310:2191–4. doi: 10.1001/jama.2013.281053

PubMed Abstract | PubMed Abstract | Crossref Full Text | Google Scholar

14. Chen S, Zhou Y, Chen Y, Gu J. fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics. (2018) 34:i884–90. doi: 10.1093/bioinformatics/bty560

PubMed Abstract | PubMed Abstract | Crossref Full Text | Google Scholar

15. Roberts A, Trapnell C, Donaghey J, Rinn JL, Pachter L. Improving RNA-Seq expression estimates by correcting for fragment bias. Genome Biol. (2011) 12:R22. doi: 10.1186/gb-2011-12-3-r22

PubMed Abstract | PubMed Abstract | Crossref Full Text | Google Scholar

16. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. (2014) 15:550. doi: 10.1186/s13059-014-0550-8

PubMed Abstract | PubMed Abstract | Crossref Full Text | Google Scholar

17. The Gene Ontology Consortium. The Gene Ontology Resource: 20 years and still GOing strong. Nucleic Acids Res. (2019) 47:D330–d338. doi: 10.1093/nar/gky1055

PubMed Abstract | PubMed Abstract | Crossref Full Text | Google Scholar

18. Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, et al. KEGG for linking genomes to life and the environment. Nucleic Acids Res. (2008) 36:D480–4. doi: 10.1093/nar/gkm882

PubMed Abstract | PubMed Abstract | Crossref Full Text | Google Scholar

19. Shadkhoo S, Mani M. The role of intracellular interactions in the collective polarization of tissues and its interplay with cellular geometry. PloS Comput Biol. (2019) 15:e1007454. doi: 10.1371/journal.pcbi.1007454

PubMed Abstract | PubMed Abstract | Crossref Full Text | Google Scholar

20. Leonardo TR, Shi J, Chen D, Trivedi HM, Chen L. Differential expression and function of bicellular tight junctions in skin and oral wound healing. Int J Mol Sci. (2020) 21(8):2966. doi: 10.3390/ijms21082966

PubMed Abstract | PubMed Abstract | Crossref Full Text | Google Scholar

21. Kakanj P, Bhide S, Moussian B, Leptin M. Autophagy-mediated plasma membrane removal promotes the formation of epithelial syncytia. EMBO J. (2022) 41:e109992. doi: 10.15252/embj.2021109992

PubMed Abstract | PubMed Abstract | Crossref Full Text | Google Scholar

22. Rastogi A, Kulkarni SA, Agarwal S, Akhtar M, Arsule S, Bhamre S, et al. Topical esmolol hydrochloride as a novel treatment modality for diabetic foot ulcers: A phase 3 randomized clinical trial. JAMA Netw Open. (2023) 6:e2311509. doi: 10.1001/jamanetworkopen.2023.11509

PubMed Abstract | PubMed Abstract | Crossref Full Text | Google Scholar

23. Theocharidis G, Baltzis D, Roustit M, Tellechea A, Dangwal S, Khetani RS, et al. Integrated skin transcriptomics and serum multiplex assays reveal novel mechanisms of wound healing in diabetic foot ulcers. Diabetes. (2020) 69:2157–69. doi: 10.2337/db20-0188

PubMed Abstract | PubMed Abstract | Crossref Full Text | Google Scholar

24. Veith AP, Henderson K, Spencer A, Sligar AD, Baker AB. Therapeutic strategies for enhancing angiogenesis in wound healing. Adv Drug Delivery Rev. (2019) 146:97–125. doi: 10.1016/j.addr.2018.09.010

PubMed Abstract | PubMed Abstract | Crossref Full Text | Google Scholar

25. Alven S, Peter S, Mbese Z, Aderibigbe BA. Polymer-based wound dressing materials loaded with bioactive agents: potential materials for the treatment of diabetic wounds. Polymers (Basel). (2022) 14(4):724. doi: 10.3390/polym14040724

PubMed Abstract | PubMed Abstract | Crossref Full Text | Google Scholar

26. Liu Y, Liu Y, Deng J, Li W, Nie X. Fibroblast growth factor in diabetic foot ulcer: progress and therapeutic prospects. Front Endocrinol (Lausanne). (2021) 12:744868. doi: 10.3389/fendo.2021.744868

PubMed Abstract | PubMed Abstract | Crossref Full Text | Google Scholar

27. Zhao X, Liu L, An T, Xian M, Luckanagul JA, Su Z, et al. A hydrogen sulfide-releasing alginate dressing for effective wound healing. Acta Biomater. (2020) 104:85–94. doi: 10.1016/j.actbio.2019.12.032

PubMed Abstract | PubMed Abstract | Crossref Full Text | Google Scholar

28. Li S, Su L, Li X, Yang L, Yang M, Zong H, et al. Reconstruction of abdominal wall with scaffolds of electrospun poly (l-lactide-co caprolactone) and porcine fibrinogen: An experimental study in the canine. Mater Sci Eng C Mater Biol Appl. (2020) 110:110644. doi: 10.1016/j.msec.2020.110644

PubMed Abstract | PubMed Abstract | Crossref Full Text | Google Scholar

29. Deng Y, Gao Y, Li T, Xiao S, Adeli M, Rodriguez RD, et al. Amorphizing Metal Selenides-Based ROS Biocatalysts at Surface Nanolayer toward Ultrafast Inflammatory Diabetic Wound Healing. ACS nano. (2023) 17:2943–57. doi: 10.1021/acsnano.2c11448

PubMed Abstract | PubMed Abstract | Crossref Full Text | Google Scholar

30. Lorden ER, Miller KJ, Bashirov L, Ibrahim MM, Hammett E, Jung Y, et al. Mitigation of hypertrophic scar contraction via an elastomeric biodegradable scaffold. Biomaterials. (2015) 43:61–70. doi: 10.1016/j.biomaterials.2014.12.003

PubMed Abstract | PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: diabetic foot ulcers, wound healing, messenger RNA, diabetes, diabetic foot ulcer (DFU)

Citation: Wang G, Wu D, Lu D, Wu H, Cai Y, Meng Q and Liu Z (2024) Sequencing of messenger RNA in the healing process of diabetes foot ulcer. Front. Endocrinol. 15:1468301. doi: 10.3389/fendo.2024.1468301

Received: 21 July 2024; Accepted: 30 October 2024;
Published: 19 November 2024.

Edited by:

Kanhaiya Singh, University of Pittsburgh, United States

Reviewed by:

Mark Dalman, Kent State University, United States
Silvia M. Becerra-Bayona, Autonomous University of Bucaramanga, Colombia

Copyright © 2024 Wang, Wu, Lu, Wu, Cai, Meng and Liu. 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: Zhaoxuan Liu, d2dsNDYxMUAxMjYuY29t

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.