Hypertrophic chondrocytes serve as a reservoir for marrow-associated skeletal stem and progenitor cells, osteoblasts, and adipocytes during skeletal development

  1. Jason T Long
  2. Abigail Leinroth
  3. Yihan Liao
  4. Yinshi Ren
  5. Anthony J Mirando
  6. Tuyet Nguyen
  7. Wendi Guo
  8. Deepika Sharma
  9. Douglas Rouse
  10. Colleen Wu
  11. Kathryn Song Eng Cheah
  12. Courtney M Karner
  13. Matthew J Hilton  Is a corresponding author
  1. Department of Cell Biology, Duke University School of Medicine, United States
  2. Department of Orthopaedic Surgery, Duke University School of Medicine, United States
  3. Department of Pharmacology and Cancer Biology, Duke University School of Medicine, United States
  4. Program of Developmental and Stem Cell Biology, Duke University School of Medicine, United States
  5. Division of Laboratory Animal Resources, Duke University School of Medicine, United States
  6. School of Biomedical Sciences, University of Hong Kong, Hong Kong
9 figures, 1 table and 2 additional files

Figures

Figure 1 with 1 supplement
Hypertrophic chondrocytes are the primary source of osteoblasts/osteocytes in trabecular bone.

Tibia sections of Col10a1CreERT2; Rosa26fs-tdTomato mice injected with tamoxifen at E13.5, 14.5, and 15.5 and subsequently sectioned at E16.5 and stained for (a) DAPI/RFP. (b–c) Higher magnification …

Figure 1—source data 1

Osteoblast and osteocyte quantification on Col10a1Cre;Rosa26fs-tdTomato 2M old mice.

https://cdn.elifesciences.org/articles/76932/elife-76932-fig1-data1-v2.xlsx
Figure 1—figure supplement 1
Contribution of hypertrophic chondrocytes to the bone marrow are observed into skeletal maturity and CRE expression is restricted to hypertrophic chondrocytes in Col10a1Cre mice.

(a) Tibia sections of Col10a1CreERT2; Rosa26fs-tdTomato mice injected with tamoxifen at 2 M of age and subsequently sectioned at 2-week chase. (b) Tibia sections of Col10a1CreERT2; Rosa26fs-tdTomato

Figure 2 with 2 supplements
Single cell transcriptomics at E16.5 captures hypertrophic chondrocytes and osteoblasts.

(a) UMAP shown in two-dimensional space produced using Seurat 4 package of R from single-cell RNA-sequencing of cleaned skeletal rudiments from Col10Cre; Rosa26fs-tdTomato mice at E16.5. (b–d) …

Figure 2—figure supplement 1
FACS plots for isolation of tdTOMATO+ cells at E16.5 and 2 M for scRNA-sequencing.

(a) E16.5 FACS plot identified tdTOMATO+ expressing cells isolated by total bone digestion and utilized for 10 X Genomics scRNA-sequencing represented in Figure 2. (b) 2 M FACS plot identified …

Figure 2—figure supplement 2
Additional genes of interest expressed in hypertrophic chondrocytes and descendants at E16.5.

(a) Feature plot of Cre expression was enriched in clusters 1–2 at E16.5. (b–c) Feature plot of genes associated with chondrocytes were observed enriched in clusters 1–7. (d-g) Feature plot of genes …

Figure 3 with 1 supplement
Single-cell transcriptomics of E16.5 hypertrophic chondrocytes and descendants reveals an intermediate SSPC upstream of osteoblasts.

(a–d) Feature plot of SSPC-associated genes identified between hypertrophic chondrocyte and osteoblast clusters (e) Violin plot representing the relative level of SSPC associated gene expression (a–d

Figure 3—source data 1

PDGFRa+ TOMATO+ quantification at e16.5 of Col10a1Cre; Rosa26-fs-tdTomato; PDGFRa-H2B-GFP.

https://cdn.elifesciences.org/articles/76932/elife-76932-fig3-data1-v2.xlsx
Figure 3—figure supplement 1
SSPC associated proteins are detected in hypertrophic chondrocyte-derived marrow-associated cells.

(a) Immunostaining of SSPC associated proteins in both bone sections from Col10a1CreERT2; Rosa26fs-tdTomato mice injected with tamoxifen at P6 and chased for 2 weeks and (b) 2 M old Col10a1Cre; …

Figure 3—figure supplement 1—source data 1

LEPR+ TOMATO+ quantification of Col10a1Cre;Rosa26fs-tdTomato at 2M.

https://cdn.elifesciences.org/articles/76932/elife-76932-fig3-figsupp1-data1-v2.xlsx
Figure 3—figure supplement 1—source data 2

Flow cytometric analysis of PDGFRa-H2B-GFP and tdTOMATO on Col10a1Cre;Rosa26fs-tdTomato at 2M.

https://cdn.elifesciences.org/articles/76932/elife-76932-fig3-figsupp1-data2-v2.xlsx
Figure 3—figure supplement 1—source data 3

Flow cytometric analysis of PDGFRb and tdTOMATO on Col10a1Cre;Rosa26fs-tdTomato at 2M.

https://cdn.elifesciences.org/articles/76932/elife-76932-fig3-figsupp1-data3-v2.xlsx
Figure 4 with 4 supplements
Single cell transcriptomics of hypertrophic chondrocyte descendants following FACS and 10 X Genomics sequencing at 2 M of age.

(a) UMAP shown in two-dimensional space produced using Seurat 3 package of R from single-cell RNA-sequencing of bone marrow digest from Col10Cre; Rosa26fs-tdTomato mice at 2 M. (b–f) Feature plots …

Figure 4—figure supplement 1
Additional osteoblast associated genes observed in clusters 1 and 2 at 2 months of age.

(a–g) Feature plots of genes associated with osteoblasts were observed in clusters 1 and 2 with minimal expression in clusters 3–5. (h) Violin plots representing relative gene expression levels in (a…

Figure 4—figure supplement 2
Differentially expressed genes mostly unique to clusters 3, 4, and 5.

(a–j) Feature plots of genes identified by differential gene expression in cluster 3 (a–d), cluster 4 (e–h), and cluster 5 (i–l). (m) Violin plots of differentially expressed genes represented in (a–…

Figure 4—figure supplement 3
Hypertrophic chondrocyte-derived cells can associate with blood vessels; however, do not express genes associated with pericytes or vascular smooth muscle cells.

(a) Immunostaining for endothelial-cell-associated protein, ENDOMUCIN, on 2-month-old Col10a1Cre; Rosa26fs-tdTomato bone sections show adjacent tdTOMATO+ cells indicated by orange arrows. Feature …

Figure 4—figure supplement 4
LeprCre; R26-tdTomato and AdipoqCre; R26-tdTomato mice exhibit reporter expression in marrow associated cells, osteoblasts, and adipocytes with age.

(a) Immunostaining of bone sections from LeprCre; Rosa26fs-tdTomato mice exhibit tdTOMATO+ marrow-associated cells; however, do not display tdTOMATO+, OCN+ osteoblasts at P0. Immunostaining at 2 …

Figure 5 with 1 supplement
Many hypertrophic chondrocyte derived cells express genes associated with SSPCs.

(a–f) Feature plots of genes previously identified as SSPC markers in genetic mouse models. (g–i) Feature plots of genes previously identified as SSPC markers for use in flow cytometry and FACS. (j) …

Figure 5—figure supplement 1
Contribution of hypertrophic chondrocyte descendants to total SSPC populations.

(a) Quantitation of immunostaining for LEPR on bone sections from 2 M old Col10a1Cre; Rosa26fs-tdTomato mice exhibit 37.4% LEPR+; tdTOMATO+ as compared to total LEPR+ marrow-associated cells (Figure …

Figure 5—figure supplement 1—source data 1

LEPR+ TOMATO+ quantification on Col10a1Cre;Rosa26fs-tdTomato at 2M.

https://cdn.elifesciences.org/articles/76932/elife-76932-fig5-figsupp1-data1-v2.xlsx
Figure 5—figure supplement 1—source data 2

Flow cytometric analysis of PDGFRa-H2B-GFP and TOMATO on Col10a1Cre;Rosa26fs-tdTomato at 2M.

https://cdn.elifesciences.org/articles/76932/elife-76932-fig5-figsupp1-data2-v2.xlsx
Hypertrophic chondrocyte derived SSPCs exhibit osteogenic and adipogenic differentiation capacities.

(a) Feature plot and violin plot of the osteoblast specification gene, Runx2. (b–c) Feature plots and violin plots of the adipogenic specification gene, Pparg, and adipogenic associated gene, …

Figure 6—source data 1

PERILIPIN+ TOMATO+ quantifications on Col10a1Cre;Rosa26fs-tdTomato at 2M.

https://cdn.elifesciences.org/articles/76932/elife-76932-fig6-data1-v2.xlsx
CFU-Fs derived from hypertrophic chondrocytes are similar to CFU-Fs derived from other cell sources, but contain SSPCs with enhanced proliferative capacities.

(a–c) CFU-Fs derived from hypertrophic chondrocytes (tdTOMATO+) (a) and those derived from other cell sources (tdTOMATO-)(b) are established and develop at similar frequencies at 2- and 8 months of …

Figure 7—source data 1

Tomato+ colony unit formation quantifications at 2 and 8M of Col10a1Cre;Rosa26fs-tdTomato.

https://cdn.elifesciences.org/articles/76932/elife-76932-fig7-data1-v2.xlsx
Figure 7—source data 2

Bulk RNA-sequencing total gene list of 3 TOMATO+ and 3 TOMATO- colonies from Col10a1Cre;Rosa26fs-tdTomato.

https://cdn.elifesciences.org/articles/76932/elife-76932-fig7-data2-v2.xlsx
Figure 7—source data 3

Gene ontology analysis of genes enriched in TOMATO+ colonies.

https://cdn.elifesciences.org/articles/76932/elife-76932-fig7-data3-v2.xlsx
Figure 7—source data 4

Gene ontology analysis of genes enriched in TOMATO- colonies.

https://cdn.elifesciences.org/articles/76932/elife-76932-fig7-data4-v2.xlsx
Figure 7—source data 5

EDU+ quantifications of TOMATO+ and TOMATO- colonies.

https://cdn.elifesciences.org/articles/76932/elife-76932-fig7-data5-v2.xlsx
Figure 7—source data 6

Colony size of TOMATO+ and TOMATO- colonies.

https://cdn.elifesciences.org/articles/76932/elife-76932-fig7-data6-v2.xlsx
Figure 8 with 1 supplement
Kidney capsule transplantations of tdTOMATO+ cells exhibit complete ossicle formation with bone, adipocytes, and bone marrow compared to tdTOMATO- resulting in only bone formation.

(a) Immunofluroescent stain for bone (osteopontin), adipocytes (perilipin), and vessels (endomucin) for tomato+ transplant (top) and tomato- transplant (bottom). (b) Hematoxylin and Eosin stain for …

Figure 8—figure supplement 1
In vitro differentiation of colonies reveals relatively similar osteogenic and adipogenic differentiation capacities.

(a–b) CFUs from Col10a1Cre; Rosa26fs-tdTomato bone marrow differentiate into osteoblasts stained with alkaline phosphatase (a) and von kossa (b) (CFU-OB). alkaline phosphatase: SD tdTOMATO+ (T) ± …

Figure 8—figure supplement 1—source data 1

Differentiation assay quantifications of TOMATO+ and TOMATO- colonies.

https://cdn.elifesciences.org/articles/76932/elife-76932-fig8-figsupp1-data1-v2.xlsx
Revised model for the generation of trabecular/endocortical bone associated osteoblasts/osteocytes from multiple sources.

(a) Perichondrial or periosteal SSPCs/osteoprogenitors (pSSPC/pOP) migrate into the marrow utilizing blood vessels and possess the ability to further differentiate into osteoblasts/osteocytes …

Tables

Appendix 1—key resources table
Reagent type (species) or resourceDesignationSource or referenceIdentifiersAdditional information
Strain, strain background (Mus musculus)C57BL/6 JJackson LaboratoryRRID:IMSR_JAX:000664
Genetic reagent (Mus musculus)Col10a1CreDOI: 10.1242/bio.201411031Dr. Klaus von der Mark
Genetic reagent (Mus musculus)Col10a1CreERT2This paperDr. Kathryn Cheah
Genetic reagent (Mus musculus)Rosa26fs-tdTomatoJackson LaboratoryRRID:IMSR_JAX:007909
Genetic reagent (Mus musculus)PdgfraH2B-GFPJackson LaboratoryRRID:IMSR_JAX:007669
Genetic reagent (Mus musculus)LeprCreJackson LaboratoryRRID:IMSR_JAX:008320
Genetic reagent (Mus musculus)AdipoqCreJackson LaboratoryRRID:IMSR_JAX:028020
Sequence-based reagentCol10a1Cre/
ERT2_F
DOI: 10.1242/bio.201411031TTTAGAGCATTATT
TCAAGGCA
GTTTCCADr. Klaus von der Mark
Sequence-based reagentCol10a1Cre/
ERT2_R
DOI: 10.1242/bio.201411031AGGCAAATCTT
GGTGTACGGDr. Klaus von der Mark
AntibodyENDOMUCIN (rat mAb)AbcamCat# ab106100, RRID:AB_10859306IF (1:100), 0.3% Triton X-100 in 1 x PBS
AntibodyPERILIPIN (rabbit mAb)Cell SignalingCat# 9349, RRID:AB_10829911IF (1:100)
AntibodyOSTEOCALCIN (rabbit pAb)MilliporeCat# AB10911, RRID:AB_1587337IF (1:200), 10 ug/mL Proteinase K
AntibodyOSTERIX(Rabbit pAb)AbcamCat# ab22552, RRID:AB_2194492IF (1:400)
AntibodyCRE(Rabbit mAb)Cell SignalingCat# 15036, RRID:AB_2798694IF (1:100), 0.3% Triton X-100 in 1 x PBS
AntibodyRFP (rabbit pAb)AbcamCat# ab62341, RRID:AB_945213IF (1:100)
AntibodyLEPR (goat pAb)R&D SystemsCat# AF497, RRID:AB_2281270IF (1:50), 0.3% Triton X-100 in 1 x PBS
AntibodyPDGFRa (goat mAb)R&D SystemsCat# AF1062, RRID:AB_2236897IF (1:50), 0.3% Triton X-100 in 1 x PBS
AntibodyPDGFRb (rat mAb)InvitrogenCat# 14-1402-81, RRID:AB_467492IF (1:50), 0.3% Triton X-100 in 1 x PBS
AntibodyOsteopontin (goat pAb)R&D SystemsCat# AF808, RRID:AB_2194992IF (1:100), 0.3% Triton X-100 in 1 x PBS
AntibodyPDGFRb (rat mAb)InvitrogenCat# 17-1402-80, RRID:AB_1548752Flow Cytometry (1:25)
AntibodyCD45-FITC (rat mAb)BiolegendCat# 103107, RRID:AB_312972Flow Cytometry (1:1000)
AntibodyCD45-APC (rat mAb)BiolegendCat# 103111, RRID:AB_312976Flow Cytometry (1:1000)
AntibodyCD31-FITC (rat mAb)BiolegendCat# 102405, RRID:AB_312900Flow Cytometry (1:1000)
AntibodyCD31-APC (rat mAb)BiolegendCat# 102409, RRID:AB_312904Flow Cytometry (1:1000)
AntibodyTER119-FITC (rat mAb)BD BiosciencesCat# 561032, RRID:AB_10563083Flow Cytometry
(1:500)
AntibodyTER119-APC (rat mAb)BiolegendCat# 116211, RRID:AB_313712Flow Cytometry
(1:250)
AntibodyDAPI (FLOW)ThermoFisherCat# D1306, RRID:AB_2629482Flow Cytometry (1:1000)
AntibodyGoat anti-Mouse Alexa Fluor 488 pAbInvitrogenCat# A-11001, RRID:AB_2534069IF (1:200)
AntibodyDonkey anti-Goat Alexa Fluor 488 pAbInvitrogenCat# A-11055, RRID:AB_2534102IF (1:200)
AntibodyGoat anti-Rat Alexa Fluor 488 pAbInvitrogenCat# A-11006, RRID:AB_2534074IF (1:200)
AntibodyGoat anti-Rabbit Alexa Fluor 594 pAbInvitrogenCat# A-11037, RRID:AB_2534095IF (1:200)
AntibodyGoat anti-Rabbit Alexa Fluor 488 pAbInvitrogenCat# A-11034, RRID:AB_2576217IF (1:200)
AntibodyGoat anti-Rat Alexa Fluor 647 pAbInvitrogenCat# A-21247, RRID:AB_141778IF (1:200)
AntibodyGoat anti-Rabbit Alexa Fluor 488 pAbInvitrogenCat# A-11008, RRID:AB_143165For LepRCre;tdTomato
stained PERILIPIN only -
IF (1:200)
Chemical compound, drugCollagenase IIGibco17101–015
Chemical compound, drugCollagenase DRoche11088866001
Chemical compound, drugPolyvinylpyrrolidoneSigmaP5288
Chemical compound, drugGlutamax 100 xGibco35050–061
Chemical compound, drugL-Ascorbic AcidSigmaA4544
Chemical compound, drugB-glycerophosphateSigmaG9422
Chemical compound, drugInsulin HumanSigmaI3536
Chemical compound, drugIBMXSigmaI7018
Chemical compound, drugDexamethasoneSigmaD2915
Chemical compound, drugAntigen Unmasking solution, Citrate Acid BasedVectorH3300
Commercial assay, kitClick-iT Plus EdU Cell
Proliferation Kit
InvitrogenC10637
Software, algorithmFIJIdoi:10.1038/
nmeth.2019
Fiji, RRID:SCR_002285
Software, algorithmSeuratCRANSeurat, RRID:SCR_007322
Software, algorithmMonocle3GithubMonocle3, RRID:SCR_018685
Software, algorithmTrim GaloreBabraham BioinformaticsTrim Galore,RRID:SCR_011847Felix Krueger
Software, algorithemCutadaptGithubcutadapt, RRID:SCR_011841
Software, algorithmSTAR RNA-seq AlignmentGithubSTAR, RRID:SCR_004463
Software, algorithmHTSeq toolGithubHTSeq, RRID:SCR_005514
Software, algorithmDESeq2BioconductorDESeq2, RRID:SCR_015687
Software, algorithmGene OntologyGeneontoloy.orgGene Ontology, RRID:SCR_002811
OtherSilver Nitrate 10%Ricca6830–4
OtherOil Red OSigmaO0625
OtherAlizarin RedSigmaA5533
OtherMayer’s
Hematoxylin
solution
Electron
Microscopy
Sciences
26043–06
OtherEosin Y with Phloxine
staining solution
Electron
Microscopy
Sciences
26051–21

Additional files

Transparent reporting form
https://cdn.elifesciences.org/articles/76932/elife-76932-transrepform1-v2.pdf
Source code 1

Code used for generation of Single Cell datasets and analysis at e16.5/2M of Col10a1Cre;Rosa26fs-tdTomato.

https://cdn.elifesciences.org/articles/76932/elife-76932-code1-v2.zip

Download links