SEARCH · Search PubMed

Results for “co-expression”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6Linked to original sources

PubMed2025

Liver cancer-specific prognostic model developed using endoplasmic reticulum stress-related LncRNAs and LINC01011 as a potential therapeutic target.

Liver cancer is a serious malignancy worldwide, and long noncoding RNAs (lncRNAs) have been implicated in its prognosis.It remains unclear how lncRNAs related to endoplasmic reticulum stress (ERS) influence liver cancer prognosis. Here, we analyzed RNA and clinical data from the Cancer Genome Atlas and sourced ERS-related genes from the Molecular Signatures Database. Co-expression analysis identified ERS-related lncRNAs, and Cox regression analysis as well as least absolute shrinkage and selection operator regression highlighted three lncRNAs for a prognostic model. Based on median risk scores, we classified patients into two risk groups. The high-risk group displayed poor prognosis, and this finding was validated in the test set. According to consistency clustering, the patients were assigned to two clusters, and tumor microenvironment scores were computed. Patients with a high mutation burden had worse outcomes. Furthermore, immune infiltration analysis indicated more immune cells and mutations in checkpoint molecules among high-risk individuals. Drug sensitivity varied between the risk groups. LINC01011 was selected for functional assays. Colony formation assay and CCK-8 assay revealed that silencing LINC01011 suppressed liver cancer cell proliferation. Transwell and scratch assays indicated that silencing LINC01011 inhibited liver cancer cell migration. Western blotting assay revealed that inhibiting LINC01011 induced apoptosis and simultaneously inhibited epithelial-mesenchymal transition. These findings confirm the validity of the prognostic model and indicate that LINC01011 could serve as a potential research target.

Humans↗

PubMed2025

Transcriptomics-based exploration of ubiquitination-related biomarkers and potential molecular mechanisms in laryngeal squamous cell carcinoma.

BACKGROUND: One of the most common and prevalent cancers is laryngeal squamous cell carcinoma (LSCC), which poses a great threat to the life and health of the patient. Nonetheless, it has been demonstrated that ubiquitination is crucial for the development and course of LSCC. Therefore, it is particularly important to identify biomarkers for ubiquitination-related genes (UbRGs) in LSCC. METHODS: Differentially expressed genes (DEGs) in the LSCC versus controls were obtained by differential expression analysis. Also, key modular genes associated with LSCC were obtained using weighted gene co-expression network analysis (WGCNA). Next, DEGs, key module genes, and UbRGs were taken to intersect to obtain candidate genes. And then machine algorithms were to screen potential biomarkers, further their diagnostic value were analyzed and validated. Then, therapeutic agents for biomarkers were predict. In addition, the regulatory networks of the biomarkers were mapped. The expression levels of biomarkers were detected in clinical samples using reverse transcription-quantitative PCR (RT-qPCR). RESULTS: A total of eight candidate genes were acquired by the overlap 1,911 DEGs, the key modular genes of WGCNA, and 1,393 UbRGs. A sum of four biomarkers (WDR54, KAT2B, NBEAL2 and LNX1) were identified by two machine learning, then these four biomarkers were validated in GSE127165 and the expression trend was consistent with TCGA-LSCC, they were recorded as biomarkers. Moreover, the accuracy of the biomarkers in predicting clinical aspects of LSCC was confirmed by the receiver operating characteristic (ROC) curves. Subsequently, cancers such as malignant neoplasms, colorectal cancers, tumors, and primary malignant neoplasms were significantly associated with the biomarkers, which further suggests that these four biomarkers were strongly associated with cancer. Meanwhile, the drugs garcinol, cocaine, and triazolam, among others, used for LSCC treatment were predicted. Finally, transcription factors (TFs) (BRD4, MYC, AR, and CTCF) were predicted to regulate the biomarkers. RT-qPCR assays illustrated that the expression trends of KAT2B, LNX1 and NBEAL2 remained consistent with the dataset. CONCLUSION: The identification of four biomarkers (WDR54, KAT2B, NBEAL2 and LNX1) associated with UbRGs could ultimately serve as a predictive clinical diagnosis of LSCC and provide insight into the molecular mechanisms of LSCC.

Humans↗

PubMed2025

Identification of potential biomarkers and mechanisms for keloid disorder based on comprehensive bioinformatics analysis and machine learning algorithms.

BACKGROUND: Keloid disorder (KD) encompasses a spectrum of fibroproliferative dermal conditions, the pathogenesis remains complex and incompletely understood. This study sought to identify biomarkers and potential therapeutic targets for KD through an integrative bioinformatics approach and machine learning analysis of RNA sequencing data. METHODS: RNA sequencing was performed on skin tissue samples from 13 patients with KD and 14 healthy controls. Using weighted gene co-expression network analysis and differential expression analysis revealed differentially expressed key module genes, and the CytoHubba plugin identified candidate genes. Subsequently analyzed using least absolute shrinkage and selection operator (LASSO) and support vector machine recursive feature elimination (SVM-RFE) methods to pinpoint feature genes associated with KD. Following this, biomarkers were determined through expression level validation, enrichment analysis, and immune infiltration analysis. RESULTS: A total of 420 differentially expressed key module genes were identified, and the top 10 genes with DMNC values were selected as candidate genes. Five feature genes were selected through LASSO and SVM-RFE, with NID2, MFAP2, COL8A1, and P4HA3 showing significant expression differences between KD and control samples, along with consistent expression patterns across datasets, identified as potential biomarkers. These four biomarkers were proved to possess high diagnostic potential, and they were found to exhibit significant positive correlations with one another. Functional enrichment analysis indicated that the primary KEGG pathways associated with these biomarkers included "steroid hormone biosynthesis" and "cytokine-cytokine receptor interaction." Moreover, immune infiltration analysis revealed that the four biomarkers were negatively correlated with type 17 T helper cells and positively correlated with 15 immune cell types, including activated B cells and central memory CD4 T cells. CONCLUSION: In conclusion, NID2, MFAP2, COL8A1, and P4HA3 were identified as key biomarkers for KD, offering new avenues for more targeted and effective diagnostic and therapeutic strategies for managing this condition.

Humans↗

PubMed2025

Identification of mitophagy-related biomarkers with immune cell infiltration in psoriasis.

BACKGROUND: Psoriasis is an inflammatory disorder characterized by scaly erythematous plaques and significant comorbidities. Recent studies have suggested that impaired mitophagy, the cellular mechanism for removing dysfunctional mitochondria, may contribute to the pathogenesis of psoriasis. METHODS: In this study, we analyzed bulk RNA sequencing data from 167 healthy individuals and 177 patients with psoriasis obtained from the Gene Expression Omnibus database (GSE30999 and GSE54456). Mitophagy-related genes were isolated using weighted gene co-expression network analysis. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed and protein-protein interaction networks were constructed for the functional enrichment of genes associated with mitophagy. The correlations between genes associated with mitophagy, signaling pathways, and immune cell infiltration were analyzed. The potential diagnostic value of genes associated with mitophagy was evaluated using receiver operating characteristic (ROC) curves, which were validated in imiquimod-induced psoriatic skin lesions in mice. RESULTS: We identified 3,839 differentially expressed genes between healthy individuals and patients with psoriasis, and 23 genes were selected as hub genes showing a high correlation with mitophagy in psoriasis. GO and KEGG analyses revealed that hub and associated genes were significantly correlated with skin functions, such as epidermal development and keratinocyte differentiation. In addition, mitophagy-related genes were negatively associated with pro-inflammatory and pro-proliferation pathways in psoriasis. Among the immune cells, CD4+ T cells were most significantly affected by mitophagy-related genes. ROC analysis demonstrated that mitophagy-related genes, especially ACER1, C1ORF68, CST6, FLG2, GJB3, GJB5, GPRIN2, KRT2, and SPRR4 were potential biomarkers of psoriasis for use in diagnosis or treatment. CONCLUSIONS: Mitophagy-related genes play crucial roles in psoriasis and have potential use as biomarkers, providing insights into disease mechanisms and therapeutic targets. Further research may lead to the development of new strategies for psoriasis management.

Psoriasis↗

PubMed2026

Hepatic metabolic adaptation to endurance exercise: temporal and sex differences by multiomics integration and validation.

BACKGROUND: Although endurance exercise benefits liver health, sex-specific adaptive trajectories remain unclear. This study mapped dynamic liver adaptation in males and females during prolonged training and identified underlying molecular programs. METHODS: Using publicly available time-resolved liver multi-omics data generated by the Molecular Transducers of Physical Activity Consortium (MoTrPAC), we established a computational pipeline for differential analysis of transcriptomic, proteomic, phosphoproteomic, and metabolomic data with FDR correction, followed by FGSEA pathway enrichment. Kinase activities were inferred through ortholog mapping and PhosphoSitePlus. Cross-omics co-expression networks were constructed using WGCNA and topological overlap to link omics features with physiological phenotypes. For experimental validation, liver tissues were collected from endurance-trained Sprague-Dawley rats, and key nodes were confirmed by Western blotting, qRT-PCR, and immunofluorescence/immunohistochemical staining. Public scRNA-seq data were further integrated to map multi-omics signals to single-cell resolution and assess functional changes in specific cell types. RESULTS: The hepatic response to exercise stress was stage-specific, shifting from early transcriptional activation to later proteomic and metabolic remodeling. Multi-omics integration revealed distinct sex-associated adaptive trajectories: males were more strongly associated with energy metabolism, redox-related programs, and amino acid/organic acid catabolism, whereas females showed prominent membrane lipid remodeling, proteostasis -related programs, and mitochondrial/ribosomal translational features. Single-cell analysis showed that tissue remodeling occurred without major lineage turnover, instead involving altered communication among pre-existing cell communities. Validation of PPP1R3G identified a protein-dominant exercise-responsive marker, supporting the contribution of post-transcriptional or protein-level regulation. CONCLUSIONS: Hepatic adaptation to endurance stress follows a cross-omics evolutionary pattern with sex-specific reprogramming of energy supply and homeostatic maintenance. This time-resolved framework clarifies how exercise improves liver function and supports sex-oriented metabolic interventions and therapeutic target discovery.

Animals↗

PubMed2026

Identification and evaluation of glutamine-related gene characteristics based on multi-omics to predict the prognosis of patients with colorectal cancer.

BACKGROUND: Colorectal cancer (CRC), a prevalent malignancy of the gastrointestinal tract, ranks among the leading causes of cancer-related morbidity and mortality. Its clinical course is marked by high fatality and poor prognosis. Elucidating the mechanisms underlying CRC initiation and recurrence is therefore critical for identifying novel therapeutic targets. METHODS: This study incorporated two datasets, TCGA-CRC and GSE17537. A total of 84 glutamine metabolism-related genes (GMRGs) were identified, and differential expression analysis was conducted using the TCGA-CRC dataset. Weighted Gene Co-expression Network Analysis (WGCNA) was applied to determine gene modules most strongly associated with GMRG scores. Single-cell RNA sequencing (scRNA-seq) was utilized to characterize key cellular clusters and to identify differentially expressed genes (DEGs) between high and low glutamine metabolism (GM) groups. Overlapping GMRGs were visualized using the ggVennDiagram package in R. A CRC risk prediction model was developed through Cox proportional hazards and LASSO regression analyses, with performance evaluated by ROC curves. Cell type enrichment across 64 immune and stromal populations was assessed via xCell, and intergroup differences were tested using the Wilcoxon rank-sum test. TIDE scores were used to estimate immunotherapy responsiveness, while oncoPredict facilitated drug sensitivity profiling. PCOLCE2 expression in CRC was validated by RT-qPCR and Western blotting. Its functional role was examined through CCK-8 assays, invasion and migration tests, flow cytometry, and glutamate quantification. RESULTS: ScRNA-seq analysis identified two key cell populations and 437 DEGs associated with GM status. WGCNA pinpointed the MEgreen module as most significantly correlated with GMRG scores, encompassing 1075 genes. Integration of DEGs, module genes, and GM-related DEGs yielded 60 candidate genes for downstream analysis. A GMRG-based prognostic model comprising six genes (SRPX, CXCL1, GPX3, PCOLCE2, CLU, SEMA3E) demonstrated strong predictive performance. Prognostic gene expression correlated with immune and stromal infiltration patterns, as indicated by Spearman correlation analysis. The high-risk group exhibited diminished predicted response to immunotherapy (TIDE scores). Drug sensitivity analysis identified four compounds—Dasatinib-51, WH-4-023-56, TWS-119-366, and LDN-193189-478—with elevated efficacy in high-risk CRC cases. PCOLCE2 expression was significantly reduced in CRC tissues. Functional assays revealed that PCOLCE2 knockdown did not substantially affect cell proliferation but significantly impaired invasion and migration in CRC cells, increased apoptosis, and suppressed both glutamine uptake and glutamate production—highlighting its oncogenic role. CONCLUSION: Six GMRGs—SRPX, CXCL1, GPX3, PCOLCE2, CLU, and SEMA3E—were identified as key components of a robust prognostic model for CRC. These findings offer valuable insights into CRC pathogenesis and potential therapeutic strategies. Notably, this study provides the first evidence implicating PCOLCE2 as a tumor-promoting factor in CRC.

Glutamine↗

PubMed2025

Integrative multi-omics analysis proposes a metabolic classification of gliomas: distinct metabolic states, immune infiltration, and prognosis.

BACKGROUND: The tumor microenvironment (TME) of glioma harbors diverse cell types; however, cell metabolic heterogeneity remains to be explored. This study aims to characterize the metabolic features of different cell types in the TME by integrating multiple datasets, including genomics, bulk and single-cell transcriptomics, and metabolomics. METHODS: Unsupervised machine learning was used to construct an energy metabolic classifier based on the metabolic pathways identified from bulk RNA-seq of gliomas in the TCGA dataset. The classifier was externally validated using multiple datasets, including genomics, bulk RNA-seq, snRNA-seq, and the metabolomics data. Furthermore, metabolic heterogeneity associated with the classifier was further characterized at single-cell resolution. RESULTS: The energy metabolism-based classifier stratified patients into two prognostic clusters: patients in cluster 1 were characterized by high pathway activity of glycolysis, the pentose phosphate pathway (PPP), and fatty acid oxidation (FAO), whereas patients in cluster 2 exhibited higher activity in glutaminolysis. This metabolic classifier revealed both intratumoral and intertumoral metabolic heterogeneity, and the complexity was further validated by the metabolomics profiling and snRNA-seq data from the CPTAC dataset. Notably, OSMR, highly expressed in cluster 1, showed significant co-expression with key glycolytic enzyme genes. The OSM/OSMR/JAK1/STAT3 axis potently drives malignant progression of glioma cells, specially enhancing their invasive and migratory capabilities. Single-cell resolution analyses demonstrated that tumor metabolic heterogeneity is primarily driven by malignant cells rather than non-malignant components, while tumor microenvironment (TME) factors were also found to modulate malignant cell metabolism. Significantly, glycolytic activity in glioma cells increased during the phenotypic transition from PN (proneural) to MES (mesenchymal), with cluster 1 metabolic phenotypes predominating in the tumor core. Compared to cluster 2, cluster 1 patients exhibited higher mRNA expression of immunosuppressive checkpoint genes, which correlated with pronounced immunosuppression in the TME. Furthermore, various immune cells demonstrated distinct metabolic preferences at single-cell resolution. CONCLUSIONS: This study developed an energy metabolic-based classifier for gliomas with prognostic and therapeutic potential. Metabolic reprogramming was linked with the PN-to-MES transition of glioma cells and immunosuppression in the tumor microenvironment. Multi-omics data, especially snRNA-seq, offered insights into metabolism heterogeneity at single-cell resolution, enabling personalized treatment strategies.

Humans↗

PubMed2026

A modular γδ TCR-T platform combining KRAS pMHC targeting with re-dosable mRNA engager redirection.

Solid tumors often evade TCR-engineered αβ T cells when antigen expression varies or when the restricting Human Leukocyte Antigen (HLA) allele is lost. γδ T cells, in contrast, detect cellular dysregulation through non-peptide/Major Histocompatibility Complex (MHC) cues, including phosphoantigens and stress ligands, and can be developed as allogeneic therapies. Although intratumoral γδ T cell signatures are associated with improved outcome across cancers, γδ recognition itself is broad and still selected within the thymus just as αβ T cell receptors (TCRs) are. It does not, however, anchor specificity to a defined driver-mutation pMHC epitope. We therefore asked whether a high-affinity, co-receptor-independent αβ TCR could graft oncogenic-driver specificity onto γδ T cells while leaving the endogenous γδ TCR intact. We knocked the KRASG12V/HLA-A*11:01 TCR A11v into primary human γδ T cells. Engineered cells co-expressed the transgenic αβ TCR and the endogenous γδ TCR and lysed KRASG12V/HLA-A*11:01+ tumor cells in vitro and in vivo. To cover potential resistance through loss of HLA-A*11:01, we delivered an mRNA lipid nanoparticle (LNP) encoding a secreted mesothelin×CD3 (M5) bispecific T cell engager (TCE). LNP-M5 produced circulating TCE that redirected γδ A11v T cells and polyclonal bystander T cells to kill mesothelin+ targets, accompanied by development of higher γδ A11v T cell counts in vivo. In humanized mice bearing mixed HLA-A*11:01+ and HLA-A*11:01 - KRASG12V tumors, γδ A11v T cells produced transient control, whereas adding LNP-M5 yielded complete responses and prolonged survival. Thus, this two-part therapy couples invariant driver targeting to tunable redirection and addresses loss of the restricting HLA allele, a central escape route for TCR-based therapy. It provides an off-the-shelf reagent to enable KRAS-anchored treatment with the ability to redeliver the reagent.

Humans↗

PubMed2025

Foxi2 and Sox3 are master transcription regulators that control ectoderm germ layer specification in Xenopus.

Germ layer specification represents a critical transition where pluripotent cells acquire lineage-specific identities. We identify the maternal transcription factors Foxi2 and Sox3 to be pivotal master regulators of ectodermal germ layer specification in Xenopus. Ectopic co-expression of Foxi2 and Sox3 in prospective endodermal tissue induces the expression of ectodermal markers while suppressing mesendodermal markers. Transcriptomic analyses reveal that Foxi2 and Sox3 jointly and independently regulate hundreds of ectodermal target genes. During early cleavage stages, Foxi2 and Sox3 pre-bind to key cis-regulatory modules (CRMs), marking sites that later recruit Ep300 and facilitate H3K27ac deposition, thereby shaping the epigenetic landscape of the ectodermal genome. These CRMs are highly enriched within ectoderm-specific super-enhancers (SEs). Our findings highlight the pivotal role of ectodermal SE-associated CRMs in precise and robust ectodermal gene activation, establishing Foxi2 and Sox3 as central architects of ectodermal lineage specification.

Animals↗

PubMed2024

Inhibiting the expression of spindle appendix cooled coil protein 1 can suppress tumor cell growth and metastasis and is associated with cancer immune cells in esophageal squamous cell carcinoma.

Inhibiting the expression of spindle appendix cooled coil protein 1 (SPDL1) can slow down disease progression and is related to poor prognosis in patients with esophageal cancer. However, the specific roles and molecular mechanisms of SPDL1 in esophageal squamous cell carcinoma (ESCC) have not been explored yet. The current study aimed to investigate the expression levels of SPDL1 in ESCC via transcriptome analysis using data from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus databases. Moreover, the biological roles, molecular mechanisms, and protein networks involved in SPDL1 were identified using machine learning and bioinformatics. The cell counting kit-8 assay, EdU staining, and transwell assay were used to investigate the effects of inhibiting SPDL1 expression on ESCC cell proliferation, migration, and invasion. Finally, the correlation between the SPDL1 expression and cancer immune infiltrating cells was evaluated by analyzing data from the TCGA database. Results showed that SPDL1 was overexpressed in the ESCC tissues. The SPDL1 expression was related to age in patients with ESCC. The SPDL1 co-expressed genes included those involved in cell division, cell cycle, DNA repair and replication, cell aging, and other processes. The high-risk scores of SPDL1-related long non-coding RNAs were significantly correlated with overall survival and cancer progression in patients with ESCC (P < 0.05). Inhibiting the SPDL1 expression was effective in suppressing the proliferation, migration, and invasion of ESCC TE-1 cells (P < 0.05). The overexpression of SPDL1 was positively correlated with the levels of Th2 and T-helper cells, and was negatively correlated with the levels of plasmacytoid dendritic cells and mast cells. In conclusion, SPDL1 was overexpressed in ESCC and was associated with immune cells. Further, inhibiting the SPDL1 expression could effectively slow down cancer cell growth and migration. SPDL1 is a promising biomarker for treating patients with ESCC.

Humans↗

PubMed2024

Transcriptome analysis under pecan scab infection reveals the molecular mechanisms of the defense response in pecans.

Pecan scab, caused by the fungal pathogen Venturia effusa, is the most devastating disease of pecan (Carya illinoinensis) in the southeastern United States. Resistance to this pathogen is determined by a complex interaction between host genetics and disease pathotype with even field-susceptible cultivars being resistant to most scab isolates. To understand the underlying molecular mechanisms of scab resistance in pecan, we performed a transcriptome analysis of the pecan cultivar, 'Desirable', in response to inoculation with a pathogenic and a non-pathogenic scab isolate at three different time points (24, 48, and 96 hrs. post-inoculation). Differential gene expression and gene ontology enrichment analyses showed contrasting gene expression patterns and pathway enrichment in response to the contrasting isolates with varying pathogenicity. The weighted gene co-expression network analysis of differentially expressed genes detected 11 gene modules. Among them, two modules had significant enrichment of genes involved with defense responses. These genes were particularly upregulated in the resistant reaction at the early stage of fungal infection (24 h) compared to the susceptible reaction. Hub genes in these modules were predominantly related to receptor-like protein kinase activity, signal reception, signal transduction, biosynthesis and transport of plant secondary metabolites, and oxidoreductase activity. Results of this study suggest that the early response of pathogen-related signal transduction and development of cellular barriers against the invading fungus are likely defense mechanisms employed by pecan cultivars against non-virulent scab isolates. The transcriptomic data generated here provide the foundation for identifying candidate resistance genes in pecan against V. effusa and for exploring the molecular mechanisms of disease resistance.

Carya↗

PubMed2026

Metabolic-cell-death gene trio predicts survival and cuproptosis sensitivity in colorectal cancer.

BACKGROUND: Metabolic cell death (MCD) modulates colorectal cancer (CRC) progression, yet its prognostic value remains unexplored. We aimed to build an MCD-centred gene signature for outcome prediction and precision therapy. METHODS: Transcriptomes of 1,174 CRC patients were integrated. Weighted gene co-expression network analysis, differential expressions and least absolute shrinkage and selection operator (LASSO) + random survival forest were successively applied to derive a three-gene (CDKN2A/MPC1/AHCY) risk model. Functional, immune-infiltration, drug-sensitivity and genomic analyses were performed, followed by validation in fresh clinical specimens and cell lines. RESULTS: Integrative metabolic-death transcriptomics identified CDKN2A, MPC1 and AHCY as the hub drivers of CRC. Their three-gene signature robustly stratified patients into high- and low-risk subsets [3-year area under the curve (AUC) 0.83-0.85, P<0.001]. High-risk tumors were enriched for extracellular matrix (ECM)-receptor-interaction pathways, displayed abundant myeloid-derived suppressor cell (MDSC) infiltration and were more vulnerable to AZD8186, AZ960 and JAK inhibitors. Guided by these in-silico findings, we functionally confirmed that CDKN2A silencing markedly repressed proliferation, invasion and migration of SW480/HCT116 cells and potentiated cuproptosis via up-regulation of lipoylated DLAT/DLST and CTR1. CONCLUSIONS: We report the first MCD-derived prognostic platform for CRC that simultaneously predicts survival and therapeutic response. Targeting CDKN2A-enhanced cuproptosis represents a promising metabolic-precision strategy for high-risk patients.

Colorectal cancer (CRC)↗

PubMed2026

Enterocutaneous Fistula-Associated Sepsis and Mortality: Development and Validation of a Multimodal Artificial Intelligence Prediction Model.

BACKGROUND: Predicting enterocutaneous fistula (ECF)-associated sepsis and mortality poses significant challenges in digital health care due to the disease's complexity and heterogeneous clinical manifestations. Current approaches that rely on single-modal data or traditional scoring systems often fail to capture the intricate immune-inflammatory dynamics and multisystem involvement in patients with ECF. OBJECTIVE: This study aims to develop an artificial intelligence (AI)-driven multimodal fusion model integrating clinical, imaging, and transcriptomic data for early prediction of ECF-associated sepsis and 28-day mortality, addressing the limitations of conventional single-dimensional models. METHODS: This study leveraged publicly available datasets (Medical Information Mart for Intensive Care III [MIMIC-III], electronic Intensive Care Unit [eICU], and The Cancer Genome Atlas) to construct a multimodal framework. Clinical parameters were processed using Extreme Gradient Boosting, abdominal imaging features were extracted via convolutional neural networks, and transcriptomic profiles were analyzed with variational autoencoders. A Transformer-based fusion network was employed for joint prediction and validated through cross-validation and external testing. Key features were identified using Shapley Additive Explanations and Local Interpretable Model-Agnostic Explanations interpretability algorithms, while immune regulatory mechanisms were explored via weighted gene co-expression network analysis. RESULTS: The multimodal model achieved an area under the curve (AUC) of 0.89 for predicting sepsis and 28-day mortality, outperforming unimodal models (clinical-only model, AUC 0.72, and imaging-only model, AUC 0.78). Critical predictors included Sequential Organ Failure Assessment score, lactate levels, intra-abdominal free fluid on imaging, and immunoregulatory genes (programmed death-ligand 1 [PD-L1] and indoleamine 2,3-dioxygenase 1 [IDO1]). Mechanistic analysis revealed distinct immune reprogramming in patients with sepsis, characterized by increased regulatory T cells and M2 macrophages, along with downregulated cluster of differentiation 8+ (CD8+) T cells. CONCLUSIONS: This multimodal AI model offers an innovative digital solution in medical informatics, enabling precise early risk stratification for ECF-associated sepsis. By integrating multisource data and providing interpretable insights into immune-inflammatory pathways, the model enhances health care quality for patients with ECF and paves the way for personalized intervention strategies.

Humans↗

PubMed2026

Gut microbiota-derived metabolites target C5AR1/KDM2A/HCAR3 axis in inflammatory bowel disease: a multi-machine learning algorithms and molecular docking study.

BACKGROUND: Inflammatory bowel disease (IBD) is a chronic recurrent disorder. Gut microbiota-derived metabolites regulate intestinal homeostasis, but their molecular mechanisms in IBD remain unclear. Current studies lack systematic "microbiota-metabolite-target" network mining with multi-method validation. This study integrates network pharmacology, three machine learning algorithms, and molecular docking to construct this regulatory network in IBD. METHODS: Transcriptome data were obtained from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) were identified using limma (p < 0.05, |log2FC| > 0.5). Weighted gene co-expression network analysis (WGCNA) with an optimal soft threshold of β = 7 was performed to identify key module genes. Candidate genes were obtained by intersecting DEGs, gut microbiota-associated genes from the gutMGene database, and WGCNA module genes. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were conducted to explore the functional roles of candidate genes. Core genes were identified using three machine learning algorithms (LASSO, Boruta, and SVM-RFE), followed by protein-protein interaction (PPI) network analysis. Molecular docking was performed to assess the binding affinities between hub proteins and gut microbiota-derived metabolites. RESULTS: A total of 885 DEGs were identified between the IBD and control groups, including 463 upregulated and 422 downregulated genes. WGCNA identified 280 key module genes from the purple and yellow modules. The intersection of DEGs, gut microbiota-associated genes, and WGCNA module genes yielded 19 core candidate genes. PPI network analysis combined with three machine learning algorithms jointly identified C5AR1, KDM2A, and HCAR3 as core hub genes. ROC curve analysis demonstrated that all three hub genes achieved AUC values greater than 0.7 in both the training and validation sets, indicating excellent diagnostic performance for IBD. Enrichment analysis revealed significant associations with the TNF, NF-κB, and IL-17 signaling pathways. Molecular docking confirmed stable binding of C5AR1 with 1,3-Diphenylpropan-2-Ol (-7.87 ± 0.83 kcal·mol-¹) and HCAR3 with 3-Indolepropionic Acid (-6.35 ± 0.70 kcal·mol-¹), both below -5.0 kcal·mol-¹. CONCLUSION: This study first constructs a "gut microbiota-metabolite-hub gene" axis in IBD, providing a computational framework for microbiota-targeted precision therapy, and identifying C5AR1/KDM2A/HCAR3 as computationally predicted diagnostic biomarkers and 1,3-Diphenylpropan-2-Ol/3-Indolepropionic Acid as candidate intervention molecules that warrant further experimental validation.

Molecular Docking Simulation↗

PubMed2026

Genome-wide cis-expression Quantitative Trait Loci (eQTL) and transcriptomic signals reveal distinct molecular regulation across correlated feed efficiency traits.

INTRODUCTION: Feed efficiency (FE) is a complex trait which determines livestock production profitability, yet the molecular mechanisms behind it remain unclear. This study investigated the blood transcriptomic profile of lambs, alongside genotype data with the aim to uncover the genetic basis of FE traits such as absolute dry matter intake (DMIabsolute), DMI adjusted for body size (DMIadjusted), average daily live weight gain (ADG), and residual feed intake (RFI). MATERIALS AND METHODS: Bulk RNA-Seq and genotype data were analysed using three complementary approaches: differential gene expression (DGE) analysis, weighted gene co-expression network analysis (WGCNA), and cis-expression Quantitative Trait Loci (cis-eQTL) mapping. These methods were used independently to identify genes and regulatory networks associated with FE traits and to investigate evidence supporting multi-trait candidate gene selection. RESULTS: DGE analysis revealed 2, 24, 85 and 4 differentially expressed genes for DMIabsolute, DMIadjusted, ADG, and RFI (Padjusted < 0.05), functionally enriched in sensory perception, ATP-dependent chromatin remodeling, Notch signaling and immune response pathways. 9 gene modules significantly associated with the FE traits (P ≤ 0.05) with correlations ranging from r = -0.56 to 0.49, were identified using WGCNA. Single nucleotide polymorphism (SNP)-level cis-eQTL analysis identified 93 eSNPs associated with 74 genes (false discovery rate (FDR) < 0.05), while permutation-derived gene level analysis identified 280 eGenes (FDR < 0.2, empirical P < 0.03). Across the three analyses, applying thresholds of DGE (Padjusted < 0.05), WGCNA (correlation, P ≤ 0.05), and cis-eQTL gene-level significance (empirical P < 0.05), multiple overlapping genes were identified including DNMT3A, KANSL1, NCOR1 for DMIadjusted, ACOX2, FANCF, CIMIP2B, LOC101115106, ARMH2, LOC132657496 for ADG, and LOC114114576 for RFI representing regulators of variations in FE. DISCUSSION: The integration of DGE, WGCNA, and cis-eQTL analyses identified key genes and regulatory mechanisms associated with variation in FE traits. These results highlight that integrated multi-trait candidate gene identification approaches can reveal key genes that lower feed intake while maintaining animal growth, supporting breeding strategies aimed at improving efficiency and long-term economic sustainability in sheep.

average daily gain (ADG)↗

PubMed2026

ZBTB16-associated NK cell alterations reveal shared immunometabolic signatures linking primary Sjögren's syndrome and type 1 diabetes mellitus.

BACKGROUND: Primary Sjögren's syndrome (pSS) and type 1 diabetes mellitus (T1DM) share immune-inflammatory features, yet conserved pathogenic signatures linking these autoimmune disorders remain incompletely understood. The present research sought to uncover common molecular markers and dissect the underlying immune-metabolic cross-talk underlying pSS and T1DM. METHODS: Gene expression profiles of patients with pSS and T1DM were retrieved from the Gene Expression Omnibus database, normalized, and corrected for batch effects prior to downstream analyses. Overlapping potential biomarkers were screened by integrating differential expression analysis, weighted gene co-expression network analysis and least absolute shrinkage and selection operator regression. Functional enrichment based on Gene Ontology and Kyoto Encyclopedia of Genes and Genomes databases was implemented to interpret gene biological properties, and a protein-protein interaction network was further established afterwards. Diagnostic performance was evaluated using receiver operating characteristic analysis. Experimental validation was conducted in non-obese diabetic (NOD) mice using quantitative PCR, immunohistochemistry, and flow cytometry. The CIBERSORT algorithm was adopted to quantify immune cell infiltration levels. RESULTS: ZBTB16 was identified as a shared hub biomarker in both pSS and T1DM and exhibited favorable diagnostic performance. Experimental validation confirmed significantly reduced ZBTB16 expression in peripheral blood mononuclear cells, salivary gland tissues, and pancreatic tissues of NOD mice. Gene Set Enrichment Analysis indicated that ZBTB16-associated signatures were enriched in mitochondrial-related processes, neuroactive ligand-receptor interactions, and ribosome-related pathways. Immune infiltration analysis revealed that resting natural killer (NK) cells were positively correlated with ZBTB16 expression in both diseases. Flow cytometric analysis further confirmed a reduced proportion of resting NK cells in peripheral blood of NOD mice, consistent with the CIBERSORT-based prediction. CONCLUSION: This study identifies ZBTB16 as a shared biomarker linking pSS and T1DM. Reduced resting NK-cell abundance was consistently observed in both computational and experimental analyses, and bioinformatic correlation analysis suggested a positive association with ZBTB16 expression. These findings provide evidence for shared molecular and immunological signatures underlying the two autoimmune disorders and support further investigation of the biological role and diagnostic value of ZBTB16 in pSS and T1DM.

Sjogren's Syndrome↗

PubMed2025

Multi-Omics and Integrative Analytics in Natural Products Discovery.

Natural products (NPs) have long been an essential source of new bioactive compounds for drug discovery; however, traditional methods for screening and isolating these compounds can be slow and often yield diminishing returns. Fortunately, advanced multi-omics and computational approaches present powerful solutions to these challenges. This review highlights innovative methodologies that integrate metabolomics, genomics, transcriptomics, and proteomics with bioinformatics and analytical chemistry to accelerate NP discovery. For instance, untargeted metabolomics platforms like high-resolution liquid chromatography-tandem mass spectrometry (LC-MS/MS) and Global Natural Products Social (GNPS) molecular networking allow for comprehensive profiling of new compounds, while targeted isotope-labeling strategies enhance this process. Additionally, genome and metagenome mining tools such as antibiotics and secondary metabolite analysis shell (antiSMASH), Deep Biosynthetic Gene Cluster (DeepBGC), and Pipeline for Reconstructing Integrated Syntheses of Metabolites (PRISM) quickly identify biosynthetic gene clusters (BGCs) in both cultured and uncultured organisms, often using heterologous expression to validate products. Transcriptomic analyses, including RNA sequencing (RNA-seq), co-expression networks, and fluxomics, help clarify how pathways are regulated, while quantitative proteomics techniques like tandem mass tags/isobaric tags for relative and absolute quantitation (TMT/iTRAQ) and label-free methods, along with chemoproteomics approaches such as cellular thermal shift assay and thermal proteome profiling (TPP), uncover molecular targets and their mechanisms of action. This review also places significant emphasis on the role of artificial intelligence (AI) and machine learning (ML) in integrating multi-omics data, spanning activities from constructing gene-metabolite correlation networks to leveraging knowledge graphs and graph neural networks for data fusion and functional prediction. Finally, this review concludes by discussing the synergistic benefits of multi-omics for natural-product discovery, addressing current technical challenges, and exploring future directions toward high-throughput, intelligent data integration for next-generation NP research.

Biological Products↗

PubMed2026

S100P as a Shared Biomarker in Inflammatory Bowel Disease, Colorectal Cancer, and Pancreatic Adenocarcinoma: An Integrated Transcriptomic Analysis.

Inflammatory bowel disease (IBD) is associated with an increased risk of colorectal cancer (CRC) and pancreatic adenocarcinoma (PAAD), yet the molecular features shared among these diseases remain incompletely understood. This study aimed to identify common genes and biological pathways associated with IBD, CRC, and PAAD through integrated transcriptomic analysis and experimental validation. Gene expression datasets for IBD, CRC, and PAAD were obtained from The Cancer Genome Atlas and Gene Expression Omnibus databases. Weighted gene co-expression network analysis and differential expression analysis were performed to identify disease-associated and shared genes. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes (analyses were used to explore enriched biological functions and pathways. Immune cell infiltration was evaluated using Cell-type Identification by Estimating Relative Subsets of RNA Transcripts. Receiver operating characteristic analysis was performed to assess the diagnostic performance of common genes. Single-cell RNA sequencing analysis was conducted to examine the cellular distribution of S100P. In addition, the effects of S100P downregulation were evaluated in lipopolysaccharide (LPS)-stimulated colonic epithelial cells. A total of 162 disease-associated genes and four common genes were identified. Functional enrichment analyses indicated significant enrichment of immune- and inflammation-related pathways, including the interleukin-17 signaling pathway. Immune infiltration analysis revealed similar trends in several immune cell populations across IBD, CRC, and PAAD. Single-cell analysis showed elevated S100P expression in epithelial cells from all three diseases. Downregulation of S100P restored the proliferative capacity of LPS-stimulated colonic epithelial cells and reduced inflammatory cytokine expression. Integrated transcriptomic analysis identified S100P as a biomarker associated with IBD, CRC, and PAAD and highlighted shared immune-related features across these diseases.

Humans↗