Integrative Analysis of Single-cell and Bulk RNA Sequencing Reveals Macrophage Heterogeneity in Lyme Arthritis

Read the full article See related articles

Listed in

This article is not in any list yet, why not save it to one of your lists.
Log in to save this article

Abstract

Introduction

Lyme 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.

Methods

We 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.

Results

We 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-β 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.

Conclusion

Our 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-β signaling networks we identified may contribute to disease pathogenesis and offer new avenues for intervention strategies.

Article activity feed