Swim-training changes the spatio-temporal dynamics of skeletogenesis in zebrafish larvae (Danio rerio)

A Zebrafish Mutant in the Extracellular Matrix Protein Gene efemp1 as a Model for Spinal Osteoarthritis

Osteoarthritis is a degenerative articular disease affecting mainly aging animals and people. The extracellular matrix protein Efemp1 was previously shown to have higher turn-over and increased secretion in the blood serum, urine, and subchondral bone of knee joints in osteoarthritic patients. Here, we use the zebrafish as a model system to investigate the function of Efemp1 in vertebrate skeletal development and homeostasis. Using in situ hybridization, we show that the efemp1 gene is expressed in the brain, the pharyngeal arches, and in the chordoblasts surrounding the notochord at 48 hours post-fertilization. We generated an efemp1 mutant line, using the CRISPR/Cas9 method, that produces a severely truncated Efemp1 protein. These mutant larvae presented a medially narrower chondrocranium at 5 days, which normalized later at day 10. At age 1.5 years, µCT analysis revealed an increased tissue mineral density and thickness of the vertebral bodies, as well as a decreased distance between individual vertebrae and ruffled borders of the vertebral centra. This novel defect, which has, to our knowledge, never been described before, suggests that the efemp1 mutant represents the first zebrafish model for spinal osteoarthritis.

From zero to ossified: Larval skeletal ontogeny of the Neotropical Cichlid fish Cichlasoma dimerus

The identification of skeletal elements, the analysis of their developmental sequence, and the time of their appearance during larval development are essential to broaden the knowledge of each fish species and to recognize skeletal abnormalities that may affect further fish performance. Therefore, this study aimed to provide a general description of the development of the entire skeleton highlighting its variability in Cichlasoma dimerus. Larvae of C. dimersus were stained with alcian blue and alizarin red from hatching to 25 days posthatching. Skeletogenesis began with the endoskeletal disk and some cartilage structures from the caudal fin and the splachnocranium, while the first bony structure observed was the cleithrum. When larvae reached the free-swimming and exogenous feeding stage, mostly bones from the jaws, the branchial arches, and the opercle series evidenced some degree of ossification, suggesting that the ossification sequence of C. dimerus adjusts to physiological demands such as feeding and ventilation. The caudal region was the most variable regarding meristic counts and evidenced higher incidence of bone deformities. In conclusion, this work provides an overview of C. dimerus skeletogenesis and lays the groundwork for further studies on diverse topics, like developmental plasticity, rearing conditions, or phylogenetic relationships.

Optimizing assays of zebrafish larvae swimming performance for drug discovery

INTRODUCTION

Zebrafish larvae are one of the few vertebrates amenable to large-scale drug discovery screens. Larval swimming behavior is often used as an outcome variable and many fields of study have developed assays for evaluating swimming performance. An unintended consequence of this wide interest is that details related to assay methodology and interpretation become scattered across the literature. The aim of this review is to consolidate this information, particularly as it relates to high-throughput approaches.

AREAS COVERED

The authors describe larval swimming behaviors as this forms the basis for understanding their experimentally evoked swimming or spontaneous activity. Next, they detail how swimming activity can serve as an outcome variable, particularly in the multi-well formats used in large-scale screening studies. They also highlight biological and technical factors that can impact the sensitivity and variability of these measurements.

EXPERT OPINION

Careful attention to animal husbandry, experimental design, data acquisition, and interpretation of results can improve screen outcomes by maximizing swimming activity while minimizing intra- and inter-larval variability. The development of more sensitive, quantitative methods of assessing swimming performance that can be incorporated into high-throughput workflows will be important in order to take full advantage of the zebrafish model.

Integration of Metabolomics and Transcriptomics to Reveal the Metabolic Characteristics of Exercise-Improved Bone Mass

(1) Background: Exercise is effective in promoting and maintaining bone mass. The aim of this study was to detect the exercise-induced metabolic changes in bone tissue of zebrafish. (2) Methods: Thirty-eight zebrafish (Danio rerio, six months old) were analyzed. The exercise group (n = 19) received 8 weeks of counter-current swimming training. The control group (n = 19) was not subjected to exercise. Mineralization was quantified, and alkaline phosphatase (Alp) and anti-tartrate acid phosphatase (Trap) activities were estimated (n = 12). The metabolomics (n = 12) and transcriptomics (n = 14) data of bone tissue were used for the integration analyses. (3) Results: The results showed that the exercise training improved the bone mineralization of zebrafish, e.g., the exercise group (5.74 × 104 ± 7.63 × 103) had a higher mean optical density than the control group (5.26 × 104 ± 8.56 × 103, p = 0.046) for the caudal vertebrae. The amount of mineralized matrix in scales of the exercised zebrafish was also higher (0.156 ± 0.012 vs. 0.102 ± 0.003, p = 0.005). Both histological staining and biochemical analysis revealed increased Alp activity (0.81 ± 0.26 vs. 0.76 ± 0.01, p = 0.002) and decreased Trap activity (1.34 ± 0.01 vs. 1.36 ± 0.01, p = 0.005) in the exercise group. A total of 103 different metabolites (DMs, VIP ≥ 1, fold change (FC) ≥ 1.20 or ≤0.83, p < 0.050) were identified. Alanine, aspartate and glutamate metabolism, β-alanine metabolism, pyrimidine metabolism, and pantothenate and CoA biosynthesis were the significantly enriched metabolic pathways (p < 0.050). A total of 35 genes (q ≤ 0.050 (BH), |Log2FC| ≥ 0.5) were coenriched with the 103 DMs in the four identified pathways. Protein–protein interaction network analysis of the 35 genes showed that entpd3, entpd1, and cmpk2 were the core genes. (4) Conclusions: The results of this study suggest that alanine, aspartate and glutamate metabolism, β-alanine metabolism, pyrimidine metabolism, and pantothenate and CoA biosynthesis contributed to exercise-induced improvements in bone mass.

Juvenile zebrafish (Danio rerio) are able to recover from lordosis

Haemal lordosis, a frequent skeletal deformity in teleost fish, has long been correlated with increased mechanical loads induced by swimming activity. In the present study, we examine whether juvenile zebrafish can recover from haemal lordosis and explore the musculoskeletal mechanisms involved. Juveniles were subjected to a swimming challenge test (SCT) that induced severe haemal lordosis in 49% of the animals and then immediately transferred them to 0.0 total body lengths (TL) per second of water velocity for a week. The recovery from lordosis was examined by means of whole mount staining, histology and gene expression analysis. Results demonstrate that 80% of the lordotic zebrafish are capable of internal and external recovery within a week after the SCT. Recovered individuals presented normal shape of the vertebral centra, maintaining though distorted internal tissue organization. Through the transcriptomic analysis of the affected haemal regions, several processes related to chromosome organization, DNA replication, circadian clock and transcription regulation were enriched within genes significantly regulated behind this musculoskeletal recovery procedure. Genes especially involved in adipogenesis, bone remodeling and muscular regeneration were regulated. A remodeling tissue-repair hypothesis behind haemal lordosis recovery is raised. Limitations and future possibilities for zebrafish as a model organism to clarify mechanically driven musculoskeletal changes are discussed.

A Baseline for Skeletal Investigations in Medaka (Oryzias latipes): The Effects of Rearing Density on the Postcranial Phenotype

Oryzias latipes is increasingly used as a model in biomedical skeletal research. The standard approach is to generate genetic variants with particular skeletal phenotypes which resemble skeletal diseases in humans. The proper diagnosis of skeletal variation is key for this type of research. However, even laboratory rearing conditions can alter skeletal phenotypes. The subject of this study is the link between skeletal phenotypes and rearing conditions. Thus, wildtype medaka were reared from hatching to an early juvenile stage at low (LD: 5 individuals/L), medium (MD: 15 individuals/L), and high (HD: 45 individuals/L) densities. The objectives of the study are: (I) provide a comprehensive overview of the postcranial skeletal elements in medaka; (II) evaluate the effects of rearing density on specific meristic counts and on the variability in type and incidence of skeletal anomalies; (III) define the best laboratory settings to obtain a skeletal reference for a sound evaluation of future experimental conditions; (IV) contribute to elucidating the structural and cellular changes related to the onset of skeletal anomalies. The results from this study reveal that rearing densities greater than 5 medaka/L reduce the animals' growth. This reduction is related to decreased mineralization of dermal (fin rays) and perichondral (fin supporting elements) bone. Furthermore, high density increases anomalies affecting the caudal fin endoskeleton and dermal rays, and the preural vertebral centra. A series of static observations on Alizarin red S whole mount-stained preural fusions provide insights into the etiology of centra fusion. The fusion of preural centra involves the ectopic formation of bony bridges over the intact intervertebral ligament. An apparent consequence is the degradation of the intervertebral ligaments and the remodeling and reshaping of the fused vertebral centra into a biconoid-shaped centrum. From this study it can be concluded that it is paramount to take into account the rearing conditions, natural variability, skeletal phenotypic plasticity, and the genetic background along with species-specific peculiarities when screening for skeletal phenotypes of mutant or wildtype medaka.

Mechanical Static Force Negatively Regulates Vitality and Early Skeletal Development in Zebrafish Embryos

Skeletal system development and remodelling is regulated by several different factors, including hormones, cytokines, and mechanical forces. It is known that gravity and pressure stimulate mechanosensors on bone cells which transduce mechanical signals to chemical ones. Nevertheless, few data have been provided about the role of mechanical forces on embryo osteogenesis in vivo. Since the zebrafish is an elective model for developmental studies, in particular on bone formation and tissue mineralization, we analyzed in vivo the effects of a static mechanical force generated by a water column on fertilized zebrafish eggs. The results have shown that an increase in the hydrostatic pressure (HP) of up to 5.9% was lethal for 100% of treated embryos at 48 h post fertilization (hpf). A small decrease in length (−2%) and 49% mortality were found in the +4.4% HP embryos compared with the controls. To analyze skeletal development, we evaluated the number of mineralized vertebral bodies in the trunk at five days post fertilization. The embryos grown under +2.4% HP showed a physiological intramembranous mineralization of vertebral bodies whereas the embryos which grew with +3.4% HP showed a significant decrease in mineralization rate (−54%). Morphological analysis of cartilage and bones in embryos at +3.4% HP revealed a delay of both intramembranous and chondrogenic mineralization, respectively, in axial and head bones, whereas the chondrogenesis appeared normal. These data suggested that developing osteoblasts and different mineralization programs are sensitive to mechanical pressure when applied to early embryogenesis.

Functional Validation of Osteoporosis Genetic Findings Using Small Fish Models

The advancement of human genomics has revolutionized our understanding of the genetic architecture of many skeletal diseases, including osteoporosis. However, interpreting results from human association studies remains a challenge, since index variants often reside in non-coding regions of the genome and do not possess an obvious regulatory function. To bridge the gap between genetic association and causality, a systematic functional investigation is necessary, such as the one offered by animal models. These models enable us to identify causal mechanisms, clarify the underlying biology, and apply interventions. Over the past several decades, small teleost fishes, mostly zebrafish and medaka, have emerged as powerful systems for modeling the genetics of human diseases. Due to their amenability to genetic intervention and the highly conserved genetic and physiological features, fish have become indispensable for skeletal genomic studies. The goal of this review is to summarize the evidence supporting the utility of Zebrafish (Danio rerio) for accelerating our understanding of human skeletal genomics and outlining the remaining gaps in knowledge. We provide an overview of zebrafish skeletal morphophysiology and gene homology, shedding light on the advantages of human skeletal genomic exploration and validation. Knowledge of the biology underlying osteoporosis through animal models will lead to the translation into new, better and more effective therapeutic approaches.

Hydrodynamic stress and phenotypic plasticity of the zebrafish regenerating fin

Understanding how extrinsic factors modulate genetically encoded information to produce a specific phenotype is of prime scientific interest. In particular, the feedback mechanism between abiotic forces and locomotory organs during morphogenesis to achieve efficient movement is a highly relevant example of such modulation. The study of this developmental process can provide unique insights on the transduction of cues at the interface between physics and biology. Here, we take advantage of the natural ability of adult zebrafish to regenerate their amputated fins to assess its morphogenic plasticity upon external modulations. Using a variety of surgical and chemical treatments, we could induce phenotypic responses to the structure of the fin. Through the ablation of specific rays in regenerating caudal fins, we generated artificially narrowed appendages in which the fin cleft depth and the positioning of rays bifurcations were perturbed compared with normal regenerates. To dissect the role of mechanotransduction in this process, we investigated the patterns of hydrodynamic forces acting on the surface of a zebrafish fin during regeneration by using particle tracking velocimetry on a range of biomimetic hydrofoils. This experimental approach enabled us to quantitatively compare hydrodynamic stress distributions over flapping fins of varying sizes and shapes. As a result, viscous shear stress acting on the distal margin of regenerating fins and the resulting internal tension are proposed as suitable signals for guiding the regulation of ray growth dynamics and branching pattern. Our findings suggest that mechanical forces are involved in the fine-tuning of the locomotory organ during fin morphogenesis.

Plasticity of the skeleton and skeletal deformities in zebrafish (Danio rerio) linked to rearing density

The teleost zebrafish (Danio rerio), an established model for human skeletal diseases, is reared under controlled conditions with defined parameters for temperature and photoperiod. Studies aimed at defining the proper rearing density have been performed with regard to behavioural and physiological stress response, sex ratio and reproduction. Studies concerning the effect of rearing density on the skeletal phenotype are lacking. This study analyses the response of the skeleton to different rearing densities and describes the skeletal deformities. Wild-type zebrafish were reared up to 30 dpf (days post-fertilization) in a common environment. From 30 to 90 dpf, animals were reared at three different densities: high density (HD), 32 fish l^-1 ; medium density (MD), 8 fish l^-1 and low density (LD), 2 fish l^-1 . Animals at 30 and 90 dpf were collected and whole-mount stained with Alizarin red S to visualize mineralized tissues. The entire skeleton was analysed for meristic counts and 172 types of deformities. The results showed that the rearing density significantly influenced the specimens' average standard length, which decreased with the increase in the rearing density. Differences in meristic counts among the three groups were not observed. Rearing density-independent malformations affected the ribs, neural arches and the spines of the abdominal region, as well as vertebrae of the caudal complex. The HD group showed the highest number of deformities per specimen, the highest number of observed types of deformities and, together with the MD group, the highest frequency of specimens affected by severe deformities. In particular, the HD group showed deformities affecting arches, spines and vertebral centra in the caudal region of the vertebral column. This study provides evidence of an effect of the rearing density on the development of different skeletal phenotypes.

Exercise-induced lordosis in zebrafish Danio rerio (Hamilton, 1822)

The anabolic effect of exercise on muscles and bones is well documented. In teleost fish, exercise has been shown to accelerate skeletogenesis, to increase bone volume, and to change the shape of vertebral bodies. Still, increased swimming has also been reported to induce malformations of the teleost vertebral column, particularly lordosis. This study examines whether zebrafish (Danio rerio) develops lordosis as a result of continuous physical exercise. Zebrafish were subjected, for 1 week, to an increased swimming exercise of 5.0, 6.5 or 8.0 total body lengths (TL) per second. Control and exercise group zebrafish were examined for the presence of vertebral abnormalities, by in vivo examination, whole mount staining for bone and cartilage and histology and micro-computed tomography (CT) scanning. Exercise zebrafish developed a significantly higher rate of lordosis in the haemal part of the vertebral column. At the end of the experiment, the frequency of lordosis in the control groups was 0.5 ± 1.3% and that in the exercise groups was 7.5 ± 10.6%, 47.5 ± 10.6% and 92.5 ± 6.0% of 5.0, 6.5 and 8.0 TL∙s^-1 , respectively. Histological analysis and CT scanning revealed abnormal vertebrae with dorsal folding of the vertebral body end plates. Possible mechanisms that trigger lordotic spine malformations are discussed. This is the first study to report a quick, reliable and welfare-compatible method of inducing skeletal abnormalities in a vertebrate model during the post-embryonic period.

A comparison of the structure, composition and mechanical properties of anosteocytic vertebrae of medaka (O. latipes) and osteocytic vertebrae of zebrafish (D. rerio)

Medaka (O. latipes) and zebrafish (D. rerio) are two teleost fish increasingly used as models to study human skeletal diseases. Although they are similar in size, swimming pattern and many other characteristics, these two species are very distant from an evolutionary point of view (by at least 100 million years). A prominent difference between the skeletons of medaka and zebrafish is the total absence of osteocytes in medaka (anosteocytic), while zebrafish bone contains numerous osteocytes (osteocytic). This fundamental difference suggests the possibility that the bony elements of their skeleton may be different in a variety of other aspects, structural, mechanical or both, particularly in heavily loaded bones like the vertebrae. Here we report on the results of a comparative study that aimed to determine the similarities and differences in medaka and zebrafish vertebrae in terms of their macro- to nanostructure, composition and mechanical properties. Our results reveal many similarities between medaka and zebrafish vertebrae, making the lack or presence of osteocytes the only major difference between the bones of these two species.

Novel model of restricted mobility induced osteopenia in zebrafish

Immobilization, such as prolonged bed rest, is a risk factor for bone loss in humans. Motivated by the emerging utility of zebrafish (Danio rerio) as an animal of choice for the study of musculoskeletal disease, here we report a model of restricted mobility induced osteopenia in adult zebrafish. Aquatic tanks with small cubical compartments to restrict the movement and locomotion of single fish were designed and fabricated for this study. Adult zebrafish were divided into two groups: a normal control (CONT) and a restricted mobility group (RMG) (18 fish/group). Six fish from each group were euthanized on days 14, 21 and 35 of the movement restriction. By using microcomputed tomography (micro-CT), we assessed bone volume/tissue volume (BV/TV) and bone density in the whole skeleton of the fish. Furthermore, we assessed skeletal shape in the vertebrae (radius, length, volume, neural and haemal arch aperture areas, neural and haemal arch angle, and thickness of the intervertebral space), single vertebra bone volume and bone density. Movement restriction significantly decreased vertebral skeletal parameters such as radius, length, volume, arch aperture areas and angles as well as the thickness of the intervertebral space in RMG. Furthermore, restricted mobility significantly (P < 0.001) decreased BV/TV and bone density as compared to the CONT group, starting as early as 14 days. By analysing zebrafish from CONT and RMG, we show that micro-CT imaging is a sensitive method to quantify distinct skeletal properties in zebrafish. We further defined the micro-CT parameters which can be used to examine the effects of restricted mobility on the skeleton of the fish. Our findings propose a rapid and effective osteopenia "stabulation" model, which could be used widely for osteoporosis drug screening.

Skeletal Biology and Disease Modeling in Zebrafish

Zebrafish are teleosts (bony fish) that share with mammals a common ancestor belonging to the phylum Osteichthyes, from which their endoskeletal systems have been inherited. Indeed, teleosts and mammals have numerous genetically conserved features in terms of skeletal elements, ossification mechanisms, and bone matrix components in common. Yet differences related to bone morphology and function need to be considered when investigating zebrafish in skeletal research. In this review, we focus on zebrafish skeletal architecture with emphasis on the morphology of the vertebral column and associated anatomical structures. We provide an overview of the different ossification types and osseous cells in zebrafish and describe bone matrix composition at the microscopic tissue level with a focus on assessing mineralization. Processes of bone formation also strongly depend on loading in zebrafish, as we elaborate here. Furthermore, we illustrate the high regenerative capacity of zebrafish bones and present some of the technological advantages of using zebrafish as a model. We highlight zebrafish axial and fin skeleton patterning mechanisms, metabolic bone disease such as after immunosuppressive glucocorticoid treatment, as well as osteogenesis imperfecta (OI) and osteopetrosis research in zebrafish. We conclude with a view of why larval zebrafish xenografts are a powerful tool to study bone metastasis. © 2021 American Society for Bone and Mineral Research (ASBMR).

Physiological Effects of Water Flow Induced Swimming Exercise in Seabream Sparus aurata

A longer on-land rearing period of Gilthead seabream Sparus aurata before transfer to sea-cages would allow the farmer to benefit from exercise-enhanced growth, resilience, and robustness as induced by increasing water flow in the tanks. In this study, the physiological effects of flow-conditioning were investigated by subjecting large groups of experimental fish to minimal flow or to flow regimes inducing swimming exercise at 1 or 2 body length (BL) s⁻¹ for a period of 8 months (February–October) in 1,500 L tanks. Fish representing the three treatment groups were then used for: (1) a stress challenge netting test and plasma cortisol measurement (baseline, peaking, and recovery levels), (2) blood plasma measurements of glucose, triglycerides, lactate, cholesterol, growth hormone (GH), and insulin-like growth factor 1 (IGF1), and (3) heart and muscle gene expression of the GH and IGF1 receptors and the muscle transcriptome by deep RNA sequencing (RNAseq). Fish size after 8 months of flow conditioning was 92 ± 27 g body weight (BW) for fish under minimal flow, 106 ± 24 g BW (+15%) at 1 BL s⁻¹, and 125 ± 27 g BW (+36%) at 2 BL s⁻¹. Flow conditioning at 1 BL s⁻¹ provided optimal conditions for growth and uniformity, but also stress (lowest baseline plasma cortisol), robustness (higher condition factor and larger hearts), and energy mobilization (increased plasma glucose). Although flow enhanced growth linearly with swimming speed, also the percentage of lordotic fish increased with exercise, particularly high for swimming at 2 BL s⁻¹. The absence of important differences in plasma GH and IGF1, and expression levels of their receptors in heart and white skeletal muscle, indicated that other factors may be involved in growth enhancement. RNAseq of the white skeletal muscle showed upregulated expression of genes involved in muscle contraction, muscle development and its molecular regulation, and immune genes that may play a role in the muscle repair mechanism. An exercise regime of swimming at 1 BL s⁻¹ can be considered as optimal for farming robust seabream although the increase of skeletal deformities should be avoided.

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: A Suitable Tool for the Study of Cell Signaling in Bone

In recent decades, many studies using the zebrafish model organism have been performed. Zebrafish, providing genetic mutants and reporter transgenic lines, enable a great number of studies aiming at the investigation of signaling pathways involved in the osteoarticular system and at the identification of therapeutic tools for bone diseases. In this review, we will discuss studies which demonstrate that many signaling pathways are highly conserved between mammals and teleost and that genes involved in mammalian bone differentiation have orthologs in zebrafish. We will also discuss as human diseases, such as osteogenesis imperfecta, osteoarthritis, osteoporosis and Gaucher disease can be investigated in the zebrafish model.

Terrestrial acclimation and exercise lead to bone functional response in Polypterus senegalus pectoral fins

The ability of bones to sense and respond to mechanical loading is a central feature of vertebrate skeletons. However, the functional demands imposed on terrestrial and aquatic animals differ vastly. The pectoral girdle of the basal actinopterygian fish Polypterus senegalus was previously shown to exhibit plasticity following terrestrial acclimation, but the pectoral fin itself has yet to be examined. We investigated skeletal plasticity in the pectoral fins of P. senegalus after exposure to terrestrial loading. Juvenile fish were divided into three groups: a control group was kept under aquatic conditions without intervention, an exercised group was also kept in water but received daily exercise on land, and a terrestrial group was kept in a chronic semi-terrestrial condition. After 5 weeks, the pectoral fins were cleared and stained with Alcian Blue and Alizarin Red to visualize cartilage and bone, allowing measurements of bone length, bone width, ossification and curvature to be taken for the endochondral radial bones. Polypterus senegalus fin bones responded most strongly to chronic loading in the terrestrial condition. Fish that were reared in a terrestrial environment had significantly longer bones compared with those of aquatic controls, wider propterygia and metapterygia, and more ossified metapterygia and medial radials, and they showed changes in propterygial curvature. Exercised fish also had longer and more ossified medial radials compared with those of controls. Polypterus senegalus fin bones exhibit plasticity in response to novel terrestrial loading. Such plasticity could be relevant for transitions between water and land on evolutionary scales, but key differences between fish and tetrapod bone make direct comparisons challenging.

Zebrafish: An Emerging Model for Orthopedic Research

Advances in next-generation sequencing have transformed our ability to identify genetic variants associated with clinical disorders of the musculoskeletal system. However, the means to functionally validate and analyze the physiological repercussions of genetic variation have lagged behind the rate of genetic discovery. The zebrafish provides an efficient model to leverage genetic analysis in an in vivo context. Its utility for orthopedic research is becoming evident in regard to both candidate gene validation as well as therapeutic discovery in tissues such as bone, tendon, muscle, and cartilage. With the development of new genetic and analytical tools to better assay aspects of skeletal tissue morphology, mineralization, composition, and biomechanics, researchers are emboldened to systematically approach how the skeleton develops and to identify the root causes, and potential treatments, of skeletal disease. © 2019 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 38:925-936, 2020.

Finite element and deformation analyses predict pattern of bone failure in loaded zebrafish spines

The spine is the central skeletal support structure in vertebrates consisting of repeated units of bone, the vertebrae, separated by intervertebral discs (IVDs) that enable the movement of the spine. Spinal pathologies such as idiopathic back pain, vertebral compression fractures and IVD failure affect millions of people worldwide. Animal models can help us to understand the disease process, and zebrafish are increasingly used as they are highly genetically tractable, their spines are axially loaded like humans, and they show similar pathologies to humans during ageing. However, biomechanical models for the zebrafish are largely lacking. Here, we describe the results of loading intact zebrafish spinal motion segments on a material testing stage within a micro-computed tomography machine. We show that vertebrae and their arches show predictable patterns of deformation prior to their ultimate failure, in a pattern dependent on their position within the segment. We further show using geometric morphometrics which regions of the vertebra deform the most during loading, and that finite-element models of the trunk subjected reflect the real patterns of deformation and strain seen during loading and can therefore be used as a predictive model for biomechanical performance.

Using zebrafish to study skeletal genomics

While genome-wide association studies (GWAS) have revolutionized our understanding of the genetic architecture of skeletal diseases, animal models are required to identify causal mechanisms and to translate underlying biology into new therapies. Despite large-scale knockout mouse phenotyping efforts, the skeletal functions of most genes residing at GWAS-identified loci remain unknown, highlighting a need for complementary model systems to accelerate gene discovery. Over the past several decades, zebrafish (Danio rerio) has emerged as a powerful system for modeling the genetics of human diseases. In this review, our goal is to outline evidence supporting the utility of zebrafish for accelerating our understanding of human skeletal genomics, as well as gaps in knowledge that need to be filled for this purpose. We do this by providing a basic foundation of the zebrafish skeletal morphophysiology and phenotypes, and surveying evidence of skeletal gene homology and the use of zebrafish for post-GWAS analysis in other tissues and organs. We also outline challenges in translating zebrafish mutant phenotypes. Finally, we conclude with recommendations of future directions and how to leverage the large body of tools and knowledge of skeletal genetics in zebrafish for the needs of human skeletal genomic exploration. Due to their amenability to rapid genetic approaches, as well as the large number of conserved genetic and phenotypic features, there is a strong rationale supporting the use of zebrafish for human skeletal genomic studies.

Gill remodelling during terrestrial acclimation in the amphibious fish Polypterus senegalus

Fishes are effectively weightless in water due to the buoyant support of the environment, but amphibious fishes must cope with increased effective weight when on land. Delicate structures such as gills are especially vulnerable to collapse and loss of surface area out of water. We tested the 'structural support' hypothesis that amphibious Polypterus senegalus solve this problem using phenotypically plastic changes that provide mechanical support and increase stiffness at the level of the gill lamellae, the filaments, and the whole arches. After 7 d in terrestrial conditions, enlargement of an inter-lamellar cell mass filled the water channels between gill lamellae, possibly to provide structural support and/or reduce evaporative water loss. Similar gill remodelling has been described in several other actinopterygian fishes, suggesting this may be an ancestral trait. There was no change in the mechanical properties or collagen composition of filaments or arches after 7 days out of water, but 8 months of terrestrial acclimation caused a reduction in gill arch length and mineralized bone volume. Thus, rather than increasing the size and stiffness of the gill skeleton, P. senegalus may instead reduce investment in supportive gill tissue while on land. These results are strikingly similar to the evolutionary trend of gill loss that occurred during the tetrapod invasion of land, raising the possibility that genetic assimilation of gill plasticity was an underlying mechanism.

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.

Age, but not short-term intensive swimming, affects chondrocyte turnover in zebrafish vertebral cartilage

Both age and intensive exercise are generally considered critical risk factors for osteoarthritis. In this work, we intend to establish zebrafish models to assess the role of these two factors on cartilage homeostasis. We designed a swimming device for zebrafish intensive exercise. The body measurements, bone mineral density (BMD) and the histology of spinal cartilages of 4- and 12-month-old zebrafish, as well the 12-month-old zebrafish before and after a 2-week exercise were compared. Our results indicate that both age and exercise affect the body length and body weight, and the micro-computed tomography reveals that both age and exercise affect the spinal BMD. However, quantitative analysis of immunohistochemistry and histochemistry indicate that short-term intensive exercise does not affect the extracellular matrix (ECM) of spinal cartilage. On the other hand, the cartilage ECM significantly grew from 4 to 12 months of age with an increase in total chondrocytes. dUTP nick end labeling staining shows that the percentages of apoptotic cells significantly increase as the zebrafish grows, whereas the BrdU labeling shows that proliferative cells dramatically decrease from 4 to 12 months of age. A 30-day chase of BrdU labeling shows some retention of labeling in cells in 4-month-old spinal cartilage but not in cartilage from 12-month-old zebrafish. Taken together, our results suggest that zebrafish chondrocytes are actively turned over, and indicate that aging is a critical factor that alters cartilage homeostasis. Zebrafish vertebral cartilage may serve as a good model to study the maturation and homeostasis of articular cartilage.

Stage-dependent effects of chlorpyrifos on medaka (Oryzias latipes) swimming behavior using a miniaturized swim flume

By considering chlorpyrifos (CPF), an organophosphorus pesticide with known mechanisms of action that affect neurobehavioral development, we assessed the validity and sensitivity of a miniaturized swim flume by investigating the effects of the insecticide on swimming behavior in medaka (Oryzias latipes) fish growing stages. Medaka in three developmental periods, namely 0, 20 and 40 day-old post-hatch (i.e. time points 0, 20 and 40, respectively), were exposed to CPF (12.5, 25, 50 and 100 μg/L) for 48 h under semi-static conditions. The CPF half-lives during exposures were evaluated and the swimming patterns in a flume section (arena) were presented on two-dimensional gradient maps of forced movement of fish against water current. A comparative numerical analysis of fish residence times between each time point control and the corresponding CPF groups was performed by dividing arenas into 15 proportional areas. The time point 0 control group gradient map showed a noticeably different swim pattern from those of the ≥12.5 μg CPF/L groups, which was statistically supported by the differences for residence times seen in ≥12 corresponding areas. The control group gradient maps for time points 20 and 40 differed from those of the respective ≥12.5 μg CPF/L groups. The comparative analysis of the residence times in the corresponding 15 areas revealed differences in ≥5 areas for time point 20 and in ≥3 areas for time point 40. The integrative analysis of the gradient maps and the numerical statistics revealed stage-specific effects and a concentration-response relationship between CPF and alterations on forced medaka swimming despite the dissipation of CPF from the water column. These results indicate the validity of the miniaturized swim flume toward a more environmentally realistic scenario for the evaluation of neurodevelopmental and behavioral toxicity in small fish models.

An open source microcontroller based flume for evaluating swimming performance of larval, juvenile, and adult zebrafish

Zebrafish are a preferred vertebrate model for delineating genotype-phenotype relationships. One of the most studied features of zebrafish is their exceptional swimming ability. By 7 days postfertilization (dpf), zebrafish spend over two-thirds of their time engaged in spontaneous swimming activity and several months later they are capable of attaining some of the fastest swimming velocities relative to body length ever recorded in the laboratory. However, laboratory-assembled flumes capable of achieving the slow flow velocities characteristics of larvae as well as the relatively fast maximal velocities of adults have not been described in sufficient detail to allow easy replication. Here we describe an easily assembled, open-source zebrafish-scaled flume for assessing swimming performance. The flume uses two independent spherical-impeller pumps modulated by a microcontroller to achieve flow velocities ranging from 1 to 70 cm s⁻¹. The microcontroller also monitors water temperature and flow velocity and sends these data to a personal computer for real-time display and storage. Incremental protocols for assessing maximal swimming speed (U_max) were developed, stored in custom software, and then uploaded to the microcontroller in order to assess performance of larval (14, 21, 28 dpf), juvenile (35, 42 dpf), and adult (8, 22 month) zebrafish. The flume had sufficient range and sensitivity to detect developmental changes in U_max of larvae and juveniles, an 18–24% faster U_max of adult males vs. females, and a 14–20% age-related reduction in U_max for the oldest zebrafish. Detailed information is provided to assemble and operate this low-cost, versatile, and reliable tool for assessing zebrafish swimming performance.

Evolution of caudal fin ray development and caudal fin hypural diastema complex in spotted gar, teleosts, and other neopterygian fishes

BACKGROUND

The caudal fin of actinopterygians transitioned from a heterocercal dorsoventrally asymmetrical fin to a homocercal externally symmetrical fin in teleosts through poorly understood evolutionary developmental mechanisms. We studied the caudal skeleton of major living actinopterygian lineages, including polypteriformes, acipenseriformes, Holostei (gars and bowfin), and teleosts, compared with reports of extinct neopterygians and basal teleosteans. We focused on the hypural diastema complex, which includes (1) a gap between hypurals 2 and 3, that (2) separates two plates of connective tissue at (3) the branching of caudal vasculature; these features had been considered as a shared, derived trait of teleosts, a synapomorphy.

RESULTS

These studies revealed that gars and teleosts share all three features of the hypural diastema complex. Absence of a complex with these features from bowfin, fossil Holostei, and stem Teleostei argues in favor of repetitive, independent emergence in several neopterygian and basal Teleostei lineages, or less likely, many independent losses. We further observed that, in gars and teleosts, the earliest developing lepidotrichia align with the horizontal adult body axis, thus participating in external symmetry.

CONCLUSIONS

These results suggest that the hypural diastema complex in teleosts and gars represents a homoplasy among neopterygians and that it emerged repeatedly by parallel evolution due to shared inherited underlying genetic and developmental programs (latent homology). Because the hypural diastema complex exists in gars with heterocercal tails, this complex is independent of homocercality. Developmental Dynamics 247:832-853, 2018. © 2018 Wiley Periodicals, Inc.

Zebrafish as a model to study bone maturation: Nanoscale structural and mechanical characterization of age-related changes in the zebrafish vertebral column

Zebrafish (Danio rerio) is a useful model for understanding biomedical properties of bone and are widely employed in developmental and genetics studies. Here, we have studied the development of zebrafish vertebral bone at the nanoscale from adolescence (6 months), early adulthood (10 months) to mid-life (14 months). Characterization of the bone was conducted using energy-dispersive X-ray spectroscopy (EDX), scanning electron microscopy (SEM) and PeakForce QNM atomic force microscopy (AFM) techniques. SEM and AFM revealed a lamellar structure with mineralized collagen fibrils. There was a significant increase in the wall thickness from 6 to 10 months (76%) and 10 months to 14 months (26%), which is positively correlated with nanomechanical behavior. An increase in the Ca/P ratio was found which was also positively correlated with nanomechanical properties. The change in mechanical properties and Ca/P are similar to those expected in humans when transitioning from adolescence to mid-life. We suggest that zebrafish serve as a suitable model for further studies on age-related changes in bone ultrastructure.

Skeletal development in the heterocercal caudal fin of spotted gar (lepisosteus oculatus) and other lepisosteiformes

BACKGROUND

The caudal fin of actinopterygians experienced substantial morphological changes during evolution. In basal actinopterygians, the caudal fin skeleton supports an asymmetrical heterocercal caudal fin, while most teleosts have a symmetrical homocercal caudal fin. The transition from the ancestral heterocercal form to the derived homocercal caudal fin remains poorly understood. Few developmental studies provide an understanding of derived and ancestral characters among basal actinopterygians. To fill this gap, we examined the development of the caudal fin of spotted gar Lepisosteus oculatus, one of only eight living species of Holostei, the sister group to the teleosts.

RESULTS

Our observations of animals from fertilization to more than a year old provide the most detailed description of the development of caudal fin skeletal elements in any Holostean species. We observed two different types of distal caudal radials replacing two transient plates of connective tissue, identifying two hypaxial ensembles separated by a space between hypurals 2 and 3. These features have not been described in any gar species, but can be observed in other gar species, and thus represent anatomical structures common to lepisosteiformes.

CONCLUSIONS

The present work highlights the power and importance of ontogenic studies and provides bases for future evolutionary and morphological investigations on actinopterygians fins. Developmental Dynamics 247:724-740, 2018. © 2018 Wiley Periodicals, Inc.

Increased mechanical loading through controlled swimming exercise induces bone formation and mineralization in adult zebrafish

Exercise promotes gain in bone mass through adaptive responses of the vertebrate skeleton. This mechanism counteracts age- and disease-related skeletal degradation, but remains to be fully understood. In life sciences, zebrafish emerged as a vertebrate model that can provide new insights into the complex mechanisms governing bone quality. To test the hypothesis that musculoskeletal exercise induces bone adaptation in adult zebrafish and to characterize bone reorganization, animals were subjected to increased physical exercise for four weeks in a swim tunnel experiment. Cellular, structural and compositional changes of loaded vertebrae were quantified using integrated high-resolution analyses. Exercise triggered rapid bone adaptation with substantial increases in bone-forming osteoblasts, bone volume and mineralization. Clearly, modeling processes in zebrafish bone resemble processes in human bone. This study highlights how exercise experiments in adult zebrafish foster in-depth insight into aging-related bone diseases and can thus catalyze the search for appropriate prevention and new treatment options.

Effects of dietary 2,2', 4,4'-tetrabromodiphenyl ether (BDE-47) exposure on medaka (Oryzias latipes) swimming behavior

A diet fortified with 2,2', 4,4'-tetrabromodiphenyl ether (BDE-47: 0, 10, 100, and 1000 ng/g) was dosed to 4-7-day-old post-hatch medaka fish for 40 days to evaluate the effects on the swimming activity of fish using a miniaturized swimming flume. Chlorpyrifos (CF)-exposed fish were selected as the positive control to assess the validity and sensitivity of the behavioral findings. After 20 and 40 days of exposure, the locomotor activity was analyzed for 6 min in a flume section (arena). The CF positive control for each time point were fish exposed to 50 ng CF/ml for 48 h. Swimming patterns, presented as two-dimensional heat maps of fish movement and positioning, were obtained by geostatistical analyses. The heat maps of the control groups at time point 20 revealed visually comparable swimming patterns to those of the BDE-47-treated groups. For the comparative fish positioning analysis, both the arenas were divided into 15 proportional areas. No statistical differences were found between residence times in the areas from the control groups and those from the BDE-47-treated groups. At time point 40, the heat map overall patterns of the control groups differed visually from that of the 100-ng BDE-47/g-treated group, but a comparative analysis of the residence times in the corresponding 15 areas did not reveal consistent differences. The relative distances traveled by the control and treated groups at time points 20 and 40 were also comparable. The heat maps of CF-treated fish at both time points showed contrasting swim patterns with respect to those of the controls. These differential patterns were statistically supported with differences in the residence times for different areas. The relative distances traveled by the CF-treated fish were also significantly shorter. These results confirm the validity of the experimental design and indicate that a dietary BDE-47 exposure does not affect forced swimming in medaka at growing stages.

MicroCT-based phenomics in the zebrafish skeleton reveals virtues of deep phenotyping in a distributed organ system

Phenomics, which ideally involves in-depth phenotyping at the whole-organism scale, may enhance our functional understanding of genetic variation. Here, we demonstrate methods to profile hundreds of phenotypic measures comprised of morphological and densitometric traits at a large number of sites within the axial skeleton of adult zebrafish. We show the potential for vertebral patterns to confer heightened sensitivity, with similar specificity, in discriminating mutant populations compared to analyzing individual vertebrae in isolation. We identify phenotypes associated with human brittle bone disease and thyroid stimulating hormone receptor hyperactivity. Finally, we develop allometric models and show their potential to aid in the discrimination of mutant phenotypes masked by alterations in growth. Our studies demonstrate virtues of deep phenotyping in a spatially distributed organ system. Analyzing phenotypic patterns may increase productivity in genetic screens, and facilitate the study of genetic variants associated with smaller effect sizes, such as those that underlie complex diseases.

Broadening the antibacterial spectrum of histidine kinase autophosphorylation inhibitors via the use of ε-poly-L-lysine capped mesoporous silica-based nanoparticles

Two-component systems (TCS) regulate diverse processes such as virulence, stress responses, metabolism and antibiotic resistance in bacteria but are absent in humans, making them promising targets for novel antibacterials. By incorporating recently described TCS histidine kinase autophosphorylation inhibitors (HKAIs) into ε-poly-L-lysine capped nanoparticles (NPs) we could overcome the Gram negative (Gr-) permeability barrier for the HKAIs. The observed bactericidal activity against Gr- bacteria was shown to be due to the enhanced delivery and internalization of the HKAIs and not an inhibitory or synergistic effect of the NPs. The NPs had no adverse effects on mammalian cell viability or the immune function of macrophages in vitro and showed no signs of toxicity to zebrafish larvae in vivo. These results show that HKAIs are promising antibacterials for both Gr- and Gr+pathogens and that NPs are a safe drug delivery technology that can enhance the selectivity and efficacy of HKAIs against bacteria.

Physiological and Behavioral Effects of Exposure to Environmentally Relevant Concentrations of Prednisolone During Zebrafish (Danio rerio) Embryogenesis

The presence of synthetic glucocorticoids within the aquatic environment has been highlighted as a potential environmental concern as they may mimic the role of endogenous glucocorticoids during vertebrate ontogeny. Prednisolone is a commonly prescribed synthetic glucocorticoid which has been repeatedly detected in the environment. This study investigated the impact of environmentally relevant concentrations of prednisolone (0.1, 1, and 10 μg/L) during zebrafish embryogenesis using physiological and behavioral end points which are known to be mediated by endogenous glucocorticoids. The frequency of spontaneous muscle contractions (24 hpf) was significantly reduced by prednisolone and 0.1 μg/L increased the distance embryos swam in response to a mechanosensory stimulus (48 hpf). The percentage of embryos hatched significantly increased following prednisolone treatment (1 and 10 μg/L), while growth and mortality were unaffected. The onset of heart contraction was differentially affected by prednisolone while heart rate and oxygen consumption both increased significantly throughout embryogenesis. No substantial effect on the axial musculature was observed. Morphological changes to the lower jaw were detected at 96 hpf in response to 1 μg/L of prednisolone. Several parameters of swim behavior were also significantly affected. Environmentally relevant concentrations of prednisolone therefore alter early zebrafish ontogeny and significantly affect embryo behavior.

A primary phosphorus-deficient skeletal phenotype in juvenile Atlantic salmon Salmo salar: the uncoupling of bone formation and mineralization

To understand the effect of low dietary phosphorus (P) intake on the vertebral column of Atlantic salmon Salmo salar, a primary P deficiency was induced in post-smolts. The dietary P provision was reduced by 50% for a period of 10 weeks under controlled conditions. The animal's skeleton was subsequently analysed by radiology, histological examination, histochemical detection of minerals in bones and scales and chemical mineral analysis. This is the first account of how a primary P deficiency affects the skeleton in S. salar at the cellular and at the micro-anatomical level. Animals that received the P-deficient diet displayed known signs of P deficiency including reduced growth and soft, pliable opercula. Bone and scale mineral content decreased by c. 50%. On radiographs, vertebral bodies appear small, undersized and with enlarged intervertebral spaces. Contrary to the X-ray-based diagnosis, the histological examination revealed that vertebral bodies had a regular size and regular internal bone structures; intervertebral spaces were not enlarged. Bone matrix formation was continuous and uninterrupted, albeit without traces of mineralization. Likewise, scale growth continues with regular annuli formation, but new scale matrix remains without minerals. The 10 week long experiment generated a homogeneous osteomalacia of vertebral bodies without apparent induction of skeletal malformations. The experiment shows that bone formation and bone mineralization are, to a large degree, independent processes in the fish examined. Therefore, a deficit in mineralization must not be the only cause of the alterations of the vertebral bone structure observed in farmed S. salar. It is discussed how the observed uncoupling of bone formation and mineralization helps to better diagnose, understand and prevent P deficiency-related malformations in farmed S. salar.

Finite element modelling predicts changes in joint shape and cell behaviour due to loss of muscle strain in jaw development

Abnormal joint morphogenesis is linked to clinical conditions such as Developmental Dysplasia of the Hip (DDH) and to osteoarthritis (OA). Muscle activity is known to be important during the developmental process of joint morphogenesis. However, less is known about how this mechanical stimulus affects the behaviour of joint cells to generate altered morphology. Using zebrafish, in which we can image all joint musculoskeletal tissues at high resolution, we show that removal of muscle activity through anaesthetisation or genetic manipulation causes a change to the shape of the joint between the Meckel's cartilage and Palatoquadrate (the jaw joint), such that the joint develops asymmetrically leading to an overlap of the cartilage elements on the medial side which inhibits normal joint function. We identify the time during which muscle activity is critical to produce a normal joint. Using Finite Element Analysis (FEA), to model the strains exerted by muscle on the skeletal elements, we identify that minimum principal strains are located at the medial region of the joint and interzone during mouth opening. Then, by studying the cells immediately proximal to the joint, we demonstrate that biomechanical strain regulates cell orientation within the developing joint, such that when muscle-induced strain is removed, cells on the medial side of the joint notably change their orientation. Together, these data show that biomechanical forces are required to establish symmetry in the joint during development.

Spinal deformity in aged zebrafish is accompanied by degenerative changes to their vertebrae that resemble osteoarthritis

Age-related degenerative changes within the vertebral column are a significant cause of morbidity with considerable socio-economic impact worldwide. An improved understanding of these changes through the development of experimental models may lead to improvements in existing clinical treatment options. The zebrafish is a well-established model for the study of skeletogenesis with significant potential in gerontological research. With advancing age, zebrafish frequently develop gross deformities of their vertebral column, previously ascribed to reduced trunk muscle tone. In this study, we assess degenerative changes specifically within the bone and cartilage of the vertebral column of zebrafish at 1, 2 and 3-years of age. We show increased frequency and severity of spinal deformities/curvatures with age. Underlying the most severe phenotypes are partial or complete vertebral dislocations and focal thickening of the vertebral bone at the joint margins. MicroCT examination demonstrates small defects, fractures and morphological evidence suggestive of bone erosion and remodeling (i.e. osteophytes) within the vertebrae during aging, but no significant change in bone density. Light and electron microscopic examination reveal striking age-related changes in cell morphology, suggestive of chondroptosis, and tissue remodelling of the vertebral cartilage, particularly within the pericellular micro-environment. Glycosaminoglycan analysis of the vertebral column by HPLC demonstrates a consistent, age-related increase in the yield of total chondroitin sulfate disaccharide, but no change in sulfation pattern, supported by immunohistochemical analysis. Immunohistochemistry strongly identifies all three chondroitin/dermatan sulphate isoforms (C-0-S, C-4-S/DS and C-6-S) within the vertebral cartilage, particularly within the pericellular micro-environment. In contrast, keratan sulfate immunolocalises specifically with the notochordal tissue of the intervertebral disc, and its labelling diminishes with age. In summary, these observations raise the prospect that zebrafish, in addition to modelling skeletal development, may have utility in modelling age-related degenerative changes that affect the skeleton during senescence.

The enigmas of bone without osteocytes

One of the hallmarks of tetrapod bone is the presence of numerous cells (osteocytes) within the matrix. Osteocytes are vital components of tetrapod bone, orchestrating the processes of bone building, reshaping and repairing (modeling and remodeling), and probably also participating in calcium-phosphorus homeostasis via both the local process of osteocytic osteolysis, and systemic effect on the kidneys. Given these critical roles of osteocytes, it is thought-provoking that the entire skeleton of many fishes consists of bone material that does not contain osteocytes. This raises the intriguing question of how the skeleton of these animals accomplishes the various essential functions attributed to osteocytes in other vertebrates, and raises the possibility that in acellular bone some of these functions are either accomplished by non-osteocytic routes or not necessary at all. In this review, we outline evidence for and against the fact that primary functions normally ascribed to osteocytes, such as mechanosensation, regulation of osteoblast/clast activity and mineral metabolism, also occur in fish bone devoid of these cells, and therefore must be carried out through alternative and perhaps ancient pathways. To enable meaningful comparisons with mammalian bone, we suggest thorough, phylogenetic examinations of regulatory pathways, studies of structure and mechanical properties and surveys of the presence/absence of bone cells in fishes. Insights gained into the micro-/nanolevel structure and architecture of fish bone, its mechanical properties and its physiology in health and disease will contribute to the discipline of fish skeletal biology, but may also help answer questions of basic bone biology.