Phenomics-Based Quantification of CRISPR-Induced Mosaicism in Zebrafish

Zebrafish Models for Skeletal and Extraskeletal Osteogenesis Imperfecta Features: Unveiling Pathophysiology and Paving the Way for Drug Discovery

In the last decades, the easy genetic manipulation, the external fertilization, the high percentage of homology with human genes and the reduced husbandry costs compared to rodents, made zebrafish a valid model for studying human diseases and for developing new therapeutical strategies. Since zebrafish shares with mammals the same bone cells and ossification types, it became widely used to dissect mechanisms and possible new therapeutic approaches in the field of common and rare bone diseases, such as osteoporosis and osteogenesis imperfecta (OI), respectively. OI is a heritable skeletal disorder caused by defects in gene encoding collagen I or proteins/enzymes necessary for collagen I synthesis and secretion. Nevertheless, OI patients can be also characterized by extraskeletal manifestations such as dentinogenesis imperfecta, muscle weakness, cardiac valve and pulmonary abnormalities and skin laxity. In this review, we provide an overview of the available zebrafish models for both dominant and recessive forms of OI. An updated description of all the main similarities and differences between zebrafish and mammal skeleton, muscle, heart and skin, will be also discussed. Finally, a list of high- and low-throughput techniques available to exploit both larvae and adult OI zebrafish models as unique tools for the discovery of new therapeutic approaches will be presented.

Loss of cped1 does not affect bone and lean tissue in zebrafish

Human genetic studies have nominated Cadherin-like and PC-esterase Domain-containing 1 (CPED1) as a candidate target gene mediating bone mineral density (BMD) and fracture risk heritability. Recent efforts to define the role of CPED1 in bone in mouse and human models have revealed complex alternative splicing and inconsistent results arising from gene targeting, making its function in bone difficult to interpret. To better understand the role of CPED1 in adult bone mass and morphology, we conducted a comprehensive genetic and phenotypic analysis of cped1 in zebrafish, an emerging model for bone and mineral research. We analyzed two different cped1 mutant lines and performed deep phenotyping to characterize more than 200 measures of adult vertebral, craniofacial, and lean tissue morphology. We also examined alternative splicing of zebrafish cped1 and gene expression in various cell/tissue types. Our studies fail to support an essential role of cped1 in adult zebrafish bone. Specifically, homozygous mutants for both cped1 mutant alleles, which are expected to result in loss-of-function and impact all cped1 isoforms, exhibited no significant differences in the measures examined when compared to their respective wildtype controls, suggesting that cped1 does not significantly contribute to these traits. We identified sequence differences in critical residues of the catalytic triad between the zebrafish and mouse orthologs of CPED1, suggesting that differences in key residues, as well as distinct alternative splicing, could underlie different functions of CPED1 orthologs in the two species. Our studies fail to support a requirement of cped1 in zebrafish bone and lean tissue, adding to evidence that variants at 7q31.31 can act independently of CPED1 to influence BMD and fracture risk.

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

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

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.

Zebrafish as a Model for Osteoporosis: Functional Validations of Genome-Wide Association Studies

PURPOSE OF REVIEW

GWAS, as a largely correlational analysis, requires in vitro or in vivo validation. Zebrafish (Danio rerio) have many advantages for studying the genetics of human diseases. Since gene editing in zebrafish has been highly valuable for studying embryonic skeletal developmental processes that are prenatally or perinatally lethal in mammalian models, we are reviewing pros and cons of this model.

RECENT FINDINGS

The true power for the use of zebrafish is the ease by which the genome can be edited, especially using the CRISPR/Cas9 system. Gene editing, followed by phenotyping, for complex traits such as BMD, is beneficial, but the major physiological differences between the fish and mammals must be considered. Like mammals, zebrafish do have main bone cells; thus, both in vivo stem cell analyses and in vivo imaging are doable. Yet, the "long" bones of fish are peculiar, and their bone cavities do not contain bone marrow. Partial duplication of the zebrafish genome should be taken into account. Overall, small fish toolkit can provide unmatched opportunities for genetic modifications and morphological investigation as a follow-up to human-first discovery.

Syntaxin 18 Defects in Human and Zebrafish Unravel Key Roles in Early Cartilage and Bone Development

SNARE proteins comprise a conserved protein family responsible for catalyzing membrane fusion during vesicle traffic. Syntaxin18 (STX18) is a poorly characterized endoplasmic reticulum (ER)-resident t-SNARE. Recently, together with TANGO1 and SLY1, its involvement was shown in ER to Golgi transport of collagen II during chondrogenesis. We report a fetus with a severe osteochondrodysplasia in whom we identified a homozygous substitution of the highly conserved p.Arg10 to Pro of STX18. CRISPR/Cas9-mediated Stx18 deficiency in zebrafish reveals a crucial role for Stx18 in cartilage and bone development. Furthermore, increased expression of multiple components of the Stx18 SNARE complex and of COPI and COPII proteins suggests that Stx18 deficiency impairs antero- and retrograde vesicular transport in the crispant stx18 zebrafish. Taken together, our studies highlight a new candidate gene for a recessive form of osteochondrodysplasia, thereby possibly broadening the SNAREopathy phenotypic spectrum and opening new doors toward future research avenues. © 2023 American Society for Bone and Mineral Research (ASBMR).

Why Animal Experiments Are Still Indispensable in Bone Research: A Statement by the European Calcified Tissue Society

Major achievements in bone research have always relied on animal models and in vitro systems derived from patient and animal material. However, the use of animals in research has drawn intense ethical debate and the complete abolition of animal experimentation is demanded by fractions of the population. This phenomenon is enhanced by the reproducibility crisis in science and the advance of in vitro and in silico techniques. 3D culture, organ-on-a-chip, and computer models have improved enormously over the last few years. Nevertheless, the overall complexity of bone tissue cross-talk and the systemic and local regulation of bone physiology can often only be addressed in entire vertebrates. Powerful genetic methods such as conditional mutagenesis, lineage tracing, and modeling of the diseases enhanced the understanding of the entire skeletal system. In this review endorsed by the European Calcified Tissue Society (ECTS), a working group of investigators from Europe and the US provides an overview of the strengths and limitations of experimental animal models, including rodents, fish, and large animals, as well the potential and shortcomings of in vitro and in silico technologies in skeletal research. We propose that the proper combination of the right animal model for a specific hypothesis and state-of-the-art in vitro and/or in silico technology is essential to solving remaining important questions in bone research. This is crucial for executing most efficiently the 3R principles to reduce, refine, and replace animal experimentation, for enhancing our knowledge of skeletal biology, and for the treatment of bone diseases that affect a large part of society. © 2023 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

Examining craniofacial variation among crispant and mutant zebrafish models of human skeletal diseases

Genetic diseases affecting the skeletal system present with a wide range of symptoms that make diagnosis and treatment difficult. Genome-wide association and sequencing studies have identified genes linked to human skeletal diseases. Gene editing of zebrafish models allows researchers to further examine the link between genotype and phenotype, with the long-term goal of improving diagnosis and treatment. While current automated tools enable rapid and in-depth phenotyping of the axial skeleton, characterizing the effects of mutations on the craniofacial skeleton has been more challenging. The objective of this study was to evaluate a semi-automated screening tool can be used to quantify craniofacial variations in zebrafish models using four genes that have been associated with human skeletal diseases (meox1, plod2, sost, and wnt16) as test cases. We used traditional landmarks to ground truth our dataset and pseudolandmarks to quantify variation across the 3D cranial skeleton between the groups (somatic crispant, germline mutant, and control fish). The proposed pipeline identified variation between the crispant or mutant fish and control fish for four genes. Variation in phenotypes parallel human craniofacial symptoms for two of the four genes tested. This study demonstrates the potential as well as the limitations of our pipeline as a screening tool to examine multi-dimensional phenotypes associated with the zebrafish craniofacial skeleton.

EmbryoNet: using deep learning to link embryonic phenotypes to signaling pathways

Evolutionarily conserved signaling pathways are essential for early embryogenesis, and reducing or abolishing their activity leads to characteristic developmental defects. Classification of phenotypic defects can identify the underlying signaling mechanisms, but this requires expert knowledge and the classification schemes have not been standardized. Here we use a machine learning approach for automated phenotyping to train a deep convolutional neural network, EmbryoNet, to accurately identify zebrafish signaling mutants in an unbiased manner. Combined with a model of time-dependent developmental trajectories, this approach identifies and classifies with high precision phenotypic defects caused by loss of function of the seven major signaling pathways relevant for vertebrate development. Our classification algorithms have wide applications in developmental biology and robustly identify signaling defects in evolutionarily distant species. Furthermore, using automated phenotyping in high-throughput drug screens, we show that EmbryoNet can resolve the mechanism of action of pharmaceutical substances. As part of this work, we freely provide more than 2 million images that were used to train and test EmbryoNet.

Genetically engineered zebrafish as models of skeletal development and regeneration

Zebrafish (Danio rerio) are aquatic vertebrates with significant homology to their terrestrial counterparts. While zebrafish have a centuries-long track record in developmental and regenerative biology, their utility has grown exponentially with the onset of modern genetics. This is exemplified in studies focused on skeletal development and repair. Herein, the numerous contributions of zebrafish to our understanding of the basic science of cartilage, bone, tendon/ligament, and other skeletal tissues are described, with a particular focus on applications to development and regeneration. We summarize the genetic strengths that have made the zebrafish a powerful model to understand skeletal biology. We also highlight the large body of existing tools and techniques available to understand skeletal development and repair in the zebrafish and introduce emerging methods that will aid in novel discoveries in skeletal biology. Finally, we review the unique contributions of zebrafish to our understanding of regeneration and highlight diverse routes of repair in different contexts of injury. We conclude that zebrafish will continue to fill a niche of increasing breadth and depth in the study of basic cellular mechanisms of skeletal biology.

Murine Animal Models in Osteogenesis Imperfecta: The Quest for Improving the Quality of Life

Osteogenesis imperfecta is a rare genetic disorder characterized by bone fragility, due to alterations in the type I collagen molecule. It is a very heterogeneous disease, both genetically and phenotypically, with a high variability of clinical phenotypes, ranging from mild to severe forms, the most extreme cases being perinatal lethal. There is no curative treatment for OI, and so great efforts are being made in order to develop effective therapies. In these attempts, the in vivo preclinical studies are of paramount importance; therefore, serious analysis is required to choose the right murine OI model able to emulate as closely as possible the disease of the target OI population. In this review, we summarize the features of OI murine models that have been used for preclinical studies until today, together with recently developed new murine models. The bone parameters that are usually evaluated in order to determine the relevance of new developing therapies are exposed, and finally, current and innovative therapeutic strategies attempts considered in murine OI models, along with their mechanism of action, are reviewed. This review aims to summarize the in vivo studies developed in murine models available in the field of OI to date, in order to help the scientific community choose the most accurate OI murine model when developing new therapeutic strategies capable of improving the quality of life.

Making sense of fragmentation and merging in lineage tracing experiments

Lineage tracing experiments give dynamic information on the functional behaviour of dividing cells. These experiments therefore have become an important tool for studying stem and progenitor cell fate behavior in vivo. When cell proliferation is high or the frequency of induced clones cannot be precisely controlled, the merging and fragmentation of clones renders the retrospective interpretation of clonal fate data highly ambiguous, potentially leading to unguarded interpretations about lineage relationships and fate behaviour. Here, we discuss and generalize statistical strategies to detect, resolve and make use of clonal fragmentation and merging. We first explain how to detect the rates of clonal fragmentation and merging using simple statistical estimates. We then discuss ways to restore the clonal provenance of labelled cells algorithmically and statistically and elaborate on how the process of clonal fragmentation can indirectly inform about cell fate. We generalize and extend results from the context of their original publication.

Phenotypic screening models for rapid diagnosis of genetic variants and discovery of personalized therapeutics

Precision medicine strives for highly individualized treatments for disease under the notion that each individual's unique genetic makeup and environmental exposures imprints upon them not only a disposition to illness, but also an optimal therapeutic approach. In the realm of rare disorders, genetic predisposition is often the predominant mechanism driving disease presentation. For such, mostly, monogenic disorders, a causal gene to phenotype association is likely. As a result, it becomes important to query the patient's genome for the presence of pathogenic variations that are likely to cause the disease. Determining whether a variant is pathogenic or not is critical to these analyses and can be challenging, as many disease-causing variants are novel and, ergo, have no available functional data to help categorize them. This problem is exacerbated by the need for rapid evaluation of pathogenicity, since many genetic diseases present in young children who will experience increased morbidity and mortality without rapid diagnosis and therapeutics. Here, we discuss the utility of animal models, with a focus mainly on C. elegans, as a contrast to tissue culture and in silico approaches, with emphasis on how these systems are used in determining pathogenicity of variants with uncertain significance and then used to screen for novel therapeutics.

wnt16 regulates spine and muscle morphogenesis through parallel signals from notochord and dermomyotome

Bone and muscle are coupled through developmental, mechanical, paracrine, and autocrine signals. Genetic variants at the CPED1-WNT16 locus are dually associated with bone- and muscle-related traits. While Wnt16 is necessary for bone mass and strength, this fails to explain pleiotropy at this locus. Here, we show wnt16 is required for spine and muscle morphogenesis in zebrafish. In embryos, wnt16 is expressed in dermomyotome and developing notochord, and contributes to larval myotome morphology and notochord elongation. Later, wnt16 is expressed at the ventral midline of the notochord sheath, and contributes to spine mineralization and osteoblast recruitment. Morphological changes in wnt16 mutant larvae are mirrored in adults, indicating that wnt16 impacts bone and muscle morphology throughout the lifespan. Finally, we show that wnt16 is a gene of major effect on lean mass at the CPED1-WNT16 locus. Our findings indicate that Wnt16 is secreted in structures adjacent to developing bone (notochord) and muscle (dermomyotome) where it affects the morphogenesis of each tissue, thereby rendering wnt16 expression into dual effects on bone and muscle morphology. This work expands our understanding of wnt16 in musculoskeletal development and supports the potential for variants to act through WNT16 to influence bone and muscle via parallel morphogenetic processes.

In humans, genetic variants (DNA sequences that vary amongst individuals) have been identified that appear to influence two tissues, bone and skeletal muscle. However, how single genes and genetic variants exert dual influence on both tissues is not well understood. In this study, we found that the wnt16 gene is necessary for specifying the size and shape of both muscle and bone during development in zebrafish. We also disentangled how wnt16 affects both tissues: distinct cellular populations adjacent to muscle and bone secrete Wnt16, where it acts as a signal guiding the size and shape of each tissue. This is important because in humans, genetic variants near the WNT16 gene have effects on both bone- and muscle-related traits. This study expands our understanding of the role of WNT16 in bone and muscle development, and helps to explain how genetic variants near WNT16 affect traits for both tissues. Moreover, WNT16 is actively being explored as a target for osteoporosis therapies, thus our study could have implications with regard to the potential of targeting WNT16 to treat bone and muscle simultaneously.

Repositioning of Lansoprazole as a Protective Agent Against Cisplatin-Induced Ototoxicity

Cisplatin (CDDP) is a well-known chemotherapeutic drug approved for various cancers. However, CDDP accumulates in the inner ear cochlea via organic cation transporter 2 (OCT2) and causes ototoxicity, which is a major clinical limitation. Since lansoprazole (LPZ), a proton pump inhibitor, is known to inhibit OCT2-mediated transport of CDDP, we hypothesized that LPZ might ameliorate CDDP-induced ototoxicity (CIO). To test this hypothesis, we utilized in vivo fluorescence imaging of zebrafish sensory hair cells. The fluorescence signals in hair cells in zebrafish treated with CDDP dose-dependently decreased. Co-treatment with LPZ significantly suppressed the decrease of fluorescence signals in zebrafish treated with CDDP. Knockout of a zebrafish homolog of OCT2 also ameliorated the reduction of fluorescence signals in hair cells in zebrafish treated with CDDP. These in vivo studies suggest that CDDP damages the hair cells of zebrafish through oct2-mediated accumulation and that LPZ protects against CIO, possibly inhibiting the entry of CDDP into the hair cells via oct2. We also evaluated the otoprotective effect of LPZ using a public database containing adverse event reports. The analysis revealed that the incidence rate of CIO was significantly decreased in patients treated with LPZ. We then retrospectively analyzed the medical records of Mie University Hospital to examine the otoprotective effect of LPZ. The incidence rate of ototoxicity was significantly lower in patients co-treated with LPZ compared to those without LPZ. These retrospective findings suggest that LPZ is also protective against CIO in humans. Taken together, co-treatment with LPZ may reduce the risk of CIO.

Partial fads2 Gene Knockout Diverts LC-PUFA Biosynthesis via an Alternative Δ8 Pathway with an Impact on the Reproduction of Female Zebrafish (Danio rerio)

The zebrafish (Danio rerio) genome contains a single gene fads2 encoding a desaturase (FADS2) with both Δ6 and Δ5 activities, the key player in the endogenous biosynthesis of long-chain polyunsaturated fatty acids (LC-PUFAs), which serve essential functions as membrane components, sources of energy and signaling molecules. LC-PUFAs include the precursors of eicosanoids and are thus predicted to be indispensable molecules for reproductive health in virtually all vertebrates. In mice, an amniotic vertebrate, fads2 deletion mutants, both males and females, have been confirmed to be sterile. In anamniotic vertebrates, such as fish, there is still no information available on the reproductive (in)ability of fads2 mutants, although zebrafish have become an increasingly important model of lipid metabolism, including some aspects of the generation of germ cells and early embryonic development. In the present study, we apply the CRISPR/Cas9 genome editing system to induce mutations in the zebrafish genome and create crispants displaying a degree of fads2 gene editing within the range of 50–80%. Focusing on adult G₀ crispant females, we investigated the LC-PUFA profiles of eggs. Our data suggest an impaired pathway of the LC-PUFA biosynthesis of the ω6 and ω3 series in the first-rate limiting steps of the conversion of linoleic acid (LA) into γ-linolenic acid (GLA), and α-linolenic acid (ALA) into stearidonic acid (SDA), respectively, finally resulting in bad-quality eggs. Our data suggest the existence of an alternative Δ8 pathway, which bypasses the first endogenous LC-PUFA biosynthetic step in zebrafish in vivo, and suggest that the zebrafish bifunctional FADS2 enzyme is actually a trifunctional Δ6/Δ5/Δ8 desaturase.

Computational anatomy and geometric shape analysis enables analysis of complex craniofacial phenotypes in zebrafish

Due to the complexity of fish skulls, previous attempts to classify craniofacial phenotypes have relied on qualitative features or sparce 2D landmarks. In this work we aim to identify previously unknown 3D craniofacial phenotypes with a semiautomated pipeline in adult zebrafish mutants. We first estimate a synthetic 'normative' zebrafish template using MicroCT scans from a sample pool of wild-type animals using the Advanced Normalization Tools (ANTs). We apply a computational anatomy (CA) approach to quantify the phenotype of zebrafish with disruptions in bmp1a, a gene implicated in later skeletal development and whose human ortholog when disrupted is associated with Osteogenesis Imperfecta. Compared to controls, the bmp1a fish have larger otoliths, larger normalized centroid sizes, and exhibit shape differences concentrated around the operculum, anterior frontal, and posterior parietal bones. Moreover, bmp1a fish differ in the degree of asymmetry. Our CA approach offers a potential pipeline for high-throughput screening of complex fish craniofacial shape to discover novel phenotypes for which traditional landmarks are too sparce to detect. The current pipeline successfully identifies areas of variation in zebrafish mutants, which are an important model system for testing genome to phenome relationships in the study of development, evolution, and human diseases. This article has an associated First Person interview with the first author of the paper.

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.

Evaluation of CRISPR gene-editing tools in zebrafish

BACKGROUND

Zebrafish have practical features that make them a useful model for higher-throughput tests of gene function using CRISPR/Cas9 editing to create 'knockout' models. In particular, the use of G0 mosaic mutants has potential to increase throughput of functional studies significantly but may suffer from transient effects of introducing Cas9 via microinjection. Further, a large number of computational and empirical tools exist to design CRISPR assays but often produce varied predictions across methods leaving uncertainty in choosing an optimal approach for zebrafish studies.

METHODS

To systematically assess accuracy of tool predictions of on- and off-target gene editing, we subjected zebrafish embryos to CRISPR/Cas9 with 50 different guide RNAs (gRNAs) targeting 14 genes. We also investigate potential confounders of G0-based CRISPR screens by assaying control embryos for spurious mutations and altered gene expression.

RESULTS

We compared our experimental in vivo editing efficiencies in mosaic G0 embryos with those predicted by eight commonly used gRNA design tools and found large discrepancies between methods. Assessing off-target mutations (predicted in silico and in vitro) found that the majority of tested loci had low in vivo frequencies (< 1%). To characterize if commonly used 'mock' CRISPR controls (larvae injected with Cas9 enzyme or mRNA with no gRNA) exhibited spurious molecular features that might exacerbate studies of G0 mosaic CRISPR knockout fish, we generated an RNA-seq dataset of various control larvae at 5 days post fertilization. While we found no evidence of spontaneous somatic mutations of injected larvae, we did identify several hundred differentially-expressed genes with high variability between injection types. Network analyses of shared differentially-expressed genes in the 'mock' injected larvae implicated a number of key regulators of common metabolic pathways, and gene-ontology analysis revealed connections with response to wounding and cytoskeleton organization, highlighting a potentially lasting effect from the microinjection process that requires further investigation.

CONCLUSION

Overall, our results provide a valuable resource for the zebrafish community for the design and execution of CRISPR/Cas9 experiments.

The genetic overlap between osteoporosis and craniosynostosis

Osteoporosis is the most prevalent bone condition in the ageing population. This systemic disease is characterized by microarchitectural deterioration of bone, leading to increased fracture risk. In the past 15 years, genome-wide association studies (GWAS), have pinpointed hundreds of loci associated with bone mineral density (BMD), helping elucidate the underlying molecular mechanisms and genetic architecture of fracture risk. However, the challenge remains in pinpointing causative genes driving GWAS signals as a pivotal step to drawing the translational therapeutic roadmap. Recently, a skull BMD-GWAS uncovered an intriguing intersection with craniosynostosis, a congenital anomaly due to premature suture fusion in the skull. Here, we recapitulate the genetic contribution to both osteoporosis and craniosynostosis, describing the biological underpinnings of this overlap and using zebrafish models to leverage the functional investigation of genes associated with skull development and systemic skeletal homeostasis.

Perspective of the GEMSTONE Consortium on Current and Future Approaches to Functional Validation for Skeletal Genetic Disease Using Cellular, Molecular and Animal-Modeling Techniques

The availability of large human datasets for genome-wide association studies (GWAS) and the advancement of sequencing technologies have boosted the identification of genetic variants in complex and rare diseases in the skeletal field. Yet, interpreting results from human association studies remains a challenge. To bridge the gap between genetic association and causality, a systematic functional investigation is necessary. Multiple unknowns exist for putative causal genes, including cellular localization of the molecular function. Intermediate traits ("endophenotypes"), e.g. molecular quantitative trait loci (molQTLs), are needed to identify mechanisms of underlying associations. Furthermore, index variants often reside in non-coding regions of the genome, therefore challenging for interpretation. Knowledge of non-coding variance (e.g. ncRNAs), repetitive sequences, and regulatory interactions between enhancers and their target genes is central for understanding causal genes in skeletal conditions. Animal models with deep skeletal phenotyping and cell culture models have already facilitated fine mapping of some association signals, elucidated gene mechanisms, and revealed disease-relevant biology. However, to accelerate research towards bridging the current gap between association and causality in skeletal diseases, alternative in vivo platforms need to be used and developed in parallel with the current -omics and traditional in vivo resources. Therefore, we argue that as a field we need to establish resource-sharing standards to collectively address complex research questions. These standards will promote data integration from various -omics technologies and functional dissection of human complex traits. In this mission statement, we review the current available resources and as a group propose a consensus to facilitate resource sharing using existing and future resources. Such coordination efforts will maximize the acquisition of knowledge from different approaches and thus reduce redundancy and duplication of resources. These measures will help to understand the pathogenesis of osteoporosis and other skeletal diseases towards defining new and more efficient therapeutic targets.

Integration of multiple imaging platforms to uncover cardiovascular defects in adult zebrafish

Aims

Mammalian models have been instrumental in investigating adult heart function and human disease. However, electrophysiological differences with human hearts and high costs motivate the need for non-mammalian models. The zebrafish is a well-established genetic model to study cardiovascular development and function; however, analysis of cardiovascular phenotypes in adult specimens is particularly challenging as they are opaque.

Methods and results

Here, we optimized and combined multiple imaging techniques including echocardiography, magnetic resonance imaging, and micro-computed tomography to identify and analyse cardiovascular phenotypes in adult zebrafish. Using alk5a/tgfbr1a mutants as a case study, we observed morphological and functional cardiovascular defects that were undetected with conventional approaches. Correlation analysis of multiple parameters revealed an association between haemodynamic defects and structural alterations of the heart, as observed clinically.

Conclusion

We report a new, comprehensive, and sensitive platform to identify otherwise indiscernible cardiovascular phenotypes in adult zebrafish.

Lrp5 Mutant and Crispant Zebrafish Faithfully Model Human Osteoporosis, Establishing the Zebrafish as a Platform for CRISPR-Based Functional Screening of Osteoporosis Candidate Genes

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 mass-spectrometry 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. © 2021 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

3D assessment of intervertebral disc degeneration in zebrafish identifies changes in bone density that prime disc disease

Back pain is a common condition with a high social impact and represents a global health burden. Intervertebral disc disease (IVDD) is one of the major causes of back pain; no therapeutics are currently available to reverse this disease. The impact of bone mineral density (BMD) on IVDD has been controversial, with some studies suggesting osteoporosis as causative for IVDD and others suggesting it as protective for IVDD. Functional studies to evaluate the influence of genetic components of BMD in IVDD could highlight opportunities for drug development and repurposing. By taking a holistic 3D approach, we established an aging zebrafish model for spontaneous IVDD. Increased BMD in aging, detected by automated computational analysis, is caused by bone deformities at the endplates. However, aged zebrafish spines showed changes in bone morphology, microstructure, mineral heterogeneity, and increased fragility that resembled osteoporosis. Elements of the discs recapitulated IVDD symptoms found in humans: the intervertebral ligament (equivalent to the annulus fibrosus) showed disorganized collagen fibers and herniation, while the disc center (nucleus pulposus equivalent) showed dehydration and cellular abnormalities. We manipulated BMD in young zebrafish by mutating sp7 and cathepsin K, leading to low and high BMD, respectively. Remarkably, we detected IVDD in both groups, demonstrating that low BMD does not protect against IVDD, and we found a strong correlation between high BMD and IVDD. Deep learning was applied to high-resolution synchrotron µCT image data to analyze osteocyte 3D lacunar distribution and morphology, revealing a role of sp7 in controlling the osteocyte lacunar 3D profile. Our findings suggest potential avenues through which bone quality can be targeted to identify beneficial therapeutics for IVDD.

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).

Opportunities and Challenges in Functional Genomics Research in Osteoporosis: Report From a Workshop Held by the Causes Working Group of the Osteoporosis and Bone Research Academy of the Royal Osteoporosis Society on October 5th 2020

The discovery that sclerostin is the defective protein underlying the rare heritable bone mass disorder, sclerosteosis, ultimately led to development of anti-sclerostin antibodies as a new treatment for osteoporosis. In the era of large scale GWAS, many additional genetic signals associated with bone mass and related traits have since been reported. However, how best to interrogate these signals in order to identify the underlying gene responsible for these genetic associations, a prerequisite for identifying drug targets for further treatments, remains a challenge. The resources available for supporting functional genomics research continues to expand, exemplified by "multi-omics" database resources, with improved availability of datasets derived from bone tissues. These databases provide information about potential molecular mediators such as mRNA expression, protein expression, and DNA methylation levels, which can be interrogated to map genetic signals to specific genes based on identification of causal pathways between the genetic signal and the phenotype being studied. Functional evaluation of potential causative genes has been facilitated by characterization of the "osteocyte signature", by broad phenotyping of knockout mice with deletions of over 7,000 genes, in which more detailed skeletal phenotyping is currently being undertaken, and by development of zebrafish as a highly efficient additional in vivo model for functional studies of the skeleton. Looking to the future, this expanding repertoire of tools offers the hope of accurately defining the major genetic signals which contribute to osteoporosis. This may in turn lead to the identification of additional therapeutic targets, and ultimately new treatments for osteoporosis.

A simple and effective F0 knockout method for rapid screening of behaviour and other complex phenotypes

Hundreds of human genes are associated with neurological diseases, but translation into tractable biological mechanisms is lagging. Larval zebrafish are an attractive model to investigate genetic contributions to neurological diseases. However, current CRISPR-Cas9 methods are difficult to apply to large genetic screens studying behavioural phenotypes. To facilitate rapid genetic screening, we developed a simple sequencing-free tool to validate gRNAs and a highly effective CRISPR-Cas9 method capable of converting >90% of injected embryos directly into F0 biallelic knockouts. We demonstrate that F0 knockouts reliably recapitulate complex mutant phenotypes, such as altered molecular rhythms of the circadian clock, escape responses to irritants, and multi-parameter day-night locomotor behaviours. The technique is sufficiently robust to knockout multiple genes in the same animal, for example to create the transparent triple knockout crystal fish for imaging. Our F0 knockout method cuts the experimental time from gene to behavioural phenotype in zebrafish from months to one week.

Rapid and efficient cataract gene evaluation in F0 zebrafish using CRISPR-Cas9 ribonucleoprotein complexes

Cataract is the leading cause of blindness worldwide. Congenital or paediatric cataract can result in permanent visual impairment or blindness even with best attempts at treatment. A significant proportion of paediatric cataract has a genetic cause. Therefore, identifying the genes that lead to cataract formation is essential for understanding the pathological process of inherited paediatric cataract as well as to the development of new therapies. Despite clear progress in genomics technologies, verification of the biological effects of newly identified candidate genes and variants is still challenging. Here, we provide a step-by-step pipeline to evaluate cataract candidate genes in F0 zebrafish using CRISPR-Cas9 ribonucleoprotein complexes (RNP). Detailed descriptions of CRISPR-Cas9 RNP design and formulation, microinjection, optimization of CRISPR-Cas9 RNP reagent dose and delivery route, editing efficacy analysis as well as cataract formation evaluation are included. Following this protocol, any cataract candidates can be readily and efficiently evaluated within 2 weeks using basic laboratory supplies.

Zebrafish: A Resourceful Vertebrate Model to Investigate Skeletal Disorders

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