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

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

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

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

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*

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

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

OBJECTIVES: Multiple myeloma (MM) is a highly heterogeneous hematologic malignancy, with disease progression driven by cytogenetic abnormalities and a complex bone marrow microenvironment. This study aims to construct a prognostic model for MM based on transcriptomic data and lipid metabolism related genes (LRGs), and to identify potential drug targets for high-risk patients to support clinical decision-making. METHODS: In this study, 2 transcriptomic datasets covering 985 newly diagnosed MM patients were retrieved from the Gene Expression Omnibus (GEO) database. Univariate Cox regression and 101 machine learning algorithms were used for gene selection. An LRG-based prognostic model was constructed using Stepwise Cox (both directions) and random survival forest (RSF) algorithms. The association between the prognostic score and clinical events was evaluated, and model performance was assessed using time-dependent receiver operating characteristic (ROC) curves and the C-index. The added predictive value of combining prognostic scores with clinical variables and staging systems was also analyzed. Differentially expressed genes between high- and low-risk groups were identified using limma and clusterProfiler and subjected to pathway enrichment analysis. Drug sensitivity analysis was conducted using the Genomics of Drug Sensitivity in Cancer (GDSC) database and oncoPredict to identify potential therapeutic targets for high-risk patients. The functional role of key LRGs in the model was validated via in vitro cell experiments. RESULTS: An LRG-based prognostic model (LRG17) was successfully developed using transcriptomic data and machine learning. The model demonstrated robust predictive performance, with area under the curve (AUC) values of 0.962, 0.912, and 0.842 for 3-, 5-, and 7-year survival, respectively. Patients were stratified into high- and low-risk groups, with high-risk patients showing significantly shorter overall survival (OS) and event-free survival (EFS) (both P<0.001) and worse clinical profiles (e.g., lower albumin, higher β2-microglobulin and lactate dehydrogenase levels). Enrichment analysis revealed that high-risk patients were significantly enriched for pathways related to chromosome segregation and mitosis, whereas low-risk patients were enriched for immune response and immune cell activation pathways. Drug screening suggested that AURKA inhibitor BMS-754807 and FGFR3 inhibitor I-BET-762 may be more effective in high-risk patients. Functional assays demonstrated that silencing of key LRG PLA2G4A significantly inhibited cell viability and induced apoptosis. CONCLUSIONS LRGs serve as promising biomarkers for prognosis prediction and risk stratification in MM. The overexpression of chromosomal instability-related and high-risk genetic event-associated genes in high-risk patients may explain their poorer outcomes. Given the observed resistance to bortezomib and lenalidomide in high-risk patients, combination therapies involving BMS-754807 or I-BET-762 may represent effective alternatives.
Humans ; Multiple Myeloma/mortality* ; Prognosis ; Lipid Metabolism/genetics* ; Transcriptome ; Machine Learning ; Male ; Female ; Gene Expression Profiling ; Algorithms