Charting human development using a multi-endodermal organ atlas and organoid models.

Organs are composed of diverse cell types that traverse transient states during organogenesis. To interrogate this diversity during human development, we generate a single-cell transcriptome atlas from multiple developing endodermal organs of the respiratory and gastrointestinal tract. We illuminate cell states, transcription factors, and organ-specific epithelial stem cell and mesenchyme interactions across lineages. We implement the atlas as a high-dimensional search space to benchmark human pluripotent stem cell (hPSC)-derived intestinal organoids (HIOs) under multiple culture conditions. We show that HIOs recapitulate reference cell states and use HIOs to reconstruct the molecular dynamics of intestinal epithelium and mesenchyme emergence. We show that the mesenchyme-derived niche cue NRG1 enhances intestinal stem cell maturation in vitro and that the homeobox transcription factor CDX2 is required for regionalization of intestinal epithelium and mesenchyme in humans. This work combines cell atlases and organoid technologies to understand how human organ development is orchestrated.

Inferring and perturbing cell fate regulomes in human brain organoids.

Self-organizing neural organoids grown from pluripotent stem cells1-3 combined with single-cell genomic technologies provide opportunities to examine gene regulatory networks underlying human brain development. Here we acquire single-cell transcriptome and accessible chromatin data over a dense time course in human organoids covering neuroepithelial formation, patterning, brain regionalization and neurogenesis, and identify temporally dynamic and brain-region-specific regulatory regions. We developed Pando-a flexible framework that incorporates multi-omic data and predictions of transcription-factor-binding sites to infer a global gene regulatory network describing organoid development. We use pooled genetic perturbation with single-cell transcriptome readout to assess transcription factor requirement for cell fate and state regulation in organoids. We find that certain factors regulate the abundance of cell fates, whereas other factors affect neuronal cell states after differentiation. We show that the transcription factor GLI3 is required for cortical fate establishment in humans, recapitulating previous research performed in mammalian model systems. We measure transcriptome and chromatin accessibility in normal or GLI3-perturbed cells and identify two distinct GLI3 regulomes that are central to telencephalic fate decisions: one regulating dorsoventral patterning with HES4/5 as direct GLI3 targets, and one controlling ganglionic eminence diversification later in development. Together, we provide a framework for how human model systems and single-cell technologies can be leveraged to reconstruct human developmental biology.

Single-cell brain organoid screening identifies developmental defects in autism.

The development of the human brain involves unique processes (not observed in many other species) that can contribute to neurodevelopmental disorders1-4. Cerebral organoids enable the study of neurodevelopmental disorders in a human context. We have developed the CRISPR-human organoids-single-cell RNA sequencing (CHOOSE) system, which uses verified pairs of guide RNAs, inducible CRISPR-Cas9-based genetic disruption and single-cell transcriptomics for pooled loss-of-function screening in mosaic organoids. Here we show that perturbation of 36 high-risk autism spectrum disorder genes related to transcriptional regulation uncovers their effects on cell fate determination. We find that dorsal intermediate progenitors, ventral progenitors and upper-layer excitatory neurons are among the most vulnerable cell types. We construct a developmental gene regulatory network of cerebral organoids from single-cell transcriptomes and chromatin modalities and identify autism spectrum disorder-associated and perturbation-enriched regulatory modules. Perturbing members of the BRG1/BRM-associated factor (BAF) chromatin remodelling complex leads to enrichment of ventral telencephalon progenitors. Specifically, mutating the BAF subunit ARID1B affects the fate transition of progenitors to oligodendrocyte and interneuron precursor cells, a phenotype that we confirmed in patient-specific induced pluripotent stem cell-derived organoids. Our study paves the way for high-throughput phenotypic characterization of disease susceptibility genes in organoid models with cell state, molecular pathway and gene regulatory network readouts.

A nomenclature consensus for nervous system organoids and assembloids.

Self-organizing three-dimensional cellular models derived from human pluripotent stem cells or primary tissue have great potential to provide insights into how the human nervous system develops, what makes it unique and how disorders of the nervous system arise, progress and could be treated. Here, to facilitate progress and improve communication with the scientific community and the public, we clarify and provide a basic framework for the nomenclature of human multicellular models of nervous system development and disease, including organoids, assembloids and transplants.

A roadmap for the Human Developmental Cell Atlas.

The Human Developmental Cell Atlas (HDCA) initiative, which is part of the Human Cell Atlas, aims to create a comprehensive reference map of cells during development. This will be critical to understanding normal organogenesis, the effect of mutations, environmental factors and infectious agents on human development, congenital and childhood disorders, and the cellular basis of ageing, cancer and regenerative medicine. Here we outline the HDCA initiative and the challenges of mapping and modelling human development using state-of-the-art technologies to create a reference atlas across gestation. Similar to the Human Genome Project, the HDCA will integrate the output from a growing community of scientists who are mapping human development into a unified atlas. We describe the early milestones that have been achieved and the use of human stem-cell-derived cultures, organoids and animal models to inform the HDCA, especially for prenatal tissues that are hard to acquire. Finally, we provide a roadmap towards a complete atlas of human development.

The road ahead in genetics and genomics.

In celebration of the 20th anniversary of Nature Reviews Genetics, we asked 12 leading researchers to reflect on the key challenges and opportunities faced by the field of genetics and genomics. Keeping their particular research area in mind, they take stock of the current state of play and emphasize the work that remains to be done over the next few years so that, ultimately, the benefits of genetic and genomic research can be felt by everyone.

Lineage recording in human cerebral organoids.

Induced pluripotent stem cell (iPSC)-derived organoids provide models to study human organ development. Single-cell transcriptomics enable highly resolved descriptions of cell states within these systems; however, approaches are needed to directly measure lineage relationships. Here we establish iTracer, a lineage recorder that combines reporter barcodes with inducible CRISPR-Cas9 scarring and is compatible with single-cell and spatial transcriptomics. We apply iTracer to explore clonality and lineage dynamics during cerebral organoid development and identify a time window of fate restriction as well as variation in neurogenic dynamics between progenitor neuron families. We also establish long-term four-dimensional light-sheet microscopy for spatial lineage recording in cerebral organoids and confirm regional clonality in the developing neuroepithelium. We incorporate gene perturbation (iTracer-perturb) and assess the effect of mosaic TSC2 mutations on cerebral organoid development. Our data shed light on how lineages and fates are established during cerebral organoid formation. More broadly, our techniques can be adapted in any iPSC-derived culture system to dissect lineage alterations during normal or perturbed development.

Organoid single-cell genomic atlas uncovers human-specific features of brain development.

The human brain has undergone substantial change since humans diverged from chimpanzees and the other great apes1,2. However, the genetic and developmental programs that underlie this divergence are not fully understood. Here we have analysed stem cell-derived cerebral organoids using single-cell transcriptomics and accessible chromatin profiling to investigate gene-regulatory changes that are specific to humans. We first analysed cell composition and reconstructed differentiation trajectories over the entire course of human cerebral organoid development from pluripotency, through neuroectoderm and neuroepithelial stages, followed by divergence into neuronal fates within the dorsal and ventral forebrain, midbrain and hindbrain regions. Brain-region composition varied in organoids from different iPSC lines, but regional gene-expression patterns remained largely reproducible across individuals. We analysed chimpanzee and macaque cerebral organoids and found that human neuronal development occurs at a slower pace relative to the other two primates. Using pseudotemporal alignment of differentiation paths, we found that human-specific gene expression resolved to distinct cell states along progenitor-to-neuron lineages in the cortex. Chromatin accessibility was dynamic during cortex development, and we identified divergence in accessibility between human and chimpanzee that correlated with human-specific gene expression and genetic change. Finally, we mapped human-specific expression in adult prefrontal cortex using single-nucleus RNA sequencing analysis and identified developmental differences that persist into adulthood, as well as cell-state-specific changes that occur exclusively in the adult brain. Our data provide a temporal cell atlas of great ape forebrain development, and illuminate dynamic gene-regulatory features that are unique to humans.

Single-cell-resolution transcriptome map of human, chimpanzee, bonobo, and macaque brains.

Identification of gene expression traits unique to the human brain sheds light on the molecular mechanisms underlying human evolution. Here, we searched for uniquely human gene expression traits by analyzing 422 brain samples from humans, chimpanzees, bonobos, and macaques representing 33 anatomical regions, as well as 88,047 cell nuclei composing three of these regions. Among 33 regions, cerebral cortex areas, hypothalamus, and cerebellar gray and white matter evolved rapidly in humans. At the cellular level, astrocytes and oligodendrocyte progenitors displayed more differences in the human evolutionary lineage than the neurons. Comparison of the bulk tissue and single-nuclei sequencing revealed that conventional RNA sequencing did not detect up to two-thirds of cell-type-specific evolutionary differences.

Multilineage communication regulates human liver bud development from pluripotency.

Conventional two-dimensional differentiation from pluripotency fails to recapitulate cell interactions occurring during organogenesis. Three-dimensional organoids generate complex organ-like tissues; however, it is unclear how heterotypic interactions affect lineage identity. Here we use single-cell RNA sequencing to reconstruct hepatocyte-like lineage progression from pluripotency in two-dimensional culture. We then derive three-dimensional liver bud organoids by reconstituting hepatic, stromal, and endothelial interactions, and deconstruct heterogeneity during liver bud development. We find that liver bud hepatoblasts diverge from the two-dimensional lineage, and express epithelial migration signatures characteristic of organ budding. We benchmark three-dimensional liver buds against fetal and adult human liver single-cell RNA sequencing data, and find a striking correspondence between the three-dimensional liver bud and fetal liver cells. We use a receptor-ligand pairing analysis and a high-throughput inhibitor assay to interrogate signalling in liver buds, and show that vascular endothelial growth factor (VEGF) crosstalk potentiates endothelial network formation and hepatoblast differentiation. Our molecular dissection reveals interlineage communication regulating organoid development, and illuminates previously inaccessible aspects of human liver development.

Dissecting direct reprogramming from fibroblast to neuron using single-cell RNA-seq.

Direct lineage reprogramming represents a remarkable conversion of cellular and transcriptome states. However, the intermediate stages through which individual cells progress during reprogramming are largely undefined. Here we use single-cell RNA sequencing at multiple time points to dissect direct reprogramming from mouse embryonic fibroblasts to induced neuronal cells. By deconstructing heterogeneity at each time point and ordering cells by transcriptome similarity, we find that the molecular reprogramming path is remarkably continuous. Overexpression of the proneural pioneer factor Ascl1 results in a well-defined initialization, causing cells to exit the cell cycle and re-focus gene expression through distinct neural transcription factors. The initial transcriptional response is relatively homogeneous among fibroblasts, suggesting that the early steps are not limiting for productive reprogramming. Instead, the later emergence of a competing myogenic program and variable transgene dynamics over time appear to be the major efficiency limits of direct reprogramming. Moreover, a transcriptional state, distinct from donor and target cell programs, is transiently induced in cells undergoing productive reprogramming. Our data provide a high-resolution approach for understanding transcriptome states during lineage differentiation.

A novel population of Hopx-dependent basal radial glial cells in the developing mouse neocortex.

A specific subpopulation of neural progenitor cells, the basal radial glial cells (bRGCs) of the outer subventricular zone (OSVZ), are thought to have a key role in the evolutionary expansion of the mammalian neocortex. In the developing lissencephalic mouse neocortex, bRGCs exist at low abundance and show significant molecular differences from bRGCs in developing gyrencephalic species. Here, we demonstrate that the developing mouse medial neocortex (medNcx), in contrast to the canonically studied lateral neocortex (latNcx), exhibits an OSVZ and an abundance of bRGCs similar to that in developing gyrencephalic neocortex. Unlike bRGCs in developing mouse latNcx, the bRGCs in medNcx exhibit human bRGC-like gene expression, including expression of Hopx, a human bRGC marker. Disruption of Hopx expression in mouse embryonic medNcx and forced Hopx expression in mouse embryonic latNcx demonstrate that Hopx is required and sufficient, respectively, for bRGC abundance as found in the developing gyrencephalic neocortex. Taken together, our data identify a novel bRGC subpopulation in developing mouse medNcx that is highly related to bRGCs of developing gyrencephalic neocortex.

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.

Lineage-negative progenitors mobilize to regenerate lung epithelium after major injury.

Broadly, tissue regeneration is achieved in two ways: by proliferation of common differentiated cells and/or by deployment of specialized stem/progenitor cells. Which of these pathways applies is both organ- and injury-specific. Current models in the lung posit that epithelial repair can be attributed to cells expressing mature lineage markers. By contrast, here we define the regenerative role of previously uncharacterized, rare lineage-negative epithelial stem/progenitor (LNEP) cells present within normal distal lung. Quiescent LNEPs activate a ΔNp63 (a p63 splice variant) and cytokeratin 5 remodelling program after influenza or bleomycin injury in mice. Activated cells proliferate and migrate widely to occupy heavily injured areas depleted of mature lineages, at which point they differentiate towards mature epithelium. Lineage tracing revealed scant contribution of pre-existing mature epithelial cells in such repair, whereas orthotopic transplantation of LNEPs, isolated by a definitive surface profile identified through single-cell sequencing, directly demonstrated the proliferative capacity and multipotency of this population. LNEPs require Notch signalling to activate the ΔNp63 and cytokeratin 5 program, and subsequent Notch blockade promotes an alveolar cell fate. Persistent Notch signalling after injury led to parenchymal 'micro-honeycombing' (alveolar cysts), indicative of failed regeneration. Lungs from patients with fibrosis show analogous honeycomb cysts with evidence of hyperactive Notch signalling. Our findings indicate that distinct stem/progenitor cell pools repopulate injured tissue depending on the extent of the injury, and the outcomes of regeneration or fibrosis may depend in part on the dynamics of LNEP Notch signalling.

Reconstructing complex lineage trees from scRNA-seq data using MERLoT.

Advances in single-cell transcriptomics techniques are revolutionizing studies of cellular differentiation and heterogeneity. It has become possible to track the trajectory of thousands of genes across the cellular lineage trees that represent the temporal emergence of cell types during dynamic processes. However, reconstruction of cellular lineage trees with more than a few cell fates has proved challenging. We present MERLoT (https://github.com/soedinglab/merlot), a flexible and user-friendly tool to reconstruct complex lineage trees from single-cell transcriptomics data. It can impute temporal gene expression profiles along the reconstructed tree. We show MERLoT's capabilities on various real cases and hundreds of simulated datasets.

Human organomics: a fresh approach to understanding human development using single-cell transcriptomics.

Innovative methods designed to recapitulate human organogenesis from pluripotent stem cells provide a means to explore human developmental biology. New technologies to sequence and analyze single-cell transcriptomes can deconstruct these 'organoids' into constituent parts, and reconstruct lineage trajectories during cell differentiation. In this Spotlight article we summarize the different approaches to performing single-cell transcriptomics on organoids, and discuss the opportunities and challenges of applying these techniques to generate organ-level, mechanistic models of human development and disease. Together, these technologies will move past characterization to the prediction of human developmental and disease-related phenomena.

Reconstructing lineage hierarchies of the distal lung epithelium using single-cell RNA-seq.

The mammalian lung is a highly branched network in which the distal regions of the bronchial tree transform during development into a densely packed honeycomb of alveolar air sacs that mediate gas exchange. Although this transformation has been studied by marker expression analysis and fate-mapping, the mechanisms that control the progression of lung progenitors along distinct lineages into mature alveolar cell types are still incompletely known, in part because of the limited number of lineage markers and the effects of ensemble averaging in conventional transcriptome analysis experiments on cell populations. Here we show that single-cell transcriptome analysis circumvents these problems and enables direct measurement of the various cell types and hierarchies in the developing lung. We used microfluidic single-cell RNA sequencing (RNA-seq) on 198 individual cells at four different stages encompassing alveolar differentiation to measure the transcriptional states which define the developmental and cellular hierarchy of the distal mouse lung epithelium. We empirically classified cells into distinct groups by using an unbiased genome-wide approach that did not require a priori knowledge of the underlying cell types or the previous purification of cell populations. The results confirmed the basic outlines of the classical model of epithelial cell-type diversity in the distal lung and led to the discovery of many previously unknown cell-type markers, including transcriptional regulators that discriminate between the different populations. We reconstructed the molecular steps during maturation of bipotential progenitors along both alveolar lineages and elucidated the full life cycle of the alveolar type 2 cell lineage. This single-cell genomics approach is applicable to any developing or mature tissue to robustly delineate molecularly distinct cell types, define progenitors and lineage hierarchies, and identify lineage-specific regulatory factors.

Dynamic architecture of a minimal RNA polymerase II open promoter complex.

The open promoter complex (OC) is a central intermediate during transcription initiation that contains a DNA bubble. Here, we employ single-molecule Förster resonance energy transfer experiments and Nano-Positioning System analysis to determine the three-dimensional architecture of a minimal OC consisting of promoter DNA, including a TATA box and an 11-nucleotide mismatched region around the transcription start site, TATA box-binding protein (TBP), RNA polymerase (Pol) II, and general transcription factor (TF)IIB and TFIIF. In this minimal OC, TATA-DNA and TBP reside above the Pol II cleft between clamp and protrusion domains. Downstream DNA is dynamically loaded into and unloaded from the Pol II cleft at a timescale of seconds. The TFIIB core domain is displaced from the Pol II wall, where it is located in the closed promoter complex. These results reveal large overall structural changes during the initiation-elongation transition, which are apparently accommodated by the intrinsic flexibility of TFIIB.