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

OBJECTIVE: To screen novel diagnostic marker or therapeutic target for multiple myeloma (MM). METHODS: Sel1L, SPAG4, KCNN3 and PARM1 were identified by bioinformatics method based on GEO database as high expression genes in MM. Their RNA and protein expression levels in bone marrow mononuclear cells from myeloma cell lines U266, NCI-H929, MM.1s, RPMI8226 and leukemia cell line THP1, as well as 31 MM patients were evaluated by RT-PCR and Western blot, respectively. Meanwhile, 5 samples of bone marrow from healthy donors for allogeneic hematopoietic stem cell transplantation were employed as controls. RESULTS: Compared with leukemia cell line THP1, the expression levels of KCNN3, PARM1 and Sel1L mRNA were significantly increased in myeloma cell lines U266, NCI-H929 and MM.1s, while PARM1 was further increased in myeloma cell lines 8226. Western blot showed that the 4 genes were all expressed in the 4 myeloma cell lines. Compared with healthy controls, the expression levels of Sel1L, SPAG4, KCNN3 and PARM1 mRNA were significantly higher in MM patients (all P < 0.05). Western blot showed that the 4 genes were all expressed in MM patients, and the protein expression level of Sel1L and KCNN3 were significantly different compared with healthy donors (all P < 0.01). CONCLUSION Sel1L, SPAG4, KCNN3 and PARM1 may be potential diagnostic markers and therapeutic targets for MM.
Humans ; Multiple Myeloma/genetics* ; Cell Line, Tumor ; Proteins/metabolism* ; Computational Biology ; RNA, Messenger/genetics*

OBJECTIVE: The aim of this study is to investigate the main active substances of Qilin pills by high performance liquid chromatogre-electrostatic field orbitrap mass spectrometry (HPLC-Q-Orbitrap /MS), and explore the mechanism of its action in the treatment of asthenozoospermia by combining network pharmacology and molecular docking. METHODS: (1) Qilin pills were quantitatively and qualitatively analyzed by HPLC-Q-Orbitrap /MS. (2) The top 100 compounds in Qilin pills were screened by content analysis and SwissADME, and their targets were predicted. The asthenozoospermia targets were searched through the database. And a "protein-protein interaction" (PPI) network was constructed. KEGG and GO analysis was performed using the DAVID database. And a "drug-target-pathway" network was constructed. (3) SailVina was used for molecular docking. RESULTS: (1) A total of 1 275 known components were found and ranked in Qilin pills by HPLC-Q-Orbitrap /MS analysis. (2) The top 100 compounds in Qilin pills predicted a total of 1 053 targets and 184 potential therapeutic targets for asthenozoospermia. KEGG pathway analysis and GO analysis showed that the treatment of asthenozoospermia by Qilin pills may be related to the steroid hormone synthesis pathway, the response to steroid hormones, the chromosomal region of cells and the activity of steroid hydroxylase. The mechanism of Qilin pills in treating asthenozoospermia may be related to regulating the synthesis, metabolism and reaction process of sex hormone in the body. (3) The molecular docking results of its key targets (CYP19A1, ESR1, HSP90AA1, p53, HIF1α and BCL2) showed that the key active ingredients M030, M039, M043, M050, M055 and M073 of Qilin pills had spontaneous binding. It had a binding energy of less than －5 kJ /mol. CONCLUSION The material basis of Qilin pills has been explored by this study. And the mechanism of action of Qilin pills in the treatment of asthenozoospermia is highly bound to the expression and response process of steroid hormones, which provides a theoretical basis for the clinical application of Qilin pills.
Asthenozoospermia/drug therapy* ; Chromatography, High Pressure Liquid ; Molecular Docking Simulation ; Drugs, Chinese Herbal/chemistry* ; Male ; Computational Biology ; Humans ; Mass Spectrometry ; Protein Interaction Maps ; Liquid Chromatography-Mass Spectrometry

OBJECTIVES: Hepatocellular carcinoma is one of the most common primary malignant tumors with the third highest mortality rate worldwide. This study aims to identify key genes associated with hepatocellular carcinoma prognosis using the Gene Expression Omnibus (GEO) database and provide a theoretical basis for discovering novel prognostic biomarkers for hepatocellular carcinoma. METHODS: Hepatocellular carcinoma-related datasets were retrieved from the GEO database. Differentially expressed genes (DEGs) were identified using the GEO2R tool. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed using the Database for Annotation, Visualization, and Integrated Discovery (DAVID). A protein-protein interaction (PPI) network was constructed using the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING), and key genes were identified using Cytoscape software. The University of Alabama at Birmingham Cancer Data Analysis Resource (UALCAN) was used to analyze the expression levels of key genes in normal and hepatocellular carcinoma tissues, as well as their associations with pathological grade, clinical stage, and patient survival. The Human Protein Atlas (THPA) was used to further validate the impact of key genes on overall survival. Expression levels of key genes in the blood of hepatocellular carcinoma patients were evaluated using the expression atlas of blood-based biomarkers in the early diagnosis of cancers (BBCancer). RESULTS: A total of 78 DEGs were identified from the GEO database. GO and KEGG analyses indicated that these genes may contribute to hepatocellular carcinoma progression by promoting cell division and regulating protein kinase activity. Sixteen key genes were screened via Cytoscape and validated using UALCAN and THPA. These genes were overexpressed in hepatocellular carcinoma tissues and were associated with disease progression and poor prognosis. Finally, BBCancer analysis showed that ASPM and NCAPG were also elevated in the blood of hepatocellular carcinoma patients. CONCLUSIONS This study identified 16 key genes as potential prognostic biomarkers for hepatocellular carcinoma, among which ASPM and NCAPG may serve as promising blood-based markers for hepatocellular carcinoma.
Humans ; Carcinoma, Hepatocellular/mortality* ; Liver Neoplasms/pathology* ; Prognosis ; Computational Biology/methods* ; Protein Interaction Maps/genetics* ; Biomarkers, Tumor/genetics* ; Gene Expression Regulation, Neoplastic ; Gene Expression Profiling ; Gene Ontology ; Databases, Genetic

OBJECTIVES: Osteoarthritis (OA) and sarcopenia are significant health concerns in the elderly, substantially impacting their daily activities and quality of life. However, the relationship between them remains poorly understood. This study aims to uncover common biomarkers and pathways associated with both OA and sarcopenia. METHODS: Gene expression profiles related to OA and sarcopenia were retrieved from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) between disease and control groups were identified using R software. Common DEGs were extracted via Venn diagram analysis. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were conducted to identify biological processes and pathways associated with shared DEGs. Protein-protein interaction (PPI) networks were constructed, and candidate hub genes were ranked using the maximal clique centrality (MCC) algorithm. Further validation of hub gene expression was performed using 2 independent datasets. Receiver operating characteristic (ROC) curve analysis was used to evaluate the predictive value of key genes for OA and sarcopenia. Mouse models of OA and sarcopenia were established. Hematoxylin-eosin and Safranin O/Fast Green staining were used to validate the OA model. The sarcopenia model was validated via rotarod testing and quadriceps muscle mass measurement. Real-time reverse transcription PCR (real-time RT-PCR) was employed to assess the mRNA expression levels of candidate key genes in both models. Gene set enrichment analysis (GSEA) was conducted to identify pathways associated with the selected shared key genes in both diseases. RESULTS: A total of 89 common DEGs were identified in the gene expression profiles of OA and sarcopenia, including 76 upregulated and 13 downregulated genes. These 89 DEGs were significantly enriched in protein digestion and absorption, the PI3K-Akt signaling pathway, and extracellular matrix-receptor interaction. PPI network analysis and MCC algorithm analysis of the 89 common DEGs identified the top 17 candidate hub genes. Based on the differential expression analysis of these 17 candidate hub genes in the validation datasets, AEBP1 and COL8A2 were ultimately selected as the common key genes for both diseases, both of which showed a significant upregulation trend in the disease groups (all P<0.05). The value of area under the curve (AUC) for AEBP1 and COL8A2 in the OA and sarcopenia datasets were all greater than 0.7, indicating that both genes have potential value in predicting OA and sarcopenia. Real-time RT-PCR results showed that the mRNA expression levels of AEBP1 and COL8A2 were significantly upregulated in the disease groups (all P<0.05), consistent with the results observed in the bioinformatics analysis. GSEA revealed that AEBP1 and COL8A2 were closely related to extracellular matrix-receptor interaction, ribosome, and oxidative phosphorylation in OA and sarcopenia. CONCLUSIONS AEBP1 and COL8A2 have the potential to serve as common biomarkers for OA and sarcopenia. The extracellular matrix-receptor interaction pathway may represent a potential target for the prevention and treatment of both OA and sarcopenia.
Sarcopenia/genetics* ; Osteoarthritis/genetics* ; Computational Biology/methods* ; Humans ; Protein Interaction Maps/genetics* ; Animals ; Mice ; Gene Expression Profiling ; Gene Ontology ; Transcriptome ; Male ; Signal Transduction/genetics* ; Gene Regulatory Networks

OBJECTIVES: Polycystic ovary syndrome (PCOS) is a common endocrine disorder that affects women's health. This study aims to investigate gene and transcription factor (TF) expression differences between PCOS patients and healthy individuals using bioinformatics approaches, and to verify the function of key transcription factors, with the goal of providing new insights into the pathogenesis of PCOS. METHODS: Differentially expressed genes (DEGs) and differentially expressed transcription factors (DETFs) between PCOS patients and controls were identified from the RNA sequencing dataset GSE168404 using bioinformatics methods. Functional enrichment analysis was performed using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. The expression and function of core transcription factors were further validated in ovarian tissues of PCOS model mice and control mice using Western blotting and reverse transcription quantitative polymerase chain reaction (RT-qPCR). RESULTS: A total of 332 DEGs were identified between PCOS patients and controls, including 259 upregulated and 73 downregulated genes in the PCOS group. 19 DETFs were further screened, of which 16 were upregulated and 3 were downregulated in PCOS. The upregulated DETFs (including TFCP2L1, DACH1, ESR2, AFF3, SMAD9, ZNF331, HOPX,ATOH8, HIF3α, DPF3, HOXC4, HES1, ID1, JDP2, SOX4, and ID3) were primarily associated with lipid metabolism, development, and cell adhesion. Protein and mRNA expression analysis in PCOS model mice revealed significantly decreased levels of hypoxia-inducible factor (HIF) 1α and HIF2α, and significantly increased expression of HIF3α compared to control mice (all P<0.001). CONCLUSIONS Significant differences in gene and TF expression exist between PCOS patients and healthy individuals. HIF-3α may play a crucial role in PCOS and could serve as a novel biomarker for diagnosis and a potential therapeutic target.
Polycystic Ovary Syndrome/metabolism* ; Female ; Humans ; Animals ; Mice ; Transcription Factors/metabolism* ; Computational Biology ; Gene Expression Profiling ; Adult