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.

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.

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.

What is a cell type?

A next step for cell atlases should be to chart perturbations in human model 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.

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.

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 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 rATPas = 0.44 [0.13-0.67], [95% CI], and rmyosin = 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.

Characterization of RNA content in individual phase-separated coacervate microdroplets.

Condensates formed by complex coacervation are hypothesized to have played a crucial part during the origin-of-life. In living cells, condensation organizes biomolecules into a wide range of membraneless compartments. Although RNA is a key component of biological condensates and the central component of the RNA world hypothesis, little is known about what determines RNA accumulation in condensates and to which extend single condensates differ in their RNA composition. To address this, we developed an approach to read the RNA content from single synthetic and protein-based condensates using high-throughput sequencing. We find that certain RNAs efficiently accumulate in condensates. These RNAs are strongly enriched in sequence motifs which show high sequence similarity to short interspersed elements (SINEs). We observe similar results for protein-derived condensates, demonstrating applicability across different in vitro reconstituted membraneless organelles. Thus, our results provide a new inroad to explore the RNA content of phase-separated droplets at single condensate resolution.

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.

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.

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.

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.