Compressive stress gradients direct mechanoregulation of anisotropic growth in the zebrafish jaw joint.

Mechanical stimuli arising from fetal movements are critical factors underlying joint growth. Abnormal fetal movements negatively affect joint shape features with important implications for joint health, but the mechanisms by which mechanical forces from fetal movements influence joint growth are still unclear. In this research, we quantify zebrafish jaw joint growth in 3D in free-to-move and immobilised fish larvae between four and five days post fertilisation. We found that the main changes in size and shape in normally moving fish were in the ventrodorsal axis, while growth anisotropy was lost in the immobilised larvae. We next sought to determine the cell level activities underlying mechanoregulated growth anisotropy by tracking individual cells in the presence or absence of jaw movements, finding that the most dramatic changes in growth rates due to jaw immobility were in the ventrodorsal axis. Finally, we implemented mechanobiological simulations of joint growth with which we tested hypotheses relating specific mechanical stimuli to mechanoregulated growth anisotropy. Different types of mechanical stimulation were incorporated into the simulation to provide the mechanoregulated component of growth, in addition to the baseline (non-mechanoregulated) growth which occurs in the immobilised animals. We found that when average tissue stress over the opening and closing cycle of the joint was used as the stimulus for mechanoregulated growth, joint morphogenesis was not accurately predicted. Predictions were improved when using the stress gradients along the rudiment axes (i.e., the variation in magnitude of compression to magnitude of tension between local regions). However, the most accurate predictions were obtained when using the compressive stress gradients (i.e., the variation in compressive stress magnitude) along the rudiment axes. We conclude therefore that the dominant biophysical stimulus contributing to growth anisotropy during early joint development is the gradient of compressive stress experienced along the growth axes under cyclical loading.

NCOA3 identified as a new candidate to explain autosomal dominant progressive hearing loss.

Hearing loss is a frequent sensory impairment in humans and genetic factors account for an elevated fraction of the cases. We have investigated a large family of five generations, with 15 reported individuals presenting non-syndromic, sensorineural, bilateral and progressive hearing loss, segregating as an autosomal dominant condition. Linkage analysis, using SNP-array and selected microsatellites, identified a region of near 13 cM in chromosome 20 as the best candidate to harbour the causative mutation. After exome sequencing and filtering of variants, only one predicted deleterious variant in the NCOA3 gene (NM_181659, c.2810C > G; p.Ser937Cys) fit in with our linkage data. RT-PCR, immunostaining and in situ hybridization showed expression of ncoa3 in the inner ear of mice and zebrafish. We generated a stable homozygous zebrafish mutant line using the CRISPR/Cas9 system. ncoa3-/- did not display any major morphological abnormalities in the ear, however, anterior macular hair cells showed altered orientation. Surprisingly, chondrocytes forming the ear cartilage showed abnormal behaviour in ncoa3-/-, detaching from their location, invading the ear canal and blocking the cristae. Adult mutants displayed accumulation of denser material wrapping the otoliths of ncoa3-/- and increased bone mineral density. Altered zebrafish swimming behaviour corroborates a potential role of ncoa3 in hearing loss. In conclusion, we identified a potential candidate gene to explain hereditary hearing loss, and our functional analyses suggest subtle and abnormal skeletal behaviour as mechanisms involved in the pathogenesis of progressive sensory function impairment.

Dominating lengthscales of zebrafish collective behaviour.

Collective behaviour in living systems is observed across many scales, from bacteria to insects, to fish shoals. Zebrafish have emerged as a model system amenable to laboratory study. Here we report a three-dimensional study of the collective dynamics of fifty zebrafish. We observed the emergence of collective behaviour changing between ordered to randomised, upon adaptation to new environmental conditions. We quantify the spatial and temporal correlation functions of the fish and identify two length scales, the persistence length and the nearest neighbour distance, that capture the essence of the behavioural changes. The ratio of the two length scales correlates robustly with the polarisation of collective motion that we explain with a reductionist model of self-propelled particles with alignment interactions.

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-3). The DEGs were also enriched for human orthologues associated with polygenetic skeletal traits, including height (P< 6× 10-4), and estimated bone mineral density (eBMD, P< 2× 10-5). Zebrafish mutants of two human orthologues that were robustly associated with height (COL11A2, P=6× 10-24) or eBMD (SPP1, P=6× 10-20) showed both exo- and endo- skeletal abnormalities as predicted by our genetic association analyses; col11a2Y228X/Y228X mutants showed exoskeletal and endoskeletal features consistent with abnormal growth, whereas spp1P160X/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.

Growth orientations, rather than heterogeneous growth rates, dominate jaw joint morphogenesis in the larval zebrafish.

In early limb embryogenesis, synovial joints acquire specific shapes which determine joint motion and function. The process by which the opposing cartilaginous joint surfaces are moulded into reciprocal and interlocking shapes, called joint morphogenesis, is one of the least understood aspects of joint formation and the cell-level dynamics underlying it are yet to be unravelled. In this research, we quantified key cellular dynamics involved in growth and morphogenesis of the zebrafish jaw joint and synthesised them in a predictive computational simulation of joint development. Cells in larval zebrafish jaw joints labelled with cartilage markers were tracked over a 48-h time window using confocal imaging. Changes in distance and angle between adjacent cell centroids resulting from cell rearrangement, volume expansion and extracellular matrix (ECM) deposition were measured and used to calculate the rate and direction of local tissue deformations. We observed spatially and temporally heterogeneous growth patterns with marked anisotropy over the developmental period assessed. There was notably elevated growth at the level of the retroarticular process of the Meckel's cartilage, a feature known to undergo pronounced shape changes during zebrafish development. Analysis of cell dynamics indicated a dominant role for cell volume expansion in growth, with minor influences from ECM volume increases and cell intercalation. Cell proliferation in the joint was minimal over the timeframe of interest. Synthesising the dynamic cell data into a finite element model of jaw joint development resulted in accurate shape predictions. Our biofidelic computational simulation demonstrated that zebrafish jaw joint growth can be reasonably approximated based on cell positional information over time, where cell positional information derives mainly from cell orientation and cell volume expansion. By modifying the input parameters of the simulation, we were able to assess the relative contributions of heterogeneous growth rates and of growth orientation. The use of uniform rather than heterogeneous growth rates only minorly impacted the shape predictions, whereas isotropic growth fields resulted in altered shape predictions. The simulation results suggest that growth anisotropy is the dominant influence on joint growth and morphogenesis. This study addresses the gap of the cellular processes underlying joint morphogenesis, with implications for understanding the aetiology of developmental joint disorders such as developmental dysplasia of the hip and arthrogryposis.

Autophagy coordinates chondrocyte development and early joint formation in zebrafish.

Autophagy is a catabolic process responsible for the removal of waste and damaged cellular components by lysosomal degradation. It plays a key role in fundamental cell processes, including ER stress mitigation, control of cell metabolism, and cell differentiation and proliferation, all of which are essential for cartilage cell (chondrocyte) development and survival, and for the formation of cartilage. Correspondingly, autophagy dysregulation has been implicated in several skeletal disorders such as osteoarthritis and osteoporosis. To test the requirement for autophagy during skeletal development in zebrafish, we generated an atg13 CRISPR knockout zebrafish line. This line showed a complete loss of atg13 expression, and restricted autophagic activity in vivo. In the absence of autophagy, chondrocyte maturation was accelerated, with chondrocytes exhibiting signs of premature hypertrophy. Focussing on the jaw element, autophagy disruption affected joint articulation causing restricted mouth opening. This gross behavioural phenotype corresponded with a failure to thrive, and death in homozygote atg13 nulls within 17 days. Taken together, our results are consistent with autophagy contributing to the timely regulation of chondrocyte maturation and for extracellular matrix formation.

Regulator of calcineurin-2 is a centriolar protein with a role in cilia length control.

Almost every cell in the human body extends a primary cilium. Defective cilia function leads to a set of disorders known as ciliopathies, which are characterised by debilitating developmental defects that affect many tissues. Here, we report a new role for regulator of calcineurin 2 (RCAN2) in primary cilia function. It localises to centrioles and the basal body and is required to maintain normal cilia length. RCAN2 was identified as the most strongly upregulated gene from a comparative RNAseq analysis of cells in which expression of the Golgi matrix protein giantin had been abolished by gene editing. In contrast to previous work where we showed that depletion of giantin by RNAi results in defects in ciliogenesis and in cilia length control, giantin knockout cells generate normal cilia after serum withdrawal. Furthermore, giantin knockout zebrafish show increased expression of RCAN2. Importantly, suppression of RCAN2 expression in giantin knockout cells results in the same defects in the control of cilia length that are seen upon RNAi of giantin itself. Together, these data define RCAN2 as a regulator of cilia function that can compensate for the loss of giantin function.

Wnt signalling controls the response to mechanical loading during zebrafish joint development.

Joint morphogenesis requires mechanical activity during development. Loss of mechanical strain causes abnormal joint development, which can impact long-term joint health. Although cell orientation and proliferation are known to shape the joint, dynamic imaging of developing joints in vivo has not been possible in other species. Using genetic labelling techniques in zebrafish we were able, for the first time, to dynamically track cell behaviours in intact moving joints. We identify that proliferation and migration, which contribute to joint morphogenesis, are mechanically controlled and are significantly reduced in immobilised larvae. By comparison with strain maps of the developing skeleton, we identify canonical Wnt signalling as a candidate for transducing mechanical forces into joint cell behaviours. We show that, in the jaw, Wnt signalling is reduced specifically in regions of high strain in response to loss of muscle activity. By pharmacological manipulation of canonical Wnt signalling, we demonstrate that Wnt acts downstream of mechanical activity and is required for joint patterning and chondrocyte maturation. Wnt16, which is also downstream of muscle activity, controls proliferation and migration, but plays no role in chondrocyte intercalation.

Potential of zebrafish as a model to characterise MicroRNA profiles in mechanically mediated joint degeneration.

Mechanically mediated joint degeneration and cartilage dyshomeostasis is implicated in highly prevalent diseases such as osteoarthritis. Increasingly, MicroRNAs are being associated with maintaining the normal state of cartilage, making them an exciting and potentially key contributor to joint health and disease onset. Here, we present a summary of current in vitro and in vivo models which can be used to study the role of mechanical load and MicroRNAs in joint degeneration, including: non-invasive murine models of PTOA, surgical models which involve ligament transection, and unloading models based around immobilisation of joints or removal of load from the joint through suspension. We also discuss how zebrafish could be used to advance this field, namely through the availability of transgenic lines relevant to cartilage homeostasis and the ability to accurately map strain through the cartilage, enabling the response of downstream MicroRNA targets to be followed dynamically at a cellular level in areas of high and low strain.

Zebrafish as a model to study autophagy and its role in skeletal development and disease.

In the last twenty years, research using zebrafish as a model organism has increased immensely. With the many advantages that zebrafish offer such as high fecundity, optical transparency, ex vivo development, and genetic tractability, they are well suited to studying developmental processes and the effect of genetic mutations. More recently, zebrafish models have been used to study autophagy. This important protein degradation pathway is needed for cell and tissue homeostasis in a variety of contexts. Correspondingly, its dysregulation has been implicated in multiple diseases including skeletal disorders. In this review, we explore how zebrafish are being used to study autophagy in the context of skeletal development and disease, and the ways these areas are intersecting to help identify potential therapeutic targets for skeletal disorders.

Scleraxis genes are required for normal musculoskeletal development and for rib growth and mineralization in zebrafish.

Tendons are an essential part of the musculoskeletal system, connecting muscle and skeletal elements to enable force generation. The transcription factor scleraxis marks vertebrate tendons from early specification. Scleraxis-null mice are viable and have a range of tendon and bone defects in the trunk and limbs but no described cranial phenotype. We report the expression of zebrafish scleraxis orthologs: scleraxis homolog (scx)-a and scxb in cranial and intramuscular tendons and in other skeletal elements. Single mutants for either scxa or scxb, generated by clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9), are viable and fertile as adult fish. Although scxb mutants show no obvious phenotype, scxa mutant embryos have defects in cranial tendon maturation and muscle misalignment. Mutation of both scleraxis genes results in more severe defects in cranial tendon differentiation, muscle and cartilage dysmorphogenesis and paralysis, and lethality by 2-5 wk, which indicates an essential function of scleraxis for craniofacial development. At juvenile and adult stages, ribs in scxa mutants fail to mineralize and/or are small and heavily fractured. Scxa mutants also have smaller muscle volume, abnormal swim movement, and defects in bone growth and composition. Scleraxis function is therefore essential for normal craniofacial form and function and vital for fish development.-Kague, E., Hughes, S. M., Lawrence, E. A., Cross, S., Martin-Silverstone, E., Hammond, C. L., Hinits, Y. Scleraxis genes are required for normal musculoskeletal development and for rib growth and mineralization in zebrafish.

Giantin-knockout models reveal a feedback loop between Golgi function and glycosyltransferase expression.

The Golgi is the cellular hub for complex glycosylation, controlling accurate processing of complex proteoglycans, receptors, ligands and glycolipids. Its structure and organisation are dependent on golgins, which tether cisternal membranes and incoming transport vesicles. Here, we show that knockout of the largest golgin, giantin, leads to substantial changes in gene expression but only limited effects on Golgi structure. Notably, 22 Golgi-resident glycosyltransferases, but not glycan-processing enzymes or the ER glycosylation machinery, are differentially expressed following giantin ablation. This includes near-complete loss of function of GALNT3 in both mammalian cell and zebrafish models. Giantin-knockout zebrafish exhibit hyperostosis and ectopic calcium deposits, recapitulating phenotypes of hyperphosphatemic familial tumoral calcinosis, a disease caused by mutations in GALNT3. These data reveal a new feature of Golgi homeostasis: the ability to regulate glycosyltransferase expression to generate a functional proteoglycome.

Glucocerebrosidase deficiency in zebrafish affects primary bone ossification through increased oxidative stress and reduced Wnt/ß-catenin signaling.

Loss of lysosomal glucocerebrosidase (GBA1) function is responsible for several organ defects, including skeletal abnormalities in type 1 Gaucher disease (GD). Enhanced bone resorption by infiltrating macrophages has been proposed to lead to major bone defects. However, while more recent evidences support the hypothesis that osteoblastic bone formation is impaired, a clear pathogenetic mechanism has not been depicted yet. Here, by combining different molecular approaches, we show that Gba1 loss of function in zebrafish is associated with defective canonical Wnt signaling, impaired osteoblast differentiation and reduced bone mineralization. We also provide evidence that increased reactive oxygen species production precedes the Wnt signaling impairment, which can be reversed upon human GBA1 overexpression. Type 1 GD patient fibroblasts similarly exhibit reduced Wnt signaling activity, as a consequence of increased ß-catenin degradation. Our results support a novel model in which a primary defect in canonical Wnt signaling antecedes bone defects in type 1 GD.

Entpd5 is essential for skeletal mineralization and regulates phosphate homeostasis in zebrafish.

Bone mineralization is an essential step during the embryonic development of vertebrates, and bone serves vital functions in human physiology. To systematically identify unique gene functions essential for osteogenesis, we performed a forward genetic screen in zebrafish and isolated a mutant, no bone (nob), that does not form any mineralized bone. Positional cloning of nob identified the causative gene to encode ectonucleoside triphosphate/diphosphohydrolase 5 (entpd5); analysis of its expression pattern demonstrates that entpd5 is specifically expressed in osteoblasts. An additional mutant, dragonfish (dgf), exhibits ectopic mineralization in the craniofacial and axial skeleton and encodes a loss-of-function allele of ectonucleotide pyrophosphatase phosphodiesterase 1 (enpp1). Intriguingly, generation of double-mutant nob/dgf embryos restored skeletal mineralization in nob mutants, indicating that mechanistically, Entpd5 and Enpp1 act as reciprocal regulators of phosphate/pyrophosphate homeostasis in vivo. Consistent with this, entpd5 mutant embryos can be rescued by high levels of inorganic phosphate, and phosphate-regulating factors, such as fgf23 and npt2a, are significantly affected in entpd5 mutant embryos. Our study demonstrates that Entpd5 represents a previously unappreciated essential player in phosphate homeostasis and skeletal mineralization.

Expression of glycosaminoglycan epitopes during zebrafish skeletogenesis.

BACKGROUND: The zebrafish is an important developmental model. Surprisingly, there are few studies that describe the glycosaminoglycan composition of its extracellular matrix during skeletogenesis. Glycosaminoglycans on proteoglycans contribute to the material properties of musculo skeletal connective tissues, and are important in regulating signalling events during morphogenesis. Sulfation motifs within the chain structure of glycosaminoglycans on cell-associated and extracellular matrix proteoglycans allow them to bind and regulate the sequestration/presentation of bioactive signalling molecules important in musculo-skeletal development. RESULTS: We describe the spatio-temporal expression of different glycosaminoglycan moieties during zebrafish skeletogenesis with antibodies recognising (1) native sulfation motifs within chondroitin and keratan sulfate chains, and (2) enzyme-generated neoepitope sequences within the chain structure of chondroitin sulfate (i.e., 0-, 4-, and 6-sulfated isoforms) and heparan sulfate glycosaminoglycans. We show that all the glycosaminoglycan moieties investigated are expressed within the developing skeletal tissues of larval zebrafish. However, subtle changes in their patterns of spatio-temporal expression over the period examined suggest that their expression is tightly and dynamically controlled during development. CONCLUSIONS: The subtle differences observed in the domains of expression between different glycosaminoglycan moieties suggest differences in their functional roles during establishment of the primitive analogues of the skeleton.

Zebrafish Scale Regeneration In Toto and Ex Vivo Culture of Scales.

Skeletal diseases are often complex in their etiology and affect millions of people worldwide. Due to the aging population, there is a need for new therapeutics that could ease the burden on healthcare systems. As these diseases are complex, it is difficult and expensive to accurately model bone pathophysiology in a lab setting. The challenge for the field is to establish a cost-effective, biologically relevant platform for modeling bone disease that can be used to test potential therapeutic compounds. Such a platform should ideally allow dynamic visualization of cell behaviors of bone-building osteoblasts and bone-degrading osteoclasts acting in their mineralized matrix environment. Zebrafish are increasingly used as models due to the availability of genetic tools, including transgenic reporter lines, and the fact that some skeletal tissues (including the scales) remain translucent to adulthood, allowing dynamic imaging options. Since zebrafish scales have both osteoblasts and osteoclasts and are highly abundant, they provide an easily accessible and abundantly available resource of independent bone units. Moreover, once removed, adult zebrafish scales fully regenerate, therefore offering a way to study the spatiotemporal growth of mineralized tissue in vivo. Here, we detail protocols for harvesting and tracking the regeneration of the scales. Lastly, a protocol for stable culture of scales ex vivo for a week and following the healing response after controlled damage to the mineralized matrix of the scale over time is also presented.

Osteomodulin downregulation is associated with osteoarthritis development.

Abnormal subchondral bone remodeling leading to sclerosis is a main feature of osteoarthritis (OA), and osteomodulin (OMD), a proteoglycan involved in extracellular matrix mineralization, is associated with the sclerotic phenotype. However, the functions of OMD remain poorly understood, specifically in vivo. We used Omd knockout and overexpressing male mice and mutant zebrafish to study its roles in bone and cartilage metabolism and in the development of OA. The expression of Omd is deeply correlated with bone and cartilage microarchitectures affecting the bone volume and the onset of subchondral bone sclerosis and spontaneous cartilage lesions. Mechanistically, OMD binds to RANKL and inhibits osteoclastogenesis, thus controlling the balance of bone remodeling. In conclusion, OMD is a key factor in subchondral bone sclerosis associated with OA. It participates in bone and cartilage homeostasis by acting on the regulation of osteoclastogenesis. Targeting OMD may be a promising new and personalized approach for OA.

Pharmacological Manipulation of Early Zebrafish Skeletal Development Shows an Important Role for Smad9 in Control of Skeletal Progenitor Populations.

Osteoporosis and other conditions associated with low bone density or quality are highly prevalent, are increasing as the population ages and with increased glucocorticoid use to treat conditions with elevated inflammation. There is an unmet need for therapeutics which can target skeletal precursors to induce osteoblast differentiation and osteogenesis. Genes associated with high bone mass represent interesting targets for manipulation, as they could offer ways to increase bone density. A damaging mutation in SMAD9 has recently been associated with high bone mass. Here we show that Smad9 labels groups of osteochondral precursor cells, which are not labelled by the other Regulatory Smads: Smad1 or Smad5. We show that Smad9+ cells are proliferative, and that the Smad9+ pocket expands following osteoblast ablation which induced osteoblast regeneration. We further show that treatment with retinoic acid, prednisolone, and dorsomorphin all alter Smad9 expression, consistent with the effects of these drugs on the skeletal system. Taken together these results demonstrate that Smad9+ cells represent an undifferentiated osteochondral precursor population, which can be manipulated by commonly used skeletal drugs. We conclude that Smad9 represents a target for future osteoanabolic therapies.

Wnt16 Elicits a Protective Effect Against Fractures and Supports Bone Repair in Zebrafish.

Bone homeostasis is a dynamic, multicellular process that is required throughout life to maintain bone integrity, prevent fracture, and respond to skeletal damage. WNT16 has been linked to bone fragility and osteoporosis in human genome wide-association studies, as well as the functional hematopoiesis of leukocytes in vivo. However, the mechanisms by which WNT16 promotes bone health and repair are not fully understood. In this study, CRISPR-Cas9 was used to generate mutant zebrafish lacking Wnt16 (wnt16 -/- ) to study its effect on bone dynamically. The wnt16 mutants displayed variable tissue mineral density (TMD) and were susceptible to spontaneous fractures and the accumulation of bone calluses at an early age. Fractures were induced in the lepidotrichia of the caudal fins of wnt16 -/- and WT zebrafish; this model was used to probe the mechanisms by which Wnt16 regulates skeletal and immune cell dynamics in vivo. In WT fins, wnt16 expression increased significantly during the early stages for bone repair. Mineralization of bone during fracture repair was significantly delayed in wnt16 mutants compared with WT zebrafish. Surprisingly, there was no evidence that the recruitment of innate immune cells to fractures or soft callus formation was altered in wnt16 mutants. However, osteoblast recruitment was significantly delayed in wnt16 mutants postfracture, coinciding with precocious activation of the canonical Wnt signaling pathway. In situ hybridization suggests that canonical Wnt-responsive cells within fractures are osteoblast progenitors, and that osteoblast differentiation during bone repair is coordinated by the dynamic expression of runx2a and wnt16. This study highlights zebrafish as an emerging model for functionally validating osteoporosis-associated genes and investigating fracture repair dynamically in vivo. Using this model, it was found that Wnt16 protects against fracture and supports bone repair, likely by modulating canonical Wnt activity via runx2a to facilitate osteoblast differentiation and bone matrix deposition. © 2021 The Authors. JBMR Plus published by Wiley Periodicals LLC. on behalf of American Society for Bone and Mineral Research.