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

ORIGINAL RESEARCH article

Front. Genet., 29 September 2025

Sec. RNA

Volume 16 - 2025 | https://doi.org/10.3389/fgene.2025.1603315

This article is part of the Research TopicRNA Viruses: From Molecular Mechanisms to Novel Therapeutic Targeting StrategiesView all articles

LncRNA and mRNA expression characteristic and bioinformatic analysis in anemic diabetic foot ulcers

Jiayu Lin&#x;Jiayu LinJingying Wang&#x;Jingying WangBo Liang
Bo Liang*
  • Department of Endocrinology, The Second Affiliated Hospital of Fujian Medical University, Quanzhou, Fujian, China

Background: Diabetic foot ulcers (DFUs) are severe complications often accompanied by hemoglobin deficiency. While transcriptomic alterations are implicated in DFU pathogenesis, the roles of long non-coding RNAs (lncRNAs) in linking hemoglobin anomalies to iron-inflammatory dysregulation remain unclear.

Objective: To identify DFU-specific lncRNA-mRNA co-expression networks associated with iron metabolism and inflammation, and evaluate their diagnostic potential.

Methods: RNA sequencing (Illumina NovaSeq 6000) was performed on skin tissues from 5 patients with type 2 diabetes and DFUs and 3 non-diabetic controls. Differentially expressed lncRNAs and mRNAs were identified using DESeq2. Functional pathway enrichment was assessed via Gene Set Enrichment Analysis (GSEA). Correlation networks integrating expression data were constructed to map lncRNA-mRNA interactions.

Results: Analysis revealed 21 differentially expressed lncRNAs (1 upregulated, 20 downregulated) and 368 differentially expressed mRNAs (139 upregulated, 229 downregulated) in DFU tissues. GSEA confirmed significant enrichment of hemoglobin-related pathways, including iron uptake/transport, oxidative stress, and neutrophil degranulation. A co-expression network integrating these findings mapped interactions between 21 lncRNAs and 21 mRNAs. Key lncRNAs showed strong correlations with mRNAs involved in iron dysregulation and chronic inflammation, suggesting their contribution to DFU pathology.

Conclusion: Our analysis of anemia in DFU identifies distinct lncRNA-mRNA co-expression networks linked to hemoglobin deficiency. Key lncRNAs correlate strongly with mRNAs involved in perturbed iron handling and sustained inflammation. This suggests that dysregulation within these networks significantly contributes to the pathology of DFU with concurrent anemia, highlighting potential regulatory roles for these lncRNAs and offering avenues for diagnostic development.

Introduction

Diabetes mellitus (DM) is a group of metabolic disorders characterized by chronic hyperglycemia, resulting from either insufficient insulin secretion or ineffective insulin action. These impairments are influenced by both genetic and environmental factors (American Diabetes Association, 2014). DM is one of the most prevalent metabolic conditions globally, especially in countries like India, China, and the United States. As of 2013, it affected approximately 382 million individuals worldwide (Zimmet et al., 2014), and projections indicate that this number will rise to around 592 million by 2035 (Kharroubi and Darwish, 2015). Among the various complications of DM, diabetic foot ulcers (DFUs) are among the most severe, affecting 15%–25% of individuals with diabetes (Lipsky et al., 2012). Recent studies emphasize the significant burden DFUs place on individuals and healthcare systems due to their high morbidity and frequent hospitalizations (Liu et al., 2023). While neuropathy, ischemia, and infection are established contributors to DFU pathogenesis (Aumiller and Dollahite, 2015), emerging evidence highlights anemia as a critical yet understudied comorbidity in DFU patients. Anemia prevalence in diabetic populations exceeds 30% and correlates with impaired wound healing and ulcer severity (Faghir-Ganji et al., 2024; Humar et al., 2022), suggesting its direct involvement in DFU progression.

Central to this association is the dysregulation of iron homeostasis and chronic inflammation. Iron deficiency disrupts oxygen delivery and collagen synthesis, while iron overload (via impaired recycling) fuels oxidative stress and inflammatory cascades (Wright et al., 2014). In DFUs, aberrant iron handling manifests as elevated hepcidin (a key iron regulator), suppressed erythropoiesis, and persistent neutrophil activation–processes implicated in tissue damage and impaired repair (Alamo et al., 2016). Despite this mechanistic link, the transcriptional drivers connecting hemoglobin anomalies to iron-inflammatory imbalance remain poorly defined.

Long non-coding RNAs (lncRNAs), a subclass of ncRNAs, are typically longer than 200 nucleotides and do not encode proteins (Bus et al., 2024). LncRNAs exert regulatory effects on protein-coding genes through epigenetic, transcriptional, and post-transcriptional mechanisms (Ferrer and Dimitrova, 2024). Long non-coding RNAs (lncRNAs) have emerged as pivotal regulators of metabolic and inflammatory pathways in diabetes (Kuai et al., 2022). Notably, lncRNAs modulate iron metabolism (e.g., LncRNA H19 suppresses hepcidin expression) (Wang et al., 2021) and inflammation (e.g., LncRNA MALAT1 regulates NF-κB signaling) (Hussain et al., 2023). However, their roles in integrating hemoglobin deficiency with iron-inflammatory dysregulation in DFUs are entirely unexplored. Current studies on DFU transcriptomics focus predominantly on protein-coding genes (Su et al., 2024), leaving a critical gap in understanding lncRNA-mediated regulatory networks specific to anemia-associated DFU pathology.

To address this, we performed RNA sequencing of DFU tissues from patients with concurrent anemia and non-diabetic controls. Our study aims to identify DFU-specific lncRNA-mRNA co-expression networks linked to hemoglobin-related pathways (iron uptake/transport, oxidative stress, neutrophil activation), offering novel targets for biomarker development and targeted interventions.

Materials and methods

Skin tissue samples in this study were obtained from two groups: the non-diabetic control (Group A) and the diabetic foot group (Group B). The samples were collected from five DFU patients and three patients with acute wounds who met the inclusion criteria and were treated at the Second Affiliated Hospital of Fujian Medical University between January 2023 and June 2023. The DFU group consisted of three male and two female patients, while the acute wound group included two male and one female patient. These samples were used for subsequent experimental research. All participants provided informed consent and were enrolled according to established criteria. All experiments were performed in accordance with relevant guidelines and regulations in the method section. All experimental procedures were approved by the Ethics Committee of the Second Affiliated Hospital of Fujian Medical University.

Upon admission, skin tissue samples were collected from all participants by trained medical professionals. Relevant demographic data, medical history, and other clinical information were also recorded for each participant.

The inclusion criteria for patients with DFU were those who met the diagnostic criteria for diabetic foot ulcers (Kuai et al., 2022). Exclusion criteria for DFU patients included individuals with any coexisting wound infection, encompassing local infections confined to the ulcer site as well as infections that had spread to become systemic.

For normal skin tissue patients, the inclusion criteria were individuals with wounds who had post-surgical skin tissue free from infection. The exclusion criteria included individuals with sepsis, autoimmune diseases, malignant tumors, or other significant systemic diseases.

Sample collection, RNA extraction, and sequencing

Tissue samples from diabetic foot ulcers and non-diabetic control wounds were collected, with surface liquid removed, and immediately placed into enzyme-free cryogenic vials. The samples were snap-frozen in liquid nitrogen and subsequently stored at −80 °C until further processing. Total RNA was isolated from approximately 30 mg of homogenized tissue using the HiPure Total RNA Mini Kit (Magen Biotechnology, Guangzhou, China, Cat. No. R4130-03) following the manufacturer’s protocol. RNA concentration and purity were determined using a NanoDrop spectrophotometer (Thermo Fisher Scientific, Waltham, MA, United States of America). RNA integrity was assessed using an Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, United States). RNA quality was examined by gel electrophoresis and with Qubit. Only samples with an RNA Integrity Number (RIN) ≥ 8.0, OD260/280 ≥ 1.8, and OD260/230 ≥ 1.0 were used for subsequent library preparation.

Sequencing libraries were constructed from 1 μg of high-quality total RNA per sample using the VAHTS Total RNA-seq (H/M/R) Library Prep Kit for Illumina (Vazyme Biotech, Nanjing, China, Cat. No. NR603). This kit utilizes a probe-based method for ribosomal RNA (rRNA) depletion and includes all steps for cDNA synthesis, end repair, adapter ligation, and PCR amplification. The final libraries were quantified using a Qubit fluorometer (Thermo Fisher Scientific) and their size distribution was confirmed using an Agilent 2100 Bioanalyzer. The qualified libraries were sequenced on an Illumina NovaSeq 6000 platform (Illumina, San Diego, CA, United States) with a paired-end 150 bp (PE150) strategy. The sequencing depth averaged 12.8 Gb per sample, with an average Q30 score of 95.0%.

Bioinformatic analysis

Raw paired-end reads were processed through a standardized pipeline. Adapter sequences and low-quality bases were trimmed using Skewer (v0.2.2). The resulting clean reads were aligned to the human reference genome GRCh38.p13 and its corresponding annotation (Homo_sapiens.GRCh38.109.gtf) from Ensembl using STAR (v2.7.10a) in two-pass mode. A reference-guided transcriptome assembly was then performed for each sample using StringTie (v2.2.1), and the resulting transcripts were merged across all samples to create a unified transcriptome using StringTie’s--merge function. Transcript abundance was quantified against this merged assembly using RSEM (v1.3.3), and expression levels were normalized and reported as Transcripts Per Million (TPM).

Identification and classification of lncRNAs

Transcripts were classified to distinguish lncRNAs from mRNAs. Known lncRNAs were directly identified from the Ensembl annotation (biotypes: lncRNA, antisense). For novel transcripts, a multi-step filtering strategy was applied. First, assembled transcripts were compared to the reference annotation using gffcompare, and unknown intergenic, intronic, and antisense transcripts were classified as preliminary lncRNA candidates. Second, candidates that were single-exon or shorter than 200 nucleotides were filtered out to remove potential artifacts. Finally, the coding potential of the remaining transcripts was assessed using four independent tools: PLEK, CPAT, CNCI, and CPC2. To maximize confidence, only transcripts predicted as non-coding by all four tools were retained as novel lncRNAs. All protein-coding transcripts were classified as mRNAs for subsequent analysis.

Differential expression and enrichment analysis

Differential expression analysis between the DFU and non-diabetic control groups was performed using DESeq2 (v1.10.1). Differentially expressed lncRNAs (DE-lncRNAs) and mRNAs (DE-mRNAs) were identified based on |log2 fold change| ≥ 0.5 and a Benjamini–Hochberg adjusted P-value <0.05. Volcano plots, generated with ggplot2 in R, provided a visual representation of the results, illustrating both the fold change and significance levels of gene expression. Gene Set Enrichment Analysis (GSEA; MSigDB v7.5.1) was applied to DE-mRNAs using pre-ranked log2FC values. The significance thresholds were set as a false discovery rate (FDR) q-value < 0.25 and |normalized enrichment score (NES)| > 1.0. Particular attention was given to pathways such as hemoglobin metabolism, iron transport, oxidative stress, and neutrophil degranulation.

Co-expression network analysis

Differentially expressed lncRNAs and mRNAs obtained by sequencing were used to construct a lncRNA - mRNA co - expression network, aimed at screening for potential lncRNA targets. Center nodes included all significant DE - lncRNAs, while peripheral nodes comprised DE - mRNAs from GSEA - enriched iron metabolism or inflammation pathways (FDR q-value <0.25). Pairwise Pearson correlations (TPM - normalized) with |r| ≥ 0.4 and P-value <0.05 were chosen as the threshold. The networks were constructed and visualized in R using igraph and ggraph, with light blue representing lncRNAs, salmon representing mRNAs, blue edges indicating positive correlations, and red edges indicating negative correlations. Each transcript corresponded to a node, and the lines between nodes represented a strong correlation between them.

Statistical analysis

Statistical analysis was performed using R version 4.4.1. Data processing, differential analysis, visualization, and functional enrichment analysis were conducted using R packages such as DESeq2, ggplot2, pheatmap, and karyoploteR. Differentially expressed genes were selected based on the criteria of P < 0.05 and |log2 fold change| ≥ 0.5. The network plot revealed gene interactions and their roles in biological processes, and was used to explore key genes associated with the disease. The Benjamini–Hochberg method was used to control the False Discovery Rate (FDR) and adjust the raw p-values for multiple hypothesis testing. The significance threshold was set at an adjusted p-value (FDR) < 0.05.

Results

Clinical characteristics

A total of 8 participants were enrolled in the study, comprising 5 patients with type 2 diabetic foot and 3 healthy controls. The mean age of the control group was 69.67 ± 11.85 years, while the mean age of the case group was 58.17 ± 6.40 years. Detailed demographic information for all participants is provided in Table 1.

Table 1
www.frontiersin.org

Table 1. Comparison of clinical characteristics between the two groups.

Principal component analysis (PCA)

PCA of skin transcriptomes revealed a clear separation between diabetic foot ulcer (DFU; n = 5, red labels: B1-B5) and non-diabetic control (n = 3, blue labels: A1, A3) groups along the primary principal component (Dim1, accounting for 21.2% variance), indicating significant global transcriptomic alterations in DFU tissues (Figure 1A).

Figure 1
Panel A shows a PCA plot with individuals categorized into control and DFU groups, distributed along Dim1 and Dim2. Panel B presents a volcano plot depicting gene expression changes, highlighting upregulated and downregulated genes. Panel C features a dot plot indicating enriched biological processes related to gene ontology, with dot size representing count and color indicating p-value significance levels.

Figure 1. Transcriptomic profiling reveals iron-inflammatory dysregulation in DFU. (A) Principal Component Analysis (PCA): Clear separation between control (red) and DFU (green) groups along principal component 1 (Dim1, explaining 16.5%–21.2% variance), indicating distinct transcriptomic signatures. Each point represents an individual tissue sample. (B) mRNA Differential Expression Volcano Plot: Significantly downregulated epithelial/neural genes (DCD, CLSTN1, SCGB2A2, KRT15) and upregulated inflammatory/proteolytic genes (SMARCD1, S100A9, RNF34) in DFU vs. controls. Threshold: |log2 (fold change)| ≥ 0.5 and FDR < 0.05 (gray: non-significant). (C) Bubble Plot of Enriched Pathways: GSEA showing co-enrichment of iron-related and inflammatory pathways. Key pathways: X-axis: Gene ratio (proportion of genes in pathway). Y-axis: Statistical significance (P-value). Bubble size: Number of genes per term (Count). Color intensity: −log10(P-value) (darker = more significant).

Differentially expressed genes

Differential gene expression analysis was conducted using DESeq2, and the results were visualized via a volcano plot. A total of 368 mRNAs were identified as differentially expressed (|log2 FC| ≥ 0.5, P < 0.05), including 139 upregulated genes (e.g., KRT14, CST4, S100A9) and 229 downregulated genes (e.g., SCGB2A2, CLSTN, EDEM3) in DFU compared to controls. These key annotated genes highlight perturbations in critical biological processes such as epidermal barrier function (KRT14), inflammation (S100A9), and stress response (EDEM3) (Figure 1B). Furthermore, 21 lncRNAs were found to be significantly dysregulated in DFU compared to non-diabetic controls, with 20 lncRNAs downregulated and only 1 lncRNA upregulated (LINC01610) (Figure 2A).

Figure 2
Graph A is a scatter plot showing gene expression changes with Log2 fold change on the x-axis and negative Log10 p-value on the y-axis. Points are color-coded for upregulation, downregulation, or stability. Graph B is a network diagram illustrating the correlation strength between various LINC genes, indicated by lines of varying thickness and color.

Figure 2. LncRNA Dysregulation and Iron-Inflammatory Gene Network in DFU with Anemia. (A) Volcano plot of differentially expressed lncRNAs: Significantly upregulated (red; e.g., LINCO1610) and downregulated (blue; e.g., LINCO1152, LINCO2989) lncRNAs in DFU vs. controls (threshold: |log2FC| ≥ 0.5, P < 0.05 by -log10(P-value)). Gray points: non-significant lncRNAs. (B) Correlation network of lncRNAs with iron/inflammation pathway genes: Central nodes (light blue): Key lncRNAs (p < 0.05). Peripheral nodes (salmon): Iron metabolism/inflammation-related mRNAs. Edges: Blue = positive correlation (r ≥ 0.4), Red = negative correlation (r ≤ −0.4).All connections: |r| ≥ 0.4 and P < 0.05 after pathway-based filtering.

Enriched pathways in DFU

Gene Set Enrichment Analysis (GSEA) using Reactome and Gene Ontology Biological Process (GOBP) gene sets revealed significant upregulation of iron-inflammatory pathways in DFU tissues (FDR q < 0.25, |NES| > 1.0). The most enriched pathway was Reactome: Iron uptake and transport (P < 0.005), indicating profound iron dysregulation. Concurrently, inflammatory and immune activation pathways were upregulated, including Reactome: Neutrophil degranulation, GOBP: Inflammatory response, and Reactome: Innate immune system. Oxidative stress processes (GOBP: Regulation of reactive oxygen species metabolic process) and immune cell proliferation (GOBP: Positive regulation of T cell proliferation) were also significantly elevated (Figure 1C).

Co-expression network reveals lncRNA-mRNA regulatory hubs in DFU

The targeted co-expression network integrated 21 differentially expressed lncRNAs (center nodes, light blue) with iron metabolism and inflammation-associated mRNAs (peripheral nodes, salmon), uncovering critical regulatory architectures in DFU pathogenesis. The primary network hub was identified as LINC02397, which exhibited strong positive correlations (r > 0.7, P < 0.01) with multiple pro-inflammatory mediators: cytokine signaling components (IL6R, IL17C), tissue remodeling enzyme MMP8, metabolic sensors (FFAR4, ABCD2), and inflammatory mediator PTGER3. Concurrently, suppressive lncRNA-mRNA axes emerged, most notably with LINC01610 showing significant negative correlations (r < −0.6, P < 0.05) to immunomodulators SEMA7A and CIQTNF12 (Figure 2B).

Discussion

Anemia is prevalent in diabetes mellitus, particularly among patients with diabetic foot ulcers (DFU), potentially impeding healing and increasing amputation/mortality risks due to impaired microcirculation. A meta-analysis of 15 studies (n = 2,895 patients) revealed weighted anemia prevalences of 69.7% (total DFU), 49.5% (mild-to-moderate DFU), and 73.0% (severe DFU) (Yammine et al., 2021). Our transcriptomic analysis provides new insights into the pathogenesis of DFU-associated anemia, highlighting lncRNAs as master regulators of an iron-inflammatory-metabolic axis. Unlike prior models that primarily emphasize hyperglycemia-induced vascular damage or nutritional deficiencies, our findings uncover a previously unrecognized lncRNA-mediated regulatory architecture that drives iron-inflammatory dysregulation in DFU tissues from patients with concomitant anemia.

The significant enrichment of Reactome: Iron uptake and transport alongside inflammatory pathways (Neutrophil degranulation, Inflammatory response) establishes a molecular framework for iron-inflammatory dysregulation in diabetic foot ulcers with hemoglobin deficiency. Our clinical data confirm a significant reduction in hemoglobin levels in diabetic foot ulcer (DFU) patients compared to non-diabetic controls, consistent with the established prevalence of anemia in chronic wounds (Li J. et al., 2023). Transcriptomically, profound downregulation of epithelial defense genes—including secretoglobin family members SCGB2A2, which modulates IL-1β signaling; SCGB1D2, an anti-protease inhibiting neutrophil elastase; and DCD, an antimicrobial peptide critical for skin barrier function—reflects multilayered barrier compromise (Onida et al., 2025). This suppression facilitates bacterial invasion and unchecked proteolytic activity, as evidenced by MMP8 upregulation in our co-expression network (Li S. et al., 2023). Resultant TLR4-driven inflammation propagates IL-6 overproduction through LINC02397-mediated IL6R activation, directly suppressing erythropoiesis and impairing hemoglobin synthesis (Rui et al., 2024). Concurrent neural-epithelial disruption arises from suppression of CLSTN1, essential for sensory nerve communication, and KRT15, maintaining epidermal stemness. These deficits collectively permit uncontrolled IL17C activation (r = 0.718 with LINC02397) and neutrophil-mediated tissue destruction, creating a self-reinforcing cycle of barrier failure, chronic inflammation, erythropoietic suppression, and impaired hemoglobin production that perpetuates DFU progression (Percival et al., 2012; Louiselle et al., 2021; Liu et al., 2025).

The profound transcriptional imbalance in diabetic foot ulcers with anemia—characterized by suppression of 20 lncRNAs regulating epithelial integrity (LINC01752, LINC03050), neural function (LINC00173, LINC01152), and inflammation control (LINC00327, LINC02939), alongside exclusive upregulation of LINC01610—reflects systemic failure of regulatory networks governing tissue homeostasis. Mechanistically, barrier collapse initiates this pathology: co-downregulation of LINC01752 correlates with DCD deficiency, impairing antimicrobial defense; LINC03050 targets the SCGB2A2 promoter, disabling anti-protease activity; and LINC02939 destabilizes KRT15 mRNA, depleting epidermal stem cells. This barrier breach facilitates bacterial invasion, triggering TLR4/MMP8-driven inflammation that induces hepcidin-mediated iron blockade (Shi et al., 2023). Concurrently, LINC01610 upregulation disrupts neuroimmune crosstalk through positive correlation with SEMA7A—amplifying neuropathic pain—and negative regulation of CLSTN1, compromising neuron-keratinocyte communication. These defects collectively disable neurogenic inflammation resolution, perpetuating IL-17 signaling. Critically, suppressed anti-inflammatory lncRNAs exacerbate this cascade: LINC00327 normally targets IL6R 3′UTR to constrain signaling, while LINC02397 binds NF-κB repressors; their loss enables constitutive IL-6/NF-κB activation that promotes hepcidin overexpression and macrophage iron sequestration, completing a vicious cycle of anemia and non-healing (Li et al., 2022; Anguiano-Hernandez et al., 2019).

As the highest-degree node within the co-expression network, LINC02397 emerges as a master regulator of cytokine-driven erythropoietic suppression through its robust positive correlations with core inflammatory effectors. IL6R, whose activation sustains JAK-STAT3 signaling to induce hepcidin—a key hormone that blocks duodenal iron absorption, as established in anemia of inflammation (Yadav et al., 2024). IL17C, where overexpression promotes neutrophil infiltration (Mobile et al., 2023), triggering reactive oxygen species release that directly damages erythrocytes. PTGER3, whose signaling escalates vascular permeability (Morimoto et al., 2014), exacerbating ulcer hypoxia and impairing erythrocyte production. This triad of LINC02397-driven inflammatory dysregulation—iron malabsorption, oxidative erythrocyte damage, and microvascular dysfunction—collectively explains the therapeutic resistance of DFU-associated anemia to iron supplementation, positioning unmitigated inflammation rather than absolute iron deficiency as the primary pathological driver.

Limitations

Our network analysis provides foundational insights into lncRNA-mediated iron-inflammatory dysregulation in DFU-anemia, yet future research must address key gaps: cohort expansion to ensure robustness, functional validation of regulatory mechanisms, resolution of cell-type-specific functions, and integration of dynamic hematological parameters. Only through such multidimensional validation can these findings mature into targeted therapies.

Conclusion

In DFU with anemia, coordinated suppression of 20 lncRNAs disables barrier protection (DCD/SCGB2A2), neural repair (CLSTN1), and inflammation control (IL6R), while upregulation of LINC01610 hijacks neuroimmune signaling (SEMA7A). This dual dysregulation creates a self-sustaining ecosystem where iron becomes biologically inaccessible. Restoring key lncRNAs (e.g., LINC01752, LINC00327) and silencing LINC01610 offers a precision medicine approach to reverse anemia and heal wounds.

Data availability statement

The data presented in the study are deposited in the NCBI repository, accession number PRJNA1234057.

Ethics statement

The studies involving human participants were reviewed and approved by the Ethics Committee of the Second Affiliated Hospital of Fujian Medical University (Approval No. 684). The studies were conducted in accordance with the institutional requirements. Written informed consent was obtained from all participants to participate in this study.

Author contributions

JL: Funding acquisition, Investigation, Writing – original draft, Writing – review and editing. JW: Data curation, Software, Writing – original draft. BL: Project administration, Validation, Writing – review and editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This study was supported by Natural Science Foundation of Fujian Province, China (Grant No.: 2023J01739), the Joint Funds for the Innovation of Science and Technology, Fujian Province (Grant No.: 2024Y9371), and Science and Technology Plan Project of Quanzhou (Grant No.: 2025QZNY082).

Acknowledgments

We thank all of the participants in the study for their cooperation 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.

References

Alamo, I. G., Kannan, K. B., Smith, M. A., Efron, P. A., and Mohr, A. M. (2016). Characterization of erythropoietin and hepcidin in the regulation of persistent injury-associated anemia. J. Trauma Acute Care Surg. 81 (4), 705–712. doi:10.1097/TA.0000000000001163

PubMed Abstract | CrossRef Full Text | Google Scholar

American Diabetes Association (2014). Diagnosis and classification of diabetes mellitus. Diabetes Care 37 (Suppl. 1), S81–S90. doi:10.2337/dc14-S081

PubMed Abstract | CrossRef Full Text | Google Scholar

Anguiano-Hernandez, Y. M., Contreras-Mendez, L., de Los Angeles Hernandez-Cueto, M., Muand Oz-Medina, J. E., Santillan-Verde, M. A., Barbosa-Cabrera, R. E., et al. (2019). Modification of HIF-1α, NF-κB, IGFBP-3, VEGF and adiponectin in diabetic foot ulcers treated with hyperbaric oxygen. Undersea Hyperb. Med. 46 (1), 35–44. doi:10.22462/01.03.2019.4

PubMed Abstract | CrossRef Full Text | Google Scholar

Aumiller, W. D., and Dollahite, H. A. (2015). Pathogenesis and management of diabetic foot ulcers. JAAPA 28 (5), 28–34. doi:10.1097/01.JAA.0000464276.44117.b1

PubMed Abstract | CrossRef Full Text | Google Scholar

Bus, S. A., Armstrong, D. G., Crews, R. T., Gooday, C., Jarl, G., Kirketerp-Moller, K., et al. (2024). Guidelines on offloading foot ulcers in persons with diabetes (IWGDF 2023 update). Diabetes Metab. Res. Rev. 40 (3), e3647. doi:10.1002/dmrr.3647

PubMed Abstract | CrossRef Full Text | Google Scholar

Faghir-Ganji, M., Abdolmohammadi, N., Nikbina, M., Amanollahi, A., Ansari-Moghaddam, A., Rozhan, R., et al. (2024). Prevalence of anemia in patients with diabetes mellitus: a systematic review and meta-analysis. Biomed. Environ. Sci. 37 (1), 96–107. doi:10.3967/bes2024.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Ferrer, J., and Dimitrova, N. (2024). Transcription regulation by long non-coding RNAs: mechanisms and disease relevance. Nat. Rev. Mol. Cell Biol. 25 (5), 396–415. doi:10.1038/s41580-023-00694-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Humar, R., Schaer, D. J., and Vallelian, F. (2022). Erythrophagocytes in hemolytic anemia, wound healing, and cancer. Trends Mol. Med. 28 (11), 906–915. doi:10.1016/j.molmed.2022.08.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Hussain, M. S., Afzal, O., Gupta, G., Altamimi, A. S. A., Almalki, W. H., Alzarea, S. I., et al. (2023). Probing the links: long non-coding RNAs and NF-κB signalling in atherosclerosis. Pathol. Res. Pract. 249, 154773. doi:10.1016/j.prp.2023.154773

PubMed Abstract | CrossRef Full Text | Google Scholar

Kharroubi, A. T., and Darwish, H. M. (2015). Diabetes mellitus: the epidemic of the century. World J. Diabetes 6 (6), 850–867. doi:10.4239/wjd.v6.i6.850

PubMed Abstract | CrossRef Full Text | Google Scholar

Kuai, L., Jiang, J. S., Li, W., Li, B., and Yin, S. Y. (2022). Long non-coding RNAs in diabetic wound healing: current research and clinical relevance. Int. Wound J. 19 (3), 583–600. doi:10.1111/iwj.13655

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, J., Zhang, J., Yang, M., Huang, X., Zhang, M., Fang, X., et al. (2022). Kirenol alleviates diabetic nephropathy via regulating TGF-β/Smads and the NF-κB signal pathway. Pharm. Biol. 60 (1), 1690–1700. doi:10.1080/13880209.2022.2112239

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, J., Zhang, Z., Wei, J., Li, Y., Cheng, C., Ma, S., et al. (2023a). Association between anemia and the risk and outcomes of diabetic foot in patients with type 2 diabetes mellitus. Exp. Ther. Med. 26 (2), 384. doi:10.3892/etm.2023.12083

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, S., Ding, X., Yan, X., Qian, J., and Tan, Q. (2023b). ceAF ameliorates diabetic wound healing by alleviating inflammation and oxidative stress via TLR4/NF-κB and Nrf2 pathways. J. Diabetes Res. 2023, 2422303. doi:10.1155/2023/2422303

PubMed Abstract | CrossRef Full Text | Google Scholar

Lipsky, B. A., Berendt, A. R., Cornia, P. B., Pile, J. C., Peters, E. J. G., Armstrong, D. G., et al. (2012). 2012 Infectious Diseases Society of America clinical practice guideline for the diagnosis and treatment of diabetic foot infections. Clin. Infect. Dis. 54 (12), e132–e173. doi:10.1093/cid/cis346

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Y., Zhang, X., Yang, L., Zhou, S., Li, Y., Shen, Y., et al. (2023). Proteomics and transcriptomics explore the effect of mixture of herbal extract on diabetic wound healing process. Phytomedicine 116, 154892. doi:10.1016/j.phymed.2023.154892

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Z., Bian, X., Luo, L., Björklund, Å. K., Li, L., Zhang, L., et al. (2025). Spatiotemporal single-cell roadmap of human skin wound healing. Cell Stem Cell 32 (3), 479–498.e8. doi:10.1016/j.stem.2024.11.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Louiselle, A. E., Niemiec, S. M., Zgheib, C., and Liechty, K. W. (2021). Macrophage polarization and diabetic wound healing. Transl. Res. 236, 109–116. doi:10.1016/j.trsl.2021.05.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Mobile, R. Z., Mendes, M. C., Machado-Souza, C., Queiroz, P. d. M., Bonfim, C. M. S., Torres-Pereira, C. C., et al. (2023). IL17A and IL17RA gene polymorphisms in Fanconi anemia. Braz Oral Res. 37, e012. doi:10.1590/1807-3107bor-2023.vol37.0012

PubMed Abstract | CrossRef Full Text | Google Scholar

Morimoto, K., Shirata, N., Taketomi, Y., Tsuchiya, S., Segi-Nishida, E., Inazumi, T., et al. (2014). Prostaglandin E2-EP3 signaling induces inflammatory swelling by mast cell activation. J. Immunol. 192 (3), 1130–1137. doi:10.4049/jimmunol.1300290

PubMed Abstract | CrossRef Full Text | Google Scholar

Onida, S., Tan, M., Balan, V., Heatley, F., Peerbux, S., Bolton-Saghdaoui, L., et al. (2025). Decellularized dermis allograft for the treatment of venous leg ulceration: the DAVE RCT. Br. J. Surg. 112 (2), znae330. doi:10.1093/bjs/znae330

PubMed Abstract | CrossRef Full Text | Google Scholar

Percival, S. L., Hill, K. E., Williams, D. W., Hooper, S. J., Thomas, D. W., and Costerton, J. W. (2012). A review of the scientific evidence for biofilms in wounds. Wound Repair Regen. 20 (5), 647–657. doi:10.1111/j.1524-475X.2012.00836.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Rui, S., Dai, L., Zhang, X., He, M., Xu, F., Wu, W., et al. (2024). Exosomal miRNA-26b-5p from PRP suppresses NETs by targeting MMP-8 to promote diabetic wound healing. J. Control Release 372, 221–233. doi:10.1016/j.jconrel.2024.06.050

PubMed Abstract | CrossRef Full Text | Google Scholar

Shi, H., Yuan, X., Liu, G., and Fan, W. (2023). Identifying and validating GSTM5 as an immunogenic gene in diabetic foot ulcer using bioinformatics and machine learning. J. Inflamm. Res. 16, 6241–6256. doi:10.2147/JIR.S442388

PubMed Abstract | CrossRef Full Text | Google Scholar

Su, N., Wang, J., Zhang, H., Jin, H., Miao, B., Zhao, J., et al. (2024). Identification and clinical validation of the role of anoikis-related genes in diabetic foot. Int. Wound J. 21 (3), e14771. doi:10.1111/iwj.14771

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Z., Chen, X., Liu, N., Shi, Y., Liu, Y., Ouyang, L., et al. (2021). A nuclear long non-coding RNA LINC00618 accelerates ferroptosis in a manner dependent upon apoptosis. Mol. Ther. 29 (1), 263–274. doi:10.1016/j.ymthe.2020.09.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Wright, J. A., Richards, T., and Srai, S. K. (2014). The role of iron in the skin and cutaneous wound healing. Front. Pharmacol. 5, 156. doi:10.3389/fphar.2014.00156

PubMed Abstract | CrossRef Full Text | Google Scholar

Yadav, J. P., Singh, A. K., Grishina, M., Pathak, P., Verma, A., Kumar, V., et al. (2024). Insights into the mechanisms of diabetic wounds: pathophysiology, molecular targets, and treatment strategies through conventional and alternative therapies. Inflammopharmacology 32 (1), 149–228. doi:10.1007/s10787-023-01407-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Yammine, K., Hayek, F., and Assi, C. (2021). Is there an association between anemia and diabetic foot ulcers? A systematic review and meta-analysis. Wound Repair Regen. 29 (3), 432–442. doi:10.1111/wrr.12902

PubMed Abstract | CrossRef Full Text | Google Scholar

Zimmet, P. Z., Magliano, D. J., Herman, W. H., and Shaw, J. E. (2014). Diabetes: a 21st century challenge. Lancet Diabetes Endocrinol. 2 (1), 56–64. doi:10.1016/S2213-8587(13)70112-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: bioinformatic analysis, diabetic foot ulcers, long non-coding RNAs, anemia, iron metabolism dysregulation, chronic inflammation

Citation: Lin J, Wang J and Liang B (2025) LncRNA and mRNA expression characteristic and bioinformatic analysis in anemic diabetic foot ulcers. Front. Genet. 16:1603315. doi: 10.3389/fgene.2025.1603315

Received: 31 March 2025; Accepted: 12 September 2025;
Published: 29 September 2025.

Edited by:

Elena Criscuolo, Vita-Salute San Raffaele University, Italy

Reviewed by:

Michael Radzieta, Western Sydney University, Australia
Fang Miao, Shanxi Provincial People’s Hospital, China

Copyright © 2025 Lin, Wang and Liang. 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: Bo Liang, NDI5NTU3MzcwQHFxLmNvbQ==

These authors have contributed equally to this work and share first authorship

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.