This code was used to analyze publicly available datasets from NCBI GEO: GSM5324643, GSM6146223 and GSM7508877. Cells from E12.5 and E13.5 wild-type craniofacial mesenchyme were subsetted for Csf1r expression and clustered by UMAP. A population of hemoglobin-expressing cells was identified and deemed to be contamination, and was removed from analysis, resulting in 3 distinct clusters. The expression of various macrophage and osteoclast markers and cytokines/chemokines by each population were characterized by violin plots.