Single-cell mapping of gene expression landscapes and lineage in the zebrafish embryo

Single-Cell RNA-Seq by Multiple Annealing and Tailing-Based Quantitative Single-Cell RNA-Seq (MATQ-Seq)

Single-cell technologies have emerged as advanced tools to study various biological processes that demand the single cell resolution. To detect subtle heterogeneity in the transcriptome, high accuracy and sensitivity are still desired for single-cell RNA-seq. We describe here multiple annealing and dC-tailing-based quantitative single-cell RNA-seq (MATQ-seq) with ~90% capture efficiency. In addition, MATQ-seq is a total RNA assay allowing for detection of nonpolyadenylated transcripts.

The flow responsive transcription factor Klf2 is required for myocardial wall integrity by modulating Fgf signaling

Complex interplay between cardiac tissues is crucial for their integrity. The flow responsive transcription factor KLF2, which is expressed in the endocardium, is vital for cardiovascular development but its exact role remains to be defined. To this end, we mutated both klf2 paralogues in zebrafish, and while single mutants exhibit no obvious phenotype, double mutants display a novel phenotype of cardiomyocyte extrusion towards the abluminal side. This extrusion requires cardiac contractility and correlates with the mislocalization of N-cadherin from the lateral to the apical side of cardiomyocytes. Transgenic rescue data show that klf2 expression in endothelium, but not myocardium, prevents this cardiomyocyte extrusion phenotype. Transcriptome analysis of klf2 mutant hearts reveals that Fgf signaling is affected, and accordingly, we find that inhibition of Fgf signaling in wild-type animals can lead to abluminal cardiomyocyte extrusion. These studies provide new insights into how Klf2 regulates cardiovascular development and specifically myocardial wall integrity.

Combined Numerical and Experimental Investigation of Localized Electroporation-Based Cell Transfection and Sampling

Localized electroporation has evolved as an effective technology for the delivery of foreign molecules into cells while preserving their viability. Consequently, this technique has potential applications in sampling the contents of live cells and the temporal assessment of cellular states at the single-cell level. Although there have been numerous experimental reports on localized electroporation-based delivery, a lack of a mechanistic understanding of the process hinders its implementation in sampling. In this work, we develop a multiphysics model that predicts the transport of molecules into and out of the cell during localized electroporation. Based on the model predictions, we optimize experimental parameters such as buffer conditions, electric field strength, cell confluency, and density of nanochannels in the substrate for successful delivery and sampling via localized electroporation. We also identify that cell membrane tension plays a crucial role in enhancing both the amount and the uniformity of molecular transport, particularly for macromolecules. We qualitatively validate the model predictions on a localized electroporation platform by delivering large molecules (bovine serum albumin and mCherry-encoding plasmid) and by sampling an exogeneous protein (tdTomato) in an engineered cell line.

An update on clonality: what smooth muscle cell type makes up the atherosclerotic plaque?

Almost 50 years ago, Earl Benditt and his son John described the clonality of the atherosclerotic plaque. This led Benditt to propose that the atherosclerotic lesion was a smooth muscle neoplasm, similar to the leiomyomata seen in the uterus of most women. Although the observation of clonality has been confirmed many times, interest in the idea that atherosclerosis might be a form of neoplasia waned because of the clinical success of treatments for hyperlipemia and because animal models have made great progress in understanding how lipid accumulates in the plaque and may lead to plaque rupture. Four advances have made it important to reconsider Benditt's observations. First, we now know that clonality is a property of normal tissue development. Second, this is even true in the vessel wall, where we now know that formation of clonal patches in that wall is part of the development of smooth muscle cells that make up the tunica media of arteries. Third, we know that the intima, the "soil" for development of the human atherosclerotic lesion, develops before the fatty lesions appear. Fourth, while the cells comprising this intima have been called "smooth muscle cells", we do not have a clear definition of cell type nor do we know if the initial accumulation is clonal. As a result, Benditt's hypothesis needs to be revisited in terms of changes in how we define smooth muscle cells and the quite distinct developmental origins of the cells that comprise the muscular coats of all arterial walls. Finally, since clonality of the lesions is real, the obvious questions are do these human tumors precede the development of atherosclerosis, how do the clones develop, what cell type gives rise to the clones, and in what ways do the clones provide the soil for development and natural history of atherosclerosis?

Single-cell mapping of lineage and identity in direct reprogramming

Direct lineage reprogramming involves the conversion of cellular identity. Single-cell technologies are useful for deconstructing the considerable heterogeneity that emerges during lineage conversion. However, lineage relationships are typically lost during cell processing, complicating trajectory reconstruction. Here we present 'CellTagging', a combinatorial cell-indexing methodology that enables parallel capture of clonal history and cell identity, in which sequential rounds of cell labelling enable the construction of multi-level lineage trees. CellTagging and longitudinal tracking of fibroblast to induced endoderm progenitor reprogramming reveals two distinct trajectories: one leading to successfully reprogrammed cells, and one leading to a 'dead-end' state, paths determined in the earliest stages of lineage conversion. We find that expression of a putative methyltransferase, Mettl7a1, is associated with the successful reprogramming trajectory; adding Mettl7a1 to the reprogramming cocktail increases the yield of induced endoderm progenitors. Together, these results demonstrate the utility of our lineage-tracing method for revealing the dynamics of direct reprogramming.

In Toto Imaging of Dynamic Osteoblast Behaviors in Regenerating Skeletal Bone

Osteoblasts are matrix-depositing cells that can divide and heal bone injuries. Their deep-tissue location and the slow progression of bone regeneration challenge attempts to capture osteoblast behaviors in live tissue at high spatiotemporal resolution. Here, we have developed an imaging platform to monitor and quantify individual and collective behaviors of osteoblasts in adult zebrafish scales, skeletal body armor discs that regenerate rapidly after loss. Using a panel of transgenic lines that visualize and manipulate osteoblasts, we find that a founder pool of osteoblasts emerges through de novo differentiation within one day of scale plucking. These osteoblasts undergo division events that are largely uniform in frequency and orientation to establish a primordium. Osteoblast proliferation dynamics diversify across the primordium by two days after injury, with cell divisions focused near, and with orientations parallel to, the scale periphery, occurring coincident with dynamic localization of fgf20a gene expression. In posterior scale regions, cell elongation events initiate in areas soon occupied by mineralized grooves called radii, beginning approximately 2 days post injury, with patterned osteoblast death events accompanying maturation of these radii. By imaging at single-cell resolution, we detail acquisition of spatiotemporally distinct cell division, motility, and death dynamics within a founder osteoblast pool as bone regenerates.

Revealing the Critical Regulators of Cell Identity in the Mouse Cell Atlas

Recent progress in single-cell technologies has enabled the identification of all major cell types in mouse. However, for most cell types, the regulatory mechanism underlying their identity remains poorly understood. By computational analysis of the recently published mouse cell atlas data, we have identified 202 regulons whose activities are highly variable across different cell types, and more importantly, predicted a small set of essential regulators for each major cell type in mouse. Systematic validation by automated literature and data mining provides strong additional support for our predictions. Thus, these predictions serve as a valuable resource that would be useful for the broad biological community. Finally, we have built a user-friendly, interactive web portal to enable users to navigate this mouse cell network atlas.

Large-scale reconstruction of cell lineages using single-cell readout of transcriptomes and CRISPR-Cas9 barcodes by scGESTALT

Lineage relationships among the large number of heterogeneous cell types generated during development are difficult to reconstruct in a high-throughput manner. We recently established a method, scGESTALT, that combines cumulative editing of a lineage barcode array by CRISPR-Cas9 with large-scale transcriptional profiling using droplet-based single-cell RNA sequencing (scRNA-seq). The technique generates edits in the barcode array over multiple timepoints using Cas9 and pools of single-guide RNAs (sgRNAs) introduced during early and late zebrafish embryonic development, which distinguishes it from similar Cas9 lineage-tracing methods. The recorded lineages are captured, along with thousands of cellular transcriptomes, to build lineage trees with hundreds of branches representing relationships among profiled cell types. Here, we provide details for (i) generating transgenic zebrafish; (ii) performing multi-timepoint barcode editing; (iii) building scRNA-seq libraries from brain tissue; and (iv) concurrently amplifying lineage barcodes from captured single cells. Generating transgenic lines takes 6 months, and performing barcode editing and generating single-cell libraries involve 7 d of hands-on time. scGESTALT provides a scalable platform to map lineage relationships between cell types in any system that permits genome editing during development, regeneration, or disease.

Cellular barcoding: lineage tracing, screening and beyond

Cellular barcoding is a technique in which individual cells are labeled with unique nucleic acid sequences, termed barcodes, so that they can be tracked through space and time. Cellular barcoding can be used to track millions of cells in parallel, and thus is an efficient approach for investigating heterogeneous populations of cells. Over the past 25 years, cellular barcoding has been used for fate mapping, lineage tracing and high-throughput screening, and has led to important insights into developmental biology and gene function. Driven by plummeting sequencing costs and the power of synthetic biology, barcoding is now expanding beyond traditional applications and into diverse fields such as neuroanatomy and the recording of cellular activity. In this review, we discuss the fundamental principles of cellular barcoding, including the underlying mathematics, and its applications in both new and established fields.

Control of Lineage-Specific Gene Expression by Functionalized gRNA Barcodes

Lineage tracking delivers essential quantitative insight into dynamic, probabilistic cellular processes, such as somatic tumor evolution and differentiation. Methods for high diversity lineage quantitation rely on sequencing a population of DNA barcodes. However, manipulation of specific individual lineages is not possible with this approach. To address this challenge, we developed a functionalized lineage tracing tool, Control of Lineages by Barcode Enabled Recombinant Transcription (COLBERT), that enables high diversity lineage tracing and lineage-specific manipulation of gene expression. This modular platform utilizes expressed barcode gRNAs to both track cell lineages and direct lineage-specific gene expression.

Ionizing radiation induces transgenerational effects of DNA methylation in zebrafish

Ionizing radiation is known to cause DNA damage, yet the mechanisms underlying potential transgenerational effects of exposure have been scarcely studied. Previously, we observed effects in offspring of zebrafish exposed to gamma radiation during gametogenesis. Here, we hypothesize that these effects are accompanied by changes of DNA methylation possibly inherited by subsequent generations. We assessed DNA methylation in F1 embryos (5.5 hours post fertilization) with whole genome bisulfite sequencing following parental exposure to 8.7 mGy/h for 27 days and found 5658 differentially methylated regions (DMRs). DMRs were predominantly located at known regulatory regions, such as gene promoters and enhancers. Pathway analysis indicated the involvement of DMRs related to similar pathways found with gene expression analysis, such as development, apoptosis and cancers, which could be linked to previous observed developmental defects and genomic instability in the offspring. Follow up of 19 F1 DMRs in F2 and F3 embryos revealed persistent effects up to the F3 generation at 5 regions. These results indicate that ionizing radiation related effects in offspring can be linked to DNA methylation changes that partly can persist over generations. Monitoring DNA methylation could serve as a biomarker to provide an indication of ancestral exposures to ionizing radiation.

Bloody Zebrafish: Novel Methods in Normal and Malignant Hematopoiesis

Hematopoiesis is an optimal system for studying stem cell maintenance and lineage differentiation under physiological and pathological conditions. In vertebrate organisms, billions of differentiated hematopoietic cells need to be continuously produced to replenish the blood cell pool. Disruptions in this process have immediate consequences for oxygen transport, responses against pathogens, maintenance of hemostasis and vascular integrity. Zebrafish is a widely used and well-established model for studying the hematopoietic system. Several new hematopoietic regulators were identified in genetic and chemical screens using the zebrafish model. Moreover, zebrafish enables in vivo imaging of hematopoietic stem cell generation and differentiation during embryogenesis, and adulthood. Finally, zebrafish has been used to model hematopoietic diseases. Recent technological advances in single-cell transcriptome analysis, epigenetic regulation, proteomics, metabolomics, and processing of large data sets promise to transform the current understanding of normal, abnormal, and malignant hematopoiesis. In this perspective, we discuss how the zebrafish model has proven beneficial for studying physiological and pathological hematopoiesis and how these novel technologies are transforming the field.

Developmental asynchrony and antagonism of sex determination pathways in a lizard with temperature-induced sex reversal

Vertebrate sex differentiation follows a conserved suite of developmental events: the bipotential gonads differentiate and shortly thereafter sex specific traits become dimorphic. However, this may not apply to squamates, a diverse vertebrate lineage comprising of many species with thermosensitive sexual development. Of the three species with data on the relative timing of gonad differentiation and genital dimorphism, the females of two (Niveoscincus ocellatus and Barisia imbricata) exhibit a phase of temporary pseudohermaphroditism or TPH (gonads have differentiated well before genital dimorphism). We report a third example of TPH in Pogona vitticeps, an agamid with temperature-induced male to female sex reversal. These findings suggest that for female squamates, genital and gonad development may not be closely synchronised, so that TPH may be common. We further observed a high frequency of ovotestes, a usually rare gonadal phenotype characterised by a mix of male and female structures, exclusively associated with temperature-induced sex reversal. We propose that ovotestes are evidence of a period of antagonism between male and female sex-determining pathways during sex reversal. Female sexual development in squamates is considerably more complex than has been appreciated, providing numerous avenues for future exploration of the genetic and hormonal cues that govern sexual development.

Regulatory cocktail for dopaminergic neurons in a protovertebrate identified by whole-embryo single-cell transcriptomics

The CNS of the protovertebrate Ciona intestinalis contains a single cluster of dopaminergic (DA) neurons, the coronet cells, which have been likened to the hypothalamus of vertebrates. Whole-embryo single-cell RNA sequencing (RNA-seq) assays identified Ptf1a as the most strongly expressed cell-specific transcription factor (TF) in DA/coronet cells. Knockdown of Ptf1a activity results in their loss, while misexpression results in the appearance of supernumerary DA/coronet cells. Photoreceptor cells and ependymal cells are the most susceptible to transformation, and both cell types express high levels of Meis Coexpression of both Ptf1a and Meis caused the wholesale transformation of the entire CNS into DA/coronet cells. We therefore suggest that the reiterative use of functional manipulations and single-cell RNA-seq assays is an effective means for the identification of regulatory cocktails underlying the specification of specific cell identities.

Elucidating the developmental trajectories of GABAergic cortical interneuron subtypes

GABAergic interneurons in the neocortex play pivotal roles in the feedforward and feedback inhibition that control higher order information processing and thus, malfunction in the inhibitory circuits often leads to neurodevelopmental disorders. Very interestingly, a large diversity of morphology, synaptic targeting specificity, electrophysiological properties and molecular expression profiles are found in cortical interneurons, which originate within the distantly located embryonic ganglionic eminences. Here, I will review the still ongoing effort to understand the developmental trajectories of GABAergic cortical interneuron subtypes.

Single-Cell Multi-omics: An Engine for New Quantitative Models of Gene Regulation

Cells in a multicellular organism fulfill specific functions by enacting cell-type-specific programs of gene regulation. Single-cell RNA sequencing technologies have provided a transformative view of cell-type-specific gene expression, the output of cell-type-specific gene regulatory programs. This review discusses new single-cell genomic technologies that complement single-cell RNA sequencing by providing additional readouts of cellular state beyond the transcriptome. We highlight regression models as a simple yet powerful approach to relate gene expression to other aspects of cellular state, and in doing so, gain insights into the biochemical mechanisms that are necessary to produce a given gene expression output.

Single-cell genomics to guide human stem cell and tissue engineering

To understand human development and disease, as well as to regenerate damaged tissues, scientists are working to engineer certain cell types in vitro and to create 3D microenvironments in which cells behave physiologically. Single-cell genomics (SCG) technologies are being applied to primary human organs and to engineered cells and tissues to generate atlases of cell diversity in these systems at unparalleled resolution. Moving beyond atlases, SCG methods are powerful tools for gaining insight into the engineering and disease process. Here we discuss how scientists can use single-cell sequencing to optimize human cell and tissue engineering by measuring precision, detecting inefficiencies, and assessing accuracy. We also provide a perspective on how emerging SCG methods can be used to reverse-engineer human cells and tissues and unravel disease mechanisms.

NeuroSystematics and Periodic System of Neurons: Model vs Reference Species at Single-Cell Resolution

There is more than one way to develop neuronal complexity, and animals frequently use different molecular toolkits to achieve similar functional outcomes (=convergent evolution). Neurons are different not only because they have different functions, but also because neurons and circuits have different genealogies, and perhaps independent origins at the broadest scale from ctenophores and cnidarians to cephalopods and primates. By combining modern phylogenomics, single-neuron sequencing (scRNA-seq), machine learning, single-cell proteomics, and metabolomic across Metazoa, it is possible to reconstruct the evolutionary histories of neurons tracing them to ancestral secretory cells. Comparative data suggest that neurons, and perhaps synapses, evolved at least 2-3 times (in ctenophore, cnidarian and bilateral lineages) during ∼600 million years of animal evolution. There were also several independent events of the nervous system centralization either from a common bilateral/cnidarian ancestor without the bona fide neurons or from the urbilaterian with diffuse, nerve-net type nervous system. From the evolutionary standpoint, (i) a neuron should be viewed as a functional rather than a genetic character, and (ii) any given neural system might be chimeric and composed of different cell lineages with distinct origins and evolutionary histories. The identification of distant neural homologies or examples of convergent evolution among 34 phyla will not only allow the reconstruction of neural systems' evolution but together with single-cell "omic" approaches the proposed synthesis would lead to the "Periodic System of Neurons" with predictive power for neuronal phenotypes and plasticity. Such a phylogenetic classification framework of Neuronal Systematics (NeuroSystematics) might be a conceptual analog of the Periodic System of Chemical Elements. scRNA-seq profiling of all neurons in an entire brain or Brain-seq is now fully achievable in many nontraditional reference species across the entire animal kingdom. Arguably, marine animals are the most suitable for the proposed tasks because the world oceans represent the greatest taxonomic and body-plan diversity.

Problems in Fish-to-Tetrapod Transition: Genetic Expeditions Into Old Specimens

The fish-to-tetrapod transition is one of the fundamental problems in evolutionary biology. A significant amount of paleontological data has revealed the morphological trajectories of skeletons, such as those of the skull, vertebrae, and appendages in vertebrate history. Shifts in bone differentiation, from dermal to endochondral bones, are key to explaining skeletal transformations during the transition from water to land. However, the genetic underpinnings underlying the evolution of dermal and endochondral bones are largely missing. Recent genetic approaches utilizing model organisms—zebrafish, frogs, chickens, and mice—reveal the molecular mechanisms underlying vertebrate skeletal development and provide new insights for how the skeletal system has evolved. Currently, our experimental horizons to test evolutionary hypotheses are being expanded to non-model organisms with state-of-the-art techniques in molecular biology and imaging. An integration of functional genomics, developmental genetics, and high-resolution CT scanning into evolutionary inquiries allows us to reevaluate our understanding of old specimens. Here, we summarize the current perspectives in genetic programs underlying the development and evolution of the dermal skull roof, shoulder girdle, and appendages. The ratio shifts of dermal and endochondral bones, and its underlying mechanisms, during the fish-to-tetrapod transition are particularly emphasized. Recent studies have suggested the novel cell origins of dermal bones, and the interchangeability between dermal and endochondral bones, obscuring the ontogenetic distinction of these two types of bones. Assimilation of ontogenetic knowledge of dermal and endochondral bones from different structures demands revisions of the prevalent consensus in the evolutionary mechanisms of vertebrate skeletal shifts.

Single-cell biology: resolving biological complexity, one cell at a time

In March 2018, over 250 researchers came together at the Wellcome Genome Campus in Hinxton, Cambridge, UK, to present their latest research in the area of single-cell biology. A highly interdisciplinary meeting, the Single Cell Biology conference covered a variety of topics, ranging from cutting-edge technological innovation, developmental biology and stem cell research to evolution and cancer. This meeting report summarises the key findings presented and the major research themes that emerged during the conference.

Hematopoietic Hierarchy - An Updated Roadmap

The classical roadmap of hematopoietic hierarchy has been proposed for nearly 20 years and has become a dogma of stem cell research for most types of adult stem cells, including hematopoietic stem cells (HSCs). However, with the development of new technologies such as omics approaches at single-cell resolution, recent studies in vitro and in vivo have suggested that heterogeneity is a common feature of HSCs and their progenies. While these findings broaden our understanding of hematopoiesis, they also challenge the well-accepted hematopoietic hierarchy roadmap. Here, we review recent advances in the hematopoiesis field and provide an updated view to incorporate these new findings as well as to reflect on the complexity of HSCs and their derivatives in development and adulthood.

The dynamics of gene expression in vertebrate embryogenesis at single-cell resolution

Time series of single-cell transcriptome measurements can reveal dynamic features of cell differentiation pathways. From measurements of whole frog embryos spanning zygotic genome activation through early organogenesis, we derived a detailed catalog of cell states in vertebrate development and a map of differentiation across all lineages over time. The inferred map recapitulates most if not all developmental relationships and associates new regulators and marker genes with each cell state. We find that many embryonic cell states appear earlier than previously appreciated. We also assess conflicting models of neural crest development. Incorporating a matched time series of zebrafish development from a companion paper, we reveal conserved and divergent features of vertebrate early developmental gene expression programs.