Integrative single-cell and bulk transcriptomics with reverse network pharmacology prioritizes macrophage-associated candidate genes in Lyme arthritis
5 affiliations
- 1Department of Pain Medicine, The Affiliated Anning First People’s Hospital of Kunming University of Science and Technology
- 2School of Basic Medical Sciences, Kunming University of Science and Technology
- 3Faculty of Life Science and Technology and the Affiliated Anning First People’s Hospital, Kunming University of Science and Technology
- 4Clinical Medical College, Kunming University
- 5First School of Clinical Medicine, Kunming Medical University
Abstract
Introduction: Lyme arthritis (LA) involves complex immune dysregulation, and macrophages contribute centrally to persistent joint inflammation. However, synovial macrophage heterogeneity and its potential relevance to therapeutic intervention remain incompletely defined.
Methods: We integrated scRNA-seq data from Borrelia burgdorferi ( Bb )-infected murine joints (GSE233850) with bulk RNA-seq data (GSE283617). The workflow included Seurat clustering, Monocle 3 trajectory inference, CellChat communication analysis, and functional enrichment. We constructed a PPI network from intersecting genes and applied five complementary topological measures to prioritize consensus hubs. Candidate gene expression was evaluated in an external transcriptomic dataset (GSE63085) and by qRT-PCR. CTD-based association analysis and molecular docking were used to prioritize candidate compounds. Virtual knockout analysis explored the potential regulatory effects of hub genes.
Results: The two pooled libraries differed descriptively in the relative representation of nine major cell populations. Among macrophages/monocytes, we resolved seven marker-defined states. Monocle 3 trajectories were rooted at Classical_Inflam_Mono, and CellChat highlighted Tgfb1-(Tgfbr1+Tgfbr2) as a prominent inferred TGF-β axis. Intersecting GSE283617 with GSE233850 yielded 121 shared genes, and GSEA showed negative enrichment of oxidative phosphorylation during infection. Mrc1 and Msr1 were the only genes ranked among the top 10 by all five measures: Degree, Bottleneck, EcCentricity, MNC, and MCC. Across these complementary analyses, Msr1 emerged as a candidate Bb -responsive gene potentially relevant to LA. A CTD-based source-compound-gene network and molecular docking prioritized curcumin and demethoxycurcumin for future experimental testing.
Conclusion: These findings provide an exploratory atlas of macrophage transcriptional states in the analyzed LA model and identify Msr1 as a candidate infection-responsive gene. Its disease specificity and clinical utility require validation in joint-derived specimens from patients with clinically confirmed LA. The computationally prioritized compounds remain hypotheses that require direct binding and functional testing.
Related
- Integrative transcriptomics and reverse network pharmacology prioritize DPP4 and Achyranthes bidentata in Barrett’s esophagus
- Integrative Multiomics Mendelian Randomization Prioritizes Endoplasmic Reticulum Stress-Related Genes Associated With Gastric Cancer Risk.
- Integrative Multi-Omics and Single-Cell Analysis Prioritizes YWHAB Within Ubiquitination-Related Transcriptional Programs in Endometrial Carcinoma
- Integrative Transcriptomics and Machine Learning Identify Temporal Stress Signatures and Astrocytic Neurocan Associated With Hypothalamic Injury.
- An evidence-governed agent harness for spatial transcriptomics prioritizes ICAM-1 as a target in lung fibrosis
BibTeX
@article{dong2026integrative,
title = {{Integrative single-cell and bulk transcriptomics with reverse network pharmacology prioritizes macrophage-associated candidate genes in Lyme arthritis}},
author = {Yan Dong and Yanshuang Luo and Yantong Chen and Meng Liu and Yingchen Zhou and Yutian Yang and Xuesong Chen and Fusong Yang and Zhi Liang and Liye Zi and Hongbing Gao and Zhichao Cao and Fucun Gao and Chao Song and Guozhong Zhou},
journal = {Frontiers in Immunology},
year = {2026},
month = oct,
doi = {10.3389/fimmu.2026.1938154},
url = {https://www.frontiersin.org/articles/10.3389/fimmu.2026.1938154/full},
note = {CC BY 4.0},
}
RIS
TY - JOUR
TI - Integrative single-cell and bulk transcriptomics with reverse network pharmacology prioritizes macrophage-associated candidate genes in Lyme arthritis
AU - Yan Dong
AU - Yanshuang Luo
AU - Yantong Chen
AU - Meng Liu
AU - Yingchen Zhou
AU - Yutian Yang
AU - Xuesong Chen
AU - Fusong Yang
AU - Zhi Liang
AU - Liye Zi
AU - Hongbing Gao
AU - Zhichao Cao
AU - Fucun Gao
AU - Chao Song
AU - Guozhong Zhou
PY - 2026
DA - 2026/10/13
JO - Frontiers in Immunology
DO - 10.3389/fimmu.2026.1938154
UR - https://www.frontiersin.org/articles/10.3389/fimmu.2026.1938154/full
AB - Introduction: Lyme arthritis (LA) involves complex immune dysregulation, and macrophages contribute centrally to persistent joint inflammation. However, synovial macrophage heterogeneity and its potential relevance to therapeutic intervention remain incompletely defined. Methods: We integrated scRNA-seq data from Borrelia burgdorferi ( Bb )-infected murine joints (GSE233850) with bulk RNA-seq data (GSE283617). The workflow included Seurat clustering, Monocle 3 trajectory inference, CellChat communication analysis, and functional enrichment. We constructed a PPI network from intersecting genes and applied five complementary topological measures to prioritize consensus hubs. Candidate gene expression was evaluated in an external transcriptomic dataset (GSE63085) and by qRT-PCR. CTD-based association analysis and molecular docking were used to prioritize candidate compounds. Virtual knockout analysis explored the potential regulatory effects of hub genes. Results: The two pooled libraries differed descriptively in the relative representation of nine major cell populations. Among macrophages/monocytes, we resolved seven marker-defined states. Monocle 3 trajectories were rooted at Classical_Inflam_Mono, and CellChat highlighted Tgfb1-(Tgfbr1+Tgfbr2) as a prominent inferred TGF-β axis. Intersecting GSE283617 with GSE233850 yielded 121 shared genes, and GSEA showed negative enrichment of oxidative phosphorylation during infection. Mrc1 and Msr1 were the only genes ranked among the top 10 by all five measures: Degree, Bottleneck, EcCentricity, MNC, and MCC. Across these complementary analyses, Msr1 emerged as a candidate Bb -responsive gene potentially relevant to LA. A CTD-based source-compound-gene network and molecular docking prioritized curcumin and demethoxycurcumin for future experimental testing. Conclusion: These findings provide an exploratory atlas of macrophage transcriptional states in the analyzed LA model and identify Msr1 as a candidate infection-responsive gene. Its disease specificity and clinical utility require validation in joint-derived specimens from patients with clinically confirmed LA. The computationally prioritized compounds remain hypotheses that require direct binding and functional testing.
SN - 1664-3224
ER -