Snake Venom Gland Organoids.

Wnt dependency and Lgr5 expression define multiple mammalian epithelial stem cell types. Under defined growth factor conditions, such adult stem cells (ASCs) grow as 3D organoids that recapitulate essential features of the pertinent epithelium. Here, we establish long-term expanding venom gland organoids from several snake species. The newly assembled transcriptome of the Cape coral snake reveals that organoids express high levels of toxin transcripts. Single-cell RNA sequencing of both organoids and primary tissue identifies distinct venom-expressing cell types as well as proliferative cells expressing homologs of known mammalian stem cell markers. A hard-wired regional heterogeneity in the expression of individual venom components is maintained in organoid cultures. Harvested venom peptides reflect crude venom composition and display biological activity. This study extends organoid technology to reptilian tissues and describes an experimentally tractable model system representing the snake venom gland.

Cell Type Purification by Single-Cell Transcriptome-Trained Sorting.

Much of current molecular and cell biology research relies on the ability to purify cell types by fluorescence-activated cell sorting (FACS). FACS typically relies on the ability to label cell types of interest with antibodies or fluorescent transgenic constructs. However, antibody availability is often limited, and genetic manipulation is labor intensive or impossible in the case of primary human tissue. To date, no systematic method exists to enrich for cell types without a priori knowledge of cell-type markers. Here, we propose GateID, a computational method that combines single-cell transcriptomics with FACS index sorting to purify cell types of choice using only native cellular properties such as cell size, granularity, and mitochondrial content. We validate GateID by purifying various cell types from zebrafish kidney marrow and the human pancreas to high purity without resorting to specific antibodies or transgenes.

Unravelling cellular relationships during development and regeneration using genetic lineage tracing.

Tracking the progeny of single cells is necessary for building lineage trees that recapitulate processes such as embryonic development and stem cell differentiation. In classical lineage tracing experiments, cells are fluorescently labelled to allow identification by microscopy of a limited number of cell clones. To track a larger number of clones in complex tissues, fluorescent proteins are now replaced by heritable DNA barcodes that are read using next-generation sequencing. In prospective lineage tracing, unique DNA barcodes are introduced into single cells through genetic manipulation (using, for example, Cre-mediated recombination or CRISPR-Cas9-mediated editing) and tracked over time. Alternatively, in retrospective lineage tracing, naturally occurring somatic mutations can be used as endogenous DNA barcodes. Finally, single-cell mRNA-sequencing datasets that capture different cell states within a developmental or differentiation trajectory can be used to recapitulate lineages. In this Review, we discuss methods for prospective or retrospective lineage tracing and demonstrate how trajectory reconstruction algorithms can be applied to single-cell mRNA-sequencing datasets to infer developmental or differentiation tracks. We discuss how these approaches are used to understand cell fate during embryogenesis, cell differentiation and tissue regeneration.

Dissecting Immune Circuits by Linking CRISPR-Pooled Screens with Single-Cell RNA-Seq.

In multicellular organisms, dedicated regulatory circuits control cell type diversity and responses. The crosstalk and redundancies within these circuits and substantial cellular heterogeneity pose a major research challenge. Here, we present CRISP-seq, an integrated method for massively parallel single-cell RNA sequencing (RNA-seq) and clustered regularly interspaced short palindromic repeats (CRISPR)-pooled screens. We show that profiling the genomic perturbation and transcriptome in the same cell enables us to simultaneously elucidate the function of multiple factors and their interactions. We applied CRISP-seq to probe regulatory circuits of innate immunity. By sampling tens of thousands of perturbed cells in vitro and in mice, we identified interactions and redundancies between developmental and signaling-dependent factors. These include opposing effects of Cebpb and Irf8 in regulating the monocyte/macrophage versus dendritic cell lineages and differential functions for Rela and Stat1/2 in monocyte versus dendritic cell responses to pathogens. This study establishes CRISP-seq as a broadly applicable, comprehensive, and unbiased approach for elucidating mammalian regulatory circuits.

Design and Analysis of Single-Cell Sequencing Experiments.

Recent advances in single-cell sequencing hold great potential for exploring biological systems with unprecedented resolution. Sequencing the genome of individual cells can reveal somatic mutations and allows the investigation of clonal dynamics. Single-cell transcriptome sequencing can elucidate the cell type composition of a sample. However, single-cell sequencing comes with major technical challenges and yields complex data output. In this Primer, we provide an overview of available methods and discuss experimental design and single-cell data analysis. We hope that these guidelines will enable a growing number of researchers to leverage the power of single-cell sequencing.

Genome-wide maps of nuclear lamina interactions in single human cells.

Mammalian interphase chromosomes interact with the nuclear lamina (NL) through hundreds of large lamina-associated domains (LADs). We report a method to map NL contacts genome-wide in single human cells. Analysis of nearly 400 maps reveals a core architecture consisting of gene-poor LADs that contact the NL with high cell-to-cell consistency, interspersed by LADs with more variable NL interactions. The variable contacts tend to be cell-type specific and are more sensitive to changes in genome ploidy than the consistent contacts. Single-cell maps indicate that NL contacts involve multivalent interactions over hundreds of kilobases. Moreover, we observe extensive intra-chromosomal coordination of NL contacts, even over tens of megabases. Such coordinated loci exhibit preferential interactions as detected by Hi-C. Finally, the consistency of NL contacts is inversely linked to gene activity in single cells and correlates positively with the heterochromatic histone modification H3K9me3. These results highlight fundamental principles of single-cell chromatin organization. VIDEO ABSTRACT.

Prospective derivation of a living organoid biobank of colorectal cancer patients.

In Rspondin-based 3D cultures, Lgr5 stem cells from multiple organs form ever-expanding epithelial organoids that retain their tissue identity. We report the establishment of tumor organoid cultures from 20 consecutive colorectal carcinoma (CRC) patients. For most, organoids were also generated from adjacent normal tissue. Organoids closely recapitulate several properties of the original tumor. The spectrum of genetic changes within the "living biobank" agrees well with previous large-scale mutational analyses of CRC. Gene expression analysis indicates that the major CRC molecular subtypes are represented. Tumor organoids are amenable to high-throughput drug screens allowing detection of gene-drug associations. As an example, a single organoid culture was exquisitely sensitive to Wnt secretion (porcupine) inhibitors and carried a mutation in the negative Wnt feedback regulator RNF43, rather than in APC. Organoid technology may fill the gap between cancer genetics and patient trials, complement cell-line- and xenograft-based drug studies, and allow personalized therapy design. PAPERCLIP.

Every cell is special: genome-wide studies add a new dimension to single-cell biology.

Single-cell analyses have provided invaluable insights into studying heterogenity, signaling, and stochastic gene expression. Recent technological advances now open the door to genome-wide single-cell studies.

Genome-wide RNA Tomography in the zebrafish embryo.

Advancing our understanding of embryonic development is heavily dependent on identification of novel pathways or regulators. Although genome-wide techniques such as RNA sequencing are ideally suited for discovering novel candidate genes, they are unable to yield spatially resolved information in embryos or tissues. Microscopy-based approaches, using in situ hybridization, for example, can provide spatial information about gene expression, but are limited to analyzing one or a few genes at a time. Here, we present a method where we combine traditional histological techniques with low-input RNA sequencing and mathematical image reconstruction to generate a high-resolution genome-wide 3D atlas of gene expression in the zebrafish embryo at three developmental stages. Importantly, our technique enables searching for genes that are expressed in specific spatial patterns without manual image annotation. We envision broad applicability of RNA tomography as an accurate and sensitive approach for spatially resolved transcriptomics in whole embryos and dissected organs.

Studying lineage decision-making in vitro: emerging concepts and novel tools.

Correct and timely lineage decisions are critical for normal embryonic development and homeostasis of adult tissues. Therefore, the search for fundamental principles that underlie lineage decision-making lies at the heart of developmental biology. Here, we review attempts to understand lineage decision-making as the interplay of single-cell heterogeneity and gene regulation. Fluctuations at the single-cell level are an important driving force behind cell-state transitions and the creation of cell-type diversity. Gene regulatory networks amplify such fluctuations and define stable cell types. They also mediate the influence of signaling inputs on the lineage decision. In this review, we focus on insights gleaned from in vitro differentiation of embryonic stem cells. We discuss emerging concepts, with an emphasis on transcriptional regulation, dynamical aspects of differentiation, and functional single-cell heterogeneity. We also highlight some novel tools to study lineage decision-making in vitro.

Feedback control of gene expression variability in the Caenorhabditis elegans Wnt pathway.

Variability in gene expression contributes to phenotypic heterogeneity even in isogenic populations. Here, we used the stereotyped, Wnt signaling-dependent development of the Caenorhabditis elegans Q neuroblast to probe endogenous mechanisms that control gene expression variability. We found that the key Hox gene that orients Q neuroblast migration exhibits increased gene expression variability in mutants in which Wnt pathway activity has been perturbed. Distinct features of the gene expression distributions prompted us on a systematic search for regulatory interactions, revealing a network of interlocked positive and negative feedback loops. Interestingly, positive feedback appeared to cooperate with negative feedback to reduce variability while keeping the Hox gene expression at elevated levels. A minimal model correctly predicts the increased gene expression variability across mutants. Our results highlight the influence of gene network architecture on expression variability and implicate feedback regulation as an effective mechanism to ensure developmental robustness.

Embryo-uterine interaction coordinates mouse embryogenesis during implantation.

Embryo implantation into the uterus marks a key transition in mammalian development. In mice, implantation is mediated by the trophoblast and is accompanied by a morphological transition from the blastocyst to the egg cylinder. However, the roles of trophoblast-uterine interactions in embryo morphogenesis during implantation are poorly understood due to inaccessibility in utero and the remaining challenges to recapitulate it ex vivo from the blastocyst. Here, we engineer a uterus-like microenvironment to recapitulate peri-implantation development of the whole mouse embryo ex vivo and reveal essential roles of the physical embryo-uterine interaction. We demonstrate that adhesion between the trophoblast and the uterine matrix is required for in utero-like transition of the blastocyst to the egg cylinder. Modeling the implanting embryo as a wetting droplet links embryo shape dynamics to the underlying changes in trophoblast adhesion and suggests that the adhesion-mediated tension release facilitates egg cylinder formation. Light-sheet live imaging and the experimental control of the engineered uterine geometry and trophoblast velocity uncovers the coordination between trophoblast motility and embryo growth, where the trophoblast delineates space for embryo morphogenesis.

Quantifying DNA replication speeds in single cells by scEdU-seq.

In a human cell, thousands of replication forks simultaneously coordinate duplication of the entire genome. The rate at which this process occurs might depend on the epigenetic state of the genome and vary between, or even within, cell types. To accurately measure DNA replication speeds, we developed single-cell 5-ethynyl-2'-deoxyuridine sequencing to detect nascent replicated DNA. We observed that the DNA replication speed is not constant but increases during S phase of the cell cycle. Using genetic and pharmacological perturbations we were able to alter this acceleration of replication and conclude that DNA damage inflicted by the process of transcription limits the speed of replication during early S phase. In late S phase, during which less-transcribed regions replicate, replication accelerates and approaches its maximum speed.

The Majority of CD45- Ter119- CD31- Bone Marrow Cell Fraction Is of Hematopoietic Origin and Contains Erythroid and Lymphoid Progenitors.

The non-hematopoietic cell fraction of the bone marrow (BM) is classically identified as CD45- Ter119- CD31- (herein referred to as triple-negative cells or TNCs). Although TNCs are believed to contain heterogeneous stromal cell populations, they remain poorly defined. Here we showed that the vast majority of TNCs (∼85%) have a hematopoietic rather than mesenchymal origin. Single cell RNA-sequencing revealed erythroid and lymphoid progenitor signatures among CD51- TNCs. Ly6D+ CD44+ CD51- TNCs phenotypically and functionally resembled CD45+ pro-B lymphoid cells, whereas Ly6D- CD44+ CD51- TNCs were enriched in previously unappreciated stromal-dependent erythroid progenitors hierarchically situated between preCFU-E and proerythroblasts. Upon adoptive transfer, CD44+ CD51- TNCs contributed to repopulate the B-lymphoid and erythroid compartments. CD44+ CD51- TNCs also expanded during phenylhydrazine-induced acute hemolysis or in a model of sickle cell anemia. These findings thus uncover physiologically relevant new classes of stromal-associated functional CD45- hematopoietic progenitors.

Optimality in the development of intestinal crypts.

Intestinal crypts in mammals are comprised of long-lived stem cells and shorter-lived progenies. These two populations are maintained in specific proportions during adult life. Here, we investigate the design principles governing the dynamics of these proportions during crypt morphogenesis. Using optimal control theory, we show that a proliferation strategy known as a "bang-bang" control minimizes the time to obtain a mature crypt. This strategy consists of a surge of symmetric stem cell divisions, establishing the entire stem cell pool first, followed by a sharp transition to strictly asymmetric stem cell divisions, producing nonstem cells with a delay. We validate these predictions using lineage tracing and single-molecule fluorescence in situ hybridization of intestinal crypts in infant mice, uncovering small crypts that are entirely composed of Lgr5-labeled stem cells, which become a minority as crypts continue to grow. Our approach can be used to uncover similar design principles in other developmental systems.

Transcription of two long noncoding RNAs mediates mating-type control of gametogenesis in budding yeast.

The cell-fate decision leading to gametogenesis is essential for sexual reproduction. In S. cerevisiae, only diploid MATa/α but not haploid MATa or MATα cells undergo gametogenesis, known as sporulation. We find that transcription of two long noncoding RNAs (lncRNAs) mediates mating-type control of sporulation. In MATa or MATα haploids, expression of IME1, the central inducer of gametogenesis, is inhibited in cis by transcription of the lncRNA IRT1, located in the IME1 promoter. IRT1 transcription recruits the Set2 histone methyltransferase and the Set3 histone deacetylase complex to establish repressive chromatin at the IME1 promoter. Inhibiting expression of IRT1 and an antisense transcript that antagonizes the expression of the meiotic regulator IME4 allows cells expressing the haploid mating type to sporulate with kinetics that are indistinguishable from that of MATa/α diploids. Conversely, expression of the two lncRNAs abolishes sporulation in MATa/α diploids. Thus, transcription of two lncRNAs governs mating-type control of gametogenesis in yeast.

Single-cell expression analyses during cellular reprogramming reveal an early stochastic and a late hierarchic phase.

During cellular reprogramming, only a small fraction of cells become induced pluripotent stem cells (iPSCs). Previous analyses of gene expression during reprogramming were based on populations of cells, impeding single-cell level identification of reprogramming events. We utilized two gene expression technologies to profile 48 genes in single cells at various stages during the reprogramming process. Analysis of early stages revealed considerable variation in gene expression between cells in contrast to late stages. Expression of Esrrb, Utf1, Lin28, and Dppa2 is a better predictor for cells to progress into iPSCs than expression of the previously suggested reprogramming markers Fbxo15, Fgf4, and Oct4. Stochastic gene expression early in reprogramming is followed by a late hierarchical phase with Sox2 being the upstream factor in a gene expression hierarchy. Finally, downstream factors derived from the late phase, which do not include Oct4, Sox2, Klf4, c-Myc, and Nanog, can activate the pluripotency circuitry.

Slug and Sox9 cooperatively determine the mammary stem cell state.

Regulatory networks orchestrated by key transcription factors (TFs) have been proposed to play a central role in the determination of stem cell states. However, the master transcriptional regulators of adult stem cells are poorly understood. We have identified two TFs, Slug and Sox9, that act cooperatively to determine the mammary stem cell (MaSC) state. Inhibition of either Slug or Sox9 blocks MaSC activity in primary mammary epithelial cells. Conversely, transient coexpression of exogenous Slug and Sox9 suffices to convert differentiated luminal cells into MaSCs with long-term mammary gland-reconstituting ability. Slug and Sox9 induce MaSCs by activating distinct autoregulatory gene expression programs. We also show that coexpression of Slug and Sox9 promotes the tumorigenic and metastasis-seeding abilities of human breast cancer cells and is associated with poor patient survival, providing direct evidence that human breast cancer stem cells are controlled by key regulators similar to those operating in normal murine MaSCs.

Single-cell Ribo-seq reveals cell cycle-dependent translational pausing.

Single-cell sequencing methods have enabled in-depth analysis of the diversity of cell types and cell states in a wide range of organisms. These tools focus predominantly on sequencing the genomes1, epigenomes2 and transcriptomes3 of single cells. However, despite recent progress in detecting proteins by mass spectrometry with single-cell resolution4, it remains a major challenge to measure translation in individual cells. Here, building on existing protocols5-7, we have substantially increased the sensitivity of these assays to enable ribosome profiling in single cells. Integrated with a machine learning approach, this technology achieves single-codon resolution. We validate this method by demonstrating that limitation for a particular amino acid causes ribosome pausing at a subset of the codons encoding the amino acid. Of note, this pausing is only observed in a sub-population of cells correlating to its cell cycle state. We further expand on this phenomenon in non-limiting conditions and detect pronounced GAA pausing during mitosis. Finally, we demonstrate the applicability of this technique to rare primary enteroendocrine cells. This technology provides a first step towards determining the contribution of the translational process to the remarkable diversity between seemingly identical cells.

Cellular decision making and biological noise: from microbes to mammals.

Cellular decision making is the process whereby cells assume different, functionally important and heritable fates without an associated genetic or environmental difference. Such stochastic cell fate decisions generate nongenetic cellular diversity, which may be critical for metazoan development as well as optimized microbial resource utilization and survival in a fluctuating, frequently stressful environment. Here, we review several examples of cellular decision making from viruses, bacteria, yeast, lower metazoans, and mammals, highlighting the role of regulatory network structure and molecular noise. We propose that cellular decision making is one of at least three key processes underlying development at various scales of biological organization.