A gut microbiota-metabolite-host gene score for subtyping and predicting the treatment response in rheumatoid arthritis: integrative analysis with in silico knockout validation.
- 1School of Management, Shanxi Medical University, Taiyuan, China.
- 2School of Management, Shanxi Medical University, Taiyuan, China. [email protected].
- 3School of Management, Shanxi Medical University, Taiyuan, China. [email protected].
Abstract
Background
The pathogenesis of rheumatoid arthritis (RA) involves a complex interplay between immune dysregulation and metabolic disturbance, with the gut microbiota and its metabolites playing pivotal roles. However, integrated biomarkers based on the "gut microbiota-metabolite-host gene" regulatory axis that can quantify its impact on RA stratification and therapeutic response are lacking.Methods
We integrated multiple RA transcriptomic datasets to screen for host genes explicitly associated with the gut microbiota and metabolites using the gutMGene database. A gut microbiota-metabolite-host score (GMS) was constructed via principal component analysis (PCA). Immune infiltration patterns were evaluated using CIBERSORT. Furthermore, scRNA-seq was employed to map the cellular origins of key genes, and in silico knockout analysis using scTenifoldKnk was performed to simulate the functional ablation of core host genes within this regulatory axis. Additionally, molecular docking was conducted to validate the binding affinities between the key host proteins and the optimal small molecule.Results
We identified seven microbiota-metabolite-host hub genes (FOS, JUN, IDO1, PEMT, EGFR, ABCC4, and FADS3) and constructed the GMS. The GMS stratified patients with RA into two distinct functional subtypes: an adaptive immune activation subtype and an inflammatory-metabolic composite subtype. Notably, the GMS demonstrated superior performance in predicting responses to anti-TNF and rituximab therapies. The scRNA-seq analysis revealed high expression of FOS and JUN in T cells and NK cells. In silico knockout analysis further revealed the central transcription factors FOS and JUN as key executors of this axis, as their ablation significantly affected osteoclast differentiation and the enrichment of TNF/MAPK inflammatory pathways. These computational findings support the core mechanism by which gut microbiota-metabolite signalling drives RA joint pathology via FOS/JUN. Molecular docking results confirmed strong binding interactions between the identified active molecules and their respective target proteins.Conclusion
This study establishes a GMS based on the gut microbiota-metabolite-host gene regulatory axis as a robust tool for molecular subtyping and predicting biological therapy responses in patients with RA. Crucially, through integrated computational modelling, including in silico knockout and molecular docking validation, we identified FOS/JUN as critical executor genes linking the gut microenvironment to RA bone destruction, providing novel targets for therapeutic strategies aimed at the gut-joint axis.
Related
- Microbiota-derived metabolite signatures in Alzheimer’s disease across the cognitive aging continuum: a human-focused systematic review and evidence-weighted knowledge graph
- Genetic-multi-omics causal inference coupled with integrative in-vitro toxicology reveals immune-metabolite mediators for aortic aneurysm risk
- Integrative population genetics and multi-omics reveal IL15RA as a shared pathogenic hub in rheumatoid arthritis-related heart failure
- Integrative Advances in Honey Bee Gut–Brain Interactions: From Microbiota Dynamics to Behavioral Regulation and Colony Health
- Radish Red Anthocyanins Extract Attenuated High-Fat Diet-Induced Aging in Association with Improving Antioxidant Defense and Gut Microbiota–Metabolite Changes in Drosophila Melanogaster
BibTeX
@article{sun2026microbiotametabolitehost,
title = {{A gut microbiota-metabolite-host gene score for subtyping and predicting the treatment response in rheumatoid arthritis: integrative analysis with in silico knockout validation.}},
author = {Xiang-Fei Sun and Yi-Jie Yang and Chen-Long Li and Tai Huang and Qi Wang and Qi Yu},
journal = {Bioresources and bioprocessing},
year = {2026},
month = oct,
doi = {10.1186/s40643-026-01149-7},
url = {https://europepmc.org/article/MED/42900335},
note = {CC BY},
}
RIS
TY - JOUR
TI - A gut microbiota-metabolite-host gene score for subtyping and predicting the treatment response in rheumatoid arthritis: integrative analysis with in silico knockout validation.
AU - Xiang-Fei Sun
AU - Yi-Jie Yang
AU - Chen-Long Li
AU - Tai Huang
AU - Qi Wang
AU - Qi Yu
PY - 2026
DA - 2026/10/20
JO - Bioresources and bioprocessing
DO - 10.1186/s40643-026-01149-7
UR - https://europepmc.org/article/MED/42900335
AB - Background The pathogenesis of rheumatoid arthritis (RA) involves a complex interplay between immune dysregulation and metabolic disturbance, with the gut microbiota and its metabolites playing pivotal roles. However, integrated biomarkers based on the "gut microbiota-metabolite-host gene" regulatory axis that can quantify its impact on RA stratification and therapeutic response are lacking.Methods We integrated multiple RA transcriptomic datasets to screen for host genes explicitly associated with the gut microbiota and metabolites using the gutMGene database. A gut microbiota-metabolite-host score (GMS) was constructed via principal component analysis (PCA). Immune infiltration patterns were evaluated using CIBERSORT. Furthermore, scRNA-seq was employed to map the cellular origins of key genes, and in silico knockout analysis using scTenifoldKnk was performed to simulate the functional ablation of core host genes within this regulatory axis. Additionally, molecular docking was conducted to validate the binding affinities between the key host proteins and the optimal small molecule.Results We identified seven microbiota-metabolite-host hub genes (FOS, JUN, IDO1, PEMT, EGFR, ABCC4, and FADS3) and constructed the GMS. The GMS stratified patients with RA into two distinct functional subtypes: an adaptive immune activation subtype and an inflammatory-metabolic composite subtype. Notably, the GMS demonstrated superior performance in predicting responses to anti-TNF and rituximab therapies. The scRNA-seq analysis revealed high expression of FOS and JUN in T cells and NK cells. In silico knockout analysis further revealed the central transcription factors FOS and JUN as key executors of this axis, as their ablation significantly affected osteoclast differentiation and the enrichment of TNF/MAPK inflammatory pathways. These computational findings support the core mechanism by which gut microbiota-metabolite signalling drives RA joint pathology via FOS/JUN. Molecular docking results confirmed strong binding interactions between the identified active molecules and their respective target proteins.Conclusion This study establishes a GMS based on the gut microbiota-metabolite-host gene regulatory axis as a robust tool for molecular subtyping and predicting biological therapy responses in patients with RA. Crucially, through integrated computational modelling, including in silico knockout and molecular docking validation, we identified FOS/JUN as critical executor genes linking the gut microenvironment to RA bone destruction, providing novel targets for therapeutic strategies aimed at the gut-joint axis.
SN - 2197-4365
ER -