MicroCT-Based Phenomics in the Zebrafish Skeleton Reveals Virtues of Deep Phenotyping in a Distributed Organ System

Quantitative Geometric Modeling of Blood Cells from X-ray Histotomograms of Whole Zebrafish Larvae

Tissue phenotyping is foundational to understanding and assessing the cellular aspects of disease in organismal context and an important adjunct to molecular studies in the dissection of gene function, chemical effects, and disease. As a first step toward computational tissue phenotyping, we explore the potential of cellular phenotyping from 3-Dimensional (3D), 0.74 µm isotropic voxel resolution, whole zebrafish larval images derived from X-ray histotomography, a form of micro-CT customized for histopathology. As proof of principle towards computational tissue phenotyping of cells, we created a semi-automated mechanism for the segmentation of blood cells in the vascular spaces of zebrafish larvae, followed by modeling and extraction of quantitative geometric parameters. Manually segmented cells were used to train a random forest classifier for blood cells, enabling the use of a generalized cellular segmentation algorithm for the accurate segmentation of blood cells. These models were used to create an automated data segmentation and analysis pipeline to guide the steps in a 3D workflow including blood cell region prediction, cell boundary extraction, and statistical characterization of 3D geometric and cytological features. We were able to distinguish blood cells at two stages in development (4- and 5-days-post-fertilization) and wild-type vs. polA2 huli hutu ( hht ) mutants. The application of geometric modeling across cell types to and across organisms and sample types may comprise a valuable foundation for computational phenotyping that is more open, informative, rapid, objective, and reproducible.

A wide-field micro-computed tomography detector: micron resolution at half-centimetre scale

Ideal three-dimensional imaging of complex samples made up of micron-scale structures extending over mm to cm, such as biological tissues, requires both wide field of view and high resolution. For existing optics and detectors used for micro-CT (computed tomography) imaging, sub-micron pixel resolution can only be achieved for fields of view of <2 mm. This article presents a unique detector system with a 6 mm field-of-view image circle and 0.5 µm pixel size that can be used in micro-CT units utilizing both synchrotron and commercial X-ray sources. A resolution-test pattern with linear microstructures and whole adult Daphnia magna were imaged at beamline 8.3.2 of the Berkeley Advanced Light Source. Volumes of 10000 × 10000 × 7096 isotropic 0.5 µm voxels were reconstructed over a 5.0 mm × 3.5 mm field of view. Measurements in the projection domain confirmed a 0.90 µm measured spatial resolution that is largely Nyquist-limited. This unprecedented combination of field of view and resolution dramatically reduces the need for sectional scans and computational stitching for large samples, ultimately offering the means to elucidate changes in tissue and cellular morphology in the context of larger, whole, intact model organisms and specimens. This system is also anticipated to benefit micro-CT imaging in materials science, microelectronics, agricultural science and biomedical engineering.

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

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.

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.

Quantitative analyses of adiposity dynamics in zebrafish

Adipose tissues often exhibit subtle, quantitative differences between individuals, leading to a graded series of adiposity phenotypes at the population level. Robust, quantitative analyses are vital for studying these differences. In this Commentary we highlight two articles from our lab that employ sensitive new methods in zebrafish capable of delineating complex and quantitative adiposity phenotypes. In the first article, we utilized in vivo imaging to systematically quantify zebrafish adipose tissues. We identified 34 regionally distinct zebrafish adipose tissues and developed statistical models to predict the size and variance of each adipose tissue over the course of zebrafish growth. We then employed these models to identify effects of strain and diet on adipose tissue growth. In the second article, we employed deep phenotyping to study complex disease-related adiposity traits. Using this methodology, we identified that adipose tissues have unique capacities to re-deposit lipid following food restriction and re-feeding. These distinct re-deposition potentials led to widespread fat distribution changes following re-feeding. We discuss how these novel findings may provide relevance to health conditions such as anorexia nervosa. Together, the strategies described in these two articles can be used as unbiased and quantitative methods to uncover new relationships between genotype, diet and adiposity.

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.

Deep phenotyping in zebrafish reveals genetic and diet-induced adiposity changes that may inform disease risk

The regional distribution of adipose tissues is implicated in a wide range of diseases. For example, proportional increases in visceral adipose tissue increase the risk for insulin resistance, diabetes, and CVD. Zebrafish offer a tractable model system by which to obtain unbiased and quantitative phenotypic information on regional adiposity, and deep phenotyping can explore complex disease-related adiposity traits. To facilitate deep phenotyping of zebrafish adiposity traits, we used pairwise correlations between 67 adiposity traits to generate stage-specific adiposity profiles that describe changing adiposity patterns and relationships during growth. Linear discriminant analysis classified individual fish according to an adiposity profile with 87.5% accuracy. Deep phenotyping of eight previously uncharacterized zebrafish mutants identified neuropilin 2b as a novel gene that alters adipose distribution. When we applied deep phenotyping to identify changes in adiposity during diet manipulations, zebrafish that underwent food restriction and refeeding had widespread adiposity changes when compared with continuously fed, equivalently sized control animals. In particular, internal adipose tissues (e.g., visceral adipose) exhibited a reduced capacity to replenish lipid following food restriction. Together, these results in zebrafish establish a new deep phenotyping technique as an unbiased and quantitative method to help uncover new relationships between genotype, diet, and adiposity.