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

OBJECTIVES: This study aimed to investigate the impact of foam macrophages (FMs) on the intracellular survival of Mycobacterium tuberculosis (MTB) and identify the molecular mechanisms influencing MTB survival. METHODS: An in vitro FM model was established using oleic acid induction. Transcriptomic and metabolomic analyses were conducted to identify the key molecular pathways involved in FM-mediated MTB survival. RESULTS: Induced FMs effectively restricted MTB survival. Transcriptomic and metabolomic profiling revealed distinct changes in gene and metabolite expression in FMs during MTB infection compared with normal macrophages. Integrated analyses identified significant alterations in the cyclic adenosine monophosphate (cAMP) signaling pathway, indicating that its activation contributes to the FM-mediated restriction of MTB survival. CONCLUSIONS FMs inhibit MTB survival. The cAMP signaling pathway is a key contributor. These findings enhance the understanding of the role of FMs in tuberculosis progression, suggest potential targets for host-directed therapies, and offer new directions for developing diagnostic and therapeutic strategies against tuberculosis.
Mycobacterium tuberculosis/physiology* ; Transcriptome ; Metabolomics ; Foam Cells/microbiology* ; Humans ; Metabolome ; Tuberculosis/microbiology* ; Gene Expression Profiling

OBJECTIVE: Patients with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection frequently develop central nervous system damage, yet the mechanisms driving this pathology remain unclear. This study investigated the primary pathways and key factors underlying brain tissue damage induced by the SARS-CoV-2 beta variant (lineage B.1.351). METHODS: K18-hACE2 and C57BL/6 mice were intranasally infected with the SARS-CoV-2 beta variant. Viral replication, pathological phenotypes, and brain transcriptomes were analyzed. Gene Ontology (GO) analysis was performed to identify altered pathways. Expression changes of host genes were verified using reverse transcription-quantitative polymerase chain reaction and Western blot. RESULTS: Pathological alterations were observed in the lungs of both mouse strains. However, only K18-hACE2 mice exhibited elevated viral RNA loads and infectious titers in the brain at 3 days post-infection, accompanied by neuropathological injury and weight loss. GO analysis of infected K18-hACE2 brain tissue revealed significant dysregulation of genes associated with innate immunity and antiviral defense responses, including type I interferons, pro-inflammatory cytokines, Toll-like receptor signaling components, and interferon-stimulated genes. Neuroinflammation was evident, alongside activation of apoptotic and pyroptotic pathways. Furthermore, altered neural cell marker expression suggested viral-induced neuroglial activation, resulting in caspase 4 and lipocalin 2 release and disruption of neuronal molecular networks. CONCLUSION These findings elucidate mechanisms of neuropathogenicity associated with the SARS-CoV-2 beta variant and highlight therapeutic targets to mitigate COVID-19-related neurological dysfunction.
Animals ; COVID-19/genetics* ; Mice ; Brain/metabolism* ; Apoptosis ; Mice, Inbred C57BL ; SARS-CoV-2/physiology* ; Pyroptosis ; Gene Expression Profiling ; Transcriptome ; Male ; Female

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