1.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
2.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
3.Conserved translational control in cardiac hypertrophy revealed by ribosome profiling.
Bao-Sen WANG ; Jian LYU ; Hong-Chao ZHAN ; Yu FANG ; Qiu-Xiao GUO ; Jun-Mei WANG ; Jia-Jie LI ; An-Qi XU ; Xiao MA ; Ning-Ning GUO ; Hong LI ; Zhi-Hua WANG
Acta Physiologica Sinica 2025;77(5):757-774
A primary hallmark of pathological cardiac hypertrophy is excess protein synthesis due to enhanced translational activity. However, regulatory mechanisms at the translational level under cardiac stress remain poorly understood. Here we examined the translational regulations in a mouse cardiac hypertrophy model induced by transaortic constriction (TAC) and explored the conservative networks versus the translatome pattern in human dilated cardiomyopathy (DCM). The results showed that the heart weight to body weight ratio was significantly elevated, and the ejection fraction and fractional shortening significantly decreased 8 weeks after TAC. Puromycin incorporation assay showed that TAC significantly increased protein synthesis rate in the left ventricle. RNA-seq revealed 1,632 differentially expressed genes showing functional enrichment in pathways including extracellular matrix remodeling, metabolic processes, and signaling cascades associated with pathological cardiomyocyte growth. When combined with ribosome profiling analysis, we revealed that translation efficiency (TE) of 1,495 genes was enhanced, while the TE of 933 genes was inhibited following TAC. In DCM patients, 1,354 genes were upregulated versus 1,213 genes were downregulated at the translation level. Although the majority of the genes were not shared between mouse and human, we identified 93 genes, including Nos3, Kcnj8, Adcy4, Itpr1, Fasn, Scd1, etc., with highly conserved translational regulations. These genes were remarkably associated with myocardial function, signal transduction, and energy metabolism, particularly related to cGMP-PKG signaling and fatty acid metabolism. Motif analysis revealed enriched regulatory elements in the 5' untranslated regions (5'UTRs) of transcripts with differential TE, which exhibited strong cross-species sequence conservation. Our study revealed novel regulatory mechanisms at the translational level in cardiac hypertrophy and identified conserved translation-sensitive targets with potential applications to treat cardiac hypertrophy and heart failure in the clinic.
Animals
;
Humans
;
Cardiomegaly/physiopathology*
;
Ribosomes/physiology*
;
Protein Biosynthesis/physiology*
;
Mice
;
Cardiomyopathy, Dilated/genetics*
;
Ribosome Profiling
4.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
5.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
6.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*
7.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
8.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
9.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
10.Single-cell transcriptomics combined with bioinformatics for comprehensive analysis of macrophage subpopulations and hub genes in ischemic stroke.
Jingyao XU ; Xiaolu WANG ; Shuai HOU ; Meng PANG ; Gang WANG ; Yanqiang WANG
Chinese Journal of Cellular and Molecular Immunology 2025;41(6):505-513
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

Result Analysis
Print
Save
E-mail