bioRxivGenomicsPreprint
Integrative Analysis of Single-cell and Bulk RNA Sequencing Reveals Macrophage Heterogeneity in Lyme Arthritis
Yanhong DongiD, Yantong CheniD, Yanshuang LuoiD, Meng LiuiD, Yingchen Zhou, Yutian Yang, Xuesong CheniD, Xiaorong LiuiD, Fusong Yang, Zhi Liang, Liye Zi, Hongbin GaoiD, Zhichao Cao, Fucun GaoiD, Chao Ling SongiD, Guozhong ZhouiD
Abstract
IntroductionLyme arthritis (LA) is a common late-stage manifestation of Lyme disease caused by Borrelia burgdorferi (Bb) infection. Macrophages play a crucial role in LA pathogenesis, yet their heterogeneity and functional states within the arthritic microenvironment remain poorly understood. MethodsWe integrated single-cell RNA sequencing data from ankle joint tissues of Bb-infected C57BL/6 mice with bulk RNA sequencing data from both infected ankle joints and bone marrow-derived macrophages. We developed a LA mouse model and validated key findings through histopathological examination and RT-qPCR analysis. ResultsWe identified nine distinct cell types in the joint microenvironment with significant compositional shifts during disease progression. Analysis of 16,617 Macrophages/Monocytes revealed seven functionally distinct subpopulations with differential responses to infection. Pseudo-time trajectory analysis demonstrated three principal differentiation pathways from monocyte-derived macrophages toward anti-inflammatory, TNF+, and IL1b+ inflammatory phenotypes. Cell-cell communication analysis revealed significantly altered TGF-{beta} signaling networks in LA, with the Tgfb1-(Tglbr1+TgIbr2) axis playing a critical role. GSVA and pathway analyses showed metabolic reprogramming from oxidative phosphorylation toward glycolysis in macrophages during infection. Through integrated analysis and LASSO regression, we identified four characteristic genes (SIRPB1C, FABP5, MMP14, and EGR1) as potential LA biomarkers and developed a diagnostic nomogram with high predictive accuracy. ConclusionOur integrative analysis provides comprehensive insights into macrophage heterogeneity and functional plasticity in LA, identifying characteristic biomarkers and potential therapeutic targets. The metabolic reprogramming and altered TGF-{beta} signaling networks we identified may contribute to disease pathogenesis and offer new avenues for intervention strategies.
- organism: mouse
- assay: scRNA-seq
- assay: bulk RNA-seq
- research focus: cell type annotation
- research focus: differential expression
- research focus: gene regulatory network inference
- research focus: benchmarking
- approach: statistical model
- approach: machine learning (classical)
- Citations
- 1
Reproducibility ledger
Nobody has reported trying to reproduce this paper yet. Did you replicate an experiment or run the code? Your report saves the next lab a week.
Sign in to file a reproducibility report.
Related papers
- Spatial Transcriptomic Profiling Reveals Microenvironment-Dependent Immune Signatures in a Lyme Arthritis model
Kaitlyn Gura, Mollie Hostetter, … Archana Dhasarathy·bioRxiv·Jun 25, 2026
- HBEGF+ macrophages identified in rheumatoid arthritis promote joint tissue invasiveness and are reshaped differentially by medications
D. Kuo, J. Ding, … L. T. Donlin·bioRxiv·Jan 31, 2019
- Secretory leukocyte protease inhibitor influences periarticular joint inflammation in B. burgdorferi-infected mice
Qian Yu, Xiao‐Tian Tang, … Erol Fikrig·bioRxiv·Nov 24, 2024
- Single cell ATAC sequencing identifies sleepy macrophages during reciprocity of cytokines in L.major infection
Shweta Khandibharad, Komal Kharat, Shailza Singh·Microbiology SpectrumPeer-reviewed·Feb 1, 2024
- Macrophage mediated recognition and clearance of Borrelia burgdorferi elicits MyD88-dependent and - independent phagosomal signals that contribute to phagocytosis and inflammation
S. J. Benjamin, K. L. Hawley, … J. C. Salazar·bioRxiv·Mar 30, 2019
- Co-infection of Ehrlichia with B. burgdorferi drives emergency myelopoiesis and promotes Lyme arthritis
J. L. Bonin, S. R. Torres, … K. C. MacNamara·bioRxiv·Aug 23, 2022