On the evolutionary relationship between chondrocytes and osteoblasts

Xenopus tropicalis osteoblast-specific open chromatin regions reveal promoters and enhancers involved in human skeletal phenotypes and shed light on early vertebrate evolution

While understanding the genetic underpinnings of osteogenesis has far-reaching implications for skeletal diseases and evolution, a comprehensive characterization of the osteoblastic regulatory landscape in non-mammalian vertebrates is still lacking. Here, we compared the ATAC-Seq profile of Xenopus tropicalis (Xt) osteoblasts to a variety of non mineralizing control tissues, and identified osteoblast-specific nucleosome free regions (NFRs) at 527 promoters and 6747 distal regions. Sequence analyses, Gene Ontology, RNA-Seq and ChIP-Seq against four key histone marks confirmed that the distal regions correspond to bona fide osteogenic transcriptional enhancers exhibiting a shared regulatory logic with mammals. We report 425 regulatory regions conserved with human and globally associated to skeletogenic genes. Of these, 35 regions have been shown to impact human skeletal phenotypes by GWAS, including one trps1 enhancer and the runx2 promoter, two genes which are respectively involved in trichorhinophalangeal syndrome type I and cleidocranial dysplasia. Intriguingly, 60 osteoblastic NFRs also align to the genome of the elephant shark, a species lacking osteoblasts and bone tissue. To tackle this paradox, we chose to focus on dlx5 because its conserved promoter, known to integrate regulatory inputs during mammalian osteogenesis, harbours an osteoblast-specific NFR in both frog and human. Hence, we show that dlx5 is expressed in Xt and elephant shark odontoblasts, supporting a common cellular and genetic origin of bone and dentine. Taken together, our work (i) unravels the Xt osteogenic regulatory landscape, (ii) illustrates how cross-species comparisons harvest data relevant to human biology and (iii) reveals that a set of genes including bnc2, dlx5, ebf3, mir199a, nfia, runx2 and zfhx4 drove the development of a primitive form of mineralized skeletal tissue deep in the vertebrate lineage.

Oral citrate supplementation mitigates age-associated pathological intervertebral disc calcification in LG/J mice

Despite the high prevalence of age-dependent intervertebral disc calcification, there is a glaring lack of treatment options for this debilitating pathology. Here, we investigate the efficacy of long-term oral K3Citrate supplementation in ameliorating disc calcification in LG/J mice, a model of spontaneous age-associated disc calcification. K3Citrate successfully reduced the incidence of disc calcification in LG/J mice without deleterious effects on vertebral bone structure, plasma chemistry, and locomotion. Notably, a positive effect on grip strength was evident in treated mice. Spectroscopic investigation of the persisting calcified nodules indicated K3Citrate did not alter the mineral composition and revealed that reactivation of an endochondral differentiation program in endplates may drive LG/J disc calcification. Importantly, K3Citrate reduced calcification incidence without altering the pathological endplate chondrocyte hypertrophy, suggesting mitigation of disc calcification primarily occurred through Ca2+ chelation, a conclusion supported by chondrogenic differentiation and Seahorse metabolic assays. Overall, this study underscores the therapeutic potential of K3Citrate as a systemic intervention strategy for disc calcification.

Estradiol-17β and bisphenol A affect growth and mineralization in early life stages of seabass

Natural and synthetic estrogens are contaminants present in aquatic ecosystems. They can have significant consequences on the estrogen-sensitive functions of organisms, including skeletal development and growth of vertebrate larvae. Synthetic polyphenols represent a group of environmental xenoestrogens capable of binding the receptors for the natural hormone estradiol-17β (E2). To better understand how (xeno-)estrogens can affect the skeleton in fish species with high ecological and commercial interest, 16 days post-hatch larvae of the seabass were experimentally exposed for 7 days to E2 and Bisphenol A (BPA), both used at the regulatory concentration of surface water quality (E2: 0.4 ng.L-1, BPA: 1.6 μg.L-1) or at a concentration 100 times higher. Skeletal mineralization levels were evaluated using Alizarin red staining, and expression of several genes playing key roles in growth, skeletogenesis and estrogen signaling pathways was assessed by qPCR. Our results show that E2 exerts an overall negative effect on skeletal mineralization at the environmental concentration of 0.4 ng.L-1, correlated with an increase in the expression of genes associated only with osteoblast bone cells. Both BPA exposures inhibited mineralization with less severe effects and modified bone homeostasis by regulating the expression of gene encoding osteoblasts and osteoclasts markers. Our results demonstrate that environmental E2 exposure inhibits larval growth and has an additional inhibitory effect on skeleton mineralization while both BPA exposures have marginal inhibitory effect on skeletal mineralization. All exposures have significant effects on transcriptional levels of genes involved in the skeletal development of seabass larvae.

Common features of cartilage maturation are not conserved in an amphibian model

BACKGROUND

Mouse, chick, and zebrafish undergo a highly conserved program of cartilage maturation during endochondral ossification (bone formation via a cartilage template). Standard histological and molecular features of cartilage maturation are chondrocyte hypertrophy, downregulation of the chondrogenic markers Sox9 and Col2a1, and upregulation of Col10a1. We tested whether cartilage maturation is conserved in an amphibian, the western clawed frog Xenopus tropicalis, using in situ hybridization for standard markers and a novel laser-capture microdissection RNAseq data set. We also functionally tested whether thyroid hormone drives cartilage maturation in X tropicalis, as it does in other vertebrates.

RESULTS

The developing frog humerus mostly followed the standard progression of cartilage maturation. Chondrocytes gradually became hypertrophic as col2a1 and sox9 were eventually down-regulated, but col10a1 was not up-regulated. However, the expression levels of several genes associated with the early formation of cartilage, such as acan, sox5, and col9a2, remained highly expressed even as humeral chondrocytes matured. Greater deviances were observed in head cartilages, including the ceratohyal, which underwent hypertrophy within hours of becoming cartilaginous, maintained relatively high levels of col2a1 and sox9, and lacked col10a1 expression. Interestingly, treating frog larvae with thyroid hormone antagonists did not specifically reduce head cartilage hypertrophy, resulting rather in a global developmental delay.

CONCLUSION

These data reveal that basic cartilage maturation features in the head, and to a lesser extent in the limb, are not conserved in X tropicalis. Future work revealing how frogs deviate from the standard cartilage maturation program might shed light on both evolutionary and health studies.

Multi-modulation of immune-inflammatory response using bioactive molecule-integrated PLGA composite for spinal fusion

Despite current developments in bone substitute technology for spinal fusion, there is a lack of adequate materials for bone regeneration in clinical applications. Recombinant human bone morphogenetic protein-2 (rhBMP-2) is commercially available, but a severe inflammatory response is a known side effect. Bone graft substitutes that enhance osteogenesis without adverse effects are needed. We developed a bioactive molecule-laden PLGA composite with multi-modulation for bone fusion. This bioresorbable composite scaffold was considered for bone tissue engineering. Among the main components, magnesium hydroxide (MH) aids in reduction of acute inflammation affecting disruption of new bone formation. Decellularized bone extracellular matrix (bECM) and demineralized bone matrix (DBM) composites were used for osteoconductive and osteoinductive activities. A bioactive molecule, polydeoxyribonucleotide (PDRN, PN), derived from trout was used for angiogenesis during bone regeneration. A nano-emulsion method that included Span 80 was used to fabricate bioactive PLGA-MH-bECM/DBM-PDRN (PME2/PN) composite to obtain a highly effective and safe scaffold. The synergistic effect provided by PME2/PN improved not only osteogenic and angiogenic gene expression for bone fusion but also improved immunosuppression and polarization of macrophages that were important for bone tissue repair, using a rat model of posterolateral spinal fusion (PLF). It thus had sufficient biocompatibility and bioactivity for spinal fusion.

Limb Mesoderm and Head Ectomesenchyme Both Express a Core Transcriptional Program During Chondrocyte Differentiation

To explain how cartilage appeared in different parts of the vertebrate body at discrete times during evolution, we hypothesize that different embryonic populations co-opted expression of a core gene regulatory network (GRN) driving chondrocyte differentiation. To test this hypothesis, laser-capture microdissection coupled with RNA-seq was used to reveal chondrocyte transcriptomes in the developing chick humerus and ceratobranchial, which are mesoderm- and neural crest-derived, respectively. During endochondral ossification, two general types of chondrocytes differentiate. Immature chondrocytes (IMM) represent the early stages of cartilage differentiation, while mature chondrocytes (MAT) undergo additional stages of differentiation, including hypertrophy and stimulating matrix mineralization and degradation. Venn diagram analyses generally revealed a high degree of conservation between chondrocyte transcriptomes of the limb and head, including SOX9, COL2A1, and ACAN expression. Typical maturation genes, such as COL10A1, IBSP, and SPP1, were upregulated in MAT compared to IMM in both limb and head chondrocytes. Gene co-expression network (GCN) analyses of limb and head chondrocyte transcriptomes estimated the core GRN governing cartilage differentiation. Two discrete portions of the GCN contained genes that were differentially expressed in limb or head chondrocytes, but these genes were enriched for biological processes related to limb/forelimb morphogenesis or neural crest-dependent processes, respectively, perhaps simply reflecting the embryonic origin of the cells. A core GRN driving cartilage differentiation in limb and head was revealed that included typical chondrocyte differentiation and maturation markers, as well as putative novel "chondrocyte" genes. Conservation of a core transcriptional program during chondrocyte differentiation in both the limb and head suggest that the same core GRN was co-opted when cartilage appeared in different regions of the skeleton during vertebrate evolution.

The impact of Drew Noden's work on our understanding of craniofacial musculoskeletal integration

The classical anatomist Drew Noden spearheaded craniofacial research, laying the foundation for our modern molecular understanding of development, evolution and disorders of the craniofacial skeleton. His work revealed the origin of cephalic musculature and the role of cranial neural crest in early formation and patterning of the head musculoskeletal structures. Much of modern cranial tendon research advances a foundation of knowledge that Noden built using classical quail-chick transplantation experiments. This elegant avian chimeric system involves grafting of donor quail cells into host chick embryos to identify the cell types they can form and their interactions with the surrounding tissues. In this review, we will give a brief background of vertebrate head formation and the impact of cranial neural crest on the patterning, development and evolution of the head musculoskeletal attachments. Using the zebrafish as a model system, we will discuss examples of modifications of craniofacial structures in evolution with a special focus on the role of tendon and ligaments. Lastly, we will discuss pathologies in craniofacial tendons and the importance of understanding the molecular and cellular dynamics during craniofacial tendon development in human disease. This article is protected by copyright. All rights reserved.

The Evolution of Biomineralization through the Co-Option of Organic Scaffold Forming Networks

Biomineralization is the process in which organisms use minerals to generate hard structures like teeth, skeletons and shells. Biomineralization is proposed to have evolved independently in different phyla through the co-option of pre-existing developmental programs. Comparing the gene regulatory networks (GRNs) that drive biomineralization in different species could illuminate the molecular evolution of biomineralization. Skeletogenesis in the sea urchin embryo was extensively studied and the underlying GRN shows high conservation within echinoderms, larval and adult skeletogenesis. The organic scaffold in which the calcite skeletal elements form in echinoderms is a tubular compartment generated by the syncytial skeletogenic cells. This is strictly different than the organic cartilaginous scaffold that vertebrates mineralize with hydroxyapatite to make their bones. Here I compare the GRNs that drive biomineralization and tubulogenesis in echinoderms and in vertebrates. The GRN that drives skeletogenesis in the sea urchin embryo shows little similarity to the GRN that drives bone formation and high resemblance to the GRN that drives vertebrates’ vascular tubulogenesis. On the other hand, vertebrates’ bone-GRNs show high similarity to the GRNs that operate in the cells that generate the cartilage-like tissues of basal chordate and invertebrates that do not produce mineralized tissue. These comparisons suggest that biomineralization in deuterostomes evolved through the phylum specific co-option of GRNs that control distinct organic scaffolds to mineralization.

Magnesium Picolinate Improves Bone Formation by Regulation of RANK/RANKL/OPG and BMP-2/Runx2 Signaling Pathways in High-Fat Fed Rats

Magnesium (Mg) deficiency may affect bone metabolism by increasing osteoclasts, decreasing osteoblasts, promoting inflammation/oxidative stress, and result in subsequent bone loss. The objective of the present study was to identify the molecular mechanism underlying the bone protective effect of different forms of Mg (inorganic magnesium oxide (MgO) versus organic magnesium picolinate (MgPic) compound) in rats fed with a high-fat diet (HFD). Forty-two Wistar albino male rats were divided into six group (n = 7): (i) control, (ii) MgO, (iii) MgPic, (iv) HFD, (v) HFD + MgO, and (vi) HFD + MgPic. Bone mineral density (BMD) increased in the Mg supplemented groups, especially MgPic, as compared with the HFD group (p < 0.001). As compared with the HFD + MgO group, the HFD + MgPic group had higher bone P (p < 0.05) and Mg levels (p < 0.001). In addition, as compared to MgO, MgPic improved bone formation by increasing the levels of osteogenetic proteins (COL1A1 (p < 0.001), BMP2 (p < 0.001), Runx2 (p < 0.001), OPG (p < 0.05), and OCN (p < 0.001), IGF-1 (p < 0.001)), while prevented bone resorption by reducing the levels of RANK and RANKL (p < 0.001). In conclusion, the present data showed that the MgPic could increase osteogenic protein levels in bone more effectively than MgO, prevent bone loss, and contribute to bone formation in HFD rats.

Comparative Analyses of Gene Co-expression Networks: Implementations and Applications in the Study of Evolution

Similarities and differences in the associations of biological entities among species can provide us with a better understanding of evolutionary relationships. Often the evolution of new phenotypes results from changes to interactions in pre-existing biological networks and comparing networks across species can identify evidence of conservation or adaptation. Gene co-expression networks (GCNs), constructed from high-throughput gene expression data, can be used to understand evolution and the rise of new phenotypes. The increasing abundance of gene expression data makes GCNs a valuable tool for the study of evolution in non-model organisms. In this paper, we cover motivations for why comparing these networks across species can be valuable for the study of evolution. We also review techniques for comparing GCNs in the context of evolution, including local and global methods of graph alignment. While some protein-protein interaction (PPI) bioinformatic methods can be used to compare co-expression networks, they often disregard highly relevant properties, including the existence of continuous and negative values for edge weights. Also, the lack of comparative datasets in non-model organisms has hindered the study of evolution using PPI networks. We also discuss limitations and challenges associated with cross-species comparison using GCNs, and provide suggestions for utilizing co-expression network alignments as an indispensable tool for evolutionary studies going forward.

Nasal Septum Deviation as the Consequence of BMP-Controlled Changes to Cartilage Properties

The nasal septum cartilage is a specialized hyaline cartilage important for normal midfacial growth. Abnormal midfacial growth is associated with midfacial hypoplasia and nasal septum deviation (NSD). However, the underlying genetics and associated functional consequences of these two anomalies are poorly understood. We have previously shown that loss of Bone Morphogenetic Protein 7 (BMP7) from neural crest (BMP7 ncko ) leads to midfacial hypoplasia and subsequent septum deviation. In this study we elucidate the cellular and molecular abnormalities underlying NSD using comparative gene expression, quantitative proteomics, and immunofluorescence analysis. We show that reduced cartilage growth and septum deviation are associated with acquisition of elastic cartilage markers and share similarities with osteoarthritis (OA) of the knee. The genetic reduction of BMP2 in BMP7 ncko mice was sufficient to rescue NSD and suppress elastic cartilage markers. To our knowledge this investigation provides the first genetic example of an in vivo cartilage fate switch showing that this is controlled by the relative balance of BMP2 and BMP7. Cellular and molecular changes similar between NSD and knee OA suggest a related etiology underlying these cartilage abnormalities.

Assessing evolutionary and developmental transcriptome dynamics in homologous cell types

Background

During development, complex organ patterns emerge through the precise temporal and spatial specification of different cell types. On an evolutionary timescale, these patterns can change, resulting in morphological diversification. It is generally believed that homologous anatomical structures are built—largely—by homologous cell types. However, whether a common evolutionary origin of such cell types is always reflected in the conservation of their intrinsic transcriptional specification programs is less clear.

Results

Here, we developed a user‐friendly bioinformatics workflow to detect gene co‐expression modules and test for their conservation across developmental stages and species boundaries. Using a paradigm of morphological diversification, the tetrapod limb, and single‐cell RNA‐sequencing data from two distantly related species, chicken and mouse, we assessed the transcriptional dynamics of homologous cell types during embryonic patterning. With mouse limb data as reference, we identified 19 gene co‐expression modules with varying tissue or cell type‐restricted activities. Testing for co‐expression conservation revealed modules with high evolutionary turnover, while others seemed maintained—to different degrees, in module make‐up, density or connectivity—over developmental and evolutionary timescales.

Conclusions

We present an approach to identify evolutionary and developmental dynamics in gene co‐expression modules during patterning‐relevant stages of homologous cell type specification using single‐cell RNA‐sequencing data.

We present an approach to identify evolutionary and developmental dynamics in gene co‐expression modules during patterning‐relevant stages of homologous cell type specification using single‐cell RNA‐sequencing data.

Somite Compartments in Amphioxus and Its Implications on the Evolution of the Vertebrate Skeletal Tissues

Mineralized skeletal tissues of vertebrates are an evolutionary novelty within the chordate lineage. While the progenitor cells that contribute to vertebrate skeletal tissues are known to have two embryonic origins, the mesoderm and neural crest, the evolutionary origin of their developmental process remains unclear. Using cephalochordate amphioxus as our model, we found that cells at the lateral wall of the amphioxus somite express SPARC (a crucial gene for tissue mineralization) and various collagen genes. During development, some of these cells expand medially to surround the axial structures, including the neural tube, notochord and gut, while others expand laterally and ventrally to underlie the epidermis. Eventually these cell populations are found closely associated with the collagenous matrix around the neural tube, notochord, and dorsal aorta, and also with the dense collagen sheets underneath the epidermis. Using known genetic markers for distinct vertebrate somite compartments, we showed that the lateral wall of amphioxus somite likely corresponds to the vertebrate dermomyotome and lateral plate mesoderm. Furthermore, we demonstrated a conserved role for BMP signaling pathway in somite patterning of both amphioxus and vertebrates. These results suggest that compartmentalized somites and their contribution to primitive skeletal tissues are ancient traits that date back to the chordate common ancestor. The finding of SPARC-expressing skeletal scaffold in amphioxus further supports previous hypothesis regarding SPARC gene family expansion in the elaboration of the vertebrate mineralized skeleton.

Histological and molecular characterization of the growing nasal septum in mice

The nasal septum is a cartilaginous structure that serves as a pacemaker for the development of the midface. The septum is a hyaline cartilage which is surrounded by a perichondrium and epithelium. It remains cartilaginous anteriorly, but posteriorly it undergoes endochondral ossification to form the perpendicular plate of the ethmoid. Understanding of hyaline cartilage differentiation stems predominantly from investigations of growth plate cartilage. It is currently unclear if the morphological and molecular properties of the differentiating nasal septum align with what is known from the growth plate. In this study, we describe growth, molecular, and cellular characteristics of the nasal septum with reference to hyaline cartilage differentiation. The nasal septum grows asynchronous across its length with phases of rapid growth interrupted by more stagnant growth. Growth appears to be driven predominantly by acquisition of chondrocyte hypertrophy. Similarly, cellular differentiation is asynchronous, and differentiation observed in the anterior part precedes posterior differentiation. Overall, the nasal septum is structurally and molecularly heterogeneous. Early and extensive chondrocyte hypertrophy but no ossification is observed in the anterior septum. Onset of hypertrophic chondrocyte differentiation coincided with collagen fiber deposition along the perichondrium. Sox9, Col2, Col10, Mmp13, Sp7, and Runx2 expression was heterogeneous and did not always follow the expected pattern established from chondrocyte differentiation in the growth plate. The presence of hypertrophic chondrocytes expressing bone-related proteins early on in regions where the nasal septum does not ossify displays incongruities with current understanding of hyaline cartilage differentiation. Runx2, Collagen II, Collagen X, and Sp7 commonly used to mark distinct stages of chondrocyte maturation and early bone formation show wider expression than expected and do not align with expected cellular characteristics. Thus, the hyaline cartilage of the nasal septum is quite distinct from growth plate hyaline cartilage, and caution should be taken before assigning cartilage properties to less well-defined cartilage structures using these commonly used markers. Beyond the structural description of the nasal cartilage, this study also provides important information for cartilage tissue engineering when using nasal septal cartilage for tissue regeneration.

Whispers from vestigial nubbins: Arrested development provokes trait loss in toads

Despite the use of acoustic communication, many species of toads (family Bufonidae) have lost parts of the tympanic middle ear, representing at least 12 independent evolutionary occurrences of trait loss. The comparative development of the tympanic middle ear in toads is poorly understood. Here, we compared middle ear development among two pairs of closely related toad species in the genera Atelopus and Rhinella that have (eared) or lack (earless) middle ear structures. We bred toads in Peru and Ecuador, preserved developmental series from tadpoles to juveniles, and examined ontogenetic timing and volume of the otic capsule, oval window, operculum, opercularis muscle, columella (stapes), and extracolumella in three-dimensional histological reconstructions. All species had similar ontogenesis of the otic capsule, oval window, operculum, and opercularis muscle. Moreover, cell clusters of primordial columella in the oval window appeared just before metamorphosis in both eared and earless lineages. However, in earless lineages, the cell clusters either remained as small nubbins or cell buds in the location of the columella footplate within the oval window or disappeared by juvenile and adult stages. Thus, columella growth began around metamorphosis in all species but was truncated and/or degenerated after metamorphosis in earless species, leaving earless adults with morphology typical of metamorphic anurans. Shifts in the timing or expression of biochemical pathways that regulate the extension or differentiation of the columella after metamorphosis may be the developmental mechanism underlying convergent trait loss among toad lineages.

Evolutionary repression of chondrogenic genes in the vertebrate osteoblast

Gene expression in extant animals might reveal how skeletal cells have evolved over the past 500 million years. The cells that make up cartilage (chondrocytes) and bone (osteoblasts) express many of the same genes, but they also have important molecular differences that allow us to distinguish them as separate cell types. For example, traditional studies of later-diverged vertebrates, such as mouse and chick, defined the genes Col2a1 and sex-determining region Y-box 9 as cartilage-specific. However, recent studies have shown that osteoblasts of earlier-diverged vertebrates, such as frog, gar, and zebrafish, express these 'chondrogenic' markers. In this review, we examine the resulting hypothesis that chondrogenic gene expression became repressed in osteoblasts over evolutionary time. The amphibian is an underexplored skeletal model that is uniquely positioned to address this hypothesis, especially given that it diverged when life transitioned from water to land. Given the relationship between phylogeny and ontogeny, a novel discovery for skeletal cell evolution might bolster our understanding of skeletal cell development.

Skeletal Mineralization in Association with Type X Collagen Expression Is an Ancestral Feature for Jawed Vertebrates

In order to characterize the molecular bases of mineralizing cell evolution, we targeted type X collagen, a nonfibrillar network forming collagen encoded by the Col10a1 gene. It is involved in the process of endochondral ossification in ray-finned fishes and tetrapods (Osteichthyes), but until now unknown in cartilaginous fishes (Chondrichthyes). We show that holocephalans and elasmobranchs have respectively five and six tandemly duplicated Col10a1 gene copies that display conserved genomic synteny with osteichthyan Col10a1 genes. All Col10a1 genes in the catshark Scyliorhinus canicula are expressed in ameloblasts and/or odontoblasts of teeth and scales, during the stages of extracellular matrix protein secretion and mineralization. Only one duplicate is expressed in the endoskeletal (vertebral) mineralizing tissues. We also show that the expression of type X collagen is present in teeth of two osteichthyans, the zebrafish Danio rerio and the western clawed frog Xenopus tropicalis, indicating an ancestral jawed vertebrate involvement of type X collagen in odontode formation. Our findings push the origin of Col10a1 gene prior to the divergence of osteichthyans and chondrichthyans, and demonstrate its ancestral association with mineralization of both the odontode skeleton and the endoskeleton.

Scaffold-free bioprinted osteogenic and chondrogenic systems to model osteochondral physiology

Three-dimensional bioprinted culture platforms mimic the native microenvironment of tissues more accurately than two-dimensional cell cultures or animal models. Scaffold-free bioprinting eliminates many complications associated with traditional scaffold-dependent printing as well as provides better cell-to-cell interactions and long-term functionality. In this study, constructs were produced from bone marrow derived mesenchymal stem cells (BM-MSCs) using a scaffold-free bioprinter. These constructs were cultured in either osteogenic, chondrogenic, a 50:50 mixture of osteogenic and chondrogenic ('osteo-chondro'), or BM-MSC growth medium. Osteogenic and chondrogenic differentiation capacity was determined over an 8-week culture period using histological and immunohistochemical staining and RT-qPCR (Phase I). After 6 weeks in culture, individual osteogenic and chondrogenic differentiated constructs were adhered to create a bone-cartilage interaction model. Adhered differentiated constructs were cultured for an additional 8 weeks in either chondrogenic or osteo-chondro medium to evaluate sustainability of lineage specification and transdifferentiation potential (Phase II). Constructs cultured in their respective osteogenic and/or chondrogenic medium differentiated directly into bone (model of intramembranous ossification) or cartilage. Positive histological and immunohistochemical staining for bone or cartilage identification was shown after 4 and 8 weeks in culture. Expression of osteogenesis and chondrogenesis associated genes increased between weeks 2 and 6. Adhered individual osteogenic and chondrogenic differentiated constructs sustained their differentiated phenotype when cultured in chondrogenic medium. However, adhered individual chondrogenic differentiated constructs cultured in osteo-chondro medium were converted to bone (model of metaplastic transformation). These bioprinted models of bone-cartilage interaction, intramembranous ossification, and metaplastic transformation of cartilage into bone offer a useful and promising approach for bone and cartilage tissue engineering research. Specifically, these models can be potentially used as functional tissue systems for studying osteochondral defect repair, drug discovery and response, and many other potential applications.

A single-cell transcriptomic atlas of the developing chicken limb

BACKGROUND

Through precise implementation of distinct cell type specification programs, differentially regulated in both space and time, complex patterns emerge during organogenesis. Thanks to its easy experimental accessibility, the developing chicken limb has long served as a paradigm to study vertebrate pattern formation. Through decades' worth of research, we now have a firm grasp on the molecular mechanisms driving limb formation at the tissue-level. However, to elucidate the dynamic interplay between transcriptional cell type specification programs and pattern formation at its relevant cellular scale, we lack appropriately resolved molecular data at the genome-wide level. Here, making use of droplet-based single-cell RNA-sequencing, we catalogue the developmental emergence of distinct tissue types and their transcriptome dynamics in the distal chicken limb, the so-called autopod, at cellular resolution.

RESULTS

Using single-cell RNA-sequencing technology, we sequenced a total of 17,628 cells coming from three key developmental stages of chicken autopod patterning. Overall, we identified 23 cell populations with distinct transcriptional profiles. Amongst them were small, albeit essential populations like the apical ectodermal ridge, demonstrating the ability to detect even rare cell types. Moreover, we uncovered the existence of molecularly distinct sub-populations within previously defined compartments of the developing limb, some of which have important signaling functions during autopod pattern formation. Finally, we inferred gene co-expression modules that coincide with distinct tissue types across developmental time, and used them to track patterning-relevant cell populations of the forming digits.

CONCLUSIONS

We provide a comprehensive functional genomics resource to study the molecular effectors of chicken limb patterning at cellular resolution. Our single-cell transcriptomic atlas captures all major cell populations of the developing autopod, and highlights the transcriptional complexity in many of its components. Finally, integrating our data-set with other single-cell transcriptomics resources will enable researchers to assess molecular similarities in orthologous cell types across the major tetrapod clades, and provide an extensive candidate gene list to functionally test cell-type-specific drivers of limb morphological diversification.

Exogenous Nitric Oxide Induced Early Mineralization in Rat Bone Marrow Mesenchymal Stem Cells via Activation of Alkaline Phosphatase

Background

Since the low concentration and short-time treatment with sodium nitroprusside (SNP), a nitric oxide (NO)–donor, cause no harm to rat bone marrow mesenchymal stem cells (MSCs), we studied the impact of SNP on MSCs differentiation.

Methods

MSCs were treated with 100 and 1000 µM of SNP for 1 hour in every 48 hours and after 5, 10, 15, and 21 days in osteogenic media. The viability and the level of mineralization were determined using MTT assay and alizarin red staining, respectively. Morphology of the cells was studied using fluorescent dye. Concentration of calcium and the activity of alanine transaminase (ALT), aspartate transaminase (AST), lactate dehydrogenase (LDH), and alkaline phosphatase (ALP) were evaluated by commercial kits.

Results

SNP with the concentration of 1000 µM significantly reduced viability from day 5 to day 20, but 100 µM did not affect the viability until the day 15. The low concentration of SNP increased matrix deposition from day 10 and reached almost its maximum (4.40 ± 2.4) at the day 15. Also, increasing the activity of ALP (419 ± 2.2), due to low concentration of SNP, started at day 10 and continued till the day 20, while LDH (2026 ± 11) and AST (25.6 ± 0.4) elevations were observed from day 5 onwards. In case of ALT, we observed a significant decrease (36%) from day 5 till day 20.

Conclusion

Based on our findings, low concentrations of SNP might be useful in the promotion of bone repair.

The Forgotten Skeletogenic Condensations: A Comparison of Early Skeletal Development Amongst Vertebrates

The development of a skeletogenic condensation is perhaps the most critical yet considerably overlooked stage of skeletogenesis. Described in this comprehensive review are the mechanisms that facilitate skeletogenic condensation formation, growth, and maintenance to allow for overt differentiation into a skeletal element. This review discusses the current knowledge of gene regulation and characterization of skeletogenic condensations in the chicken, mouse, zebrafish, and other developmental models. We limited our scope to condensations that give rise to the bones and cartilages of the vertebrate skeleton, with a particular focus on craniofacial and limb bud regions. While many of the skeletogenic processes are similar among vertebrate lineages, differences are apparent in the site and timing of the initial epithelial⁻mesenchymal interactions as well as in whether the condensation has an osteogenic or chondrogenic fate, both within and among species. Further comparative studies are needed to clarify and broaden the existing knowledge of this intricate phenomenon.

Skeletal development of the chondrocranium in the tongue sole Cynoglossus semilaevis (Pleuronectiformes: Cynoglossidae)

This study provides a comprehensive description of chondrocranial development before, during and after larval metamorphosis in the tongue sole Cynoglossus semilaevis, a commercially valuable flatfish in China. Samples were collected at regular intervals ranging from 1 to 23 days post hatching (dph). Based on observations of cleared and double-stained specimens and images from sections stained with safranin O-fast green, major morphological events during early development were described. No cartilaginous structure was visible at hatching. From 2 dph onwards, cartilaginous structures such as the trabecular bar and some elements of the mandibular, hyoid and branchial arches appeared. At this time also, cartilaginous structures of the neurocranium started to form. Hypertrophic chondrocytes could be observed in many splanchnocranium elements at 5 dph. The start of ossification was indicated by alizarin red stain visible at 14 dph. At 17 dph, most of the cartilaginous skeleton was ossified. Soon after, the right eye started to migrate and pass through a slit beneath the dorsal-fin base and above the skull. Metamorphosis was complete at 20 dph, at which time the dorsal-fin base cartilage extended onto the anterior region of the head. Meanwhile, extremities of the hyoid and branchial arch elements remained cartilaginous. At 23 dph, endochondral ossification of the splanchnocranium was nearly complete. Unlike previous observations of other Pleuronectiformes, our study indicates that endochondral ossification of C. semilaevis skull cartilage occurs before metamorphosis.

Enzymatically Cross-Linked Silk Fibroin-Based Hierarchical Scaffolds for Osteochondral Regeneration

Osteochondral (OC) regeneration faces several limitations in orthopedic surgery, owing to the complexity of the OC tissue that simultaneously entails the restoration of articular cartilage and subchondral bone diseases. In this study, novel biofunctional hierarchical scaffolds composed of a horseradish peroxidase (HRP)-cross-linked silk fibroin (SF) cartilage-like layer (HRP-SF layer) fully integrated into a HRP-SF/ZnSr-doped β-tricalcium phosphate (β-TCP) subchondral bone-like layer (HRP-SF/dTCP layer) were proposed as a promising strategy for OC tissue regeneration. For comparative purposes, a similar bilayered structure produced with no ion incorporation (HRP-SF/TCP layer) was used. A homogeneous porosity distribution was achieved throughout the scaffolds, as shown by micro-computed tomography analysis. The ion-doped bilayered scaffolds presented a wet compressive modulus (226.56 ± 60.34 kPa) and dynamic mechanical properties (ranging from 403.56 ± 111.62 to 593.56 ± 206.90 kPa) superior to that of the control bilayered scaffolds (189.18 ± 90.80 kPa and ranging from 262.72 ± 59.92 to 347.68 ± 93.37 kPa, respectively). Apatite crystal formation, after immersion in simulated body fluid (SBF), was observed in the subchondral bone-like layers for the scaffolds incorporating TCP powders. Human osteoblasts (hOBs) and human articular chondrocytes (hACs) were co-cultured onto the bilayered structures and monocultured in the respective cartilage and subchondral bone half of the partitioned scaffolds. Both cell types showed good adhesion and proliferation in the scaffold compartments, as well as adequate integration of the interface regions. Osteoblasts produced a mineralized extracellular matrix (ECM) in the subchondral bone-like layers, and chondrocytes showed GAG deposition. The gene expression profile was different in the distinct zones of the bilayered constructs, and the intermediate regions showed pre-hypertrophic chondrocyte gene expression, especially on the BdTCP constructs. Immunofluorescence analysis supported these observations. This study showed that the proposed bilayered scaffolds allowed a specific stimulation of the chondrogenic and osteogenic cells in the co-culture system together with the formation of an osteochondral-like tissue interface. Hence, the structural adaptability, suitable mechanical properties, and biological performance of the hierarchical scaffolds make these constructs a desired strategy for OC defect regeneration.

Bone-like features in skate suggest a novel elasmobranch synapomorphy and deep homology of trabecular mineralization patterns

Bone is a defining characteristic of the vertebrate skeleton, and while chondrichthyans (sharks, skates, and other cartilaginous fishes) are vertebrates, they are hypothesized to have lost the ability to make bone during their evolution. Multiple descriptions of a bone-like tissue in neural arches of vertebrae in various shark species (selachians), however, challenge this hypothesis. Here, we extend this argument by analyzing vertebrae of two members of the batoids (the little skate Leucoraja erinacea and Eaton's skate Bathyraja eatonii), the sister group to selachians within elasmobranchs. Micro-CT images showed a bone-like mineralization pattern in neural arches of each skate species, and histological analyses confirmed that this bone-like tissue surrounded a cartilage core, exactly as described in sharks. Another mineralization pattern identified in skate vertebrae was distinct from the polygonal tesseral and areolar patterns that classically are associated with the chondrichthyan endoskeleton. Many regions of the vertebrae, including the neural spine and transverse processes, showed this perichondral mineralization pattern, termed here trabecular tesseral. Other than the cartilage core of the neural arch, all mineralized tissues in skate vertebrae had flattened cells surrounded by matrix with bone-like histology. Analyses of quantitative microstructural parameters revealed that, compared to rat vertebrae, the bone-like mineralization pattern in the neural arches of skate vertebrae was more similar to compact bone than trabecular bone. In contrast, the thickness of the trabecular tesseral pattern was more similar to trabecular bone than compact bone of rat vertebrae. In conclusion, a bone-like tissue in neural arches of skate vertebrae appears to be a novel elasmobranch synapomorphy. We propose that the trabecular tesseral mineralization pattern in the skate might have deep homology to the mineralization pattern utilized in trabecular bone. STATEMENT OF SIGNIFICANCE: Mineralization patterns of skeletal tissues have not been investigated thoroughly in all vertebrate clades. Despite their designation as 'cartilaginous fish', chondrichthyans clearly evolved from ancestral vertebrates that made bone. The consensus that chondrichthyans lost the ability to make bone during their evolution, however, is challenged by reports of bone and bone-like tissues in the neural arches of vertebrae in extant sharks (selachians). Here, we provide evidence from micro-CT imaging and histological analyses to support our hypothesis that a bone-like tissue is present in the neural arches of batoids (the sister group to selachians within elasmobranchs). These results argue strongly that the neural arch bone-like tissue is a previously unknown synapomorphy of elasmobranchs. In addition to the bone-like mineralization pattern identified in the neural arches, micro-CT images also showed a novel mineralization pattern which we described as trabecular tesseral. Quantitative microstructural features shared between trabecular tesseral pattern and trabecular bone (from homologous rat vertebrae) suggest that both patterns might derive from an ancestral gene network driving trabecular mineralization (i.e., deep homology).

Chasing Chimeras - The elusive stable chondrogenic phenotype

The choice of the best-suited cell population for the regeneration of damaged or diseased cartilage depends on the effectiveness of culture conditions (e.g. media supplements, three-dimensional scaffolds, mechanical stimulation, oxygen tension, co-culture systems) to induce stable chondrogenic phenotype. Herein, advances and shortfalls in in vitro, preclinical and clinical setting of various in vitro microenvironment modulators on maintaining chondrocyte phenotype or directing stem cells towards chondrogenic lineage are critically discussed. Chondrocytes possess low isolation efficiency, limited proliferative potential and rapid phenotypic drift in culture. Mesenchymal stem cells are relatively readily available, possess high proliferation potential, exhibit great chondrogenic differentiation capacity, but they tend to acquire a hypertrophic phenotype when exposed to chondrogenic stimuli. Embryonic and induced pluripotent stem cells, despite their promising in vitro and preclinical data, are still under-investigated. Although a stable chondrogenic phenotype remains elusive, recent advances in in vitro microenvironment modulators are likely to develop clinically- and commercially-relevant therapies in the years to come.

Simulation enabled search for explanatory mechanisms of the fracture healing process

A significant portion of bone fractures fail to heal properly, increasing healthcare costs. Advances in fracture management have slowed because translation barriers have limited generation of mechanism-based explanations for the healing process. When uncertainties are numerous, analogical modeling can be an effective strategy for developing plausible explanations of complex phenomena. We demonstrate the feasibility of engineering analogical models in software to facilitate discovery of biomimetic explanations for how fracture healing may progress. Concrete analogical models—Callus Analogs—were created using the MASON simulation toolkit. We designated a Target Region initial state within a characteristic tissue section of mouse tibia fracture at day-7 and posited a corresponding day-10 Target Region final state. The goal was to discover a coarse-grain analog mechanism that would enable the discretized initial state to transform itself into the corresponding Target Region final state, thereby providing an alternative way to study the healing process. One of nine quasi-autonomous Tissue Unit types is assigned to each grid space, which maps to an 80×80 μm region of the tissue section. All Tissue Units have an opportunity each time step to act based on individualized logic, probabilities, and information about adjacent neighbors. Action causes transition from one Tissue Unit type to another, and simulation through several thousand time steps generates a coarse-grain analog—a theory—of the healing process. We prespecified a minimum measure of success: simulated and actual Target Region states achieve ≥ 70% Similarity. We used an iterative refinement protocol to explore many combinations of Tissue Unit logic and action constraints. Workflows progressed through four stages of analog mechanisms. Similarities of 73–90% were achieved for Mechanisms 2–4. The range of Upper-Level similarities increased to 83–94% when we allowed for uncertainty about two Tissue Unit designations. We have demonstrated how Callus Analog experiments provide domain experts with a fresh medium and tools for thinking about and understanding the fracture healing process.

Translation barriers have limited the generation of mechanism-based explanations of fracture healing processes. Those barriers help explain why, to date, biological therapeutics have had only a minor impact on fracture management. Alternative approaches are needed, and we present one that is intended to help develop incrementally better mechanism-based explanations of fracture healing phenomena. We created virtual Callus Analogs to simulate how the histologic appearance of a mouse fracture callus may transition from day-7 to day-10. Callus Analogs use software-based model mechanisms, and simulation experiments enable challenging and improving those model mechanisms. During execution, model mechanism operation provides a coarse-grain explanation (a theory) of a four-day portion of the healing process. Simulated day-10 callus histologic images achieved 73–94% Similarity to a corresponding day-10 fracture callus image, thus demonstrating feasibility. Simulated healing provides an alternative perspective on the actual healing process and an alternative way of thinking about plausible fracture healing mechanisms. Our working hypothesis is that the approach can be extended to cover more of the healing process while making features of simulated and actual fracture healing increasingly analogous. The methods presented are intended to be extensible to other research areas that use histologic analysis to investigate and explain tissue level phenomena.

Posttranscriptional Interaction Between miR-450a-5p and miR-28-5p and STAT1 mRNA Triggers Osteoblastic Differentiation of Human Mesenchymal Stem Cells

We demonstrate that the interaction between miR-450a-5p and miR-28-5p and signal transducer and activator of transcription 1 (STAT1) mRNA correlates with the osteoblastic differentiation of mesenchymal stem cells from human exfoliated deciduous teeth (shed cells). STAT1 negatively regulates runx-related transcription factor 2 (RUNX2), which is an essential transcription factor in this process. However, the elements that trigger osteoblastic differentiation and therefore pause the inhibitory effect of STAT1 need investigation. Usually, STAT1 can be posttranscriptionally regulated by miRNAs. To test this, we used an in vitro model system in which shed cells were chemically induced toward osteoblastic differentiation and temporally analyzed, comparing undifferentiated cells with their counterparts in the early (2 days) or late (7 or 21 days) periods of induction. The definition of the entire functional genome expression signature demonstrated that the transcriptional activity of a large set of mRNAs and miRNAs changes during this process. Interestingly, STAT1 and RUNX2 mRNAs feature contrasting expression levels during the course of differentiation. While undifferentiated or early differentiating cells express high levels of STAT1 mRNA, which was gradually downregulated, RUNX2 mRNA was upregulated toward differentiation. The reconstruction of miRNA-mRNA interaction networks allowed the identification of six miRNAs (miR-17-3p, miR-28-5p, miR-29b, miR-29c-5p, miR-145-3p, and miR-450a-5p), and we predicted their respective targets, from which we focused on miR-450a-5p and miR-28-5p STAT1 mRNA interactions, whose intracellular occurrence was validated through the luciferase assay. Transfections of undifferentiated shed cells with miR-450a-5p or miR-28-5p mimics or with miR-450a-5p or miR-28-5p antagonists demonstrated that these miRNAs might play a role as posttranscriptional controllers of STAT1 mRNA during osteoblastic differentiation. J. Cell. Biochem. 118: 4045-4062, 2017. © 2017 Wiley Periodicals, Inc.

Evolutionary origin of endochondral ossification: the transdifferentiation hypothesis

The vertebrate endoskeleton results from the piecemeal assembly of bone and cartilage as well as additional types of calcified extracellular matrices produced by seemingly hybrid cell types of intermediate phenotypes between osteoblasts and chondrocytes. Hence, shedding light on the emergence and subsequent diversification of skeletal tissues represents a major challenge in vertebrate evolutionary developmental biology. A 150-year-old debate in the field was recently solved by lineage tracing experiments demonstrating that, during mouse endochondral bone development, a subset of chondrocytes evades apoptosis and transdifferentiates into osteoblasts at the chondro-osseous junction. Here, we interpret these new data from a broad phylogenetic perspective, integrating fossil, histological, cellular, and genetic evidence from a wide range of vertebrates. We propose a testable scenario according to which transdifferentiation played a fundamental role in the emergence of endochondral ossification, an osteichthyan-specific evolutionary novelty. The remarkable skeletal cell plasticity might be contingent on the similar architectures of the osteoblastic and chondrocytic gene regulatory networks, thereby providing a unifying mechanism underlying both complete transdifferentiation as well as partial cell type conversions observed in intermediate skeletal tissues such as the chondroid bone or globuli ossei.

Lessons on skeletal cell plasticity from studying jawbone regeneration in zebrafish

Three major mesenchymal cell types have important roles in determining the shapes of vertebrate animals: bone-producing osteoblasts, cartilage-producing chondrocytes, and fat-producing adipocytes. Although often considered discrete cell types, accumulating evidence is revealing mesenchymal cells of intermediate identities and interconversion of cell types. Such plasticity is particularly evident during adult skeletal repair. In this Review, we highlight recent work in zebrafish showing a role for hybrid cartilage-bone cells in large-scale regeneration of the adult jawbone, as well as their origins in the periosteum. An emerging theme is that the unique mechanical and signaling environment of the adult wound causes skeletal cell differentiation to diverge from the discrete lineages seen during development, which may aid in rapid and extensive regeneration of bone.

Signaling pathways regulating cartilage growth plate formation and activity

The growth plate is a highly specialized and dynamic cartilage structure that serves many essential functions in skeleton patterning, growth and endochondral ossification in developing vertebrates. Major signaling pathways initiated by classical morphogens and by other systemic and tissue-specific factors are intimately involved in key aspects of growth plate development. As a corollary of these essential functions, disturbances in these pathways due to mutations or environmental factors lead to severe skeleton disorders. Here, we review these pathways and the most recent progress made in understanding their roles in chondrocyte differentiation in growth plate development and activity. Furthermore, we discuss newly uncovered pathways involved in growth plate formation, including mTOR, the circadian clock, and the COP9 signalosome.