Organogenetic transcriptomes of the Drosophila embryo at single cell resolution.

To gain insight into the transcription programs activated during the formation of Drosophila larval structures, we carried out single cell RNA sequencing during two periods of Drosophila embryogenesis: stages 10-12, when most organs are first specified and initiate morphological and physiological specialization; and stages 13-16, when organs achieve their final mature architectures and begin to function. Our data confirm previous findings with regards to functional specialization of some organs - the salivary gland and trachea - and clarify the embryonic functions of another - the plasmatocytes. We also identify two early developmental trajectories in germ cells and uncover a potential role for proteolysis during germline stem cell specialization. We identify the likely cell type of origin for key components of the Drosophila matrisome and several commonly used Drosophila embryonic cell culture lines. Finally, we compare our findings with other recent related studies and with other modalities for identifying tissue-specific gene expression patterns. These data provide a useful community resource for identifying many new players in tissue-specific morphogenesis and functional specialization of developing organs.

Generation and molecular characterization of human pluripotent stem cell-derived pharyngeal foregut endoderm.

Approaches to study human pharyngeal foregut endoderm-a developmental intermediate that is linked to various human syndromes involving pharynx development and organogenesis of tissues such as thymus, parathyroid, and thyroid-have been hampered by scarcity of tissue access and cellular models. We present an efficient stepwise differentiation method to generate human pharyngeal foregut endoderm from pluripotent stem cells. We determine dose and temporal requirements of signaling pathway engagement for optimized differentiation and characterize the differentiation products on cellular and integrated molecular level. We present a computational classification tool, "CellMatch," and transcriptomic classification of differentiation products on an integrated mouse scRNA-seq developmental roadmap confirms cellular maturation. Integrated transcriptomic and chromatin analyses infer differentiation stage-specific gene regulatory networks. Our work provides the method and integrated multiomic resource for the investigation of disease-relevant loci and gene regulatory networks and their role in developmental defects affecting the pharyngeal endoderm and its derivatives.

A systematic comparison of computational methods for expression forecasting.

Due to the abundance of single cell RNA-seq data, a number of methods for predicting expression after perturbation have recently been published. Expression prediction methods are enticing because they promise to answer pressing questions in fields ranging from developmental genetics to cell fate engineering and because they are faster, cheaper, and higher-throughput than their experimental counterparts. However, the absolute and relative accuracy of these methods is poorly characterized, limiting their informed use, their improvement, and the interpretation of their predictions. To address these issues, we created a benchmarking platform that combines a panel of large-scale perturbation datasets with an expression forecasting software engine that encompasses or interfaces to current methods. We used our platform to systematically assess methods, parameters, and sources of auxiliary data. We found that uninformed baseline predictions, which were not always included in prior evaluations, yielded the same or better mean absolute error than benchmarked methods in all test cases. These results cast doubt on the ability of current expression forecasting methods to provide mechanistic insights or to rank hypotheses for experimental follow-up. However, given the rapid pace of innovation in the field, new approaches may yield more accurate expression predictions. Our platform will serve as a neutral benchmark to improve methods and to identify contexts in which expression prediction can succeed.

Integration of single-cell transcriptomes and chromatin landscapes reveals regulatory programs driving pharyngeal organ development.

Maldevelopment of the pharyngeal endoderm, an embryonic tissue critical for patterning of the pharyngeal region and ensuing organogenesis, ultimately contributes to several classes of human developmental syndromes and disorders. Such syndromes are characterized by a spectrum of phenotypes that currently cannot be fully explained by known mutations or genetic variants due to gaps in characterization of critical drivers of normal and dysfunctional development. Despite the disease-relevance of pharyngeal endoderm, we still lack a comprehensive and integrative view of the molecular basis and gene regulatory networks driving pharyngeal endoderm development. To close this gap, we apply transcriptomic and chromatin accessibility single-cell sequencing technologies to generate a multi-omic developmental resource spanning pharyngeal endoderm patterning to the emergence of organ-specific epithelia in the developing mouse embryo. We identify cell-type specific gene regulation, distill GRN models that define developing organ domains, and characterize the role of an immunodeficiency-associated forkhead box transcription factor.

Single-Cell RNA-Sequencing-Based CRISPRi Screening Resolves Molecular Drivers of Early Human Endoderm Development.

Studies in vertebrates have outlined conserved molecular control of definitive endoderm (END) development. However, recent work also shows that key molecular aspects of human END regulation differ even from rodents. Differentiation of human embryonic stem cells (ESCs) to END offers a tractable system to study the molecular basis of normal and defective human-specific END development. Here, we interrogated dynamics in chromatin accessibility during differentiation of ESCs to END, predicting DNA-binding proteins that may drive this cell fate transition. We then combined single-cell RNA-seq with parallel CRISPR perturbations to comprehensively define the loss-of-function phenotype of those factors in END development. Following a few candidates, we revealed distinct impairments in the differentiation trajectories for mediators of TGFß signaling and expose a role for the FOXA2 transcription factor in priming human END competence for human foregut and hepatic END specification. Together, this single-cell functional genomics study provides high-resolution insight on human END development.

Inferring population dynamics from single-cell RNA-sequencing time series data.

Recent single-cell RNA-sequencing studies have suggested that cells follow continuous transcriptomic trajectories in an asynchronous fashion during development. However, observations of cell flux along trajectories are confounded with population size effects in snapshot experiments and are therefore hard to interpret. In particular, changes in proliferation and death rates can be mistaken for cell flux. Here we present pseudodynamics, a mathematical framework that reconciles population dynamics with the concepts underlying developmental trajectories inferred from time-series single-cell data. Pseudodynamics models population distribution shifts across trajectories to quantify selection pressure, population expansion, and developmental potentials. Applying this model to time-resolved single-cell RNA-sequencing of T-cell and pancreatic beta cell maturation, we characterize proliferation and apoptosis rates and identify key developmental checkpoints, data inaccessible to existing approaches.

A Single-Cell Transcriptomic Atlas of Thymus Organogenesis Resolves Cell Types and Developmental Maturation.

Thymus development is critical to the adaptive immune system, yet a comprehensive transcriptional framework capturing thymus organogenesis at single-cell resolution is still needed. We applied single-cell RNA sequencing (RNA-seq) to capture 8 days of thymus development, perturbations of T cell receptor rearrangement, and in vitro organ cultures, producing profiles of 24,279 cells. We resolved transcriptional heterogeneity of developing lymphocytes, and genetic perturbation confirmed T cell identity of conventional and non-conventional lymphocytes. We characterized maturation dynamics of thymic epithelial cells in vivo, classified cell maturation state in a thymic organ culture, and revealed the intrinsic capacity of thymic epithelium to preserve transcriptional regularity despite exposure to exogenous retinoic acid. Finally, by integrating the cell atlas with human genome-wide association study (GWAS) data and autoimmune-disease-related genes, we implicated embryonic thymus-resident cells as possible participants in autoimmune disease etiologies. This resource provides a single-cell transcriptional framework for biological discovery and molecular analysis of thymus organogenesis.