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.

Postnatal Cardiac Gene Editing Using CRISPR/Cas9 With AAV9-Mediated Delivery of Short Guide RNAs Results in Mosaic Gene Disruption.

RATIONALE: CRISPR/Cas9 (clustered regularly interspaced palindromic repeats/CRISPR-associated protein 9)-based DNA editing has rapidly evolved as an attractive tool to modify the genome. Although CRISPR/Cas9 has been extensively used to manipulate the germline in zygotes, its application in postnatal gene editing remains incompletely characterized. OBJECTIVE: To evaluate the feasibility of CRISPR/Cas9-based cardiac genome editing in vivo in postnatal mice. METHODS AND RESULTS: We generated cardiomyocyte-specific Cas9 mice and demonstrated that Cas9 expression does not affect cardiac function or gene expression. As a proof-of-concept, we delivered short guide RNAs targeting 3 genes critical for cardiac physiology, Myh6, Sav1, and Tbx20, using a cardiotropic adeno-associated viral vector 9. Despite a similar degree of DNA disruption and subsequent mRNA downregulation, only disruption of Myh6 was sufficient to induce a cardiac phenotype, irrespective of short guide RNA exposure or the level of Cas9 expression. DNA sequencing analysis revealed target-dependent mutations that were highly reproducible across mice resulting in differential rates of in- and out-of-frame mutations. Finally, we applied a dual short guide RNA approach to effectively delete an important coding region of Sav1, which increased the editing efficiency. CONCLUSIONS: Our results indicate that the effect of postnatal CRISPR/Cas9-based cardiac gene editing using adeno-associated virus serotype 9 to deliver a single short guide RNA is target dependent. We demonstrate a mosaic pattern of gene disruption, which hinders the application of the technology to study gene function. Further studies are required to expand the versatility of CRISPR/Cas9 as a robust tool to study novel cardiac gene functions in vivo.

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.

Passenger Gene Coamplifications Create Collateral Therapeutic Vulnerabilities in Cancer.

DNA amplifications in cancer do not only harbor oncogenes. We sought to determine whether passenger coamplifications could create collateral therapeutic vulnerabilities. Through an analysis of >3,000 cancer genomes followed by the interrogation of CRISPR-Cas9 loss-of-function screens across >700 cancer cell lines, we determined that passenger coamplifications are accompanied by distinct dependency profiles. In a proof-of-principle study, we demonstrate that the coamplification of the bona fide passenger gene DEAD-Box Helicase 1 (DDX1) creates an increased dependency on the mTOR pathway. Interaction proteomics identified tricarboxylic acid (TCA) cycle components as previously unrecognized DDX1 interaction partners. Live-cell metabolomics highlighted that this interaction could impair TCA activity, which in turn resulted in enhanced mTORC1 activity. Consequently, genetic and pharmacologic disruption of mTORC1 resulted in pronounced cell death in vitro and in vivo. Thus, structurally linked coamplification of a passenger gene and an oncogene can result in collateral vulnerabilities. SIGNIFICANCE: We demonstrate that coamplification of passenger genes, which were largely neglected in cancer biology in the past, can create distinct cancer dependencies. Because passenger coamplifications are frequent in cancer, this principle has the potential to expand target discovery in oncology. This article is featured in Selected Articles from This Issue, p. 384.

Origin and function of activated fibroblast states during zebrafish heart regeneration.

The adult zebrafish heart has a high capacity for regeneration following injury. However, the composition of the regenerative niche has remained largely elusive. Here, we dissected the diversity of activated cell states in the regenerating zebrafish heart based on single-cell transcriptomics and spatiotemporal analysis. We observed the emergence of several transient cell states with fibroblast characteristics following injury, and we outlined the proregenerative function of collagen-12-expressing fibroblasts. To understand the cascade of events leading to heart regeneration, we determined the origin of these cell states by high-throughput lineage tracing. We found that activated fibroblasts were derived from two separate sources: the epicardium and the endocardium. Mechanistically, we determined Wnt signalling as a regulator of the endocardial fibroblast response. In summary, our work identifies specialized activated fibroblast cell states that contribute to heart regeneration, thereby opening up possible approaches to modulating the regenerative capacity of the vertebrate heart.

Spatio-temporal mRNA tracking in the early zebrafish embryo.

Early stages of embryogenesis depend on subcellular localization and transport of maternal mRNA. However, systematic analysis of these processes is hindered by a lack of spatio-temporal information in single-cell RNA sequencing. Here, we combine spatially-resolved transcriptomics and single-cell RNA labeling to perform a spatio-temporal analysis of the transcriptome during early zebrafish development. We measure spatial localization of mRNA molecules within the one-cell stage embryo, which allows us to identify a class of mRNAs that are specifically localized at an extraembryonic position, the vegetal pole. Furthermore, we establish a method for high-throughput single-cell RNA labeling in early zebrafish embryos, which enables us to follow the fate of individual maternal transcripts until gastrulation. This approach reveals that many localized transcripts are specifically transported to the primordial germ cells. Finally, we acquire spatial transcriptomes of two xenopus species and compare evolutionary conservation of localized genes as well as enriched sequence motifs.

GPSeq reveals the radial organization of chromatin in the cell nucleus.

With the exception of lamina-associated domains, the radial organization of chromatin in mammalian cells remains largely unexplored. Here we describe genomic loci positioning by sequencing (GPSeq), a genome-wide method for inferring distances to the nuclear lamina all along the nuclear radius. GPSeq relies on gradual restriction digestion of chromatin from the nuclear lamina toward the nucleus center, followed by sequencing of the generated cut sites. Using GPSeq, we mapped the radial organization of the human genome at 100-kb resolution, which revealed radial patterns of genomic and epigenomic features and gene expression, as well as A and B subcompartments. By combining radial information with chromosome contact frequencies measured by Hi-C, we substantially improved the accuracy of whole-genome structure modeling. Finally, we charted the radial topography of DNA double-strand breaks, germline variants and cancer mutations and found that they have distinctive radial arrangements in A and B subcompartments. We conclude that GPSeq can reveal fundamental aspects of genome architecture.

Simultaneous lineage tracing and cell-type identification using CRISPR-Cas9-induced genetic scars.

A key goal of developmental biology is to understand how a single cell is transformed into a full-grown organism comprising many different cell types. Single-cell RNA-sequencing (scRNA-seq) is commonly used to identify cell types in a tissue or organ. However, organizing the resulting taxonomy of cell types into lineage trees to understand the developmental origin of cells remains challenging. Here we present LINNAEUS (lineage tracing by nuclease-activated editing of ubiquitous sequences)-a strategy for simultaneous lineage tracing and transcriptome profiling in thousands of single cells. By combining scRNA-seq with computational analysis of lineage barcodes, generated by genome editing of transgenic reporter genes, we reconstruct developmental lineage trees in zebrafish larvae, and in heart, liver, pancreas, and telencephalon of adult fish. LINNAEUS provides a systematic approach for tracing the origin of novel cell types, or known cell types under different conditions.

Spatial Transcriptomics of C. elegans Males and Hermaphrodites Identifies Sex-Specific Differences in Gene Expression Patterns.

To advance our understanding of the genetic programs that drive cell and tissue specialization, it is necessary to obtain a comprehensive overview of gene expression patterns. Here, we have used spatial transcriptomics to generate high-resolution, anteroposterior gene expression maps of C. elegans males and hermaphrodites. To explore these maps, we have developed computational methods for discovering region- and tissue-specific genes. We have found extensive sex-specific gene expression differences in the germline and sperm and discovered genes that are specifically expressed in the male reproductive tract. These include a group of uncharacterized genes that encode small secreted proteins that are required for male fertility. We conclude that spatial gene expression maps provide a powerful resource for identifying tissue-specific gene functions in C. elegans. Importantly, we found that expression maps from different animals can be precisely aligned, enabling transcriptome-wide comparisons of gene expression patterns.

Methods for lineage tracing on the organism-wide level.

Determining the lineage origin of cell types is a major goal in developmental biology. Furthermore, lineage tracing is a powerful approach for understanding the origin of developmental defects as well as the origin of diseases such as cancer. There is now a variety of complementary approaches for identifying lineage relationships, ranging from direct observation of cell divisions by light microscopy to genetic labeling of cells using inducible recombinases and fluorescent reporters. A recent development, and the main topic of this review article, is the use of high-throughput sequencing data for lineage analysis. This emerging approach holds the promise of increased multiplexing capacity, allowing lineage analysis of large cell numbers up to the organism-wide level combined with simultaneous transcription profiling by single cell RNA sequencing.

Integrated genome and transcriptome sequencing of the same cell.

Single-cell genomics and single-cell transcriptomics have emerged as powerful tools to study the biology of single cells at a genome-wide scale. However, a major challenge is to sequence both genomic DNA and mRNA from the same cell, which would allow direct comparison of genomic variation and transcriptome heterogeneity. We describe a quasilinear amplification strategy to quantify genomic DNA and mRNA from the same cell without physically separating the nucleic acids before amplification. We show that the efficiency of our integrated approach is similar to existing methods for single-cell sequencing of either genomic DNA or mRNA. Further, we find that genes with high cell-to-cell variability in transcript numbers generally have lower genomic copy numbers, and vice versa, suggesting that copy number variations may drive variability in gene expression among individual cells. Applications of our integrated sequencing approach could range from gaining insights into cancer evolution and heterogeneity to understanding the transcriptional consequences of copy number variations in healthy and diseased tissues.