BACKGROUND: Diabetic foot is a complex condition with high incidence, recurrence, mortality, and disability rates. Current treatments for diabetic foot ulcers are often insufficient. This study was conducted to identify potential therapeutic targets for diabetic foot. METHODS: Datasets related to diabetic foot and diabetic skin were retrieved from the Gene Expression Omnibus database. Differentially expressed genes (DEGs) were identified using R software. Enrichment analysis was conducted to screen for critical gene functions and pathways. A protein interaction network was constructed to identify node genes corresponding to key proteins. The DEGs and node genes were overlapped to pinpoint target genes. Plasma and chronic ulcer samples from diabetic and non-diabetic individuals were collected. Western blotting, immunohistochemistry, and enzyme-linked immunosorbent assays were performed to verify the S100 calcium binding protein A9 (S100A9), inflammatory cytokine, and related pathway protein levels. Hematoxylin and eosin staining was used to measure epidermal layer thickness. RESULTS: In total, 283 common DEGs and 42 node genes in diabetic foot ulcers were identified. Forty-three genes were differentially expressed in the skin of diabetic and non-diabetic individuals. The overlapping of the most significant DEGs and node genes led to the identification of S100A9 as a target gene. The S100A9 level was significantly higher in diabetic than in non-diabetic plasma (178.40 ± 44.65 ng/mL vs. 40.84 ± 18.86 ng/mL) and in chronic ulcers, and the wound healing time correlated positively with the plasma S100A9 level. The levels of inflammatory cytokines (tumor necrosis factor-α, interleukin [IL]-1, and IL-6) and related pathway proteins (phospho-extracellular signal regulated kinase [ERK], phospho-p38, phospho-p65, and p-protein kinase B [Akt]) were also elevated. The epidermal layer was notably thinner in chronic diabetic ulcers than in non-diabetic skin (24.17 ± 25.60 μm vs. 412.00 ± 181.60 μm). CONCLUSIONS S100A9 was significantly upregulated in diabetic foot and was associated with prolonged wound healing. S100A9 may impair diabetic wound healing by disrupting local inflammatory responses and skin re-epithelialization.
Calgranulin B/therapeutic use* ; Diabetic Foot/metabolism* ; Humans ; Datasets as Topic ; Computational Biology ; Mice, Inbred C57BL ; Animals ; Mice ; Protein Interaction Maps ; Immunohistochemistry

BACKGROUND: Gluconeogenesis is a critical metabolic pathway for maintaining glucose homeostasis, and its dysregulation can lead to glycometabolic disorders. This study aimed to identify hub biomarkers of these disorders to provide a theoretical foundation for enhancing diagnosis and treatment. METHODS: Gene expression profiles from liver tissues of three well-characterized gluconeogenesis mouse models were analyzed to identify commonly differentially expressed genes (DEGs). Weighted gene co-expression network analysis (WGCNA), machine learning techniques, and diagnostic tests on transcriptome data from publicly available datasets of type 2 diabetes mellitus (T2DM) patients were employed to assess the clinical relevance of these DEGs. Subsequently, we identified hub biomarkers associated with gluconeogenesis-related glycometabolic disorders, investigated potential correlations with immune cell types, and validated expression using quantitative polymerase chain reaction in the mouse models. RESULTS: Only a few common DEGs were observed in gluconeogenesis-related glycometabolic disorders across different contributing factors. However, these DEGs were consistently associated with cytokine regulation and oxidative stress (OS). Enrichment analysis highlighted significant alterations in terms related to cytokines and OS. Importantly, osteomodulin ( OMD ), apolipoprotein A4 ( APOA4 ), and insulin like growth factor binding protein 6 ( IGFBP6 ) were identified with potential clinical significance in T2DM patients. These genes demonstrated robust diagnostic performance in T2DM cohorts and were positively correlated with resting dendritic cells. CONCLUSIONS Gluconeogenesis-related glycometabolic disorders exhibit considerable heterogeneity, yet changes in cytokine regulation and OS are universally present. OMD , APOA4 , and IGFBP6  may serve as hub biomarkers for gluconeogenesis-related glycometabolic disorders.
Machine Learning ; Humans ; Computational Biology/methods* ; Biomarkers/metabolism* ; Diabetes Mellitus, Type 2/genetics* ; Animals ; Mice ; Gluconeogenesis/physiology* ; Gene Expression Profiling ; Transcriptome/genetics* ; Gene Regulatory Networks/genetics* ; Clinical Relevance

With the rapid development of artificial intelligence, computational pathology has been seamlessly integrated into the entire clinical workflow, which encompasses diagnosis, treatment, prognosis, and biomarker discovery. This integration has significantly enhanced clinical accuracy and efficiency while reducing the workload for clinicians. Traditionally, research in this field has depended on the collection and labeling of large datasets for specific tasks, followed by the development of task-specific computational pathology models. However, this approach is labor intensive and does not scale efficiently for open-set identification or rare diseases. Given the diversity of clinical tasks, training individual models from scratch to address the whole spectrum of clinical tasks in the pathology workflow is impractical, which highlights the urgent need to transition from task-specific models to foundation models (FMs). In recent years, pathological FMs have proliferated. These FMs can be classified into three categories, namely, pathology image FMs, pathology image-text FMs, and pathology image-gene FMs, each of which results in distinct functionalities and application scenarios. This review provides an overview of the latest research advancements in pathological FMs, with a particular emphasis on their applications in oncology. The key challenges and opportunities presented by pathological FMs in precision oncology are also explored.
Humans ; Precision Medicine/methods* ; Medical Oncology/methods* ; Artificial Intelligence ; Neoplasms/pathology* ; Computational Biology/methods*

This study aims to investigate the therapeutic effect and mechanism of Daotan Xixin Decoction on APP/PS1 mice. Twelve APP/PS1 male mice were randomized into four groups: APP/PS1 and low-, medium-, and high-dose Daotan Xixin Decoction. Three C57BL/6 wild-type mice were used as the control group. The learning and memory abilities of mice in each group were examined by the Morris water maze test. The pathological changes of hippocampal nerve cells were observed by hematoxylin-eosin staining and Nissl staining. Immunohistochemistry was employed to detect the expression of β-amyloid(Aβ)_(1-42) in the hippocampal tissue. The high-dose Daotan Xixin Decoction group with significant therapeutic effects and the model group were selected for high-throughput sequencing. The differentially expressed gene(DEG) analysis, Gene Ontology(GO) analysis, Kyoto Encyclopedia of Genes and Genomes(KEGG) pathway enrichment analysis, and Gene Set Variation Analysis(GSVA) were performed on the sequencing results. RT-qPCR and Western blot were conducted to determine the mRNA and protein levels, respectively, of some DEGs. Compared with the APP/PS1 group, Daotan Xixin Decoction at different doses significantly improved the learning and memory abilities of APP/PS1 mice, ameliorated the neuropathological damage in the CA1 region of the hippocampus, increased the number of neurons, and decreased the deposition of Aβ_(1-42) in the brain. A total of 1 240 DEGs were screened out, including 634 genes with up-regulated expression and 606 genes with down-regulated expression. The GO analysis predicted the biological processes including RNA splicing and protein folding, the cellular components including spliceosome complexes and nuclear spots, and the molecular functions including unfolded protein binding and heat shock protein binding. The KEGG pathway enrichment analysis revealed the involvement of neurodegenerative disease pathways, amyotrophic lateral sclerosis, and splicing complexes. Further GSVA pathway enrichment analysis showed that the down-regulated pathways involved nuclear factor-κB(NF-κB)-mediated tumor necrosis factor-α(TNF-α) signaling pathway, UV response, and unfolded protein response, while the up-regulated pathways involved the Wnt/β-catenin signaling pathway. The results of RT-qPCR and Western blot showed that compared with the APP/PS1 group, Daotan Xixin Decoction at different doses down-regulated the mRNA and protein levels of signal transducer and activator of transcription 3(STAT3), NF-κB, and interleukin-6(IL-6) in the hippocampus. In conclusion, Daotan Xixin Decoction can improve the learning and memory abilities of APP/PS1 mice by regulating the STAT3/NF-κB/IL-6 signaling pathway.
Animals ; Drugs, Chinese Herbal/administration & dosage* ; Mice ; Male ; Alzheimer Disease/metabolism* ; Computational Biology ; Mice, Inbred C57BL ; High-Throughput Nucleotide Sequencing ; Amyloid beta-Protein Precursor/metabolism* ; Hippocampus/metabolism* ; Mice, Transgenic ; Presenilin-1/metabolism* ; Humans ; Memory/drug effects* ; Maze Learning/drug effects* ; Amyloid beta-Peptides/genetics* ; Disease Models, Animal

With the advancement in big data and artificial intelligence technologies, the extensive application of omics technologies in traditional Chinese medicine(TCM) research has generated large experimental datasets, enabling the exploration of cross-scale correlations among massive data and thereby resulting in the shift toward a data-intensive research paradigm. The emerging approach of multi-omics data fusion analysis, emphasizing technical and computational tools, presents a potential breakthrough in this field. The holistic perspective of TCM aligns with the concept of multi-omics data fusion, yet the data types encountered exhibit high dimensionality with small sample sizes, necessitating data processing techniques such as dimensionality reduction. The current challenge lies in selecting suitable analytical methods for these data to enhance the systematic understanding of physiological functions and disease diagnosis/treatment processes. This paper explores the theories and frameworks of multi-omics data fusion, analyzes methods for fusing high-dimensional, small-sample multi-omics data in TCM, and aims to provide insights for advancing TCM research.
Medicine, Chinese Traditional/methods* ; Humans ; Computational Biology/methods* ; Genomics/methods* ; Sample Size ; Artificial Intelligence ; Multiomics

This study employed bioinformatics to screen the feature genes related to efferocytosis in diabetic kidney disease(DKD) and explores traditional Chinese medicine(TCM) regulating these feature genes. The GSE96804 and GSE30528 datasets were integrated as the training set, and the intersection of differentially expressed genes and efferocytosis-related genes(ERGs) was identified as DKD-ERGs. Subsequently, correlation analysis, protein-protein interaction(PPI) network construction, enrichment analysis, and immune infiltration analysis were performed. Consensus clustering was conducted on DKD patients based on the expression levels of DKD-ERGs, and the expression levels, immune infiltration characteristics, and gene set variations between different subtypes were explored. Eight machine learning models were constructed and their prediction performance was evaluated. The best-performing model was evaluated by nomograms, calibration curves, and external datasets, followed by the identification of efferocytosis-related feature genes associated with DKD. Finally, potential TCMs that can regulate these feature genes were predicted. The results showed that the training set contained 640 differentially expressed genes, and after intersecting with ERGs, 12 DKD-ERGs were obtained, which demonstrated mutual regulation and immune modulation effects. Consensus clustering divided DKD into two subtypes, C1 and C2. The support vector machine(SVM) model had the best performance, predicting that growth arrest-specific protein 6(GAS6), S100 calcium-binding protein A9(S100A9), C-X3-C motif chemokine ligand 1(CX3CL1), 5'-nucleotidase(NT5E), and interleukin 33(IL33) were the feature genes of DKD. Potential TCMs with therapeutic effects included Astragali Radix, Trionycis Carapax, Sargassum, Rhei Radix et Rhizoma, Curcumae Radix, and Alismatis Rhizoma, which mainly function to clear heat, replenish deficiency, activate blood, resolve stasis, and promote urination and drain dampness. Molecular docking revealed that the key components of these TCMs, including β-sitosterol, quercetin, and sitosterol, exhibited good binding activity with the five target genes. These results indicated that efferocytosis played a crucial role in the development and progression of DKD. The feature genes closely related to both DKD and efferocytosis, such as GAS6, S100A9, CX3CL1, NT5E, and IL33, were identified. TCMs such as Astragali Radix, Trionycis Carapa, Sargassum, Rhei Radix et Rhizoma, Curcumae Radix, and Alismatis Rhizoma may provide a new therapeutic strategy for DKD by regulating efferocytosis.
Humans ; Computational Biology ; Diabetic Nephropathies/physiopathology* ; Protein Interaction Maps ; Medicine, Chinese Traditional ; Drugs, Chinese Herbal ; Phagocytosis/genetics* ; Efferocytosis

Chronic cerebral hypoperfusion (CCH) plays an important role in the occurrence and development of vascular dementia (VD). Recent studies have indicated that multiple stages of immune-inflammatory response are involved in the process of cerebral ischemia, drawing increasing attention to immune therapies for cerebral ischemia. This study aims to identify potential immune therapeutic targets for CCH using bioinformatics methods from an immunological perspective. We identified a total of 823 differentially expressed genes associated with CCH, and further screened for 9 core immune-related genes, namely RASGRP1, FGF12, SEMA7A, PAK6, EDN3, BPHL, FCGRT, HSPA1B and MLNR. Gene enrichment analysis showed that core genes were mainly involved in biological functions such as cell growth, neural projection extension, and mesenchymal stem cell migration. Biological signaling pathway analysis indicated that core genes were mainly involved in the regulation of T cell receptor, Ras and MAPK signaling pathways. Through LASSO regression, we identified RASGRP1 and BPHL as key immune-related core genes. Additionally, by integrating differential miRNAs and the miRwalk database, we identified miR-216b-5p as a key immune-related miRNA that regulates RASGRP1. In summary, the predicted miR-216b-5p/ RASGRP1 signaling pathway plays a significant role in immune regulation during CCH, which may provide new targets for immune therapy in CCH.
Humans ; Computational Biology/methods* ; Brain Ischemia/therapy* ; Immunotherapy ; MicroRNAs/genetics* ; Signal Transduction ; Dementia, Vascular/genetics* ; Chronic Disease

Objective To investigate the common molecular features of immunothrombosis, thus enhancing the comprehension of thrombosis triggered by immune and inflammatory responses and offering crucial insights for identifying potential diagnostic and therapeutic targets. Methods Differential gene expression analysis and functional enrichment analysis were conducted on datasets of systemic lupus erythematosus (SLE) and venous thromboembolism (VTE). The intersection of differentially expressed genes in SLE and VTE with those of neutrophil extracellular traps (NET) yielded cross-talk genes (CG) for SLE-NET and VTE-NET interaction. Further analysis included functional enrichment and protein-protein interaction (PPI) network assessments of these CG to identify hub genes. Venn diagrams and receiver operating characteristic (ROC) curve analysis were employed to pinpoint the most effective shared diagnostic CG, which were validated using a graft-versus-host disease (GVHD) dataset. Results Differential expression genes in SLE and VTE were associated with distinct biological processes, whereas SLE-NET-CG and VTE-NET-CG were implicated in pathways related to leukocyte migration, inflammatory response, and immune response. Through PPI network analysis, several hub genes were identified, with matrix metalloproteinase 9 (MMP9) and S100 calcium-binding protein A12 (S100A12) emerging as the best shared diagnostic CG for SLE (AUC: 0.936 and 0.832) and VTE (AUC: 0.719 and 0.759). Notably, MMP9 exhibited good diagnostic performance in the GVHD dataset (AUC: 0.696). Conclusion This study unveils the common molecular features of SLE, VTE, and NET, emphasizing MMP9 and S100A12 as the optimal shared diagnostic CG, thus providing valuable evidence for the diagnosis and therapeutic strategies related to immunothrombosis. Additionally, the expression of MMP9 in GVHD highlights its critical role in the risk of VTE associated with immune system disorders.
Humans ; Computational Biology/methods* ; Lupus Erythematosus, Systemic/immunology* ; Protein Interaction Maps/genetics* ; Venous Thromboembolism/therapy* ; Matrix Metalloproteinase 9/genetics* ; Extracellular Traps/metabolism* ; Gene Regulatory Networks ; Thrombosis/immunology* ; Graft vs Host Disease/genetics* ; Gene Expression Profiling

Objective To explore macrophage subpopulations in ischemic stroke (IS) by using single-cell RNA sequencing (scRNA-seq) data analysis and High-Dimensional Weighted Gene Co-Expression Network Analysis (hdWGCNA). Methods Based on single-cell sequencing data, transcriptomic information for different cell types was obtained, and macrophages were selected for subpopulation identification. hdWGCNA, cell-cell communication, and pseudotime trajectory analysis were used to explore the characteristics of macrophage subpopulations following IS. Key genes related to IS were identified using microarray data and validated for diagnostic potential through Receiver Operating Characteristic (ROC) analysis. Gene Set Enrichment Analysis (GSEA) was conducted to investigate the potential functions of these genes. Results The scRNA-seq data analysis revealed significant changes in macrophage subpopulation composition after IS. A specific macrophage subpopulation enriched in the stroke group was identified and designated as MCAO-specific macrophages (MSM). Pseudotime trajectory analysis indicated that MSM cells were in an intermediate stage of macrophage differentiation. Cell-cell communication analysis uncovered complex interactions between MSM cells and other cells, with the CCL6-CCR1 signaling axis potentially playing a crucial role in neuroinflammation. Two gene modules associated with MSM were identified via hdWGCNA, significantly enriched in pathways related to NOD-like receptors and antigen processing. By integrating differentially expressed MSM genes with conventional transcriptomic data, three IS-related hub genes were identified: Arg1, CLEC4D, and CLEC4E. Conclusion This study reveals the characteristics and functions of macrophage subpopulations following IS and identifies three hub genes with potential diagnostic value, providing novel insights into the pathological mechanisms of IS.
Macrophages/metabolism* ; Computational Biology/methods* ; Single-Cell Analysis/methods* ; Transcriptome ; Ischemic Stroke/metabolism* ; Animals ; Gene Regulatory Networks ; Gene Expression Profiling ; Humans ; Male

Objective Through integrative bioinformatics analysis of multi-source transcriptomic data, potential biomarkers to asthma epithelial cells were identified. The expression of these candidate target was subsequently validated in lung tissues and epithelial cells from asthma models. Methods The gene expression profile data of epithelial cells from three asthma patient cohorts and corresponding healthy controls were integrated from the Gene Expression Omnibus (GEO) database. Differential expression analysis and gene co-expression network analysis were performed to identify key genes and biological pathways associated with asthma. The key genes were validated in lung tissues and epithelial cells in asthma animal models. Results Differential gene expression analysis revealed 1121 upregulated and 1484 downregulated genes in epithelial cells from asthma patients compared with healthy controls. The biological pathway enrichment analysis revealed that the upregulated genes were mainly involved in glycosylation processes, whereas the downregulated genes were mainly associated with immune cell differentiation process. The gene co-expression network analysis revealed that module 9, enriched in glycosylation-related pathways, was significantly positively correlated with asthma, whereas module 17, associated with insulin and other signaling pathways, showed a significant negative correlation with asthma. We identified the genes of polypeptide N-acetylgalactosaminyltransferase 5 (GALNT5), pyrroline-5-carboxylate reductase 1 (PYCR1), and carcinoembryonic antigen-related cell adhesion molecule 5 (CEACAM5) as key genes within module 9, all of which were significantly upregulated in asthma. Finally, we validated that the expression levels of GALNT5, PYCR1, and CEACAM5 were significantly upregulated in epithelial cells from asthmatic lung tissue. Additionally, using a rat asthma model, we further confirmed that the protein levels of these three genes were significantly upregulated in lung tissues of the model group. Conclusion Through data integration and experimental validation, this study identified key genes and biological pathways closely associated with asthma pathogenesis. These findings provide a novel theoretical basis and potential targets for the diagnosis and treatment of asthma.
Asthma/metabolism* ; Humans ; Epithelial Cells/metabolism* ; Animals ; Biomarkers/metabolism* ; Gene Expression Profiling ; Transcriptome ; Gene Regulatory Networks ; Rats ; Computational Biology