Search

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 mechanisms that regulate stripe patterning and scale morphogenesis may lead to discovery of fundamental mechanisms that govern 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 discovery of mechanisms that govern skin development and homeostasis.

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 cell and non-sensory supporting cells. We identify a putative progenitor population for hair cells and supporting cells, as well as distinct hair cells 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 two major macular organs, the utricle and saccule, 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 also support using zebrafish as a relevant model for understanding inner ear-specific hair cell function and regeneration.

The maturation of single cell transcriptomic technologies has facilitated the generation of comprehensive cellular atlases from whole embryos. A majority of this data, however, has been collected from wild type embryos without an appreciation for latent variation present in development. Here we present single cell transcriptomic data from 1812 individually resolved developing zebrafish embryos, encompassing 19 time points, 23 genetic perturbations, and totaling 3.2M cells. The high degree of replication in our study (8 or more embryos per condition) allows 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 the 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 the genetic dependencies of zebrafish cell types, while also addressing long-standing challenges in developmental genetics, including the cellular and transcriptional plasticity underlying phenotypic diversity across individuals.

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.

The genetic program of embryonic development is remarkably robust, but temperature stress can degrade its ability to generate animals with invariant anatomy. While the stereotyped, consistent phenotypes associated with environmental stress during vertebrate development suggest that some cell types are more sensitive to stress than others, the basis of this sensitivity is unknown. Here, we characterize hundreds of individual zebrafish embryos under temperature stress using whole-animal single cell RNA-seq to identify cell types and molecular programs within them that drive phenotypic variability. We find that temperature perturbs the normal proportions and gene expression programs of numerous cell types and also introduces asynchrony in their development. The notochord is particularly sensitive to temperature stress, which we show is due to a specialized cell type, sheath cells. Further analyses show that sheath cells accumulate misfolded protein at elevated temperature, leading to a cascading structural failure of the notochord and irreversible anatomic defects in the embryo. Our study demonstrates that whole-animal single cell RNA-seq can characterize mechanisms important for developmental robustness and pinpoint molecular programs within specific cell types that constitute key failure points.

Single-cell chemical transcriptomics Single-cell transcriptomic technologies have emerged as powerful tools to explore cellular heterogeneity at the resolution of individual cells. Srivatsan et al. now add another layer of information and complexity by combining single-cell transcriptomics with oligo hashing and small molecule screening in a method called sci-Plex. Because screening many chemical compounds requires the ability to profile many cells, and because screens perturb cells in many different ways, the authors demonstrate the effects of 188 compounds in three cancer lines. The sci-Plex method can capture gene expression profiles from thousands of experimental conditions in a single experiment. Science , this issue p. 45

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.HighlightsThyroid hormone (TH) is necessary for normal squamation patterning in zebrafish.Stratified dermis develops by migration of primary hypodermal cells.Dermis stratifies in an invariant wave.TH regulates the rates of multiple aspects of dermis development.Scale size and density are sensitive to skin size at onset of squamation.

Single cell transcriptomic profiling can map developmental trajectories and assess cell-type–specific changes during heat shock at single cell resolution., Single cell RNA sequencing can yield high-resolution cell-type–specific expression signatures that reveal new cell types and the developmental trajectories of cell lineages. Here, we apply this approach to Arabidopsis (Arabidopsis thaliana) root cells to capture gene expression in 3,121 root cells. We analyze these data with Monocle 3, which orders single cell transcriptomes in an unsupervised manner and uses machine learning to reconstruct single cell developmental trajectories along pseudotime. We identify hundreds of genes with cell-type–specific expression, with pseudotime analysis of several cell lineages revealing both known and novel genes that are expressed along a developmental trajectory. We identify transcription factor motifs that are enriched in early and late cells, together with the corresponding candidate transcription factors that likely drive the observed expression patterns. We assess and interpret changes in total RNA expression along developmental trajectories and show that trajectory branch points mark developmental decisions. Finally, by applying heat stress to whole seedlings, we address the longstanding question of possible heterogeneity among cell types in the response to an abiotic stress. Although the response of canonical heat-shock genes dominates expression across cell types, subtle but significant differences in other genes can be detected among cell types. Taken together, our results demonstrate that single cell transcriptomics holds promise for studying plant development and plant physiology with unprecedented resolution.

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 A 1 into vitamin A 2 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 A 1 - to A 2 -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.

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 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., Dimensionality reduction is often used to visualize expression profiling data in order to find relationships among cells. Here, the authors use Uniform Manifold Approximation and Projection (UMAP) on published expression data of gene deletions of S. cerevisiae to find novel protein interactions.

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 of biologists to generate hypotheses for genetic (mutant) or functional analysis, to launch an effort to define the diversity of cell-types during zebrafish organogenesis, and to examine the transcriptional profiles that produce each cell type over developmental time.

Significance A foundational question in biology is how phenotypically similar traits arise. Here we identify two distinct white pigment cell populations, leucophores, that arise from independent progenitors in zebrafish and its relatives. Remarkably, one of these leucophore populations develops from previously differentiated, black melanophores. These different leucophores exhibited distinct pigment chemistries, cytological features, gene expression profiles, and genetic requirements, and whole-animal experiments implicated these cells in behavioral interactions. Our several approaches provide insights into pigment cell complements and origins in zebrafish and contribute to our understanding of form and function in the spectacular pigment patterns of teleost fishes., 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.

Single-cell RNA-seq can yield high-resolution cell-type-specific expression signatures that reveal new cell types and the developmental trajectories of cell lineages. Here, we apply this approach toA. thalianaroot cells to capture gene expression in 3,121 root cells. We analyze these data with Monocle 3, which orders single cell transcriptomes in an unsupervised manner and uses machine learning to reconstruct single-cell developmental trajectories along pseudotime. We identify hundreds of genes with cell-type-specific expression, with pseudotime analysis of several cell lineages revealing both known and novel genes that are expressed along a developmental trajectory. We identify transcription factor motifs that are enriched in early and late cells, together with the corresponding candidate transcription factors that likely drive the observed expression patterns. We assess and interpret changes in total RNA expression along developmental trajectories and show that trajectory branch points mark developmental decisions. Finally, by applying heat stress to whole seedlings, we address the longstanding question of possible heterogeneity among cell types in the response to an abiotic stress. Although the response of canonical heat shock genes dominates expression across cell types, subtle but significant differences in other genes can be detected among cell types. Taken together, our results demonstrate that single-cell transcriptomics holds promise for studying plant development and plant physiology with unprecedented resolution.

Single-cell transcriptome sequencing now routinely samples thousands of cells, potentially providing enough data to reconstruct causal gene regulatory networks from observational data. Here, we present Scribe, a toolkit for detecting and visualizing causal regulatory interactions between genes and explore the potential for single-cell experiments to power network reconstruction. Scribe employs Restricted Directed Information to determine causality by estimating the strength of information transferred from a potential regulator to its downstream target. We apply Scribe and other leading approaches for causal network reconstruction to several types of single-cell measurements and show that there is a dramatic drop in performance for "pseudotime” ordered single-cell data compared to true time series data. We demonstrate that performing causal inference requires temporal coupling between measurements. We show that methods such as “RNA velocity” restore some degree of coupling through an analysis of chromaffin cell fate commitment. These analyses therefore highlight an important shortcoming in experimental and computational methods for analyzing gene regulation at single-cell resolution and point the way towards overcoming it.

Quantitative PCR is a diagnostic pillar for clinical virology testing, and reference materials are necessary for accurate, comparable quantitation between clinical laboratories. Accurate quantitation of human herpesvirus 6A/B (HHV-6A/B) is important for detection of viral reactivation and inherited chromosomally integrated HHV-6A/B in immunocompromised patients. Reference materials in clinical virology commonly consist of laboratory-adapted viral strains that may be affected by the culture process. We performed next-generation sequencing to make relative copy number measurements at single nucleotide resolution of eight candidate HHV-6A and seven HHV-6B reference strains and DNA materials from the HHV-6 Foundation and Advanced Biotechnologies Inc. Eleven of 17 (65%) HHV-6A/B candidate reference materials showed multiple copies of the origin of replication upstream of the U41 gene by next-generation sequencing. These large tandem repeats arose independently in culture-adapted HHV-6A and HHV-6B strains, measuring 1,254 bp and 983 bp, respectively. The average copy number measured was between 5 and 10 times the number of copies of the rest of the genome. We also report the first interspecies recombinant HHV-6A/B strain with a HHV-6A backbone and a >5.5-kb region from HHV-6B, from U41 to U43, that covered the origin tandem repeat. Specific HHV-6A reference strains demonstrated duplication of regions at U1/U2, U87, and U89, as well as deletion in the U12-to-U24 region and the U94/U95 genes. HHV-6A/B strains derived from cord blood mononuclear cells from different laboratories on different continents with fewer passages revealed no copy number differences throughout the viral genome. These data indicate that large origin tandem duplications are an adaptation of both HHV-6A and HHV-6B in culture and show interspecies recombination is possible within the Betaherpesvirinae. IMPORTANCE Anything in science that needs to be quantitated requires a standard unit of measurement. This includes viruses, for which quantitation increasingly determines definitions of pathology and guidelines for treatment. However, the act of making standard or reference material in virology can alter its very accuracy through genomic duplications, insertions, and rearrangements. We used deep sequencing to examine candidate reference strains for HHV-6, a ubiquitous human virus that can reactivate in the immunocompromised population and is integrated into the human genome in every cell of the body for 1% of people worldwide. We found large tandem repeats in the origin of replication for both HHV-6A and HHV-6B that are selected for in culture. We also found the first interspecies recombinant between HHV-6A and HHV-6B, a phenomenon that is well known in alphaherpesviruses but to date has not been seen in betaherpesviruses. These data critically inform HHV-6A/B biology and the standard selection process.

Quantitative PCR is the diagnostic pillar for clinical virology testing, and reference materials are necessary for accurate, comparable quantitation between clinical laboratories. Accurate quantitation of HHV-6 is important for detection of viral reactivation and inherited chromosomally integrated HHV-6 in immunocompromised patients. Reference materials in clinical virology commonly consist of laboratory-adapted viral strains that may be affected by the culture process. We performed next-generation sequencing to make relative copy number measurements at single nucleotide resolution of eight candidate HHV-6A and seven HHV-6B reference strains and DNA materials from the HHV-6 Foundation and Advanced Biotechnologies. 11 of 17 (65%) HHV6 candidate reference materials showed multiple copies of the origin of replication upstream of the U41 gene by next-generation sequencing. These large tandem repeats arose independently in culture-adapted HHV-6A and HHV-6B strains, measuring 1254 bp and 983 bp, respectively. Copy number measured between 4-10X copies relative to the rest of the genome. We also report the first interspecies recombinant HHV-6 strain with a HHV-6A GS backbone and >5.5kb region from HHV-6B Z29 from U41-U43 that covered the origin tandem repeat. Specific HHV-6A reference strains demonstrated duplication of regions at UL1/UL2, U87, and U89, as well as deletion in the U12-U24 region and U94/95 genes. HHV-6 strains derived from cord blood mononuclear cells from different labs on different continents revealed no copy number differences throughout the viral genome. These data indicate large origin tandem duplications are an adaptation of both HHV-6A and HHV-6B in culture and show interspecies recombination is possible within theBetaherpesvirinae.ImportanceAnything in science that needs to be quantitated requires a standard unit of measurement. This includes viruses, for which quantitation increasingly determines definitions of pathology and guidelines for treatment. However, the act of making standard or reference material in virology can alter its very usefulness through genomic duplications, insertions, and rearrangements. We used deep sequencing to examine candidate reference strains for HHV-6, a ubiquitous human virus that can reactivate in the immunocompromised population and is integrated into the human genome in every cell of the body for 1% of people worldwide. We found large tandem repeats in the origin of replication for both HHV-6A and HHV-6B that are selected for in culture. We also found the first interspecies recombinant between HHV-6A and HHV-6B, a phenomenon that is well-known in alphaherpesviruses but to date has not been seen in betaherpesviruses. These data critically inform HHV-6 biology and the standard selection process.

Here, we present Scribe (https://github.com/aristoteleo/Scribe-py), a toolkit for detecting and visualizing causal regulatory interactions between genes and explore the potential for single-cell experiments to power network reconstruction. Scribe employs restricted directed information to determine causality by estimating the strength of information transferred from a potential regulator to its downstream target. We apply Scribe and other leading approaches for causal network reconstruction to several types of single-cell measurements and show that there is a dramatic drop in performance for "pseudotime"-ordered single-cell data compared with true time-series data. We demonstrate that performing causal inference requires temporal coupling between measurements. We show that methods such as "RNA velocity" restore some degree of coupling through an analysis of chromaffin cell fate commitment. These analyses highlight a shortcoming in experimental and computational methods for analyzing gene regulation at single-cell resolution and suggest ways of overcoming it.

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., The mechanism underlying endothelial cell responses to BMP signals is unknown. Here, the authors show that the endothelial response to pro-angiogenic BMP ligands is regulated by Notch via its effect on SMAD6, a known inhibitor of BMP intracellular signaling cascade.