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.

A single-cell atlas of de novo ß-cell regeneration reveals the contribution of hybrid ß/δ-cells to diabetes recovery in zebrafish.

Regeneration-competent species possess the ability to reverse the progression of severe diseases by restoring the function of the damaged tissue. However, the cellular dynamics underlying this capability remain unexplored. Here, we have used single-cell transcriptomics to map de novo ß-cell regeneration during induction and recovery from diabetes in zebrafish. We show that the zebrafish has evolved two distinct types of somatostatin-producing δ-cells, which we term δ1- and δ2-cells. Moreover, we characterize a small population of glucose-responsive islet cells, which share the hormones and fate-determinants of both ß- and δ1-cells. The transcriptomic analysis of ß-cell regeneration reveals that ß/δ hybrid cells provide a prominent source of insulin expression during diabetes recovery. Using in vivo calcium imaging and cell tracking, we further show that the hybrid cells form de novo and acquire glucose-responsiveness in the course of regeneration. The overexpression of dkk3, a gene enriched in hybrid cells, increases their formation in the absence of ß-cell injury. Finally, interspecies comparison shows that plastic δ1-cells are partially related to PP cells in the human pancreas. Our work provides an atlas of ß-cell regeneration and indicates that the rapid formation of glucose-responsive hybrid cells contributes to the resolution of diabetes in zebrafish.

Dissecting the spatiotemporal diversity of adult neural stem cells.

Adult stem cells are important for tissue turnover and regeneration. However, in most adult systems it remains elusive how stem cells assume different functional states and support spatially patterned tissue architecture. Here, we dissected the diversity of neural stem cells in the adult zebrafish brain, an organ that is characterized by pronounced zonation and high regenerative capacity. We combined single-cell transcriptomics of dissected brain regions with massively parallel lineage tracing and in vivo RNA metabolic labeling to analyze the regulation of neural stem cells in space and time. We detected a large diversity of neural stem cells, with some subtypes being restricted to a single brain region, while others were found globally across the brain. Global stem cell states are linked to neurogenic differentiation, with different states being involved in proliferative and non-proliferative differentiation. Our work reveals principles of adult stem cell organization and establishes a resource for the functional manipulation of neural stem cell subtypes.

A single-cell RNA labeling strategy for measuring stress response upon tissue dissociation.

Tissue dissociation, a crucial step in single-cell sample preparation, can alter the transcriptional state of a sample through the intrinsic cellular stress response. Here we demonstrate a general approach for measuring transcriptional response during sample preparation. In our method, transcripts made during dissociation are labeled for later identification upon sequencing. We found general as well as cell-type-specific dissociation response programs in zebrafish larvae, and we observed sample-to-sample variation in the dissociation response of mouse cardiomyocytes despite well-controlled experimental conditions. Finally, we showed that dissociation of the mouse hippocampus can lead to the artificial activation of microglia. In summary, our approach facilitates experimental optimization of dissociation procedures as well as computational removal of transcriptional perturbation response.

Single-cell transcriptomics enters the age of mass production.

Two publications in the current issue of Cell introduce novel methods for high-throughput single-cell transcriptomics by using droplet microfluidics and sophisticated barcoding schemes for transcriptional profiling of thousands of individual cells.

Systematic identification of A-to-I RNA editing in zebrafish development and adult organs.

A-to-I RNA editing is a common post transcriptional mechanism, mediated by the Adenosine deaminase that acts on RNA (ADAR) enzymes, that increases transcript and protein diversity. The study of RNA editing is limited by the absence of editing maps for most model organisms, hindering the understanding of its impact on various physiological conditions. Here, we mapped the vertebrate developmental landscape of A-to-I RNA editing, and generated the first comprehensive atlas of editing sites in zebrafish. Tens of thousands unique editing events and 149 coding sites were identified with high-accuracy. Some of these edited sites are conserved between zebrafish and humans. Sequence analysis of RNA over seven developmental stages revealed high levels of editing activity in early stages of embryogenesis, when embryos rely on maternal mRNAs and proteins. In contrast to the other organisms studied so far, the highest levels of editing were detected in the zebrafish ovary and testes. This resource can serve as the basis for understanding of the role of editing during zebrafish development and maturity.

Long-term imaging reveals behavioral plasticity during C. elegans dauer exit.

BACKGROUND: During their lifetime, animals must adapt their behavior to survive in changing environments. This ability requires the nervous system to undergo adjustments at distinct temporal scales, from short-term dynamic changes in expression of neurotransmitters and receptors to longer-term growth, spatial and connectivity reorganization, while integrating external stimuli. The nematode Caenorhabditis elegans provides a model of nervous system plasticity, in particular its dauer exit decision. Under unfavorable conditions, larvae will enter the non-feeding and non-reproductive stress-resistant dauer stage and adapt their behavior to cope with the harsh new environment, with active reversal under improved conditions leading to resumption of reproductive development. However, how different environmental stimuli regulate the exit decision mechanism and thereby drive the larva's behavioral change is unknown. To fill this gap and provide insights on behavioral changes over extended periods of time, we developed a new open hardware method for long-term imaging (12h) of C. elegans larvae. RESULTS: Our WormObserver platform comprises open hardware and software components for video acquisition, automated processing of large image data (> 80k images/experiment) and data analysis. We identified dauer-specific behavioral motifs and characterized the behavioral trajectory of dauer exit in different environments and genetic backgrounds to identify key decision points and stimuli promoting dauer exit. Combining long-term behavioral imaging with transcriptomics data, we find that bacterial ingestion triggers a change in neuropeptide gene expression to establish post-dauer behavior. CONCLUSIONS: Taken together, we show how a developing nervous system can robustly integrate environmental changes activate a developmental switch and adapt the organism's behavior to a new environment. WormObserver is generally applicable to other research questions within and beyond the C. elegans field, having a modular and customizable character and allowing assessment of behavioral plasticity over longer periods.

Multispecies RNA tomography reveals regulators of hematopoietic stem cell birth in the embryonic aorta.

The defined location of a stem cell within a niche regulates its fate, behavior, and molecular identity via a complex extrinsic regulation that is far from being fully elucidated. To explore the molecular characteristics and key components of the aortic microenvironment, where the first hematopoietic stem cells are generated during development, we performed genome-wide RNA tomography sequencing on zebrafish, chicken, mouse, and human embryos. The resulting anterior-posterior and dorsal-ventral transcriptional maps provided a powerful resource for exploring genes and regulatory pathways active in the aortic microenvironment. By performing interspecies comparative RNA sequencing analyses and functional assays, we explored the complexity of the aortic microenvironment landscape and the fine-tuning of various factors interacting to control hematopoietic stem cell generation, both in time and space in vivo, including the ligand-receptor couple ADM-RAMP2 and SVEP1. Understanding the regulatory function of the local environment will pave the way for improved stem cell production in vitro and clinical cell therapy.

The transcriptome dynamics of single cells during the cell cycle.

The cell cycle is among the most basic phenomena in biology. Despite advances in single-cell analysis, dynamics and topology of the cell cycle in high-dimensional gene expression space remain largely unknown. We developed a linear analysis of transcriptome data which reveals that cells move along a planar circular trajectory in transcriptome space during the cycle. Non-cycling gene expression adds a third dimension causing helical motion on a cylinder. We find in immortalized cell lines that cell cycle transcriptome dynamics occur largely independently from other cellular processes. We offer a simple method ("Revelio") to order unsynchronized cells in time. Precise removal of cell cycle effects from the data becomes a straightforward operation. The shape of the trajectory implies that each gene is upregulated only once during the cycle, and only two dynamic components represented by groups of genes drive transcriptome dynamics. It indicates that the cell cycle has evolved to minimize changes of transcriptional activity and the related regulatory effort. This design principle of the cell cycle may be of relevance to many other cellular differentiation processes.

When noisy neighbors are a blessing: analysis of gene expression noise identifies coregulated genes.

In this issue of Molecular Cell, Stewart-Ornstein et al. (2012) use systematic pair-wise correlation analysis of expression noise in a large number of yeast genes to identify clusters of functionally related genes and signaling pathways responsible for elevated noise.

Neural-specific Sox2 input and differential Gli-binding affinity provide context and positional information in Shh-directed neural patterning.

In the vertebrate neural tube, regional Sonic hedgehog (Shh) signaling invokes a time- and concentration-dependent induction of six different cell populations mediated through Gli transcriptional regulators. Elsewhere in the embryo, Shh/Gli responses invoke different tissue-appropriate regulatory programs. A genome-scale analysis of DNA binding by Gli1 and Sox2, a pan-neural determinant, identified a set of shared regulatory regions associated with key factors central to cell fate determination and neural tube patterning. Functional analysis in transgenic mice validates core enhancers for each of these factors and demonstrates the dual requirement for Gli1 and Sox2 inputs for neural enhancer activity. Furthermore, through an unbiased determination of Gli-binding site preferences and analysis of binding site variants in the developing mammalian CNS, we demonstrate that differential Gli-binding affinity underlies threshold-level activator responses to Shh input. In summary, our results highlight Sox2 input as a context-specific determinant of the neural-specific Shh response and differential Gli-binding site affinity as an important cis-regulatory property critical for interpreting Shh morphogen action in the mammalian neural tube.

Tomo-Seq Identifies SOX9 as a Key Regulator of Cardiac Fibrosis During Ischemic Injury.

BACKGROUND: Cardiac ischemic injury induces a pathological remodeling response, which can ultimately lead to heart failure. Detailed mechanistic insights into molecular signaling pathways relevant for different aspects of cardiac remodeling will support the identification of novel therapeutic targets. METHODS: Although genome-wide transcriptome analysis on diseased tissues has greatly advanced our understanding of the regulatory networks that drive pathological changes in the heart, this approach has been disadvantaged by the fact that the signals are derived from tissue homogenates. Here we used tomo-seq to obtain a genome-wide gene expression signature with high spatial resolution spanning from the infarcted area to the remote to identify new regulators of cardiac remodeling. Cardiac tissue samples from patients suffering from ischemic heart disease were used to validate our findings. RESULTS: Tracing transcriptional differences with a high spatial resolution across the infarcted heart enabled us to identify gene clusters that share a comparable expression profile. The spatial distribution patterns indicated a separation of expressional changes for genes involved in specific aspects of cardiac remodeling, such as fibrosis, cardiomyocyte hypertrophy, and calcium handling (Col1a2, Nppa, and Serca2). Subsequent correlation analysis allowed for the identification of novel factors that share a comparable transcriptional regulation pattern across the infarcted tissue. The strong correlation between the expression levels of these known marker genes and the expression of the coregulated genes could be confirmed in human ischemic cardiac tissue samples. Follow-up analysis identified SOX9 as common transcriptional regulator of a large portion of the fibrosis-related genes that become activated under conditions of ischemic injury. Lineage-tracing experiments indicated that the majority of COL1-positive fibroblasts stem from a pool of SOX9-expressing cells, and in vivo loss of Sox9 blunted the cardiac fibrotic response on ischemic injury. The colocalization between SOX9 and COL1 could also be confirmed in patients suffering from ischemic heart disease. CONCLUSIONS: Based on the exact local expression cues, tomo-seq can serve to reveal novel genes and key transcription factors involved in specific aspects of cardiac remodeling. Using tomo-seq, we were able to unveil the unknown relevance of SOX9 as a key regulator of cardiac fibrosis, pointing to SOX9 as a potential therapeutic target for cardiac fibrosis.

A Gene Regulatory Program for Meiotic Prophase in the Fetal Ovary.

The chromosomal program of meiotic prophase, comprising events such as laying down of meiotic cohesins, synapsis between homologs, and homologous recombination, must be preceded and enabled by the regulated induction of meiotic prophase genes. This gene regulatory program is poorly understood, particularly in organisms with a segregated germline. We characterized the gene regulatory program of meiotic prophase as it occurs in the mouse fetal ovary. By profiling gene expression in the mouse fetal ovary in mutants with whole tissue and single-cell techniques, we identified 104 genes expressed specifically in pre-meiotic to pachytene germ cells. We characterized the regulation of these genes by 1) retinoic acid (RA), which induces meiosis, 2) Dazl, which is required for germ cell competence to respond to RA, and 3) Stra8, a downstream target of RA required for the chromosomal program of meiotic prophase. Initial induction of practically all identified meiotic prophase genes requires Dazl. In the presence of Dazl, RA induces at least two pathways: one Stra8-independent, and one Stra8-dependent. Genes vary in their induction by Stra8, spanning fully Stra8-independent, partially Stra8-independent, and fully Stra8-dependent. Thus, Stra8 regulates the entirety of the chromosomal program but plays a more nuanced role in governing the gene expression program. We propose that Stra8-independent gene expression enables the stockpiling of selected meiotic structural proteins prior to the commencement of the chromosomal program. Unexpectedly, we discovered that Stra8 is required for prompt down-regulation of itself and Rec8. Germ cells that have expressed and down-regulated Stra8 are refractory to further Stra8 expression. Negative feedback of Stra8, and subsequent resistance to further Stra8 expression, may ensure a single, restricted pulse of Stra8 expression. Collectively, our findings reveal a gene regulatory logic by which germ cells prepare for the chromosomal program of meiotic prophase, and ensure that it is induced only once.

Licensing of primordial germ cells for gametogenesis depends on genital ridge signaling.

In mouse embryos at mid-gestation, primordial germ cells (PGCs) undergo licensing to become gametogenesis-competent cells (GCCs), gaining the capacity for meiotic initiation and sexual differentiation. GCCs then initiate either oogenesis or spermatogenesis in response to gonadal cues. Germ cell licensing has been considered to be a cell-autonomous and gonad-independent event, based on observations that some PGCs, having migrated not to the gonad but to the adrenal gland, nonetheless enter meiosis in a time frame parallel to ovarian germ cells -- and do so regardless of the sex of the embryo. Here we test the hypothesis that germ cell licensing is cell-autonomous by examining the fate of PGCs in Gata4 conditional mutant (Gata4 cKO) mouse embryos. Gata4, which is expressed only in somatic cells, is known to be required for genital ridge initiation. PGCs in Gata4 cKO mutants migrated to the area where the genital ridge, the precursor of the gonad, would ordinarily be formed. However, these germ cells did not undergo licensing and instead retained characteristics of PGCs. Our results indicate that licensing is not purely cell-autonomous but is induced by the somatic genital ridge.

DAZL regulates Tet1 translation in murine embryonic stem cells.

Embryonic stem cell (ESC) cultures display a heterogeneous gene expression profile, ranging from a pristine naïve pluripotent state to a primed epiblast state. Addition of inhibitors of GSK3ß and MEK (so-called 2i conditions) pushes ESC cultures toward a more homogeneous naïve pluripotent state, but the molecular underpinnings of this naïve transition are not completely understood. Here, we demonstrate that DAZL, an RNA-binding protein known to play a key role in germ-cell development, marks a subpopulation of ESCs that is actively transitioning toward naïve pluripotency. Moreover, DAZL plays an essential role in the active reprogramming of cytosine methylation. We demonstrate that DAZL associates with mRNA of Tet1, a catalyst of 5-hydroxylation of methyl-cytosine, and enhances Tet1 mRNA translation. Overexpression of DAZL in heterogeneous ESC cultures results in elevated TET1 protein levels as well as increased global hydroxymethylation. Conversely, null mutation of Dazl severely stunts 2i-mediated TET1 induction and hydroxymethylation. Our results provide insight into the regulation of the acquisition of naïve pluripotency and demonstrate that DAZL enhances TET1-mediated cytosine hydroxymethylation in ESCs that are actively reprogramming to a pluripotent ground state.

Inference of differentiation trajectories by transfer learning across biological processes.

Stem cells differentiate into distinct fates by transitioning through a series of transcriptional states. Current computational approaches allow reconstruction of differentiation trajectories from single-cell transcriptomics data, but it remains unknown to what degree differentiation can be predicted across biological processes. Here, we use transfer learning to infer differentiation processes and quantify predictability in early embryonic development and adult hematopoiesis. Overall, we find that non-linear methods outperform linear approaches, and we achieved the best predictions with a custom variational autoencoder that explicitly models changes in transcriptional variance. We observed a high accuracy of predictions in embryonic development, but we found somewhat lower agreement with the real data in adult hematopoiesis. We demonstrate that this discrepancy can be explained by a higher degree of concordant transcriptional processes along embryonic differentiation compared with adult homeostasis. In summary, we establish a framework for quantifying and exploiting predictability of cellular differentiation trajectories.

Antigen presentation plays positive roles in the regenerative response to cardiac injury in zebrafish.

In contrast to adult mammals, adult zebrafish can fully regenerate injured cardiac tissue, and this regeneration process requires an adequate and tightly controlled immune response. However, which components of the immune response are required during regeneration is unclear. Here, we report positive roles for the antigen presentation-adaptive immunity axis during zebrafish cardiac regeneration. We find that following the initial innate immune response, activated endocardial cells (EdCs), as well as immune cells, start expressing antigen presentation genes. We also observe that T helper cells, a.k.a. Cd4+ T cells, lie in close physical proximity to these antigen-presenting EdCs. We targeted Major Histocompatibility Complex (MHC) class II antigen presentation by generating cd74a; cd74b mutants, which display a defective immune response. In these mutants, Cd4+ T cells and activated EdCs fail to efficiently populate the injured tissue and EdC proliferation is significantly decreased. cd74a; cd74b mutants exhibit additional defects in cardiac regeneration including reduced cardiomyocyte dedifferentiation and proliferation. Notably, Cd74 also becomes activated in neonatal mouse EdCs following cardiac injury. Altogether, these findings point to positive roles for antigen presentation during cardiac regeneration, potentially involving interactions between activated EdCs, classical antigen-presenting cells, and Cd4+ T cells.

Spatial transcriptomics unveils ZBTB11 as a regulator of cardiomyocyte degeneration in arrhythmogenic cardiomyopathy.

AIMS: Arrhythmogenic cardiomyopathy (ACM) is an inherited cardiac disorder that is characterized by progressive loss of myocardium that is replaced by fibro-fatty cells, arrhythmias, and sudden cardiac death. While myocardial degeneration and fibro-fatty replacement occur in specific locations, the underlying molecular changes remain poorly characterized. Here, we aim to delineate local changes in gene expression to identify new genes and pathways that are relevant for specific remodelling processes occurring during ACM. METHODS AND RESULTS: Using Tomo-Seq, genome-wide transcriptional profiling with high spatial resolution, we created transmural epicardial-to-endocardial gene expression atlases of explanted ACM hearts to gain molecular insights into disease-driving processes. This enabled us to link gene expression profiles to the different regional remodelling responses and allowed us to identify genes that are potentially relevant for disease progression. In doing so, we identified distinct gene expression profiles marking regions of cardiomyocyte degeneration and fibro-fatty remodelling and revealed Zinc finger and BTB domain-containing protein 11 (ZBTB11) to be specifically enriched at sites of active fibro-fatty replacement of myocardium. Immunohistochemistry indicated ZBTB11 to be induced in cardiomyocytes flanking fibro-fatty areas, which could be confirmed in multiple cardiomyopathy patients. Forced overexpression of ZBTB11 induced autophagy and cell death-related gene programmes in human cardiomyocytes, leading to increased apoptosis. CONCLUSION: Our study shows the power of Tomo-Seq to unveil new molecular mechanisms in human cardiomyopathy and uncovers ZBTB11 as a novel driver of cardiomyocyte loss.