Heterogeneity of murine periosteum progenitors involved in fracture healing

  1. Brya G Matthews  Is a corresponding author
  2. Sanja Novak
  3. Francesca V Sbrana
  4. Jessica L Funnell
  5. Ye Cao
  6. Emma J Buckels
  7. Danka Grcevic
  8. Ivo Kalajzic
  1. Department of Molecular Medicine and Pathology, University of Auckland, New Zealand
  2. Department of Reconstructive Sciences, UConn Health, United States
  3. Department of Physiology and Immunology, University of Zagreb, Croatia
  4. Croatian Intitute for Brain Research, University of Zagreb, Croatia

Abstract

The periosteum is the major source of cells involved in fracture healing. We sought to characterize progenitor cells and their contribution to bone fracture healing. The periosteum is highly enriched with progenitor cells, including Sca1+ cells, fibroblast colony-forming units, and label-retaining cells compared to the endosteum and bone marrow. Using lineage tracing, we demonstrate that alpha smooth muscle actin (αSMA) identifies long-term, slow-cycling, self-renewing osteochondroprogenitors in the adult periosteum that are functionally important for bone formation during fracture healing. In addition, Col2.3CreER-labeled osteoblast cells contribute around 10% of osteoblasts but no chondrocytes in fracture calluses. Most periosteal osteochondroprogenitors following fracture can be targeted by αSMACreER. Previously identified skeletal stem cell populations were common in periosteum but contained high proportions of mature osteoblasts. We have demonstrated that the periosteum is highly enriched with skeletal progenitor cells, and there is heterogeneity in the populations of cells that contribute to mature lineages during periosteal fracture healing.

Introduction

Bone tissue retains the ability to heal and regenerate throughout life. This process relies on tissue-resident stem and progenitor cells capable of generating new matrix. Until recently, most studies focused on the bone marrow compartment as a source of osteoprogenitors, with respect to their contribution to growth and remodeling, and their potential as a source of cells for regenerative applications. In the bone marrow compartment, it is well established that the endosteal region, as well as regions near trabecular bone surfaces, are highly enriched for stem and progenitor cells, including cells of the mesenchymal or skeletal lineage (Siclari et al., 2013; Morikawa et al., 2009). Skeletal stem cells (SSCs) can be isolated from various bone compartments based on their plastic adherence; however, more recently, numerous markers have been proposed to identify or isolate skeletal progenitors prospectively (Cao et al., 2020).

The skeleton has different requirements for the functionality of tissue-resident progenitors at different life stages. During development, the majority of the skeleton forms through endochondral ossification via a cartilage template, and this process continues during postnatal life at the growth plates. Hypertrophic chondrocytes can give rise to osteoblasts, particularly those in the primary spongiosa and trabecular region, likely via a progenitor intermediate (Wolff and Hartmann, 2019; Yang et al., 2014a; Yang et al., 2014b; Zhou et al., 2014a; Mizuhashi et al., 2018). However, this process is growth-related, and SSC pools must become established in adult tissues by the completion of growth when the growth plate activity reduces dramatically (as in mice) or fuses completely (in humans). Many stem and progenitor markers have been studied primarily during early life. For example, skeletal progenitors expressing Gli1, Gremlin1, and CTGF are prevalent during development and adolescent growth but appear to diminish or disappear in mice over 1 month of age (Shi et al., 2017; Worthley et al., 2015; Wang et al., 2015). Chan and colleagues have reported flow cytometry-based methods for identifying SSCs and downstream bone, cartilage, and stromal progenitors (BCSPs), as well as various other lineage-restricted populations in mice and humans (Chan et al., 2018; Chan et al., 2015). The murine cell surface antigen combinations were validated only in cells isolated from whole neonatal mouse long bones, but have since been applied to various settings in adult systems, including fracture healing (Marecic et al., 2015; Tevlin et al., 2017). Similarly, markers for putative SSCs in humans were identified using fetal tissue, with a focus on the hypertrophic cartilage zone, then applied to adult tissue. Another recent study used similar marker combinations to identify stem cell populations in the periosteum, with the majority of studies performed using very young animals (late embryonic up to P32, many at P6) (Debnath et al., 2018). PDGFRα+Sca1+ (PαS) cells were characterized in adult animals, and leptin receptor Cre-labeled cells that ultimately give rise to osteoblasts and bone marrow adipocytes do not become established before 2 months of age (Morikawa et al., 2009; Zhou et al., 2014b). These results indicate that stem and progenitor cell populations and markers change at different life stages, meaning that markers identified in neonates or juveniles may not apply to the adult setting.

The periosteum is the tissue surrounding the outer surface of the bone. It is physically separated from the bone marrow compartment and appears to have a much greater regenerative capacity, playing a critical role in fracture healing (Colnot, 2009; Roberts et al., 2015). Periosteum-based healing involves a process similar to endochondral ossification. Following an inflammatory phase, the periosteum expands and forms a callus initially comprised of fibrocartilage in the central area, combined with direct bone formation at the periphery. It has been challenging to identify markers that are specific to the periosteum; however, recent studies have characterized cathepsin K and periostin as markers that are highly enriched in periosteal mesenchymal cells, and alpha smooth muscle actin (αSMA) was shown to enrich Mx1-labeled periosteal cells for progenitors (Debnath et al., 2018; Duchamp de Lageneste et al., 2018; Ortinau et al., 2019). We have previously shown that αSMA identifies a population of growth-related osteoprogenitors in the trabecular region, and cells that contribute to fracture healing in the adult periosteum (Matthews et al., 2014; Grcevic et al., 2012).

In the current study, we have used murine models to characterize the periosteal progenitor populations. We have demonstrated that αSMACreER identifies a subset of long-term osteochondroprogenitors involved in fracture callus formation. In addition, we have evaluated different bone compartments for the expression of stem/progenitor markers, including cell surface profile, SSC, and BCSP populations, and lifespan by label retention assays in different bone compartments, and evaluated growth potential of different periosteal populations ex vivo.

Results

In vivo identification of long-term osteochondroprogenitors

In order to identify murine periosteal progenitor populations in vivo, we initially utilized a lineage tracing approach with a tamoxifen-inducible Cre combined with the Ai9 Tomato (Tom) reporter. We have demonstrated that αSMACreER identifies osteoprogenitors in both the periosteum and the bone marrow compartment (Matthews et al., 2014; Grcevic et al., 2012; Matthews et al., 2020). Previously, we delivered tamoxifen at the time of fracture, which potentially labels cells that activate αSMA expression at the time of or shortly after fracture (Mori et al., 2016). In this study, we delivered tamoxifen up to 90 days before fracture to evaluate whether the reporter was identifying long-term progenitor cells. First, we confirmed that αSMA-labeled cells were retained in the periosteum over the experimental time course by both histology (Figure 1A–B) and flow cytometry (Figure 1C). αSMA-labeled cells remained in the periosteum at a similar frequency over the timeframe evaluated and also remained in the endosteum (Figure 1D). In juvenile mice (4–5 weeks old), we see contribution of αSMA+ cells from the bone marrow compartment, particularly in the trabecular region, to osteoblasts and osteocytes after 17 days (Figure 1—figure supplement 1A–BGrcevic et al., 2012). However, 70 days after tamoxifen delivery, osteoblasts are no longer labeled, but some Tomosteocytes remain (Figure 1—figure supplement 1C–F). In contrast, when tamoxifen is delivered to sexually mature mice, there is no contribution to osteoblasts in trabecular bone up to 100 days after tamoxifen (Figure 1—figure supplement 2).

Figure 1 with 3 supplements see all
Alpha smooth muscle actin (αSMA) identifies long-term progenitor cells in the periosteum.

Representative images of tibia mid-shaft periosteum: (A) 14 days or (B) 56 days after tamoxifen with the periosteal surface near the top and endosteal surface near the bottom of the image. Percentage of αSMA-labeled (Tom+) cells within the CD45- population at different time points in (C) periosteum and (D) endosteum as measured by flow cytometry (n = 2–5). Fracture histology of representative day 7 tibial fracture calluses in different groups. Diagrams indicate the experimental design. (E) αSMA-labeled with tamoxifen (Tam) at day −1, 0; (F) αSMA-labeled with tamoxifen at day −91, –90. Quantification of Col2.3GFP+ osteoblasts derived from αSMA-labeled cells in the fracture callus based on image analysis (Tom+GFP+ cells/GFP+ cells) in (G) day 7 fracture callus and (H) day 14 fracture callus at different times after tamoxifen administration. (I) Quantification of Sox9+ chondrocytes derived from αSMA-labeled cells in day 7 fracture callus. All images show mineralized bone by dark field (gray) and DAPI+ nuclei in blue. Scale bars are 100 µm (A, B) or 250 µm (E–F). In F, the boxed area is magnified and αSMA-labeled osteoblasts as quantified in G-H are indicated with arrowheads. Approximate callus area used for analysis is indicated by the dashed lines in E-F, cortical bone surface defines the inner boundary. *p<0.05 compared to no tamoxifen control with ANOVA followed by Dunnett’s posttest. #p<0.05 compared to tamoxifen day 0 group.

Figure 1—source data 1

Image analysis processed data and statistical analysis tables (from Graphpad Prism) for graphs in Figure 1G–I and supplement 3B-C.

https://cdn.elifesciences.org/articles/58534/elife-58534-fig1-data1-v2.xlsx

Following tibial fracture, αSMA-labeled cells made substantial contributions to fracture callus osteoblasts and chondrocytes (Figure 1E–I). Very few Tom+ cells were present or contributed to mature cell types in the absence of tamoxifen (Figure 1G–I, Figure 1—figure supplement 3A). Most fracture callus osteoblasts are derived from cells that expressed αSMA at some point in their lineage, with 81.5% derived from αSMA+ cells labeled just after fracture. A smaller population of osteoblasts were derived from long-term αSMA-labeled progenitors (around 20%, Figure 1G–H). The contribution of αSMA-labeled cells was similar in both the initial intramembranous bone formation in the day 7 callus and the more fully mineralized day 14 callus. Contribution to chondrocytes showed a very similar trend based on the timing of tamoxifen; however, the proportion of Tom+ chondrocytes was always lower (Figure 1I). Following the pin insertion and fracture, bone formation also occurs inside the marrow compartment, particularly in the proximal diaphysis surrounding the pin over 1 mm from the fracture site (Figure 1—figure supplement 3). αSMA-labeled cells labeled at the time of fracture make a major contribution to newly formed osteoblasts in the marrow compartment; however, there was lower contribution at the later time points suggesting that αSMA identifies long-term injury-responsive osteochondroprogenitors in the periosteum but does not consistently label long-term injury-responsive cells in the bone marrow compartment.

Putative mesenchymal stem cell markers are enriched in periosteum

We performed further characterization of cell populations in the periosteum, first, to evaluate surface profile using putative mesenchymal stem cell or SSC markers and their corresponding growth potential and, second, to characterize αSMA-labeled cell marker expression. There are numerous markers for SSCs and progenitor cells reported in the literature; however, many have not been evaluated systematically in different tissue compartments. We evaluated the expression of several cell surface markers in periosteum, endosteum, and bone marrow of adult mice (Figure 2A, Figure 2—figure supplement 1). All stains included a CD45/CD31/Ter119 cocktail to exclude hematopoietic and endothelial cells and select for mesenchymal lineages (termed CD45- population). We confirmed that the freshly isolated CD45+ population was unable to form fibroblast colony-forming units (CFU-F) or any type of colonies in the CFU-F assay. The CD45- population in periosteum was highly enriched with most of the markers evaluated compared to both endosteum and bone marrow (Figure 2A). Sca1, PDGFRα, and PDGFRβ were abundant in periosteum, but expressed on <2% CD45- endosteal and bone marrow cells. The only marker that did not follow this trend was CD105, which was consistently present on the majority of endosteal and bone marrow CD45- cells. Various marker combinations that have previously been used to identify stem cells such as Sca1/PDGFRα or Sca1/CD51 were easily identifiable in periosteum but very rare in the endosteum or bone marrow (Figure 2—figure supplement 1). The CD45- periosteal population formed around 20× more CFU-F than the CD45- endosteal population, suggesting functional enrichment of stem and progenitor cells (Figure 2B). This is in agreement with Duchamp de Lageneste et al. (Duchamp de Lageneste et al., 2018) who used a different method to isolate periosteal cells. We went on to evaluate the expansion and differentiation potential of selected periosteal populations independently of αSMA expression using in vitro assays. CFU-F formation was restricted to CD51+ cells, and most enriched in Sca1+ CD51+ cells (Figure 2—figure supplement 2A). CD90+ cells also showed high CFU-F potential, whereas CD105+ cells, which are almost exclusively CD90-, formed few or no colonies. Sca1+ CD51+ cells showed the ability to differentiate toward osteogenic (ALP+) and adipogenic (Oil Red O+) lineages, in many cases within a single colony (Figure 2—figure supplement 2B-C). Similar results were obtained for the CD90+ population. Sca1- CD51+ cells and CD90- CD105- cells efficiently activated ALP but never differentiated into adipocytes. This suggests that in addition to including some mature osteoblasts, these populations contain committed osteoprogenitors that remain capable of density-independent growth. All these populations contained a small proportion of αSMA-labeled cells, with the highest frequency of αSMA-labeled cells in the CD90+ population (Figure 2—figure supplement 2D).

Figure 2 with 2 supplements see all
Mesenchymal stem cell markers are enriched in periosteum.

Flow cytometry and cell sorting was used to evaluate marker expression in tissue from adult long bones. (A) Expression of cell surface markers in periosteum, endosteum, and bone marrow samples as a percentage of the CD45/Ter119/CD31- fraction, n = 4–50. Expression in periosteum is significantly different from both bone marrow and endosteum for all markers (two-way ANOVA with posttest comparing cell types). (B) Colony-forming unit fibroblast (CFU-F) formation in the total CD45/Ter119/CD31- fraction. (C) Representative dot plot showing gating of alpha smooth muscle actin (αSMA)-labeled cells. (D) Proportion of αSMA-labeled cells within the CD45/Ter119/CD31- gate 2–3 days after tamoxifen in different tissue compartments (n = 17). Cell surface marker expression in the Tom+ population 2–3 days after tamoxifen in (E) periosteum and (F) endosteum (n = 4–17). (G) Cell surface profile of αSMA-labeled periosteum cells at different time points (n = 2–7). *p<0.05 compared to CD45- population.

Figure 2—source data 1

Data tables and statistical analysis tables (from Graphpad Prism) for graphs shown in Figure 2 and figure supplement 2.

https://cdn.elifesciences.org/articles/58534/elife-58534-fig2-data1-v2.xlsx

We next performed further characterization of αSMA-labeled cells in the absence of injury. By flow cytometry, we confirmed that αSMA+ cells were much more prevalent in periosteum compared to endosteum and bone marrow in adult mice (Figure 2C–D), although they were rarer than the other markers we analyzed. We evaluated expression of various markers in the αSMA+ population. The majority of αSMA-labeled cells in the periosteum and endosteum express CD51 and CD90 (Figure 2E–F). The other markers tested were expressed in smaller subsets of αSMA-labeled cells, indicating there is heterogeneity within the labeled population. We next evaluated whether the prevalence of any of these markers in the αSMA-labeled population changed over time since labeling. We did not detect differences in the level of expression of any of the markers evaluated individually or in various combinations at different times after tamoxifen, suggesting that additional markers may be required to define the heterogeneity and phenotype of these cells (Figure 2G).

To further evaluate heterogeneity, we performed plate-based single cell RNAseq analysis on αSMA-labeled periosteal cells isolated at different time points after tamoxifen administration (Figure 3). We obtained data from 185 cells from 2 days after tamoxifen, 197 cells from 6 weeks, and 147 cells from 13 weeks after tamoxifen. The median unique molecular identifier (UMI) number was 697, and the median number of expressed genes was 449 (Figure 3—figure supplement 1). Since expression levels for the majority of genes were low in this study, we restricted our analysis to basic clustering studies. Clustering using the Leiden community detection algorithm generated four clusters from the αSMA-labeled cells, confirming heterogeneity in this population (Figure 3A, Figure 3—figure supplement 2Traag et al., 2019). Clusters 1–3 are similar and appear to represent mainly periosteal populations. Cluster 4 is quite separate from the rest and appears to be derived from muscle contamination, expressing markers such as skeletal muscle actin, Acta1, and myoglobin, Mb (Figure 3—figure supplement 2B). We have previously demonstrated that αSMA labels a subset of muscle satellite cells that go on to generate myofibers in both juvenile and adult animals (Matthews et al., 2016). Notably, of the 26 cells in this cluster, only one was from the day 2 chase group, consistent with the differentiation of αSMA+ satellite cells into mature muscle lineage cells over time.

Figure 3 with 2 supplements see all
Single cell RNAseq analysis of αSMA-Tom+ cells.

Plate-based RNAseq was performed on alpha smooth muscle actin (αSMA)-Tom+ cells at up to 13 weeks after tamoxifen delivery. (A) UMAP cluster plot indicating four clusters. (B) Violin plots of top markers for clusters 1–3. (C) Distribution of cells in clusters 1–3 at different time points following tamoxifen. (D) Heat map of expression of putative skeletal stem and periosteal progenitor markers.

Figure 3—source data 1

scRNAseq cluster makers with AUROC values >0.7.

https://cdn.elifesciences.org/articles/58534/elife-58534-fig3-data1-v2.xlsx

Cluster 3 is characterized by expression of a number of pericyte markers, including Rgs5, and frequent expression of PDGFRβ and αSMA. They also frequently express Cxcl12, which is a marker of bone marrow resident reticular cells, and Fabp4, which is known as an adipocyte lineage marker (Figure 3B, Figure 3—figure supplement 2A). Other adipocyte markers like Adipoq and Pparg were undetectable in most cells. The proportion of cluster 3 cells increased slightly over time, suggesting they might represent a population of mature perivascular cells (Figure 3C). Clusters 1 and 2 show some overlap, and cluster 1 had few selective markers, with only the matrix protein gelsolin, Gsn, passing our marker criteria. Interestingly, cluster 1 shows the most frequent expression of the genes for Sca1, PDGFRα, and CD90, suggesting it may contain multipotent stem/progenitor populations (Figure 3D). The proportion of cells in cluster 1 was reduced following longer chase periods, in line with the reduction in αSMA-lineage contribution to fracture healing at later time points. Markers of cluster 2 are primarily matrix genes that can be expressed by the osteoblast lineage, including Col1a2, Tnmd, Fmod, and Ogn. Unpublished studies we performed in parallel with differentiated osteoblasts show consistently high levels of Bglap in differentiated osteoblasts, but Bglap is absent or present at only low levels in the αSMA-labeled cells, indicating that cluster 2 does not represent mature osteoblast lineage cells. It is feasible that at least some of these cells are from the fibrotendonous lineage, as the periosteum contains fascia and fibrous tendon attachments, which appear to be targeted by αSMA histologically. Overall, our data indicate that αSMA+ cells are a fairly rare, but heterogeneous cell population within the periosteum that can be separated into at least three different cell types.

Label-retaining cells in the periosteum are enriched with progenitor markers

One feature of long-term stem cells is a slow cycling rate, and this can be evaluated experimentally using label retention assays. We used a doxycycline (dox)-inducible GFP-tagged histone system, H2B-GFP, which is incorporated into the DNA during dox feeding, then diluted during subsequent cell division. This system has been used for enrichment of highly self-renewing populations of hematopoietic stem cells and muscle satellite cells (Foudi et al., 2009; Chakkalakal et al., 2014). The experimental design for generating label-retaining cells (LRCs) is shown in Figure 4A. H2B-GFP was incorporated into the majority of periosteal cells (Figure 4B, Figure 4—figure supplement 1A). Quiescent terminally differentiated cells such as osteocytes also effectively retain the label (Figure 4—figure supplement 1B). Flow cytometry was used to track GFP expression in the CD45- fraction of the periosteum (Figure 4B–C). Our periosteal isolation procedure excludes bone tissue and therefore osteocytes, and we have previously demonstrated that our endosteal preparations contain almost no osteocytes (Matic et al., 2016). After the chase phase, the periosteum contained a small population of true LRCs (GFPhi), as well as a slightly larger GFPint population which would be expected to contain some LRCs. Histologically, the LRCs were mostly in the inner cambium layer of the periosteum (Figure 4—figure supplement 1B). We confirmed that most cells capable of forming CFU-F initially incorporated a GFP label (Figure 4D). After the chase, both GFP+ populations showed enrichment of CFU-F compared to the GFP- fraction (Figure 4E). We also evaluated cell surface marker expression in LRCs and found strong enrichment of Sca1 and CD51, more frequent expression of CD90, and fewer CD105+ cells in periosteum (Figure 4F) and endosteum (Figure 4G) compared to GFP- cells. Since label retention alone was not very effective at enriching CFU-F, we combined this approach with αSMA-labeling. Initially, we labeled 73.4 ± 7.6% of αSMA+ cells, and after at least 13 weeks chase, 22.6 ± 7.7% were LRCs, indicating a large proportion of these cells are quiescent (Figure 4—figure supplement 2A). This is a substantially higher proportion of LRCs than we found in the overall CD45- population of the same samples (3.2 ± 1.5%, p<0.0001). Ex vivo, αSMA-labeled LRCs (GFP+Tom+) were enriched for CFU-F compared to the GFP-Tom+ population (Figure 4—figure supplement 2B). Interestingly, αSMA-labeled LRCs appeared to be lineage-restricted, with evidence of differentiation along the osteogenic lineage but absence of adipogenic differentiation (Figure 4—figure supplement 2C). αSMA-labeled LRCs were more frequently CD51+, but fewer were CD90+ that Tom+ non-LRCs (Figure 4—figure supplement 2D).

Figure 4 with 2 supplements see all
Label-retaining cells are enriched with most progenitor markers.

(A) Experimental design for label-retaining cell experiments. (B) Representative gating of the CD45- population in different periosteum samples to define three populations based on GFP expression. (C) Quantification of different GFP populations in a representative experiment for periosteum and endosteum (n = 1, no doxycycline [dox]; n = 4, no chase; n = 6, 13-week chase). (D) Colony-forming unit fibroblast (CFU-F) frequency in different CD45- GFP populations in periosteal no chase samples (n = 3) and (E) following chase (n = 6). Cell surface marker expression is shown in different CD45- GFP subsets following chase (populations equivalent to right-hand label-retaining cell plot in B) in (F) periosteum cells and (G) endosteal cells. Data is pooled from three separate experiments involving 13–22 weeks chase (n = 19). *p<0.05, repeated measures ANOVA with Dunnett’s posttest compared to GFP- population.

Figure 4—source data 1

Data tables and statistical analysis tables (from Graphpad Prism) for graphs shown in Figure 4F–G.

https://cdn.elifesciences.org/articles/58534/elife-58534-fig4-data1-v2.xlsx

Evaluation of SSC populations in adult mice

In order to systematically evaluate stem cell populations in adult animals, we evaluated the utility of the stain described by Chan et al. to identify SSC and BCSP (Chan et al., 2015). Stem/progenitor markers and gating strategies should, by definition, exclude mature cell types of the lineage. There are no validated methodologies available to identify osteoblasts using cell surface markers, although some authors report Sca1+CD51+ as stem cells and Sca1-CD51+ cells as committed osteoprogenitors or osteoblasts (Arthur et al., 2019; Lundberg et al., 2007). We utilized the Col2.3GFP transgene to identify mature osteoblasts (Kalajzic et al., 2002; Dacquin et al., 2002). Col2.3GFP+ cells were always present in endosteal samples (5.4 ± 5.0% of CD45-) and were detectable in numbers sufficient for analysis in some periosteum samples (1.3 ± 1.3% of CD45-). CD200 expression was present on over 90% of Col2.3GFP osteoblasts in both endosteum and periosteum, suggesting that CD200 may be useful for positive selection of osteoblasts. Around 35% of endosteal osteoblasts were CD51+ (Figure 5AMatic et al., 2016). A subset of osteoblasts in both tissues were CD105+, while the other markers including Sca1, CD90, and Ly51 (also known as 6C3) were rare or absent in endosteal osteoblasts. The gating strategy to identify SSCs and BCSPs in both periosteum and endosteum is shown in Figure 5C–D. We found that Col2.3GFP+ osteoblasts constituted almost 40% of endosteal cells defined as SSC (CD45-CD51+CD90-Ly51-CD105-CD200+) and 20% defined as BCSP (CD45-CD51+CD90-Ly51-CD105+) (Figure 5D). The proportion of mature osteoblasts in these populations was lower in periosteal samples but still represented 7–8% of SSC and BCSP (Figure 5C). The pre-BCSP population (CD45-CD51+CD90-Ly51-CD105-CD200-) contains very few osteoblasts. In addition, when we evaluate the proportion of SSC and BCSP in different tissue compartments, we find that the GFP+ osteoblast populations, particularly in the endosteum, have the highest frequency of these ‘stem cells’ (Figure 5B). Therefore, these marker combinations used for SSC and BCSP do not exclude mature cells within the osteoblast lineage in adult mice, so do not specifically identify stem cells in this setting.

Bone, cartilage, and stromal progenitor (BCSP) and skeletal stem cell (SSC) populations in adult bone contain mature osteoblasts.

(A) Marker expression on Col2.3GFP+ osteoblasts, n = 4–32. (B) Frequency of BCSP and SSC populations in different tissue compartments and populations (n = 8–12). (C–D) Gating to identify BCSPs, pre-BCSPs, SSCs, and presence of Col2.3GFP+ osteoblasts as a major component of these populations in periosteum (C) and endosteum (D). Mean ± standard deviation for % GFP+ cells are shown.

αSMA identifies transplantable, self-renewing osteoprogenitors

In ex vivo assays, αSMA+ periosteal cells were capable of CFU-F formation (Figure 6A). Following transplantation into a critical-sized calvarial defect, αSMA+ cells were readily identified in all sections evaluated, indicating the ability to engraft and expand in the defect (Figure 6B). Some cells activated Col2.3GFP, indicating limited differentiation into osteoblasts, similar to a recent study (Ortinau et al., 2019). We did not see evidence of cartilage formation in these defects at the time point evaluated. We tested the in vivo self-renewal potential of αSMA-labeled cells by performing a fracture, allowing it to heal for 8 weeks, then refracturing the same bone in a similar location. When tamoxifen was given at the time of initial fracture, around 20% of both osteoblasts and chondrocytes in the secondary fracture callus are derived from the original αSMA-labeled population (Figure 6C–J). This is 35–50% of the contribution in the first fracture but indicates self-renewal of a proportion of αSMA-labeled progenitors. Contribution to the secondary fracture when tamoxifen was given 14 days prior to the first fracture was not statistically significantly different to the day 0 group, confirming identification of self-renewing osteochondroprogenitors.

Alpha smooth muscle actin (αSMA)-labeled cells are transplantable self-renewing osteochondroprogenitors.

(A) Colony-forming unit fibroblast (CFU-F) frequency of periosteal αSMA-labeled cells sorted 2 days after tamoxifen. (B) Engraftment of transplanted αSMA-labeled cells in a calvarial defect 5 weeks after transplantation. Arrowheads indicate cells that are Col2.3GFP+. Quantification of αSMA-labeled cell contribution to (C) osteoblasts and (D) chondrocytes in day 7 secondary fracture calluses is shown. (E) Experimental design and fracture callus histology following a second tibia fracture in αSMACreER/Ai9/Col2.3GFP mice showing contribution of Tom+ cells to callus tissue, including osteoblasts (F) and Sox9+ chondrocytes (white) (G). (H–J) Histology of a similar secondary fracture where tamoxifen was given 14 days prior to the primary fracture. White arrowheads indicate αSMA-labeled Col2.3GFP+ osteoblasts, orange arrowheads indicate αSMA-labeled Sox9+ chondrocytes. DAPI stain is shown in blue. Approximate callus area used for analysis is indicated by the dashed lines in E and H. Scale bars indicate 250 µm (low magnification) and 100 µm (high magnification).

Figure 6—source data 1

Data tables and statistical analysis tables (from Graphpad Prism) for graphs shown in Figure 6C-D.

https://cdn.elifesciences.org/articles/58534/elife-58534-fig6-data1-v2.xlsx

Functional contribution of αSMA+ cells to fracture healing

We evaluated the functional importance of αSMA+ cells during fracture healing using the Rosa-DTA model of cellular ablation. Diphtheria toxin (DTA) expression, which causes death of the cells where it is expressed, was driven by αSMACreER and induced by tamoxifen delivery around the time of fracture. When ablation was performed at the time of fracture (days 0, 1), we saw a 50% reduction in αSMA-labeled cells in 4 days post-fracture (DPF) callus tissue by flow cytometry (Figure 7A). This is consistent with the 50% efficiency of recombination we demonstrated for αSMACreER in vitro (Sinder et al., 2020). Histologically, there is a clear reduction in fracture size and αSMA-labeled cells at 4 DPF (Figure 7B). Histological analysis at 14 DPF revealed a significant reduction in total callus size as well as mineralized area (Figure 7C–D). Micro-computed tomography (microCT) results at day 14 confirmed the histological observations (Figure 7E). At 21 DPF, microCT showed no change in total callus volume but a significant decrease in callus bone mass (Figure 7F). Similar impairments in healing were seen whether ablation was performed at the time of fracture or throughout the first 8 days of healing, suggesting that early loss of αSMA+ cells was most important for this phenotype. Overall, our data demonstrate that partial ablation of αSMA+ progenitors reduced callus formation and delayed mineralization. The less severe phenotype at 21 DPF suggests this may be a delay in healing rather than complete disruption of the process, although the fractures from mice with αSMA-ablation probably never reach the size of control calluses. Healing in the DTA mice may still progress due to compensation of αSMA+ progenitors that escaped ablation or larger contribution from αSMA-negative progenitors in this setting.

Ablation of alpha smooth muscle actin (αSMA)+ cells impairs fracture healing.

(A) Flow cytometry of day 4 fracture callus from control αSMACreER/Ai9 mice and DTA animals (αSMACreER/Rosa-DTA/Ai9) indicating significantly fewer (p<0.01) αSMA-labeled cells in the ablation model (mean ± SD shown, n = 4). (B) Representative histological sections from day 4 fracture callus in controls or with DTA-mediated ablation. The outline of the callus is indicated. (C) Histological analysis of callus area and mineralized area at day 14, and (D) representative images of von-Kossa staining (n = 4–7). (E) Micro-computed tomography (MicroCT) analysis of day 14 fracture callus in DTA- and DTA+ mice that had undergone two different tamoxifen regimens to ablate αSMA+ cells, day 0, 1 of fracture or day 0, 2, 4, 6, 8. Callus tissue volume and bone mass were measured. Cre-: n = 14, 8 male (M), 6 female (F); DTA D0-1: n = 10, 4M 6F; DTA D0-8: n = 8, 4M 4F. (F) MicroCT analysis of day 21 fracture callus in mice treated with tamoxifen on day 0, 2, and 4 (n = 14–15).

Figure 7—source data 1

Data tables and statistical analysis tables (from Graphpad Prism) for graphs shown in Figure 7C,E and F.

https://cdn.elifesciences.org/articles/58534/elife-58534-fig7-data1-v2.xlsx

Contribution of other periosteal populations to fracture healing

We have demonstrated that αSMA identifies some, but not all, of the cells contributing to osteoblasts and chondrocytes in the fracture callus. We therefore examined the potential contribution of other populations. We tested the previously unexplored hypothesis that mature osteoblast lineage cells may contribute to osteoblast formation in the fracture callus. Col2.3CreER-labeled cells made a small but significant contribution to osteoblast formation in the fracture callus, but only when tamoxifen was given close to the time of fracture (Figure 8.A-B). This is in line with previous studies that estimate a lifespan of around 2 weeks for osteoblasts (Matic et al., 2016; Jilka et al., 1998). We did not detect any contribution of Col2.3-labeled cells to chondrocytes. This suggests that committed osteoblast lineage cells present in the periosteum can contribute to fracture osteoblasts but do not have the plasticity to contribute to chondrocytes.

Contribution of mature osteoblasts to fracture healing and PDGFRα distribution.

(A) Fracture histology of a representative day 7 tibial fracture callus from a Col2.3CreER/Ai9 mouse treated with tamoxifen at days −4 and −3. Mineralized bone is shown by dark field (gray) and DAPI+ nuclei in blue. Scale bar is 250 µm. (B) Quantification of Col2.3-labeled cell contribution to Osx+ osteoblasts. *p<0.05 based on ANOVA with Dunnett’s posttest compared to no tamoxifen group. (C) Tibia histology of PdgfraCreER/Ai9 mice given three daily doses of tamoxifen and euthanized 1 day later. Scale bar is 200 µm. Magnified areas of (D) trabecular and (E) cortical bone are shown. (F) Tibial histology of a PdgfraCreER/Ai9 mouse that had not been exposed to tamoxifen showing minimal leakiness of the CreER. Representative images of n = 3 aged 5–8 weeks of age are shown. (G) PDGFRα immunostaining of bone. Scale bar is 50 µm. bm, bone marrow; gp, growth plate; p, periosteum.

Figure 8—source data 1

Data table and statistical analysis table (from Graphpad Prism) for graph shown in Figure 8B.

https://cdn.elifesciences.org/articles/58534/elife-58534-fig8-data1-v2.xlsx

Finally, we further evaluated PDGFRα as an alternate marker of periosteal progenitor cells as we hypothesized that it may target a broader periosteal progenitor population than αSMA. PDGFRaCreER targeted the majority of periosteal cells 1 day after completion of tamoxifen administration, as well as the majority of osteoblasts on trabecular and endocortical bone surfaces (Figure 8C–E). It also targeted a population of cells in the bone marrow. This labeling was tamoxifen-dependent (Figure 8F). These results were confirmed with immunostaining for PDGFRα that indicated expression in the periosteum and in endosteal osteoblasts (Figure 8G). There is some discrepancy between our flow data and this result (see Figure 2A and Figure 5A). PDGFRα staining by flow cytometry tends to be weak, and there is rarely a clear demarcation between positive and negative cells (Figure 2—figure supplement 1). In contrast, PDGFRα-driven Cre in combination with a strong reporter gene such as Ai9 amplifies even low transgene activity. It is therefore likely that we are only gating for cells highly expressing PDGFRα in flow cytometry, but this makes it difficult to use as a reliable marker. This discrepancy is consistent with the results of Ambrosi et al. who report much wider expression of PDGFRα using a GFP reporter mouse compared to previous results with antibody staining and flow cytometry (Morikawa et al., 2009; Ambrosi et al., 2017). Their data also suggests expression of this reporter in osteoblasts. Another recent study reported much higher PDGFRα expression in bone marrow compared to periosteum in contrast with our flow cytometry results (Tournaire et al., 2020). We conclude that PDGFRα is not an optimal marker of progenitor cells within the periosteum due to both technical variability depending on the method used and expression in mature osteoblasts.

Discussion

Our results indicate that the periosteum of adult mice is highly enriched for cells with markers and phenotype of stem and progenitor cells. Together with other recent studies, our results highlight the presence of unique cell populations within the periosteum that are rare or absent in the bone marrow compartment, and demonstrate that stem cell markers can change with age (Debnath et al., 2018; Duchamp de Lageneste et al., 2018; Ortinau et al., 2019; Tournaire et al., 2020). There is abundant evidence suggesting that fracture healing is impaired if periosteum is removed, or that otherwise critical-sized defects in a variety of bones can regrow if periosteum is retained (Mashimo et al., 2013; Kuwahara et al., 2019; Ozaki et al., 2000; Zhang et al., 2005). There are also numerous studies reporting superior outcomes with cultured periosteal cells compared to other cell sources like bone marrow in terms of bone formation both in vitro and in vivo (Duchamp de Lageneste et al., 2018; van Gastel et al., 2014). Of the markers we examined, we found that most were much more abundant in the periosteum than endosteum. Sca1, CD51, CD90, αSMA, and PDGFRα (Wang et al., 2019a) all enrich for CFU-F in periosteum, similar to what has been reported in other cell types. Although we did not thoroughly characterize PαS cells in the periosteum, we note that these cells were very rare in the endosteal compartment. Morikawa et al. report large enrichment of PαS cells following collagenase digestion of crushed bones (Morikawa et al., 2009). Our results raise the possibility that some of the cells isolated in their studies were actually derived from the periosteum. We confirmed by immunostaining and use of an inducible Cre that PDGFRα expression was very abundant in the periosteum. This result and abundant PDGFRα expression shown histologically in osteoblasts were not completely replicated in our flow cytometry data. Ambrosi et al. reported that all Sca1+ cells in the bone marrow co-express PDGFRα using a PDGFRα-GFP line which contrasted with the original study characterizing PαS cells (Ambrosi et al., 2017). This highlights the difficulty in using PDGFRα consistently as a stem and progenitor cell marker. With respect to CD51, our results agree with the studies using neonatal bones, indicating that progenitor cells are solely contained within the CD51+ population (Chan et al., 2009); however, a significant proportion of mature osteoblasts are CD51+, so additional markers that exclude differentiated cells are required. Finally, our results suggest that strong Sca1 expression is specific to periosteal mesenchymal cells and enriches multipotent stem or progenitor cells in combination with CD51, at least in vitro. Sca1 expression was also highly enriched in LRCs in all bone tissue compartments. A recent study using scRNAseq in periosteal cells suggested Sca1 as a marker of differentiated cells in the periosteum, although it was unclear what cell type these mature cells were (Debnath et al., 2018). Inducible Cre models driven by expression of the Sca1 gene are now available and would potentially help to confirm the in vivo function of Sca1+ cells. Notably, they target endothelial progenitors which will complicate any analysis with these models (Vagnozzi et al., 2018; Tang et al., 2018), and one of these models was reported to target very few cells within the bone (Vesprey et al., 2020).

LRCs in the periosteum showed large differences in cell surface marker profile compared to non-LRCs and enrichment of CFU-F formation. However, the experimental system presented some limitations, as we did not achieve uniform labeling initially, and some cells had low GFP expression even in the absence of dox. This ‘leaky’ GFP expression may be related to positioning of the tetO-H2B-GFP cassette in the Col1a1 locus. Following chase, we found that CFU-F formation from the CD45- GFPint population was similar, or in some cases higher than the true GFPhi LRCs. The GFPint population potentially contains a mix of cells including those with dox-independent GFP expression, those who have lost some but not all label, and LRCs that did not achieve bright labeling initially. We are unable to distinguish between these populations. It is also possible that the true LRC population contains terminally differentiated quiescent cells that are present in the periosteum. Finally, CFU-F formation appears to be a feature of progenitor cells as well as stem cells, and progenitors may be enriched in the GFPint population. The use of label retention is likely to perform better in combination with other markers that already enrich stem cell populations (Foudi et al., 2009; Chakkalakal et al., 2014). Within the αSMA-labeled population, we saw fourfold higher CFU-F frequency than αSMA-labeled non-LRCs. A previous study also reported the presence of heterogeneous LRCs in the periosteum; however, they were unable to perform functional characterization, while Nes-GFP+ periosteal cells were not LRCs (Tournaire et al., 2020; Cherry et al., 2014). Interestingly, the previous study indicated that αSMA+ cells were not LRCs; however, many antibodies to αSMA stain perivascular cells in larger blood vessels but do not detect αSMA+ cells in the cambium layer of the periosteum (Cherry et al., 2014). The αSMA+ cells in the cambium layer appear to be the cell subset involved in fracture healing (Ortinau et al., 2019; Wang et al., 2019b). Label retention may be a useful tool to enrich periosteal stem cells in ex vivo assays, although more robust assays for testing self-renewal such as serial transplantation are required.

One of the key findings of this study is that αSMA identifies long-term, slow-cycling, self-renewing osteoprogenitors in the adult periosteum that are functionally important for bone and cartilage formation during fracture healing. αSMA+ periosteal cells also contribute to osteoblasts in response to loading (Matthews et al., 2020). The lineage tracing results suggest that αSMA identifies at least two distinct osteochondroprogenitor populations in the periosteum that contribute to fracture healing. First, there is a long-term quiescent tissue-resident progenitor population that contributes to a subset of osteoblasts and chondrocytes for at least 3 months after labeling. This population is presumably re-established following fracture and contributes to healing of a second fracture in the same bone. Second, there are cells that activate αSMA after injury, as visualized in studies using real-time αSMA reporters (Mori et al., 2016). The majority of osteoblasts are derived from αSMA-labeled cells when tamoxifen is given at the time of fracture. The results from the fractures performed after a 2-week chase, which show slightly lower contribution than at the time of fracture but a trend toward higher contribution than after 6 or more weeks chase, suggest that there may be a third population of tissue-resident osteochondroprogenitor cells that have a limited lifespan or lack self-renewal capacity. RNAseq analysis of αSMA periosteal cells in intact bone suggested there are at least three cell populations present. While we cannot be certain which of these contains stem and progenitor cells, marker gene expression suggests that one cluster is perivascular cells, another is more fibroblast-like, and the intermediate population is likely to contain stem and progenitor cells. Unfortunately, few of the cluster marker genes were cell surface markers suitable for flow cytometry-based separation of populations.

We consistently saw lower contribution of αSMA to chondrocytes than to osteoblasts. This may be due to contribution of muscle-derived cells to chondrocytes (Abou-Khalil et al., 2015; Liu et al., 2011). Notably, αSMA expression is present in some muscle cells, both myogenic satellite cells and a separate population of stromal cells that can contribute to heterotopic bone formation (Matthews et al., 2016; Kolind et al., 2015). It is also possible that there is a periosteal population that is restricted to chondrogenic and perhaps fibroblastic lineages, although we are not aware of any cell populations that have shown these properties in vivo. Gli1+ cells contributed to more chondrocytes than osteoblasts in a similar fracture model (Shi et al., 2017). Gli1+ cells also formed fibrous tissue in fractures that failed to heal due to radiation exposure, while αSMA+ cells did not (Wang et al., 2019b). The contribution of Gli1+ cells in the fibrous layer of the periosteum or from non-periosteal sources could explain these observations. Cxcl12-CreER also appears to identify cells that form mainly chondrocytes in a periosteal fracture callus; however, the presence of labeled cells outside the bone marrow compartment was not reported, making it difficult to identify the tissue where these cells originated (Matsushita et al., 2020).

Lineage tracing with Col2.3CreER indicated that there are cells the periosteum already committed to the osteogenic lineage that participate in osteoblast expansion during fracture healing. The Col2.3CreER+ cells only contribute significantly when tamoxifen is given near the time of fracture, indicating that it does not target long-lived self-renewing progenitors. Cells labeled by Prx1, Col2a1, Cathepsin K, LepR, and Mx1-driven Cres all make a major contribution to fracture callus osteoblasts; however, all these models involve either constitutive Cre expression or long-term labeling that includes differentiated cells, meaning that a combination of stem/progenitor cells and pre-osteoblasts/osteoblasts could be contributing to new osteoblast formation (Debnath et al., 2018; Zhou et al., 2014b; Duchamp de Lageneste et al., 2018; Ortinau et al., 2019; Wang et al., 2019b). Prx1-CreER and Sox9-CreER also label cells in the periosteum that contribute to fracture healing (He et al., 2017; Xiao et al., 2020; Kawanami et al., 2009). It is unclear if these populations overlap with αSMA+ cells or represent separate progenitor pools. However, our results suggest that there is heterogeneity in the cell types contributing to both osteoblasts and chondrocytes in the fracture callus, but by administering tamoxifen at the time of fracture, the majority of periosteal progenitors, particularly those that form osteoblasts, can be genetically targeted by αSMACreER.

αSMA+ cells are present in the endosteum but are less frequent than in the periosteum. αSMA is expressed in a growth-related transiently amplifying progenitor population in the trabecular region. It is also activated in cells in situations where there is a greater need for osteoblasts, such as injury, or genetic ablation of osteoblasts (Matic et al., 2016; Kalajzic et al., 2008). However, it does not appear to label cells involved in bone turnover during adulthood or consistently label long-term injury-responsive cells in the bone marrow compartment. αSMA also labels cells involved transiently in tendon growth as well as identifying progenitors of cementoblasts and periodontal fibroblasts in the periodontium (Dyment et al., 2014; Roguljic et al., 2013). αSMA+ cells are also a major contributor to healing in both these tissues as well as during reparative dentinogenesis following injury in murine molars (Vidovic et al., 2017). There is increasing evidence in bone tissue that different cell types could contribute to homeostasis compared to regeneration. In calvarial bones, the suture provides a niche for cells involved in both growth and healing. Axin2+ and Gli1+ cells contribute to both growth and healing; however, Prx1CreER labels a more restricted population only involved in healing (Zhao et al., 2015; Maruyama et al., 2016; Wilk et al., 2017). Cxcl12+ cells in the central bone marrow do not contribute to osteoblast formation under homeostatic conditions, but are activated and recruited following injury involving the bone marrow (Matsushita et al., 2020). αSMA+ cells in the periosteum appear to represent a population involved primarily in response to injury or increased need for bone formation.

A number of the markers and stains we evaluated did not specifically identify stem cells in the adult periosteum. CD105 was the only marker that was more prevalent in endosteum than periosteum (Tournaire et al., 2020) and was mostly absent from LRCs in all cell types analyzed. CD105+ periosteal cells also failed to form CFU-F. CD105 is best known as one of the markers used to identify mesenchymal stem or stromal cells in culture but has also been used as part of stains to identify or enrich progenitor populations in the periosteum and other tissues (Chan et al., 2015; Debnath et al., 2018; Zhou et al., 2014b; Ortinau et al., 2019). We demonstrated that the BCSP and SSC marker combinations, but not the pre-BSCP population, which were originally developed using neonatal bones, contain a significant proportion of mature osteoblasts when applied to adult bone tissue. This likely influences the observation that these populations expand in settings such as BMP2 treatment or fracture healing, which means some results should be revisited in future when more specific markers for adult bone-residing progenitors are developed (Chan et al., 2015; Marecic et al., 2015). Interestingly, a similar stain was used to enrich periosteal stem cells in a recent study (Debnath et al., 2018). Finally, CD200 was almost universally expressed on mature osteoblasts. A number of studies suggest that CD200 is present in various stem cell populations, but it definitely does not exclude all mature cells, so should be used with caution in settings where mature osteoblasts are present (Chan et al., 2015; Debnath et al., 2018; Lukač et al., 2020).

In conclusion, we have demonstrated that the periosteum is highly enriched with skeletal stem and progenitor cells compared to endosteum and bone marrow. Our data demonstrate the presence of a long-term periosteum residing osteochondroprogenitor cell population that is identified by expression of αSMA, along with injury-activated αSMA+ progenitors and mature Col2.3+ osteoblasts that also contribute to the healing process. Overall, the periosteum contains quite different mesenchymal populations to those present within the bone marrow compartment, and there are multiple populations that directly contribute to new tissue formation during fracture healing.

Materials and methods

Key resources table
Reagent type
(species) or resource
DesignationSource or referenceIdentifiersAdditional information
Genetic reagent (Mus musculus)αSMACreER (Tg(Acta2-cre/ERT2)1Ikal)PMID:22083974MGI:5461154
Genetic reagent (Mus musculus)Col2.3CreER (B6.Cg-Tg(Col1a1-cre/ERT2)1Crm/J)Jackson LabsJax: 016241
Genetic reagent (Mus musculus)Col2.3GFP (B6.Cg-Tg(Col1a1*2.3-GFP)1Rowe/J)PMID:11771662Jax: 013134
Genetic reagent (Mus musculus)Ai9 (B6.Cg-Gt(ROSA)26Sortm9(CAG-tdTomato)Hze/J)Jackson LabsJax: 007909
Genetic reagent (Mus musculus)Rosa-DTA (B6.129P2-Gt(ROSA)26Sortm1(DTA)Lky/J)Jackson LabsJax: 009669
Genetic reagent (Mus musculus)R26-M2rtTA; TetOP-H2B-GFP (B6;129S4-Gt(ROSA)26Sortm1(rtTA*M2)JaeCol1a1tm7(tetO-HIST1H2BJ/GFP)Jae/J)Jackson LabsJax: 016836
Genetic reagent (Mus musculus)NSG (NOD.Cg-Prkdcscid Il2rgtm1Wjl/SzJ)Jackson LabsJax: 005557
Genetic reagent (Mus musculus)PdgfraCreER (B6.129S-Pdgfratm1.1(cre/ERT2)Blh/J)Jackson LabsJax: 032770
AntibodyCD45 eFluor450 (rat monoclonal, clone 30-F11)eBioscience48–0451Flow (1:400)
AntibodyTer119 eFluor450 (rat monoclonal, clone TER-119)eBioscience48–5921Flow (1:200)
AntibodyCD31 eFluor450 (rat monoclonal, clone 390)eBioscience48–0311Flow (1:400)
AntibodyCD45 APC (rat monoclonal, clone 30-F11)eBioscience17–0451Flow (1:400)
AntibodyTer119 APC (rat monoclonal, clone TER-119)eBioscience17–5921Flow (1:200)
AntibodyCD31 APC (rat monoclonal, clone 30-F11)eBioscience17–0311Flow (1:400)
AntibodyCD140a APC (rat monoclonal, clone APA5)eBioscience17–1401Flow (1:100)
AntibodyCD140a PE-Cy7 (rat monoclonal, clone APA5)eBioscience25–1401Flow (1:400)
AntibodyCD140b APC (rat monoclonal, clone APB5)eBioscience17–1402Flow (1:100)
AntibodyCD105 APC (rat monoclonal, clone MJ7/18)eBioscience17–1051Flow (1:200)
AntibodyCD90.2 APC eFluor 780 (rat monoclonal, clone 53–2.1)eBioscience47–0902Flow (1:200)
AntibodyCD90.2 BV605 (rat monoclonal, clone 30-H12)BD Bioscience740334Flow (1:100)
AntibodyCD51 Biotin (rat monoclonal, clone RMV-7)eBioscience13–0512Flow (1:100)
AntibodySca1 AlexaFluor700 (rat monoclonal, clone D7)eBioscience56–5981Flow (1:100)
AntibodyLy51 BV711 (rat monoclonal, clone 6C3)BD Bioscience740691Flow (1:100)
AntibodyCD200 PE (rat monoclonal, clone OX-2)Biolegend123807Flow (1:100)
AntibodySox9 (rabbit polyclonal)EMD MilliporeABE2868IHC (1:200)
AntibodyOsterix (rabbit monoclonal)Abcamab209484IHC (1:400)
AntibodyPDGFRα (rabbit polyclonal)R and D SystemsAF1062IHC (1:80)
AntibodyGoat anti rabbit AlexaFluor647 (polyclonal)ThermoFisherA21244IHC (1:300)
Chemical compound, drugStreptavidin APC eFluor 780eBioscience47–4317Flow (1:400)
Chemical compound, drugStreptavidin PE-Cy7eBioscience25–4317Flow (1:400)
Chemical compound, drugStreptavidin APCeBioscience17–4317Flow (1:400)
Chemical compound, drugTamoxifenSigma AldrichT564875 mg/kg i.p.
Chemical compound, drugDoxycyclineEnvigoTD.01306625 mg/kg in Teklad Custom Diet
Sequence-based reagentBarcode oligo primerIDT5’-AAGCAGTGGTATCAACGC AGAGTACJJJJJJJJJJJJNNNNN NNNTTTTTTTTTTTTTTTTTTT TTTTTTTTTTTVN-3’
Sequence-based reagentCustom primer sequenceIDT5'-AATGATACGGCGACCACC GAGATCTACACGCCTGTCCG CGGAAGCAGTGGTATCAACG CAGAGT*A*C-3’
Sequence-based reagentCustom read one sequenceIDT5’-CGGAAGCAGTGGTATCAA CGCAGAGTAC-3’

Mice

All animal studies were approved by the institutional animal care and use committee at UConn Health or the University of Auckland. Mouse lines are listed in the Key Resources Table. αSMACreER and Col2.3GFP mice have been described (Grcevic et al., 2012; Kalajzic et al., 2002). R26-M2rtTA; TetOP-H2B-GFP dual homozygous mice (Foudi et al., 2009), Ai9 Tomato reporter mice (Madisen et al., 2010), Rosa-DTA (Voehringer et al., 2008), Col2.3CreER (Kim et al., 2004), and PdgfraCreER (Chung et al., 2018) were purchased from Jackson Labs. All these strains were maintained on a C57Bl/6J background. NSG immunodeficient mice (NOD.Cg-Prkdcscid Il2rgtm1Wjl/SzJ, stock # 005557) were a gift from David Rowe. All transgenic mice for experiments were bred at UConn Health and were group housed in individually ventilated cages maintained at 22 ± 2°C, 55 ± 5% humidity, and a 12 hr light-dark cycle and allowed ad libitum access to food and water. Male C57Bl/6J mice were used for experiments at the University of Auckland and were bred at the onsite facility under similar conditions. CreER was activated by tamoxifen in corn oil (75 µg/g i.p.), and two doses of tamoxifen were administered 24 hr apart unless otherwise stated. Animals were assigned to experimental groups in blocks, and cage-mates were generally exposed to the same tamoxifen regimen to avoid inadvertent exposure of low-dose tamoxifen.

Tibia fracture

Request a detailed protocol

Closed tibia fractures stabilized by an intramedullary pin were generated mid-shaft as previously described (Matthews et al., 2014). All mice undergoing fracture lineage tracing were female and were sexually mature young adults (8–16 weeks of age at tamoxifen administration, 8–24 weeks at fracture). The presence of fracture was confirmed by X-ray (Faxitron LX 60). Animals were euthanized and limbs collected for histology 7 or 14 DPF. For secondary fractures, tibia fractures were performed and allowed to heal for 8 weeks with the pin in place. Then a secondary fracture was performed in a similar location, confirmed by X-ray, and collected 7 DPF for histological analysis.

Histology

Request a detailed protocol

Bones were fixed for 5–7 days in 4% paraformaldehyde, incubated in 30% sucrose overnight and embedded in cryomatrix (Fisher Scientific). Cryosections (7 µm) were obtained using a cryostat and tape transfer system (Section-lab, Hiroshima, Japan) as previously described (Dyment et al., 2016). Imaging settings were kept consistent for each dataset. Fluorescent images were obtained with an Axioscan slide scanner (Zeiss). Immunostaining antibodies are listed in the Key Resources Table. For Sox9 and PDGFRα staining, sections were permeabilized in 0.1% Triton X for 15 min, blocked in 2% bovine serum albumin (BSA)/phosphate buffered saline (PBS) for 1 hr, and stained with primary antibody overnight at 4°C. For osterix staining, sections were permeabilized in 0.1% Tween-20 for 15 min, blocked in 10% goat serum, 2% BSA/PBS for 1 hr, stained with anti-Osterix antibody overnight at 4°C. Both sets of slides were washed in 0.1% Tween-20/PBS and stained with secondary antibody for 1 hr at RT, then washed thoroughly. Slides were coverslipped in 50% glycerol/PBS containing 4′,6-diamidino-2-phenylindole (DAPI). Image analysis was performed as described previously using ImageJ (Matthews et al., 2020). Samples were excluded from all image analysis for the following reasons: fracture was not evident or callus was absent. For osteoblast analysis, sections were excluded if they did not contain a clear fracture site. Endosteal analysis was performed only in sections that had a clear central marrow cavity present, and new bone formation evident in this region. Chondrocyte analysis was performed at day 7 only, as most cartilage is gone by day 14, and some fractures were excluded due to lack of cartilage tissue for analysis. Briefly, to evaluate Tom+ contribution to osteoblasts, regions of interest (ROIs) were drawn encompassing the fracture callus. Adjacent ROIs encompassing the central area of the callus (1 mm either side of the fracture site) as well as more distal regions were drawn (up to 6 ROIs/section), then data were pooled to determine total callus values. The DAPI channel was thresholded and separated using the watershed algorithm, then fluorescent signal for each channel was measured in nuclear regions. Standardized thresholds were set manually for each channel. We calculated the number of cells that were GFP+ or Osx+ (osteoblasts), the number that were Tom+, and the percentage of osteoblasts that were Tom+ (dual positive cells/osteoblasts). In most cases, one section/bone was analyzed. Where multiple sections were analyzed, results from different sections were consistently very similar and were pooled. A minimum of 3000 nuclei (median 9262) were evaluated for each sample. Similar analysis was performed on the bone marrow compartment when sections were suitable. For chondrocyte analysis, since there was some non-specific staining with the Sox9 antibody, only regions that demonstrated chondrocytic morphology were incorporated in the ROI analysis.

Cell isolation

Request a detailed protocol

For studies comparing different tissue compartments, bone marrow, endosteum, and periosteum were isolated from hindlimbs of the same mice as previously described (Matic et al., 2016; Wang et al., 2019a). In most cases, two or three animals of the same sex were pooled to generate a sample. Briefly, after removal of muscle tissue, bone marrow was flushed and underwent red blood cell lysis (Matthews et al., 2017). Periosteum was scraped off the bone surface and digested in PBS containing 0.05% Collagenase P and 0.2% hyaluronidase (Sigma Aldrich, St Louis, MO, USA) for 1 hr at 37°C. To isolate endosteal cells, bones were then opened, crushed, loosely attached cells removed by washing, then digested in PBS containing 0.05% Collagenase P for 1 hr at 37°C. Cells were washed and resuspended in staining medium (SM, 2% fetal bovine serum (FBS), 1 mM ethylenediaminetetraacetic acid in PBS).

Flow cytometry and cell sorting

Request a detailed protocol

Cells were stained using antibody cocktails in SM. Antibodies and secondary reagents are listed in the Key Resources Table. Analysis was performed on an LSRII or Aria II instrument (BD Biosciences), and sorting was performed on an Aria II. All experiments included unstained controls for establishing gates derived from non-fluorescent mice, and single color controls generated using periosteal or bone marrow cells. Some experiments included fluorescence minus one (FMO) controls to assist with gating. Dead cells were excluded in most experiments using DAPI staining (~50 ng/ml final concentration). Four-way sorting was performed on selected populations. Cells were sorted into 1.5 ml tubes containing 500 µl collection medium (αMEM 10% FBS). Acquisition and sorting was performed using FACSDiva software and analysis was performed in FlowJo (BD Biosciences). Experiments involved at least two biological replicates per group, and most were performed on more than one occasion to ensure consistent outcomes.

In vitro assays

Request a detailed protocol

CFU-F assays were performed using freshly sorted cells. Live cell numbers were assumed to be 50% of the cells sorted by the machine. Sorted cells were resuspended in αMEM 20% FBS and seeded in six-well plates, generally at 50 cells/cm2, although this varied depending on yield. In most cases, three wells were seeded per population, and at least three independent samples per population were evaluated. Cells were cultured in a 37°C humidified incubator, 5% O2, 5% CO2. Medium was changed on day 4, and cultures were terminated for staining on day 7–8. Cells were fixed in 10% formalin for 10 min, washed in water, then underwent crystal violet staining. For differentiation, after 7 days, medium was changed to an osteogenic/adipogenic combined cocktail containing αMEM 10% FBS, 50 µg/ml ascorbic acid, 4 mM β-glycerophosphate, 1 µM insulin, 0.5 µM rosiglitazone, and transferred to a normoxic incubator. Medium was changed after 3 days and cultures were terminated after 5 days of differentiation. Plates underwent concurrent ALP staining (86R or B5655, Sigma Aldrich) and Oil Red O staining, were imaged, then underwent crystal violet staining.

Single cell RNAseq analysis

Request a detailed protocol

Periosteum was prepared as described above but with a 30-min enzymatic digestion. We sorted Tom+ cells from periosteal isolations from two to three female αSMACreER/Ai9 mice. Cells were isolated from αSMACreER/Ai9/Col2.3GFP mice at 2 days and 6 weeks after tamoxifen and from H2B-GFP/αSMACreER/Ai9 mice 13 weeks after tamoxifen and withdrawal of dox. No staining was performed and dead cells were excluded by DAPI staining. Single Tom+ cells were sorted using the index sort function on the BD Aria II directly into individual wells of a Bio-Rad hard shell 384-well plate. The plate was immediately stored in a −80°C freezer. Approximately half a plate per sample was processed for sequencing (192 cells for 2 days and 13 weeks, 216 cells for 6 weeks). Custom designed Drop-Seq barcoded primers from Integrated DNA Technologies (IDT) were delivered into each well of the 384-well plate, with primers in one well sharing the same unique cell barcode and different UMIs. An Echo 525 liquid handler was used to dispense 1 µl of primers and reaction reagents into each well for cell lysis and cDNA synthesis. Following cDNA synthesis, the contents of each well were collected and pooled into one tube using a Caliper SciClone Liquid Handler. After treatment with exonuclease I to remove unextended primers, the cDNA was polymerase chain reaction (PCR) amplified for 13 cycles. The cDNA was fragmented and amplified for sequencing with the Nextera XT DNA sample prep kit (Illumina) using custom primers that enabled the specific amplification of only the 3′ ends (see Key Resources Table). The libraries were purified, quantified, and sequenced on an Illumina MiSeq. The data preprocessing was performed using Drop-seq tools (version 1.13), with picard (version 2.17.8) for fastq to bam conversion and STAR (version 2.5.4a) for alignment to the reference genome, mm10 (version 1.2.0). Gene expression matrices were generated, and analysis was done using scanpy (version 1.3.7) (Wolf et al., 2018). The data is deposited in GEO under accession GSE165846. The data was filtered to remove cells with fewer than 200 genes, normalized and log transformed. The top 900 highly variable genes as measured by dispersion were used for neighborhood graph generation (nearest neighbor = 20) and dimensionality reduction with uniform manifold approximation and projection (UMAP) (Satija et al., 2015; Becht et al., 2019). Leiden community detection algorithm was used for clustering (Traag et al., 2019). Data were visualized and plots and heat maps generated by CellView (Bolisetty et al., 2017).

Label retention studies

Request a detailed protocol

R26-M2rtTA; TetOP-H2B-GFP (H2B-GFP) males were crossed with αSMACreER/Ai9/Ai9 females to generate experimental mice. Mice received doxin their diet (625 mg/kg) from P10 for 5 weeks. Males were removed from breeding cages prior to birth of pups to ensure mothers did not ingest dox during subsequent pregnancies. Some animals received two doses of i.p. tamoxifen at the time of dox withdrawal. The chase period was a minimum of 13 weeks. Animals that were never exposed to dox were included in all experiments as controls.

Calvarial defect transplantation

Request a detailed protocol

Mice were anesthetized with ketamine/xylazine (135 and 15 mg/kg) and bilateral critical-sized circular defects (3.5 mm diameter) created in the parietal bone as previously described (Gohil et al., 2014). About 2500–5000 sorted αSMA-labeled cells were mixed with 2 × 105 BMSCs, resuspended in 7 µl PBS, loaded on a healos scaffold, and placed into calvaria defects of immunodeficient NSG mice (n = 6). After 5 weeks, mice were euthanized for histological analysis.

Ablation studies

Request a detailed protocol

Animals were generated by crossing Rosa-DTA homozygous mice with αSMACreER, or where lineage tracing was required, crossing αSMACreER/Rosa-DTA with Ai9 and using only Cre+ animals. Fractures were performed in 7–10-week-old animals as previously described (Novak et al., 2020). Tamoxifen was given 0 and 1 DPF in most cases. Briefly, closed transverse diaphyseal femoral fractures were created using a blunt guillotine, after inserting a 24G needle into the intramedullary canal. To evaluate efficacy of the ablation, flow cytometry and histology were performed 4 DPF. For flow, the callus tissue was dissected and single cell suspension from individual male mice (n = 4) was generated, stained with CD45/CD31/Ter119 cocktail and DAPI. For microCT, fractured bones were collected 14 and 21 DPF, fixed in 10% formalin, and scanned (µCT40, Scanco Medical AG, Bassersdorf, Switzerland) with a voxel size of 12 µm, 55 kV, and intensity of 145 µA. We evaluated 150 slices above and below the fracture site to measure density and callus volume, and calculate bone mass. For histological analysis on day 14, frozen sections were stained using von Kossa and callus size and mineralized areas were quantified using ImageJ software as previously described (Novak et al., 2020). An additional group given tamoxifen 0, 2, 4, 6, and 8 DPF were included. The 21 DPF analysis group received tamoxifen 0, 2, 4 DPF. Controls with different tamoxifen regimens showed no difference between groups, so were pooled. Both male and female mice were included (day 14: Cre-: n = 14, 8M 6F; DTA D0,1: n = 10, 4M 6F; DTA D0-8: n = 8, 4M 4F).

Statistics

All experiments include at least three biological replicates. For lineage tracing fracture experiments, we planned to include at least n = 6 per group, although long-term tracing groups in particular had more animals included to account for possible exclusions. We planned to include n = 10 for DTA studies based on our previous experience with quantitative fracture studies, but numbers and sex distribution were limited by the genotypes of available animals. Exact n values are listed in figures or figure legends, or where ranges are stated, in the source data files. Values represent the number of biological replicates (biological replicates represent one mouse/one fracture for injury studies, and one pool of two to three mice for most flow cytometry studies; most of the flow and CFU-F data involves analysis or sorts performed over multiple days or experiments). Graphs show mean ± standard error of the mean. Flow cytometry data on dot plots and in the text show mean ± standard deviation. Statistical analysis was performed in GraphPad Prism, using t-tests or one-way or two-way ANOVAs with appropriate post-hoc tests. For many of the flow cytometry datasets where different populations or tissues from one sample are evaluated, paired tests are used. Details are included in figure legends, and analysis or p values are supplied in the source data files. p<0.05 was used as the threshold for statistical significance.

Data availability

RNAseq data have been deposited in GEO under accession GSE165846. Source data files are provided for all figures (1–8).

The following data sets were generated

References

Decision letter

  1. Subburaman Mohan
    Reviewing Editor; Loma Linda University, United States
  2. Clifford J Rosen
    Senior Editor; Maine Medical Center Research Institute, United States

In the interests of transparency, eLife publishes the most substantive revision requests and the accompanying author responses.

Acceptance summary:

The major findings of the study that fracture healing is contributed by not a single but distinct stem populations, and, that periosteum contains unique progenitor population (αSMA+) that develop into osteoblasts in response to injury, are relevant in terms of identifying effective therapeutic approaches to promote healing of delayed or nonunion fractures.

Decision letter after peer review:

Thank you for submitting your article "Heterogeneity of murine periosteum osteochondroprogenitors involved in fracture healing" for consideration by eLife. Your article has been reviewed by three peer reviewers, one of whom is a member of our Board of Reviewing Editors, and the evaluation has been overseen by Clifford Rosen as the Senior Editor. The reviewers have opted to remain anonymous.

The reviewers have discussed the reviews with one another and the Reviewing Editor has drafted this decision to help you prepare a revised submission. If you have questions about the review and the concerns please contact me directly. I think the biggest issue remains novelty (see below).

As the editors have judged that your manuscript is of interest, but as described below that additional edits are required before it is published, we would like to draw your attention to changes in our revision policy that we have made in response to COVID-19 (https://elifesciences.org/articles/57162). First, because many researchers have temporarily lost access to the labs, we will give authors as much time as they need to submit revised manuscripts. We are also offering, if you choose, to post the manuscript to bioRxiv (if it is not already there) along with this decision letter and a formal designation that the manuscript is "in revision at eLife". Please let us know if you would like to pursue this option. (If your work is more suitable for medRxiv, you will need to post the preprint yourself, as the mechanisms for us to do so are still in development.)

All three reviewers agree that the authors have performed extensive amount of work to support the major conclusion that different osteoprogenitors in the periosteum are involved in fracture healing. However, they all have raised the same concern regarding the limited novelty of the study which is a weakness of this paper. There is now a general agreement in the literature that different cells in the periosteum contribute to the skeletal development. Since the αSMA-labeled cells characterized in this study is considered heterogeneous, the delayed healing in the α-SAM depleted mice is not surprising. It is recommended that you consider performing additional experiments based on the comments of the reviewers provided below and rewrite the manuscript focused on the novel aspects of the paper during the revision.

Reviewer #1:

Much attention has been focused recently in the skeletal biology research in characterizing the skeletal progenitors in the periosteum using various stem cell markers and evaluating their role in fracture healing. These studies have led to identification of a number of self-renewing skeletal stem cell populations (Nestin+; Mx1+αSMA+; Cathepsin K +) in the periosteum (Deveza et al., 2018; Ortinau et al., 2019; Debnath et al., 2018; Tournaire et al., 2020). Furthermore, these studies have identified stem cells in the periosteum to contribute significantly to fracture healing. In this study, Mathews et al. first characterized the putative mesenchymal stem cell markers in the periosteum, endosteum and bone marrow of adult mice. They further used lineage tracing strategy to identify αSMA cells in the periosteum to contribute to approximately 80% of osteoblasts and 40% of chondrocytes following fracture. Based on the finding that ablation of SMA+ cells reduced fracture callus volume and bone mass, it is concluded long-term, slow cycling SMA+ periosteal stem cells are functionally important for bone and cartilage formation during fracture healing. Since identifying cells that contribute to repair of fractured bone is important to the goal of developing therapies to promote healing of nonunion fractures, research in this area is of translational value. Overall, the manuscript is well written and the findings are clearly described. However, there are a number of issues that need the authors attention during the revision of this manuscript.

1) The first part of the studies is focused on characterizing the stem cell surface markers in the different regions of bone, i.e. periosteum, endosteum and bone marrow. These studies seem to clearly show evidence for presence of different populations of stem cells with different osteogenic potential in the different bone compartments. It would have been worthwhile to characterize which of these different populations, for example, Sca1+, CD51+; CD90+; Sca1-, CD51+, CD90-CD105-) are positive for SMA or other known markers of periosteal stem cells (nestin, cathepsin-k). Thus, there seems to be no connection between the first part of the studies and the subsequent studies that are focused on cells expressing SMA.

2) Fracture healing involves both intramembranous and endochondral ossification. Cartilage formation occurs early by chondrocytes which is subsequently converted to bone via endochondral ossification. Periosteal stem cells have been shown to contribute to both formation of chondrocyte-mediated endochondral bone formation as well as osteoblast-mediated intramembranous bone formation. In a recent study, Debnath et al. (Nature 562:133, 2018) showed that cathepsin-k positive periosteal cells contributed to intramembranous bone formation. In addition, recent studies also demonstrate that chondrocytes can directly transdifferentiate into bone forming osteoblasts during postnatal bone formation and fracture repair. Thus, the issues of whether there are different sub populations of stem cells in the periosteum that contribute to chondrocytes and osteoblasts during endochondral bone formation and osteoblasts during intramembranous bone formation and the extent to which the SMA+ stem cells in the periosteum contribute to the different cell types that contribute to the fracture healing remain unresolved.

3) The data in Figure 3 showing that SMA identifies long term progenitor cells in the periosteum seem very interesting. The contribution of SMA+ cells to osteoblasts was evaluated at day 7 after fracture healing. In the closed femur mid-shaft fracture model, cartilage callus is more predominant at day 7. One would, therefore, expect much more chondrocytes at day 7 fracture healing while at later stages around day 10-14, there are more osteoblasts. Accordingly, Figures 3E-G show that while there abundant SMA+ cells in the fracture callus area, the number of osteoblasts seem limited which could be due to the early time point used for histology when one would anticipate higher chondrocytes in the callus area. However, the quantitative data in Figure 3H show that 80% of Tom+ cells are osteoblasts which seem to be inconsistent with the data seen in Figure 3E.

4) To investigate the contribution of the SMA+ cells to fracture healing, the authors have used an established ROSA-DTA mouse model to ablate SMA expressing cells. The fracture phenotype in this model is not adequately characterized. The phenotype is evaluated only at one time point which is day 14 which seems early. Three and four week time points are more appropriate to evaluate if fracture healing is impaired. Time course studies of histological evaluation of cartilage and bone matrix would provide some mechanistic information on whether the reduced callus bone mass is due to reduced cartilage formation, impaired endochondral ossification or defective periosteal bridging.

Reviewer #2:

The authors investigate heterogeneity of periosteal progenitor cells using a multitude of approaches, such as markers for skeletal stem cells, label-retaining assays and lineage tracing studies of cells expressing αSMA-CreER and Col2.3-CreER. They found that SSC markers are particularly enriched in αSMA-labeled label-retaining periosteal cells that contribute to fracture healing. This strength of this study is an extensive application of cutting-edge mouse genetic and flow cytometry technologies to characterize periosteal cells.

1) The novelty of this study is not clear from the manuscript. I see the value of a wealth of important information presented here, including periosteal SSC marker expression, label-retaining cells and SMA-labeled cells in fracture healing. However, most of the major findings have been already reported elsewhere by the authors' group and others. For example, the enrichment of SSCs in the periosteum has been reported by Duchamp de Lageneste et al., 2018, and the role of αSMA+ cells in fracture healing has been reported by the authors' group (Matthews et al., 2014; Grevic et al., 2012). The way it is presented now does not really convey the progress that they made in this study. This is the critical point that the authors need to carefully address in any revision of the manuscript.

2) Heterogeneity of periosteal cells is not experimentally examined in this study, although this point is emphasized in the title. They present many findings without much cohesiveness, though each of which is quite interesting. The authors suggest that αSMA+ cells are composed of short-term and long-term cells that differentially contribute to fracture healing. Single cell RNA-seq analyses of αSMA+ cells would be essential to support this central theme.

3) The authors should present additional histological data to support their important findings. They stated that LRCs were mostly found in the inner cambium layer of the periosteum. However, no high magnification image is provided to support this claim, with appropriate counterstaining. In addition, the authors raised a very interesting theory that PaS cells might be located in the periosteum. They should perform immunostaining for PDGFRa and Sca1, and provide evidence that these cells are indeed present in the periosteum. This will enhance the novelty of this study.

Reviewer #3:

Current studies have evaluated that αSMA+ progenitor cells are long-term label retaining cells and majority of them are present in particularly on periosteum. These cells are activated with injury and become bone forming osteoblasts to contribute to the healing process. Authors provided evidences that αSMA-CreER-labeled cells express skeletal stem cell (SSC) markers and these cells maintain the expression of SSC markers for long period of time. Further, αSMA+ cells labeled at the time of injury make major contribute to bone healing process by turning into osteoblasts. However, the notion of the periosteal stem cell heterogeneity is not new and authors did not show new models, new functions and/or mechanism of αSMA+ population as compared to the authors previous study before. Further, there are many important issues to be addressed.

1) As authors have acknowledged, αSMA+ cells are heterogeneous. In fact, another group already demonstrated that a subset of αSMA-GFP+ cells overlap with Mx1+ cells and these double labeled cells are long-term repopulating stem/progenitor cells on periosteum and contribute to bone healing process with serial transplantation capability. Thus, the study with single αSMA-CreER mouse model may not further specify long-term progenitors within highly heterogenous periosteal cells and in bone healing process following the fracture.

2) LRCs in Figure 2C arrows are not convincing. There are only a few cells in periosteum but majority of GFP+ LRCs are cortical osteocytes in Figure 2C. Therefore, there are high possibility of mixture of these population in GFP high gating. Progenitor marker analysis in Figure 2H and 2I is not clear. What are the cell source for CFU-F assay from Figure 2F and 2G? periosteum? Endosteum?

3) Authors mentioned that αSMA+ cells from the bone marrow in juvenile mice became osteoblasts and osteocytes after 17 days. In Figure S3, Osteocytes are quite obvious to tell from histological section, while osteoblasts are not easy to make this conclusion without additional markers. it could be helpful to use the αSMA-CreER/Col2.3-GFP mice to investigate what proportion of αSMA+ cells turn into osteoblasts under normal bone homeostasis.

4) In previous literature, αSMA-GFP+ cells in periosteum are heterogeneous and only a subset (~25%) of them showed periosteal SSC characteristics (Ortinau et al., 2019). In Figure 3, authors claimed that αSMACreER-labeled cells are long-term osteochondroprogenitor cells. Are all αSMA-labeled cells long-term progenitors? In fact in supplemental figure, authors showed that αSMACreER-labeled cells are weak or no H2B-GFP (LRCs). Please explain this discrepancy.

5) From injury model of Figure 3, authors claimed that αSMA-labeled cells at the time of fracture make a major contribution to newly formed osteoblasts in the marrow compartment. However, it is hard to define the source of these cells as bone marrow αSMA+ cells. First, Figure S3 shows there are some αSMA+ cells in bone marrow regions (Figure S3E). Second, fracture model from current study have potential of having mixture of periosteum and bone marrow cells in both periosteum and bone marrow regions at the fracture site. Therefore, the source of those αSMA+Col2.3GFP+ cells in bone marrow compartment shown in Figure 3E is not clear.

6) In Figure 3F and H, authors showed that αSMA-labeled cells contribute to most fracture callus osteoblasts (~80%) at D+1 tam but only ~20% of osteoblasts at D+90 tam. If the labeled cells are long-term repopulating stem cells, their osteoblast contribution should be stable over the time. These results likely implicate that αSMA-labeled cells are not long-term repopulating progenitors. It will be more convincing if authors trace αSMA-labeled cells even longer time (8-10 months).

7) Figure 4G shows the transplantation of αSMA-labeled cells in a calvarial defect with a few osteoblast differentiation from αSMA+ cells. However, there is no evidence showing that there were no Tom+GFP+ cells from the beginning. Authors should include Tom+GFP+ cell transplantation as a control. Serial transplantation will be more informative.

https://doi.org/10.7554/eLife.58534.sa1

Author response

All three reviewers agree that the authors have performed extensive amount of work to support the major conclusion that different osteoprogenitors in the periosteum are involved in fracture healing. However, they all have raised the same concern regarding the limited novelty of the study which is a weakness of this paper. There is now a general agreement in the literature that different cells in the periosteum contribute to the skeletal development. Since the αSMA-labeled cells characterized in this study is considered heterogeneous, the delayed healing in the αSMA depleted mice is not surprising. It is recommended that you consider performing additional experiments based on the comments of the reviewers provided below and rewrite the manuscript focused on the novel aspects of the paper during the revision.

We have revised the manuscript extensively. While the idea of heterogeneity of progenitors is becoming more commonplace, this is a recent development, and a lot of the literature, even in the past 10 years, appears to assume that there is one stem cell population that gives rise to osteoblasts in animals at various ages, and in settings including homeostasis and injury. Our study presents novel in vivo experimental data that proves there are separate populations contributing to fracture healing, and strongly suggests that the periosteum contains unique progenitor populations. We therefore believe it is an important contribution to the field.

Reviewer #1:

Much attention has been focused recently in the skeletal biology research in characterizing the skeletal progenitors in the periosteum using various stem cell markers and evaluating their role in fracture healing. These studies have led to identification of a number of self-renewing skeletal stem cell populations (Nestin+; Mx1+αSMA+; Cathepsin K +) in the periosteum (Deveza et al., 2018; Ortinau et al., 2019; Debnath et al., 2018; Tournaire et al., 2020). Furthermore, these studies have identified stem cells in the periosteum to contribute significantly to fracture healing. In this study, Mathews et al. first characterized the putative mesenchymal stem cell markers in the periosteum, endosteum and bone marrow of adult mice. They further used lineage tracing strategy to identify αSMA cells in the periosteum to contribute to approximately 80% of osteoblasts and 40% of chondrocytes following fracture. Based on the finding that ablation of SMA+ cells reduced fracture callus volume and bone mass, it is concluded long-term, slow cycling SMA+ periosteal stem cells are functionally important for bone and cartilage formation during fracture healing. Since identifying cells that contribute to repair of fractured bone is important to the goal of developing therapies to promote healing of nonunion fractures, research in this area is of translational value. Overall, the manuscript is well written and the findings are clearly described. However, there are a number of issues that need the authors attention during the revision of this manuscript.

1) The first part of the studies is focused on characterizing the stem cell surface markers in the different regions of bone, i.e. periosteum, endosteum and bone marrow. These studies seem to clearly show evidence for presence of different populations of stem cells with different osteogenic potential in the different bone compartments. It would have been worthwhile to characterize which of these different populations, for example, Sca1+, CD51+; CD90+; Sca1-, CD51+, CD90-CD105-) are positive for SMA or other known markers of periosteal stem cells (nestin, cathepsin-k). Thus, there seems to be no connection between the first part of the studies and the subsequent studies that are focused on cells expressing SMA.

We have reorganized the presentation of the data to improve the flow and the integration of the αSMA data with the rest of the results focusing on characterizing periosteal progenitor populations. We have added quantification of αSMA+ cells within the Sca1/CD51 populations, and CD90/CD105 populations in Figure 2—figure supplement 2D, Results subsection “Putative mesenchymal stem cell markers are enriched in periosteum”. Since Nestin and CatK are not cell surface markers that can be evaluated without transgenic lines, we were not able to assess the overlap of these populations. We have, however, included all these markers, and various others used in recent studies, in the evaluation of αSMA+ cells’ gene expression (Figure 3D). CatK expression was detected in 79/529 αSMA+ cells. Most of these cells (90%) were within the progenitor and fibroblast-like clusters (#1 and 2). Nestin was detected in 15/529 cells, the majority (67%) of which were in the perivascular-like cluster 3.

2) Fracture healing involves both intramembranous and endochondral ossification. Cartilage formation occurs early by chondrocytes which is subsequently converted to bone via endochondral ossification. Periosteal stem cells have been shown to contribute to both formation of chondrocyte-mediated endochondral bone formation as well as osteoblast-mediated intramembranous bone formation. In a recent study, Debnath et al. (Nature 562:133, 2018) showed that cathepsin-k positive periosteal cells contributed to intramembranous bone formation. In addition, recent studies also demonstrate that chondrocytes can directly transdifferentiate into bone forming osteoblasts during postnatal bone formation and fracture repair. Thus, the issues of whether there are different sub populations of stem cells in the periosteum that contribute to chondrocytes and osteoblasts during endochondral bone formation and osteoblasts during intramembranous bone formation and the extent to which the SMA+ stem cells in the periosteum contribute to the different cell types that contribute to the fracture healing remain unresolved.

We have performed a much more thorough characterization and quantification of αSMA+ cell contribution to fracture healing than other studies in the literature that have characterized other progenitor markers that we are aware of. While Debnath et al. report that CatK+ progenitors show intramembranous bone formation upon transplantation, they demonstrate some contribution to fibrocartilage following fracture in vivo (see extended data Figure 8 in their publication) (Debnath et al., 2018). The CatK model is problematic for lineage tracing since it is a constitutive Cre that doesn’t enable temporal understanding of the CatK expression. Reaching a consensus on mesenchymal/osteoprogenitor cell identity will require the research community to compile data from numerous lineage tracing models to confirm contribution of different cell populations. There is a great deal more to understand, but we believe that our data using αSMACreER and Col2.3CreER clearly demonstrate that there are different populations contributing to healing, and some of these populations are lineage-restricted. We therefore think this is an important contribution to the field that can be considered in the context of other results.

3) The data in Figure 3 showing that SMA identifies long term progenitor cells in the periosteum seem very interesting. The contribution of SMA+ cells to osteoblasts was evaluated at day 7 after fracture healing. In the closed femur mid-shaft fracture model, cartilage callus is more predominant at day 7. One would, therefore, expect much more chondrocytes at day 7 fracture healing while at later stages around day 10-14, there are more osteoblasts. Accordingly, Figures 3E-G show that while there abundant SMA+ cells in the fracture callus area, the number of osteoblasts seem limited which could be due to the early time point used for histology when one would anticipate higher chondrocytes in the callus area. However, the quantitative data in Figure 3H show that 80% of Tom+ cells are osteoblasts which seem to be inconsistent with the data seen in Figure 3E.

The data the reviewer is referring to (now Figure 1E-I) actually present the percent of Col2.3GFP+ osteoblasts that are Tom+ (dual positive cells/GFP+) rather than the proportion of SMA cells differentiating into osteoblasts. The reviewer’s comment about the dynamics of osteoblast and chondrocyte numbers in the callus are correct – there are generally fewer osteoblasts at day 7, and almost no chondrocytes by day 14, which is why this cell type was only analyzed in the day 7 callus. We analyzed both time points for osteoblasts as we hypothesized that there could be different cells contributing to osteoblasts at different phases of healing. Our data do not support this hypothesis, and we restricted certain analyses such as the 90 day chase to day 7 only. The different time points are now mentioned in the Results section: “The contribution of αSMA-labeled cells was similar in both the initial intramembranous bone formation in the day 7 callus, and the more fully mineralized day 14 callus.” We did not present additional data about the fracture composition as it is well-established, and we felt it made it more difficult to clearly present our findings. However, we clarified the parameter that is presented in both the Materials and methods: We calculated the number of cells that were GFP+ or Osx+ (osteoblasts), the number that were Tom+, and the percentage of osteoblasts that were Tom+ (dual positive cells/osteoblasts).”; and the figure legend for Figure 1. In addition, we have added much more detailed source data for Figure 1G-I which shows all the key outputs from the image analysis. By comparing the GFP+ values in column N for day 7 and day 14, it is clear that the proportion of osteoblasts in the callus tends to increase over the time course. We hope this will make it clearer to readers what analysis is presented.

4) To investigate the contribution of the SMA+ cells to fracture healing, the authors have used an established ROSA-DTA mouse model to ablate SMA expressing cells. The fracture phenotype in this model is not adequately characterized. The phenotype is evaluated only at one time point which is day 14 which seems early. Three and four week time points are more appropriate to evaluate if fracture healing is impaired. Time course studies of histological evaluation of cartilage and bone matrix would provide some mechanistic information on whether the reduced callus bone mass is due to reduced cartilage formation, impaired endochondral ossification or defective periosteal bridging.

We performed the DTA ablation studies as a proof of principle that αSMA+ cells have a functional role in the fracture healing process, but limited the evaluation as this is not a physiologically relevant model, and there are potentially other effects of ablating αSMA+ cells that affect the health of the animals. We have added additional data to what is now Figure 7, and the accompanying Results section. We have confirmed that the callus area and mineralized area are reduced histologically at 14 DPF, and now present microCT data from 21 DPF. The phenotype at day 21 is more subtle, but still indicates a deficiency in fracture healing. The following statement has been added to the results to describe this process: “Overall, our data demonstrate that partial ablation of αSMA+ progenitors reduced callus formation and delayed mineralization. The less severe phenotype at 21 DPF suggests this may be a delay in healing rather than complete disruption of the process, although the fractures from mice with αSMA-ablation probably never reach the size of control calluses. Healing in the DTA mice may still progress due to compensation of αSMA+ progenitors that escaped ablation, or larger contribution from αSMA-negative progenitors in this setting.”

Reviewer #2:

The authors investigate heterogeneity of periosteal progenitor cells using a multitude of approaches, such as markers for skeletal stem cells, label-retaining assays and lineage tracing studies of cells expressing αSMA-CreER and Col2.3-CreER. They found that SSC markers are particularly enriched in αSMA-labeled label-retaining periosteal cells that contribute to fracture healing. This strength of this study is an extensive application of cutting-edge mouse genetic and flow cytometry technologies to characterize periosteal cells.

1) The novelty of this study is not clear from the manuscript. I see the value of a wealth of important information presented here, including periosteal SSC marker expression, label-retaining cells and SMA-labeled cells in fracture healing. However, most of the major findings have been already reported elsewhere by the authors' group and others. For example, the enrichment of SSCs in the periosteum has been reported by Duchamp de Lageneste et al., 2018, and the role of αSMA+ cells in fracture healing has been reported by the authors' group (Matthews et al., 2014; Grevic et al., 2012). The way it is presented now does not really convey the progress that they made in this study. This is the critical point that the authors need to carefully address in any revision of the manuscript.

We have reorganized the manuscript in a way that better highlights the important findings and improves the cohesiveness. While some of the observations have been reported in other studies, we believe that it is important to have multiple lines of evidence from a variety of sources to support the presence and definition of different populations of stem/progenitor cells within the periosteum. More importantly, while other studies imply heterogeneity of healing-related progenitors, our study shows direct evidence of this through in vivo lineage tracing of different populations.

2) Heterogeneity of periosteal cells is not experimentally examined in this study, although this point is emphasized in the title. They present many findings without much cohesiveness, though each of which is quite interesting. The authors suggest that αSMA+ cells are composed of short-term and long-term cells that differentially contribute to fracture healing. Single cell RNA-seq analyses of αSMA+ cells would be essential to support this central theme.

Flow cytometry data indicates heterogeneity with respect to surface marker expression, however we have extended this by adding single cell RNAseq data for αSMA+ cells to the revised manuscript (Figure 3 and supplements). This data has indicated heterogeneity in the population, and shifts in relative abundance of the populations over time. Further studies will be required, potentially with the use of new genetic tools, to confirm the roles of these subpopulations, but we believe this new data strengthens our findings.

3) The authors should present additional histological data to support their important findings. They stated that LRCs were mostly found in the inner cambium layer of the periosteum. However, no high magnification image is provided to support this claim, with appropriate counterstaining. In addition, the authors raised a very interesting theory that PaS cells might be located in the periosteum. They should perform immunostaining for PDGFRa and Sca1, and provide evidence that these cells are indeed present in the periosteum. This will enhance the novelty of this study.

High magnification of H2B-GFP expression in the periosteum is shown in what is now Figure 4—figure supplement 1. We have added additional annotation to this figure to indicate labeled cells within the inner layer of the periosteum. Unfortunately, we don’t have counterstained images available, but we believe that the DAPI staining clearly demonstrates the location within the periosteum.

We have added additional data addressing the use of PDGFRα, and this is presented in Figure 8 and the accompanying Results section. Both immunostaining and Cre-mediated labeling indicated that PDGFRα expression is very prevalent in the periosteum, likely more so than our flow data suggested. Unfortunately we were unable to achieve dual Sca1/PDGFRα staining of suitable quality for publication in the current circumstances. However, the new data we present indicates serious caveats of using PDGFRα as an osteoprogenitor marker, primarily due to clear expression in the osteoblast layer, that we thought were important to include, and we have altered the Results and Discussion accordingly.

Reviewer #3:

Current studies have evaluated that αSMA+ progenitor cells are long-term label retaining cells and majority of them are present in particularly on periosteum. These cells are activated with injury and become bone forming osteoblasts to contribute to the healing process. Authors provided evidences that αSMA-CreER-labeled cells express skeletal stem cell (SSC) markers and these cells maintain the expression of SSC markers for long period of time. Further, αSMA+ cells labeled at the time of injury make major contribute to bone healing process by turning into osteoblasts. However, the notion of the periosteal stem cell heterogeneity is not new and authors did not show new models, new functions and/or mechanism of αSMA+ population as compared to the authors previous study before. Further, there are many important issues to be addressed.

1) As authors have acknowledged, αSMA+ cells are heterogeneous. In fact, another group already demonstrated that a subset of αSMA-GFP+ cells overlap with Mx1+ cells and these double labeled cells are long-term repopulating stem/progenitor cells on periosteum and contribute to bone healing process with serial transplantation capability. Thus, the study with single αSMA-CreER mouse model may not further specify long-term progenitors within highly heterogenous periosteal cells and in bone healing process following the fracture.

We have taken advantage of the inducible Cre model that allows us to characterize cells ex vivo and track their potential in vivo in a physiological setting. The data from other studies including the Ortinau study, and Debnath et al. 2018 (Debnath et al., 2018; Ortinau et al., 2019), as well as our data, suggest that transplantation results do not completely recapitulate in vivo differentiation potential as we generally see osteogenic differentiation only in following transplantation, while observing osteochondrogenic potential in the fracture setting. In the revised manuscript, we have changed the order to emphasize novel findings, added new data to confirm self-renewal, as well as single cell RNAseq data to further characterize the αSMA+ population.

2) LRCs in Figure 2C arrows are not convincing. There are only a few cells in periosteum but majority of GFP+ LRCs are cortical osteocytes in Figure 2C. Therefore, there are high possibility of mixture of these population in GFP high gating. Progenitor marker analysis in Figure 2H and 2I is not clear. What are the cell source for CFU-F assay from Figure 2F and 2G? periosteum? Endosteum?

We have annotated the H2B-GFP histology more clearly in what is now Figure 4—figure supplement 1.

We are confident that our cell preparation method for the periosteum does not result in osteocyte contamination as we scrape the periosteum surface and do not digest cortical bone. Our endosteal preparation involves digestion of bone pieces, however we have previously demonstrated using Dmp1CreER/Ai9 mice, which show reporter expression only in osteocytes in the absence of tamoxifen treatment (Tam-), that we get almost no Tom+ cells in these preparations (see Matic, Matthews et al. 2016, Supplemental Figure S6C)

The following statement has been added to the Results: “Our periosteal isolation procedure excludes bone tissue and therefore osteocytes, and we have previously demonstrated that our endosteal preparations contain almost no osteocytes.”

We have added additional annotation to the figure (now Figure 4), to clarify that the CFU-F data is from periosteal cells, and the populations analyzed in 4F and G are spelt out more clearly in the figure legend.

3) Authors mentioned that αSMA+ cells from the bone marrow in juvenile mice became osteoblasts and osteocytes after 17 days. In Figure S3, Osteocytes are quite obvious to tell from histological section, while osteoblasts are not easy to make this conclusion without additional markers. it could be helpful to use the αSMA-CreER/Col2.3-GFP mice to investigate what proportion of αSMA+ cells turn into osteoblasts under normal bone homeostasis.

This data is replication of data we published in 2012 (Grcevic et al., 2012) that we included for completeness, and as we had not previously published the longer term tracing data. We have therefore not completed this study as our data clearly indicate contribution to osteoblasts and osteocytes, and it is outside the major scope of this study.

4) In previous literature, αSMA-GFP+ cells in periosteum are heterogeneous and only a subset (~25%) of them showed periosteal SSC characteristics (Ortinau et al., 2019). In Figure 3, authors claimed that αSMACreER-labeled cells are long-term osteochondroprogenitor cells. Are all αSMA-labeled cells long-term progenitors? In fact in supplemental Figure, authors showed that αSMACreER-labeled cells are weak or no H2B-GFP (LRCs). Please explain this discrepancy.

We have noted at various points of the manuscript that αSMA-labeled cells are heterogeneous, and the fracture lineage tracing data (now Figure 1) suggests that only a subset have long-term progenitor characteristics. This is now more clearly demonstrated using single cell RNAseq data presented in the new Figure 3. As mentioned in the results, over 20% of αSMA-labeled cells were LRCs, and we have added additional annotation to what is now Figure 4—figure supplement 1 to show this histologically as well.

5) From injury model of Figure 3, authors claimed that αSMA-labeled cells at the time of fracture make a major contribution to newly formed osteoblasts in the marrow compartment. However, it is hard to define the source of these cells as bone marrow αSMA+ cells. First, Figure S3 shows there are some αSMA+ cells in bone marrow regions (Figure S3E). Second, fracture model from current study have potential of having mixture of periosteum and bone marrow cells in both periosteum and bone marrow regions at the fracture site. Therefore, the source of those αSMA+Col2.3GFP+ cells in bone marrow compartment shown in Figure 3E is not clear.

It is true that we cannot completely exclude some bone marrow contribution in the fracture model, however the weight of evidence in the literature indicates that when present, the periosteum is the major contributor to fracture callus formation. We have also performed preliminary studies with a periosteal scratch model where the cortex is not injured, and demonstrated formation of a callus containing αSMA-labeled cells (data not shown). In addition, there is clear evidence indicating that bone marrow/endosteum-derived cells do not form fibrocartilage in vivo, and we see fibrocartilage near the fracture site in most day 7 samples (Colnot et al., 2009). On this basis, combined with our data indicating αSMA+ cells are present but much rarer in the endosteum (new Figure 2D), we have persisted with our interpretation that the cells contributing to the fracture callus osteoblasts and chondrocytes are primarily derived from the periosteum. The majority of the bone formation in the bone marrow compartment is around the pin and is usually in the proximal area of the bone at least 1mm from the fracture site. We also never see evidence of chondrocytes. On this basis, we are fairly confident that bone marrow-resident cells are responsible for this bone formation. The text has been modified to clarify this:

“Following the pin insertion and fracture, bone formation also occurs inside the marrow compartment, particularly in the proximal diaphysis surrounding the pin over 1mm from the fracture site (Figure 1—figure supplement 3).”

6) In Figure 3F and H, authors showed that αSMA-labeled cells contribute to most fracture callus osteoblasts (~80%) at D+1 tam but only ~20% of osteoblasts at D+90 tam. If the labeled cells are long-term repopulating stem cells, their osteoblast contribution should be stable over the time. These results likely implicate that αSMA-labeled cells are not long-term repopulating progenitors. It will be more convincing if authors trace αSMA-labeled cells even longer time (8-10 months).

We agree with the reviewer than it would be interesting to do even longer tracing, but given the amount of time this study would take, it is outside the scope of revisions for this journal. As explained in the manuscript, the difference between the labeling at the time of fracture, and 90 days before is likely to be primarily due to cells that activate αSMA expression soon after injury. It is also notable that contribution stabilizes between 42 and 90 days chase.

7) Figure 4G shows the transplantation of αSMA-labeled cells in a calvarial defect with a few osteoblast differentiation from αSMA+ cells. However, there is no evidence showing that there were no Tom+GFP+ cells from the beginning. Authors should include Tom+GFP+ cell transplantation as a control. Serial transplantation will be more informative.

Almost all the αSMA-Tom+ cells isolated for transplantation were GFP-. This is shown in the flow plots from the cell sorting in Author response image 1. On this basis, it is not feasible to sort and transplant Tom+GFP+ cells as it is such a rare population. In other studies, Col2.3GFP+ cells show minimal growth following sorting, at least in vitro (Siclari et al., 2013).

Author response image 1
Dot plot of periosteum sample used for sorting αSMA-labeled cells for calvarial defect transplantation.

Note that almost all Tom+ cells are GFP-.

We agree with the reviewer that self-renewal was not adequately assessed. We have included a different experiment to address self-renewal of these cells in the new Figure 6 where two cycles of fractures were performed (tamoxifen labeling, initial fracture, healing for 8 weeks, secondary fracture at same site). This shows significant contributions of initially αSMA-labeled cells to secondary fracture healing, even if the initial tamoxifen labeling was 14 days before fracture. This strongly suggests identification of at least a subset of self-renewing progenitors. We believe this experimental design is more physiologically relevant, particularly as there are many technical challenges involved in serial transplantation experiments.

https://doi.org/10.7554/eLife.58534.sa2

Article and author information

Author details

  1. Brya G Matthews

    1. Department of Molecular Medicine and Pathology, University of Auckland, Auckland, New Zealand
    2. Department of Reconstructive Sciences, UConn Health, Farmington, United States
    Contribution
    Conceptualization, Formal analysis, Supervision, Funding acquisition, Validation, Investigation, Visualization, Methodology, Writing - original draft, Project administration
    For correspondence
    brya.matthews@auckland.ac.nz
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-4145-4696
  2. Sanja Novak

    Department of Reconstructive Sciences, UConn Health, Farmington, United States
    Contribution
    Formal analysis, Investigation, Methodology, Writing - review and editing
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-8042-932X
  3. Francesca V Sbrana

    Department of Reconstructive Sciences, UConn Health, Farmington, United States
    Contribution
    Formal analysis, Investigation
    Competing interests
    No competing interests declared
  4. Jessica L Funnell

    Department of Reconstructive Sciences, UConn Health, Farmington, United States
    Contribution
    Formal analysis
    Competing interests
    No competing interests declared
  5. Ye Cao

    Department of Molecular Medicine and Pathology, University of Auckland, Auckland, New Zealand
    Contribution
    Validation, Investigation
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-1654-9060
  6. Emma J Buckels

    Department of Molecular Medicine and Pathology, University of Auckland, Auckland, New Zealand
    Contribution
    Formal analysis
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-9256-5014
  7. Danka Grcevic

    1. Department of Physiology and Immunology, University of Zagreb, Zagreb, Croatia
    2. Croatian Intitute for Brain Research, University of Zagreb, Zagreb, Croatia
    Contribution
    Investigation, Methodology, Writing - review and editing
    Competing interests
    No competing interests declared
  8. Ivo Kalajzic

    Department of Reconstructive Sciences, UConn Health, Farmington, United States
    Contribution
    Conceptualization, Supervision, Funding acquisition, Investigation, Methodology, Writing - review and editing
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-6752-3442

Funding

Connecticut Innovations (14-SCA-UCHC-02)

  • Brya G Matthews

Health Research Council of New Zealand (Sir Charles Hercus Fellowship)

  • Brya G Matthews

American Society for Bone and Mineral Research (Rising Star Award)

  • Brya G Matthews

National Institute of Arthritis and Musculoskeletal and Skin Diseases (AR055607)

  • Ivo Kalajzic

National Institute of Arthritis and Musculoskeletal and Skin Diseases (AR070813)

  • Ivo Kalajzic

Connecticut Innovations (16-RMB-UCHC-10)

  • Ivo Kalajzic

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

This work has been supported by Connecticut Stem Cell grant 14-SCA-UCHC-02, the Health Research Council of New Zealand Sir Charles Hercus Fellowship, and the American Society for Bone and Mineral Research Rising Star Award to BGM; Connecticut Regenerative Medicine Research Fund grant 16-RMB-UCHC-10 and NIH/NIAMS grants AR055607 and AR070813 to IK. We gratefully acknowledge the contribution of the Single Cell Biology Lab at The Jackson Laboratory for Genomic Medicine for their expert assistance with the work described in this publication. We thank Prof. David Rowe for providing NSG mice. We thank Drs Emilie Roeder and India Azevedo for assistance with fracture studies, and Dr Liping Wang for performing some of the calvarial defect surgeries.

Ethics

Animal experimentation: The majority of the study was performed at UConn Health in an AAALAC accredited facility in accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health. Studies were approved by the UConn Health institutional animal care and use committee (IACUC) under protocol numbers 100490-0815, 101095-0518, 101757-0221, Hz#-Dox0322e-101058 and Hz#-MCh0331e-101086. Experiments at the University of Auckland were performed in accordance with the University of Auckland Code of Ethical Conduct (CEC) and the Animal Welfare Act 1999, under Animal Ethical Committee approval 001940.

Senior Editor

  1. Clifford J Rosen, Maine Medical Center Research Institute, United States

Reviewing Editor

  1. Subburaman Mohan, Loma Linda University, United States

Version history

  1. Received: May 4, 2020
  2. Accepted: February 8, 2021
  3. Accepted Manuscript published: February 9, 2021 (version 1)
  4. Version of Record published: February 25, 2021 (version 2)

Copyright

© 2021, Matthews et al.

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 2,819
    Page views
  • 551
    Downloads
  • 34
    Citations

Article citation count generated by polling the highest count across the following sources: Crossref, Scopus, PubMed Central.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Brya G Matthews
  2. Sanja Novak
  3. Francesca V Sbrana
  4. Jessica L Funnell
  5. Ye Cao
  6. Emma J Buckels
  7. Danka Grcevic
  8. Ivo Kalajzic
(2021)
Heterogeneity of murine periosteum progenitors involved in fracture healing
eLife 10:e58534.
https://doi.org/10.7554/eLife.58534

Further reading

    1. Medicine
    2. Stem Cells and Regenerative Medicine
    Thiruma V Arumugam, Asfa Alli-Shaik ... Jayantha Gunaratne
    Research Article

    Intermittent fasting (IF) has been shown to reduce cardiovascular risk factors in both animals and humans, and can protect the heart against ischemic injury in models of myocardial infarction. However, the underlying molecular mechanisms behind these effects remain unclear. To shed light on the molecular and cellular adaptations of the heart to IF, we conducted comprehensive system-wide analyses of the proteome, phosphoproteome, and transcriptome, followed by functional analysis. Using advanced mass spectrometry, we profiled the proteome and phosphoproteome of heart tissues obtained from mice that were maintained on daily 12- or 16 hr fasting, every-other-day fasting, or ad libitum control feeding regimens for 6 months. We also performed RNA sequencing to evaluate whether the observed molecular responses to IF occur at the transcriptional or post-transcriptional levels. Our analyses revealed that IF significantly affected pathways that regulate cyclic GMP signaling, lipid and amino acid metabolism, cell adhesion, cell death, and inflammation. Furthermore, we found that the impact of IF on different metabolic processes varied depending on the length of the fasting regimen. Short IF regimens showed a higher correlation of pathway alteration, while longer IF regimens had an inverse correlation of metabolic processes such as fatty acid oxidation and immune processes. Additionally, functional echocardiographic analyses demonstrated that IF enhances stress-induced cardiac performance. Our systematic multi-omics study provides a molecular framework for understanding how IF impacts the heart’s function and its vulnerability to injury and disease.

    1. Developmental Biology
    2. Stem Cells and Regenerative Medicine
    Tao Zhang, Liyang Wan ... Hongbin Lu
    Research Article Updated

    The attachment site of the rotator cuff (RC) is a classic fibrocartilaginous enthesis, which is the junction between bone and tendon with typical characteristics of a fibrocartilage transition zone. Enthesis development has historically been studied with lineage tracing of individual genes selected a priori, which does not allow for the determination of single-cell landscapes yielding mature cell types and tissues. Here, in together with open-source GSE182997 datasets (three samples) provided by Fang et al., we applied Single-cell RNA sequencing (scRNA-seq) to delineate the comprehensive postnatal RC enthesis growth and the temporal atlas from as early as postnatal day 1 up to postnatal week 8. And, we furtherly performed single-cell spatial transcriptomic sequencing on postnatal day 1 mouse enthesis, in order to deconvolute bone-tendon junction (BTJ) chondrocytes onto spatial spots. In summary, we deciphered the cellular heterogeneity and the molecular dynamics during fibrocartilage differentiation. Combined with current spatial transcriptomic data, our results provide a transcriptional resource that will support future investigations of enthesis development at the mechanistic level and may shed light on the strategies for enhanced RC healing outcomes.