Spatial localization is a key determinant of cellular fate and behavior, but spatial RNA assays traditionally rely on staining for a limited number of RNA species. In contrast, single-cell RNA-seq allows for deep profiling of cellular gene expression, but established methods separate cells from their native spatial context. Here we present Seurat, a computational strategy to infer cellular localization by integrating single-cell RNA-seq data with in situ RNA patterns. We applied Seurat to spatially map 851 single cells from dissociated zebrafish (Danio rerio) embryos, inferring a transcriptome-wide map of spatial patterning. We confirmed Seurat’s accuracy using several experimental approaches, and used it to identify a set of archetypal expression patterns and spatial markers. Additionally, Seurat correctly localizes rare subpopulations, accurately mapping both spatially restricted and scattered groups. Seurat will be applicable to mapping cellular localization within complex patterned tissues in diverse systems.

MicroRNAs (miRNAs) comprise 1 to 3% of all vertebrate genes, but their in vivo functions and mechanisms of action remain largely unknown. Zebrafish miR-430 is expressed at the onset of zygotic transcription and regulates morphogenesis during early development. By using a microarray approach and in vivo target validation, we find that miR-430 directly regulates several hundred target messenger RNA molecules (mRNAs). Most targets are maternally expressed mRNAs that accumulate in the absence of miR-430. We also show that miR-430 accelerates the deadenylation of target mRNAs. These results suggest that miR-430 facilitates the deadenylation and clearance of maternal mRNAs during early embryogenesis.

MicroRNAs (miRNAs) are small RNAs that regulate gene expression posttranscriptionally. To block all miRNA formation in zebrafish, we generated maternal-zygotic dicer (MZdicer) mutants that disrupt the Dicer ribonuclease III and double-stranded RNA-binding domains. Mutant embryos do not process precursor miRNAs into mature miRNAs, but injection of preprocessed miRNAs restores gene silencing, indicating that the disrupted domains are dispensable for later steps in silencing. MZdicer mutants undergo axis formation and differentiate multiple cell types but display abnormal morphogenesis during gastrulation, brain formation, somitogenesis, and heart development. Injection of miR-430 miRNAs rescues the brain defects in MZdicer mutants, revealing essential roles for miRNAs during morphogenesis.

Long noncoding RNAs (lncRNAs) comprise a diverse class of transcripts that structurally resemble mRNAs but do not encode proteins. Recent genome-wide studies in humans and the mouse have annotated lncRNAs expressed in cell lines and adult tissues, but a systematic analysis of lncRNAs expressed during vertebrate embryogenesis has been elusive. To identify lncRNAs with potential functions in vertebrate embryogenesis, we performed a time-series of RNA-seq experiments at eight stages during early zebrafish development. We reconstructed 56,535 high-confidence transcripts in 28,912 loci, recovering the vast majority of expressed RefSeq transcripts while identifying thousands of novel isoforms and expressed loci. We defined a stringent set of 1133 noncoding multi-exonic transcripts expressed during embryogenesis. These include long intergenic ncRNAs (lincRNAs), intronic overlapping lncRNAs, exonic antisense overlapping lncRNAs, and precursors for small RNAs (sRNAs). Zebrafish lncRNAs share many of the characteristics of their mammalian counterparts: relatively short length, low exon number, low expression, and conservation levels comparable to that of introns. Subsets of lncRNAs carry chromatin signatures characteristic of genes with developmental functions. The temporal expression profile of lncRNAs revealed two novel properties: lncRNAs are expressed in narrower time windows than are protein-coding genes and are specifically enriched in early-stage embryos. In addition, several lncRNAs show tissue-specific expression and distinct subcellular localization patterns. Integrative computational analyses associated individual lncRNAs with specific pathways and functions, ranging from cell cycle regulation to morphogenesis. Our study provides the first systematic identification of lncRNAs in a vertebrate embryo and forms the foundation for future genetic, genomic, and evolutionary studies.

The CRISPR/Cas9 system has been implemented in a variety of model organisms to mediate site-directed mutagenesis. A wide range of mutation rates has been reported, but at a limited number of genomic target sites. To uncover the rules that govern effective Cas9-mediated mutagenesis in zebrafish, we targeted over a hundred genomic loci for mutagenesis using a streamlined and cloning-free method. We generated mutations in 85% of target genes with mutation rates varying across several orders of magnitude, and identified sequence composition rules that influence mutagenesis. We increased rates of mutagenesis by implementing several novel approaches. The activities of poor or unsuccessful single-guide RNAs (sgRNAs) initiating with a 5′ adenine were improved by rescuing 5′ end homogeneity of the sgRNA. In some cases, direct injection of Cas9 protein/sgRNA complex further increased mutagenic activity. We also observed that low diversity of mutant alleles led to repeated failure to obtain frame-shift mutations. This limitation was overcome by knock-in of a stop codon cassette that ensured coding frame truncation. Our improved methods and detailed protocols make Cas9-mediated mutagenesis an attractive approach for labs of all sizes.

During embryogenesis, cells acquire distinct fates by transitioning through transcriptional states. To uncover these transcriptional trajectories during zebrafish embryogenesis, we sequenced 38,731 cells and developed URD, a simulated diffusion-based computational reconstruction method. URD identified the trajectories of 25 cell types through early somitogenesis, gene expression along them, and their spatial origin in the blastula. Analysis of Nodal signaling mutants revealed that their transcriptomes were canalized into a subset of wild-type transcriptional trajectories. Some wild-type developmental branchpoints contained cells expressing genes characteristic of multiple fates. These cells appeared to trans-specify from one fate to another. These findings reconstruct the transcriptional trajectories of a vertebrate embryo, highlight the concurrent canalization and plasticity of embryonic specification, and provide a framework to reconstruct complex developmental trees from single-cell transcriptomes.

Cilia, as motile and sensory organelles, have been implicated in normal development, as well as diseases including cystic kidney disease, hydrocephalus and situs inversus. In kidney epithelia, cilia are proposed to be nonmotile sensory organelles, while in the mouse node, two cilia populations, motile and non-motile have been proposed to regulate situs. We show that cilia in the zebrafish larval kidney, the spinal cord and Kupffer's vesicle are motile, suggesting that fluid flow is a common feature of each of these organs. Disruption of cilia structure or motility resulted in pronephric cyst formation, hydrocephalus and left-right asymmetry defects. The data show that loss of fluid flow leads to fluid accumulation, which can account for organ distension pathologies in the kidney and brain. In Kupffer's vesicle, loss of flow is associated with loss of left-right patterning, indicating that the 'nodal flow' mechanism of generating situs is conserved in non-mammalian vertebrates.

The zebrafish EGF-CFC gene one-eyed pinhead (oep) is required zygotically for the formation of the ventral neuroectoderm, endoderm, and prechordal plate. Here we report that embryos lacking both maternal and zygotic Oep activity are defective in germ layer formation, organizer development, and the positioning of the anterior-posterior axis. An identical phenotype is displayed by double mutants for the nodal-related genes squint and cyclops. Mutations in oep eliminate the response to Squint and Cyclops overexpression but are suppressed by expression of Activin and activated forms of the type I receptor ActRIB and Smad2. Expression of the murine EGF-CFC gene cripto rescues oep mutants. These results suggest a conserved role for EGF-CFC proteins as essential extracellular cofactors for Nodal signaling during vertebrate development.