Integrative Analysis of Single-cell and Bulk RNA Sequencing Reveals Macrophage Heterogeneity in Lyme Arthritis
Zhou, G.; Chen, X.; Dong, Y.; Luo, Y.; Liu, M.; Chen, Y.; Song, C.; Liu, X.; Yang, F.; Liang, Z.; Zi, L.; Gao, H.; Cao, Z.; Gao, F.
Show 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.
Matching journals
The top 6 journals account for 50% of the predicted probability mass.
Similar papers in this journal
- Transcriptome analysis of PBMCs reveals distinct immune response in the asymptomatic and re-detectable positive COVID-19 patients 94%
- A gene expression map of host immune response in human brucellosis 93%
- Cross-Tissue Transcriptomic Analysis Leveraging Machine Learning Approaches Identifies New Biomarkers for Rheumatoid Arthritis 93%
Similar papers in this journal
Similar papers in this journal
- Identification of a novel transcriptome signature for predicting the response to anti-TNF-α treatment in rheumatoid arthritis patients 92%
- Shotgun metagenomics reveals an enrichment of potentially cross-reactive bacterial epitopes in ankylosing spondylitis patients, as well as the effects of TNFi therapy and the host’s genotype upon microbiome composition 91%
- Risk factors for hospital admission related to COVID-19 in inflammatory rheumatic diseases 91%
Similar papers in this journal
- Ginkgolic acid attenuates e chinococcus granulosus infection-induced hepatic fibrosis by inhibiting Smad4 SUMOylation 93%
- Polarized Lung Inflammation and Tie2/Angiopoietin-Mediated Endothelial Dysfunction during Severe Orientia tsutsugamushi Infection 92%
- Annexin A2 depletion exacerbates the intracerebral microhemorrhage induced by acute rickettsia and Ebola virus infections 92%
Similar papers in this journal
- Multi-omics profiling reveals phenotypic and functional heterogeneity of neutrophils in COVID-19 94%
- Berberine delays onset of collagen induced arthritis through T cell suppression. 93%
- RIPK2 is crucial for the microglial inflammatory response to bacterial muramyl dipeptide but not to lipopolysaccharide 92%
"Similar papers" are the closest papers from that journal in the model's embedding space. They show what the match is built on, but the ranking comes mostly from a classifier over the whole training set, not from these examples alone.