Evolution of the osteoblast: skeletogenesis in gar and zebrafish

CHD7 regulates craniofacial cartilage development via controlling HTR2B expression

Mutations in the Chromodomain helicase DNA-binding protein 7 - coding gene (CHD7) cause CHARGE syndrome (CS). Although craniofacial and skeletal abnormalities are major features of CS patients, the role of CHD7 in bone and cartilage development remain largely unexplored. Here, using a zebrafish (Danio rerio) CS model, we show that chd7-/- larvae display abnormal craniofacial cartilage development and spinal deformities. The craniofacial and spine defects are accompanied by a marked reduction of bone mineralization. At the molecular level, we show that these phenotypes are associated with significant reduction in the expression levels of osteoblast differentiation markers. Additionally, we detected a marked depletion of collagen 2α1 in the cartilage of craniofacial regions and vertebrae, along with significantly reduced number of chondrocytes. Chondrogenesis defects are at least in part due to downregulation of htr2b, which we found to be also dysregulated in human cells derived from an individual with CHD7 mutation-positive CS. Overall, this study thus unveils an essential role for CHD7 in cartilage and bone development, with potential clinical relevance for the craniofacial defects associated with CS.

A Collagen10a1 mutation disrupts cell polarity in a medaka model for metaphyseal chondrodysplasia type Schmid

Heterozygous mutations in COL10A1 lead to metaphyseal chondrodysplasia type Schmid (MCDS), a skeletal disorder characterized by epiphyseal abnormalities. Prior analysis revealed impaired trimerization and intracellular retention of mutant collagen type X alpha 1 chains as cause for elevated endoplasmic reticulum (ER) stress. However, how ER stress translates into structural defects remained unclear. We generated a medaka (Oryzias latipes) MCDS model harboring a 5 base pair deletion in col10a1, which led to a frameshift and disruption of 11 amino acids in the conserved trimerization domain. col10a1Δ633a heterozygotes recapitulated key features of MCDS and revealed early cell polarity defects as cause for dysregulated matrix secretion and deformed skeletal structures. Carbamazepine, an ER stress-reducing drug, rescued this polarity impairment and alleviated skeletal defects in col10a1Δ633a heterozygotes. Our data imply cell polarity dysregulation as a potential contributor to MCDS and suggest the col10a1Δ633a medaka mutant as an attractive MCDS animal model for drug screening.

The Osteoblast Transcriptome in Developing Zebrafish Reveals Key Roles for Extracellular Matrix Proteins Col10a1a and Fbln1 in Skeletal Development and Homeostasis

Zebrafish are now widely used to study skeletal development and bone-related diseases. To that end, understanding osteoblast differentiation and function, the expression of essential transcription factors, signaling molecules, and extracellular matrix proteins is crucial. We isolated Sp7-expressing osteoblasts from 4-day-old larvae using a fluorescent reporter. We identified two distinct subpopulations and characterized their specific transcriptome as well as their structural, regulatory, and signaling profile. Based on their differential expression in these subpopulations, we generated mutants for the extracellular matrix protein genes col10a1a and fbln1 to study their functions. The col10a1a-/- mutant larvae display reduced chondrocranium size and decreased bone mineralization, while in adults a reduced vertebral thickness and tissue mineral density, and fusion of the caudal fin vertebrae were observed. In contrast, fbln1-/- mutants showed an increased mineralization of cranial elements and a reduced ceratohyal angle in larvae, while in adults a significantly increased vertebral centra thickness, length, volume, surface area, and tissue mineral density was observed. In addition, absence of the opercle specifically on the right side was observed. Transcriptomic analysis reveals up-regulation of genes involved in collagen biosynthesis and down-regulation of Fgf8 signaling in fbln1-/- mutants. Taken together, our results highlight the importance of bone extracellular matrix protein genes col10a1a and fbln1 in skeletal development and homeostasis.

Contribution of Signal Transducer and Activator of Transcription 3 (STAT3) to Bone Development and Repair

Signal transducer and activator of transcription 3 (STAT3) is a transcription factor activated canonically by numerous cytokines and other factors, with significant roles in immunity, immune diseases, and cancer. It has also been implicated in several human skeletal disorders, with loss-of-function (LOF) mutations associated with aberrant skeletal development. To gain further insights, two zebrafish STAT3 lines were investigated: a complete LOF knockout (KO) mutant and a partial LOF mutant with the transactivation domain truncated (ΔTAD). Consistent with other studies, the KO mutants were smaller, with reduced length in early embryos exacerbated by a decreased growth rate from 5 days postfertilization (dpf). They displayed skeletal deformities that approached 80% incidence by 30 dpf, with a significant reduction in early bone but not cartilage formation. Further analysis additionally identified considerable abrogation of caudal fin regeneration, concomitant with a paucity of infiltrating macrophages and neutrophils, which may be responsible for this. Most of these phenotypes were also observed in the ΔTAD mutants, indicating that loss of canonical STAT3 signaling was the likely cause. However, the impacts on early bone formation and regeneration were muted in the ΔTAD mutant, suggesting the potential involvement of noncanonical functions in these processes.

Identification and expression analysis of transcripts involved in taurine biosynthesis during early ontogeny of tropical gar Atractosteus tropicus

In fishes, the availability of taurine is regulated during ontogenetic development, where its endogenous synthesis capacity is species dependent. Thus, different pathways and involved enzymes have been described: pathway I (cysteine sulfinate-dependent pathway), cysteine dioxygenase type 1 (cdo1) and cysteine sulfinic acid decarboxylase (csad); pathway II (cysteic acid pathway), cdo1 and glutamic acid decarboxylase (gad); and pathway III (cysteamine pathway), 2-aminoethanethiol dioxygenase (ado); whereas taurine transporter (taut) is responsible for taurine entry into cells on the cell membrane and the mitochondria. This study determined if the tropical gar (Atractosteus tropicus), an ancient holostean fish model, has the molecular mechanism to synthesize taurine through the identification and analysis expression of transcripts coding for proteins involved in its biosynthesis and transportation, at different embryo-larvae stages and in different organs of juveniles (31 dah). We observed a fluctuating expression of all transcripts involved in the three pathways at all analyzed stages. All transcripts are expressed during the beginning of larval development; however, ado and taut show a peak expression at 9 dah, and all transcripts but csad decreased at 23 dah, when the organism ended the larval period. Furthermore, at 31 dah, we observed taut expression in all examined organs. The transcripts involved in pathways I and III are expressed differently across all organs, whereas pathway II was only observed in the brain, eye, and skin. The results suggested that taurine biosynthesis in tropical gar is regulated during its early development before first feeding, and the pathway might also be organ-type dependent.

T, Mishoe Hernandez L, DeLaurier A. Exposure to valproic acid (VPA) reproduces hdac1 loss of function phenotypes in zebrafish

Histone deacetylases are enzymes that remove acetyl groups from histone tails and are understood to act as repressors of transcriptional activity. Hdac1 has been previously shown to function in eye, pectoral fin, heart, liver, and pharyngeal skeletal development. We show that high doses of Valproic Acid (VPA) reproduce the hdac1 phenotype. We identify tbx5 genes as potential targets of Hdac1 in eye, pectoral fin, and heart development. Using timed exposures, we show that skeletal structures in the pharyngeal arches are impacted by VPA between 24-36 hours post-fertilization, indicating a role for Hdac1 during post-migration patterning, differentiation, or proliferation of cranial neural crest cells.

Disruption of the foxe1 gene in zebrafish reveals conserved functions in development of the craniofacial skeleton and the thyroid

Introduction: Mutations in the FOXE1 gene are implicated in cleft palate and thyroid dysgenesis in humans.

Methods: To investigate whether zebrafish could provide meaningful insights into the etiology of developmental defects in humans related to FOXE1, we generated a zebrafish mutant that has a disruption in the nuclear localization signal in the foxe1 gene, thereby restraining nuclear access of the transcription factor. We characterized skeletal development and thyroidogenesis in these mutants, focusing on embryonic and larval stages.

Results: Mutant larvae showed aberrant skeletal phenotypes in the ceratohyal cartilage and had reduced whole body levels of Ca, Mg and P, indicating a critical role for foxe1 in early skeletal development. Markers of bone and cartilage (precursor) cells were differentially expressed in mutants in post-migratory cranial neural crest cells in the pharyngeal arch at 1 dpf, at induction of chondrogenesis at 3 dpf and at the start of endochondral bone formation at 6 dpf. Foxe1 protein was detected in differentiated thyroid follicles, suggesting a role for the transcription factor in thyroidogenesis, but thyroid follicle morphology or differentiation were unaffected in mutants.

Discussion: Taken together, our findings highlight the conserved role of Foxe1 in skeletal development and thyroidogenesis, and show differential signaling of osteogenic and chondrogenic genes related to foxe1 mutation.

Meis2 controls skeletal formation in the hyoid region

A vertebrate skull is composed of many skeletal elements which display enormous diversity of shapes. Cranial bone formation embodies a multitude of processes, i.e., epithelial-mesenchymal induction, mesenchymal condensation, and endochondral or intramembranous ossification. Molecular pathways determining complex architecture and growth of the cranial skeleton during embryogenesis are poorly understood. Here, we present a model of the hyoid apparatus development in Wnt1-Cre2-induced Meis2 conditional knock-out (cKO) mice. Meis2 cKO embryos develop an aberrant hyoid apparatus—a complete skeletal chain from the base of the neurocranium to lesser horns of the hyoid, resembling extreme human pathologies of the hyoid-larynx region. We examined key stages of hyoid skeletogenesis to obtain a complex image of the hyoid apparatus formation. Lack of Meis2 resulted in ectopic loci of mesenchymal condensations, ectopic cartilage and bone formation, disinhibition of skeletogenesis, and elevated proliferation of cartilage precursors. We presume that all these mechanisms contribute to formation of the aberrant skeletal chain in the hyoid region. Moreover, Meis2 cKO embryos exhibit severely reduced expression of PBX1 and HAND2 in the hyoid region. Altogether, MEIS2 in conjunction with PBX1 and HAND2 affects mesenchymal condensation, specification and proliferation of cartilage precursors to ensure development of the anatomically correct hyoid apparatus.

The Shape of the Jaw-Zebrafish Col11a1a Regulates Meckel's Cartilage Morphogenesis and Mineralization

The expression of the col11a1a gene is essential for normal skeletal development, affecting both cartilage and bone. Loss of function mutations have been shown to cause abnormalities in the growth plate of long bones, as well as in craniofacial development. However, the specific effects on Meckel's cartilage have not been well studied. To further understand the effect of col11a1a gene function, we analyzed the developing jaw in zebrafish using gene knockdown by the injection of an antisense morpholino oligonucleotide using transgenic Tg(sp7:EGFP) and Tg(Fli1a:EGFP) EGFP reporter fish, as well as wildtype AB zebrafish. Our results demonstrate that zebrafish col11a1a knockdown impairs the cellular organization of Meckel's cartilage in the developing jaw and alters the bone formation that occurs adjacent to the Meckel's cartilage. These results suggest roles for Col11a1a protein in cartilage intermediates of bone development, the subsequent mineralization of the bony collar of long bones, and that which occurs adjacent to Meckel's cartilage in the developing jaw.

A non-disruptive and efficient knock-in approach allows fate tracing of resident osteoblast progenitors during repair of vertebral lesions in medaka

During bone development and repair, osteoblasts are recruited to bone deposition sites. To identify the origin of recruited osteoblasts, cell lineage tracing using Cre/loxP recombination is commonly used. However, a confounding factor is the use of transgenic Cre drivers that do not accurately recapitulate endogenous gene expression or the use of knock-in Cre drivers that alter endogenous protein activity or levels. Here, we describe a CRISPR/Cas9 homology-directed repair knock-in approach that allows efficient generation of Cre drivers controlled by the endogenous gene promoter. In addition, a self-cleaving peptide preserves the reading frame of the endogenous protein. Using this approach, we generated col10a1p2a-CreERT2 knock-in medaka and show that tamoxifen-inducible CreERT2 efficiently recombined loxP sites in col10a1 cells. Similar knock-in efficiencies were obtained when two unrelated loci (osr1 and col2a1a) were targeted. Using live imaging, we traced the fate of col10a1 osteoblast progenitors during bone lesion repair in the medaka vertebral column. We show that col10a1 cells at neural arches represent a mobilizable cellular source for bone repair. Together, our study describes a previously unreported strategy for precise cell lineage tracing via efficient and non-disruptive knock-in of Cre.

A Comprehensive Analysis of Fibrillar Collagens in Lamprey Suggests a Conserved Role in Vertebrate Musculoskeletal Evolution

Vertebrates have distinct tissues which are not present in invertebrate chordates nor other metazoans. The rise of these tissues also coincided with at least one round of whole-genome duplication as well as a suite of lineage-specific segmental duplications. Understanding whether novel genes lead to the origin and diversification of novel cell types, therefore, is of great importance in vertebrate evolution. Here we were particularly interested in the evolution of the vertebrate musculoskeletal system, the muscles and connective tissues that support a diversity of body plans. A major component of the musculoskeletal extracellular matrix (ECM) is fibrillar collagens, a gene family which has been greatly expanded upon in vertebrates. We thus asked whether the repertoire of fibrillar collagens in vertebrates reflects differences in the musculoskeletal system. To test this, we explored the diversity of fibrillar collagens in lamprey, a jawless vertebrate which diverged from jawed vertebrates (gnathostomes) more than five hundred million years ago and has undergone its own gene duplications. Some of the principal components of vertebrate hyaline cartilage are the fibrillar collagens type II and XI, but their presence in cartilage development across all vertebrate taxa has been disputed. We particularly emphasized the characterization of genes in the lamprey hyaline cartilage, testing if its collagen repertoire was similar to that in gnathostomes. Overall, we discovered thirteen fibrillar collagens from all known gene subfamilies in lamprey and were able to identify several lineage-specific duplications. We found that, while the collagen loci have undergone rearrangement, the Clade A genes have remained linked with the hox clusters, a phenomenon also seen in gnathostomes. While the lamprey muscular tissue was largely similar to that seen in gnathostomes, we saw considerable differences in the larval lamprey skeletal tissue, with distinct collagen combinations pertaining to different cartilage types. Our gene expression analyses were unable to identify type II collagen in the sea lamprey hyaline cartilage nor any other fibrillar collagen during chondrogenesis at the stages observed, meaning that sea lamprey likely no longer require these genes during early cartilage development. Our findings suggest that fibrillar collagens were multifunctional across the musculoskeletal system in the last common ancestor of vertebrates and have been largely conserved, but these genes alone cannot explain the origin of novel cell types.

Regenerating zebrafish scales express a subset of evolutionary conserved genes involved in human skeletal disease

Background

Scales are mineralised exoskeletal structures that are part of the dermal skeleton. Scales have been mostly lost during evolution of terrestrial vertebrates whilst bony fish have retained a mineralised dermal skeleton in the form of fin rays and scales. Each scale is a mineralised collagen plate that is decorated with both matrix-building and resorbing cells. When removed, an ontogenetic scale is quickly replaced following differentiation of the scale pocket-lining cells that regenerate a scale. Processes promoting de novo matrix formation and mineralisation initiated during scale regeneration are poorly understood. Therefore, we performed transcriptomic analysis to determine gene networks and their pathways involved in dermal scale regeneration.

Results

We defined the transcriptomic profiles of ontogenetic and regenerating scales of zebrafish and identified 604 differentially expressed genes (DEGs). These were enriched for extracellular matrix, ossification, and cell adhesion pathways, but not in enamel or dentin formation processes indicating that scales are reminiscent to bone. Hypergeometric tests involving monogenetic skeletal disorders showed that DEGs were strongly enriched for human orthologues that are mutated in low bone mass and abnormal bone mineralisation diseases (P< 2× 10⁻³). The DEGs were also enriched for human orthologues associated with polygenetic skeletal traits, including height (P< 6× 10⁻⁴), and estimated bone mineral density (eBMD, P< 2× 10⁻⁵). Zebrafish mutants of two human orthologues that were robustly associated with height (COL11A2, P=6× 10⁻²⁴) or eBMD (SPP1, P=6× 10⁻²⁰) showed both exo- and endo- skeletal abnormalities as predicted by our genetic association analyses; col11a2^Y228X/Y228X mutants showed exoskeletal and endoskeletal features consistent with abnormal growth, whereas spp1^P160X/P160X mutants predominantly showed mineralisation defects.

Conclusion

We show that scales have a strong osteogenic expression profile comparable to other elements of the dermal skeleton, enriched in genes that favour collagen matrix growth. Despite the many differences between scale and endoskeletal developmental processes, we also show that zebrafish scales express an evolutionarily conserved sub-population of genes that are relevant to human skeletal disease.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01209-8.

Larval Development in Tropical Gar (Atractosteus tropicus) Is Dependent on the Embryonic Thermal Regime: Ecological Implications under a Climate Change Context

In ectotherm species, environmental temperature plays a key role in development, growth, and survival. Thus, determining how temperature affects fish populations is of utmost importance to accurately predict the risk of climate change over fisheries and aquaculture, critical to warrant nutrition and food security in the coming years. Here, the potential effects of abnormal thermal regimes (24, 28 and 32 °C; TR24, TR28, and TR32, respectively) exclusively applied during embryogenesis in tropical gar (Atractosteus tropicus) has been explored to decipher the potential consequences on hatching and growth from fertilization to 16 days post-fertilization (dpf), while effects on skeletal development and body morphology were explored at fertilization and 16 dpf. Egg incubation at higher temperatures induced an early hatching and mouth opening. A higher hatching rate was obtained in eggs incubated at 28 °C when compared to those at 24 °C. No differences were found in fish survival at 16 dpf, with values ranging from 84.89 to 88.86%, but increased wet body weight and standard length were found in larvae from TR24 and TR32 groups. Thermal regime during embryogenesis also altered the rate at which the skeletal development occurs. Larvae from the TR32 group showed an advanced skeletal development, with a higher development of cartilaginous structures at hatching but reduced at 16 dpf when compared with the TR24 and TR28 groups. Furthermore, this advanced skeletal development seemed to determine the fish body morphology. Based on biometric measures, a principal component analysis showed how along development, larvae from each thermal regime were clustered together, but with each population remaining clearly separated from each other. The current study shows how changes in temperature may induce craniofacial and morphological alterations in fish during early stages and contribute to understanding the possible effects of global warming in early development of fish and its ecological implications.

Mineralized Cartilage and Bone-Like Tissues in Chondrichthyans Offer Potential Insights Into the Evolution and Development of Mineralized Tissues in the Vertebrate Endoskeleton

The impregnation of biominerals into the extracellular matrix of living organisms, a process termed biomineralization, gives rise to diverse mineralized (or calcified) tissues in vertebrates. Preservation of mineralized tissues in the fossil record has provided insights into the evolutionary history of vertebrates and their skeletons. However, current understanding of the vertebrate skeleton and of the processes underlying its formation is biased towards biomedical models such as the tetrapods mouse and chick. Chondrichthyans (sharks, skates, rays, and chimaeras) and osteichthyans are the only vertebrate groups with extant (living) representatives that have a mineralized skeleton, but the basal phylogenetic position of chondrichthyans could potentially offer unique insights into skeletal evolution. For example, bone is a vertebrate novelty, but the internal supporting skeleton (endoskeleton) of extant chondrichthyans is commonly described as lacking bone. The molecular and developmental basis for this assertion is yet to be tested. Subperichondral tissues in the endoskeleton of some chondrichthyans display mineralization patterns and histological and molecular features of bone, thereby challenging the notion that extant chondrichthyans lack endoskeletal bone. Additionally, the chondrichthyan endoskeleton demonstrates some unique features and others that are potentially homologous with other vertebrates, including a polygonal mineralization pattern, a trabecular mineralization pattern, and an unconstricted perichordal sheath. Because of the basal phylogenetic position of chondrichthyans among all other extant vertebrates with a mineralized skeleton, developmental and molecular studies of chondrichthyans are critical to flesh out the evolution of vertebrate skeletal tissues, but only a handful of such studies have been carried out to date. This review discusses morphological and molecular features of chondrichthyan endoskeletal tissues and cell types, ultimately emphasizing how comparative embryology and transcriptomics can reveal homology of mineralized skeletal tissues (and their cell types) between chondrichthyans and other vertebrates.

Detection of selection signatures in the genome of a farmed population of anadromous rainbow trout (Oncorhynchus mykiss)

Domestication processes and artificial selection are likely to leave signatures that can be detected at a molecular level in farmed rainbow trout (Oncorhynchus mykiss). These signatures of selection are genomic regions that contain functional genetic variants conferring a higher fitness to their bearers. We genotyped 749 rainbow trout from a commercial population using a rainbow trout Axiom 57 K SNP array panel and identified putative genomic regions under selection using the pcadapt, Composite Likelihood Ratio (CLR) and Integrated Haplotype Score (iHS) methods. After applying quality-control pipelines and statistical analyses, we detected 12, 96 and 16 SNPs putatively under selection, associated with 96, 781 and 115 candidate genes, respectively. Several of these candidate genes were associated with growth, early development, reproduction, behavior and immune system traits. In addition, some of the SNPs were found in interesting regions located in autosomal inversions on Omy05 and Omy20. These findings could represent a genome-wide map of selection signatures in farmed rainbow trout and could be important in explaining domestication and selection for genetic traits of commercial interest.

Making and shaping endochondral and intramembranous bones

Skeletal elements have a diverse range of shapes and sizes specialized to their various roles including protecting internal organs, locomotion, feeding, hearing, and vocalization. The precise positioning, size, and shape of skeletal elements is therefore critical for their function. During embryonic development, bone forms by endochondral or intramembranous ossification and can arise from the paraxial and lateral plate mesoderm or neural crest. This review describes inductive mechanisms to position and pattern bones within the developing embryo, compares and contrasts the intrinsic vs extrinsic mechanisms of endochondral and intramembranous skeletal development, and details known cellular processes that precisely determine skeletal shape and size. Key cellular mechanisms are employed at distinct stages of ossification, many of which occur in response to mechanical cues (eg, joint formation) or preempting future load‐bearing requirements. Rapid shape changes occur during cellular condensation and template establishment. Specialized cellular behaviors, such as chondrocyte hypertrophy in endochondral bone and secondary cartilage on intramembranous bones, also dramatically change template shape. Once ossification is complete, bone shape undergoes functional adaptation through (re)modeling. We also highlight how alterations in these cellular processes contribute to evolutionary change and how differences in the embryonic origin of bones can influence postnatal bone repair.

Compares and contrasts Endochondral and intramembranous bone development

Reviews embryonic origins of different bones

Describes the cellular and molecular mechanisms of positioning skeletal elements.

Describes mechanisms of skeletal growth with a focus on the generation of skeletal shape

Zebrafish models of skeletal dysplasia induced by cholesterol biosynthesis deficiency

Human disorders of the post-squalene cholesterol biosynthesis pathway frequently result in skeletal abnormalities, yet our understanding of the mechanisms involved is limited. In a forward-genetic approach, we have found that a late-onset skeletal mutant, named kolibernu7 , is the result of a cis-acting regulatory mutation leading to loss of methylsterol monooxygenase 1 (msmo1) expression within pre-hypertrophic chondrocytes. Generated msmo1nu81 knockdown mutation resulted in lethality at larval stage. We demonstrated that this is a result of both cholesterol deprivation and sterol intermediate accumulation by creating a mutation eliminating activity of Lanosterol synthase (Lss). Our results indicate that double lssnu60;msmo1nu81 and single lssnu60 mutants survive significantly longer than msmo1nu81 homozygotes. Liver-specific restoration of either Msmo1 or Lss in corresponding mutant backgrounds suppresses larval lethality. Rescued mutants develop dramatic skeletal abnormalities, with a loss of Msmo1 activity resulting in a more-severe patterning defect of a near-complete loss of hypertrophic chondrocytes marked by col10a1a expression. Our analysis suggests that hypertrophic chondrocytes depend on endogenous cholesterol synthesis, and blocking C4 demethylation exacerbates the cholesterol deficiency phenotype. Our findings offer new insight into the genetic control of bone development and provide new zebrafish models for human disorders of the cholesterol biosynthesis pathway.

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.

Zebrafish: A Resourceful Vertebrate Model to Investigate Skeletal Disorders

Animal models are essential tools for addressing fundamental scientific questions about skeletal diseases and for the development of new therapeutic approaches. Traditionally, mice have been the most common model organism in biomedical research, but their use is hampered by several limitations including complex generation, demanding investigation of early developmental stages, regulatory restrictions on breeding, and high maintenance cost. The zebrafish has been used as an efficient alternative vertebrate model for the study of human skeletal diseases, thanks to its easy genetic manipulation, high fecundity, external fertilization, transparency of rapidly developing embryos, and low maintenance cost. Furthermore, zebrafish share similar skeletal cells and ossification types with mammals. In the last decades, the use of both forward and new reverse genetics techniques has resulted in the generation of many mutant lines carrying skeletal phenotypes associated with human diseases. In addition, transgenic lines expressing fluorescent proteins under bone cell- or pathway- specific promoters enable in vivo imaging of differentiation and signaling at the cellular level. Despite the small size of the zebrafish, many traditional techniques for skeletal phenotyping, such as x-ray and microCT imaging and histological approaches, can be applied using the appropriate equipment and custom protocols. The ability of adult zebrafish to remodel skeletal tissues can be exploited as a unique tool to investigate bone formation and repair. Finally, the permeability of embryos to chemicals dissolved in water, together with the availability of large numbers of small-sized animals makes zebrafish a perfect model for high-throughput bone anabolic drug screening. This review aims to discuss the techniques that make zebrafish a powerful model to investigate the molecular and physiological basis of skeletal disorders.

Warfarin-exposed zebrafish embryos resembles human warfarin embryopathy in a dose and developmental-time dependent manner - From molecular mechanisms to environmental concerns

Warfarin is the most worldwide used anticoagulant drug and rodenticide. Since it crosses placental barrier it can induce warfarin embryopathy (WE), a fetal mortality in neonates characterized by skeletal deformities in addition to brain hemorrhages. Although the effects of warfarin exposure in aquatic off target species were already described, the particular molecular toxicological mechanisms during early development are still unclear. Here, we used zebrafish (Danio rerio) to describe and compare the developmental effects of warfarin exposure (0, 15.13, 75.68 and 378.43 mM) on two distinct early developmental phases (embryos and eleuthero-embryos). Although exposure to both developmental phases induced fish mortality, only embryos exposed to the highest warfarin level exhibited features mimicking mammalian WE, e.g. high mortality, higher incidence of hemorrhages and altered skeletal development, among other effects. To gain insights into the toxic mechanisms underlying warfarin exposure, the transcriptome of embryos exposed to warfarin was explored through RNA-Seq and compared to that of control embryos. 766 differentially expressed (564 up- and 202 down-regulated) genes were identified. Gene Ontology analysis revealed particular cellular components (cytoplasm, extracellular matrix, lysosome and vacuole), biological processes (mainly amino acid and lipid metabolism and response to stimulus) and pathways (oxidative stress response and apoptosis signaling pathways) being significantly overrepresented in zebrafish embryos upon warfarin exposure. Protein-protein interaction further evidenced an altered redox system, blood coagulation and vasculogenesis, visual phototransduction and collagen formation upon warfarin exposure. The present study not only describes for the first time the WE in zebrafish, it provides new insights for a better risk assessment, and highlights the need for programming the rat eradication actions outside the fish spawning season to avoid an impact on off target fish community. The urge for the development of more species-specific anticoagulants for rodent pest control is also highlighted.

Analysis of transgenic zebrafish expressing the Lenz-Majewski syndrome gene PTDSS1 in skeletal cell lineages

Background: Lenz-Majewski syndrome (LMS) is characterized by osteosclerosis and hyperostosis of skull, vertebrae and tubular bones as well as craniofacial, dental, cutaneous, and digit abnormalities. We previously found that LMS is caused by de novo dominant missense mutations in the  PTDSS1 gene, which encodes phosphatidylserine synthase 1 (PSS1), an enzyme that catalyses the conversion of phosphatidylcholine to phosphatidylserine. The mutations causing LMS result in a gain-of-function, leading to increased enzyme activity and blocking end-product inhibition of PSS1. Methods: Here, we have used transpose-mediated transgenesis to attempt to stably express wild-type and mutant forms of human PTDSS1 ubiquitously or specifically in chondrocytes, osteoblasts or osteoclasts in zebrafish. Results: We report multiple genomic integration sites for each of 8 different transgenes. While we confirmed that the ubiquitously driven transgene constructs were functional in terms of driving gene expression following transient transfection in HeLa cells, and that all lines exhibited expression of a heart-specific cistron within the transgene, we failed to detect PTDSS1 gene expression at either the RNA or protein levels in zebrafish. All wild-type and mutant transgenic lines of zebrafish exhibited mild scoliosis with variable incomplete penetrance which was never observed in non-transgenic animals. Conclusions: Collectively the data suggest that the transgenes are silenced, that animals with integrations that escape silencing are not viable, or that other technical factors prevent transgene expression. In conclusion, the incomplete penetrance of the phenotype and the lack of a matched transgenic control model precludes further meaningful investigations of these transgenic lines.

Programmed conversion of hypertrophic chondrocytes into osteoblasts and marrow adipocytes within zebrafish bones

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.

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.

A novel nonosteocytic regulatory mechanism of bone modeling

Osteocytes, cells forming an elaborate network within the bones of most vertebrate taxa, are thought to be the master regulators of bone modeling, a process of coordinated, local bone-tissue deposition and removal that keeps bone strains at safe levels throughout life. Neoteleost fish, however, lack osteocytes and yet are known to be capable of bone modeling, although no osteocyte-independent modeling regulatory mechanism has so far been described. Here, we characterize a novel, to our knowledge, bone-modeling regulatory mechanism in a fish species (medaka), showing that although lacking osteocytes (i.e., internal mechanosensors), when loaded, medaka bones model in mechanically directed ways, successfully reducing high tissue strains. We establish that as in mammals, modeling in medaka is regulated by the SOST gene, demonstrating a mechanistic link between skeletal loading, SOST down-regulation, and intense bone deposition. However, whereas mammalian SOST is expressed almost exclusively by osteocytes, in both medaka and zebrafish (a species with osteocytic bones), SOST is expressed by a variety of nonosteocytic cells, none of which reside within the bone bulk. These findings argue that in fishes (and perhaps other vertebrates), nonosteocytic skeletal cells are both sensors and responders, shouldering duties believed exclusive to osteocytes. This previously unrecognized, SOST-dependent, osteocyte-independent mechanism challenges current paradigms of osteocyte exclusivity in bone-modeling regulation, suggesting the existence of multivariate feedback networks in bone modeling—perhaps also in mammalian bones—and thus arguing for the possibility of untapped potential for cell targets in bone therapeutics.

Bone’s ability to change its morphology in response to load is widely attributed to osteocytes. A study of fish shows that bone can respond to load even in the absence of osteocytes, using a molecular mechanism that is conserved across vertebrates, albeit with different cellular effectors.

Bone is a “smart” tissue, able to sense loads within its bulk and change its morphology when needed by a process named bone modeling. This process is carried out by bone-depositing cells (osteoblasts) and bone-resorbing cells (osteoclasts) and is regulated by osteocytes—cells that reside in small cavities within the bone tissue. Osteocytes are considered to function as mechanosensors, detecting areas of high loads that require modeling, and master regulators of osteoblasts and osteoclasts. Curiously, evolutionarily advanced fish do not have osteocytes in their bones, although more basal fish and all other bony vertebrates have them. In this paper, we show how the bones of advanced fish can respond to load in a mechanically efficient way despite the absence of osteocytes. We describe the molecular mechanism, which we found to be the same as in all other vertebrates; however, we show that the cellular effectors are different. The protein sclerostin, which is produced by osteocytes in mammals and is a potent suppressor of bone building by osteoblasts, is produced by a variety of nonosteocytic cells in medaka and zebrafish, and nonosteocytic skeletal cells serve as sensors and responders in these species. These results challenge current paradigms of osteocyte exclusivity in the regulation of bone modeling.

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.

The mechanical impact of col11a2 loss on joints; col11a2 mutant zebrafish show changes to joint development and function, which leads to early-onset osteoarthritis

Collagen is the major structural component of cartilage, and mutations in the genes encoding type XI collagen are associated with severe skeletal dysplasias (fibrochondrogenesis and Stickler syndrome) and early-onset osteoarthritis (OA). The impact of the lack of type XI collagen on cell behaviour and mechanical performance during skeleton development is unknown. We studied a zebrafish mutant for col11a2 and evaluated cartilage, bone development and mechanical properties to address this. We show that in col11a2 mutants, type II collagen is made but is prematurely degraded in maturing cartilage and ectopically expressed in the joint. These changes are correlated with increased stiffness of both bone and cartilage; quantified using atomic force microscopy. In the mutants, the skeletal rudiment terminal region in the jaw joint is broader and the interzone smaller. These differences in shape and material properties impact on joint function and mechanical performance, which we modelled using finite element analyses. Finally, we show that col11a2 heterozygous carriers reach adulthood but show signs of severe early-onset OA. Taken together, our data demonstrate a key role for type XI collagen in maintaining the properties of cartilage matrix; which when lost leads to alterations to cell behaviour that give rise to joint pathologies.This article is part of the Theo Murphy meeting issue 'Mechanics of development'.

Skeletal stiffening in an amphibious fish out of water is a response to increased body weight

Terrestrial animals must support their bodies against gravity, while aquatic animals are effectively weightless because of buoyant support from water. Given this evolutionary history of minimal gravitational loading of fishes in water, it has been hypothesized that weight-responsive musculoskeletal systems evolved during the tetrapod invasion of land and are thus absent in fishes. Amphibious fishes, however, experience increased effective weight when out of water - are these fishes responsive to gravitational loading? Contrary to the tetrapod-origin hypothesis, we found that terrestrial acclimation reversibly increased gill arch stiffness (∼60% increase) in the amphibious fish Kryptolebias marmoratus when loaded normally by gravity, but not under simulated microgravity. Quantitative proteomics analysis revealed that this change in mechanical properties occurred via increased abundance of proteins responsible for bone mineralization in other fishes as well as in tetrapods. Type X collagen, associated with endochondral bone growth, increased in abundance almost ninefold after terrestrial acclimation. Collagen isoforms known to promote extracellular matrix cross-linking and cause tissue stiffening, such as types IX and XII collagen, also increased in abundance. Finally, more densely packed collagen fibrils in both gill arches and filaments were observed microscopically in terrestrially acclimated fish. Our results demonstrate that the mechanical properties of the fish musculoskeletal system can be fine-tuned in response to changes in effective body weight using biochemical pathways similar to those in mammals, suggesting that weight sensing is an ancestral vertebrate trait rather than a tetrapod innovation.

Calcified cartilage or bone? Collagens in the tessellated endoskeletons of cartilaginous fish (sharks and rays)

The primary skeletal tissue in elasmobranchs -sharks, rays and relatives- is cartilage, forming both embryonic and adult endoskeletons. Only the skeletal surface calcifies, exhibiting mineralized tiles (tesserae) sandwiched between a cartilage core and overlying fibrous perichondrium. These two tissues are based on different collagens (Coll II and I, respectively), fueling a long-standing debate as to whether tesserae are more like calcified cartilage or bone (Coll 1-based) in their matrix composition. We demonstrate that stingray (Urobatis halleri) tesserae are bipartite, having an upper Coll I-based 'cap' that merges into a lower Coll II-based 'body' zone, although tesserae are surrounded by cartilage. We identify a 'supratesseral' unmineralized cartilage layer, between tesserae and perichondrium, distinguished from the cartilage core in containing Coll I and X (a common marker for mammalian mineralization), in addition to Coll II. Chondrocytes within tesserae appear intact and sit in lacunae filled with Coll II-based matrix, suggesting tesserae originate in cartilage, despite comprising a diversity of collagens. Intertesseral joints are also complex in their collagenous composition, being similar to supratesseral cartilage closer to the perichondrium, but containing unidentified fibrils nearer the cartilage core. Our results indicate a unique potential for tessellated cartilage in skeletal biology research, since it lacks features believed diagnostic for vertebrate cartilage mineralization (e.g. hypertrophic and apoptotic chondrocytes), while offering morphologies amenable for investigating the regulation of complex mineralized ultrastructure and tissues patterned on multiple collagens.

Ancestral benzo[a]pyrene exposure affects bone integrity in F3 adult fish (Oryzias latipes)

Benzo[a]pyrene (BaP) at an environmentally relevant concentration (1μg/L) has previously been shown to affect bone development in a transgenerational manner in F3 medaka (Oryzias latipes) larvae (17dph). Here, we provide novel histomorphometric data demonstrating that the impaired bone formation at an early life stage is not recoverable and can result in a persistent transgenerational impairment of bone metabolism in F3 adult fish. A decrease in bone thickness and the occurrence of microcracks in ancestrally BaP-treated adult male fish (F3) were revealed by MicroCt measurement and histopathological analysis. The expression of twenty conserved bone miRNAs were screened in medaka and their relative expression (in the F3 ancestral BaP treatment vs the F3 control fish) were determined by quantitative real-time PCR. Attempt was made to link bone miRNA expression with the potential target bone mRNA expression in medaka. Five functional pairs of mRNA/miRNA were identified (Osx/miR-214, Col2a1b/miR-29b, Runx2/miR-204, Sox9b/miR-199a-3p, APC/miR-27b). Unique knowledge of bone-related miRNA expression in medaka in response to ancestral BaP-exposure in the F3 generation is presented. From the ecological risk assessment perspective, BaP needs to be regarded as a transgenerational skeletal toxicant which exerts a far-reaching impact on fish survival and fitness. Given that the underlying mechanisms of cartilage/bone formation are conserved between medaka and mammals, the results may also shed light on the potential transgenerational effect of BaP on skeletal disorders in mammals/humans.

Distinct requirements of wls, wnt9a, wnt5b and gpc4 in regulating chondrocyte maturation and timing of endochondral ossification

Formation of the mandible requires progressive morphologic change, proliferation, differentiation and organization of chondrocytes preceding osteogenesis. The Wnt signaling pathway is involved in regulating bone development and maintenance. Chondrocytes that are fated to become bone require Wnt to polarize and orientate appropriately to initiate the endochondral ossification program. Although the canonical Wnt signaling has been well studied in the context of bone development, the effects of non-canonical Wnt signaling in regulating the timing of cartilage maturation and subsequent bone formation in shaping ventral craniofacial structure is not fully understood.. Here we examined the role of the non-canonical Wnt signaling pathway (wls, gpc4, wnt5b and wnt9a) in regulating zebrafish Meckel's cartilage maturation to the onset of osteogenic differentiation. We found that disruption of wls resulted in a significant loss of craniofacial bone, whereas lack of gpc4, wnt5b and wnt9a resulted in severely delayed endochondral ossification. This study demonstrates the importance of the non-canonical Wnt pathway in regulating coordinated ventral cartilage morphogenesis and ossification.

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.

The Morphogenesis of Cranial Sutures in Zebrafish

Using morphological, histological, and TEM analyses of the cranium, we provide a detailed description of bone and suture growth in zebrafish. Based on expression patterns and localization, we identified osteoblasts at different degrees of maturation. Our data confirm that, unlike in humans, zebrafish cranial sutures maintain lifelong patency to sustain skull growth. The cranial vault develops in a coordinated manner resulting in a structure that protects the brain. The zebrafish cranial roof parallels that of higher vertebrates and contains five major bones: one pair of frontal bones, one pair of parietal bones, and the supraoccipital bone. Parietal and frontal bones are formed by intramembranous ossification within a layer of mesenchyme positioned between the dermal mesenchyme and meninges surrounding the brain. The supraoccipital bone has an endochondral origin. Cranial bones are separated by connective tissue with a distinctive architecture of osteogenic cells and collagen fibrils. Here we show RNA in situ hybridization for col1a1a, col2a1a, col10a1, bglap/osteocalcin, fgfr1a, fgfr1b, fgfr2, fgfr3, foxq1, twist2, twist3, runx2a, runx2b, sp7/osterix, and spp1/ osteopontin, indicating that the expression of genes involved in suture development in mammals is preserved in zebrafish. We also present methods for examining the cranium and its sutures, which permit the study of the mechanisms involved in suture patency as well as their pathological obliteration. The model we develop has implications for the study of human disorders, including craniosynostosis, which affects 1 in 2,500 live births.

Citalopram and sertraline exposure compromises embryonic bone development

Selective serotonin reuptake inhibitors (SSRIs) are the most commonly prescribed treatments for depression and, as a class of drugs, are among the most used medications in the world. Concern regarding possible effects of SSRI treatment on fetal development has arisen recently as studies have suggested a link between maternal SSRI use and an increase in birth defects such as persistent pulmonary hypertension, seizures and craniosynostosis. Furthermore, SSRI exposure in adults is associated with decreased bone mineral density and increased fracture risk, and serotonin receptors are expressed in human osteoblasts and osteoclasts. To determine possible effects of SSRI exposure on developing bone, we treated both zebrafish, during embryonic development, and human mesenchymal stem cells (MSCs), during differentiation into osteoblasts, with the two most prescribed SSRIs, citalopram and sertraline. SSRI treatment in zebrafish decreased bone mineralization, visualized by alizarin red staining and decreased the expression of mature osteoblast-specific markers during embryogenesis. Furthermore, we showed that this inhibition was not associated with increased apoptosis. In differentiating human MSCs, we observed a decrease in osteoblast activity that was associated with a decrease in expression of the osteoblast-specific genes Runx2, Sparc and Spp1, measured with quantitative real-time PCR (qRT-PCR). Similar to the developing zebrafish, no increase in expression of the apoptotic marker Caspase 3 was observed. Therefore, we propose that SSRIs inhibit bone development by affecting osteoblast maturation during embryonic development and MSC differentiation. These results highlight the need to further investigate the risks of SSRI use during pregnancy in exposing unborn babies to potential skeletal abnormalities.

Analyzing Biological Performance of 3D-Printed, Cell-Impregnated Hybrid Constructs for Cartilage Tissue Engineering

Three-dimensional (3D) bioprinting of hybrid constructs is a promising biofabrication method for cartilage tissue engineering because a synthetic polymer framework and cell-impregnated hydrogel provide structural and biological features of cartilage, respectively. During bioprinting, impregnated cells may be subjected to high temperatures (caused by the adjacent melted polymer) and process-induced mechanical forces, potentially compromising cell function. This study addresses these biofabrication issues, evaluating the heat distribution of printed polycaprolactone (PCL) strands and the rheological property and structural stability of alginate hydrogels at various temperatures and concentrations. The biocompatibility of parameters from these studies was tested by culturing 3D hybrid constructs bioprinted with primary cells from embryonic chick cartilage. During initial two-dimensional culture expansion of these primary cells, two morphologically and molecularly distinct cell populations ("rounded" and "fibroblastic") were isolated. The biological performance of each population was evaluated in 3D hybrid constructs separately. The cell viability, proliferation, and cartilage differentiation were observed at high levels in hybrid constructs of both cell populations, confirming the validity of these 3D bioprinting parameters for effective cartilage tissue engineering. Statistically significant performance variations were observed, however, between the rounded and fibroblastic cell populations. Molecular and morphological data support the notion that such performance differences may be attributed to the relative differentiation state of rounded versus fibroblastic cells (i.e., differentiated chondrocytes vs. chondroprogenitors, respectively), which is a relevant issue for cell-based tissue engineering strategies. Taken together, our study demonstrates that bioprinting 3D hybrid constructs of PCL and cell-impregnated alginate hydrogel is a promising approach for cartilage tissue engineering.

Transcriptome analysis of vertebral bone in the flounder, Paralichthys olivaceus (Teleostei, Pleuronectiformes), using Illumina sequencing

The processes underlying vertebral development in teleosts and tetrapods differ markedly in a variety of ways. At present, the molecular basis of teleost vertebral development and growth is poorly understood. Understanding vertebral development at the molecular level is important for aquaculture to prevent vertebral anomalies that can arise from a variety of factors, including excess vitamin A (all-trans retinol, VA) in the diet. To facilitate studies on teloest vertebral development, we performed transcriptome analysis of four month old flounder, Paralichthys olivaceus, vertebrae using next-generation sequencing. Expression profile obtained demonstrates that some members of the hh, bmp, fgf, wnt gene families, and their receptors, hox, pax, sox, dlx and tbx gene families and ntl, which are known to function in notochord and somite development in embryos, are expressed in the vertebrae. It was also showed that in addition to the retinoic acid receptor (Rar), the vertebrae express alcohol dehydrogenase 1 and retinal dehydrogenase 2 which convert VA to all-trans-retinoic acid (RA). The assembled contigs also included cytochrome p450 family members, which inactivate RA, as well as phosphatidylcholine-retinol O-acetyltransferase, which converts VA to all-trans-retinyl ester, a stock form of VA. These data suggest that in teleost vertebrae, expression of various signals and transcription factors which function in the notochord and somite development is maintained until adult stage, and RA metabolism and signaling are active to regulate transcription of RA-responsible genes, such as hedgehog and hox genes. This is the first transcriptome analysis of teleost fish vertebrae.

Danio rerio: the Janus of the bone from embryo to scale

Danio rerio (zebrafish), like the Roman god Janus, is an old animal model which is recently emerged and looks to the future with an increasing scientific success. Unlike other traditional animal models, zebrafish represents a versatile way to approach the study of the skeleton. Transparency of the larval stage, genetic manipulability and unique anatomical structures (scales) makes zebrafish a powerful and versatile instrument to investigate the bone tissue in terms of structure and function. Like Janus, zebrafish offers two different faces, or better, two models in one animal: larval and adult stage. The embryo can be used to isolate new genes involved in osteogenesis by large-scale mutagenesis screenings. The behavior of bone cells and genes in osteogenesis can be investigate by using transgenic lines, vital dyes, mutants and traditional molecular biology techniques. The adult zebrafish represents an important resource to study the pathways related to the bone metabolism and turnover. In particular, the properties of the caudal fin allow to study mechanisms of bone regeneration and reparation whereas the elasmoid scale represents an unique tool to investigate the bone metabolism under physiological or pathological conditions.

On the evolutionary relationship between chondrocytes and osteoblasts

Vertebrates are the only animals that produce bone, but the molecular genetic basis for this evolutionary novelty remains obscure. Here, we synthesize information from traditional evolutionary and modern molecular genetic studies in order to generate a working hypothesis on the evolution of the gene regulatory network (GRN) underlying bone formation. Since transcription factors are often core components of GRNs (i.e., kernels), we focus our analyses on Sox9 and Runx2. Our argument centers on three skeletal tissues that comprise the majority of the vertebrate skeleton: immature cartilage, mature cartilage, and bone. Immature cartilage is produced during early stages of cartilage differentiation and can persist into adulthood, whereas mature cartilage undergoes additional stages of differentiation, including hypertrophy and mineralization. Functionally, histologically, and embryologically, these three skeletal tissues are very similar, yet unique, suggesting that one might have evolved from another. Traditional studies of the fossil record, comparative anatomy and embryology demonstrate clearly that immature cartilage evolved before mature cartilage or bone. Modern molecular approaches show that the GRNs regulating differentiation of these three skeletal cell fates are similar, yet unique, just like the functional and histological features of the tissues themselves. Intriguingly, the Sox9 GRN driving cartilage formation appears to be dominant to the Runx2 GRN of bone. Emphasizing an embryological and evolutionary transcriptomic view, we hypothesize that the Runx2 GRN underlying bone formation was co-opted from mature cartilage. We discuss how modern molecular genetic experiments, such as comparative transcriptomics, can test this hypothesis directly, meanwhile permitting levels of constraint and adaptation to be evaluated quantitatively. Therefore, comparative transcriptomics may revolutionize understanding of not only the clade-specific evolution of skeletal cells, but also the generation of evolutionary novelties, providing a modern paradigm for the evolutionary process.

Molecular footprinting of skeletal tissues in the catshark Scyliorhinus canicula and the clawed frog Xenopus tropicalis identifies conserved and derived features of vertebrate calcification

Understanding the evolutionary emergence and subsequent diversification of the vertebrate skeleton requires a comprehensive view of the diverse skeletal cell types found in distinct developmental contexts, tissues, and species. To date, our knowledge of the molecular nature of the shark calcified extracellular matrix, and its relationships with osteichthyan skeletal tissues, remain scarce. Here, based on specific combinations of expression patterns of the Col1a1, Col1a2, and Col2a1 fibrillar collagen genes, we compare the molecular footprint of endoskeletal elements from the chondrichthyan Scyliorhinus canicula and the tetrapod Xenopus tropicalis. We find that, depending on the anatomical location, Scyliorhinus skeletal calcification is associated to cell types expressing different subsets of fibrillar collagen genes, such as high levels of Col1a1 and Col1a2 in the neural arches, high levels of Col2a1 in the tesserae, or associated to a drastic Col2a1 downregulation in the centrum. We detect low Col2a1 levels in Xenopus osteoblasts, thereby revealing that the osteoblastic expression of this gene was significantly reduced in the tetrapod lineage. Finally, we uncover a striking parallel, from a molecular and histological perspective, between the vertebral cartilage calcification of both species and discuss the evolutionary origin of endochondral ossification.

Structural divergence in vertebrate phylogeny of a duplicated prototype galectin

Prototype galectins, endogenously expressed animal lectins with a single carbohydrate recognition domain, are well-known regulators of tissue properties such as growth and adhesion. The earliest discovered and best studied of the prototype galectins is Galectin-1 (Gal-1). In the Gallus gallus (chicken) genome, Gal-1 is represented by two homologs: Gal-1A and Gal-1B, with distinct biochemical properties, tissue expression, and developmental functions. We investigated the origin of the Gal-1A/Gal-1B divergence to gain insight into when their developmental functions originated and how they could have contributed to vertebrate phenotypic evolution. Sequence alignment and phylogenetic tree construction showed that the Gal-1A/Gal-1B divergence can be traced back to the origin of the sauropsid lineage (consisting of extinct and extant reptiles and birds) although lineage-specific duplications also occurred in the amphibian and actinopterygian genomes. Gene synteny analysis showed that sauropsid gal-1b (the gene for Gal-1B) and its frog and actinopterygian gal-1 homologs share a similar chromosomal location, whereas sauropsid gal-1a has translocated to a new position. Surprisingly, we found that chicken Gal-1A, encoded by the translocated gal-1a, was more similar in its tertiary folding pattern than Gal-1B, encoded by the untranslocated gal-1b, to experimentally determined and predicted folds of nonsauropsid Gal-1s. This inference is consistent with our finding of a lower proportion of conserved residues in sauropsid Gal-1Bs, and evidence for positive selection of sauropsid gal-1b, but not gal-1a genes. We propose that the duplication and structural divergence of Gal-1B away from Gal-1A led to specialization in both expression and function in the sauropsid lineage.

A new model army: Emerging fish models to study the genomics of vertebrate Evo-Devo

Many fields of biology--including vertebrate Evo-Devo research--are facing an explosion of genomic and transcriptomic sequence information and a multitude of fish species are now swimming in this "genomic tsunami." Here, we first give an overview of recent developments in sequencing fish genomes and transcriptomes that identify properties of fish genomes requiring particular attention and propose strategies to overcome common challenges in fish genomics. We suggest that the generation of chromosome-level genome assemblies--for which we introduce the term "chromonome"--should be a key component of genomic investigations in fish because they enable large-scale conserved synteny analyses that inform orthology detection, a process critical for connectivity of genomes. Orthology calls in vertebrates, especially in teleost fish, are complicated by divergent evolution of gene repertoires and functions following two rounds of genome duplication in the ancestor of vertebrates and a third round at the base of teleost fish. Second, using examples of spotted gar, basal teleosts, zebrafish-related cyprinids, cavefish, livebearers, icefish, and lobefin fish, we illustrate how next generation sequencing technologies liberate emerging fish systems from genomic ignorance and transform them into a new model army to answer longstanding questions on the genomic and developmental basis of their biodiversity. Finally, we discuss recent progress in the genetic toolbox for the major fish models for functional analysis, zebrafish, and medaka, that can be transferred to many other fish species to study in vivo the functional effect of evolutionary genomic change as Evo-Devo research enters the postgenomic era.

An RNA-seq protocol to identify mRNA expression changes in mouse diaphyseal bone: applications in mice with bone property altering Lrp5 mutations

Loss-of-function and certain missense mutations in the Wnt coreceptor low-density lipoprotein receptor-related protein 5 (LRP5) significantly decrease or increase bone mass, respectively. These human skeletal phenotypes have been recapitulated in mice harboring Lrp5 knockout and knock-in mutations. We hypothesized that measuring mRNA expression in diaphyseal bone from mice with Lrp5 wild-type (Lrp5(+/+) ), knockout (Lrp5(-/-) ), and high bone mass (HBM)-causing (Lrp5(p.A214V/+) ) knock-in alleles could identify genes and pathways that regulate or are regulated by LRP5 activity. We performed RNA-seq on pairs of tibial diaphyseal bones from four 16-week-old mice with each of the aforementioned genotypes. We then evaluated different methods for controlling for contaminating nonskeletal tissue (ie, blood, bone marrow, and skeletal muscle) in our data. These methods included predigestion of diaphyseal bone with collagenase and separate transcriptional profiling of blood, skeletal muscle, and bone marrow. We found that collagenase digestion reduced contamination, but also altered gene expression in the remaining cells. In contrast, in silico filtering of the diaphyseal bone RNA-seq data for highly expressed blood, skeletal muscle, and bone marrow transcripts significantly increased the correlation between RNA-seq data from an animal's right and left tibias and from animals with the same Lrp5 genotype. We conclude that reliable and reproducible RNA-seq data can be obtained from mouse diaphyseal bone and that lack of LRP5 has a more pronounced effect on gene expression than the HBM-causing LRP5 missense mutation. We identified 84 differentially expressed protein-coding transcripts between LRP5 "sufficient" (ie, Lrp5(+/+) and Lrp5(p.A214V/+) ) and "insufficient" (Lrp5(-/-) ) diaphyseal bone, and far fewer differentially expressed genes between Lrp5(p.A214V/+) and Lrp5(+/+) diaphyseal bone.

Comparative analysis of a teleost skeleton transcriptome provides insight into its regulation

An articulated endoskeleton that is calcified is a unifying innovation of the vertebrates, however the molecular basis of the structural divergence between terrestrial and aquatic vertebrates, such as teleost fish, has not been determined. In the present study long-read next generation sequencing (NGS, Roche 454 platform) was used to characterize acellular perichondral bone (vertebrae) and chondroid bone (gill arch) in the gilthead sea bream (Sparus auratus). A total of 15.97 and 14.53Mb were produced, respectively from vertebrae and gill arch cDNA libraries and yielded 32,374 and 28,371 contigs (consensus sequences) respectively. 10,455 contigs from vertebrae and 10,625 contigs from gill arches were annotated with gene ontology terms. Comparative analysis of the global transcriptome revealed 4249 unique transcripts in vertebrae, 4201 unique transcripts in the gill arches and 3700 common transcripts. Several core gene networks were conserved between the gilthead sea bream and mammalian skeleton. Transcripts for putative endocrine factors were identified in acellular gilthead sea bream bone suggesting that in common with mammalian bone it can act as an endocrine tissue. The acellular bone of the vertebra, in contrast to current opinion based on histological analysis, was responsive to a short fast and significant (p<0.05) down-regulation of several transcripts identified by NGS, osteonectin, osteocalcin, cathepsin K and IGFI occurred. In gill arches fasting caused a significant (p<0.05) down-regulation of osteocalcin and up-regulation of MMP9.

Establishment of a bone-specific col10a1:GFP transgenic zebrafish

During skeletal development, both osteogenic and chondrogenic programs are initiated from multipotent mesenchymal cells, requiring a number of signaling molecules, transcription factors, and downstream effectors to orchestrate the sophisticated process. Col10a1, an important downstream effector gene, has been identified as a marker for maturing chondrocytes in higher vertebrates, such as mammals and birds. In zebrafish, this gene has been shown to be expressed in both osteoblasts and chondrocytes, but no study has reported its role in osteoblast development. To initially delineate the osteogenic program from chondrogenic lineage development, we used the zebrafish col10a1 promoter to establish a transgenic zebrafish expressing a GFP reporter specifically in osteoblast-specific bone structures that do not involve cartilaginous programs. A construct harboring a -2.2-kb promoter region was found to be sufficient to drive the reporter gene in osteoblast-specific bone structures within the endogenous col10a1 expression domain, confirming that separable cis-acting elements exist for distinct cell type-specific expression of col10a1 during zebrafish skeletal development. The -2.2-kb col10a1:GFP transgenic zebrafish marking only bone structures derived from osteoblasts will undoubtedly be an invaluable tool for identifying and characterizing molecular events driving osteoblast development in zebrafish, which may further provide a differential mechanism where col10a1 is involved in the development of chondrocytes undergoing maturation in other vertebrate systems.

Evolution of the vertebrate bone matrix: an expression analysis of the network forming collagen paralogues in amphibian osteoblasts

The emergence of vertebrates is closely associated to the evolution of mineralized bone tissue. However, the molecular basis underlying the origin and subsequent diversification of the skeletal mineralized matrix is still poorly understood. One efficient way to tackle this issue is to compare the expression, between vertebrate species, of osteoblastic genes coding for bone matrix proteins. In this work, we have focused on the evolution of the network forming collagen family which contains the Col8a1, Col8a2, and Col10a1 genes. Both phylogeny and synteny reveal that these three paralogues are vertebrate-specific and derive from two independent duplications in the vertebrate lineage. To shed light on the evolution of this family, we have analyzed the osteoblastic expression of the network forming collagens in endochondral and intramembraneous skeletal elements of the amphibian Xenopus tropicalis. Remarkably, we find that amphibian osteoblasts express Col10a1, a gene strongly expressed in osteoblasts in actinopterygians but not in amniotes. In addition, while Col8a1 is known to be robustly expressed in mammalian osteoblasts, the expression levels of its amphibian orthologue are dramatically reduced. Our work reveals that while a skeletal expression of network forming collagen members is widespread throughout vertebrates, osteoblasts from divergent vertebrate lineages express different combinations of network forming collagen paralogues.