1.Construction of a Disulfidptosis-Related Prediction Model for Acute Myocardial Infarction Based on Transcriptome Data.
Qiu-Rong TANG ; Yang FENG ; Yao ZHAO ; Yun-Fei BIAN
Acta Academiae Medicinae Sinicae 2025;47(3):354-365
Objective To identify disulfidptosis-related gene(DRG)in acute myocardial infarction(AMI)by bioinformatics,analyze the molecular pattern of DRGs in AMI,and construct a DRGs-related prediction model.Methods AMI-related datasets were downloaded from the Gene Expression Omnibus database,and DRGs with differential expression were screened in AMI.CIBERSORT method was used to analyze the immune infiltration.Based on the differentially expressed DRGs,the AMI patients were classified into distinct subtypes via consensus clustering,followed by immune infiltration analysis,differential expression analysis,gene ontology and Kyoto encyclopedia of genes and genomes enrichment analysis,and gene set variation analysis.Weighted gene co-expression network analysis(WGCNA)was then performed to construct subtype-associated modules and identify hub genes.Finally,least absolute shrinkage and selection operator,random forest,and support vector machine-recursive feature elimination were used to screen feature genes to construct a DRGs-related prediction model.The model's diagnostic efficacy was evaluated by nomogram and receiver operating characteristic(ROC)curve analysis,followed by external validation.Results Nine differentially expressed DRGs were identified between AMI patients and controls.Based on the expression levels of these nine DRGs,AMI patients were divided into two DRGs subtypes,C1 and C2.Increased infiltration of monocytes,M0 macrophages,and neutrophils was observed in AMI patients and C1 subtype(all P<0.05),indicating a close correlation between DRGs and immune cells.There were 257 differentially expressed genes between the C1 and C2 subtypes,which were related to biological processes such as myeloid leukocyte activation and positive regulation of cytokines.Fcγ receptor-mediated phagocytosis and NOD-like receptor signaling pathway activity were enhanced in C1 subtype.WGCNA analysis suggested that the brown module exhibited the strongest correlation with DRG subtypes(r=0.67),from which 23 differentially expressed genes were identified.The feature genes screened by three machine learning methods were interpolated to obtain a DRGs-related prediction model consisting of three genes(AQP9,F5 and PYGL).Nomogram and ROC curves(AUCtrain=0.891,AUCtest=0.840)showed good diagnostic efficacy.Conclusions DRGs were closely related to the occurrence and progression of AMI.The DRGs-related prediction model consisting of AQP9,F5 and PYGL may provide targets for the diagnosis and personalized treatment of AMI.
Humans
;
Myocardial Infarction/diagnosis*
;
Transcriptome
;
Computational Biology
;
Gene Expression Profiling
;
ROC Curve
;
Gene Regulatory Networks
;
Nomograms
;
Disulfidptosis
2.Single-cell transcriptomic analysis reveals immune dysregula-tion and macrophage reprogramming in diabetic foot ulcers.
Chunli HUANG ; Yu JIANG ; Wei JIAO ; Ying SUI ; Chunlei WANG ; Yongtao SU
Journal of Zhejiang University. Medical sciences 2025;54(5):602-610
OBJECTIVES:
To elucidate the underlying mechanisms of macrophage-mediated inflammation and tissue injury in diabetic foot ulcer (DFU).
METHODS:
Skin tissue samples were collected from patients with DFU and with non-DFU. A total of 79 272 high-quality cell transcriptomes were obtained using single-cell RNA sequencing. An unbiased clustering approach was employed to identify cell subpopulations. Seurat functions were used to identify differentially expressed genes between DFU and non-DFU groups, and gene ontology (GO) enrichment analysis was used to reveal gene function. Furthermore, cell-cell communication network construction and ligand-receptor interaction analysis were performed to reveal the mechanisms underlying cellular interactions and signaling regulation in the DFU microenvironment from multiple perspectives.
RESULTS:
The results revealed a significant expansion of myeloid cells in DFU tissues, alongside a marked reduction in structural cells such as endothelial cells, epithelial cells, and smooth muscle cells. Major cell types underwent functional reprogramming, characterized by immune activation and impaired tissue remodeling. Specifically, macrophages in DFU skin tissues exhibited a shift toward a pro-inflammatory M1 phenotype, with upregulation of genes associated with inflammation and oxidative stress. Cell communication analysis further demonstrated that M1 macrophages served as both primary signal receivers and influencers in the COMPLEMENT pathway mediated communication network, and as key signal senders and mediators in the secreted phosphoprotein 1 (SPP1) pathway mediated communication network, actively shaping the inflammatory microenvironment. Key ligand-receptor interactions driving macrophage signaling were identified, including C3-(ITGAM+ITGB2) and SPP1-CD44.
CONCLUSIONS
This study establishes a comprehensive single-cell atlas of DFU, revealing the role of macrophage-driven cellular networks in chronic inflammation and impaired healing. These findings may offer potential novel therapeutic targets for DFU treatment.
Humans
;
Macrophages/immunology*
;
Diabetic Foot/pathology*
;
Single-Cell Analysis
;
Transcriptome
;
Gene Expression Profiling
;
Inflammation
;
Skin
;
Cell Communication
;
Signal Transduction
;
Cellular Reprogramming
3.Hub biomarkers and their clinical relevance in glycometabolic disorders: A comprehensive bioinformatics and machine learning approach.
Liping XIANG ; Bing ZHOU ; Yunchen LUO ; Hanqi BI ; Yan LU ; Jian ZHOU
Chinese Medical Journal 2025;138(16):2016-2027
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
4.Transcriptomic analysis of key genes involved in sex differences in intellectual development.
Jia-Wei ZHANG ; Xiao-Li ZHENG ; Hai-Qian ZHOU ; Zhen ZHU ; Wei HAN ; Dong-Min YIN
Acta Physiologica Sinica 2025;77(2):211-221
Intelligence encompasses various abilities, including logical reasoning, comprehension, self-awareness, learning, planning, creativity, and problem-solving. Extensive research and practical experience suggest that there are sex differences in intellectual development, with females typically maturing earlier than males. However, the key genes and molecular network mechanisms underlying these sex differences in intellectual development remain unclear. To date, Genome-Wide Association Studies (GWAS) have identified 507 genes that are significantly associated with intelligence. This study first analyzed RNA sequencing data from different stages of brain development (from BrainSpan), revealing that during the late embryonic stage, the average expression levels of intelligence-related genes are higher in males than in females, while the opposite is observed during puberty. This study further constructed interaction networks of intelligence-related genes with sex-differential expression in the brain, including the prenatal male network (HELP-M: intelligence genes with higher expression levels in prenatal males) and the pubertal female network (HELP-F: intelligence genes with higher expression levels in pubertal females). The findings indicate that the key genes in both networks are Ep300 and Ctnnb1. Specifically, Ep300 regulates the transcription of 53 genes in both HELP-M and HELP-F, while Ctnnb1 regulates the transcription of 45 genes. Ctnnb1 plays a more prominent role in HELP-M, while Ep300 is more crucial in HELP-F. Finally, this study conducted sequencing validation on rats at different developmental stages, and the results indicated that in the prefrontal cortex of female rats during adolescence, the expression levels of the intelligence genes in HELP-F, as well as key genes Ep300 and Ctnnb1, were higher than those in male rats. These genes were also involved in neurodevelopment-related biological processes. The findings reveal a sex-differentiated intelligence gene network and its key genes, which exhibit varying expression levels during the neurodevelopmental process.
Female
;
Intelligence/physiology*
;
Male
;
Sex Characteristics
;
Animals
;
Brain/growth & development*
;
E1A-Associated p300 Protein/physiology*
;
beta Catenin/physiology*
;
Transcriptome
;
Rats
;
Gene Expression Profiling
;
Genome-Wide Association Study
5.Transcriptome sequencing reveals molecular mechanism of seed dormancy release of Zanthoxylum nitidum.
Chang-Qian QUAN ; Dan-Feng TANG ; Jian-Ping JIANG ; Yan-Xia ZHU
China Journal of Chinese Materia Medica 2025;50(1):102-110
The transcriptome sequencing based on Illumina Novaseq 6000 Platform was performed with the untreated seed embryo(DS), stratified seed embryo(SS), and germinated seed embryo(GS) of Zanthoxylum nitidum, aiming to explore the molecular mechanism regulating the seed dormancy and germination of Z. nitidum and uncover key differentially expressed genes(DEGs). A total of 61.41 Gb clean data was obtained, and 86 386 unigenes with an average length of 773.49 bp were assembled. A total of 29 290 DEGs were screened from three comparison groups(SS vs DS, GS vs SS, and GS vs DS), and these genes were annotated on 134 Kyoto Encyclopedia of Genes and Genomes(KEGG) pathways. KEGG enrichment analysis revealed that the plant hormone signal transduction pathway is the richest pathway, containing 226 DEGs. Among all DEGs, 894 transcription factors were identified, which were distributed across 34 transcription factor families. These transcription factors were also mainly concentrated in plant hormone signal transduction and mitogen-activated protein kinase(MAPK) signaling pathways. Further real-time quantitative polymerase chain reaction(RT-qPCR) validation of 12 DEGs showed that the transcriptome data is reliable. During the process of seed dormancy release and germination, a large number of DEGs involved in polysaccharide degradation, protein synthesis, lipid metabolism, and hormone signal transduction were expressed. These genes were involved in multiple metabolic pathways, forming a complex regulatory network for dormancy and germination. This study lays a solid foundation for analyzing the molecular mechanisms of seed dormancy and germination of Z. nitidum.
Zanthoxylum/metabolism*
;
Plant Dormancy/genetics*
;
Seeds/metabolism*
;
Gene Expression Regulation, Plant
;
Plant Proteins/metabolism*
;
Transcriptome
;
Gene Expression Profiling
;
Germination
;
Transcription Factors/metabolism*
;
Plant Growth Regulators/genetics*
;
Signal Transduction
6.Selection and validation of reference genes for quantitative real-time PCR analysis in Tujia medicine Xuetong.
Qian XIAO ; Chen-Si TAN ; Jiang ZENG ; Yuan-Shu XU ; Tian-Hao FU ; Lu-Yun NING ; Wei WANG
China Journal of Chinese Materia Medica 2025;50(3):682-692
Tujia ethnic group medicine Xuetong is derived from Kadsura heteroclita, the stem of which has the medicinal value for anti-rheumatoid arthritis, liver protection, anti-tumor, anti-oxidation effects, and has been widely used in Hunan and Guangdong in China. The selection of reliable and stable reference genes is the basis for subsequent molecular research on K. heteroclita. In this study, GAPDH, TUA, Actin, UBQ, EF-1α, 18S-rRNA, CYP, UBC, TUB, H2A, and RPL were selected as candidate reference genes in Kadsura heteroclita. The gene expression levels of the 11 candidate reference genes of K. heteroclita in its 6 different parts(stem-inside of the cambium, stem-outside of the cambium, fruit, flower, root, and leaf) and under different intervention conditions [drought stress, salt stress, and methyl jasmonate(MeJA) treatment] were detected by quantitative real-time polymerase chain reaction(qRT-PCR). The expression stability of the 11 candidate reference genes was comprehensively analyzed and evaluated by geNorm, NormFinder, ΔCT algorithm, and RefFinder software. The results showed that the expression of UBC and RPL was relatively stable in 6 different parts, and UBC and GAPDH genes were relatively stable under different intervention conditions. To verify the reliability of reference genes for K. heteroclita, this study further examined the relative expression levels of KhFPS, KhIDI, KhCAS, KhSQE, KhSQS, KhSQS-2, KhHMGS, KhHMGR, KhMVD, KhMVK, KhDXR, KhDXS, KhPMVK, and KhGGPS in different parts and under different intervention conditions, which might relate to the synthesis of the main component(Xuetongsu) of K. heteroclita. The results showed that with UBC and RPL or UBC and GAPDH as the reference genes, the expression trends of these 14 genes were basically consistent in different parts or under different intervention conditions for K. heteroclita. In conclusion, UBC can be used as a reference gene of K. heteroclita for its different parts and different intervention conditions, which lays a foundation for further research on the biosynthetic pathway of main components in K. heteroclita.
Real-Time Polymerase Chain Reaction/methods*
;
Reference Standards
;
Gene Expression Regulation, Plant
;
Gene Expression Profiling
;
Plant Proteins/metabolism*
;
Drugs, Chinese Herbal
7.Transcriptome analysis and catechin synthesis genes in different organs of Spatholobus suberectus.
Wei-Qi QIN ; Quan LIN ; Ying LIANG ; Fan WEI ; Gui-Li WEI ; Qi GAO ; Shuang-Shuang QIN
China Journal of Chinese Materia Medica 2025;50(12):3297-3306
To study the differences in transcript levels among different organs of Spatholobus suberectus and to explore the genes encoding enzymes related to the catechin biosynthesis pathway, this study utilized the genome and full-length transcriptome data of S. suberectus as references. Transcriptome sequencing and bioinformatics analysis were performed on five different organs of S. suberectus-roots, stems, leaves, flowers, and fruits-using the Illumina NovaSeq 6000 platform. A total of 115.28 Gb of clean data were obtained, with GC content values ranging from 45.19% to 47.54%, Q20 bases at 94.17% and above, and an overall comparison rate with the reference genome around 90%. In comparisons between the stem and root, stem and leaf, stem and flower, and stem and fruit, 10 666, 9 674, 9 320, and 5 896 differentially expressed genes(DEGs) were identified, respectively. The lowest number of DEGs was found in the stem and root comparison group. KEGG enrichment analysis revealed that the DEGs were mainly concentrated in the pathways of phytohormone signaling, phenylalanine biosynthesis, etc. A total of 39 genes were annotated in the catechin biosynthesis pathway, with at least one highly expressed gene found in all organs. Among these, PAL1, PAL2, C4H1, C4H3, 4CL1, 4CL2, and DFR2 showed high expression in the stems, suggesting that they may play important roles in the biosynthesis of flavonoids in S. suberectus. This study aims to provide important information for the in-depth exploration of the regulation of catechin biosynthesis in S. suberectus through transcriptome analysis of its different organs and to provide a reference for the further realization of S. suberectus varietal improvement and molecular breeding.
Catechin/biosynthesis*
;
Gene Expression Profiling
;
Gene Expression Regulation, Plant
;
Plant Proteins/metabolism*
;
Fabaceae/metabolism*
;
Transcriptome
;
Flowers/metabolism*
;
Plant Stems/metabolism*
;
Plant Leaves/metabolism*
;
Plant Roots/metabolism*
;
Fruit/metabolism*
8.Identification and expression analysis of AP2/ERF family members in Lonicera macranthoides.
Si-Min ZHOU ; Mei-Ling QU ; Juan ZENG ; Jia-Wei HE ; Jing-Yu ZHANG ; Zhi-Hui WANG ; Qiao-Zhen TONG ; Ri-Bao ZHOU ; Xiang-Dan LIU
China Journal of Chinese Materia Medica 2025;50(15):4248-4262
The AP2/ERF transcription factor family is a class of transcription factors widely present in plants, playing a crucial role in regulating flowering, flower development, flower opening, and flower senescence. Based on transcriptome data from flower, leaf, and stem samples of two Lonicera macranthoides varieties, 117 L. macranthoides AP2/ERF family members were identified, including 14 AP2 subfamily members, 61 ERF subfamily members, 40 DREB subfamily members, and 2 RAV subfamily members. Bioinformatics and differential gene expression analyses were performed using NCBI, ExPASy, SOMPA, and other platforms, and the expression patterns of L. macranthoides AP2/ERF transcription factors were validated via qRT-PCR. The results indicated that the 117 LmAP2/ERF members exhibited both similarities and variations in protein physicochemical properties, AP2 domains, family evolution, and protein functions. Differential gene expression analysis revealed that AP2/ERF transcription factors were primarily differentially expressed in the flowers of the two L. macranthoides varieties, with the differentially expressed genes mainly belonging to the ERF and DREB subfamilies. Further analysis identified three AP2 subfamily genes and two ERF subfamily genes as potential regulators of flower development, two ERF subfamily genes involved in flower opening, and two ERF subfamily genes along with one DREB subfamily gene involved in flower senescence. Based on family evolution and expression analyses, it is speculated that AP2/ERF transcription factors can regulate flower development, opening, and senescence in L. macranthoides, with ERF subfamily genes potentially serving as key regulators of flowering duration. These findings provide a theoretical foundation for further research into the specific functions of the AP2/ERF transcription factor family in L. macranthoides and offer important theoretical insights into the molecular mechanisms underlying floral phenotypic differences among its varieties.
Plant Proteins/chemistry*
;
Gene Expression Regulation, Plant
;
Transcription Factors/chemistry*
;
Lonicera/classification*
;
Flowers/metabolism*
;
Phylogeny
;
Gene Expression Profiling
;
Multigene Family
9.A novel glycolysis-related prognostic risk model for colorectal cancer patients based on single-cell and bulk transcriptomic data.
Kai YAO ; Jingyi XIA ; Shuo ZHANG ; Yun SUN ; Junjie MA ; Bo ZHU ; Li REN ; Congli ZHANG
Chinese Journal of Cellular and Molecular Immunology 2025;41(2):105-115
Objective To explore the prognostic value of glycolysis-related genes in colorectal cancer (CRC) patients and formulate a novel glycolysis-related prognostic risk model. Methods Single-cell and bulk transcriptomic data of CRC patients, along with clinical information, were obtained from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. Glycolysis scores for each sample were calculated using single-sample Gene Set Enrichment Analysis (ssGSEA). Kaplan-Meier survival curves were generated to analyze the relationship between glycolysis scores and overall survival. Novel glycolysis-related subgroups were defined among the cell type with the highest glycolysis scores. Gene enrichment analysis, metabolic activity assessment, and univariate Cox regression were performed to explore the biological functions and prognostic impact of these subgroups. A prognostic risk model was built and validated based on genes significantly affecting the prognosis. Gene Set Enrichment Analysis (GSEA) was conducted to explore differences in biological processes between high- and low-risk groups. Differences in immune microenvironment and drug sensitivity between these groups were assessed using R packages. Potential targeted agents for prognostic risk genes were predicted using the Enrichr database. Results Tumor tissues showed significantly higher glycolysis scores than normal tissues, which was associated with a poor prognosis in CRC patients. The highest glycolysis score was observed in epithelial cells, within which we defined eight novel glycolysis-related cell subpopulations. Specifically, the P4HA1+ epithelial cell subpopulation was associated with a poor prognosis. Based on signature genes of this subpopulation, a six-gene prognostic risk model was formulated. GSEA revealed significant biological differences between high- and low-risk groups. Immune microenvironment analysis demonstrated that the high-risk group had increased infiltration of macrophages and tumor-associated fibroblasts, along with evident immune exclusion and suppression, while the low-risk group exhibited higher levels of B cell and T cell infiltration. Drug sensitivity analysis indicated that high-risk patients were more sensitive to Abiraterone, while low-risk patients responded to Cisplatin. Additionally, Valproic acid was predicted as a potential targeted agent. Conclusion High glycolytic activity is associated with a poor prognosis in CRC patients. The novel glycolysis-related prognostic risk model formulated in this study offers significant potential for enhancing the diagnosis and treatment of CRC.
Humans
;
Colorectal Neoplasms/pathology*
;
Glycolysis/genetics*
;
Prognosis
;
Transcriptome
;
Tumor Microenvironment/genetics*
;
Gene Expression Profiling
;
Single-Cell Analysis
;
Gene Expression Regulation, Neoplastic
;
Male
;
Female
;
Kaplan-Meier Estimate
10.Unveiling the molecular features and diagnosis and treatment prospects of immunothrombosis via integrated bioinformatics analysis.
Yafen WANG ; Xiaoshuang WU ; Zhixin LIU ; Xinlei LI ; Yaozhen CHEN ; Ning AN ; Xingbin HU
Chinese Journal of Cellular and Molecular Immunology 2025;41(3):228-235
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

Result Analysis
Print
Save
E-mail