A vitamin D enriched diet attenuates sex-specific behavioral deficits, increases the lifespan, but does not rescue bone abnormalities in a mouse model of cortical dysplasia.

Individuals who experience recurrent spontaneous seizures often show behavioral and physiological comorbidities. Those with epilepsy are at a high risk of bone fractures (independent of seizure-related falls) and show a higher rate of a diagnosis of Autism Spectrum Disorder. The neural subset-specific (NS) Pten knockout (KO) mouse has an epilepsy phenotype, has been characterized to show autistic-like deficits, and has an osteoporosis phenotype. The current study examined the effect of a vitamin D enriched diet (20,000 IU VD) in the NS-Pten KO and wildtype mice. Mice were placed onto a vitamin D enriched diet at 4 weeks of age and maintained on that diet throughout testing. Behavioral testing began at 6 weeks of age and included tests for general activity, anxiety, repetitive behaviors, social behaviors, and memory. Results indicated that a vitamin D diet attenuated hypoactivity levels in male KO mice (p < 0.05). In a social partition task, vitamin D increased sociability in male wildtype mice, (p < 0.05). Most significantly, vitamin D fortified diet increased percent survival in KO animals and decreased the level of microglia marker IBA-1 and mTOR (mammalian target of rapamycin) downstream targets pS6 and pAKT. A high vitamin D diet did not reverse bone deficits in male or female KO mice. Overall, these findings suggest that a vitamin D enriched diet had a significant impact on the behavioral phenotype of NS-Pten KO mice, suggesting that dietary manipulations could be a potential therapeutic option for autistic-like behavior.

Zebrafish type I collagen mutants faithfully recapitulate human type I collagenopathies.

The type I collagenopathies are a group of heterogeneous connective tissue disorders, that are caused by mutations in the genes encoding type I collagen and include specific forms of osteogenesis imperfecta (OI) and the Ehlers-Danlos syndrome (EDS). These disorders present with a broad disease spectrum and large clinical variability of which the underlying genetic basis is still poorly understood. In this study, we systematically analyzed skeletal phenotypes in a large set of zebrafish, with diverse mutations in the genes encoding type I collagen, representing different genetic forms of human OI, and a zebrafish model resembling human EDS, which harbors a number of soft connective tissues defects, typical of EDS. Furthermore, we provide insight into how zebrafish and human type I collagen are compositionally and functionally related, which is relevant in the interpretation of human type I collagen-related disease models. Our studies reveal a high degree of intergenotype variability in phenotypic expressivity that closely correlates with associated OI severity. Furthermore, we demonstrate the potential for select mutations to give rise to phenotypic variability, mirroring the clinical variability associated with human disease pathology. Therefore, our work suggests the future potential for zebrafish to aid in identifying unknown genetic modifiers and mechanisms underlying the phenotypic variability in OI and related disorders. This will improve diagnostic strategies and enable the discovery of new targetable pathways for pharmacological intervention.

Neuromuscular dysfunction, independent of gait dysfunction, modulates trabecular bone homeostasis in mice.

OBJECTIVES: To clarify the effects of neuromuscular dysfunction on hindlimb loading, muscle atrophy, and bone homeostasis. METHODS: We quantified changes to hindlimb loading, muscle atrophy, and bone morphology following either Botulinum toxin A (BTxA) induced muscle paralysis or peripheral nerve injury (PNI) in mice; two in vivo models that we anticipated would differently alter gait and mechanical loading patterns due to their distinct effects on neuromuscular signaling. To confirm the expected behavioral effects of PNI, we assessed mechanical allodynia of the ipsilateral hindlimb using von Frey testing and activity (distance traveled and speed) was monitored in both groups using open field testing. Peak vertical ground reaction forces (GRF) and ankle and knee kinematics during normal locomotion were quantified and used to estimate peak mid-diaphyseal normal strains. Muscle atrophy and trabecular and cortical bone morphology were assessed via high-resolution microCT imaging. RESULTS: BTxA-induced calf paralysis caused severe muscle atrophy and altered gait kinetics and kinematics and reduced gait-induced normal strains. PNI increased mechanical allodynia but did not alter gait, nor did it cause muscle atrophy. We observed that muscle paralysis and PNI both led to severe trabecular bone loss but only BTxA-induced paralysis increased cortical bone resorption. CONCLUSIONS: While mechanical stimuli clearly have essential functions in bone development and adaptation, these data emphasize that neuromuscular signaling, independent of load-induced mechanical strains, may modulate trabecular bone homeostasis in normal and disease states.

Muscle paralysis induces bone marrow inflammation and predisposition to formation of giant osteoclasts.

Transient muscle paralysis engendered by a single injection of botulinum toxin A (BTxA) rapidly induces profound focal bone resorption within the medullary cavity of adjacent bones. While initially conceived as a model of mechanical disuse, osteoclastic resorption in this model is disproportionately severe compared with the modest gait defect that is created. Preliminary studies of bone marrow following muscle paralysis suggested acute upregulation of inflammatory cytokines, including TNF-α and IL-1. We therefore hypothesized that BTxA-induced muscle paralysis would rapidly alter the inflammatory microenvironment and the osteoclastic potential of bone marrow. We tested this hypothesis by defining the time course of inflammatory cell infiltration, osteoinflammatory cytokine expression, and alteration in osteoclastogenic potential in the tibia bone marrow following transient muscle paralysis of the calf muscles. Our findings identified inflammatory cell infiltration within 24 h of muscle paralysis. By 72 h, osteoclast fusion and pro-osteoclastic inflammatory gene expression were upregulated in tibia bone marrow. These alterations coincided with bone marrow becoming permissive to the formation of osteoclasts of greater size and greater nuclei numbers. Taken together, our data are consistent with the thesis that transient calf muscle paralysis induces acute inflammation within the marrow of the adjacent tibia and that these alterations are temporally consistent with a role in mediating muscle paralysis-induced bone resorption.

Human microglia show unique transcriptional changes in Alzheimer's disease.

Microglia, the innate immune cells of the brain, influence Alzheimer's disease (AD) progression and are potential therapeutic targets. However, microglia exhibit diverse functions, the regulation of which is not fully understood, complicating therapeutics development. To better define the transcriptomic phenotypes and gene regulatory networks associated with AD, we enriched for microglia nuclei from 12 AD and 10 control human dorsolateral prefrontal cortices (7 males and 15 females, all aged >60 years) before single-nucleus RNA sequencing. Here we describe both established and previously unrecognized microglial molecular phenotypes, the inferred gene networks driving observed transcriptomic change, and apply trajectory analysis to reveal the putative relationships between microglial phenotypes. We identify microglial phenotypes more prevalent in AD cases compared with controls. Further, we describe the heterogeneity in microglia subclusters expressing homeostatic markers. Our study demonstrates that deep profiling of microglia in human AD brain can provide insight into microglial transcriptional changes associated with AD.

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.

Single-cell resolution of MET- and EMT-like programs in osteoblasts during zebrafish fin regeneration.

Zebrafish regenerate fin rays following amputation through epimorphic regeneration, a process that has been proposed to involve the epithelial-to-mesenchymal transition (EMT). We performed single-cell RNA sequencing (scRNA-seq) to elucidate osteoblastic transcriptional programs during zebrafish caudal fin regeneration. We show that osteoprogenitors are enriched with components associated with EMT and its reverse, mesenchymal-to-epithelial transition (MET), and provide evidence that the EMT markers cdh11 and twist2 are co-expressed in dedifferentiating cells at the amputation stump at 1 dpa, and in differentiating osteoblastic cells in the regenerate, the latter of which are enriched in EMT signatures. We also show that esrp1, a regulator of alternative splicing in epithelial cells that is associated with MET, is expressed in a subset of osteoprogenitors during outgrowth. This study provides a single cell resource for the study of osteoblastic cells during zebrafish fin regeneration, and supports the contribution of MET- and EMT-associated components to this process.

Primary cilium-dependent mechanosensing is mediated by adenylyl cyclase 6 and cyclic AMP in bone cells.

Primary cilia are chemosensing and mechanosensing organelles that regulate remarkably diverse processes in a variety of cells. We previously showed that primary cilia play a role in mediating mechanosensing in bone cells through an unknown mechanism that does not involve extracellular Ca(2+)-dependent intracellular Ca(2+) release, which has been implicated in all other cells that transduce mechanical signals via the cilium. Here, we identify a molecular mechanism linking primary cilia and bone cell mechanotransduction that involves adenylyl cyclase 6 (AC6) and cAMP. Intracellular cAMP was quantified in MLO-Y4 cells exposed to dynamic flow, and AC6 and primary cilia were inhibited using RNA interference. When exposed to flow, cells rapidly (<2 min) and transiently decreased cAMP production in a primary cilium-dependent manner. RT-PCR revealed differential expression of the membrane-bound isoforms of adenylyl cyclase, while immunostaining revealed one, AC6, preferentially localized to the cilium. Further studies showed that decreases in cAMP in response to flow were dependent on AC6 and Gd(3+)-sensitive channels but not intracellular Ca(2+) release and that this response mediated flow-induced COX-2 gene expression. The signaling events identified provide important details of a novel early mechanosensing mechanism in bone and advances our understanding of how signal transduction occurs at the primary cilium.

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.

The Reissner Fiber Is Highly Dynamic In Vivo and Controls Morphogenesis of the Spine.

Cerebrospinal fluid (CSF) physiology is important for the development and homeostasis of the central nervous system, and its disruption has been linked to scoliosis in zebrafish [1, 2]. Suspended in the CSF is an extracellular structure called the Reissner fiber, which extends from the brain through the central canal of the spinal cord. Zebrafish scospondin-null mutants are unable to assemble a Reissner fiber and fail to form a straight body axis during embryonic development [3]. Here, we describe hypomorphic missense mutations of scospondin, which allow Reissner fiber assembly and extension of a straight axis. However, during larval development, these mutants display progressive Reissner fiber disassembly, which is concomitant with the emergence of axial curvatures and scoliosis in adult animals. Using a scospondin-GFP knockin zebrafish line, we demonstrate several dynamic properties of the Reissner fiber in vivo, including embryonic fiber assembly, the continuous rostral to caudal movement of the fiber within the brain and central canal, and subcommissural organ (SCO)-spondin-GFP protein secretion from the floor plate to merge with the fiber. Finally, we show that disassembly of the Reissner fiber is also associated with the progression of axial curvatures in distinct scoliosis mutant zebrafish models. Together, these data demonstrate a critical role for the Reissner fiber for the maintenance of a straight body axis and spine morphogenesis in adult zebrafish. Our study establishes a framework for future investigations to address the cellular effectors responsible for Reissner-fiber-dependent regulation of axial morphology. VIDEO ABSTRACT.

Phenomics-Based Quantification of CRISPR-Induced Mosaicism in Zebrafish.

Genetic mosaicism can manifest as spatially variable phenotypes that vary from site to site within an organism. Here, we use imaging-based phenomics to quantitate phenotypes at many sites within the axial skeleton of CRISPR-edited G0 zebrafish. Through characterization of loss-of-function cell clusters in the developing skeleton, we identify a distinctive size distribution shown to arise from clonal fragmentation and merger events. We quantitate the phenotypic mosaicism produced by somatic mutations of two genes, plod2 and bmp1a, implicated in human osteogenesis imperfecta. Comparison of somatic, CRISPR-generated G0 mutants to homozygous germline mutants reveals phenotypic convergence, suggesting that CRISPR screens of G0 animals can faithfully recapitulate the biology of inbred disease models. We describe statistical frameworks for phenomic analysis of spatial phenotypic variation present in somatic G0 mutants. In sum, this study defines an approach for decoding spatially variable phenotypes generated during CRISPR-based screens.

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.

Developmental tuning of mineralization drives morphological diversity of gill cover bones in sculpins and their relatives.

The role of osteoblast placement in skeletal morphological variation is relatively well understood, but alternative developmental mechanisms affecting bone shape remain largely unknown. Specifically, very little attention has been paid to variation in later mineralization stages of intramembranous ossification as a driver of morphological diversity. We discover the occurrence of specific, sometimes large, regions of nonmineralized osteoid within bones that also contain mineralized tissue. We show through a variety of histological, molecular, and tomographic tests that this "extended" osteoid material is most likely nonmineralized bone matrix. This tissue type is a significant determinant of gill cover bone shape in the teleostean suborder Cottoidei. We demonstrate repeated evolution of extended osteoid in Cottoidei through ancestral state reconstruction and test for an association between extended osteoid variation and habitat differences among species. Through measurement of extended osteoid at various stages of gill cover development in species across the phylogeny, we gain insight into possible evolutionary developmental origins of the trait. We conclude that this fine-tuned developmental regulation of bone matrix mineralization reflects heterochrony at multiple biological levels and is a novel mechanism for the evolution of diversity in skeletal morphology. This research lays the groundwork for a new model in which to study bone mineralization and evolutionary developmental processes, particularly as they may relate to adaptation during a prominent evolutionary radiation of fishes.

Osteocytes as mechanosensors in the inhibition of bone resorption due to mechanical loading.

Bone has the ability to adjust its structure to meet its mechanical environment. The prevailing view of bone mechanobiology is that osteocytes are responsible for detecting and responding to mechanical loading and initiating the bone adaptation process. However, how osteocytes signal effector cells and initiate bone turnover is not well understood. Recent in vitro studies have shown that osteocytes support osteoclast formation and activation when co-cultured with osteoclast precursors. In this study, we examined the osteocytes' role in the mechanical regulation of osteoclast formation and activation. We demonstrated here that (1) mechanical stimulation of MLO-Y4 osteocyte-like cells decreases their osteoclastogenic-support potential when co-cultured with RAW264.7 monocyte osteoclast precursors; (2) soluble factors released by these mechanically stimulated MLO-Y4 cells inhibit osteoclastogenesis induced by ST2 bone marrow stromal cells or MLO-Y4 cells; and (3) soluble RANKL and OPG were released by MLO-Y4 cells, and the expressions of both were found to be mechanically regulated. Our data suggest that mechanical loading decreases the osteocyte's potential to induce osteoclast formation by direct cell-cell contact. However, it is not clear that osteocytes in vivo are able to form contacts with osteoclast precursors. Our data also demonstrate that mechanically stimulated osteocytes release soluble factors that can inhibit osteoclastogenesis induced by other supporting cells including bone marrow stromal cells. In summary, we conclude that osteocytes may function as mechanotransducers by regulating local osteoclastogenesis via soluble signals.

ScreenCube: A 3D Printed System for Rapid and Cost-Effective Chemical Screening in Adult Zebrafish.

Phenotype-based small molecule screens in zebrafish embryos and larvae have been successful in accelerating pathway and therapeutic discovery for diverse biological processes. Yet, the application of chemical screens to adult physiologies has been relatively limited due to additional demands on cost, space, and labor associated with screens in adult animals. In this study, we present a 3D printed system and methods for intermittent drug dosing that enable rapid and cost-effective chemical administration in adult zebrafish. Using prefilled screening plates, the system enables dosing of 96 fish in ∼3 min, with a 10-fold reduction in drug quantity compared to that used in previous chemical screens in adult zebrafish. We characterize water quality kinetics during immersion in the system and use these kinetics to rationally design intermittent dosing regimens that result in 100% fish survival. As a demonstration of system fidelity, we show the potential to identify two known chemical inhibitors of adult tail fin regeneration, cyclopamine and dorsomorphin. By developing methods for rapid and cost-effective chemical administration in adult zebrafish, this study expands the potential for small molecule discovery in postembryonic models of development, disease, and regeneration.

Botulinum toxin A-induced muscle paralysis stimulates Hdac4 and differential miRNA expression.

At sufficient dose, intramuscular injection of Botulinum toxin A causes muscle wasting that is physiologically consistent with surgical denervation and other types of neuromuscular dysfunction. The aim of this study was to clarify early molecular and micro-RNA alterations in skeletal muscle following Botulinum toxin A-induced muscle paralysis. Quadriceps were analyzed for changes in expression of micro- and messenger RNA and protein levels after a single injection of 0.4, 2 or 4U Botulinum toxin A (/100g body weight). After injection with 2.0U Botulinum toxin A, quadriceps exhibited significant reduction in muscle weight and increased levels of ubiquitin ligase proteins at 7, 14 and 28 days. Muscle miR-1 and miR-133a/b levels were decreased at these time points, whereas a dose-responsive increase in miR-206 expression at day 14 was observed. Expression of the miR-133a/b target genes RhoA, Tgfb1 and Ctfg, and the miR-1/206 target genes Igf-1 and Hdac4, were upregulated at 28 days after Botulinum toxin A injection. Increased levels of Hdac4 protein were observed after injection, consistent with anticipated expression changes in direct and indirect Hdac4 target genes, such as Myog. Our results suggest Botulinum toxin A-induced denervation of muscle shares molecular characteristics with surgical denervation and other types of neuromuscular dysfunction, and implicates miR-133/Tgf-ß1/Ctfg and miR-1/Hdac4/Myog signaling during the resultant muscle atrophy.

Time-dependent deformations in bone cells exposed to fluid flow in vitro: investigating the role of cellular deformation in fluid flow-induced signaling.

Numerous experiments have shown fluid flow to be a potent stimulator of bone cells in vitro, suggesting that fluid flow is an important physical signal in bone mechanotransduction. In fluid flow experiments, bone cells are exposed to both time-dependent (e.g., oscillating or pulsing) and time-independent (e.g., steady) flow profiles. Interestingly, the signaling response of bone cells shows dependence on loading frequency and/or rate that has been postulated to be due to viscoelastic behavior. Thus, the objective of this study was to investigate the time-dependent deformations of bone cells exposed to fluid flow in vitro. Specifically, our goal was to characterize the mechanical response of bone cells exposed to oscillatory flow from 0.5 to 2.0 Hz and steady flow, since these flow profiles have previously been shown to induce different morphological and biochemical responses in vitro. By tracking cell-bound sulfate and collagen coated fluorescent beads of varying sizes, we quantified the normalized peak deformation (peak displacement normalized by the maximum peak displacement observed for all frequencies) and phase lag in bone cells exposed to 1.0 Pa oscillating flow at frequencies of 0.5-2.0 Hz. The phase lag was small (3-10 degrees ) and frequency dependent, while the normalized peak displacements decreased as a weak power law of frequency ( approximately f(-0.2)). During steady flow, the cells exhibited a nearly instantaneous deformation, followed by creep. Our results suggest that while substantial viscous deformation may occur during steady flow (compared to oscillating flow at approximately 1 Hz), bone cells behave primarily as elastic bodies when exposed to flow at frequencies associated with habitual loading.

A Dynamic Anesthesia System for Long-Term Imaging in Adult Zebrafish.

Long-term in vivo imaging in adult zebrafish (i.e., 1-24 h) has been limited by the fact that regimens for long-term anesthesia in embryos and larvae are ineffective in adults. Here, we examined the potential for dynamic administration of benzocaine to enable long-term anesthesia in adult zebrafish. We developed a computer-controlled perfusion system comprised of programmable peristaltic pumps that enabled automatic exchange between anesthetic and system water. Continuous administration of benzocaine in adult zebrafish resulted in a mean time to respiratory arrest of 5.0 h and 8-h survival of 14.3%. We measured characteristic sedation and recovery times in response to benzocaine, and used them to devise an intermittent dosing regimen consisting of 14.5 min of benzocaine followed by 5.5 min of system water. Intermittent benzocaine administration in adult zebrafish resulted in a mean time to respiratory arrest of 7.6 h and 8-h survival of 71.4%. Finally, we performed a single 24-h trial and found that intermittent dosing maintained anesthesia in an adult zebrafish over the entire 24-h period. In summary, our studies demonstrate the potential for dynamic administration of benzocaine to enable prolonged anesthesia in adult zebrafish, expanding the potential for imaging in adult physiologies that unfold over 1-24 h.

Neuron subset-specific Pten deletion induces abnormal skeletal activity in mice.

Individuals with a history of epilepsy are at higher risk for bone fractures compared to the general population. Although clinical studies support an association between low bone mineral density (BMD) and anti-seizure medications, little is known on whether a history of seizures is linked to altered bone health. Therefore, in this study we tested the hypothesis that bone mass, morphology, and bone mineralization are altered by seizures in genetically epileptic animals and in animals subjected to an episode of status epilepticus. In this study, we used NS-Pten conditional knockout mice (a well-studied genetic model of epilepsy). We used microCT analysis to measure BMD, morphology, and mineralization in NS-Pten+/+ (wildtype) and NS-Pten-/- (knockout) mice at 4 and 8weeks, as well as adult Kv4.2+/+ and Kv4.2-/- mice. We measured BMD, bone morphology, and mineralization in adult NS-Pten+/+ mice that received status epilepticus through kainic acid (20mg/kg intraperitoneal). Further, we measured locomotion for NS-Pten+/+ and NS-Pten-/- mice at 4 and 6weeks. We found that NS-Pten-/- mice exhibited low BMD in the tibial metaphysis and midshaft compared to non-epileptic mice. Morphologically, NS-Pten-/- mice exhibited decreased trabecular volume fraction, and endocortical expansion in both the metaphyeal and diaphyseal compartments. In the midshaft, NS-Pten-/- mice exhibited reduced tissue mineral density, indicating impaired mineralization in addition to morphological deficits. NS-Pten-/- mice exhibited hyperactivity in open field testing, suggesting low bone mass in NS-Pten-/- mice was not attributable to hypoactivity. Differences in BMD were not observed following kainate-induced seizures or in the Kv4.2-/- model of seizure susceptibility. Our findings suggest that deletion of Pten in the brain results in impaired bone mass and mineralization, which may contribute to weaker bones and thereby a higher fracture risk.

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.