Transgene-mediated skeletal phenotypic variation in zebrafish

Loss of developmentally derived Irf8+ macrophages promotes hyperinnervation and arrhythmia in the adult zebrafish heart

Recent developments in cardiac macrophage biology have broadened our understanding of the critical functions of macrophages in the heart. As a result, there is further interest in understanding the independent contributions of distinct subsets of macrophage to cardiac development and function. Here, we demonstrate that genetic loss of interferon regulatory factor 8 (Irf8)-positive embryonic-derived macrophages significantly disrupts cardiac conduction, chamber function, and innervation in adult zebrafish. At 4 months post-fertilization (mpf), homozygous irf8st96/st96 mutants have significantly shortened atrial action potential duration and significant differential expression of genes involved in cardiac contraction. Functional in vivo assessments via electro- and echocardiograms at 12 mpf reveal that irf8 mutants are arrhythmogenic and exhibit diastolic dysfunction and ventricular stiffening. To identify the molecular drivers of the functional disturbances in irf8 null zebrafish, we perform single cell RNA sequencing and immunohistochemistry, which reveal increased leukocyte infiltration, epicardial activation, mesenchymal gene expression, and fibrosis. Irf8 null hearts are also hyperinnervated and have aberrant axonal patterning, a phenotype not previously assessed in the context of cardiac macrophage loss. Gene ontology analysis supports a novel role for activated epicardial-derived cells (EPDCs) in promoting neurogenesis and neuronal remodeling in vivo. Together, these data uncover significant cardiac abnormalities following embryonic macrophage loss and expand our knowledge of critical macrophage functions in heart physiology and governing homeostatic heart health.

The Fraser complex interconnects tissue layers to support basal epidermis and osteoblast integrated morphogenesis underlying fin skeletal patterning

Fraser Syndrome is a rare, multisystemic autosomal recessive disorder characterized by disrupted epithelial-mesenchymal associations upon loss of Fraser Complex genes. Disease manifestation and affected organs are highly variable. Digit malformations such as syndactyly are common but of unclear developmental origins. We explored if zebrafish fraser extracellular matrix complex subunit 1 (fras1) mutants model Fraser Syndrome-associated appendicular skeleton patterning defects. Approximately 10% of fras1 mutants survive to adulthood, displaying striking and varied fin abnormalities, including endochondral bone fusions, ectopic cartilage, and disrupted caudal fin symmetry. The fins of surviving fras1 mutants frequently have fewer and unbranched bony rays. fras1 mutant fins regenerate to their original size but with exacerbated ray branching and fin symmetry defects. Single cell RNA-Seq analysis, in situ hybridizations, and antibody staining show specific Fraser complex expression in the basal epidermis during regenerative outgrowth. Fras1 and Fraser Complex component Frem2 accumulate along the basal side of distal-most basal epidermal cells. Greatly reduced and mislocalized Frem2 accompanies loss of Fras1 in fras1 mutants. The Sonic hedgehog signaling between distal basal epidermis and adjacent mesenchymal pre-osteoblasts that promotes ray branching persists upon Fraser Complex loss. However, fras1 mutant regenerating fins exhibit extensive sub-epidermal blistering associated with a disorganized basal epidermis and adjacent pre-osteoblasts. We propose Fraser Complex-supported tissue layer adhesion enables robust integrated tissue morphogenesis involving the basal epidermis and osteoblasts. Further, we establish zebrafish fin development and regeneration as an accessible model to explore mechanisms of Fraser Syndrome-associated digit defects and Fraser Complex function at epithelial-mesenchymal interfaces.

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

The alx3 gene shapes the zebrafish neurocranium by regulating frontonasal neural crest cell differentiation timing

During craniofacial development, different populations of cartilage- and bone-forming cells develop in precise locations in the head. Most of these cells are derived from pluripotent cranial neural crest cells and differentiate with distinct developmental timing and cellular morphologies. The mechanisms that divide neural crest cells into discrete populations are not fully understood. Here, we use single-cell RNA sequencing to transcriptomically define different populations of cranial neural crest cells. We discovered that the gene family encoding the Alx transcription factors is enriched in the frontonasal population of neural crest cells. Genetic mutant analyses indicate that alx3 functions to regulate the distinct differentiation timing and cellular morphologies among frontonasal neural crest cell subpopulations. This study furthers our understanding of how genes controlling developmental timing shape craniofacial skeletal elements.