Single-cell, whole-embryo phenotyping of mammalian developmental disorders.

Mouse models are a critical tool for studying human diseases, particularly developmental disorders1. However, conventional approaches for phenotyping may fail to detect subtle defects throughout the developing mouse2. Here we set out to establish single-cell RNA sequencing of the whole embryo as a scalable platform for the systematic phenotyping of mouse genetic models. We applied combinatorial indexing-based single-cell RNA sequencing3 to profile 101 embryos of 22 mutant and 4 wild-type genotypes at embryonic day 13.5, altogether profiling more than 1.6 million nuclei. The 22 mutants represent a range of anticipated phenotypic severities, from established multisystem disorders to deletions of individual regulatory regions4,5. We developed and applied several analytical frameworks for detecting differences in composition and/or gene expression across 52 cell types or trajectories. Some mutants exhibit changes in dozens of trajectories whereas others exhibit changes in only a few cell types. We also identify differences between widely used wild-type strains, compare phenotyping of gain- versus loss-of-function mutants and characterize deletions of topological associating domain boundaries. Notably, some changes are shared among mutants, suggesting that developmental pleiotropy might be 'decomposable' through further scaling of this approach. Overall, our findings show how single-cell profiling of whole embryos can enable the systematic molecular and cellular phenotypic characterization of mouse mutants with unprecedented breadth and resolution.

Embryo-scale reverse genetics at single-cell resolution.

The maturation of single-cell transcriptomic technologies has facilitated the generation of comprehensive cellular atlases from whole embryos1-4. A majority of these data, however, has been collected from wild-type embryos without an appreciation for the latent variation that is present in development. Here we present the 'zebrafish single-cell atlas of perturbed embryos': single-cell transcriptomic data from 1,812 individually resolved developing zebrafish embryos, encompassing 19 timepoints, 23 genetic perturbations and a total of 3.2 million cells. The high degree of replication in our study (eight or more embryos per condition) enables us to estimate the variance in cell type abundance organism-wide and to detect perturbation-dependent deviance in cell type composition relative to wild-type embryos. Our approach is sensitive to rare cell types, resolving developmental trajectories and genetic dependencies in the cranial ganglia neurons, a cell population that comprises less than 1% of the embryo. Additionally, time-series profiling of individual mutants identified a group of brachyury-independent cells with strikingly similar transcriptomes to notochord sheath cells, leading to new hypotheses about early origins of the skull. We anticipate that standardized collection of high-resolution, organism-scale single-cell data from large numbers of individual embryos will enable mapping of the genetic dependencies of zebrafish cell types, while also addressing longstanding challenges in developmental genetics, including the cellular and transcriptional plasticity underlying phenotypic diversity across individuals.

Proteostasis governs differential temperature sensitivity across embryonic cell types.

Embryonic development is remarkably robust, but temperature stress can degrade its ability to generate animals with invariant anatomy. Phenotypes associated with environmental stress suggest that some cell types are more sensitive to stress than others, but the basis of this sensitivity is unknown. Here, we characterize hundreds of individual zebrafish embryos under temperature stress using whole-animal single-cell RNA sequencing (RNA-seq) to identify cell types and molecular programs driving phenotypic variability. We find that temperature perturbs the normal proportions and gene expression programs of numerous cell types and also introduces asynchrony in developmental timing. The notochord is particularly sensitive to temperature, which we map to a specialized cell type: sheath cells. These cells accumulate misfolded protein at elevated temperature, leading to a cascading structural failure of the notochord and anatomic defects. Our study demonstrates that whole-animal single-cell RNA-seq can identify mechanisms for developmental robustness and pinpoint cell types that constitute key failure points.

Transcriptomic profiling of tissue environments critical for post-embryonic patterning and morphogenesis of zebrafish skin.

Pigment patterns and skin appendages are prominent features of vertebrate skin. In zebrafish, regularly patterned pigment stripes and an array of calcified scales form simultaneously in the skin during post-embryonic development. Understanding the mechanisms that regulate stripe patterning and scale morphogenesis may lead to the discovery of fundamental mechanisms that govern the development of animal form. To learn about cell types and signaling interactions that govern skin patterning and morphogenesis, we generated and analyzed single-cell transcriptomes of skin from wild-type fish as well as fish having genetic or transgenically induced defects in squamation or pigmentation. These data reveal a previously undescribed population of epidermal cells that express transcripts encoding enamel matrix proteins, suggest hormonal control of epithelial-mesenchymal signaling, clarify the signaling network that governs scale papillae development, and identify a critical role for the hypodermis in supporting pigment cell development. Additionally, these comprehensive single-cell transcriptomic data representing skin phenotypes of biomedical relevance should provide a useful resource for accelerating the discovery of mechanisms that govern skin development and homeostasis.

Single-cell transcriptomic profiling of the zebrafish inner ear reveals molecularly distinct hair cell and supporting cell subtypes.

A major cause of human deafness and vestibular dysfunction is permanent loss of the mechanosensory hair cells of the inner ear. In non-mammalian vertebrates such as zebrafish, regeneration of missing hair cells can occur throughout life. While a comparative approach has the potential to reveal the basis of such differential regenerative ability, the degree to which the inner ears of fish and mammals share common hair cells and supporting cell types remains unresolved. Here, we perform single-cell RNA sequencing of the zebrafish inner ear at embryonic through adult stages to catalog the diversity of hair cells and non-sensory supporting cells. We identify a putative progenitor population for hair cells and supporting cells, as well as distinct hair and supporting cell types in the maculae versus cristae. The hair cell and supporting cell types differ from those described for the lateral line system, a distributed mechanosensory organ in zebrafish in which most studies of hair cell regeneration have been conducted. In the maculae, we identify two subtypes of hair cells that share gene expression with mammalian striolar or extrastriolar hair cells. In situ hybridization reveals that these hair cell subtypes occupy distinct spatial domains within the three macular organs, the utricle, saccule, and lagena, consistent with the reported distinct electrophysiological properties of hair cells within these domains. These findings suggest that primitive specialization of spatially distinct striolar and extrastriolar hair cells likely arose in the last common ancestor of fish and mammals. The similarities of inner ear cell type composition between fish and mammals validate zebrafish as a relevant model for understanding inner ear-specific hair cell function and regeneration.

Thyroid hormone regulates abrupt skin morphogenesis during zebrafish postembryonic development.

Thyroid hormone is a key regulator of post-embryonic vertebrate development. Skin is a biomedically important thyroid hormone target organ, but the cellular and molecular mechanisms underlying skin pathologies associated with thyroid dysfunction remain obscure. The transparent skin of zebrafish is an accessible model system for studying vertebrate skin development. During post-embryonic development of the zebrafish, scales emerge in the skin from a hexagonally patterned array of dermal papillae, like other vertebrate skin appendages such as feathers and hair follicles. We show here that thyroid hormone regulates the rate of post-embryonic dermal development through interaction with nuclear hormone receptors. This couples skin development with body growth to generate a well ordered array of correctly proportioned scales. This work extends our knowledge of thyroid hormone actions on skin by providing in-vivo evidence that thyroid hormone regulates multiple aspects of dermal development.

A-kinase-anchoring protein 1 (dAKAP1)-based signaling complexes coordinate local protein synthesis at the mitochondrial surface.

Compartmentalization of macromolecules is a ubiquitous molecular mechanism that drives numerous cellular functions. The appropriate organization of enzymes in space and time enables the precise transmission and integration of intracellular signals. Molecular scaffolds constrain signaling enzymes to influence the regional modulation of these physiological processes. Mitochondrial targeting of protein kinases and protein phosphatases provides a means to locally control the phosphorylation status and action of proteins on the surface of this organelle. Dual-specificity protein kinase A anchoring protein 1 (dAKAP1) is a multivalent binding protein that targets protein kinase A (PKA), RNAs, and other signaling enzymes to the outer mitochondrial membrane. Many AKAPs recruit a diverse set of binding partners that coordinate a broad range of cellular processes. Here, results of MS and biochemical analyses reveal that dAKAP1 anchors additional components, including the ribonucleoprotein granule components La-related protein 4 (LARP4) and polyadenylate-binding protein 1 (PABPC1). Local translation of mRNAs at organelles is a means to spatially control the synthesis of proteins. RNA-Seq data demonstrate that dAKAP1 binds mRNAs encoding proteins required for mitochondrial metabolism, including succinate dehydrogenase. Functional studies suggest that the loss of dAKAP1-RNA interactions reduces mitochondrial electron transport chain activity. Hence, dAKAP1 plays a previously unappreciated role as a molecular interface between second messenger signaling and local protein synthesis machinery.

Thyroid hormone receptors mediate two distinct mechanisms of long-wavelength vision.

Thyroid hormone (TH) signaling plays an important role in the regulation of long-wavelength vision in vertebrates. In the retina, thyroid hormone receptor ß (thrb) is required for expression of long-wavelength-sensitive opsin (lws) in red cone photoreceptors, while in retinal pigment epithelium (RPE), TH regulates expression of a cytochrome P450 enzyme, cyp27c1, that converts vitamin A1 into vitamin A2 to produce a red-shifted chromophore. To better understand how TH controls these processes, we analyzed the phenotype of zebrafish with mutations in the three known TH nuclear receptor transcription factors (thraa, thrab, and thrb). We found that no single TH nuclear receptor is required for TH-mediated induction of cyp27c1 but that deletion of all three (thraa-/-;thrab-/-;thrb-/- ) completely abrogates its induction and the resulting conversion of A1- to A2-based retinoids. In the retina, loss of thrb resulted in an absence of red cones at both larval and adult stages without disruption of the underlying cone mosaic. RNA-sequencing analysis revealed significant down-regulation of only five genes in adult thrb-/- retina, of which three (lws1, lws2, and miR-726) occur in a single syntenic cluster. In the thrb-/- retina, retinal progenitors destined to become red cones were transfated into ultraviolet (UV) cones and horizontal cells. Taken together, our findings demonstrate cooperative regulation of cyp27c1 by TH receptors and a requirement for thrb in red cone fate determination. Thus, TH signaling coordinately regulates both spectral sensitivity and sensory plasticity.

Chromatin accessibility dynamics and single cell RNA-Seq reveal new regulators of regeneration in neural progenitors.

Vertebrate appendage regeneration requires precisely coordinated remodeling of the transcriptional landscape to enable the growth and differentiation of new tissue, a process executed over multiple days and across dozens of cell types. The heterogeneity of tissues and temporally-sensitive fate decisions involved has made it difficult to articulate the gene regulatory programs enabling regeneration of individual cell types. To better understand how a regenerative program is fulfilled by neural progenitor cells (NPCs) of the spinal cord, we analyzed pax6-expressing NPCs isolated from regenerating Xenopus tropicalis tails. By intersecting chromatin accessibility data with single-cell transcriptomics, we find that NPCs place an early priority on neuronal differentiation. Late in regeneration, the priority returns to proliferation. Our analyses identify Pbx3 and Meis1 as critical regulators of tail regeneration and axon organization. Overall, we use transcriptional regulatory dynamics to present a new model for cell fate decisions and their regulators in NPCs during regeneration.

Dimensionality reduction by UMAP to visualize physical and genetic interactions.

Dimensionality reduction is often used to visualize complex expression profiling data. Here, we use the Uniform Manifold Approximation and Projection (UMAP) method on published transcript profiles of 1484 single gene deletions of Saccharomyces cerevisiae. Proximity in low-dimensional UMAP space identifies groups of genes that correspond to protein complexes and pathways, and finds novel protein interactions, even within well-characterized complexes. This approach is more sensitive than previous methods and should be broadly useful as additional transcriptome datasets become available for other organisms.

A single-cell transcriptome atlas for zebrafish development.

The ability to define cell types and how they change during organogenesis is central to our understanding of animal development and human disease. Despite the crucial nature of this knowledge, we have yet to fully characterize all distinct cell types and the gene expression differences that generate cell types during development. To address this knowledge gap, we produced an atlas using single-cell RNA-sequencing methods to investigate gene expression from the pharyngula to early larval stages in developing zebrafish. Our single-cell transcriptome atlas encompasses transcriptional profiles from 44,102 â€‹cells across four days of development using duplicate experiments that confirmed high reproducibility. We annotated 220 identified clusters and highlighted several strategies for interrogating changes in gene expression associated with the development of zebrafish embryos at single-cell resolution. Furthermore, we highlight the power of this analysis to assign new cell-type or developmental stage-specific expression information to many genes, including those that are currently known only by sequence and/or that lack expression information altogether. The resulting atlas is a resource for biologists to generate hypotheses for functional analysis, which we hope integrates with existing efforts to define the diversity of cell-types during zebrafish organogenesis, and to examine the transcriptional profiles that produce each cell type over developmental time.

Fate plasticity and reprogramming in genetically distinct populations of Danio leucophores.

Understanding genetic and cellular bases of adult form remains a fundamental goal at the intersection of developmental and evolutionary biology. The skin pigment cells of vertebrates, derived from embryonic neural crest, are a useful system for elucidating mechanisms of fate specification, pattern formation, and how particular phenotypes impact organismal behavior and ecology. In a survey of Danio fishes, including the zebrafish Danio rerio, we identified two populations of white pigment cells-leucophores-one of which arises by transdifferentiation of adult melanophores and another of which develops from a yellow-orange xanthophore or xanthophore-like progenitor. Single-cell transcriptomic, mutational, chemical, and ultrastructural analyses of zebrafish leucophores revealed cell-type-specific chemical compositions, organelle configurations, and genetic requirements. At the organismal level, we identified distinct physiological responses of leucophores during environmental background matching, and we showed that leucophore complement influences behavior. Together, our studies reveal independently arisen pigment cell types and mechanisms of fate acquisition in zebrafish and illustrate how concerted analyses across hierarchical levels can provide insights into phenotypes and their evolution.

Thyroid hormone regulates distinct paths to maturation in pigment cell lineages.

Thyroid hormone (TH) regulates diverse developmental events and can drive disparate cellular outcomes. In zebrafish, TH has opposite effects on neural crest derived pigment cells of the adult stripe pattern, limiting melanophore population expansion, yet increasing yellow/orange xanthophore numbers. To learn how TH elicits seemingly opposite responses in cells having a common embryological origin, we analyzed individual transcriptomes from thousands of neural crest-derived cells, reconstructed developmental trajectories, identified pigment cell-lineage specific responses to TH, and assessed roles for TH receptors. We show that TH promotes maturation of both cell types but in distinct ways. In melanophores, TH drives terminal differentiation, limiting final cell numbers. In xanthophores, TH promotes accumulation of orange carotenoids, making the cells visible. TH receptors act primarily to repress these programs when TH is limiting. Our findings show how a single endocrine factor integrates very different cellular activities during the generation of adult form.

Notch regulates BMP responsiveness and lateral branching in vessel networks via SMAD6.

Functional blood vessel growth depends on generation of distinct but coordinated responses from endothelial cells. Bone morphogenetic proteins (BMP), part of the TGFß superfamily, bind receptors to induce phosphorylation and nuclear translocation of SMAD transcription factors (R-SMAD1/5/8) and regulate vessel growth. However, SMAD1/5/8 signalling results in both pro- and anti-angiogenic outputs, highlighting a poor understanding of the complexities of BMP signalling in the vasculature. Here we show that BMP6 and BMP2 ligands are pro-angiogenic in vitro and in vivo, and that lateral vessel branching requires threshold levels of R-SMAD phosphorylation. Endothelial cell responsiveness to these pro-angiogenic BMP ligands is regulated by Notch status and Notch sets responsiveness by regulating a cell-intrinsic BMP inhibitor, SMAD6, which affects BMP responses upstream of target gene expression. Thus, we reveal a paradigm for Notch-dependent regulation of angiogenesis: Notch regulates SMAD6 expression to affect BMP responsiveness of endothelial cells and new vessel branch formation.