Programmed conversion of hypertrophic chondrocytes into osteoblasts and marrow adipocytes within zebrafish bones
Abstract
Much of the vertebrate skeleton develops from cartilage templates that are progressively remodeled into bone. Lineage tracing studies in mouse suggest that chondrocytes within these templates persist and become osteoblasts, yet the underlying mechanisms of this process and whether chondrocytes can generate other derivatives remain unclear. We find that zebrafish cartilages undergo extensive remodeling and vascularization during juvenile stages to generate fat-filled bones. Growth plate chondrocytes marked by sox10 and col2a1a contribute to osteoblasts, marrow adipocytes, and mesenchymal cells within adult bones. At the edge of the hypertrophic zone, chondrocytes re-enter the cell cycle and express leptin receptor (lepr), suggesting conversion into progenitors. Further, mutation of matrix metalloproteinase 9 (mmp9) results in delayed growth plate remodeling and fewer marrow adipocytes. Our data support Mmp9-dependent growth plate remodeling and conversion of chondrocytes into osteoblasts and marrow adipocytes as conserved features of bony vertebrates.
https://doi.org/10.7554/eLife.42736.001eLife digest
Our adult bones are made of a fatty tissue, called marrow, wrapped inside a hard outer layer produced by bone cells. They may appear stiff and unyielding, but our bones are actually dynamic structures. Early in life, most bones start as small ‘templates’ made of another, flexible tissue called cartilage. As the templates grow into adult bones, the cartilage is gradually replaced by bone and fat, but this process is still poorly understood. For example, it is not clear whether cartilage cells simply die and make way for new cells, or instead if they turn into bone and fat cells. To investigate this question, Giovannone, Paul et al. set out to follow the fate of early cartilage cells in zebrafish, and to compare this with what happens in mammals. Zebrafish were chosen because their skeleton and ours develop in similar ways; yet, these animals are much easier to study, in particular because their embryos are transparent.
Young cartilage cells were ‘tagged’ with a long-lasting fluorescent protein in genetically engineered zebrafish embryos, and then followed over time. As the embryos started to form bones, the cartilage cells gave rise to bone cells, fat cells, and also potentially adult stem cells within the marrow, which can become other types of cells. This process required a protein called Mmp9, which also helps shape bone development in other organisms, including humans.
Defects in how early cartilage templates morph into bone and fat may contribute to dwarfism and other severe conditions. Fully grasping the molecular mechanisms that preside over this complex transition may one day help design drugs to treat skeletal disorders.
https://doi.org/10.7554/eLife.42736.002Introduction
Vertebrate bones develop via two largely distinct processes. Intramembranous (i.e. dermal) bone, which makes up a large portion of the skull, arises through the direct differentiation of mesenchymal precursors into osteoblasts and then osteocytes. In contrast, endochondral bone, which comprises the majority of the axial and limb skeletons, arises through the progressive remodeling of an embryonic cartilage template. On the outside of developing endochondral bone, perichondral cells mature into periosteal progenitors that contribute to the bone collar. The cartilage templates of endochondral bone are organized into distinct zones of chondrocytes: resting, proliferative, pre-hypertrophic, and hypertrophic. In mammals, chondrocytes at the edge of the developing hypertrophic zone largely disappear as the cartilage matrix is degraded, a process concurrent with the invasion of blood vessels, hematopoietic cells, and progenitors for osteoblasts and marrow adipocytes (Maes et al., 2010). This growth plate remodeling contributes to the establishment of trabecular bone, complementing the cortical bone largely derived from the periosteum, and the marrow cavity supports continued hematopoiesis. As with mammals, zebrafish also have intramembranous and endochondral bones. Their endochondral bones are hollow and filled predominantly with fat yet do not support hematopoiesis as in mammals (Witten and Huysseune, 2009; Weigele and Franz-Odendaal, 2016). It has remained unclear, however, whether zebrafish bone arises solely through osteoblast differentiation in the periosteum, or also through invasion of the vasculature and conversion of growth plate cartilage to bone as in mammals. The source of marrow adipocytes also remains unclear in either fish or mammals.
It has long been appreciated that many hypertrophic chondrocytes undergo cell death during endochondral ossification, with osteoblasts forming from periosteal cells brought into the bone along with the vasculature (Maes et al., 2010). At the same time, there are numerous studies showing that cultured chondrocytes can dedifferentiate into mesenchymal progenitors and/or transdifferentiate into osteoblasts (Shimomura et al., 1975; Mayne et al., 1976; von der Mark and von der Mark, 1977). Recent lineage tracing studies of hypertrophic chondrocytes, using constitutive and inducible Col10a1-Cre- and Aggrecan-Cre-based transgenes in mice, has revealed that such transdifferentiation may also occur in vivo, with chondrocytes making a major contribution to osteoblasts within trabecular bone and potentially also the bone collar (Yang et al., 2014; Kobayashi et al., 2014; Jing et al., 2015; Park et al., 2015). A limitation of these studies is the use of population-based labeling by Cre recombination, which cannot exclude low-level and/or leaky labeling of other cell types. It is also unclear whether hypertrophic chondrocytes can give rise to other cell types, such as marrow adipocytes, and whether hypertrophic chondrocytes directly transform into osteoblasts or do so through a stem cell intermediate. Finally, it is unknown whether the ability of chondrocytes to generate osteoblasts and other cell types is specific to mammals or a more broadly shared feature of vertebrates.
In this study, we address the long-term fate of growth plate chondrocytes in zebrafish, as well as potential mechanisms of their fate plasticity. We use the ceratohyal (Ch) bone of the lower face as a model. This long bone, which is derived from cranial neural crest cells, exhibits properties in common with the long bones of mammalian limbs, including two prominent growth plates at either end and a marrow cavity (Paul et al., 2016). Here, we describe remodeling of the Ch from a cartilage template to a fat-filled bone in juvenile stages, which coincides with extensive vascularization. Using inducible Cre and long-lived histone-mCherry fusion proteins, driven by regulatory regions of the chondrocyte genes sox10 and col2a1a, we reveal contribution of chondrocytes to osteoblasts, adipocytes, and mesenchymal cells within the adult Ch. In mouse, LepR expression marks bone marrow cells that contribute to osteoblasts and adipocytes primarily after birth (Zhou et al., 2014b). In zebrafish, we find that growth plate chondrocytes express lepr and re-enter the cell cycle during the late hypertrophic phase, raising the possibility that Lepr +skeletal stem cells may derive from growth plate chondrocytes. Further, we find that delayed remodeling of the hypertrophic cartilage zone in zebrafish mmp9 mutants correlates with a paucity of marrow adipocytes. Unlike in mouse where Mmp9 functions in hematopoietic cells for timely growth plate remodeling (Vu et al., 1998), we find that Mmp9 is sufficient in neural crest-derived chondrocytes of zebrafish for growth plate remodeling. Our studies reveal that growth plate chondrocytes generate osteocytes and adipocytes in zebrafish bones, potentially by transitioning through a proliferative intermediate.
Results
Remodeling of the Ch bone in juvenile zebrafish
In order to characterize the progressive remodeling of an endochondral bone in zebrafish, we performed pentachrome staining on sections of the Ch bone from juvenile through adult stages (Figure 1). The Ch bone is shaped like a flattened barbell, and here we sectioned it to reveal the thin plane of the bone (see Figure 1—figure supplement 1A) for a view along the thicker perpendicular plane). Unlike the unidirectional growth plates in the mouse limb, the two growth plates of Ch are bidirectional with a central zone of compact, proliferative chondrocytes flanked by hypertrophic chondrocytes on either side (Paul et al., 2016). Unlike in many other fish species, the Ch bone, as with other bones in zebrafish, also contains embedded osteocytes (Witten and Huysseune, 2009). At 11 mm standard length (SL) (approx. 4.5 weeks post-fertilization (wpf)), the Ch contains chondrocytes throughout its length with the exception of a small marrow space at the anterior tip. The Ch is surrounded by a thin layer of cortical bone that has been shown to derive from osteoblasts located on the outside of the cartilage template (i.e. periosteum) (Paul et al., 2016). By 12 mm SL (approx. five wpf), both tips of the Ch contain marrow spaces, and on the central sides of the growth plates we begin to observe small fissures in the cortical bone and disruption of the hypertrophic zone. By 13 mm SL (approx. 5.5 wpf), breaks in the cortical bone become more prominent and are accompanied by further degradation of the cartilage matrix. At later stages (16 and 19 mm SL) (approx. 7 and 9 wpf), cortical bone regains integrity and increases in thickness, and marrow adipocytes containing LipidTOX +lipid vesicles are seen throughout Ch (Figure 1—figure supplement 1B). By adulthood (one year of age), the marrow cavity is filled with large fat cells and the growth plates appear largely mineralized. While we focus on the Ch for this study, a number of other cartilage-derived bones in the face and fins have been reported to have a similar structure in zebrafish, including growth plates and prominent marrow fat (Weigele and Franz-Odendaal, 2016).
Vascularization of the Ch bone in juvenile zebrafish
Given the transient breakdown of cortical bone, we examined whether this coincides with vascularization of Ch. To do so, we performed confocal imaging of the dissected Ch from fish carrying both a chondrocyte-specific col2a1a:mCherry-NTR transgene and fli1a:GFP (Figure 2A). We used fli1a:GFP to label endothelial cells of the vasculature, but we also noticed that presumptive resting chondrocytes in the middle of the growth plate express this transgene. At 10 and 11 mm SL, vessels expressing fli1a:GFP are found largely on the outside of the Ch bone. By 13, 16, and 20 mm SL, we observe increasing numbers of capillaries within the Ch, coinciding with the replacement of col2a1a:mCherry-NTR +chondrocytes with adipocytes. High-magnification confocal sections at 18 mm SL clearly show fli1a:GFP+ vessels in intimate association with adipocytes and within the Calcein Blue+ bone collar (Figure 1—figure supplement 1C). Given the more complicated expression pattern of fli1a:GFP, we also independently confirmed blood vessel identity with kdrl:GFP. At 28 mm SL (approx. 26 wpf), the Ch is heavily supplied with both kdrl:GFP+blood vessels and lyve1:DsRed+lymphatic vessels, which abut each side of the growth plate (Figure 2B). Hence, as in mammalian bones, remodeling of the Ch bone is accompanied by extensive vascular invasion of the cartilage template.
Contribution of sox10-lineage cells to osteoblasts, adipocytes, and mesenchymal cells
Given the extensive remodeling and vascularization of Ch, we next investigated the long-term fate of growth plate chondrocytes by multiple, independent methods. First, we constructed an inducible sox10:CreERT2 line and crossed it to a ubiquitous bactin2:loxP-tagBFP-stop-loxP-DsRed reporter (Blue to Red conversion: B > R). The zebrafish sox10 promoter drives expression in early cranial neural crest cells from 10 to 16 hpf, followed by a second wave of expression in all chondrocytes from two dpf onwards (Dutton et al., 2008). Here, we took advantage of this second wave of expression to label developmental chondrocytes. Upon addition of 4-hydroxytamoxifen (4-OHT) at 15 dpf, we observed extensive labeling of chondrocytes within 5 days, as well as some cells in the perichondrium surrounding Ch and other cartilages (Figure 3A). We did not observe leaky conversion in the absence of 4-OHT at either embryonic or adult stages (Figure 3—figure supplement 1A). We then converted sox10/B > R fish by 4-OHT treatment at 14 dpf and raised these to adulthood (27 mm SL) for analysis, with inclusion of an ocn:GFP transgene allowing us to label osteoblasts (Figure 3B). Confocal maximum intensity projections through Ch revealed extensive labeling of growth plate chondrocytes. We also observed numerous cells throughout the Ch, including in and around the marrow cavity. In sections through the middle of Ch, we observed labeling of a number of large diameter cells of adipocyte morphology (Figure 3C). In superficial sections through the cortical bone, we also observed labeled cells that were positive for the osteoblast marker ocn:GFP, as well as some labeled cells negative for ocn:GFP that may represent bone progenitors or other cell types (Figure 3D). As a comparison, we used a constitutive Cre driven by a human SOX10 enhancer that drives expression throughout the neural crest lineage (note that this human enhancer lacks the second wave of chondrocyte expression seen with the zebrafish regulatory region used for the sox10:CreERT2 line) (Kague et al., 2012). When crossed to the B > R line, this neural crest-specific SOX10:Cre line drives broader conversion in the five dpf head than the later conversion of sox10:CreERT2 (Figure 3—figure supplement 2A). Analysis of the adult Ch in SOX10:Cre fish shows labeling of all growth plate cartilage, as well as most if not all adipocytes and numerous smaller cells throughout the marrow and cortical bone surface (Figure 3—figure supplement 2B,C). This is consistent with previous studies showing that the Ch bone is neural crest-derived (Schilling and Kimmel, 1994). However, we also detect unconverted cells in the marrow, consistent with contribution of non-neural crest-derived cells such as the mesoderm-derived vasculature (Figure 2).
As sox10:CreERT2-mediated conversion at 14 dpf also labels cells outside the cartilage, which could also contribute to osteoblasts and adipocytes, we next examined animals with lower conversion efficiency to follow discrete growth plate clones. When analyzed at 30 mm SL, growth plate clones could be quite large, consistent with clonal selection as described in the zebrafish heart and skeletal muscle (Gupta and Poss, 2012; Nguyen et al., 2017). In one example with three discrete clones, we observed a clone that contributed to a narrow column of growth plate chondrocytes in the middle of Ch that was contiguous with mesenchymal cells and then adipocytes toward the central marrow cavity (Figure 3E). To more definitively follow growth plate clones, we also examined sox10:CreERT2; ubb:Zebrabow fish in which 4OH-T treatment at 14 dpf resulted in conversion of RFP (red) to various color combinations of CFP, YFP, and RFP at 23 mm SL. Analysis of uniquely colored clones in the growth plate revealed those that contained both chondrocytes and adjacent adipocytes in the marrow (Figure 3F). We also observed clones that appeared to contain chondrocytes and weakly labeled cells embedded in cortical bone, though the clonal contribution to these and other mesenchymal lineages will require further analysis. Together, our data are consistent with chondrocytes giving rise to adipocytes and mesenchymal cells, and potentially osteoblasts, after growth plate remodeling in zebrafish.
Contribution of col2a1a+ chondrocytes to osteoblasts, adipocytes, and mesenchymal cells
As sox10-CreERT2-mediated conversion was broader than just the cartilage, we also generated an inducible col2a1a-CreERT2 line to more precisely trace chondrocytes and their derivatives. In mice, Col2a1:CreERT2-mediated conversion at embryonic and early postnatal stages broadly labels not only chondrocytes but also osteochondroprogenitors, such as those in the perichondrium and periosteum (Ono et al., 2014). In zebrafish, col2a1a is similarly expressed at high levels in chondrocytes and in weaker levels in osteoblasts and perichondrium, although direct evidence for col2a1a marking osteochondroprogenitors in zebrafish is lacking (Eames et al., 2012). In order to restrict expression to chondrocytes, thus avoiding potential complications of labeling osteoblasts and putative perichondral progenitors, we utilized a chondrocyte-specific ‘R2’ enhancer of the zebrafish col2a1a gene that we had previously characterized (Dale and Topczewski, 2011; Askary et al., 2015). Treatment of col2a1a/B > R fish with a single dose of 4-OHT at five dpf resulted in extensive labeling of col2a1a-BAC:GFP+ chondrocytes at 12 dpf, but no labeling of the perichondrium, periosteum, and osteoblasts as marked by the sp7:GFP transgene (Figure 4A,B; Figure 3—figure supplement 1B). We also observed labeling of the notochord in larval fish, but no labeling of the vasculature, blood, or other tissues examined. After conversion of col2a1a/B > R chondrocytes at five dpf and examination at adulthood (30 mm SL), maximal intensity projections through Ch revealed extensive labeling of growth plate chondrocytes, as well as cells throughout the bone and marrow cavity (Figure 4C). Thus, the majority of adult growth plate chondrocytes in Ch appear to derive from embryonic chondrocytes. Moreover, optical sections revealed cytoplasmic DsRed staining in lipid-filled adipocytes, presumptive osteoblasts lining the inner surface of bone, and mesenchymal cells within the marrow cavity. In some animals displaying lower conversion efficiency, we observed apparent clones of cells containing growth plate chondrocytes, large adipocytes, and osteoblasts embedded in Calcein +mineralized matrix (Figure 4D). Analysis of individual sections at higher magnification revealed contribution of col2a1a-lineage cells to a subset of osteoblasts within both the endosteal and periosteal surfaces of bone, as well as embedded osteocytes with characteristic cellular processes (Figure 4E–G).
Long-lived Histone2A-mCherry protein reveals contribution of col2a1a+ chondrocytes to osteoblasts
A limitation of CreER-mediated lineage tracing is that it can be difficult to rule out contributions from rare converted cells outside the population of interest, for example in the perichondrium and periosteum. We therefore employed a Cre-independent approach to independently assess the fate of growth plate chondrocytes. To do so, we expressed a Histone2A-mCherry fusion protein in col2a1a+ cells. An advantage of this type of lineage approach is that the levels of Histone2A-mCherry protein, which is stably incorporated into chromatin, reflect those of endogenous col2a1a in chondrocytes provided continued cell division does not dilute out the fusion protein. This is in contrast to Cre-mediated approaches in which a strong ubiquitous promoter determines the level of a reporter protein; hence, even low levels of Cre recombinase activity outside of chondrocytes (e.g. in osteochondroprogenitors) can result in strong reporter expression. In zebrafish, col2a1a is expressed at high levels in the proliferative zone of the Ch growth plate and then downregulated in the hypertrophic zone (Paul et al., 2016). In a newly generated col2a1a:Histone2A-mCherry-T2A-GFP-CAAX transgenic line, membrane-localized GFP-CAAX, which is rapidly turned over, is seen primarily in the proliferative zone, whereas Histone2A-mCherry is seen uniformly throughout the Ch cartilage, confirming the long-lived nature of this fusion protein (Figure 5A and Figure 5—figure supplement 1A,B). At 7–8 mm SL (approx. three wpf), which is well before the start of growth plate remodeling at 11–12 mm SL, all Ch chondrocytes are Histone2A-mCherry+ and we do not detect Histone2A-mCherry+ cells associated with Calcein Blue+ bone or co-expressing the osteoblast transgene ocn:GFP (Figure 5A and Figure 5—figure supplement 1C,D). In contrast, at post-remodeling stages (12 and 18 mm SL), we observe extensive overlap of Histone2A-mCherry with ocn:GFP+ osteoblasts associated with cortical bone (Figure 5B and Figure 5—figure supplement 1E). We also observe numerous Histone2A-mCherry+ cells embedded in the endosteal and periosteal surfaces of the Ch bone (labeled by Calcein Blue), with several of these cells co-expressing the pre-osteoblast transgene RUNX2:GFP or the early osteoblast transgene sp7:GFP (Figure 5C–E). Note that ocn:GFP, RUNX2:GFP, and sp7:GFP can all be readily distinguished from membrane GFP-CAAX by their much stronger and cytoplasmic expression. In addition, we observed that sp7:GFP+ osteoblasts further from the growth plate tended to have weaker Histone2A-mCherry signal than those more closely associated with the edge of the hypertrophic zone, suggesting that hypertrophic chondrocytes and/or their osteoblast derivatives undergo cell division to dilute out the Histone2A-mCherry signal. These findings independently confirm the conclusions of our CreER lineage tracing studies that col2a1a+ chondrocytes generate osteoblasts in zebrafish.
Hypertrophic chondrocytes re-enter the cell cycle and express lepr
The conversion of hypertrophic chondrocytes to osteoblasts and adipocytes could occur in the absence of cell division (i.e. ‘transdifferentiation’) and/or through partial dedifferentiation into a proliferative progenitor. To test these possibilities, we first used incorporation of bromodeoxyuridine (BrdU) to detect proliferative cells in the Ch during remodeling stages (Figure 6A,B). At the beginning of remodeling (11 mm SL), we detected BrdU+ cells in the central zone of chondroblasts in the growth plate, as well as in the perichondrium and periosteum. In addition, we observed BrdU+ cells at the edge of the hypertrophic zone where the cartilage matrix is being actively degraded, similar to what has been reported in mouse (Park et al., 2015). We observed BrdU incorporation in similarly positioned hypertrophic chondrocytes at 15 and 19 mm SL, with BrdU+ cells becoming fewer in the perichondrium and periosteum by 19 mm.
We next tested whether hypertrophic chondrocytes that re-enter the cell cycle also express known skeletal stem cell markers. In mice, LepR expression marks a heterogeneous population of cells in endochondral bone, including a putative postnatal skeletal stem cell population (Zhou et al., 2014b). First, we examined expression of lepr and found dynamic expression in the zebrafish Ch endochondral bone from juvenile through adult stages. A comparison with the hypertrophic chondrocyte and early osteoblast marker runx2b shows higher lepr expression in proliferative versus hypertrophic chondrocytes at 8, 12, and 20 mm SL stages, with chondrocyte expression decreasing by 27 mm SL (Figure 6—figure supplement 1A). During remodeling of zebrafish Ch (15 mm SL), we also observe lepr+ cells in the marrow cavity that are derived from chondrocytes, based on labeling by col2a1a/B > R (Figure 6C), and we continue to observe lepr+ cells in the Ch marrow at later stages (20 and 27 mm SL) (Figure 6—figure supplement 1A). These results are consistent with lepr+ cells in the bone marrow deriving from growth plate chondrocytes in zebrafish, although direct evidence will be needed to determine if any of these chondrocyte-derived lepr+ marrow cells behave as skeletal stem cells in zebrafish.
Requirement of mmp9 for growth plate remodeling and marrow adipocyte formation
In mice, Mmp9 has been reported to function in hematopoietic lineage cells for growth plate remodeling, as bone marrow transplants can rescue the delay in growth plate remodeling seen in Mmp9 mutants (Vu et al., 1998). Here, we tested whether mmp9 might have a conserved requirement for growth plate remodeling in zebrafish, including the generation of marrow adipocytes. At the beginning of remodeling (11 mm SL), we observe expression of mmp9 at the edge of the hypertrophic zone, with this restricted expression in late-stage hypertrophic chondrocytes continuing through 17 mm SL stages (Figure 7A). Higher magnification views show that mmp9 expression is prominent in hypertrophic chondrocytes that appear to be exiting their lacunae. Next, we used CRISPR/Cas9 mutagenesis to create an early frame-shift mutation in the mmp9 gene that is predicted to abolish most if not all protein function (Figure 7B). mmp9 homozygous mutants are adult viable and do not display obvious larval craniofacial defects. Whereas trichrome staining revealed no significant differences in the mutant Ch growth plates at 17 mm SL, by 21 mm we observed that the hypertrophic zone was significantly larger, compared to the proliferative zone, in mmp9 mutants versus controls, indicating a delay in growth plate remodeling (Figure 7C,D). The defect in growth plate remodeling was still evident at 27 mm SL, with mutants displaying a wider Ch growth plate (Figure 7E,G). Strikingly, mmp9 mutants also had fewer adipocytes in the central marrow cavity compared to stage-matched controls (Figure 7E,G), consistent with adipocytes deriving from hypertrophic chondrocytes that are released from the cartilage matrix by Mmp9 activity.
-
Figure 7—source data 1
- https://doi.org/10.7554/eLife.42736.015
Given mmp9 expression in late-stage hypertrophic chondrocytes, we next tested whether Mmp9 might function in chondrocytes as opposed to hematopoietic lineage cells for growth plate remodeling and marrow adipocyte generation in zebrafish. As the Ch bone is generated from neural crest-derived cells, we used transplantation of ubiquitously labeled wild-type ectodermal cells into the neural crest precursor domain of unlabeled mmp9-/- shield-stage hosts to generate wild-type Ch bones in otherwise mutant hosts. At adult stages, we were able to recover eight mutant recipients with contribution of wild-type cells to chondrocytes of the Ch growth plate. In these animals, we observed a rescue of growth plate width, with a trend toward better rescue in wild-type versus mutant regions of the growth plate (p = 0.06), as well as a trend toward rescue of adipocyte number (p = 0.06) (Figure 7E–G). As a control, transplantation of wild-type neural crest cells into wild-type animals had no effect on Ch growth plate width and adipocyte number. These results indicate that mmp9 is required in the neural crest lineage, and potentially chondrocytes themselves, for efficient remodeling of the growth plate and the generation of marrow adipocytes from chondrocytes.
Discussion
Despite anatomical differences between zebrafish and mammalian bones, we find that growth plate remodeling is remarkably well conserved and thus likely ancestral to bony vertebrates. The zebrafish Ch undergoes a transient breakdown of cortical bone near the growth plates, which coincides with extensive vascularization and an Mmp9-dependent replacement of hypertrophic chondrocytes with fat and bone. Using multiple methods of lineage tracing, including a Cre-independent technique, we show that late-stage hypertrophic chondrocytes generate not only osteoblasts but also marrow adipocytes in zebrafish. Further support for the ability of chondrocytes to generate adipocytes is that delayed growth plate remodeling in mmp9 mutants results in a paucity of marrow adipocytes in adults. Lastly, we show that hypertrophic chondrocytes re-enter the cell cycle and contribute to lepr+ mesenchymal cells, raising the possibility that partial dedifferentiation into proliferative progenitors underlies chondrocyte fate transitions inside endochondral bones.
As zebrafish have hollow bones that lack hematopoiesis, one possibility was that most if not all of the cortical bone of adult zebrafish would be simply derived from the periosteum. However, we find significant contribution of chondrocyte-derived cells to both the endosteal and periosteal surfaces of the Ch bone, similar to what has been described in the mouse (Yang et al., 2014; Zhou et al., 2014a; Jing et al., 2015). We also reveal chondrocytes to be a significant source of marrow adipocytes in zebrafish. Although the incomplete conversion efficiency of the CreER lines, as well as the expression of sox10:CreER outside of chondrocytes, made it difficult to precisely quantify what proportion of osteoblasts and marrow adipocytes derive from chondrocytes, there are likely other sources for these cells, in particular the periosteum which houses several types of skeletal stem cells (Debnath et al., 2018). Our findings raise the question of whether marrow adipocytes also derive in part from chondrocytes in mammals. Indeed, older studies have demonstrated the ability of murine chondrocytes to differentiate into adipocytes in vitro (Heermeier et al., 1994; Hegert et al., 2002). Further, Col2a1:CreER cells converted at postnatal day three in mouse were found to give rise to adipocytes in the metaphyseal bone marrow, although a caveat is that Col2a1:CreER marks both chondrocytes and progenitors at early postnatal stages (Ono et al., 2014).
The finding that late-stage hypertrophic chondrocytes can re-enter the cell cycle and express lepr suggests that at least some of these cells may dedifferentiate into stem-like cells, which subsequently expand in number and differentiate into osteoblasts and adipocytes. Future molecular profiling of these chondrocyte-derived marrow mesenchymal cells will be needed to better characterize their relationship to previously identified skeletal stem cells. We also cannot rule out that some hypertrophic chondrocytes directly change into adipocytes and/or osteoblasts in the absence of cell division (i.e. ‘transdifferentiation’). Indeed, the detection of osteoblasts with strong col2a1a:Histone2A-mCherry signal suggests that in some cases chondrocytes can form osteoblasts with little to no cell division, as otherwise the Histone2A-mCherry signal would have been diluted out with successive cell divisions. On the other hand, osteoblasts farther from the growth plate tended to have weaker Histone2A-mCherry signal, suggesting proliferation of a progenitor intermediate, and we directly observed late-stage hypertrophic chondrocytes undergoing cell division. Whereas it has been suggested in mouse that only chondrocytes near the periosteal surface, that is ‘borderline chondrocytes’, may be capable of lineage plasticity (Bianco et al., 1998; Maes et al., 2010), we detected lineage clones containing growth plate chondrocytes, mesenchymal cells, and adipocytes in the central part of the growth plate (Figure 3E,F), arguing against such a model.
A notable feature of LepR-lineage cells in mice is that they contribute to osteoblasts and adipocytes primarily in postnatal phases, that is when growth plate remodeling is already underway (Yang et al., 2014). It would therefore be interesting to test whether LepR-expressing skeletal stem cells have a similar origin from hypertrophic chondrocytes in mammals. One caveat is that we detect endogenous lepr expression in both marrow cells and growth plate chondrocytes in zebrafish. However, the more specific labeling of marrow cells by the LepR-Cre in mouse likely reflects the Cre insertion being in the long LepR isoform (containing exon 18b) that displays more restricted expression than the short isoform (Zhou et al., 2014b). Indeed, LepR mRNA and protein has also been reported in chondrocytes of mouse (Hoggard et al., 1997), rat, and human (Morroni et al., 2004), a finding we confirmed in the postnatal mouse femur, including the same higher expression in immature versus hypertrophic chondrocytes that we observe in zebrafish (Figure 6—figure supplement 1B). Without a comparable long-isoform lepr-Cre line in zebrafish, we cannot therefore conclude whether zebrafish lepr+ marrow cells are comparable to those described in mouse. The future generation of Cre lines to specifically mark lepr+ marrow cells in zebrafish will be needed to determine whether these chondrocyte-derived marrow cells also act as stem cells for osteoblasts and adipocytes in post-embryonic fish.
A similar delay in the remodeling of the hypertrophic zone in mouse Mmp9 and zebrafish mmp9 mutants reveals genetic conservation of growth plate remodeling from fish to mammals. In contrast to mice, we find that Mmp9 in zebrafish appears to function primarily in chondrocytes for growth plate remodeling. We also observed rescue of marrow adipocyte number by wild-type chondrocytes, though this had moderate statistical significance (p = 0.06), potentially owing to low sample size (n = 8), the mosaic contribution of wild-type cells to the growth plate, and/or roles of cells outside the growth plate to marrow adipocyte generation. In mice, loss of Mmp13 in chondrocytes does result in a growth plate remodeling delay, with global Mmp13 deletion enhancing the remodeling defects of Mmp9 mutants (Inada et al., 2004; Stickens et al., 2004). It may be that MMPs are derived from both chondrocytes and invading hematopoietic cells (e.g. osteoclasts), with the relative importance of each cell source varying between zebrafish and mammals. Compensation by mmp13 might also explain why we detected growth remodeling defects in mmp9 zebrafish mutants at 21 and 27–31 mm SL stages, but not at 17 mm SL. Mmp9 and Mmp13 have known roles in degrading components of the cartilage extracellular matrix (Page-McCaw et al., 2007), consistent with our observed loss of collagen-rich matrix in hypertrophic chondrocytes at the edges of the zebrafish growth plate. Secreted Mmp9 may therefore function simply to degrade cartilage matrix and facilitate release of dedifferentiating chondrocytes into the marrow. Intriguingly, Mmp9 has recently been shown to have an additional, non-canonical function in the nucleus for histone H3 tail cleavage. Whereas this has only been demonstrated so far for osteoclasts (Kim et al., 2016), it remains possible that Mmp9 could have a similar non-canonical function in altering the chromatin structure of hypertrophic chondrocytes to allow them to acquire new potential. In conclusion, our data support conservation of mammalian-like growth plate remodeling in zebrafish, which provides new opportunities for better understanding the molecular and cellular mechanisms by which hypertrophic chondrocytes transform into osteoblasts, marrow adipocytes, and potentially adult skeletal stem cells within endochondral bones.
Materials and methods
Zebrafish transgenic lines and mmp9 mutants
Request a detailed protocolAll procedures were approved by the University of Southern California Institutional Animal Care and Use Committee. Published Danio rerio lines include Tg(Has.RUNX2:EGFP)zf259 and Tg(Mmu.Sox10-Mmu.Fos:Cre)zf384 (Kague et al., 2012), Tg(sp7:EGFP)b1212 (DeLaurier et al., 2010), Tg(Ola.Osteocalcin.1:EGFP)hu4008 (Knopf et al., 2011), Tg(col2a1aBAC:GFP)el483 and Tg(col2a1aBAC:mCherry-NTR)el559 (Askary et al., 2015), Tg(bactin2:loxP-BFP-loxP-DsRed)sd27 (Kobayashi et al., 2014), Tg(fli1a:eGFP)y1 (Lawson and Weinstein, 2002), Tg(kdrl:eGFP)s843 (Jin et al., 2005), Tg(−5.2lyve1b:DsRed)nz101 (Okuda et al., 2012), and Zebrabow - Tg(ubb:LOX2272-LOXP-RFP-LOX2272-CFP-LOXP-YFP)a131 (Pan et al., 2013). The sox10:CreERT2 transgene was generated with Gateway Cloning (Invitrogen) and the Tol2 kit (Kwan et al., 2007) by combining p5E-sox10 (Das and Crump, 2012), pME-CreERT2, p3E-pA, and pDestTol2pA2. The col2a1a-R2-E1b:CreERT2 and col2a1a-R2-E1b:H2A.F/Z-mCherry-P2A-GFPCAAX transgenes utilize a zebrafish col2a1a R2 enhancer element with a minimal E1B promoter sequence (Dale and Topczewski, 2011). For col2a1a-R2-E1b:CreERT2, CreERT2 was amplified using pENTR/D-CreERT2 as template and primers 5’- TTCTTGTACAAAGTGGCCACCGGCCACCATGTCCAATTTACTGACCGTACAC-3’ and 5’-TAGAGGCTCGAGAGGCCTTGTCAAGCTGTGGCAGGGAAACCCTC-3’. The amplified PCR product was combined with an NcoI/EcoRI fragment of pDestTol2-col2a1aR2-E1B-eGFPpA (Askary et al., 2015) using Gibson Assembly (New England Biolabs). For col2a1a-R2-E1b:H2A.F/Z-mCherry-P2A-GFPCAAX, H2A.F/Z-mCherry was amplified using template pME-H2A.F/Z-mCherry and primers 5’-TTTCTTGTACAAAGTGGCCAAAGCTTGGATCCCGGCCACCATGGCAGGTGGAAAAGCAGG-3’ and 5’-AAGTTGGTTGCTCCCGACCCCTTGTACAGCTCGTCCATGCCGCCGGTG-3’. P2A-GFPCAAX was synthesized as a gBlock (IDT). H2A.F/Z-mCherry and P2A-GFPCAAX were combined with an NcoI/EcoRI fragment of pDestTol2-col2a1aR2-E1B-eGFPpA using Gibson Assembly. All transgenes were injected at 30 ng/μL along with 50 ng/μL Tol2 mRNA into one-cell-stage embryos, with these animals raised and outcrossed to identify stable germline founders. CreERT2 transgenes were injected into Tg(bactin2:loxP-BFP-loxP-DsRed)sd27 embryos, followed by overnight treatment with 10 μM (Z)−4-Hydroxytamoxifen (4-OHT) (Sigma-Aldrich H7904) at two dpf and screening for DsRed+ chondrocytes at six dpf using a Leica fluorescent stereomicroscope. Embryos with fluorescent cartilages were raised to adulthood and outcrossed to identify founders. We used Tg(sox10:CreERT2)el777, Tg(col2a1a-R2-E1b:CreERT2)el712, and Tg(col2a1a-R2-E1b:H2A.F/Z-mCherry-P2A-GFPCAAX)el695 lines for our experiments. Two additional alleles of Tg(col2a1a-R2-E1b:CreERT2) and one additional allele of Tg(col2a1a-R2-E1b:H2A.F/Z-mCherry-P2A-GFPCAAX) gave similar cartilage-specific expression. To generate mmp9el734 we used CRISPR/Cas9 mutagenesis to target the second exon. gRNA targeting the sequence 5’-TTGATGCCATGAAGCAGCCC-3’ was injected at 25 ng/μl with 50 ng/μl Cas9 RNA into one-cell-stage embryos as described (Hwang et al., 2013). The mmp9el734 allele is an 8 bp deletion that results in the incorporation of 13 additional amino acids after amino acid 98 (P), followed by a premature stop codon that is predicted to completely abolish the catalytic metalloproteinase domain.
Histology and LipidTOX staining
Request a detailed protocolLive bone staining of dissected ceratohyal bones was performed by treating with 50 µg/ml Alizarin Red (Sigma Aldrich, cat. no. A5533), Calcein Green (Thermofisher Scientific, cat. no. C481), or Calcein Blue, AM (Thermofisher Scientific, cat. no. C1429) for 5 min and repeatedly rinsing in embryo medium as described (Paul et al., 2016). We performed adipocyte labeling of dissected ceratohyal bones by incubating in a 1:200 solution of HCS LipidTOX Deep Red (Life Technologies, cat. no. H34477) for 15 min and rinsing in embryo medium as described (Minchin and Rawls, 2017). Paraffin embedding and histology were performed as described (Paul et al., 2016). We cut blocks into 5 μm sections on a Shandon Finesse Me+ microtome (cat. no. 77500102) and collected sections on Apex Superior Adhesive slides (Leica Microsystems, cat. no. 3800080). Pentachrome and Trichrome staining were performed according to manufacturer's instructions (Movat-Russell modified pentachrome stain kit, Newcomer Supply cat. no. 9150A; Gomori One-Step, Aniline Blue, trichrome stain kit, Newcomer Supply cat. no. 9176A).
Immunohistochemistry and in situ hybridization
Request a detailed protocolFor cell proliferations assays, fish were immersed in system water containing 4.5 mg/ml BrdU (Sigma Aldrich, cat. no. B5002) for 1 hr, followed by two washes in system water, fixation in 4% paraformaldehyde, and paraffin embedding. Immunohistochemistry was performed as described except that we blocked with 2% normal goat serum (Jackson ImmunoResearch, cat. no. 005-000-121). After cutting thin sections, we performed antigen retrieval by treating slides with citrate buffer (pH 6.0) in a steamer set (IHC World, cat. no. IW-1102) for 35 min. Primary antibodies include rat anti-BrdU (1:100, Bio-Rad, cat. no. MCA2060GA) and rabbit anti-mCherry (1:250, Novus Biologicals, cat. no. NBP2-25157). We used AlexaFluor secondary antibodies and Hoechst to visualize nuclei. Colorimetric in situ hybridization was performed as described (Paul et al., 2016). The mmp9 riboprobe was generated by PCR amplification of zebrafish genomic DNA with primers 5’-GTCTCCAATACTAAAGCTCTGAAGAAG-3’ and 5-‘TAGGATGTCGAAGGTCTATAGAGAATG-3’ and cloning into pCR-BluntII-TOPO (Life Technologies). RNA probe was synthesized with T7 polymerase (Roche) after linearizing the plasmid with BamHI restriction. RNAscope probes for leptin receptor (lepr) and runx2b were synthesized by Advanced Cell Diagnostics in Channel 1 and Channel 2, respectively, and for mouse Lepr in channel 1. Paraformaldehyde-fixed paraffin-embedded sections were deparaffinized, and the RNAscope fluorescent multiplex v2 assay combined with immunofluorescence was performed according to manufacturer’s protocols and with the ACD HybEZ Hybridization oven.
Neural crest transplantations
Request a detailed protocolAt six hpf, donor ectoderm from the animal cap of bactin2:loxP-tagBFP-loxP-DsRed embryos was transplanted into the neural crest precursor domain of mmp9-/- hosts as described (Crump et al., 2004). Embryos displaying unilateral tagBFP fluorescence in the face at three dpf were raised in the nursery and then genotyped at 14 dpf for the mmp9 mutant allele. Homozygous mutant and wild-type siblings were raised at similar density until they reached the indicated sizes for analysis.
Imaging
Request a detailed protocolBrightfield images of pentachrome and trichrome stains, and colorimetric in situ hybridizations, were acquired with a Zeiss AxioImager.A1 microscope and a Zeiss slide scanner AxioScan.Z1. Focus stacking of multiple images was done in Adobe Photoshop CS5. Fluorescence images were captured on a Zeiss LSM800 confocal microscope, with representative sections or maximum intensity projections shown as specified. Tiling was performed using ZEN software or manually stitched together using Fiji. Brightness and contrast were adjusted in Adobe Photoshop CS5 with similar settings for experimental and control samples.
Quantification and statistical analyses
Request a detailed protocolWe stage-matched control and experimental zebrafish by measuring standard body length (SL) from the tip of the snout to the edge of the hypuralia. Adipocyte counts and area/width measurements of growth plates were calculated using Fiji. The proliferative zone of the growth plate was defined as the central region of compact chondrocytes. All measurements were performed blinded to genotype. Statistical significance was determined by a student’s t-test for pair-wise comparisons or Tukey-Kramer HSD tests for multiple comparisons, using GraphPad’s Prism software.
Data availability
All data generated or analysed during this study are included in the manuscript and supporting files.
References
-
Evolution of the osteoblast: skeletogenesis in gar and zebrafishBMC Evolutionary Biology 12:27.https://doi.org/10.1186/1471-2148-12-27
-
Adipose differentiation of cartilage in vitroDifferentiation 56:45–53.https://doi.org/10.1046/j.1432-0436.1994.56120045.x
-
Differentiation plasticity of chondrocytes derived from mouse embryonic stem cellsJournal of Cell Science 115:4617–4628.https://doi.org/10.1242/jcs.00171
-
Efficient genome editing in zebrafish using a CRISPR-Cas systemNature Biotechnology 31:227–229.https://doi.org/10.1038/nbt.2501
-
Chondrocytes directly transform into bone cells in mandibular condyle growthJournal of Dental Research 94:1668–1675.https://doi.org/10.1177/0022034515598135
-
The Tol2kit: a multisite gateway-based construction kit for Tol2 transposon transgenesis constructsDevelopmental Dynamics 236:3088–3099.https://doi.org/10.1002/dvdy.21343
-
In vivo imaging of embryonic vascular development using transgenic zebrafishDevelopmental Biology 248:307–318.https://doi.org/10.1006/dbio.2002.0711
-
A classification system for zebrafish adipose tissuesDisease Models & Mechanisms 10:797–809.https://doi.org/10.1242/dmm.025759
-
A subset of chondrogenic cells provides early mesenchymal progenitors in growing bonesNature Cell Biology 16:1157–1167.https://doi.org/10.1038/ncb3067
-
Matrix metalloproteinases and the regulation of tissue remodellingNature Reviews Molecular Cell Biology 8:221–233.https://doi.org/10.1038/nrm2125
-
Segment and cell type lineage restrictions during pharyngeal arch development in the zebrafish embryoDevelopment 120:483–494.
-
Osteogenesis by chondrocytes from growth cartilage of rat ribCalcified Tissue Research 19:179–187.https://doi.org/10.1007/BF02564002
-
The role of three genetically distinct collagen types in endochondral ossification and calcification of cartilageThe Journal of Bone and Joint Surgery. British Volume 59-B:458–464.https://doi.org/10.1302/0301-620X.59B4.72756
Article and author information
Author details
Funding
National Institute of Dental and Craniofacial Research (R35 DE027550)
- J Gage Crump
Burroughs Wellcome Fund (Postdoctoral Fellowship)
- D'Juan T Farmer
National Institute of Dental and Craniofacial Research (F31 025549)
- Dion Giovannone
National Institute of Dental and Craniofacial Research (K99 DE027218)
- Joanna Smeeton
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Megan Matsutani and Jennifer DeKoeyer Crump for fish care, Francesca Mariani for insightful conversations, David Traver for the bactin2:loxP-tagBFP-stop-loxP-DsRed line, and Ken Poss, Leonard Zon, and Rodney Dale for plasmids.
Ethics
Animal experimentation: This study was performed in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health. All of the animals were handled according to approved institutional animal care and use committee (IACUC) protocols (#20771) of the University of Southern California.
Copyright
© 2019, Giovannone 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
-
- 4,838
- views
-
- 642
- downloads
-
- 57
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
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)
Further reading
-
- Cell Biology
- Developmental Biology
How cells regulate the size of their organelles remains a fundamental question in cell biology. Cilia, with their simple structure and surface localization, provide an ideal model for investigating organelle size control. However, most studies on cilia length regulation are primarily performed on several single-celled organisms. In contrast, the mechanism of length regulation in cilia across diverse cell types within multicellular organisms remains a mystery. Similar to humans, zebrafish contain diverse types of cilia with variable lengths. Taking advantage of the transparency of zebrafish embryos, we conducted a comprehensive investigation into intraflagellar transport (IFT), an essential process for ciliogenesis. By generating a transgenic line carrying Ift88-GFP transgene, we observed IFT in multiple types of cilia with varying lengths. Remarkably, cilia exhibited variable IFT speeds in different cell types, with longer cilia exhibiting faster IFT speeds. This increased IFT speed in longer cilia is likely not due to changes in common factors that regulate IFT, such as motor selection, BBSome proteins, or tubulin modification. Interestingly, longer cilia in the ear cristae tend to form larger IFT compared to shorter spinal cord cilia. Reducing the size of IFT particles by knocking down Ift88 slowed IFT speed and resulted in the formation of shorter cilia. Our study proposes an intriguing model of cilia length regulation via controlling IFT speed through the modulation of the size of the IFT complex. This discovery may provide further insights into our understanding of how organelle size is regulated in higher vertebrates.
-
- Developmental Biology
- Genetics and Genomics
Repurposing of pleiotropic factors during execution of diverse cellular processes has emerged as a regulatory paradigm. Embryonic development in metazoans is controlled by maternal factors deposited in the egg during oogenesis. Here, we explore maternal role(s) of Caspar (Casp), the Drosophila orthologue of human Fas-associated factor-1 (FAF1) originally implicated in host-defense as a negative regulator of NF-κB signaling. Maternal loss of either Casp or it’s protein partner, transitional endoplasmic reticulum 94 (TER94) leads to partial embryonic lethality correlated with aberrant centrosome behavior, cytoskeletal abnormalities, and defective gastrulation. Although ubiquitously distributed, both proteins are enriched in the primordial germ cells (PGCs), and in keeping with the centrosome problems, mutant embryos display a significant reduction in the PGC count. Moreover, the total number of pole buds is directly proportional to the level of Casp. Consistently, it’s ‘loss’ and ‘gain’ results in respective reduction and increase in the Oskar protein levels, the master determinant of PGC fate. To elucidate this regulatory loop, we analyzed several known components of mid-blastula transition and identify the translational repressor Smaug, a zygotic regulator of germ cell specification, as a potential critical target. We present a detailed structure-function analysis of Casp aimed at understanding its novel involvement during PGC development.