Ophiopogon japonicus, a precious medicinal plant endemic to Zhejiang Province. Its tuberous roots are rich in bioactive components such as flavonoids, possessing anti-inflammatory, antioxidant, and immunomodulatory properties. To elucidate the impact of cadmium (Cd) stress on the accumulation and biosynthetic pathway of flavonoids in O. japonicus, this study exposed O. japonicus to different concentrations of Cd stress and explored the changes through integrated transcriptomics and metabolomics analysis. The results demonstrated that Cd stress (1 mg/L and 10 mg/L) significantly increased the content of flavonoids in O. japonicus in a concentration-dependent manner. The metabolomics analysis revealed a total of 110 flavonoids including flavones, flavanols, flavonols, flavone and flavonol derivatives, flavanones, isoflavonoids, chalcones and dihydrochalcones, and anthocyanins in O. japonicus, among which flavones, flavonols, flavone and flavonol derivatives, and anthocyanins increased under Cd stress. The transcriptomics analysis identified several key flavonoid biosynthesis-associated genes with up-regulated expression under Cd stress, including 14 genes encoding 4-coumarate CoA ligase (4CL), 2 genes encoding chalcone isomerase (CHI), and 14 genes encoding phenylalanine ammonia lyase (PAL). The gene-metabolite regulatory network indicated significant positive correlations of 4CL (Cluster-21637.5012, Cluster-21637.90648, and Cluster-21637.62637) and CHI (Cluster-21637.111909 and Cluster-21637.123300) with flavonoid metabolites, suggesting that these genes promoted the synthesis of specific flavonoid metabolites, which led to the accumulation of total flavonoids under Cd stress. These findings provide theoretical support for the cultivation and utilization of medicinal plants in Cd-contaminated environments and offered new perspectives for studying plant responses to heavy metal stress.
Cadmium/toxicity* ; Flavonoids/biosynthesis* ; Metabolomics ; Ophiopogon/drug effects* ; Stress, Physiological ; Transcriptome ; Gene Expression Profiling ; Gene Expression Regulation, Plant

The culture of suspended Vero cells is facing difficulties such as low cell viability and long doubling time. To investigate the main reasons for the slow growth and low viability of suspended Vero cells, this study conducted transcriptomic analysis of suspended Vero cells (Vero-XF) and adherent Vero cells (Vero-AD) to screen the differentially expressed genes (DEGs) affecting the growth of suspended cells. In addition, epidermal growth factor (EGF) was supplemented to the culture system to improve the growth of Vero-XF. The results showed that compared with the Vero-AD group, the Vero-XF group had 7 376 significant DEGs. Kyoto encyclopedia of genes and genomes enrichment analysis revealed that the DEGs were mainly enriched in the autophagy and mitophagy pathways. Eleven DEGs were selected and verified by quantitative real-time PCR, which showed up-regulated expression of ATG9B, WIPI2, LAMP2, OPTN, Rab7a, and DEPTOR and down-regulated expression of ATG4D, being consistent with the results of transcriptomic analysis. In addition, the Vero-XF group showed significantly up-regulated expression of ATG101, ATG2A, and STX17 and insignificant change in the expression of NBR1, compared with the Vero-AD group. The protein levels of LC3 and P62 in Vero-XF and Vero-AD were determined by Western blotting, which showed up-regulated expression of LC3Ⅱ/Ⅰ and down-regulated expression of P62 in Vero-XF, indicating a higher level of autophagy. Finally, the exogenous supplementation of EGF at 10, 20, and 30 μg/L in the culture system reduced the autophagy level of Vero-XF by 22.35%, 48.15%, and 71.29%, increased the specific growth rate by 15.48%, 33.33%, and 57.14%, and decreased the apoptosis rate by 2.84%, 15.46%, and 16.23%, respectively. The results of this study preliminarily reveal that the activation of autophagy is one of the reasons for the slow growth of Vero-XF, which provides reference for the subsequent culture of suspended Vero cells.
Animals ; Vero Cells ; Autophagy/genetics* ; Chlorocebus aethiops ; Epidermal Growth Factor/pharmacology* ; Gene Expression Profiling ; Transcriptome ; Cell Survival

Bacillus velezensis DJ1 exhibits broad-spectrum antagonistic activity against diverse phytopathogenic fungi, while its biocontrol mechanisms against Fusarium graminearum, the causal agent of maize stalk rot, remain poorly characterized. In this study, we integrated genomics and transcriptomics to elucidate the antifungal mechanisms of strain DJ1. The results demonstrated that DJ1 inhibited F. graminearum with the efficacy of 64.4%, while its polyketide crude extract achieved the control efficacy of 55% in pot experiments against this disease. Whole-genome sequencing revealed a single circular chromosome (3 929 792 bp, GC content of 47%) harboring 12 biosynthetic gene clusters for secondary metabolites, six of which encoded known antimicrobial compounds (macrolactin H, bacillaene, difficidin, surfactin, fengycin, and bacilysin). Transcriptomic analysis identified 243 differentially expressed genes (152 upregulated and 91 downregulated, P < 0.05), which were potentially associated with the antagonistic activity against F. graminearum. KEGG enrichment analysis highlighted activation (P < 0.05) of cysteine/methionine metabolism, pentose phosphate pathway, and polyketide biosynthesis pathways, indicating that DJ1 employed synergistic strategies involving antimicrobial compound synthesis, energy metabolism enhancement, and nutrient competition to suppress pathogens. This study provides a theoretical foundation for developing novel microbial resources and application technologies to combat phytopathogenic fungi.
Fusarium/drug effects* ; Bacillus/metabolism* ; Plant Diseases/prevention & control* ; Antifungal Agents/pharmacology* ; Genomics ; Zea mays/microbiology* ; Transcriptome ; Gene Expression Profiling ; Antibiosis ; Multigene Family ; Multiomics

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(AUC_train=0.891,AUC_test=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

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

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

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