Synovial macrophage diversity and activation of M-CSF signaling in post-traumatic osteoarthritis

  1. Alexander J Knights
  2. Easton C Farrell
  3. Olivia M Ellis
  4. Michelle J Song
  5. C Thomas Appleton
  6. Tristan Maerz  Is a corresponding author
  1. Department of Orthopaedic Surgery, University of Michigan, United States
  2. Department of Biomedical Engineering, University of Michigan, United States
  3. Department of Physiology and Pharmacology, Western University, Canada
  4. Bone and Joint Institute, Western University, Canada
  5. Department of Medicine, Schulich School of Medicine and Dentistry, Western University, Canada
  6. Department of Internal Medicine – Division of Rheumatology, University of Michigan, United States
6 figures, 1 table and 8 additional files

Figures

Figure 1 with 1 supplement
Diverse immune cell types are present in healthy and injured synovium.

(A) Uniform manifold approximation and projection (UMAP) plot of all immune cells by single-cell RNA sequencing (scRNA-seq) of synovium from mice subjected to Sham, 7 days anterior cruciate ligament …

Figure 1—figure supplement 1
Computational subsetting of immune cells from whole synovium single-cell RNA sequencing (scRNA-seq).

(A) Immune cell clusters (2, 5, 8) were computationally subset out of all synovial cells from previously published dataset available at NCBI GEO Accession number GSE211584 (Knights et al., 2023b). (B

Figure 2 with 5 supplements
Integration of joint immune cells from mouse arthritis single-cell RNA sequencing (scRNA-seq) datasets.

(A) Integrated uniform manifold approximation and projection (UMAP) plot of all immune cells from GSE134420 (Culemann), GSE200843 (Sebastian), GSE211584 (Knights), and GSE184609 (Muench). (B) Violin …

Figure 2—figure supplement 1
Recapitulation of published single-cell RNA sequencing (scRNA-seq) datasets: Culemann.

(A) Uniform manifold approximation and projection (UMAP) plot of all cell clusters from GSE134420 (Culemann et al., 2019). Cells were sorted CD45+ CD11b+ Ly6G- hindpaw synovial cells from mice …

Figure 2—figure supplement 2
Recapitulation of published single-cell RNA sequencing (scRNA-seq) datasets: Sebastian.

(A) Uniform manifold approximation and projection (UMAP) plot of all cell clusters from GSE200843 (Sebastian et al., 2022). Cells were sorted CD45+ cells obtained from a whole knee joint digestion …

Figure 2—figure supplement 3
Recapitulation of published single-cell RNA sequencing (scRNA-seq) datasets: Muench.

(A) Uniform manifold approximation and projection (UMAP) plot of all cell clusters from GSE184609 (Muench et al., 2022). Cells were sorted live cells obtained from hindpaw synovia of mice with …

Figure 2—figure supplement 4
Integration of immune single-cell RNA sequencing (scRNA-seq) datasets.

(A) Heatmap of top 5 gene markers for each cell cluster in the integrated immune object of all scRNA-seq datasets. (B) Total abundance of each major immune cell grouping for each dataset. (C) …

Figure 2—figure supplement 5
Assessment of neutrophils in the whole joint.

(A) Healthy mouse knee joints were subjected to two-step digestion described by Leale et al., 2022, resulting in a soft and hard tissue fraction. Joints were taken from naïve mice (unperfused) or …

Figure 3 with 2 supplements
Comparison of macrophages between arthritis disease states.

(A) Monocytes and MΦ from all post-traumatic osteoarthritis (PTOA) and rheumatoid arthritis (RA) immune cell datasets were computationally isolated, resulting in three clusters: MΦ, monocytes and …

Figure 3—figure supplement 1
Macrophages in osteoarthritis (OA) and rheumatoid arthritis (RA).

(A) Uniform manifold approximation and projection (UMAP) plots showing clusters of MΦ, monocytes, and osteoclast-like cells across datasets.

Figure 3—figure supplement 2
Outgoing signaling patterns from macrophages in post-traumatic osteoarthritis (PTOA) and rheumatoid arthritis (RA).

(A–B) River plots of outgoing communication patterns of MΦ, monocytes (Mono), and osteoclasts (Osteo) in PTOA (A) or RA (B), generated in CellChat. *Same signals from the same cell type in PTOA and …

Figure 4 with 1 supplement
Synovial macrophage subsets and trajectories in post-traumatic osteoarthritis (PTOA).

(A) Uniform manifold approximation and projection (UMAP) plot of monocytes and MΦ from synovium of Sham, 7 days anterior cruciate ligament rupture (ACLR) and 28 days ACLR mice, or split by condition …

Figure 4—figure supplement 1
Synovial macrophage subsets and trajectories.

(A) To generate an object containing only monocytes and MΦ from our dataset (GSE211584), we first subset all myeloid cell clusters from Seurat into Monocle3, then removed dendritic cell (DC), mast …

Figure 5 with 2 supplements
Stromal-immune crosstalk via M-CSF signaling.

(A) Cartoon depicting multi-directional and uni-directional crosstalk between stromal and immune cells in synovium, used for calculation of crosstalk communication probability scores. (B) Induction …

Figure 5—figure supplement 1
Stromal-immune crosstalk communication probability ratios in Sham and 7 days anterior cruciate ligament rupture (ACLR).

(A–C) Crosstalk communication probability ratios for multi-directional (stromal ←→ immune) or uni-directional (stromal → immune or immune → stromal) signaling in Sham synovium. A higher ratio …

Figure 5—figure supplement 2
M-CSF signaling hierarchy to stromal cells.

(A–C) CellChat hierarchy plots for the M-CSF signaling pathway in Sham, 7 days anterior cruciate ligament rupture (ACLR) and 28 days ACLR, with stromal cells as the receiving cell type in each case. …

Figure 6 with 1 supplement
Transcriptional control of monocyte differentiation in synovium.

(A–B) Pseudotime trajectory from monocytes to infiltrating MΦ showing directionality (arrowheads), starting point (root, star), branching points (nodes, black circles), and endpoints (termini, gray …

Figure 6—figure supplement 1
Expression of key transcription factors across pseudotime from monocytes to macrophages.

(A) RcisTarget analysis was performed to screen promoters of genes enriched in modules 1, 4, and 6 of the monocyte-to-infiltrating MΦ trajectory. Pseudotime regression plots for the top four …

Tables

Key resources table
Reagent type (species) or resourceDesignationSource or referenceIdentifiersAdditional information
Biological sample (Mus musculus)C57Bl/6 miceJackson LaboratoryStrain #000664Male and female
AntibodyTruStain FcX PLUSBioLegendRRID:AB_2783137Clone S17011E (1:1000)
AntibodyAnti-mouse CD3-APC/Fire750 (rat monoclonal)BioLegendRRID:AB_2572117Clone 17A2 (1:100)
AntibodyAnti-mouse CD11b-BV605 (rat monoclonal)BioLegendRRID:AB_11126744Clone M170 (1:400)
AntibodyAnti-mouse CD11c-PE/Dazzle594 (Armenian hamster monoclonal)BioLegendRRID:AB_2563654Clone N418 (1:200)
AntibodyAnti-mouse CD19-FITC (rat monoclonal)BioLegendRRID:AB_2629813Clone 1D3/CD19 (1:100)
AntibodyAnti-mouse CD45-BV650 (rat monoclonal)BioLegendRRID:AB_2565884Clone 30-F11 (1:400)
AntibodyAnti-mouse F4/80-APC/R700 (rat monoclonal)BD HorizonRRID:AB_2869711Clone T45-2342 (1:400)
AntibodyAnti-mouse FceRIa-PE/Cy7 (Armenian hamster monoclonal)BioLegendRRID:AB_10640122Clone MAR-1 (1:100)
AntibodyAnti-mouse Ly6G-BV421 (rat monoclonal)BioLegendRRID:AB_10897944Clone 1A8 (1:200)
AntibodyAnti-mouse MHC Class II-BV421 (rat monoclonal)BioLegendRRID:AB_2650896Clone M5/114.15.2 (1:200)
AntibodyAnti-mouse NK1.1-BV510 (mouse monoclonal)BioLegendRRID:AB_2562216Clone PK136 (1:200)
Peptide, recombinant proteinType IV CollagenaseSigmaCAS 9001-12-1
Peptide, recombinant proteinLiberaseSigma5401119001
Peptide, recombinant proteinDNaseISigmaCAS 9003-98-9
Chemical compound, drugTO-PRO-3 iodideInvitrogenT3605
Software, algorithmFlowJo v10BD/TreestarNA
Software, algorithmSeurat v4.1Hao et al., 2021https://doi.org/10.1016/j.cell.2021.04.048

Software, algorithmCluster Identity PRedictor (CIPR)Ekiz et al., 2020https://doi.org/10.1186/s12859-020-3538-2
Software, algorithmMonocle 3Cao et al., 2019https://doi.org/10.1038/s41586-019-0969-x
Software, algorithmRCisTarget (SCENIC)Aibar et al., 2017https://doi.org/10.1038/nmeth.4463
Software, algorithmCellChatJin et al., 2021

https://doi.org/10.1038/nmeth.4463
OtherscRNA-seq data of synovial cells from ACLRNCBI GEOGSE211584Knights et al., 2023a
OtherscRNA-seq data of synovial macrophages from RANCBI GEOGSE134420Culemann et al., 2019
OtherscRNA-seq data of joint immune cells from ACLRNCBI GEOGSE200843Sebastian et al., 2022
OtherscRNA-seq data of hindpaw joint cells from arthritisNCBI GEOGSE184609Muench et al., 2022
OtherBulk RNA-seq data of synoviumNCBI GEOGSE271903Bergman et al., 2024

Additional files

Supplementary file 1

Flow cytometry antibodies.

https://cdn.elifesciences.org/articles/93283/elife-93283-supp1-v1.docx
Supplementary file 2

Top 10 Cluster Identity PRedictor (CIPR) hits for each neutrophil cluster from the integrated object of all immune cells, all datasets.

https://cdn.elifesciences.org/articles/93283/elife-93283-supp2-v1.docx
Supplementary file 3

Ctrl vs post-traumatic osteoarthritis (PTOA) mac differentially expressed genes (DEGs).

https://cdn.elifesciences.org/articles/93283/elife-93283-supp3-v1.xlsx
Supplementary file 4

Ctrl vs rheumatoid arthritis (RA) mac differentially expressed genes (DEGs).

https://cdn.elifesciences.org/articles/93283/elife-93283-supp4-v1.xlsx
Supplementary file 5

Gene ontology outputs for post-traumatic osteoarthritis (PTOA) and rheumatoid arthritis (RA) macrophages.

https://cdn.elifesciences.org/articles/93283/elife-93283-supp5-v1.xlsx
Supplementary file 6

Gene ontology outputs for resident-like macrophages.

https://cdn.elifesciences.org/articles/93283/elife-93283-supp6-v1.xlsx
Supplementary file 7

Monocyte-to-infiltrating macrophage trajectory gene modules.

https://cdn.elifesciences.org/articles/93283/elife-93283-supp7-v1.xlsx
MDAR checklist
https://cdn.elifesciences.org/articles/93283/elife-93283-mdarchecklist1-v1.docx

Download links