Single-Cell RNA-Seq Mapping of Human Thymopoiesis Reveals Lineage Specification Trajectories and a Commitment Spectrum in T Cell Development

KLRG1 Cell Depletion as a Novel Therapeutic Strategy in Patients with Mature T-Cell Lymphoma Subtypes

PURPOSE

Develop a novel therapeutic strategy for patients with subtypes of mature T-cell and NK-cell neoplasms.

EXPERIMENTAL DESIGN

Primary specimens, cell lines, patient-derived xenograft models, commercially available, and proprietary anti-KLRG1 antibodies were used for screening, target, and functional validation.

RESULTS

Here we demonstrate that surface KLRG1 is highly expressed on tumor cells in subsets of patients with extranodal NK/T-cell lymphoma (ENKTCL), T-prolymphocytic leukemia (T-PLL), and gamma/delta T-cell lymphoma (G/D TCL). The majority of the CD8+/CD57+ or CD3-/CD56+ leukemic cells derived from patients with T- and NK-large granular lymphocytic leukemia (T-LGLL and NK-LGLL), respectively, expressed surface KLRG1. The humanized afucosylated anti-KLRG1 monoclonal antibody (mAb208) optimized for mouse in vivo use depleted KLRG1+ TCL cells by mechanisms of ADCC, ADCP, and CDC rather than apoptosis. mAb208 induced ADCC and ADCP of T-LGLL patient-derived CD8+/CD57+ cells ex vivo. mAb208 effected ADCC of subsets of healthy donor-derived KLRG1+ NK, CD4+, CD8+ Tem, and TemRA cells while sparing KLRG1- naïve and CD8+ Tcm cells. Treatment of cell line and TCL patient-derived xenografts with mAb208 or anti-CD47 mAb alone and in combination with the PI3K-δ/γ inhibitor duvelisib extended survival. The depletion of macrophages in vivo antagonized mAb208 efficacy.

CONCLUSIONS

Our findings suggest the potential benefit of a broader treatment strategy combining therapeutic antibodies with PI3Ki for the treatment of patients with mature T-cell and NK-cell neoplasms. See related commentary by Varma and Diefenbach, p. 2300.

Dynamic Changes in Lymphocyte Populations Establish Zebrafish as a Thymic Involution Model

The thymus is the site of T lymphocyte development and T cell education to recognize foreign, but not self, Ags. B cells also reside and develop in the thymus, although their functions are less clear. During "thymic involution," a process of lymphoid atrophy and adipose replacement linked to sexual maturation, thymocytes decline. However, thymic B cells decrease far less than T cells, such that B cells comprise ∼1% of human neonatal thymocytes but up to ∼10% in adults. All jawed vertebrates possess a thymus, and we and others have shown zebrafish (Danio rerio) also have thymic B cells. In this article, we investigated the precise identities of zebrafish thymic T and B cells and how they change with involution. We assessed the timing and specific details of zebrafish thymic involution using multiple lymphocyte-specific, fluorophore-labeled transgenic lines, quantifying the changes in thymic T- and B-lymphocytes pre- versus postinvolution. Our results prove that, as in humans, zebrafish thymic B cells increase relative to T cells postinvolution. We also performed RNA sequencing on D. rerio thymic and marrow lymphocytes of four novel double-transgenic lines, identifying distinct populations of immature T and B cells. Collectively, this is, to our knowledge, the first comprehensive analysis of zebrafish thymic involution, demonstrating its similarity to human involution and establishing the highly genetically manipulatable zebrafish model as a template for involution studies.

Transposable elements regulate thymus development and function

Transposable elements (TEs) are repetitive sequences representing ~45% of the human and mouse genomes and are highly expressed by medullary thymic epithelial cells (mTECs). In this study, we investigated the role of TEs on T-cell development in the thymus. We performed multiomic analyses of TEs in human and mouse thymic cells to elucidate their role in T-cell development. We report that TE expression in the human thymus is high and shows extensive age- and cell lineage-related variations. TE expression correlates with multiple transcription factors in all cell types of the human thymus. Two cell types express particularly broad TE repertoires: mTECs and plasmacytoid dendritic cells (pDCs). In mTECs, transcriptomic data suggest that TEs interact with transcription factors essential for mTEC development and function (e.g., PAX1 and REL), and immunopeptidomic data showed that TEs generate MHC-I-associated peptides implicated in thymocyte education. Notably, AIRE, FEZF2, and CHD4 regulate small yet non-redundant sets of TEs in murine mTECs. Human thymic pDCs homogenously express large numbers of TEs that likely form dsRNA, which can activate innate immune receptors, potentially explaining why thymic pDCs constitutively secrete IFN ɑ/β. This study highlights the diversity of interactions between TEs and the adaptive immune system. TEs are genetic parasites, and the two thymic cell types most affected by TEs (mTEcs and pDCs) are essential to establishing central T-cell tolerance. Therefore, we propose that orchestrating TE expression in thymic cells is critical to prevent autoimmunity in vertebrates.

PAGER-scFGA: unveiling cell functions and molecular mechanisms in cell trajectories through single-cell functional genomics analysis

Background: Understanding how cells and tissues respond to stress factors and perturbations during disease processes is crucial for developing effective prevention, diagnosis, and treatment strategies. Single-cell RNA sequencing (scRNA-seq) enables high-resolution identification of cells and exploration of cell heterogeneity, shedding light on cell differentiation/maturation and functional differences. Recent advancements in multimodal sequencing technologies have focused on improving access to cell-specific subgroups for functional genomics analysis. To facilitate the functional annotation of cell groups and characterization of molecular mechanisms underlying cell trajectories, we introduce the Pathways, Annotated Gene Lists, and Gene Signatures Electronic Repository for Single-Cell Functional Genomics Analysis (PAGER-scFGA). Results: We have developed PAGER-scFGA, which integrates cell functional annotations and gene-set enrichment analysis into popular single-cell analysis pipelines such as Scanpy. Using differentially expressed genes (DEGs) from pairwise cell clusters, PAGER-scFGA infers cell functions through the enrichment of potential cell-marker genesets. Moreover, PAGER-scFGA provides pathways, annotated gene lists, and gene signatures (PAGs) enriched in specific cell subsets with tissue compositions and continuous transitions along cell trajectories. Additionally, PAGER-scFGA enables the construction of a gene subcellular map based on DEGs and allows examination of the gene functional compartments (GFCs) underlying cell maturation/differentiation. In a real-world case study of mouse natural killer (mNK) cells, PAGER-scFGA revealed two major stages of natural killer (NK) cells and three trajectories from the precursor stage to NK T-like mature stage within blood, spleen, and bone marrow tissues. As the trajectories progress to later stages, the DEGs exhibit greater divergence and variability. However, the DEGs in different trajectories still interact within a network during NK cell maturation. Notably, PAGER-scFGA unveiled cell cytotoxicity, exocytosis, and the response to interleukin (IL) signaling pathways and associated network models during the progression from precursor NK cells to mature NK cells. Conclusion: PAGER-scFGA enables in-depth exploration of functional insights and presents a comprehensive knowledge map of gene networks and GFCs, which can be utilized for future studies and hypothesis generation. It is expected to become an indispensable tool for inferring cell functions and detecting molecular mechanisms within cell trajectories in single-cell studies. The web app (accessible at https://au-singlecell.streamlit.app/) is publicly available.

Mechanism study of ubiquitination in T cell development and autoimmune disease

T cells play critical role in multiple immune processes including antigen response, tumor immunity, inflammation, self-tolerance maintenance and autoimmune diseases et. Fetal liver or bone marrow-derived thymus-seeding progenitors (TSPs) settle in thymus and undergo T cell-lineage commitment, proliferation, T cell receptor (TCR) rearrangement, and thymic selections driven by microenvironment composed of thymic epithelial cells (TEC), dendritic cells (DC), macrophage and B cells, thus generating T cells with diverse TCR repertoire immunocompetent but not self-reactive. Additionally, some self-reactive thymocytes give rise to Treg with the help of TEC and DC, serving for immune tolerance. The sequential proliferation, cell fate decision, and selection during T cell development and self-tolerance establishment are tightly regulated to ensure the proper immune response without autoimmune reaction. There are remarkable progresses in understanding of the regulatory mechanisms regarding ubiquitination in T cell development and the establishment of self-tolerance in the past few years, which holds great potential for further therapeutic interventions in immune-related diseases.

Transcriptional profile of human thymus reveals IGFBP5 is correlated with age-related thymic involution

Thymus is the main immune organ which is responsible for the production of self-tolerant and functional T cells, but it shrinks rapidly with age after birth. Although studies have researched thymus development and involution in mouse, the critical regulators that arise with age in human thymus remain unclear. We collected public human single-cell transcriptomic sequencing (scRNA-seq) datasets containing 350,678 cells from 36 samples, integrated them as a cell atlas of human thymus. Clinical samples were collected and experiments were performed for validation. We found early thymocyte-specific signaling and regulons which played roles in thymocyte migration, proliferation, apoptosis and differentiation. Nevertheless, signaling patterns including number, strength and path completely changed during aging, Transcription factors (FOXC1, MXI1, KLF9, NFIL3) and their target gene, IGFBP5, were resolved and up-regulated in aging thymus and involved in promoting epithelial-mesenchymal transition (EMT), responding to steroid and adipogenesis process of thymic epithelial cell (TECs). Furthermore, we validated that IGFBP5 protein increased at TECs and Hassall's corpuscle in both human and mouse aging thymus and knockdown of IGFBP5 significantly increased the expression of proliferation-related genes in thymocytes. Collectively, we systematically explored cell-cell communications and regulons of early thymocytes as well as age-related differences in human thymus by using both bioinformatic and experimental verification, indicating IGFBP5 as a functional marker of thymic involution and providing new insights into the mechanisms of thymus involution.

Lineage tracing of T cell differentiation from T-iPSC by 2D feeder-free culture and 3D organoid culture

Introduction

T cells induced from induced pluripotent stem cells(iPSCs) derived from antigen-specific T cells (T-iPS-T cells) are an attractive tool for T cell immunotherapy. The induction of cytotoxic T-iPS-T cells is well established in feeder-free condition for the aim of off-the-shelf production, however, the induction of helper T-iPS-T cells remains challenging.

Methods

We analyzed T-iPS-T cells matured in 3D organoid culture at different steps in the culture process at the single-cell level. T-iPS-T cell datasets were merged with an available human thymocyte dataset based in single-cell RNA sequencing (scRNA-seq). Particularly, we searched for genes crucial for generation CD4+ T-iPS-T cells by comparing T-iPS-T cells established in 2D feeder-free or 3D organoid culture.

Results

The scRNA-seq data indicated that T-iPS-T cells are similar to T cells transitioning to human thymocytes, with SELENOW, GIMAP4, 7, SATB1, SALMF1, IL7R, SYTL2, S100A11, STAT1, IFITM1, LZTFL1 and SOX4 identified as candidate genes for the 2D feeder-free induction of CD4+ T-iPS-T cells.

Discussion

This study provides single cell transcriptome datasets of iPS-T cells and leads to further analysis for CD4+ T cell generation from T-iPSCs.

A spatial human thymus cell atlas mapped to a continuous tissue axis

T cells develop from circulating precursors, which enter the thymus and migrate throughout specialised sub-compartments to support maturation and selection. This process starts already in early fetal development and is highly active until the involution of the thymus in adolescence. To map the micro-anatomical underpinnings of this process in pre- vs. post-natal states, we undertook a spatially resolved analysis and established a new quantitative morphological framework for the thymus, the Cortico-Medullary Axis. Using this axis in conjunction with the curation of a multimodal single-cell, spatial transcriptomics and high-resolution multiplex imaging atlas, we show that canonical thymocyte trajectories and thymic epithelial cells are highly organised and fully established by post-conception week 12, pinpoint TEC progenitor states, find that TEC subsets and peripheral tissue genes are associated with Hassall's Corpuscles and uncover divergence in the pace and drivers of medullary entry between CD4 vs. CD8 T cell lineages. These findings are complemented with a holistic toolkit for spatial analysis and annotation, providing a basis for a detailed understanding of T lymphocyte development.

Divergent molecular events underlying initial T-cell commitment in human prenatal and postnatal thymus

Introduction

Intrathymic T-cell development is a coordinated process accompanied by dynamic changes in gene expression. Although the transcriptome characteristics of developing T cells in both human fetal and postnatal thymus at single-cell resolution have been revealed recently, the differences between human prenatal and postnatal thymocytes regarding the ontogeny and early events of T-cell development still remain obscure. Moreover, the transcriptional heterogeneity and posttranscriptional gene expression regulation such as alternative polyadenylation at different stages are also unknown.

Method

In this study, we performed integrative single-cell analyses of thymocytes at distinct developmental stages.

Results

The subsets of prenatal CD4-CD8- double-negative (DN) cells, the most immature thymocytes responsible for T-cell lineage commitment, were characterized. By comprehensively comparing prenatal and postnatal DN cells, we revealed significant differences in some key gene expressions. Specifically, prenatal DN subpopulations exhibited distinct biological processes and markedly activated several metabolic programs that may be coordinated to meet the required bioenergetic demands. Although showing similar gene expression patterns along the developmental path, prenatal and postnatal thymocytes were remarkably varied regarding the expression dynamics of some pivotal genes for cell cycle, metabolism, signaling pathway, thymus homing, and T-cell commitment. Finally, we quantified the transcriptome-wide changes in alternative polyadenylation across T-cell development and found diverse preferences of polyadenylation site usage in divergent populations along the T-cell commitment trajectory.

Discussion

In summary, our results revealed transcriptional heterogeneity and a dynamic landscape of alternative polyadenylation during T-cell development in both human prenatal and postnatal thymus, providing a comprehensive resource for understanding T lymphopoiesis in human thymus.

Circulating IRF8-expressing CD123+CD127+ lymphoid progenitors: key players in human hematopoiesis

Lymphopoiesis is the process in which B and T cells, and innate lymphoid cells (ILCs) develop from hematopoietic progenitors that exhibit early lymphoid priming. The branching points where lymphoid-primed human progenitors are further specified to B/T/ILC differentiation trajectories remain unclear. Here, we discuss the emerging role of interferon regulatory factor (IRF)8 as a key factor to bridge human lymphoid and dendritic cell (DC) differentiation, and the current evidence for the existence of circulating and tissue-resident CD123+CD127+ lymphoid progenitors. We propose a model whereby DC/B/T/ILC lineage programs in circulating CD123+CD127+ lymphoid progenitors are expressed in balance. Upon tissue seeding, the tissue microenvironment tilts this molecular balance towards a specific lineage, thereby determining in vivo lineage fates. Finally, we discuss the translational implication of these lymphoid precursors.

Dynamic Changes in Lymphocyte Populations Establish Zebrafish as a Thymic Involution Model

The thymus is the site of T lymphocyte development and T cell education to recognize foreign, but not self, antigens. B cells also reside and develop in the thymus, although their functions are less clear. During 'thymic involution,' a process of lymphoid atrophy and adipose replacement linked to sexual maturation, thymocytes decline. However, thymic B cells decrease far less than T cells, such that B cells comprise ~1% of human neonatal thymocytes, but up to ~10% in adults. All jawed vertebrates possess a thymus, and we and others have shown zebrafish (Danio rerio) also have thymic B cells. Here, we investigated the precise identities of zebrafish thymic T and B cells and how they change with involution. We assessed the timing and specific details of zebrafish thymic involution using multiple lymphocyte-specific, fluorophore-labeled transgenic lines, quantifying the changes in thymic T- and B-lymphocytes pre- vs. post-involution. Our results prove that, as in humans, zebrafish thymic B cells increase relative to T cells post-involution. We also performed RNA sequencing (RNA-seq) on D. rerio thymic and marrow lymphocytes of four novel double-transgenic lines, identifying distinct populations of immature T and B cells. Collectively, this is the first comprehensive analysis of zebrafish thymic involution, demonstrating its similarity to human involution, and establishing the highly genetically-manipulatable zebrafish model as a template for involution studies.

Malignant A-to-I RNA editing by ADAR1 drives T-cell acute lymphoblastic leukemia relapse via attenuating dsRNA sensing

Leukemia initiating cells (LICs) are regarded as the origin of leukemia relapse and therapeutic resistance. Identifying direct stemness determinants that fuel LIC self-renewal is critical for developing targeted approaches to eliminate LICs and prevent relapse. Here, we show that the RNA editing enzyme ADAR1 is a crucial stemness factor that promotes LIC self-renewal by attenuating aberrant double-stranded RNA (dsRNA) sensing. Elevated adenosine-to-inosine (A-to-I) editing is a common attribute of relapsed T-ALL regardless of molecular subtypes. Consequently, knockdown of ADAR1 severely inhibits LIC self-renewal capacity and prolongs survival in T-ALL PDX models. Mechanistically, ADAR1 directs hyper-editing of immunogenic dsRNA and retains unedited nuclear dsRNA to avoid detection by the innate immune sensor MDA5. Moreover, we uncovered that the cell intrinsic level of MDA5 dictates the dependency on ADAR1-MDA5 axis in T-ALL. Collectively, our results show that ADAR1 functions as a self-renewal factor that limits the sensing of endogenous dsRNA. Thus, targeting ADAR1 presents a safe and effective therapeutic strategy for eliminating T-ALL LICs.

Systems biology approaches to unravel lymphocyte subsets and function

Single-cell technologies have revealed the extensive heterogeneity and complexity of the immune system. Systems biology approaches in immunology have taken advantage of the high-parameter, high-throughput data and analyzed immune cell types in a 'bottom-up' data-driven method. This approach has discovered previously unrecognized cell types and functions. Especially for human immunology, in which experimental manipulations are challenging, systems approach has become a successful means to investigate physiologically relevant contexts. This review focuses on the recent findings in lymphocyte biology, from their development, differentiation into subsets, and heterogeneity in their functions, enabled by these systems approaches. Furthermore, we review examples of the application of findings from systems approach studies and discuss how now to leave the rich dataset in the curse of high dimensionality.

Speed and navigation control of thymocyte development by the fetal T-cell gene regulatory network

T-cell differentiation is a tightly regulated developmental program governed by interactions between transcription factors (TFs) and chromatin landscapes and affected by signals received from the thymic stroma. This process is marked by a series of checkpoints: T-lineage commitment, T-cell receptor (TCR)β selection, and positive and negative selection. Dynamically changing combinations of TFs drive differentiation along the T-lineage trajectory, through mechanisms that have been most extensively dissected in adult mouse T-lineage cells. However, fetal T-cell development differs from adult in ways that suggest that these TF mechanisms are not fully deterministic. The first wave of fetal T-cell differentiation occurs during a unique developmental window during thymic morphogenesis, shows more rapid kinetics of differentiation with fewer rounds of cell division, and gives rise to unique populations of innate lymphoid cells (ILCs) and invariant γδT cells that are not generated in the adult thymus. As the characteristic kinetics and progeny biases are cell-intrinsic properties of thymic progenitors, the differences could be based on distinct TF network circuitry within the progenitors themselves. Here, we review recent single-cell transcriptome data that illuminate the TF networks involved in T-cell differentiation in the fetal and adult mouse thymus.

Intrathymic dendritic cell-biased precursors promote human T cell lineage specification through IRF8-driven transmembrane TNF

The cross-talk between thymocytes and thymic stromal cells is fundamental for T cell development. In humans, intrathymic development of dendritic cells (DCs) is evident but its physiological significance is unknown. Here we showed that DC-biased precursors depended on the expression of the transcription factor IRF8 to express the membrane-bound precursor form of the cytokine TNF (tmTNF) to promote differentiation of thymus seeding hematopoietic progenitors into T-lineage specified precursors through activation of the TNF receptor (TNFR)-2 instead of TNFR1. In vitro recapitulation of TNFR2 signaling by providing low-density tmTNF or a selective TNFR2 agonist enhanced the generation of human T cell precursors. Our study shows that, in addition to mediating thymocyte selection and maturation, DCs function as hematopoietic stromal support for the early stages of human T cell development and provide proof of concept that selective targeting of TNFR2 can enhance the in vitro generation of T cell precursors for clinical application.

Integrated role of human thymic stromal cells in hematopoietic stem cell extravasation

The human thymus is the site of T-cell maturation and induction of central tolerance. Hematopoietic stem cell (HSC)-derived progenitors are recruited to the thymus from the fetal liver during early prenatal development and from bone marrow at later stages and postnatal life. The mechanism by which HSCs are recruited to the thymus is poorly understood in humans, though mouse models have indicated the critical role of thymic stromal cells (TSC). Here, we developed a 3D microfluidic assay based on human cells to model HSC extravasation across the endothelium into the extracellular matrix. We found that the presence of human TSC consisting of cultured thymic epithelial cells (TEC) and interstitial cells (TIC) increases the HSC extravasation rates by 3-fold. Strikingly, incorporating TEC or TIC alone is insufficient to perturb HSC extravasation rates. Furthermore, we identified complex gene expressions from interactions between endothelial cells, TEC and TIC modulates the HSCs extravasation. Our results suggest that comprehensive signaling from the complex thymic microenvironment is crucial for thymus seeding and that our system will allow manipulation of these signals with the potential to increase thymocyte migration in a therapeutic setting.

Endothelial cell-leukemia interactions remodel drug responses, uncovering T-ALL vulnerabilities

T-cell acute lymphoblastic leukemia (T-ALL) is an aggressive and often incurable disease. To uncover therapeutic vulnerabilities, we first developed T-ALL patient-derived tumor xenografts (PDXs) and exposed PDX cells to a library of 433 clinical-stage compounds in vitro. We identified 39 broadly active drugs with antileukemia activity. Because endothelial cells (ECs) can alter drug responses in T-ALL, we developed an EC/T-ALL coculture system. We found that ECs provide protumorigenic signals and mitigate drug responses in T-ALL PDXs. Whereas ECs broadly rescued several compounds in most models, for some drugs the rescue was restricted to individual PDXs, suggesting unique crosstalk interactions and/or intrinsic tumor features. Mechanistically, cocultured T-ALL cells and ECs underwent bidirectional transcriptomic changes at the single-cell level, highlighting distinct "education signatures." These changes were linked to bidirectional regulation of multiple pathways in T-ALL cells as well as in ECs. Remarkably, in vitro EC-educated T-ALL cells transcriptionally mirrored ex vivo splenic T-ALL at single-cell resolution. Last, 5 effective drugs from the 2 drug screenings were tested in vivo and shown to effectively delay tumor growth and dissemination thus prolonging overall survival. In sum, we developed a T-ALL/EC platform that elucidated leukemia-microenvironment interactions and identified effective compounds and therapeutic vulnerabilities.

Multi-modular structure of the gene regulatory network for specification and commitment of murine T cells

T cells develop from multipotent progenitors by a gradual process dependent on intrathymic Notch signaling and coupled with extensive proliferation. The stages leading them to T-cell lineage commitment are well characterized by single-cell and bulk RNA analyses of sorted populations and by direct measurements of precursor-product relationships. This process depends not only on Notch signaling but also on multiple transcription factors, some associated with stemness and multipotency, some with alternative lineages, and others associated with T-cell fate. These factors interact in opposing or semi-independent T cell gene regulatory network (GRN) subcircuits that are increasingly well defined. A newly comprehensive picture of this network has emerged. Importantly, because key factors in the GRN can bind to markedly different genomic sites at one stage than they do at other stages, the genes they significantly regulate are also stage-specific. Global transcriptome analyses of perturbations have revealed an underlying modular structure to the T-cell commitment GRN, separating decisions to lose "stem-ness" from decisions to block alternative fates. Finally, the updated network sheds light on the intimate relationship between the T-cell program, which depends on the thymus, and the innate lymphoid cell (ILC) program, which does not.

Multimodal human thymic profiling reveals trajectories and cellular milieu for T agonist selection

To prevent autoimmunity, thymocytes expressing self-reactive T cell receptors (TCRs) are negatively selected, however, divergence into tolerogenic, agonist selected lineages represent an alternative fate. As thymocyte development, selection, and lineage choices are dependent on spatial context and cell-to-cell interactions, we have performed Cellular Indexing of Transcriptomes and Epitopes by sequencing (CITE-seq) and spatial transcriptomics on paediatric human thymu​​s. Thymocytes expressing markers of strong TCR signalling diverged from the conventional developmental trajectory prior to CD4⁺ or CD8⁺ lineage commitment, while markers of different agonist selected T cell populations (CD8αα(I), CD8αα(II), T_(agonist), T_reg(diff), and T_reg) exhibited variable timing of induction. Expression profiles of chemokines and co-stimulatory molecules, together with spatial localisation, supported that dendritic cells, B cells, and stromal cells contribute to agonist selection, with different subsets influencing thymocytes at specific developmental stages within distinct spatial niches. Understanding factors influencing agonist T cells is needed to benefit from their immunoregulatory effects in clinical use.

Multi-omic analyses in immune cell development with lessons learned from T cell development

Traditionally, flow cytometry has been the preferred method to characterize immune cells at the single-cell level. Flow cytometry is used in immunology mostly to measure the expression of identifying markers on the cell surface, but-with good antibodies-can also be used to assess the expression of intracellular proteins. The advent of single-cell RNA-sequencing has paved the road to study immune development at an unprecedented resolution. Single-cell RNA-sequencing studies have not only allowed us to efficiently chart the make-up of heterogeneous tissues, including their most rare cell populations, it also increasingly contributes to our understanding how different omics modalities interplay at a single cell resolution. Particularly for investigating the immune system, this means that these single-cell techniques can be integrated to combine and correlate RNA and protein data at the single-cell level. While RNA data usually reveals a large heterogeneity of a given population identified solely by a combination of surface protein markers, the integration of different omics modalities at a single cell resolution is expected to greatly contribute to our understanding of the immune system.

Genetic Strategies to Study T Cell Development

Genetics approaches have been instrumental to deciphering T cell development in the thymus, including gene disruption by homologous recombination and more recently Crispr-based gene editing and transgenic gene expression, especially of specific T cell antigen receptors (TCR). This brief chapter describes commonly used tools and strategies to modify the genome of thymocytes, including mouse strains with lineage- and stage-specific expression of the Cre recombinase used for conditional allele inactivation or expressing unique antigen receptor specificities.

From pseudo to real-time dynamics of T cell thymic differentiation

Numerous methods have recently emerged for ordering single cells along developmental trajectories. However, accurate depiction of developmental dynamics can only be achieved after rescaling the trajectory according to the relative time spent at each developmental point. We formulate a model which estimates local cell densities and fluxes, and incorporates cell division and apoptosis rates, to infer the real-time dimension of the developmental trajectory. We validate the model using mathematical simulations and apply it to experimental high dimensional cytometry data obtained from the mouse thymus to construct the true time profile of the thymocyte developmental process. Our method can easily be implemented in any of the existing tools for trajectory inference.

Single-cell immune profiling reveals thymus-seeding populations, T cell commitment, and multilineage development in the human thymus

T cell development in the mouse thymus has been studied extensively, but less is known regarding T cell development in the human thymus. We used a combination of single-cell techniques and functional assays to perform deep immune profiling of human T cell development, focusing on the initial stages of prelineage commitment. We identified three thymus-seeding progenitor populations that also have counterparts in the bone marrow. In addition, we found that the human thymus physiologically supports the development of monocytes, dendritic cells, and NK cells, as well as limited development of B cells. These results are an important step toward monitoring and guiding regenerative therapies in patients after hematopoietic stem cell transplantation.

Transcriptional dynamics and epigenetic regulation of E and ID protein encoding genes during human T cell development

T cells are generated from hematopoietic stem cells through a highly organized developmental process, in which stage-specific molecular events drive maturation towards αβ and γδ T cells. Although many of the mechanisms that control αβ- and γδ-lineage differentiation are shared between human and mouse, important differences have also been observed. Here, we studied the regulatory dynamics of the E and ID protein encoding genes during pediatric human T cell development by evaluating changes in chromatin accessibility, histone modifications and bulk and single cell gene expression. We profiled patterns of ID/E protein activity and identified up- and downstream regulators and targets, respectively. In addition, we compared transcription of E and ID protein encoding genes in human versus mouse to predict both shared and unique activities in these species, and in prenatal versus pediatric human T cell differentiation to identify regulatory changes during development. This analysis showed a putative involvement of TCF3/E2A in the development of γδ T cells. In contrast, in αβ T cell precursors a pivotal pre-TCR-driven population with high ID gene expression and low predicted E protein activity was identified. Finally, in prenatal but not postnatal thymocytes, high HEB/TCF12 levels were found to counteract high ID levels to sustain thymic development. In summary, we uncovered novel insights in the regulation of E and ID proteins on a cross-species and cross-developmental level.

Function of normal oral mucosa revealed by single-cell RNA sequencing

The functions of oral mucosa include barrier, sensation, and secretion. The barrier protection function is particularly important, which includes physical barrier and immunological barrier. Few studies have revealed the function of oral mucosa by displaying the map of normal oral mucosal cells from the perspective of single cells. Here, single-cell transcriptome sequencing was used to bring a relatively comprehensive map of the normal oral mucosal cells. In total, 26,398 cells from three cases of normal oral mucosa were analyzed by single-cell RNA-sequencing and 14 distinct cell groups were defined, 7 of which were immune cells. We performed subgroup classification and heterogeneity analysis of epithelial cells, T cells, and macrophagocytes, which found a subpopulation of epithelial cells with high expression of major histocompatibility complex class II molecules, a subpopulation CD8⁺ GZMK⁺ T cells, and two kinds of active macrophagocytes. Meanwhile, we identified ligand-receptor pairs among the major cell types to explore the interactions and how they maintain the homeostasis of normal oral mucosa. Based on these results, the epithelial barrier function, immunological barrier function, and potential maintenance function of stromal cells in the oral mucosa were described at the single-cell level, which provides basic data resources for further studies of oral mucosal diseases.

A single-cell analysis of thymopoiesis and thymic iNKT cell development in pigs

Many aspects of the porcine immune system remain poorly characterized, which poses a barrier to improving swine health and utilizing pigs as preclinical models. Here, we employ single-cell RNA sequencing (scRNA-seq) to create a cell atlas of the early-adolescent pig thymus. Our data show conserved features as well as species-specific differences in cell states and cell types compared with human thymocytes. We also describe several unconventional T cell types with gene expression profiles associated with innate effector functions. This includes a cell census of more than 11,000 differentiating invariant natural killer T (iNKT) cells, which reveals that the functional diversity of pig iNKT cells differs substantially from the iNKT0/1/2/17 subset differentiation paradigm established in mice. Our data characterize key differentiation events in porcine thymopoiesis and iNKT cell maturation and provide important insights into pig T cell development.

Single-cell deletion analyses show control of pro-T cell developmental speed and pathways by Tcf7, Spi1, Gata3, Bcl11a, Erg, and Bcl11b

As early T cell precursors transition from multipotentiality to T lineage commitment, they change expression of multiple transcription factors. It is unclear whether individual transcription factors directly control choices between T cell identity and some alternative fate or whether these factors mostly affect proliferation or survival during the normal commitment process. Here, we unraveled the impacts of deleting individual transcription factors at two stages in early T cell development, using synchronized in vitro differentiation systems, single-cell RNA-seq with batch indexing, and controlled gene-disruption strategies. First, using a customized method for single-cell CRISPR disruption, we defined how the early-acting transcription factors Bcl11a, Erg, Spi1 (PU.1), Gata3, and Tcf7 (TCF1) function before commitment. The results revealed a kinetic tug of war within individual cells between T cell factors Tcf7 and Gata3 and progenitor factors Spi1 and Bcl11a, with an unexpected guidance role for Erg. Second, we tested how activation of transcription factor Bcl11b during commitment altered ongoing cellular programs. In knockout cells where Bcl11b expression was prevented, the cells did not undergo developmental arrest, instead following an alternative path as T lineage commitment was blocked. A stepwise, time-dependent regulatory cascade began with immediate-early transcription factor activation and E protein inhibition, finally leading Bcl11b knockout cells toward exit from the T cell pathway. Last, gene regulatory networks of transcription factor cross-regulation were extracted from the single-cell transcriptome results, characterizing the specification network operating before T lineage commitment and revealing its links to both the Bcl11b knockout alternative network and the network consolidating T cell identity during commitment.

acorde unravels functionally interpretable networks of isoform co-usage from single cell data

Alternative splicing (AS) is a highly-regulated post-transcriptional mechanism known to modulate isoform expression within genes and contribute to cell-type identity. However, the extent to which alternative isoforms establish co-expression networks that may be relevant in cellular function has not been explored yet. Here, we present acorde, a pipeline that successfully leverages bulk long reads and single-cell data to confidently detect alternative isoform co-expression relationships. To achieve this, we develop and validate percentile correlations, an innovative approach that overcomes data sparsity and yields accurate co-expression estimates from single-cell data. Next, acorde uses correlations to cluster co-expressed isoforms into a network, unraveling cell type-specific alternative isoform usage patterns. By selecting same-gene isoforms between these clusters, we subsequently detect and characterize genes with co-differential isoform usage (coDIU) across cell types. Finally, we predict functional elements from long read-defined isoforms and provide insight into biological processes, motifs, and domains potentially controlled by the coordination of post-transcriptional regulation. The code for acorde is available at https://github.com/ConesaLab/acorde .

The transcription factor LRF promotes integrin β7 expression by and gut homing of CD8αα+ intraepithelial lymphocyte precursors

While T cell receptor (TCR) αβ+CD8α+CD8β- intraepithelial lymphocytes (CD8αα+ IELs) differentiate from thymic IEL precursors (IELps) and contribute to gut homeostasis, the transcriptional control of their development remains poorly understood. In the present study we showed that mouse thymocytes deficient for the transcription factor leukemia/lymphoma-related factor (LRF) failed to generate TCRαβ+CD8αα+ IELs and their CD8β-expressing counterparts, despite giving rise to thymus and spleen CD8αβ+ T cells. LRF-deficient IELps failed to migrate to the intestine and to protect against T cell-induced colitis, and had impaired expression of the gut-homing integrin α4β7. Single-cell RNA-sequencing found that LRF was necessary for the expression of genes characteristic of the most mature IELps, including Itgb7, encoding the β7 subunit of α4β7. Chromatin immunoprecipitation and gene-regulatory network analyses both defined Itgb7 as an LRF target. Our study identifies LRF as an essential transcriptional regulator of IELp maturation in the thymus and subsequent migration to the intestinal epithelium.

Patterns of Differentially Expressed circRNAs in Human Thymocytes

Circular RNAs (circRNAs) are suggested to play a discriminative role between some stages of thymocyte differentiation. However, differential aspects of the stage of mature single-positive thymocytes remain to be explored. The purpose of this study is to investigate the differential expression pattern of circRNAs in three different development stages of human thymocytes, including mature single-positive cells, and perform predictions in silico regarding the ability of specific circRNAs when controlling the expression of genes involved in thymocyte differentiation. We isolate human thymocytes at three different stages of intrathymic differentiation and determine the expression of circRNAs and mRNA by RNASeq. We show that the differential expression pattern of 50 specific circRNAs serves to discriminate between the three human thymocyte populations. Interestingly, the downregulation of RAG2, a gene involved in T-cell differentiation in the thymus, could be simultaneously controlled by the downregulation of two circRNASs (hsa_circ_0031584 and hsa_circ_0019079) through the hypothetical liberation of hsa-miR-609. Our study provides, for the first time, significant insights into the usefulness of circRNAs in discriminating between different stages of thymocyte differentiation and provides new potential circRNA–miRNA–mRNA networks capable of controlling the expression of genes involved in T-cell differentiation in the thymus.

Multi-objective optimization reveals time- and dose-dependent inflammatory cytokine-mediated regulation of human stem cell derived T-cell development

The generation of T-cells from stem cells in vitro could provide an alternative source of cells for immunotherapies. T-cell development from hematopoietic stem and progenitor cells (HSPCs) is tightly regulated through Notch pathway activation by Delta-like (DL) ligands 1 and 4. Other molecules, such as stem cell factor (SCF) and interleukin (IL)-7, play a supportive role in regulating the survival, differentiation, and proliferation of developing T-cells. Numerous other signaling molecules influence T-lineage development in vivo, but little work has been done to understand and optimize their use for T-cell production. Using a defined engineered thymic niche system, we undertook a multi-stage statistical learning-based optimization campaign and identified IL-3 and tumor necrosis factor α (TNFα) as a stage- and dose-specific enhancers of cell proliferation and T-lineage differentiation. We used this information to construct an efficient three-stage process for generating conventional TCRαβ⁺CD8⁺ T-cells expressing a diverse TCR repertoire from blood stem cells. Our work provides new insight into T-cell development and a robust system for generating T-cells to enable clinical therapies for treating cancer and immune disorders.

Natural Barcodes for Longitudinal Single Cell Tracking of Leukemic and Immune Cell Dynamics

Blood malignancies provide unique opportunities for longitudinal tracking of disease evolution following therapeutic bottlenecks and for the monitoring of changes in anti-tumor immunity. The expanding development of multi-modal single-cell sequencing technologies affords newer platforms to elucidate the mechanisms underlying these processes at unprecedented resolution. Furthermore, the identification of molecular events that can serve as in-vivo barcodes now facilitate the tracking of the trajectories of malignant and of immune cell populations over time within primary human samples, as these permit unambiguous identification of the clonal lineage of cell populations within heterogeneous phenotypes. Here, we provide an overview of the potential for chromosomal copy number changes, somatic nuclear and mitochondrial DNA mutations, single nucleotide polymorphisms, and T and B cell receptor sequences to serve as personal natural barcodes and review technical implementations in single-cell analysis workflows. Applications of these methodologies include the study of acquired therapeutic resistance and the dissection of donor- and host cellular interactions in the context of allogeneic hematopoietic stem cell transplantation.

Azvudine is a thymus-homing anti-SARS-CoV-2 drug effective in treating COVID-19 patients

Azvudine (FNC) is a nucleoside analog that inhibits HIV-1 RNA-dependent RNA polymerase (RdRp). Recently, we discovered FNC an agent against SARS-CoV-2, and have taken it into Phase III trial for COVID-19 patients. FNC monophosphate analog inhibited SARS-CoV-2 and HCoV-OC43 coronavirus with an EC₅₀ between 1.2 and 4.3 μM, depending on viruses or cells, and selective index (SI) in 15-83 range. Oral administration of FNC in rats revealed a substantial thymus-homing feature, with FNC triphosphate (the active form) concentrated in the thymus and peripheral blood mononuclear cells (PBMC). Treating SARS-CoV-2 infected rhesus macaques with FNC (0.07 mg/kg, qd, orally) reduced viral load, recuperated the thymus, improved lymphocyte profiles, alleviated inflammation and organ damage, and lessened ground-glass opacities in chest X-ray. Single-cell sequencing suggested the promotion of thymus function by FNC. A randomized, single-arm clinical trial of FNC on compassionate use (n = 31) showed that oral FNC (5 mg, qd) cured all COVID-19 patients, with 100% viral ribonucleic acid negative conversion in 3.29 ± 2.22 days (range: 1-9 days) and 100% hospital discharge rate in 9.00 ± 4.93 days (range: 2-25 days). The side-effect of FNC is minor and transient dizziness and nausea in 16.12% (5/31) patients. Thus, FNC might cure COVID-19 through its anti-SARS-CoV-2 activity concentrated in the thymus, followed by promoted immunity.

Recent advances in T-cell lymphoid neoplasms

T Cells comprise many subtypes of specified lymphocytes, and their differentiation and function take place in different tissues. This cellular diversity is also observed in the multiple ways T-cell transformation gives rise to a variety of T-cell neoplasms. This review covers the main types of T-cell malignancies and their specific characteristics, emphasizing recent advances at the cellular and molecular levels as well as differences and commonalities among them.

Delineating spatiotemporal and hierarchical development of human fetal innate lymphoid cells

Whereas the critical roles of innate lymphoid cells (ILCs) in adult are increasingly appreciated, their developmental hierarchy in early human fetus remains largely elusive. In this study, we sorted human hematopoietic stem/progenitor cells, lymphoid progenitors, putative ILC progenitor/precursors and mature ILCs in the fetal hematopoietic, lymphoid and non-lymphoid tissues, from 8 to 12 post-conception weeks, for single-cell RNA-sequencing, followed by computational analysis and functional validation at bulk and single-cell levels. We delineated the early phase of ILC lineage commitment from hematopoietic stem/progenitor cells, which mainly occurred in fetal liver and intestine. We further unveiled interleukin-3 receptor as a surface marker for the lymphoid progenitors in fetal liver with T, B, ILC and myeloid potentials, while IL-3RA- lymphoid progenitors were predominantly B-lineage committed. Notably, we determined the heterogeneity and tissue distribution of each ILC subpopulation, revealing the proliferating characteristics shared by the precursors of each ILC subtype. Additionally, a novel unconventional ILC2 subpopulation (CRTH2- CCR9+ ILC2) was identified in fetal thymus. Taken together, our study illuminates the precise cellular and molecular features underlying the stepwise formation of human fetal ILC hierarchy with remarkable spatiotemporal heterogeneity.

The Metabolic Landscape of Thymic T Cell Development In Vivo and In Vitro

Although metabolic pathways have been shown to control differentiation and activation in peripheral T cells, metabolic studies on thymic T cell development are still lacking, especially in human tissue. In this study, we use transcriptomics and extracellular flux analyses to investigate the metabolic profiles of primary thymic and in vitro-derived mouse and human thymocytes. Core metabolic pathways, specifically glycolysis and oxidative phosphorylation, undergo dramatic changes between the double-negative (DN), double-positive (DP), and mature single-positive (SP) stages in murine and human thymus. Remarkably, despite the absence of the complex multicellular thymic microenvironment, in vitro murine and human T cell development recapitulated the coordinated decrease in glycolytic and oxidative phosphorylation activity between the DN and DP stages seen in primary thymus. Moreover, by inducing in vitro T cell differentiation from Rag1^-/- mouse bone marrow, we show that reduced metabolic activity at the DP stage is independent of TCR rearrangement. Thus, our findings suggest that highly conserved metabolic transitions are critical for thymic T cell development.

Single-Cell RNA Sequencing Approaches for Tracing T Cell Development

T cell development occurs in the thymus, where uncommitted progenitors are directed into a range of sublineages with distinct functions. The goal is to generate a TCR repertoire diverse enough to recognize potential pathogens while remaining tolerant of self. Decades of intensive research have characterized the transcriptional programs controlling critical differentiation checkpoints at the population level. However, greater precision regarding how and when these programs orchestrate differentiation at the single-cell level is required. Single-cell RNA sequencing approaches are now being brought to bear on this question, to track the identity of cells and analyze their gene expression programs at a resolution not previously possible. In this review, we discuss recent advances in the application of these technologies that have the potential to yield unprecedented insight to T cell development.

Integrative Single-Cell Transcriptomic Analysis of Human Fetal Thymocyte Development

Intrathymic differentiation of T lymphocytes begins as early as intrauterine stage, yet the T cell lineage decisions of human fetal thymocytes at different gestational ages are not currently understood. Here, we performed integrative single-cell analyses of thymocytes across gestational ages. We identified conserved candidates underlying the selection of T cell receptor (TCR) lineages in different human fetal stages. The trajectory of early thymocyte commitment during fetal growth was also characterized. Comparisons with mouse data revealed conserved and species-specific transcriptional dynamics of thymocyte proliferation, apoptosis and selection. Genome-wide association study (GWAS) data associated with multiple autoimmune disorders were analyzed to characterize susceptibility genes that are highly expressed at specific stages during fetal thymocyte development. In summary, our integrative map describes previously underappreciated aspects of human thymocyte development, and provides a comprehensive reference for understanding T cell lymphopoiesis in a self-tolerant and functional adaptive immune system.

Single-cell technologies and analyses in hematopoiesis and hematological malignancies

In recent years, single-cell technologies have emerged as breakthrough techniques that enable the characterization of hematopoietic cell populations of normal and malignant tissue samples and will be combined in the near future with bulk technologies, currently used in clinical practice, to improve diagnosis, prognosis, and the search for novel molecular targets. These single-cell methods have the advantage of not masking cell-to-cell variation features and involve the study of genetic, epigenetic, transcriptional, and proteomic landscapes from a single-cell perspective. Latest advances in this field have enabled the development of novel strategies that significantly increase both sensitivity and high throughput. In this review, we emphasize emerging techniques aimed at assessing individual or multiomic parameters at single-cell resolution and analyze how these technologies have helped us understand hematopoietic variability and identify unknown and/or rare subpopulations. We also summarize the impact of these single-cell profiling strategies on the characterization of cell diversity within the tumor and the clonal evolution of multiple hematological malignancies in samples from untreated and treated patients, which provide valuable information for diagnosis, prognosis, and future treatments and explain why current therapies may fail. However, despite these improvements, new challenges lie ahead.

Epigenetic Dynamics in the Function of T-Lineage Regulatory Factor Bcl11b

The transcription factor Bcl11b is critically required to support the development of diverse cell types, including T lymphocytes, type 2 innate lymphoid cells, neurons, craniofacial mesenchyme and keratinocytes. Although in T cell development its onset of expression is tightly linked to T-lymphoid lineage commitment, the Bcl11b protein in fact regulates substantially different sets of genes in different lymphocyte populations, playing strongly context-dependent roles. Somewhat unusually for lineage-defining transcription factors with site-specific DNA binding activity, much of the reported chromatin binding of Bcl11b appears to be indirect, or guided in large part by interactions with other transcription factors. We describe evidence suggesting that a further way in which Bcl11b exerts such distinct stage-dependent functions is by nucleating changes in regional suites of epigenetic modifications through recruitment of multiple families of chromatin-modifying enzyme complexes. Herein we explore what is - and what remains to be - understood of the roles of Bcl11b, its cofactors, and how it modifies the epigenetic state of the cell to enforce its diverse set of context-specific transcriptional and developmental programs.

Development of double-positive thymocytes at single-cell resolution

BACKGROUND

T cells generated from thymopoiesis are essential for the immune system, and recent single-cell studies have contributed to our understanding of the development of thymocytes at the genetic and epigenetic levels. However, the development of double-positive (DP) T cells, which comprise the majority of thymocytes, has not been well investigated.

METHODS

We applied single-cell sequencing to mouse thymocytes and analyzed the transcriptome data using Seurat. By applying unsupervised clustering, we defined thymocyte subtypes and validated DP cell subtypes by flow cytometry. We classified the cell cycle phases of each cell according to expression of cell cycle phase-specific genes. For immune synapse detection, we used immunofluorescent staining and ImageStream-based flow cytometry. We studied and integrated human thymocyte data to verify the conservation of our findings and also performed cross-species comparisons to examine species-specific gene regulation.

RESULTS

We classified blast, rearrangement, and selection subtypes of DP thymocytes and used the surface markers CD2 and Ly6d to identify these subtypes by flow cytometry. Based on this new classification, we found that the proliferation of blast DP cells is quite different from that of double-positive cells and other cell types, which tend to exit the cell cycle after a single round. At the DP cell selection stage, we observed that CD8-associated immune synapses formed between thymocytes, indicating that CD8sp selection occurred among thymocytes themselves. Moreover, cross-species comparison revealed species-specific transcription factors (TFs) that contribute to the transcriptional differences of thymocytes from humans and mice.

CONCLUSIONS

Our study classified DP thymocyte subtypes of different developmental stages and provided new insight into the development of DP thymocytes at single-cell resolution, furthering our knowledge of the fundamental immunological process of thymopoiesis.

Single-cell insights into the hematopoietic generation of T-lymphocyte precursors in mouse and human

T-Cell development is a major branch of lymphoid development and a key output of hematopoiesis, especially in early life, but the molecular requirements for T-cell potential have remained obscure. Considerable advances have now been made toward solving this problem through single-cell transcriptome studies, interfaced with in vitro differentiation assays that monitor potential efficiently at the single-cell level. This review focuses on a series of recent reports studying mouse and human early T-cell precursors, both in the developing fetus and in stringently purified postnatal samples of intrathymic and prethymic T-lineage precursors. Cross-comparison of results reveals a robustly conserved core program in mouse and human, but with some informative and provocative variations between species and between ontogenic states. Repeated findings are the multipotent progenitor regulatory signature of thymus-seeding cells and the proximity of the T-cell program to dendritic cell programs, especially to plasmacytoid dendritic cells in humans.

How transcription factors drive choice of the T cell fate

Recent evidence has elucidated how multipotent blood progenitors transform their identities in the thymus and undergo commitment to become T cells. Together with environmental signals, a core group of transcription factors have essential roles in this process by directly activating and repressing specific genes. Many of these transcription factors also function in later T cell development, but control different genes. Here, we review how these transcription factors work to change the activities of specific genomic loci during early intrathymic development to establish T cell lineage identity. We introduce the key regulators and highlight newly emergent insights into the rules that govern their actions. Whole-genome deep sequencing-based analysis has revealed unexpectedly rich relationships between inherited epigenetic states, transcription factor-DNA binding affinity thresholds and influences of given transcription factors on the activities of other factors in the same cells. Together, these mechanisms determine T cell identity and make the lineage choice irreversible.

T Cell Development: Old Tales Retold By Single-Cell RNA Sequencing

Mammalian T cell development initiates from the migration of hematopoietic progenitors to the thymus, which undergo cell proliferation, T-lineage specification and commitment, as well as positive and negative selection. These processes are precisely controlled at multiple levels and have been intensively studied using gene-modified animal models and in vitro coculture systems. However, several long-standing questions, including the characterization of the rare but crucial progenitors/precursors and the molecular mechanisms underlying their fate decision, have been dampened because of cell scarcity and lack of appropriate techniques. Single-cell RNA sequencing (scRNA-seq) makes it possible to investigate and resolve some of these questions, leading to new remarkable progress in identifying and characterizing early thymic progenitors and delineating the refined developmental trajectories of conventional and unconventional T cells.

An Integrated Epigenomic and Transcriptomic Map of Mouse and Human αβ T Cell Development

αβ lineage T cells, most of which are CD4+ or CD8+ and recognize MHC I- or MHC II-presented antigens, are essential for immune responses and develop from CD4+CD8+ thymocytes. The absence of in vitro models and the heterogeneity of αβ thymocytes have hampered analyses of their intrathymic differentiation. Here, combining single-cell RNA and ATAC (chromatin accessibility) sequencing, we identified mouse and human αβ thymocyte developmental trajectories. We demonstrated asymmetric emergence of CD4+ and CD8+ lineages, matched differentiation programs of agonist-signaled cells to their MHC specificity, and identified correspondences between mouse and human transcriptomic and epigenomic patterns. Through computational analysis of single-cell data and binding sites for the CD4+-lineage transcription factor Thpok, we inferred transcriptional networks associated with CD4+- or CD8+-lineage differentiation, and with expression of Thpok or of the CD8+-lineage factor Runx3. Our findings provide insight into the mechanisms of CD4+ and CD8+ T cell differentiation and a foundation for mechanistic investigations of αβ T cell development.

Processing Human Thymic Tissue for Single Cell RNA-Seq

Single cell RNA sequencing of human thymic cells is dependent on isolation of highly pure and viable cell populations. This protocol describes the isolation of CD34⁺ progenitor and more differentiated CD34^– fractions from post-natal thymic tissue to study thymopoiesis. CD34⁺ cells represent <1% of thymic cells, so this protocol uses magnetic- followed by fluorescence-activated cell separation to isolate highly enriched CD34⁺ cells.

For complete details on the use and execution of this protocol, please refer to Le et al. (2020).

•

Protocol for processing of human thymus for single cell RNA-seq

•

Thymus dissection and density gradient centrifugation isolate mononuclear cells

•

Magnetic-activated cell separation isolates CD34⁺/CD34^– cells

•

Fluorescence-activated cell separation isolates populations for RNA-seq