Hes1 marks peri-condensation mesenchymal cells that generate both chondrocytes and perichondrial cells in early bone development

Bone development starts with condensations of undifferentiated mesenchymal cells that set a framework for future bones within the primordium. In the endochondral pathway, mesenchymal cells inside the condensation differentiate into chondrocytes and perichondrial cells in a SOX9-dependent mechanism. However, the identity of mesenchymal cells outside the condensation and how they participate in developing bones remain undefined. Here we show that mesenchymal cells surrounding the condensation contribute to both cartilage and perichondrium, robustly generating chondrocytes, osteoblasts, and marrow stromal cells in developing bones. Single-cell RNA-seq analysis of Prrx1-cre-marked limb bud mesenchymal cells at E11.5 reveals that Notch effector Hes1 is expressed in a mutually exclusive manner with Sox9 that is expressed in pre-cartilaginous condensations. Analysis of a Notch signaling reporter CBF1:H2B-Venus reveals that peri-condensation mesenchymal cells are active for Notch signaling. In vivo lineage-tracing analysis using Hes1-creER identifies that Hes1+ early mesenchymal cells surrounding the SOX9+ condensation at E10.5 contribute to both cartilage and perichondrium at E13.5, subsequently becoming growth plate chondrocytes, osteoblasts of trabecular and cortical bones, and marrow stromal cells in postnatal bones. In contrast, Hes1+ cells in the perichondrium at E12.5 or E14.5 do not generate chondrocytes within cartilage, contributing to osteoblasts and marrow stromal cells only through the perichondrial route. Therefore, Hes1+ peri-condensation mesenchymal cells give rise to cells of the skeletal lineage through cartilage-dependent and independent pathways, supporting the theory that early mesenchymal cells outside the condensation also play important roles in early bone development.

Bone development starts with condensations of undifferentiated mesenchymal cells that set a framework for future bones within the primordium. In the endochondral pathway, mesenchymal cells inside the condensation differentiate into chondrocytes and perichondrial cells in a SOX9-dependent mechanism. However, the identity of mesenchymal cells outside the condensation and how they participate in developing bones remain undefined. Here we show that mesenchymal cells surrounding the condensation contribute to both cartilage and perichondrium, robustly generating chondrocytes, osteoblasts, and marrow stromal cells in developing bones. Single-cell RNA-seq analysis of Prrx1-cre-marked limb bud mesenchymal cells at E11.5 reveals that Notch effector Hes1 is expressed in a mutually exclusive manner with Sox9 that is expressed in pre-cartilaginous condensations. Analysis of a Notch signaling reporter CBF1:H2B-Venus reveals that peri-condensation mesenchymal cells are active for Notch signaling. In vivo lineage-tracing analysis using Hes1-creER identifies that Hes1 + early mesenchymal cells surrounding the SOX9 + condensation at E10.5 contribute to both cartilage and perichondrium at E13.5, subsequently becoming growth plate chondrocytes, osteoblasts of trabecular and cortical bones, and marrow stromal cells in postnatal bones. In contrast, Hes1 + cells in the perichondrium at E12.5 or E14.5 do not generate chondrocytes within cartilage, contributing to osteoblasts and marrow stromal cells only through the perichondrial route. Therefore, Hes1 + peri-condensation mesenchymal cells give rise to cells of the skeletal lineage through cartilage-dependent and independent pathways, supporting the theory that early mesenchymal cells outside the condensation also play important roles in early bone development.
Bone development involves highly sequential steps of mesenchymal cell proliferation and differentiation, thus represents a prime example of organogenesis deliberately executed through heterotypic cellular interactions. Condensation of undifferentiated mesenchymal cells is the first step of bone development, which sets a framework for future bones within the primordium (1,2). In this process, a previously dispersed population of undifferentiated mesenchymal cells migrate to the site of the future skeleton in the limb bud and aggregate to form a vasculature-free element termed mesenchymal condensation (3). Mesenchymal condensations serve as an important modality to generate skeletal precursor cells on a large scale in a highly organized environment. Particularly in the endochondral pathway by which most of mammalian bones are formed, mesenchymal condensations give rise to both cartilage template and perichondrium in the subsequent phase. Chondrocytes within the cartilage template and perichondrial cells provide distinct cellular sources of growth plate chondrocytes, trabecular and cortical osteoblasts, and marrow stromal cells (4)(5)(6)(7). Therefore, mesenchymal condensations provide the foundation of developing bones.
The current concept holds that SOX9, a master transcription factor for chondrogenesis, is essential for condensing mesenchymal cells to differentiate into chondrocytes and perichondrial cells (1,(8)(9)(10). Indeed, SOX9 is absolutely required for these mesenchymal cells to stay organized within the condensation (11). Fate-mapping studies using Sox9-cre/creER demonstrate that Sox9 + cells within the condensation function as osteo-chondro-progenitor cells as they give rise to chondrocytes and perichondrial cells in the subsequent stage (8,9,12). In the limb, Prrx1-cre is used to mark the entire spectrum of cells of the skeletal lineage including "skeletal stem cells" and their derivatives, as Prrx1 is expressed not only by condensing mesenchymal cells but also by other cells that originate from the lateral plate mesoderm (13,14). Interestingly, Prrx1 is predominantly expressed in the perichondrium at E13.5 when the cartilage template is established (4,15,16), suggesting that cells surrounding the condensation may also robustly contribute to the generation of skeletal lineage cells. However, the identity of early mesenchymal cells outside the condensation and how they participate in developing bones remain undefined.
We recently demonstrated that Dlx5 + early perichondrial cells located in the outer layer of the perichondrium at E12.5 sustainably contribute to cortical bone and marrow stromal compartments in developing bones (4). These Dlx5 + perichondrial cells remain in the perichondrium in subsequent stages, unlike osterix (Osx + ) osteogenic perichodrial cells that are in the inner layer of the perichondrium (17). Osx + cell descendants contribute to skeletogenesis only transiently and eventually disappear from the perichondrium and the skeletal element (8,18). Interestingly, however, Dlx5 + early perichondrial cells do not contribute to cartilage. It has not been determined if peri-condensation mesenchymal cells preceding Dlx5 + perichondrial cells can generate chondrocytes within the cartilage template.
In this study, we hypothesize that early mesenchymal cells surrounding the condensation can contribute to chondrocytes of the cartilage template while also remaining outside the condensation and contributing to perichondrial cells. To test this hypothesis, we utilized in vivo lineage-tracing approaches to define the cell fates of peri-condensation mesenchymal cells. Using single-cell RNA-seq analyses, we identified Notch effector Hes1 as a potential marker of peri-condensation mesenchymal cells that are located outside the SOX9 + condensation and subsequently revealed their cell fates using a tamoxifen-inducible Hes1-creER allele. Our findings demonstrate that Hes1 + peri-condensation mesenchymal cells give rise to cells of the skeletal cell lineage through cartilage-dependent and independent pathways in growing bones and critically support early bone development.
In search for a marker for peri-condensation mesenchymal cells, we noticed that a Notch effector gene Hes1 exhibited a pattern negatively correlated with Sox9, identified as cell-type specific markers for Cluster 3,5,6 ( Fig. 1E). Notch signaling is essential for maintaining skeletal progenitor cells, and its target gene Hes1 is most abundantly expressed in the fetal perichondrium (22,23). RNAscope analysis revealed that Hes1 was expressed in a manner surrounding the Sox9-expressing condensation (Fig. 1F). To further validate Notch-responsive status of peri-condensation mesenchymal cells in situ, we utilized a Notch signaling reporter strain expressing a histone 2B (H2B)-bound Venus protein under a C promoter binding factor 1 (CBF1) promoter (CBF1:H2B-Venus) (24). Interestingly, Notch-responsive Venus bright cells were mostly located outside the SOX9 + domain of the condensation (Fig. S1B). Therefore, Notch signaling is active in mesenchymal cells surrounding the SOX9 + condensation. In subsequent in vivo lineage-tracing studies, we focused on Notch effector Hes1 as a marker of peri-condensation mesenchymal cells.

Hes1-creER + peri-condensation mesenchymal cells can generate chondrocytes
To study the cell fates of these Hes1 + cells surrounding Sox9 + condensation, we performed lineage-tracing experiments using a Hes1-creER knock-in allele (25) that activates an R26R-tdTomato reporter in a tamoxifen-dependent manner. First, we pulsed Hes1-creER; R26R tdTomato mice at E10.5 and analyzed these mice after 24 h at E11.5 to define the identities of Hes1-creER + cells (Hes1 CE -E10.5 cells). Hes1 CE -E10.5 cells were located throughout the limb bud at E11.5 in a manner excluding the SOX9 + pre-cartilaginous condensation ( Fig. 2A). Hes1 CE -E10.5 cells were also located in the area immediately adjacent to SOX9 + cells associated with CBF1:H2B-Venus reporter activities (Fig. 2B). We also noticed that Hes1-creER simultaneously marked other types of cells outside the skeletal elements of the limb bud, including myosin heavy chain 3 (MYH3 + ) skeletal muscle cells and endomucin (EMCN + ) endothelial cells (Fig. 2, C and D). Therefore, Hes1-creER can mark mesenchymal cells in the SOX9-negative domain of the condensation upon tamoxifen injection.
Subsequently, we traced the fate of Hes1-creER + cells in endochondral bone development (see Fig. S2, A and B for H&E and Alcian Blue staining of developing endochondral bones). After 3 days of chase at E13.5, Hes1 CE -E10.5 cells contributed to both chondrocytes and perichondrial cells of the cartilage template at E13.5 (Fig. 2E). Hes1 CE -E10.5 cells became SOX9 + chondrocytes within the cartilage template while contributing to a majority of perichondrial cells including those expressing CBF1:H2B-Venus (Fig. 2E, right panel). In fact, Hes1 CE -E10.5 cells contributed to all layers of the perichondrium including OSX + osteoblast precursors, in addition to MYH3 + skeletal muscle cells outside the perichondrium and EMCN + endothelial cells (Figs. 2F and S3A). After 9 days of the chase at postnatal day (P) 0, Hes1 CE -E10.5 cells contributed to many columnar chondrocytes of the growth plate, Col1a1-GFP + osteoblasts on the cortical and trabecular bone, and stromal cells throughout the marrow space (Fig. 2G). We also examined the cell fate of Hes1 + cells at an earlier stage. For this purpose, we pulsed Hes1-creER; R26R tdTomato mice at E8.5 and analyzed these mice at E13.5 and E18.5. Hes1 CE -E8.5 cells contributed to essentially all skeletal cells, robustly generating chondrocytes, perichondrial cells, osteoblasts, and marrow stromal cells (Fig. S3, C and D). Notably, Hes1 CE -E8.5 cells contributed more substantially to the skeletal lineage than Hes1 CE -E10.5 cells did, indicating that Hes1 + cells at a precondensation stage might possess a strong chondrogenic potential.
Importantly, skeletal muscle cells or endothelial cells that were initially marked by Hes1-creER at E10.5 did not contribute to these skeletal lineage cells, as cells marked by skeletal muscle-specific cre lines (Acta1-cre, Myl1-cre, and Mck-cre) or an endothelial cell-specific cre line (Tie2-cre) did not contribute to chondrocytes, osteoblasts, or marrow stromal cells at any time points (Fig. S4, A and B). This allows us to exclude the contribution of skeletal muscle cells and endothelial cells initially labelled by Hes1-creER at E10.5 to skeletal lineage cells.
Therefore, Hes1 + mesenchymal cells surrounding the condensation can differentiate not only into perichondrial cells but also into chondrocytes within the cartilage template, both of which robustly contribute to chondrocytes, osteoblasts, and marrow stromal cells in the later stages of endochondral bone development (Fig. 2H).

Hes1-creER + perichondrial cells possess osteogenic but not chondrogenic potential
Subsequently, we sought to define the cell fate of Hes1 + cells in the later stages of endochondral bone development, at E12.5 when the perichondrium is established adjacent to the cartilage template, as well as at E14.5 when the formation of the marrow space starts due to vascular invasion into the cartilage template.

Hes1 + early mesenchymal cells in endochondral pathway
We subsequently traced the fate of Hes1 CE -E12.5 perichondrial cells during the formation of the primary ossification center. After 3 days of chase at E15.5, Hes1 CE -E12.5 cells contributed to cells within the primary ossification center as well as Col1a1-GFP + osteoblasts in the bone collar (Fig. 3D). After 6 days of chase at E18.5, Hes1 CE -E12.5 cells contributed to Cxcl12-GFP high marrow stromal cells and Col1a1-GFP + osteoblasts in the trabecular bone (Fig. 3, E and F). Importantly, unlike those pulsed at an earlier stage, Hes1 CE -E12.5 cells did not contribute to SOX9 + chondrocytes in the growth plate but maintained themselves within the perichondrium, and robustly contributed to the bone collar and Col1a1-GFP + osteoblasts in the cortical bone (Fig. 3E). Therefore, Hes1-creER marks perichondrial cells with robust capability to contribute to cortical and trabecular bone compartments as well as to marrow stromal compartments, akin to those marked by Dlx5-creER that we recently reported (4).
We also traced the fate of Hes1-creER + cells at E14.5 (Hes1 CE -E14.5 cells). At E15.5, Hes1 CE -E14.5 cells were located among the perichondrium, the periosteum, and the primary ossification center, although most of these cells appeared to be non-skeletal with minimal overlap with Col1a1-GFP + osteoblasts (Fig. 3G). After 5 days of chase at P0, Hes1 CE -E14.5 cells did not robustly contribute to Col1a1-GFP + osteoblasts in the cortical and trabecular bones (Fig. 3H), indicating that Hes1 + cells lose robust osteogenic potential at E14.5. Postnatally, Hes1-creER + cells at P3 (Hes1 CE -P3 cells) were localized in the groove of Ranvier, the metaphyseal primary spongiosa and the superficial layer of the articular cartilage at P5 (Fig. S5). Hes1 CE -P3 cells did not overtly overlap with SOX9 + cells within the growth plate, whereas a few Hes1 CE -P3 cells in a deeper layer of the articular cartilage expressed SOX9. Therefore, Hes1 + cells do not overlap with SOX9 + cells at a later stage of perichondrial development.
Second, we quantified the contribution of Hes1 CE -E10.5, Hes1 CE -E12.5 or Hes1 CE -E14.5 cells to the E15.5 perichondrium and bone collar, among outer Col1a1-GFP neg perichondrial cells (top) or inner Col1a1-GFP + osteogenic bone collar cells (bottom). Hes1-creER did not effectively mark perichondrial cells at E14.5, as Hes1 CE -E14.5 cells contributed to a significantly fewer fraction of outer perichondrial cells than Hes1 CE -E10.5 or Hes1 CE -E12.5 cells did (Fig. 3I, second panel from the left in the upper panels). Hes1 CE -E10.5 cells contributed to a majority of Col1a1-GFP + osteogenic cells in the bone collar. However, the osteogenic potential of Hes1 + cells declined in later stages, as Hes1 CE -E12.5 and Hes1 CE -E14.5 cells contributed to progressively fewer fractions of inner Col1a1-GFP + osteogenic cells in the bone collar (Fig. 3I, second panel from the left in the lower panels). These trends continued onto P0, as the contribution of Hes1 + cells to cortical and trabecular bone osteoblasts progressively declined in later stages (Fig. 3I, third panel from the left in the lower panels).
Third, we quantified the contribution of Hes1 CE -E10.5, Hes1 CE -E12.5, or Hes1 CE -E14.5 cells to SOX9 + chondrocytes within the cartilage template. The chondrogenic potential was unique to Hes1 CE -E10.5 cells. These cells contributed progressively to a larger fraction of SOX9 + chondrocytes within the growth plate during the chase (Fig. 3I, upper right panel).
Therefore, Hes1 CE -E10.5 peri-condensation mesenchymal cells have both chondrogenic and osteogenic potentials, while Hes1 CE -E12.5 perichondrial cells have only osteogenic potentials. Importantly, Hes1 + descendants appeared to continue expressing Hes1 in the perichondrium, as the relative fraction of Hes1 CE -E10.5 and Hes1 CE -E12.5 cells in the perichondrium remained comparable at E13.5. Interestingly Hes1 CE -E14.5 cells lose their robust potential to contribute neither chondrocyte or osteoblast, indicating that the ability of Hes1 + cells to generate chondrocytes and osteoblasts is restricted to an early stage of endochondral bone development.

Discussion
Here, we demonstrate that a population of early mesenchymal cells defined by Notch effector Hes1 expression outside the SOX9 + condensation can provide an important source of chondrocytes and perichondrial cells in early endochondral bone development. These Hes1 + peri-condensation mesenchymal cells represent precursors of Dlx5 + early perichondrial cells that we reported recently (4) and robustly participate in chondrogenesis, osteogenesis, and subsequently the formation of the postnatal bone marrow stromal compartment through two routes of cartilage-dependent and independent pathways. Our Hes1-creER-based in vivo lineage-tracing findings establish the concept that SOX9-negative Hes1-expressing mesenchymal cells surrounding mesenchymal condensations play important roles in early endochondral bone development. Importantly, the robust chondrogenic and osteogenic potential of Hes1 + cells is unique to the mesenchymal condensation stage, whereas these cells progressively lose such potential in later stages.
We believe that the two lineages of chondrocytes and perichondrial cells diverge prior to E12.5, as Hes1 CE -E10.5 cells, but not Hes1 CE -E12.5 cells, contribute to chondrocytes within the cartilage template. Hes1 + cells at E10.5 and earlier time are likely to provide a common progenitor population for at least some of the chondrocytes and osteoblasts. We acknowledge, however, that our approach based on a single-color reporter cannot determine if an individual peri-condensation Hes1 + cell can contribute clonally both to chondrocytes and osteoblasts. An in vivo clonal analysis with a multi-color lineage reporter would facilitate the identification of a common progenitor population for chondrocytes and osteoblasts in the pericondensation region.
Hes1 is one of the canonical Notch effector genes that is expressed by cells in which Notch signaling is activated. In fact, we observed that a Notch signaling reporter CBF1:H2B-Venus was specifically active in peri-condensation mesenchymal cells surrounding the SOX9 + pre-cartilaginous condensation. Notch signaling, which is primarily mediated by the CBF1dependent target gene Hes1, regulates both Sox9 and Runx2 expression during skeletal cell differentiation (22). Hes1 inhibits Sox9 expression, and Hes1 inhibition promotes chondrogenesis by upregulating Sox9. Hes1 also inhibits osteoblast differentiation through the inhibition of Runx2 activities (23). Notch signaling regulates asymmetric cell division and lineage decision through lateral inhibition. It is possible that Hes1 + peri-condensation mesenchymal cells divide into Hes1positive and negative cells through asymmetric cell division; the former Hes1-positive cells may stay in the perichondrium as undifferentiated skeletal progenitors, while the latter Hes1negative cells enter into the cartilage template and differentiate into chondrocytes.
We show that Hes1-creER allows the marking of early mesenchymal cells in a SOX9-negative domain of the condensation. Historically, Prrx1-cre/creER has been widely Figure 5. Early developmental Hes1-creER + cell-derived marrow stromal cells behave as skeletal progenitor cells. A, Colony-forming unit fibroblast (CFU-F) assay of Hes1-creER; R26R tdTomato at P21, pulsed at E10.5 (top), E12.5 (middle), and E14.5 (bottom). Scale bar: 5 mm. B, Percentage of tdTomato + colonies among total CFU-Fs. n = 4 (Hes1-E12.5), n = 5 (Hes1-E12.5), n = 3 (Hes1-E14.5) mice. *p < 0.05, **p < 0.01. Two-tailed, Oneway ANOVA followed by Tukey's post-hoc test. Data are presented as mean ± SD. MB, Methylene blue staining. utilized as a tool to mark "skeletal stem/progenitor cells" at the condensation stage (13,26). However, the caveat is that Prrx1cre simultaneously marks SOX9 + cells within the condensation, making it impossible to discern the contribution of SOX9-negative peri-condensation mesenchymal cells. Our findings from Hes1-creER-based lineage-tracing experiments revise the concept regarding the origin of osteochondroprogenitor cells within the limb bud, which is currently considered to rest solely upon SOX9 + (Sox9-cre + ) cells (9). The emerging concept is that both SOX9 + cells within the condensation and SOX9-negative cells outside the condensation play equally important roles in early endochondral bone development. Identifying the unique role of each respective cell-of-origin of the limb bud mesenchyme remain as an important agenda for future studies.
In conclusion, we propose a new concept that pericondensation mesenchymal cells surrounding the SOX9 + condensation provide an important source of skeletal progenitor cells in early endochondral bone development (see the proposed diagram in Fig. 6).
Tamoxifen and induction of cre-loxP recombination Tamoxifen (Sigma T5648) was mixed with 100% ethanol until completely dissolved. Subsequently, a proper volume of sunflower seed oil (Sigma S5007) was added to the tamoxifenethanol mixture and rigorously mixed. The tamoxifenethanol-oil mixture was incubated at 60 C in a chemical hood until the ethanol evaporated completely. The tamoxifenoil mixture was stored at room temperature until use. Tamoxifen was injected at a dose of 3 mg into pregnant mice intraperitoneally using a 26 to 1/2-gauge needle (BD309597).

RNAscope in situ hybridization
Samples were fixed in 4% paraformaldehyde overnight at 4 C and then cryoprotected. Frozen sections at 14 μm were prepared on positively charged glass slides. In situ hybridization was performed with RNAscope 2.5 HD Reagent kit Brown (Advanced Cell Diagnostics 322,300) using the following probes: Hes1 (417701) and Sox9 (custom-designed) according to the manufacturer's protocol.

Imaging and cell quantification
Images were captured by an automated inverted fluorescence microscope with a structured illumination system (Zeiss Axio Observer Z1 with ApoTome.2 system) and Zen 2 (blue edition) software. C for 60 min on a shaking incubator. After cell dissociation, cells were mechanically triturated using an 18-gauge needle with a 1 ml Luer-Lok syringe (BD) and a pestle with a mortar (Coors Tek), and subsequently filtered through a 70 μm cell strainer (BD) into a 50 ml tube on ice to prepare single cell suspension. These steps were repeated for five times, and dissociated cells were collected in the same tube. Cells were pelleted and resuspended in an appropriate medium for subsequent purposes. For cell culture experiments, cells were resuspended in 10 ml culture medium and counted on a hemocytometer.
Single-cell RNA-seq analysis of fluorescence-activated cell sorting-isolated cells Cell sorting was performed using a four-laser BD fluorescence-activated cell sorting (FACS) Aria III (E x .407/ Figure 6. Hes1-marked early mesenchymal cells surrounding the condensation behave as skeletal progenitors in endochondral bone development. Early perichondrial population of skeletal progenitor cells in endochondral bone development. Undifferentiated mesenchymal cells surrounding condensation or outer perichondrial cells outside the cartilage template provide an important source during fetal endochondral bone development. During the condensation stage, Hes1 + undifferentiated mesenchymal cells take two distinct routes to participate in endochondral bone development. First, these cells translocate into the cartilage template and directly differentiate into Sox9 + chondrocytes and then contribute to bone formation (chondrocyte-dependent pathway). Second, these cells contribute to the outer layer of the perichondrium and become Hes1+ perichondrial cells. These Hes1+ perichondrial cells translocate into the nascent marrow space and directly differentiate into marrow mesenchymal cells by bypassing a Sox9+ state (chondrocyte-independent pathway). These perichondrial skeletal progenitor cells keep providing osteoblasts and bone marrow stromal cells even during the postnatal stage. As a result, postnatal bone marrow is characterized by transitional mosaicism composed of both chondrocyte-derived and perichondrium-derived stromal cells.
Hes1 + early mesenchymal cells in endochondral pathway 488/561/640 nm) high-speed cell sorter with a 100 μm nozzle. tdTomato + cells were directly sorted into ice-cold DPBS/1% BSA, pelleted by centrifugation and resuspended in appropriate amount of DPBS/1% BSA (1000 cells/μl). Cell numbers were quantified by Countless II automated Cell Counter (ThermoFisher) before loading onto the Chromium Single Cell 3 0 v2 microfluidics chip (10x Genomics Inc). cDNA libraries were sequenced by Illumina HiSeq 4000 using two lanes and 50 cycle paired-end read, generating a total of 770 million reads. The sequencing data was first pre-processed using the 10X Genomics software Cell Ranger. For alignment purposes, we generated and used a custom genome fasta and index file by including the sequences of tdTomato-WPRE to the mouse genome (mm10). Further downstream analysis steps were performed using the Seurat (20) R package. We filtered out cells with less than 500 genes per cell and with more than 20% mitochondrial read content. The downstream analysis steps include normalization, identification of highly variable genes across the single cells, scaling based on number of UMI, dimensionality reduction (PCA, CCA, and t-SNE), unsupervised clustering, and the discovery of differentially expressed cell-type specific markers. Differential gene expression to identify cell-type specific genes was performed using the nonparametric Wilcoxon rank sum test.

Colony-forming assay and subcloning
Nucleated bone marrow cells were plated into tissue culture 6-well plates (BD Falcon) at a density of <10 5 cells/cm 2 and cultured in low-glucose DMEM with GlutaMAX supplement (Gibco 10567022) and 10% mesenchymal stem cell-qualified FBS (Gibco 12662029) containing penicillin-streptomycin (Sigma P0781) for 1014 days. Cell cultures were maintained at 37 C in a 5% CO 2 incubator. Representative images of at least three independent biological samples are shown in the figures. For CFU-Fs, cells were fixed with 70% Ethanol for 5 min and stained for 2% methylene blue.

Statistical analysis
Results are presented as mean values ± SD. Statistical evaluation was conducted using the Mann-Whitney's U-test or one-way ANOVA. A p value of < 0.05 was considered significant. No statistical method was used to predetermine sample size. Sample size was determined on the basis of previous literature and our previous experience to give sufficient standard deviations of the mean so as not to miss a biologically important difference between groups. The experiments were not randomized. All of the available mice of the desired genotypes were used for experiments. The investigators were not blinded during experiments and outcome assessment. One femur from each mouse was arbitrarily chosen for histological analysis. Genotypes were not particularly highlighted during quantification.

Data availability
The data generated during and/or analyzed during the current study are available from the corresponding author upon reasonable request. The single-cell RNA-seq data presented herein has been deposited in the National Center for Biotechnology Information (NCBI)'s Gene Expression Omnibus (GEO) and are accessible through GEO Series accession number GSE144411 (https://www.ncbi.nlm.nih.gov/ geo/query/acc.cgi?acc=GSE144411).
Supporting information-This article contains supporting information.