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

Single-cell transcriptional diversity is a hallmark of developmental potential

Single-cell RNA sequencing (scRNA-seq) is a powerful approach for reconstructing cellular differentiation trajectories. However, inferring both the state and direction of differentiation is challenging. Here, we demonstrate a simple, yet robust, determinant of developmental potential-the number of expressed genes per cell-and leverage this measure of transcriptional diversity to develop a computational framework (CytoTRACE) for predicting differentiation states from scRNA-seq data. When applied to diverse tissue types and organisms, CytoTRACE outperformed previous methods and nearly 19,000 annotated gene sets for resolving 52 experimentally determined developmental trajectories. Additionally, it facilitated the identification of quiescent stem cells and revealed genes that contribute to breast tumorigenesis. This study thus establishes a key RNA-based feature of developmental potential and a platform for delineation of cellular hierarchies.

A cell atlas of human thymic development defines T cell repertoire formation

The thymus provides a nurturing environment for the differentiation and selection of T cells, a process orchestrated by their interaction with multiple thymic cell types. We used single-cell RNA sequencing to create a cell census of the human thymus across the life span and to reconstruct T cell differentiation trajectories and T cell receptor (TCR) recombination kinetics. Using this approach, we identified and located in situ CD8αα⁺ T cell populations, thymic fibroblast subtypes, and activated dendritic cell states. In addition, we reveal a bias in TCR recombination and selection, which is attributed to genomic position and the kinetics of lineage commitment. Taken together, our data provide a comprehensive atlas of the human thymus across the life span with new insights into human T cell development.

The Molecular Anatomy of Mouse Skin during Hair Growth and Rest

Skin homeostasis is orchestrated by dozens of cell types that together direct stem cell renewal, lineage commitment, and differentiation. Here, we use single-cell RNA sequencing and single-molecule RNA FISH to provide a systematic molecular atlas of full-thickness skin, determining gene expression profiles and spatial locations that define 56 cell types and states during hair growth and rest. These findings reveal how the outer root sheath (ORS) and inner hair follicle layers coordinate hair production. We found that the ORS is composed of two intermingling but transcriptionally distinct cell types with differing capacities for interactions with stromal cell types. Inner layer cells branch from transcriptionally uncommitted progenitors, and each lineage differentiation passes through an intermediate state. We also provide an online tool to explore this comprehensive skin cell atlas, including epithelial and stromal cells such as fibroblasts, vascular, and immune cells, to spur further discoveries in skin biology.

Thymic Epithelial Cells Contribute to Thymopoiesis and T Cell Development

The thymus is the primary lymphoid organ responsible for the generation and maturation of T cells. Thymic epithelial cells (TECs) account for the majority of thymic stromal components. They are further divided into cortical and medullary TECs based on their localization within the thymus and are involved in positive and negative selection, respectively. Establishment of self-tolerance in the thymus depends on promiscuous gene expression (pGE) of tissue-restricted antigens (TRAs) by TECs. Such pGE is co-controlled by the autoimmune regulator (Aire) and forebrain embryonic zinc fingerlike protein 2 (Fezf2). Over the past two decades, research has found that TECs contribute greatly to thymopoiesis and T cell development. In turn, signals from T cells regulate the differentiation and maturation of TECs. Several signaling pathways essential for the development and maturation of TECs have been discovered. New technology and animal models have provided important observations on TEC differentiation, development, and thymopoiesis. In this review, we will discuss recent advances in classification, development, and maintenance of TECs and mechanisms that control TEC functions during thymic involution and central tolerance.

Qualitative Changes in Cortical Thymic Epithelial Cells Drive Postpartum Thymic Regeneration

During gestation, sex hormones cause a significant thymic involution which enhances fertility. This thymic involution is rapidly corrected following parturition. As thymic epithelial cells (TECs) are responsible for the regulation of thymopoiesis, we analyzed the sequential phenotypic and transcriptomic changes in TECs during the postpartum period in order to identify mechanisms triggering postpartum thymic regeneration. In particular, we performed flow cytometry analyses and deep RNA-sequencing on purified TEC subsets at several time points before and after parturition. We report that pregnancy-induced involution is not caused by loss of TECs since their number does not change during or after pregnancy. However, during pregnancy, we observed a significant depletion of all thymocyte subsets downstream of the double-negative 1 (DN1) differentiation stage. Variations in thymocyte numbers correlated with conspicuous changes in the transcriptome of cortical TECs (cTECs). The transcriptomic changes affected predominantly cTEC expression of Foxn1, its targets and several genes that are essential for thymopoiesis. By contrast, medullary TECs (mTECs) showed very little transcriptomic changes in the early postpartum regenerative phase, but seemed to respond to the expansion of single-positive (SP) thymocytes in the late phase of regeneration. Together, these results show that postpartum thymic regeneration is orchestrated by variations in expression of a well-defined subset of cTEC genes, that occur very early after parturition.

Single-Cell Expression Variability Implies Cell Function

As single-cell RNA sequencing (scRNA-seq) data becomes widely available, cell-to-cell variability in gene expression, or single-cell expression variability (scEV), has been increasingly appreciated. However, it remains unclear whether this variability is functionally important and, if so, what are its implications for multi-cellular organisms. Here, we analyzed multiple scRNA-seq data sets from lymphoblastoid cell lines (LCLs), lung airway epithelial cells (LAECs), and dermal fibroblasts (DFs) and, for each cell type, selected a group of homogenous cells with highly similar expression profiles. We estimated the scEV levels for genes after correcting the mean-variance dependency in that data and identified 465, 466, and 364 highly variable genes (HVGs) in LCLs, LAECs, and DFs, respectively. Functions of these HVGs were found to be enriched with those biological processes precisely relevant to the corresponding cell type’s function, from which the scRNA-seq data used to identify HVGs were generated—e.g., cytokine signaling pathways were enriched in HVGs identified in LCLs, collagen formation in LAECs, and keratinization in DFs. We repeated the same analysis with scRNA-seq data from induced pluripotent stem cells (iPSCs) and identified only 79 HVGs with no statistically significant enriched functions; the overall scEV in iPSCs was of negligible magnitude. Our results support the “variation is function” hypothesis, arguing that scEV is required for cell type-specific, higher-level system function. Thus, quantifying and characterizing scEV are of importance for our understating of normal and pathological cellular processes.

Thymic epithelial cell heterogeneity: TEC by TEC

The generation of a functional T cell repertoire in the thymus is mainly orchestrated by thymic epithelial cells (TECs), which provide developing T cells with cues for their navigation, proliferation, differentiation and survival. The TEC compartment has been segregated historically into two major populations of medullary TECs and cortical TECs, which differ in their anatomical localization, molecular characteristics and functional roles. However, recent studies have shown that TECs are highly heterogeneous and comprise multiple subpopulations with distinct molecular and functional characteristics, including tuft cell-like or corneocyte-like phenotypes. Here, we review the most recent advances in our understanding of TEC heterogeneity from a molecular, functional and developmental perspective. In particular, we highlight the key insights that were recently provided by single-cell genomic technologies and in vivo fate mapping and discuss them in the context of previously published data.

Myc controls a distinct transcriptional program in fetal thymic epithelial cells that determines thymus growth

Interactions between thymic epithelial cells (TEC) and developing thymocytes are essential for T cell development, but molecular insights on TEC and thymus homeostasis are still lacking. Here we identify distinct transcriptional programs of TEC that account for their age-specific properties, including proliferation rates, engraftability and function. Further analyses identify Myc as a regulator of fetal thymus development to support the rapid increase of thymus size during fetal life. Enforced Myc expression in TEC induces the prolonged maintenance of a fetal-specific transcriptional program, which in turn extends the growth phase of the thymus and enhances thymic output; meanwhile, inducible expression of Myc in adult TEC similarly promotes thymic growth. Mechanistically, this Myc function is associated with enhanced ribosomal biogenesis in TEC. Our study thus identifies age-specific transcriptional programs in TEC, and establishes that Myc controls thymus size.

Building a Human Thymus: A Pointillist View

In this issue of Immunity, Zeng et al. use single-cell RNA sequencing analyses of rare samples to shed light on the emergence of thymic stromal cell types, the first developing T lymphocytes, and their possible pre-thymic precursors in the early human fetus.

Single-Cell RNA Sequencing Resolves Spatiotemporal Development of Pre-thymic Lymphoid Progenitors and Thymus Organogenesis in Human Embryos

Generation of the first T lymphocytes in the human embryo involves the emergence, migration, and thymus seeding of lymphoid progenitors together with concomitant thymus organogenesis, which is the initial step to establish the entire adaptive immune system. However, the cellular and molecular programs regulating this process remain unclear. We constructed a single-cell transcriptional landscape of human early T lymphopoiesis by using cells from multiple hemogenic and hematopoietic sites spanning embryonic and fetal stages. Among heterogenous early thymic progenitors, one subtype shared common features with a subset of lymphoid progenitors in fetal liver that are known as thymus-seeding progenitors. Unbiased bioinformatics analysis identified a distinct type of pre-thymic lymphoid progenitors in the aorta-gonad-mesonephros (AGM) region. In parallel, we investigated thymic epithelial cell development and potential cell-cell interactions during thymus organogenesis. Together, our data provide insights into human early T lymphopoiesis that prospectively direct T lymphocyte regeneration, which might lead to development of clinical applications.

BBKNN: fast batch alignment of single cell transcriptomes

MOTIVATION

Increasing numbers of large scale single cell RNA-Seq projects are leading to a data explosion, which can only be fully exploited through data integration. A number of methods have been developed to combine diverse datasets by removing technical batch effects, but most are computationally intensive. To overcome the challenge of enormous datasets, we have developed BBKNN, an extremely fast graph-based data integration algorithm. We illustrate the power of BBKNN on large scale mouse atlasing data, and favourably benchmark its run time against a number of competing methods.

AVAILABILITY AND IMPLEMENTATION

BBKNN is available at https://github.com/Teichlab/bbknn, along with documentation and multiple example notebooks, and can be installed from pip.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

BBKNN: fast batch alignment of single cell transcriptomes

Motivation

Increasing numbers of large scale single cell RNA-Seq projects are leading to a data explosion, which can only be fully exploited through data integration. A number of methods have been developed to combine diverse datasets by removing technical batch effects, but most are computationally intensive. To overcome the challenge of enormous datasets, we have developed BBKNN, an extremely fast graph-based data integration algorithm. We illustrate the power of BBKNN on large scale mouse atlasing data, and favourably benchmark its run time against a number of competing methods.

Availability and implementation

BBKNN is available at https://github.com/Teichlab/bbknn, along with documentation and multiple example notebooks, and can be installed from pip.

Supplementary information

Supplementary data are available at Bioinformatics online.

Cell-type-specific role of lamin-B1 in thymus development and its inflammation-driven reduction in thymus aging

Cellular architectural proteins often participate in organ development and maintenance. Although functional decay of some of these proteins during aging is known, the cell-type-specific developmental role and the cause and consequence of their subsequent decay remain to be established especially in mammals. By studying lamins, the nuclear structural proteins, we demonstrate that lamin-B1 functions specifically in the thymic epithelial cells (TECs) for proper thymus organogenesis. An up-regulation of proinflammatory cytokines in the intra-thymic myeloid immune cells during aging accompanies a gradual reduction of lamin-B1 in adult TECs. We show that these cytokines can cause senescence and lamin-B1 reduction of the young adult TECs. Lamin-B1 supports the expression of TEC genes that can help maintain the adult TEC subtypes we identified by single-cell RNA-sequencing, thymic architecture, and function. Thus, structural proteins involved in organ building and maintenance can undergo inflammation-driven decay which can in turn contribute to age-associated organ degeneration.

Forkhead box K2 promotes human colorectal cancer metastasis by upregulating ZEB1 and EGFR

Background: Metastasis is the major reason for high recurrence rates and poor survival among patients with colorectal cancer (CRC). However, the underlying molecular mechanism of CRC metastasis is unclear. This study aimed to investigate the role of forkhead box K2 (FOXK2), one of the most markedly increased FOX genes in CRC, and the mechanism by which it is deregulated in CRC metastasis. Methods: FOXK2 levels were analyzed in two independent human CRC cohorts (cohort I, n = 363; cohort II, n = 390). In vitro Transwell assays and in vivo lung and liver metastasis models were used to examine CRC cell migration, invasion and metastasis. Chromatin immunoprecipitation and luciferase reporter assays were used to measure the binding of transcription factors to the promoters of FOXK2, zinc finger E-box binding homeobox 1 (ZEB1) and epidermal growth factor receptor (EGFR). Cetuximab was utilized to treat FOXK2-mediated metastatic CRC. Results: FOXK2 was significantly upregulated in human CRC tissues, was correlated with more aggressive features and indicated a poor prognosis. FOXK2 overexpression promoted CRC migration, invasion and metastasis, while FOXK2 downregulation had the opposite effects. ZEB1 and EGFR were determined to be direct transcriptional targets of FOXK2 and were essential for FOXK2-mediated CRC metastasis. Moreover, activation of EGFR signaling by EGF enhanced FOXK2 expression via the extracellular regulated protein kinase (ERK) and nuclear factor (NF)-κB pathways. The EGFR monoclonal antibody cetuximab significantly inhibited FOXK2-promoted CRC metastasis. In clinical CRC tissues, FOXK2 expression was positively correlated with the expression of p65, ZEB1 and EGFR. CRC patients who coexpressed p65/FOXK2, FOXK2/ZEB1 and FOXK2/EGFR had poorer prognosis. Conclusions: FOXK2 serves as a prognostic biomarker in CRC. Cetuximab can block the EGF-NF-κB-FOXK2-EGFR feedback loop and suppress CRC metastasis.

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.

Suppression of ILC2 differentiation from committed T cell precursors by E protein transcription factors

Qian et al. shows that ILC2s can be generated from not only thymic multipotent progenitors but also committed T cell precursors. These processes are greatly suppressed by E protein transcription factors. Thymic ILC2s show functional differences from those made elsewhere.

Current models propose that group 2 innate lymphoid cells (ILC2s) are generated in the bone marrow. Here, we demonstrate that subsets of these cells can differentiate from multipotent progenitors and committed T cell precursors in the thymus, both in vivo and in vitro. These thymic ILC2s exit the thymus, circulate in the blood, and home to peripheral tissues. Ablation of E protein transcription factors greatly promotes the ILC fate while impairing B and T cell development. Consistently, a transcriptional network centered on the ZBTB16 transcription factor and IL-4 signaling pathway is highly up-regulated due to E protein deficiency. Our results show that ILC2 can still arise from what are normally considered to be committed T cell precursors, and that this alternative cell fate is restrained by high levels of E protein activity in these cells. Thymus-derived lung ILC2s of E protein–deficient mice show different transcriptomes, proliferative properties, and cytokine responses from wild-type counterparts, suggesting potentially distinct functions.

Staged developmental mapping and X chromosome transcriptional dynamics during mouse spermatogenesis

Male gametes are generated through a specialised differentiation pathway involving a series of developmental transitions that are poorly characterised at the molecular level. Here, we use droplet-based single-cell RNA-Sequencing to profile spermatogenesis in adult animals and at multiple stages during juvenile development. By exploiting the first wave of spermatogenesis, we both precisely stage germ cell development and enrich for rare somatic cell-types and spermatogonia. To capture the full complexity of spermatogenesis including cells that have low transcriptional activity, we apply a statistical tool that identifies previously uncharacterised populations of leptotene and zygotene spermatocytes. Focusing on post-meiotic events, we characterise the temporal dynamics of X chromosome re-activation and profile the associated chromatin state using CUT&RUN. This identifies a set of genes strongly repressed by H3K9me3 in spermatocytes, which then undergo extensive chromatin remodelling post-meiosis, thus acquiring an active chromatin state and spermatid-specific expression.

The Role of the Thymus in the Immune Response

The thymus is a primary lymphoid organ essential for the development of T lymphocytes, which orchestrate adaptive immune responses. T-cell development in the thymus is spatially regulated; key checkpoints in T-cell maturation and selection occur in cortical and medullary regions to eliminate self-reactive T cells, establish central tolerance, and export naïve T cells to the periphery with the potential to recognize diverse pathogens. Thymic output is also temporally regulated due to age-related involution of the thymus accompanied by loss of epithelial cells. This review discusses the structural and age-related control of thymus function in humans.

Quantitative imaging of lipid droplets in single cells

The combination of next generation sequencing (NGS) and automated liquid handling platforms has led to a revolution in single-cell genomic studies. However, many molecules that are critical to understanding the functional roles of cells in a complex tissue or organs, are not directly encoded in the genome, and therefore cannot be profiled with NGS. Lipids, for example, play a critical role in many metabolic processes but cannot be detected by sequencing. Recent developments in quantitative imaging, particularly coherent Raman scattering (CRS) techniques, have produced a suite of tools for studying lipid content in single cells. This article reviews CRS imaging and computational image processing techniques for non-destructive profiling of dynamic changes in lipid composition and spatial distribution at the single-cell level. As quantitative CRS imaging progresses synergistically with microfluidic and microscopic platforms for single-cell genomic analysis, we anticipate that these techniques will bring researchers closer towards combined lipidomic and genomic analysis.

Single Cell Gene Expression to Understand the Dynamic Architecture of the Heart

The recent development of single cell gene expression technologies, and especially single cell transcriptomics, have revolutionized the way biologists and clinicians investigate organs and organisms, allowing an unprecedented level of resolution to the description of cell demographics in both healthy and diseased states. Single cell transcriptomics provide information on prevalence, heterogeneity, and gene co-expression at the individual cell level. This enables a cell-centric outlook to define intracellular gene regulatory networks and to bridge toward the definition of intercellular pathways otherwise masked in bulk analysis. The technologies have developed at a fast pace producing a multitude of different approaches, with several alternatives to choose from at any step, including single cell isolation and capturing, lysis, RNA reverse transcription and cDNA amplification, library preparation, sequencing, and computational analyses. Here, we provide guidelines for the experimental design of single cell RNA sequencing experiments, exploring the current options for the crucial steps. Furthermore, we provide a complete overview of the typical data analysis workflow, from handling the raw sequencing data to making biological inferences. Significantly, advancements in single cell transcriptomics have already contributed to outstanding exploratory and functional studies of cardiac development and disease models, as summarized in this review. In conclusion, we discuss achievable outcomes of single cell transcriptomics' applications in addressing unanswered questions and influencing future cardiac clinical applications.

Single-Cell Analysis Identifies Thymic Maturation Delay in Growth-Restricted Neonatal Mice

Fetal growth restriction (FGR) causes a wide variety of defects in the neonate which can lead to increased risk of heart disease, diabetes, anxiety and other disorders later in life. However, the effect of FGR on the immune system, is poorly understood. We used a well-characterized mouse model of FGR in which placental Igf-2 production is lost due to deletion of the placental specific Igf-2 P₀ promotor. The thymi in such animals were reduced in mass with a ~70% reduction in cellularity. We used single cell RNA sequencing (Drop-Seq) to analyze 7,264 thymus cells collected at postnatal day 6. We identified considerable heterogeneity among the Cd8/Cd4 double positive cells with one subcluster showing marked upregulation of transcripts encoding a sub-set of proteins that contribute to the surface of the ribosome. The cells from the FGR animals were underrepresented in this cluster. Furthermore, the distribution of cells from the FGR animals was skewed with a higher proportion of immature double negative cells and fewer mature T-cells. Cell cycle regulator transcripts also varied across clusters. The T-cell deficit in FGR mice persisted into adulthood, even when body and organ weights approached normal levels due to catch-up growth. This finding complements the altered immunity found in growth restricted human infants. This reduction in T-cellularity may have implications for adult immunity, adding to the list of adult conditions in which the in utero environment is a contributory factor.