No abstract

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.

SUMMARY 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 fore-gut and hepatic END specification. Together, this single-cell functional genomics study provides high-resolution insight on human END development.

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 β-cell maturation, we characterize proliferation and apoptosis rates and identify key developmental checkpoints, inaccessible to existing approaches.

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.