Engineering fibronectin-templated multi-component fibrillar extracellular matrices to modulate tissue-specific cell response

Cells assemble fibronectin, the major extracellular matrix (ECM) protein, into fibrillar matrices, which serve as 3D architectural scaffolds to provide, together with other ECM proteins tissue-specific environments. Although recent approaches enable to bioengineer 3D fibrillar fibronectin matrices in vitro, it remains elusive how fibronectin can be co-assembled with other ECM proteins into complex 3D fibrillar matrices that recapitulate tissue-specific compositions and cellular responses. Here, we introduce the engineering of fibrillar fibronectin-templated 3D matrices that can be complemented with other ECM proteins, including vitronectin, collagen, and laminin to resemble ECM architectures observed in vivo. For the co-assembly of different ECM proteins, we employed their innate fibrillogenic mechanisms including shear forces, pH-dependent electrostatic interactions, or specific binding domains. Through recapitulating various tissue-specific ECM compositions and morphologies, the large scale multi-composite 3D fibrillar ECM matrices can guide fibroblast adhesion, 3D fibroblast tissue formation, or tissue morphogenesis of epithelial cells. In other examples, we customize multi-composite 3D fibrillar matrices to support the growth of signal propagating neuronal networks and of human brain organoids. We envision that these 3D fibrillar ECM matrices can be tailored in scale and composition to modulate tissue-specific responses across various biological length scales and systems, and thus to advance manyfold studies of cell biological systems.

Multimodal spatiotemporal phenotyping of human retinal organoid development

Organoids generated from human pluripotent stem cells provide experimental systems to study development and disease, but quantitative measurements across different spatial scales and molecular modalities are lacking. In this study, we generated multiplexed protein maps over a retinal organoid time course and primary adult human retinal tissue. We developed a toolkit to visualize progenitor and neuron location, the spatial arrangements of extracellular and subcellular components and global patterning in each organoid and primary tissue. In addition, we generated a single-cell transcriptome and chromatin accessibility timecourse dataset and inferred a gene regulatory network underlying organoid development. We integrated genomic data with spatially segmented nuclei into a multimodal atlas to explore organoid patterning and retinal ganglion cell (RGC) spatial neighborhoods, highlighting pathways involved in RGC cell death and showing that mosaic genetic perturbations in retinal organoids provide insight into cell fate regulation.

Metformin rescues migratory deficits of cells derived from patients with periventricular heterotopia

Periventricular neuronal heterotopia (PH) is one of the most common forms of cortical malformation in the human cortex. We show that human neuronal progenitor cells (hNPCs) derived from PH patients with a DCHS1 or FAT4 mutation as well as isogenic lines had altered migratory dynamics when grafted in the mouse brain. The affected migration was linked to altered autophagy as observed in vivo with an electron microscopic analysis of grafted hNPCs, a Western blot analysis of cortical organoids, and time-lapse imaging of hNPCs in the presence of bafilomycin A1. We further show that deficits in autophagy resulted in the accumulation of paxillin, a focal adhesion protein involved in cell migration. Strikingly, a single-cell RNA-seq analysis of hNPCs revealed similar expression levels of autophagy-related genes. Bolstering AMPK-dependent autophagy by metformin, an FDA-approved drug, promoted migration of PH patients-derived hNPCs. Our data indicate that transcription-independent homeostatic modifications in autophagy contributed to the defective migratory behavior of hNPCs in vivo and suggest that modulating autophagy in hNPCs might rescue neuronal migration deficits in some forms of PH.

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.

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.

What is a cell type?

A next step for cell atlases should be to chart perturbations in human model systems.

Phospho-seq: Integrated, multi-modal profiling of intracellular protein dynamics in single cells

Cell signaling plays a critical role in regulating cellular behavior and fate. While multimodal single-cell sequencing technologies are rapidly advancing, scalable and flexible profiling of cell signaling states alongside other molecular modalities remains challenging. Here we present Phospho-seq, an integrated approach that aims to quantify phosphorylated intracellular and intranuclear proteins, and to connect their activity with cis-regulatory elements and transcriptional targets. We utilize a simplified benchtop antibody conjugation method to create large custom antibody panels for simultaneous protein and scATAC-seq profiling on whole cells, and integrate this information with scRNA-seq datasets via bridge integration. We apply our workflow to cell lines, induced pluripotent stem cells, and 3-month-old brain organoids to demonstrate its broad applicability. We demonstrate that Phospho-seq can define cellular states and trajectories, reconstruct gene regulatory relationships, and characterize the causes and consequences of heterogeneous cell signaling in neurodevelopment.

Single-cell analyses of axolotl telencephalon organization, neurogenesis, and regeneration

Salamanders are tetrapod models to study brain organization and regeneration; however, the identity and evolutionary conservation of brain cell types are largely unknown. We delineated the cell populations in the axolotl telencephalon during homeostasis and regeneration using single-cell genomic profiling. We identified glutamatergic neurons with similarities to amniote neurons of hippocampus, dorsal and lateral cortex, and conserved γ-aminobutyric acid-releasing (GABAergic) neuron classes. We inferred transcriptional dynamics and gene regulatory relationships of postembryonic, region-specific neurogenesis and unraveled conserved differentiation signatures. After brain injury, ependymoglia activate an injury-specific state before reestablishing lost neuron populations and axonal connections. Together, our analyses yield insights into the organization, evolution, and regeneration of a tetrapod nervous system.

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.

High-throughput muscle fiber typing from RNA sequencing data

Background

Skeletal muscle fiber type distribution has implications for human health, muscle function, and performance. This knowledge has been gathered using labor-intensive and costly methodology that limited these studies. Here, we present a method based on muscle tissue RNA sequencing data (totRNAseq) to estimate the distribution of skeletal muscle fiber types from frozen human samples, allowing for a larger number of individuals to be tested.

Methods

By using single-nuclei RNA sequencing (snRNAseq) data as a reference, cluster expression signatures were produced by averaging gene expression of cluster gene markers and then applying these to totRNAseq data and inferring muscle fiber nuclei type via linear matrix decomposition. This estimate was then compared with fiber type distribution measured by ATPase staining or myosin heavy chain protein isoform distribution of 62 muscle samples in two independent cohorts (n = 39 and 22).

Results

The correlation between the sequencing-based method and the other two were r_ATPas = 0.44 [0.13–0.67], [95% CI], and r_myosin = 0.83 [0.61–0.93], with p = 5.70 × 10^–3 and 2.00 × 10^–6, respectively. The deconvolution inference of fiber type composition was accurate even for very low totRNAseq sequencing depths, i.e., down to an average of ~ 10,000 paired-end reads.

Conclusions

This new method (https://github.com/OlaHanssonLab/PredictFiberType) consequently allows for measurement of fiber type distribution of a larger number of samples using totRNAseq in a cost and labor-efficient way. It is now feasible to study the association between fiber type distribution and e.g. health outcomes in large well-powered studies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13395-022-00299-4.

High temporal resolution proteome and phosphoproteome profiling of stem cell-derived hepatocyte development

Primary human hepatocytes are widely used to evaluate liver toxicity of drugs, but they are scarce and demanding to culture. Stem cell-derived hepatocytes are increasingly discussed as alternatives. To obtain a better appreciation of the molecular processes during the differentiation of induced pluripotent stem cells into hepatocytes, we employ a quantitative proteomic approach to follow the expression of 9,000 proteins, 12,000 phosphorylation sites, and 800 acetylation sites over time. The analysis reveals stage-specific markers, a major molecular switch between hepatic endoderm versus immature hepatocyte-like cells impacting, e.g., metabolism, the cell cycle, kinase activity, and the expression of drug transporters. Comparing the proteomes of two- (2D) and three-dimensional (3D)-derived hepatocytes with fetal and adult liver indicates a fetal-like status of the in vitro models and lower expression of important ADME/Tox proteins. The collective data enable constructing a molecular roadmap of hepatocyte development that serves as a valuable resource for future research.

Spatial transcriptomic and single-nucleus analysis reveals heterogeneity in a gigantic single-celled syncytium

In multicellular organisms, the specification, coordination, and compartmentalization of cell types enable the formation of complex body plans. However, some eukaryotic protists such as slime molds generate diverse and complex structures while remaining in a multinucleate syncytial state. It is unknown if different regions of these giant syncytial cells have distinct transcriptional responses to environmental encounters and if nuclei within the cell diversify into heterogeneous states. Here, we performed spatial transcriptome analysis of the slime mold Physarum polycephalum in the plasmodium state under different environmental conditions and used single-nucleus RNA-sequencing to dissect gene expression heterogeneity among nuclei. Our data identifies transcriptome regionality in the organism that associates with proliferation, syncytial substructures, and localized environmental conditions. Further, we find that nuclei are heterogenous in their transcriptional profile and may process local signals within the plasmodium to coordinate cell growth, metabolism, and reproduction. To understand how nuclei variation within the syncytium compares to heterogeneity in single-nucleus cells, we analyzed states in single Physarum amoebal cells. We observed amoebal cell states at different stages of mitosis and meiosis, and identified cytokinetic features that are specific to nuclei divisions within the syncytium. Notably, we do not find evidence for predefined transcriptomic states in the amoebae that are observed in the syncytium. Our data shows that a single-celled slime mold can control its gene expression in a region-specific manner while lacking cellular compartmentalization and suggests that nuclei are mobile processors facilitating local specialized functions. More broadly, slime molds offer the extraordinary opportunity to explore how organisms can evolve regulatory mechanisms to divide labor, specialize, balance competition with cooperation, and perform other foundational principles that govern the logic of life.

Extracellular LGALS3BP regulates neural progenitor position and relates to human cortical complexity

Basal progenitors (BPs), including intermediate progenitors and basal radial glia, are generated from apical radial glia and are enriched in gyrencephalic species like humans, contributing to neuronal expansion. Shortly after generation, BPs delaminate towards the subventricular zone, where they further proliferate before differentiation. Gene expression alterations involved in BP delamination and function in humans are poorly understood. Here, we study the role of LGALS3BP, so far known as a cancer biomarker, which is a secreted protein enriched in human neural progenitors (NPCs). We show that individuals with LGALS3BP de novo variants exhibit altered local gyrification, sulcal depth, surface area and thickness in their cortex. Additionally, using cerebral organoids, human fetal tissues and mice, we show that LGALS3BP regulates the position of NPCs. Single-cell RNA-sequencing and proteomics reveal that LGALS3BP-mediated mechanisms involve the extracellular matrix in NPCs' anchoring and migration within the human brain. We propose that its temporal expression influences NPCs' delamination, corticogenesis and gyrification extrinsically.

Engineering of fully humanized and vascularized 3D bone marrow niches sustaining undifferentiated human cord blood hematopoietic stem and progenitor cells

Hematopoietic stem and progenitor cells (HSPCs) are frequently located around the bone marrow (BM) vasculature. These so-called perivascular niches regulate HSC function both in health and disease, but they have been poorly studied in humans due to the scarcity of models integrating complete human vascular structures. Herein, we propose the stromal vascular fraction (SVF) derived from human adipose tissue as a cell source to vascularize 3D osteoblastic BM niches engineered in perfusion bioreactors. We show that SVF cells form self-assembled capillary structures, composed by endothelial and perivascular cells, that add to the osteogenic matrix secreted by BM mesenchymal stromal cells in these engineered niches. In comparison to avascular osteoblastic niches, vascularized BM niches better maintain immunophenotypically-defined cord blood (CB) HSCs without affecting cell proliferation. In contrast, HSPCs cultured in vascularized BM niches showed increased CFU-granulocyte-erythrocyte-monocyte-megakaryocyte (CFU-GEMM) numbers. The vascularization also contributed to better preserve osteogenic gene expression in the niche, demonstrating that niche vascularization has an influence on both hematopoietic and stromal compartments. In summary, we have engineered a fully humanized and vascularized 3D BM tissue to model native human endosteal perivascular niches and revealed functional implications of this vascularization in sustaining undifferentiated CB HSPCs. This system provides a unique modular platform to explore hemato-vascular interactions in human healthy/pathological hematopoiesis.

NGN2 induces diverse neuron types from human pluripotency

Human neurons engineered from induced pluripotent stem cells (iPSCs) through neurogenin 2 (NGN2) overexpression are widely used to study neuronal differentiation mechanisms and to model neurological diseases. However, the differentiation paths and heterogeneity of emerged neurons have not been fully explored. Here, we used single-cell transcriptomics to dissect the cell states that emerge during NGN2 overexpression across a time course from pluripotency to neuron functional maturation. We find a substantial molecular heterogeneity in the neuron types generated, with at least two populations that express genes associated with neurons of the peripheral nervous system. Neuron heterogeneity is observed across multiple iPSC clones and lines from different individuals. We find that neuron fate acquisition is sensitive to NGN2 expression level and the duration of NGN2-forced expression. Our data reveal that NGN2 dosage can regulate neuron fate acquisition, and that NGN2-iN heterogeneity can confound results that are sensitive to neuron type.

•

NGN2-iNs are molecularly heterogeneous

•

NGN2-iNs subtypes have signatures of central and peripheral nervous system

•

Neural fate acquisition is sensitive to the level and duration of NGN2 expression

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.

Mapping Development of the Human Intestinal Niche at Single-Cell Resolution

The human intestinal stem cell niche supports self-renewal and epithelial function, but little is known about its development. We used single-cell mRNA sequencing with in situ validation approaches to interrogate human intestinal development from 7-21 weeks post conception, assigning molecular identities and spatial locations to cells and factors that comprise the niche. Smooth muscle cells of the muscularis mucosa, in close proximity to proliferative crypts, are a source of WNT and RSPONDIN ligands, whereas EGF is expressed far from crypts in the villus epithelium. Instead, an PDGFRAHI/F3HI/DLL1HI mesenchymal population lines the crypt-villus axis and is the source of the epidermal growth factor (EGF) family member NEUREGULIN1 (NRG1). In developing intestine enteroid cultures, NRG1, but not EGF, permitted increased cellular diversity via differentiation of secretory lineages. This work highlights the complexities of intestinal EGF/ERBB signaling and delineates key niche cells and signals of the developing intestine.

Comparison of induced neurons reveals slower structural and functional maturation in humans than in apes

We generated induced excitatory neurons (iNeurons, iNs) from chimpanzee, bonobo, and human stem cells by expressing the transcription factor neurogenin-2 (NGN2). Single-cell RNA sequencing showed that genes involved in dendrite and synapse development are expressed earlier during iNs maturation in the chimpanzee and bonobo than the human cells. In accordance, during the first 2 weeks of differentiation, chimpanzee and bonobo iNs showed repetitive action potentials and more spontaneous excitatory activity than human iNs, and extended neurites of higher total length. However, the axons of human iNs were slightly longer at 5 weeks of differentiation. The timing of the establishment of neuronal polarity did not differ between the species. Chimpanzee, bonobo, and human neurites eventually reached the same level of structural complexity. Thus, human iNs develop slower than chimpanzee and bonobo iNs, and this difference in timing likely depends on functions downstream of NGN2.

CSS: cluster similarity spectrum integration of single-cell genomics data

It is a major challenge to integrate single-cell sequencing data across experiments, conditions, batches, time points, and other technical considerations. New computational methods are required that can integrate samples while simultaneously preserving biological information. Here, we propose an unsupervised reference-free data representation, cluster similarity spectrum (CSS), where each cell is represented by its similarities to clusters independently identified across samples. We show that CSS can be used to assess cellular heterogeneity and enable reconstruction of differentiation trajectories from cerebral organoid and other single-cell transcriptomic data, and to integrate data across experimental conditions and human individuals.

Human Stem Cell Resources Are an Inroad to Neandertal DNA Functions

Induced pluripotent stem cells (iPSCs) from diverse humans offer the potential to study human functional variation in controlled culture environments. A portion of this variation originates from an ancient admixture between modern humans and Neandertals, which introduced alleles that left a phenotypic legacy on individual humans today. Here, we show that a large iPSC repository harbors extensive Neandertal DNA, including alleles that contribute to human phenotypes and diseases, encode hundreds of amino acid changes, and alter gene expression in specific tissues. We provide a database of the inferred introgressed Neandertal alleles for each individual iPSC line, together with the annotation of the predicted functional variants. We also show that transcriptomic data from organoids generated from iPSCs can be used to track Neandertal-derived RNA over developmental processes. Human iPSC resources provide an opportunity to experimentally explore Neandertal DNA function and its contribution to present-day phenotypes, and potentially study Neandertal traits.

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.

Single-cell genomic analysis of human cerebral organoids

Investigating early brain development has previously relied on using primary developing brain tissue or two-dimensional cell culture models. Recently, stem cell-derived three-dimensional cell culture systems, collectively called brain organoids, have been developed that can faithfully recapitulate many aspects of early brain development. Together with the ability to reprogram fibroblast or blood cells into induced pluripotent stem cells from humans with neurodevelopmental disorders, this opens new inroads to study patient-specific brain development in a personalized cell culture model. Studying the transcriptomes and regulatory landscape of single cells within brain organoids presents a major advance to understand cell-type specific features and transient states during development, and to link these states to their underlying regulatory logic at high resolution. In this protocol, we describe how to generate single-cell RNA-seq and ATAC-seq (Assay for Transposase-Accessible Chromatin using sequencing) data from the same suspension of organoid cells and focus on reducing batch effects by multiplexing multiple individuals in one experiment. Moreover, we outline basic data processing, analysis, and strategies to correct for batch effects, to account for variability in organoids and for integrating gene expression and open chromatin data.

In Vitro and In Vivo Development of the Human Airway at Single-Cell Resolution

Bud tip progenitor cells give rise to all murine lung epithelial lineages and have been described in the developing human lung; however, the mechanisms controlling human bud tip differentiation into specific lineages are unclear. Here, we used homogeneous human bud tip organoid cultures and identified SMAD signaling as a key regulator of the bud tip-to-airway transition. SMAD induction led to the differentiation of airway-like organoids possessing functional basal cells capable of clonal expansion and multilineage differentiation. To benchmark in vitro-derived organoids, we developed a single-cell mRNA sequencing atlas of the human lung from 11.5 to 21 weeks of development, which revealed high degrees of similarity between the in vitro-derived and in vivo airway. Together, this work sheds light on human airway differentiation in vitro and provides a single-cell atlas of the developing human lung.

Mapping human cell phenotypes to genotypes with single-cell genomics

The cumulative activity of all of the body's cells, with their myriad interactions, life histories, and environmental experiences, gives rise to a condition that is distinctly human and specific to each individual. It is an enduring goal to catalog our human cell types, to understand how they develop, how they vary between individuals, and how they fail in disease. Single-cell genomics has revolutionized this endeavor because sequencing-based methods provide a means to quantitatively annotate cell states on the basis of high-information content and high-throughput measurements. Together with advances in stem cell biology and gene editing, we are in the midst of a fascinating journey to understand the cellular phenotypes that compose human bodies and how the human genome is used to build and maintain each cell. Here, we will review recent advances into how single-cell genomics is being used to develop personalized phenotyping strategies that cross subcellular, cellular, and tissue scales to link our genome to our cumulative cellular phenotypes.

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.

Single Cell Transcriptomes Derived from Human Cervical and Uterine Tissue during Pregnancy

This work presents the workflow for generating single cell transcriptomes derived from primary human uterine and cervical tissue obtained during planned cesarean hysterectomies. In total, a catalogue of 310 single cell transcriptomes are obtained, cell types present in these biopsies are inferred, and specific genes defining each of the cellular types present in the tissue are identified. Further validation of the inferred cell identity is also demonstrated via meta-analysis of independent repositories in literature generated by bulk sequenced data of fluorescence-activated cell sorting sorted cells.

Evaluation of Wu et al.: Strengthening Single-Cell Bioinformatic Comparisons for Improving Organoid Differentiation

An example of the peer review process, for "Comparative Analysis and Refinement of Human PSC-Derived Kidney Organoid Differentiation with Single-Cell Transcriptomics" (Wu et al., 2018), is presented herein.

High-throughput single-cell transcriptomics on organoids

Three-dimensional (3D) tissues grown in culture from human stem cells offer the incredible opportunity to analyze and manipulate human development, and to generate patient-specific models of disease. Methods to sequence DNA and RNA in single cells are being used to analyze these so-called 'organoid' systems in high-resolution. Single-cell transcriptomics has been used to quantitate the similarity of organoid cells to primary tissue counterparts in the brain, intestine, liver, and kidney, as well as identify cell-specific responses to environmental variables and disease conditions. The merging of these two technologies, single-cell genomics and organoids, will have profound impact on personalized medicine in the near future.

Single-cell analysis uncovers convergence of cell identities during axolotl limb regeneration

Amputation of the axolotl forelimb results in the formation of a blastema, a transient tissue where progenitor cells accumulate prior to limb regeneration. However, the molecular understanding of blastema formation had previously been hampered by the inability to identify and isolate blastema precursor cells in the adult tissue. We have used a combination of Cre-loxP reporter lineage tracking and single-cell messenger RNA sequencing (scRNA-seq) to molecularly track mature connective tissue (CT) cell heterogeneity and its transition to a limb blastema state. We have uncovered a multiphasic molecular program where CT cell types found in the uninjured adult limb revert to a relatively homogenous progenitor state that recapitulates an embryonic limb bud-like phenotype including multipotency within the CT lineage. Together, our data illuminate molecular and cellular reprogramming during complex organ regeneration in a vertebrate.

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.

Optimal Hypoxia Regulates Human iPSC-Derived Liver Bud Differentiation through Intercellular TGFB Signaling

Timely controlled oxygen (O2) delivery is crucial for the developing liver. However, the influence of O2 on intercellular communication during hepatogenesis is unclear. Using a human induced pluripotent stem cell-derived liver bud (hiPSC-LB) model, we found hypoxia induced with an O2-permeable plate promoted hepatic differentiation accompanied by TGFB1 and TGFB3 suppression. Conversely, extensive hypoxia generated with an O2-non-permeable plate elevated TGFBs and cholangiocyte marker expression. Single-cell RNA sequencing revealed that TGFB1 and TGFB3 are primarily expressed in the human liver mesenchyme and endothelium similar to in the hiPSC-LBs. Stromal cell-specific RNA interferences indicated the importance of TGFB signaling for hepatocytic differentiation in hiPSC-LB. Consistently, during mouse liver development, the Hif1a-mediated developmental hypoxic response is positively correlated with TGFB1 expression. These data provide insights into the mechanism that hypoxia-stimulated signals in mesenchyme and endothelium, likely through TGFB1, promote hepatoblast differentiation prior to fetal circulation establishment.

ShinyCortex: Exploring Single-Cell Transcriptome Data From the Developing Human Cortex

Single-cell mRNA sequencing (scRNA-seq) is a powerful method to identify and classify cell types and reconstruct differentiation trajectories within complex tissues, such as the developing human cortex. scRNA-seq data also enables the discovery of cell type-specific marker genes and genes that regulate developmental transitions. Here we provide a brief overview of how scRNA-seq has been shaping the study of human cortex development, and present ShinyCortex, a resource that brings together data from recent scRNA-seq studies of the developing cortex for further analysis. ShinyCortex is based in R and displays recently published scRNA-seq data from the human and mouse cortex in a comprehensible, dynamic and accessible way, suitable for data exploration by biologists.

Mapping heterogeneity in patient-derived melanoma cultures by single-cell RNA-seq

Recent technological advances in single-cell genomics make it possible to analyze cellular heterogeneity of tumor samples. Here, we applied single-cell RNA-seq to measure the transcriptomes of 307 single cells cultured from three biopsies of three different patients with a BRAF/NRAS wild type, BRAF mutant/NRAS wild type and BRAF wild type/NRAS mutant melanoma metastasis, respectively. Analysis based on self-organizing maps identified sub-populations defined by multiple gene expression modules involved in proliferation, oxidative phosphorylation, pigmentation and cellular stroma. Gene expression modules had prognostic relevance when compared with gene expression data from published melanoma samples and patient survival data. We surveyed kinome expression patterns across sub-populations of the BRAF/NRAS wild type sample and found that CDK4 and CDK2 were consistently highly expressed in the majority of cells, suggesting that these kinases might be involved in melanoma progression. Treatment of cells with the CDK4 inhibitor palbociclib restricted cell proliferation to a similar, and in some cases greater, extent than MAPK inhibitors. Finally, we identified a low abundant sub-population in this sample that highly expressed a module containing ABC transporter ABCB5, surface markers CD271 and CD133, and multiple aldehyde dehydrogenases (ALDHs). Patient-derived cultures of the BRAF mutant/NRAS wild type and BRAF wild type/NRAS mutant metastases showed more homogeneous single-cell gene expression patterns with gene expression modules for proliferation and ABC transporters. Taken together, our results describe an intertumor and intratumor heterogeneity in melanoma short-term cultures which might be relevant for patient survival, and suggest promising targets for new treatment approaches in melanoma therapy.

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.

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.

Structural reorganization of the chromatin remodeling enzyme Chd1 upon engagement with nucleosomes

The yeast Chd1 protein acts to position nucleosomes across genomes. Here, we model the structure of the Chd1 protein in solution and when bound to nucleosomes. In the apo state, the DNA-binding domain contacts the edge of the nucleosome while in the presence of the non-hydrolyzable ATP analog, ADP-beryllium fluoride, we observe additional interactions between the ATPase domain and the adjacent DNA gyre 1.5 helical turns from the dyad axis of symmetry. Binding in this conformation involves unravelling the outer turn of nucleosomal DNA and requires substantial reorientation of the DNA-binding domain with respect to the ATPase domains. The orientation of the DNA-binding domain is mediated by sequences in the N-terminus and mutations to this part of the protein have positive and negative effects on Chd1 activity. These observations indicate that the unfavorable alignment of C-terminal DNA-binding region in solution contributes to an auto-inhibited state.

DOI:http://dx.doi.org/10.7554/eLife.22510.001

Differences and similarities between human and chimpanzee neural progenitors during cerebral cortex development

Human neocortex expansion likely contributed to the remarkable cognitive abilities of humans. This expansion is thought to primarily reflect differences in proliferation versus differentiation of neural progenitors during cortical development. Here, we have searched for such differences by analysing cerebral organoids from human and chimpanzees using immunohistofluorescence, live imaging, and single-cell transcriptomics. We find that the cytoarchitecture, cell type composition, and neurogenic gene expression programs of humans and chimpanzees are remarkably similar. Notably, however, live imaging of apical progenitor mitosis uncovered a lengthening of prometaphase-metaphase in humans compared to chimpanzees that is specific to proliferating progenitors and not observed in non-neural cells. Consistent with this, the small set of genes more highly expressed in human apical progenitors points to increased proliferative capacity, and the proportion of neurogenic basal progenitors is lower in humans. These subtle differences in cortical progenitors between humans and chimpanzees may have consequences for human neocortex evolution.

Cellular Taxonomy of the Mouse Striatum as Revealed by Single-Cell RNA-Seq

The striatum contributes to many cognitive processes and disorders, but its cell types are incompletely characterized. We show that microfluidic and FACS-based single-cell RNA sequencing of mouse striatum provides a well-resolved classification of striatal cell type diversity. Transcriptome analysis revealed ten differentiated, distinct cell types, including neurons, astrocytes, oligodendrocytes, ependymal, immune, and vascular cells, and enabled the discovery of numerous marker genes. Furthermore, we identified two discrete subtypes of medium spiny neurons (MSNs) that have specific markers and that overexpress genes linked to cognitive disorders and addiction. We also describe continuous cellular identities, which increase heterogeneity within discrete cell types. Finally, we identified cell type-specific transcription and splicing factors that shape cellular identities by regulating splicing and expression patterns. Our findings suggest that functional diversity within a complex tissue arises from a small number of discrete cell types, which can exist in a continuous spectrum of functional states.

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.