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

SignificanceType I collagenopathies are a heterogenous group of connective tissue disorders, caused by genetic defects in type I collagen. Inherent to these disorders is a large clinical variability, of which the underlying molecular basis remains undefined. By systematically analyzing skeletal phenotypes in a large set of type I collagen zebrafish mutants, we show that zebrafish models are able to both genocopy and phenocopy different forms of human type I collagenopathies, arguing for a similar pathogenetic basis. This study illustrates the future potential of zebrafish as a tool to further dissect the molecular basis of phenotypic variability in human type I collagenopathies, to improve diagnostic strategies as well as promote the discovery of new targetable pathways for pharmacological intervention of these disorders.

High-risk neuroblastoma, a pediatric tumor originating from the sympathetic nervous system, has a low mutation load but highly recurrent somatic DNA copy number variants. Previously, segmental gains and/or amplifications allowed identification of drivers for neuroblastoma development. Using this approach, combined with gene dosage impact on expression and survival, we identified ribonucleotide reductase subunit M2 (RRM2) as a candidate dependency factor further supported by growth inhibition upon in vitro knockdown and accelerated tumor formation in a neuroblastoma zebrafish model coexpressing human RRM2 with MYCN. Forced RRM2 induction alleviates excessive replicative stress induced by CHK1 inhibition, while high RRM2 expression in human neuroblastomas correlates with high CHK1 activity. MYCN-driven zebrafish tumors with RRM2 co-overexpression exhibit differentially expressed DNA repair genes in keeping with enhanced ATR-CHK1 signaling activity. In vitro, RRM2 inhibition enhances intrinsic replication stress checkpoint addiction. Last, combinatorial RRM2-CHK1 inhibition acts synergistic in high-risk neuroblastoma cell lines and patient-derived xenograft models, illustrating the therapeutic potential.

Genomewide association studies (GWAS) have improved our understanding of the genetic architecture of common complex diseases such as osteoporosis. Nevertheless, to attribute functional skeletal contributions of candidate genes to osteoporosis-related traits, there is a need for efficient and cost-effective in vivo functional testing. This can be achieved through CRISPR-based reverse genetic screens, where phenotyping is traditionally performed in stable germline knockout (KO) mutants. Recently it was shown that first-generation (F0) mosaic mutant zebrafish (so-called crispants) recapitulate the phenotype of germline KOs. To demonstrate feasibility of functional validation of osteoporosis candidate genes through crispant screening, we compared a crispant to a stable KO zebrafish model for the lrp5 gene. In humans, recessive loss-of-function mutations in LRP5, a co-receptor in the Wnt signaling pathway, cause osteoporosis-pseudoglioma syndrome. In addition, several GWAS studies identified LRP5 as a major risk locus for osteoporosis-related phenotypes. In this study, we showed that early stage lrp5 KO larvae display decreased notochord mineralization and malformations of the head cartilage. Quantitative micro-computed tomography (micro-CT) scanning and massspectrometry element analysis of the adult skeleton revealed decreased vertebral bone volume and bone mineralization, hallmark features of osteoporosis. Furthermore, regenerating fin tissue displayed reduced Wnt signaling activity in lrp5 KO adults. We next compared lrp5 mutants with crispants. Next-generation sequencing analysis of adult crispant tissue revealed a mean out-of-frame mutation rate of 76%, resulting in strongly reduced levels of Lrp5 protein. These crispants generally showed a milder but nonetheless highly comparable skeletal phenotype and a similarly reduced Wnt pathway response compared with lrp5 KO mutants. In conclusion, we show through faithful modeling of LRP5-related primary osteoporosis that crispant screening in zebrafish is a promising approach for rapid functional screening of osteoporosis candidate genes.