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

Comparative transcriptomic analyses of thymocytes using 10x Genomics and Parse scRNA-seq technologies

BACKGROUND

Single-cell RNA sequencing experiments commonly use 10x Genomics (10x) kits due to their high-throughput capacity and standardized protocols. Recently, Parse Biosciences (Parse) introduced an alternative technology that uses multiple in-situ barcoding rounds within standard 96-well plates. Parse enables the analysis of more cells from multiple samples in a single run without the need for additional reagents or specialized microfluidics equipment. To evaluate the performance of both platforms, we conducted a benchmark study using biological and technical replicates of mouse thymus as a complex immune tissue.

RESULTS

We found that Parse detected nearly twice the number of genes compared to 10x, with each platform detecting a distinct set of genes. The comparison of multiplexed samples generated from 10x and Parse techniques showed 10x data to have lower technical variability and more precise annotation of biological states in the thymus compared to Parse.

CONCLUSION

Our results provide a comprehensive comparison of the suitability of both single-cell platforms for immunological studies.

Transcriptional network dynamics in early T cell development

The rate at which cells enter the T cell pathway depends not only on the immigration of hematopoietic precursors into the strong Notch signaling environment of the thymus but also on the kinetics with which each individual precursor cell reaches T-lineage commitment once it arrives. Notch triggers a complex, multistep gene regulatory network in the cells in which the steps are stereotyped but the transition speeds between steps are variable. Progenitor-associated transcription factors delay T-lineage differentiation even while Notch-induced transcription factors within the same cells push differentiation forward. Progress depends on regulator cross-repression, on breaching chromatin barriers, and on shifting, competitive collaborations between stage-specific and stably expressed transcription factors, as reviewed here.

Rediscovering the human thymus through cutting-edge technologies

Recent technological advances have transformed our understanding of the human thymus. Innovations such as high-resolution imaging, single-cell omics, and organoid cultures, including thymic epithelial cell (TEC) differentiation and culture, and improvements in biomaterials, have further elucidated the thymus architecture, cellular dynamics, and molecular mechanisms underlying T cell development, and have unraveled previously unrecognized levels of stromal cell heterogeneity. These advancements offer unprecedented insights into thymic biology and hold promise for the development of novel therapeutic strategies for immune-related disorders.

Single-cell sequencing has revealed a more complex array of thymic epithelial cells

Thymic epithelial cells participate in the maturation and selection of T lymphocytes. This review explores recent insights from single-cell sequencing regarding classifying thymic epithelial cells in both normal and neoplastic thymus. Cortical thymic epithelial cells facilitate thymocyte differentiation and contribute to positive selection. Medullary epithelial cells are distinguished by their expression of AIRE. Cells progress from a pre-AIRE state, containing precursors with cortical and medullary characteristics, termed junctional cells. Mature medullary epithelial cells exhibit promiscuous gene expression and after that downregulate AIRE mRNA. Post-AIRE cells can adopt a Hassall corpuscle-like phenotype or exhibit distinctive differentiation characteristics including tuft cells, ionocytes, neuroendocrine cells, and myoid cells.

Genomic analysis and functional properties of Lactobacillus johnsonii GJ231 isolated from healthy beagles

Probiotics are one of the management tools to improve the host's healthy microbiota. The positive effects of probiotics on host health are species-specific, so probiotics isolated from host's own gut may be most beneficial. Many of the metabolites (e.g., short-chain fatty acids, bacteriocins, and hydrogen peroxide) produced by Lactobacillus johnsonii have specific inhibitory profiles against invading pathogens. In this study, we isolated L. johnsonii GJ231 from the intestinal tract of healthy female beagles. The genome size of 1.763 M encoded a total of 1,691 predicted genes. Many carbohydrate-active enzymes responsible for carbohydrate degradation and the production of short-chain fatty acids were also predicted. The metabolic profile of short-chain fatty acids in L. johnsonii GJ231 was determined using LC-MS/MS. The bacteriocin-producing gene bacteriocin (lactacin F) in L. johnsonii GJ231 was also predicted. In vitro, experiments demonstrated that GJ231 can thrive in weak acids, 0.3% bile salts, and artificial gastrointestinal fluid models. It was tolerant of to high temperatures up to 70°C, was non- hemolytic, inhibited pathogenic bacteria, and had a high antioxidant capacity. In vivo safety experiments conducted in mice revealed that oral administration of GJ231 not only had no toxic side effect but also increased their antioxidant capacity. In conclusion, combining the above test results, which collectively demonstrate that canine-derived L. johnsonii GJ231 was a non-pathogenic, acid-tolerant and bile-salt-tolerant probiotic strain that inhibits pathogenic bacteria and improves host antioxidant function. This may make it a promising candidate for the development of innovative functional foods for pets.

Unraveling the spatial organization and development of human thymocytes through integration of spatial transcriptomics and single-cell multi-omics profiling

The structural components of the thymus are essential for guiding T cell development, but a thorough spatial view is still absent. Here we develop the TSO-his tool, designed to integrate multimodal data from single-cell and spatial transcriptomics to decipher the intricate structure of human thymus. Specifically, we characterize dynamic changes in cell types and critical markers, identifying ELOVL4 as a mediator of CD4+ T cell positive selection in the cortex. Utilizing the mapping function of TSO-his, we reconstruct thymic spatial architecture at single-cell resolution and recapitulates classical cell types and their essential co-localization for T cell development; additionally, previously unknown co-localization relationships such as that of CD8αα with memory B cells and monocytes are identified. Incorporating VDJ sequencing data, we also delineate distinct intermediate thymocyte states during αβ T cell development. Overall, these insights enhance our understanding of thymic biology and may inform therapeutic interventions targeting T cell-mediated immune responses.

Thymic epithelial organoids mediate T-cell development

Although the advent of organoids has opened unprecedented perspectives for basic and translational research, immune system-related organoids remain largely underdeveloped. Here, we established organoids from the thymus, the lymphoid organ responsible for T-cell development. We identified conditions enabling mouse thymic epithelial progenitor cell proliferation and development into organoids with diverse cell populations and transcriptional profiles resembling in vivo thymic epithelial cells (TECs) more closely than traditional TEC cultures. In contrast to these two-dimensional cultures, thymic epithelial organoids maintained thymus functionality in vitro and mediated physiological T-cell development upon reaggregation with T-cell progenitors. The reaggregates showed in vivo-like epithelial diversity and the ability to attract T-cell progenitors. Thymic epithelial organoids are the first organoids originating from the stromal compartment of a lymphoid organ. They provide new opportunities to study TEC biology and T-cell development in vitro, paving the way for future thymic regeneration strategies in ageing or acute injuries.

Age-related epithelial defects limit thymic function and regeneration

The thymus is essential for establishing adaptive immunity yet undergoes age-related involution that leads to compromised immune responsiveness. The thymus is also extremely sensitive to acute insult and although capable of regeneration, this capacity declines with age for unknown reasons. We applied single-cell and spatial transcriptomics, lineage-tracing and advanced imaging to define age-related changes in nonhematopoietic stromal cells and discovered the emergence of two atypical thymic epithelial cell (TEC) states. These age-associated TECs (aaTECs) formed high-density peri-medullary epithelial clusters that were devoid of thymocytes; an accretion of nonproductive thymic tissue that worsened with age, exhibited features of epithelial-to-mesenchymal transition and was associated with downregulation of FOXN1. Interaction analysis revealed that the emergence of aaTECs drew tonic signals from other functional TEC populations at baseline acting as a sink for TEC growth factors. Following acute injury, aaTECs expanded substantially, further perturbing trophic regeneration pathways and correlating with defective repair of the involuted thymus. These findings therefore define a unique feature of thymic involution linked to immune aging and could have implications for developing immune-boosting therapies in older individuals.

Direct presentation of inflammation-associated self-antigens by thymic innate-like T cells induces elimination of autoreactive CD8+ thymocytes

Upregulation of diverse self-antigens that constitute components of the inflammatory response overlaps spatially and temporally with the emergence of pathogen-derived foreign antigens. Therefore, discrimination between these inflammation-associated self-antigens and pathogen-derived molecules represents a unique challenge for the adaptive immune system. Here, we demonstrate that CD8+ T cell tolerance to T cell-derived inflammation-associated self-antigens is efficiently induced in the thymus and supported by redundancy in cell types expressing these molecules. In addition to thymic epithelial cells, this included thymic eosinophils and innate-like T cells, a population that expressed molecules characteristic for all major activated T cell subsets. We show that direct T cell-to-T cell antigen presentation by minute numbers of innate-like T cells was sufficient to eliminate autoreactive CD8+ thymocytes. Tolerance to such effector molecules was of critical importance, as its breach caused by decreased thymic abundance of a single model inflammation-associated self-antigen resulted in autoimmune elimination of an entire class of effector T cells.

Identifying cell states in single-cell RNA-seq data at statistically maximal resolution

Single-cell RNA sequencing (scRNA-seq) has become a popular experimental method to study variation of gene expression within a population of cells. However, obtaining an accurate picture of the diversity of distinct gene expression states that are present in a given dataset is highly challenging because of the sparsity of the scRNA-seq data and its inhomogeneous measurement noise properties. Although a vast number of different methods is applied in the literature for clustering cells into subsets with 'similar' expression profiles, these methods generally lack rigorously specified objectives, involve multiple complex layers of normalization, filtering, feature selection, dimensionality-reduction, employ ad hoc measures of distance or similarity between cells, often ignore the known measurement noise properties of scRNA-seq measurements, and include a large number of tunable parameters. Consequently, it is virtually impossible to assign concrete biophysical meaning to the clusterings that result from these methods. Here we address the following problem: Given raw unique molecule identifier (UMI) counts of an scRNA-seq dataset, partition the cells into subsets such that the gene expression states of the cells in each subset are statistically indistinguishable, and each subset corresponds to a distinct gene expression state. That is, we aim to partition cells so as to maximally reduce the complexity of the dataset without removing any of its meaningful structure. We show that, given the known measurement noise structure of scRNA-seq data, this problem is mathematically well-defined and derive its unique solution from first principles. We have implemented this solution in a tool called Cellstates which operates directly on the raw data and automatically determines the optimal partition and cluster number, with zero tunable parameters. We show that, on synthetic datasets, Cellstates almost perfectly recovers optimal partitions. On real data, Cellstates robustly identifies subtle substructure within groups of cells that are traditionally annotated as a common cell type. Moreover, we show that the diversity of gene expression states that Cellstates identifies systematically depends on the tissue of origin and not on technical features of the experiments such as the total number of cells and total UMI count per cell. In addition to the Cellstates tool we also provide a small toolbox of software to place the identified cellstates into a hierarchical tree of higher-order clusters, to identify the most important differentially expressed genes at each branch of this hierarchy, and to visualize these results.

Gene Regulatory Programs that Specify Age-Related Differences during Thymocyte Development

T cell development is fundamental to immune system establishment, yet how this development changes with age remains poorly understood. Here, we construct a transcriptional and epigenetic atlas of T cell developmental programs in neonatal and adult mice, revealing the ontogeny of divergent gene regulatory programs and their link to age-related differences in phenotype and function. Specifically, we identify a gene module that diverges with age from the earliest stages of genesis and includes programs that govern effector response and cell cycle regulation. Moreover, we reveal that neonates possess more accessible chromatin during early thymocyte development, likely establishing poised gene expression programs that manifest later in thymocyte development. Finally, we leverage this atlas, employing a CRISPR-based perturbation approach coupled with single-cell RNA sequencing as a readout to uncover a conserved transcriptional regulator, Zbtb20, that contributes to age-dependent differences in T cell development. Altogether, our study defines transcriptional and epigenetic programs that regulate age-specific differences in T cell development.

Transcriptomic profiling of thymic dysregulation and viral tropism after neonatal roseolovirus infection

Introduction

Herpesviruses, including the roseoloviruses, have been linked to autoimmune disease. The ubiquitous and chronic nature of these infections have made it difficult to establish a causal relationship between acute infection and subsequent development of autoimmunity. We have shown that murine roseolovirus (MRV), which is highly related to human roseoloviruses, induces thymic atrophy and disruption of central tolerance after neonatal infection. Moreover, neonatal MRV infection results in development of autoimmunity in adult mice, long after resolution of acute infection. This suggests that MRV induces durable immune dysregulation.

Methods

In the current studies, we utilized single-cell RNA sequencing (scRNAseq) to study the tropism of MRV in the thymus and determine cellular processes in the thymus that were disrupted by neonatal MRV infection. We then utilized tropism data to establish a cell culture system.

Results

Herein, we describe how MRV alters the thymic transcriptome during acute neonatal infection. We found that MRV infection resulted in major shifts in inflammatory, differentiation and cell cycle pathways in the infected thymus. We also observed shifts in the relative number of specific cell populations. Moreover, utilizing expression of late viral transcripts as a proxy of viral replication, we identified the cellular tropism of MRV in the thymus. This approach demonstrated that double negative, double positive, and CD4 single positive thymocytes, as well as medullary thymic epithelial cells were infected by MRV in vivo. Finally, by applying pseudotime analysis to viral transcripts, which we refer to as "pseudokinetics," we identified viral gene transcription patterns associated with specific cell types and infection status. We utilized this information to establish the first cell culture systems susceptible to MRV infection in vitro.

Conclusion

Our research provides the first complete picture of roseolovirus tropism in the thymus after neonatal infection. Additionally, we identified major transcriptomic alterations in cell populations in the thymus during acute neonatal MRV infection. These studies offer important insight into the early events that occur after neonatal MRV infection that disrupt central tolerance and promote autoimmune disease.

Single-Cell Sequencing Illuminates Thymic Development: An Updated Framework for Understanding Thymic Epithelial Tumors

Thymic epithelial tumors (TETs) are rare tumors for which treatment options are limited. The ongoing need for improved systemic therapies reflects a limited understanding of tumor biology as well as the normal thymus. The essential role of the thymus in adaptive immunity is largely effected by its epithelial compartment, which directs thymocyte (T-cell) differentiation and immunologic self-tolerance. With aging, the thymus undergoes involution whereby epithelial tissue is replaced by adipose and other connective tissue, decreasing immature T-cell production. Against this natural drive toward involution, a fraction of thymuses will instead undergo oncologic transformation, leading to the formation of TETs, including thymoma and thymic carcinoma. The rarity of these tumors restricts investigation of the mechanisms of tumorigenesis and development of rational treatment options. To this end, the development of technologies which allow deep molecular profiling of individual tumor cells permits a new window through which to view normal thymic development and contrast the malignant changes that result in oncogenic transformation. In this review, we describe the findings of recent illuminating studies on the diversity of cell types within the epithelial compartment through thymic differentiation and aging. We contextualize these findings around important unanswered questions regarding the spectrum of known somatic tumor alterations, cell of origin, and tumor heterogeneity. The perspectives informed by single-cell molecular profiling offer new approaches to clinical and basic investigation of thymic epithelial tumors, with the potential to accelerate development of improved therapeutic strategies to address ongoing unmet needs in these rare tumors.

The Proteostasis of Thymic Stromal Cells in Health and Diseases

The thymus is the key immune organ for the development of T cells. Different populations of thymic stromal cells interact with T cells, thereby controlling the dynamic development of T cells through their differentiation and function. Proteostasis represents a balance between protein expression, folding, and modification and protein clearance, and its fluctuation usually depends at least partially on related protein regulatory systems for further survival and effects. However, in terms of the substantial requirement for self-antigens and their processing burden, increasing evidence highlights that protein regulation contributes to the physiological effects of thymic stromal cells. Impaired proteostasis may expedite the progression of thymic involution and dysfunction, accompanied by the development of autoimmune diseases or thymoma. Hence, in this review, we summarize the regulation of proteostasis within different types of thymic stromal cells under physiological and pathological conditions to identify potential targets for thymic regeneration and immunotherapy.

Derivation of functional thymic epithelial organoid lines from adult murine thymus

Thymic epithelial cells (TECs) orchestrate T cell development by imposing positive and negative selection on thymocytes. Current studies on TEC biology are hampered by the absence of long-term ex vivo culture platforms, while the cells driving TEC self-renewal remain to be identified. Here, we generate long-term (>2 years) expandable 3D TEC organoids from the adult mouse thymus. For further analysis, we generated single and double FoxN1-P2A-Clover, Aire-P2A-tdTomato, and Cldn4-P2A-tdTomato reporter lines by CRISPR knockin. Single-cell analyses of expanding clonal organoids reveal cells with bipotent stem/progenitor phenotypes. These clonal organoids can be induced to express Foxn1 and to generate functional cortical- and Aire-expressing medullary-like TECs upon RANK ligand + retinoic acid treatment. TEC organoids support T cell development from immature thymocytes in vitro as well as in vivo upon transplantation into athymic nude mice. This organoid-based platform allows in vitro study of TEC biology and offers a potential strategy for ex vivo T cell development.

Immune tolerance and the prevention of autoimmune diseases essentially depend on thymic tissue homeostasis

The intricate balance of immune reactions towards invading pathogens and immune tolerance towards self is pivotal in preventing autoimmune diseases, with the thymus playing a central role in establishing and maintaining this equilibrium. The induction of central immune tolerance in the thymus involves the elimination of self-reactive T cells, a mechanism essential for averting autoimmunity. Disruption of the thymic T cell selection mechanisms can lead to the development of autoimmune diseases. In the dynamic microenvironment of the thymus, T cell migration and interactions with thymic stromal cells are critical for the selection processes that ensure self-tolerance. Thymic epithelial cells are particularly significant in this context, presenting self-antigens and inducing the negative selection of autoreactive T cells. Further, the synergistic roles of thymic fibroblasts, B cells, and dendritic cells in antigen presentation, selection and the development of regulatory T cells are pivotal in maintaining immune responses tightly regulated. This review article collates these insights, offering a comprehensive examination of the multifaceted role of thymic tissue homeostasis in the establishment of immune tolerance and its implications in the prevention of autoimmune diseases. Additionally, the developmental pathways of the thymus are explored, highlighting how genetic aberrations can disrupt thymic architecture and function, leading to autoimmune conditions. The impact of infections on immune tolerance is another critical area, with pathogens potentially triggering autoimmunity by altering thymic homeostasis. Overall, this review underscores the integral role of thymic tissue homeostasis in the prevention of autoimmune diseases, discussing insights into potential therapeutic strategies and examining putative avenues for future research on developing thymic-based therapies in treating and preventing autoimmune conditions.

Ehf and Fezf2 regulate late medullary thymic epithelial cell and thymic tuft cell development

Thymic epithelial cells are indispensable for T cell maturation and selection and the induction of central immune tolerance. The self-peptide repertoire expressed by medullary thymic epithelial cells is in part regulated by the transcriptional regulator Aire (Autoimmune regulator) and the transcription factor Fezf2. Due to the high complexity of mTEC maturation stages (i.e., post-Aire, Krt10+ mTECs, and Dclk1+ Tuft mTECs) and the heterogeneity in their gene expression profiles (i.e., mosaic expression patterns), it has been challenging to identify the additional factors complementing the transcriptional regulation. We aimed to identify the transcriptional regulators involved in the regulation of mTEC development and self-peptide expression in an unbiased and genome-wide manner. We used ATAC footprinting analysis as an indirect approach to identify transcription factors involved in the gene expression regulation in mTECs, which we validated by ChIP sequencing. This study identifies Fezf2 as a regulator of the recently described thymic Tuft cells (i.e., Tuft mTECs). Furthermore, we identify that transcriptional regulators of the ELF, ESE, ERF, and PEA3 subfamily of the ETS transcription factor family and members of the Krüppel-like family of transcription factors play a role in the transcriptional regulation of genes involved in late mTEC development and promiscuous gene expression.

CDCA7-associated global aberrant DNA hypomethylation translates to localized, tissue-specific transcriptional responses

Disruption of cell division cycle associated 7 (CDCA7) has been linked to aberrant DNA hypomethylation, but the impact of DNA methylation loss on transcription has not been investigated. Here, we show that CDCA7 is critical for maintaining global DNA methylation levels across multiple tissues in vivo. A pathogenic Cdca7 missense variant leads to the formation of large, aberrantly hypomethylated domains overlapping with the B genomic compartment but without affecting the deposition of H3K9 trimethylation (H3K9me3). CDCA7-associated aberrant DNA hypomethylation translated to localized, tissue-specific transcriptional dysregulation that affected large gene clusters. In the brain, we identify CDCA7 as a transcriptional repressor and epigenetic regulator of clustered protocadherin isoform choice. Increased protocadherin isoform expression frequency is accompanied by DNA methylation loss, gain of H3K4 trimethylation (H3K4me3), and increased binding of the transcriptional regulator CCCTC-binding factor (CTCF). Overall, our in vivo work identifies a key role for CDCA7 in safeguarding tissue-specific expression of gene clusters via the DNA methylation pathway.

HOXA3 functions as the on-off switch to regulate the development of hESC-derived third pharyngeal pouch endoderm through EPHB2-mediated Wnt pathway

Objectives

Normal commitment of the endoderm of the third pharyngeal pouch (3PP) is essential for the development and differentiation of the thymus. The aim of this study was to investigate the role of transcription factor HOXA3 in the development and differentiation of 3PP endoderm (3PPE) from human embryonic stem cells (hESCs).

Methods

The 3PPE was differentiated from hESC-derived definitive endoderm (DE) by mimicking developmental queues with Activin A, WNT3A, retinoic acid and BMP4. The function of 3PPE was assessed by further differentiating into functional thymic epithelial cells (TECs). The effect of HOXA3 inhibition on cells of 3PPE was subsequently investigated.

Results

A highly efficient approach for differentiating 3PPE cells was developed and these cells expressed 3PPE related genes HOXA3, SIX1, PAX9 as well as EpCAM. 3PPE cells had a strong potential to develop into TECs which expressed both cortical TEC markers K8 and CD205, and medullary TEC markers K5 and AIRE, and also promoted the development and maturation of T cells. More importantly, transcription factor HOXA3 not only regulated the differentiation of 3PPE, but also had a crucial role for the proliferation and migration of 3PPE cells. Our further investigation revealed that HOXA3 controlled the commitment and function of 3PPE through the regulation of Wnt signaling pathway by activating EPHB2.

Conclusion

Our results demonstrated that HOXA3 functioned as the on-off switch to regulate the development of hESC-derived 3PPE through EPHB2-mediated Wnt pathway, and our findings will provide new insights into studying the development of 3PP and thymic organ in vitro and in vivo.

Expression of the transcription factor Klf6 by thymic epithelial cells is required for thymus development

Thymic epithelial cells (TEC) control T cell development and play essential roles in establishing self-tolerance. By using Foxn1-Cre-driven ablation of Klf6 gene in TEC, we identified Klf6 as a critical factor in TEC development. Klf6 deficiency resulted in a hypoplastic thymus-evident from fetal stages into adulthood-in which a dramatic increase in the frequency of apoptotic TEC was observed. Among cortical TEC (cTEC), a previously unreported cTEC population expressing the transcription factor Sox10 was relatively expanded. Within medullary TEC (mTEC), mTEC I and Tuft-like mTEC IV were disproportionately decreased. Klf6 deficiency altered chromatin accessibility and affected TEC chromatin configuration. Consistent with these defects, naïve conventional T cells and invariant natural killer T cells were reduced in the spleen. Late stages of T cell receptor-dependent selection of thymocytes were affected, and mice exhibited autoimmunity. Thus, Klf6 has a prosurvival role and affects the development of specific TEC subsets contributing to thymic function.

MAGNETO: Cell type marker panel generator from single-cell transcriptomic data

Single-cell RNA sequencing experiments produce data useful to identify different cell types, including uncharacterized and rare ones. This enables us to study the specific functional roles of these cells in different microenvironments and contexts. After identifying a (novel) cell type of interest, it is essential to build succinct marker panels, composed of a few genes referring to cell surface proteins and clusters of differentiation molecules, able to discriminate the desired cells from the other cell populations. In this work, we propose a fully-automatic framework called MAGNETO, which can help construct optimal marker panels starting from a single-cell gene expression matrix and a cell type identity for each cell. MAGNETO builds effective marker panels solving a tailored bi-objective optimization problem, where the first objective regards the identification of the genes able to isolate a specific cell type, while the second conflicting objective concerns the minimization of the total number of genes included in the panel. Our results on three public datasets show that MAGNETO can identify marker panels that identify the cell populations of interest better than state-of-the-art approaches. Finally, by fine-tuning MAGNETO, our results demonstrate that it is possible to obtain marker panels with different specificity levels.

Protein-level mutant p53 reporters identify druggable rare precancerous clones in noncancerous tissues

Detecting and targeting precancerous cells in noncancerous tissues is a major challenge for cancer prevention. Massive stabilization of mutant p53 (mutp53) proteins is a cancer-specific event that could potentially mark precancerous cells, yet in vivo protein-level mutp53 reporters are lacking. Here we developed two transgenic protein-level mutp53 reporters, p53R172H-Akaluc and p53-mCherry, that faithfully mimic the dynamics and function of mutp53 proteins in vivo. Using these reporters, we identified and traced rare precancerous clones in deep noncancerous tissues in various cancer models. In classic mutp53-driven thymic lymphoma models, we found that precancerous clones exhibit broad chromosome number variations, upregulate precancerous stage-specific genes such as Ybx3 and enhance amino acid transport and metabolism. Inhibiting amino acid transporters downstream of Ybx3 at the early but not late stage effectively suppresses tumorigenesis and prolongs survival. Together, these protein-level mutp53 reporters reveal undercharacterized features and vulnerabilities of precancerous cells during early tumorigenesis, paving the way for precision cancer prevention.

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.

Investigating Thymic Epithelial Cell Diversity Using Systems Biology

The thymus is an intricate organ consisting of a diverse population of thymic epithelial cells (TECs). Cortical and medullary TECs and their subpopulations have distinct roles in coordinating the development and selection of functionally competent and self-tolerant T cells. Recent advances made in technologies such as single-cell RNA sequencing have made it possible to investigate and resolve the heterogeneity in TECs. These findings have provided further understanding of the molecular mechanisms regulating TEC function and expression of tissue-restricted Ags. In this brief review, we focus on the newly characterized subsets of TECs and their diversity in relation to their functions in supporting T cell development. We also discuss recent discoveries in expression of self-antigens in the context of TEC development as well as the cellular and molecular changes occurring during embryonic development to thymic involution.

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.

Revisiting Aire and tissue-restricted antigens at single-cell resolution

The thymus is a highly specialized organ that plays an indispensable role in the establishment of self-tolerance, a process characterized by the "education" of developing T-cells. To provide competent T-cells tolerant to self-antigens, medullary thymic epithelial cells (mTECs) orchestrate negative selection by ectopically expressing a wide range of genes, including various tissue-restricted antigens (TRAs). Notably, recent advancements in the high-throughput single-cell analysis have revealed remarkable heterogeneity in mTECs, giving us important clues for dissecting the mechanisms underlying TRA expression. We overview how recent single-cell studies have furthered our understanding of mTECs, with a focus on the role of Aire in inducing mTEC heterogeneity to encompass TRAs.

Revelations in Thymic Epithelial Cell Biology and Heterogeneity from Single-Cell RNA Sequencing and Lineage Tracing Methodologies

Thymic epithelial cells (TECs) make up the thymic microenvironments that support the generation of a functionally competent and self-tolerant T-cell repertoire. Cortical (c)TECs, present in the cortex, are essential for early thymocyte development including selection of thymocytes expressing functional TCRs (positive selection). Medullary (m)TECs, located in the medulla, play a key role in late thymocyte development, including depletion of self-reactive T cells (negative selection) and selection of regulatory T cells. In recent years, transcriptomic analysis by single-cell (sc)RNA sequencing (Seq) has revealed TEC heterogeneity previously masked by population-level RNA-Seq or phenotypic studies. We summarize the discoveries made possible by scRNA-Seq, including the identification of novel mTEC subsets, advances in understanding mTEC promiscuous gene expression, and TEC alterations from embryonic to adult stages. Whereas pseudotime analyses of scRNA-Seq data can suggest relationships between TEC subsets, experimental methods such as lineage tracing and reaggregate thymic organ culture (RTOC) are required to test these hypotheses. Lineage tracing - namely, of β5t or Aire expressing cells - has exposed progenitor and parent-daughter cellular relationships within TEC.

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.

In Silico Stage-Matching of Human, Marmoset, Mouse, and Pig Embryos to Enhance Organ Development Through Interspecies Chimerism

Currently, there is a significant shortage of transplantable organs for patients in need. Interspecies chimerism and blastocyst complementation are alternatives for generating transplantable human organs in host animals such as pigs to meet this shortage. While successful interspecies chimerism and organ generation have been observed between evolutionarily close species such as rat and mouse, barriers still exist for more distant species pairs such as human-mouse, marmoset-mouse, human-pig, and others. One of the proposed barriers to chimerism is the difference in developmental stages between the donor cells and the host embryo at the time the cells are introduced into the host embryo. Hence, there is a logical effort to stage-match the donor cells with the host embryos for enhancing interspecies chimerism. In this study, we used an in silico approach to simultaneously stage-match the early developing embryos of four species, including human, marmoset, mouse, and pig based on transcriptome similarities. We used an unsupervised clustering algorithm to simultaneously stage-match all four species as well as Spearman's correlation analyses to stage-match pairs of donor-host species. From our stage-matching analyses, we found that the four stages that best matched with each other are the human blastocyst (E6/E7), the gastrulating mouse embryo (E6-E6.75), the marmoset late inner cell mass, and the pig late blastocyst. We further demonstrated that human pluripotent stem cells best matched with the mouse post-implantation stages. We also performed ontology analysis of the genes upregulated and commonly expressed between donor-host species pairs at their best matched stages. The stage-matching results predicted by this study will inform in vivo and in vitro interspecies chimerism and blastocyst complementation studies and can be used to match donor cells with host embryos between multiple species pairs to enhance chimerism for organogenesis.

Embryonic ILC-poiesis across tissues

The family of innate lymphoid cells (ILCs), consisting of Group 1 ILCs (natural killer cells and ILC1), ILC2, and ILC3, are critical effectors of innate immunity, inflammation, and homeostasis post-natally, but also exert essential functions before birth. Recent studies during critical developmental periods in the embryo have hinted at complex waves of tissue colonization, and highlighted the breadth of multipotent and committed ILC progenitors from both classic fetal hematopoietic organs such as the liver, as well as tissue sites such as the lung, thymus, and intestine. Assessment of the mechanisms driving cell fate and function of the ILC family in the embryo will be vital to the understanding ILC biology throughout fetal life and beyond.

Integration of scATAC-Seq with scRNA-Seq Data

Single-cell studies are enabling our understanding of the molecular processes of normal cell development and the onset of several pathologies. For instance, single-cell RNA sequencing (scRNA-Seq) measures the transcriptome-wide gene expression at a single-cell resolution, allowing for studying the heterogeneity among the cells of the same population and revealing complex and rare cell populations. On the other hand, single-cell Assay for Transposase-Accessible Chromatin using sequencing (scATAC-Seq) can be used to define transcriptional and epigenetic changes by analyzing the chromatin accessibility at the single-cell level. However, the integration of multi-omics data still remains one of the most difficult tasks in bioinformatics. In this chapter, we focus on the combination of scRNA-Seq and scATACSeq data to perform an integrative analysis of the single-cell transcriptome and chromatin accessibility of human fetal progenitors.

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.

Mesenchymal cell replacement corrects thymic hypoplasia in murine models of 22q11.2 deletion syndrome

22q11.2 deletion syndrome (22q11.2DS) is the most common human chromosomal microdeletion, causing developmentally linked congenital malformations, thymic hypoplasia, hypoparathyroidism, and/or cardiac defects. Thymic hypoplasia leads to T cell lymphopenia, which most often results in mild SCID. Despite decades of research, the molecular underpinnings leading to thymic hypoplasia in 22q11.2DS remain unknown. Comparison of embryonic thymuses from mouse models of 22q11.2DS (Tbx1neo2/neo2) revealed proportions of mesenchymal, epithelial, and hematopoietic cell types similar to those of control thymuses. Yet, the small thymuses were growth restricted in fetal organ cultures. Replacement of Tbx1neo2/neo2 thymic mesenchymal cells with normal ones restored tissue growth. Comparative single-cell RNA-Seq of embryonic thymuses uncovered 17 distinct cell subsets, with transcriptome differences predominant in the 5 mesenchymal subsets from the Tbx1neo2/neo2 cell line. The transcripts affected included those for extracellular matrix proteins, consistent with the increased collagen deposition we observed in the small thymuses. Attenuating collagen cross-links with minoxidil restored thymic tissue expansion for hypoplastic lobes. In colony-forming assays, the Tbx1neo2/neo2-derived mesenchymal cells had reduced expansion potential, in contrast to the normal growth of thymic epithelial cells. These findings suggest that mesenchymal cells were causal to the small embryonic thymuses in the 22q11.2DS mouse models, which was correctable by substitution with normal mesenchyme.

Mesenchymal stromal cells in the thymus

The microenvironment of the thymus is composed of a group of stromal cells that include endoderm-derived thymic epithelial cells (TECs) and mesenchymal stromal cells such as fibroblasts and serves as a site for the development of T cells. TECs are known to play an essential role in T cell differentiation and selection. Mesenchymal stromal cells have been less studied in terms of their immunological significance compared to TECs. Recently, new technologies have made it possible to identify and characterize mesenchymal stromal cells in the thymus, revealing their unique functions in thymic organogenesis and T cell development. This review outlines the current views on mesenchymal stromal cells in the thymus, particularly highlighting the newly discovered function of thymic fibroblasts in T cell repertoire selection.

TOX deficiency facilitates the differentiation of IL-17A-producing γδ T cells to drive autoimmune hepatitis

The specification of the αβ/γδ lineage and the maturation of medullary thymic epithelial cells (mTECs) coordinate central tolerance to self-antigens. However, the mechanisms underlying this biological process remain poorly clarified. Here, we report that dual-stage loss of TOX in thymocytes hierarchically impaired mTEC maturation, promoted thymic IL-17A-producing γδ T-cell (Tγδ17) lineage commitment, and led to the development of fatal autoimmune hepatitis (AIH) via different mechanisms. Transfer of γδ T cells from TOX-deficient mice reproduced AIH. TOX interacted with and stabilized the TCF1 protein to maintain the balance of γδ T-cell development in thymic progenitors, and overexpression of TCF1 normalized αβ/γδ lineage specification and activation. In addition, TOX expression was downregulated in γδ T cells from AIH patients and was inversely correlated with the AIH diagnostic score. Our findings suggest multifaceted roles of TOX in autoimmune control involving mTEC and Tγδ17 development and provide a potential diagnostic marker for AIH.

Thymocyte regulatory variant alters transcription factor binding and protects from type 1 diabetes in infants

We recently mapped a genetic susceptibility locus on chromosome 6q22.33 for type 1 diabetes (T1D) diagnosed below the age of 7 years between the PTPRK and thymocyte-selection-associated (THEMIS) genes. As the thymus plays a central role in shaping the T cell repertoire, we aimed to identify the most likely causal genetic factors behind this association using thymocyte genomic data. In four thymocyte populations, we identified 253 DNA sequence motifs underlying histone modifications. The G insertion allele of rs138300818, associated with protection from diabetes, created thymocyte motifs for multiple histone modifications and thymocyte types. In a parallel approach to identifying variants that alter transcription factor binding motifs, the same variant disrupted a predicted motif for Rfx7, which is abundantly expressed in the thymus. Chromatin state and RNA sequencing data suggested strong transcription overlapping rs138300818 in fetal thymus, while expression quantitative trait locus and chromatin conformation data associate the insertion with lower THEMIS expression. Extending the analysis to other T1D loci further highlighted rs66733041 affecting the GATA3 transcription factor binding in the AFF3 locus. Taken together, our results support a role for thymic THEMIS gene expression and the rs138300818 variant in promoting the development of early-onset T1D.

The divergence between T cell and innate lymphoid cell fates controlled by E and Id proteins

T cells develop in the thymus from lymphoid primed multipotent progenitors or common lymphoid progenitors into αβ and γδ subsets. The basic helix-loop-helix transcription factors, E proteins, play pivotal roles at multiple stages from T cell commitment to maturation. Inhibitors of E proteins, Id2 and Id3, also regulate T cell development while promoting ILC differentiation. Recent findings suggest that the thymus can also produce innate lymphoid cells (ILCs). In this review, we present current findings that suggest the balance between E and Id proteins is likely to be critical for controlling the bifurcation of T cell and ILC fates at early stages of T cell development.

Type 2 cytokines in the thymus activate Sirpα+ dendritic cells to promote clonal deletion

The thymus contains a diversity of dendritic cells (DCs) that exist in defined locations and have different antigen-processing and -presenting features. This suggests that they play nonredundant roles in mediating thymocyte selection. In an effort to eliminate SIRPα+ classic DC2 subsets, we discovered that a substantial proportion expresses the surface lectin, CD301b, in the thymus. These cells resemble the CD301b+ type 2 immune response promoting DCs that are present in the skin-draining lymph nodes. Transcriptional and phenotypic comparison to other DC subsets in the thymus revealed that thymic CD301b+ cDCs represent an activated state that exhibits enhanced antigen processing and presentation. Furthermore, a CD301b+ cDC2 subset demonstrated a type 2 cytokine signature and required steady-state interleukin-4 receptor signaling. Selective ablation of CD301b+ cDC2 subsets impaired clonal deletion without affecting regulatory T cells (Treg cells). The T cell receptor α repertoire sequencing confirmed that a cDC2 subset promotes deletion of conventional T cells with minimal effect on Treg cell selection. Together, these findings suggest that cytokine-induced activation of DCs in the thymus substantially enforces central tolerance.

TCF-1 promotes chromatin interactions across topologically associating domains in T cell progenitors

The high mobility group (HMG) transcription factor TCF-1 is essential for early T cell development. Although in vitro biochemical assays suggest that HMG proteins can serve as architectural elements in the assembly of higher-order nuclear organization, the contribution of TCF-1 on the control of three-dimensional (3D) genome structures during T cell development remains unknown. Here, we investigated the role of TCF-1 in 3D genome reconfiguration. Using gain- and loss-of-function experiments, we discovered that the co-occupancy of TCF-1 and the architectural protein CTCF altered the structure of topologically associating domains in T cell progenitors, leading to interactions between previously insulated regulatory elements and target genes at late stages of T cell development. The TCF-1-dependent gain in long-range interactions was linked to deposition of active enhancer mark H3K27ac and recruitment of the cohesin-loading factor NIPBL at active enhancers. These data indicate that TCF-1 has a role in controlling global genome organization during T cell development.

hECA: The cell-centric assembly of a cell atlas

The accumulation of massive single-cell omics data provides growing resources for building biomolecular atlases of all cells of human organs or the whole body. The true assembly of a cell atlas should be cell-centric rather than file-centric. We developed a unified informatics framework for seamless cell-centric data assembly and built the human Ensemble Cell Atlas (hECA) from scattered data. hECA v1.0 assembled 1,093,299 labeled human cells from 116 published datasets, covering 38 organs and 11 systems. We invented three new methods of atlas applications based on the cell-centric assembly: “in data” cell sorting for targeted data retrieval with customizable logic expressions, “quantitative portraiture” for multi-view representations of biological entities, and customizable reference creation for generating references for automatic annotations. Case studies on agile construction of user-defined sub-atlases and “in data” investigation of CAR-T off-targets in multiple organs showed the great potential enabled by the cell-centric ensemble atlas.

•

A unified informatics framework for seamless cell-centric assembly of massive single-cell data

•

Built the general-purpose human Ensemble Cell Atlas (hECA) V1.0 from scattered data

•

Three new methods of applications enabling “in data” cell experiments and portraiture

•

Case studies of agile atlas reconstruction and target therapies side-effect discovery

Integrative analysis of scRNA-seq and scATAC-seq revealed transit-amplifying thymic epithelial cells expressing autoimmune regulator

Medullary thymic epithelial cells (mTECs) are critical for self-tolerance induction in T cells via promiscuous expression of tissue-specific antigens (TSAs), which are controlled by the transcriptional regulator, AIRE. Whereas AIRE-expressing (Aire⁺) mTECs undergo constant turnover in the adult thymus, mechanisms underlying differentiation of postnatal mTECs remain to be discovered. Integrative analysis of single-cell assays for transposase-accessible chromatin (scATAC-seq) and single-cell RNA sequencing (scRNA-seq) suggested the presence of proliferating mTECs with a specific chromatin structure, which express high levels of Aire and co-stimulatory molecules, CD80 (Aire⁺CD80^hi). Proliferating Aire⁺CD80^hi mTECs detected using Fucci technology express a minimal number of Aire-dependent TSAs and are converted into quiescent Aire⁺CD80^hi mTECs expressing high levels of TSAs after a transit amplification. These data provide evidence for the existence of transit-amplifying Aire⁺mTEC precursors during the Aire⁺mTEC differentiation process of the postnatal thymus.

Generation and clinical potential of functional T lymphocytes from gene-edited pluripotent stem cells

Engineered T cells have been shown to be highly effective in cancer immunotherapy, although T cell exhaustion presents a challenge for their long-term function. Additional T-cell sources must be exploited to broaden the application of engineered T cells for immune defense and reconstitution. Unlimited sources of pluripotent stem cells (PSCs) have provided a potential opportunity to generate precise-engineered therapeutic induced T (iT) cells. Single-cell transcriptome analysis of PSC-derived induced hematopoietic stem and progenitor cells (iHSPC)/iT identified the developmental pathways and possibilities of generating functional T cell from PSCs. To date, the PSC-to-iT platforms encounter several problems, including low efficiency of conventional T subset specification, limited functional potential, and restrictions on large-scale application, because of the absence of a thymus-like organized microenvironment. The updated PSC-to-iT platforms, such as the three-dimensional (3D) artificial thymic organoid (ATO) co-culture system and Runx1/Hoxa9-enforced iT lymphopoiesis, provide fresh perspectives for coordinating culture conditions and transcription factors, which may greatly improve the efficiency of T-cell generation greatly. In addition, the improved PSC-to-iT platform coordinating gene editing technologies will provide various functional engineered unconventional or conventional T cells. Furthermore, the clinical applications of PSC-derived immune cells are accelerating from bench to bedside.

The E-Id Axis Instructs Adaptive Versus Innate Lineage Cell Fate Choice and Instructs Regulatory T Cell Differentiation

Immune responses are primarily mediated by adaptive and innate immune cells. Adaptive immune cells, such as T and B cells, evoke antigen-specific responses through the recognition of specific antigens. This antigen-specific recognition relies on the V(D)J recombination of immunoglobulin (Ig) and T cell receptor (TCR) genes mediated by recombination-activating gene (Rag)1 and Rag2 (Rag1/2). In addition, T and B cells employ cell type-specific developmental pathways during their activation processes, and the regulation of these processes is strictly regulated by the transcription factor network. Among these factors, members of the basic helix-loop-helix (bHLH) transcription factor mammalian E protein family, including E12, E47, E2-2, and HEB, orchestrate multiple adaptive immune cell development, while their antagonists, Id proteins (Id1-4), function as negative regulators. It is well established that a majority of T and B cell developmental trajectories are regulated by the transcriptional balance between E and Id proteins (the E-Id axis). E2A is critically required not only for B cell but also for T cell lineage commitment, whereas Id2 and Id3 enforce the maintenance of naïve T cells and naïve regulatory T (Treg) cells. Here, we review the current knowledge of E- and Id-protein function in T cell lineage commitment and Treg cell differentiation.

The Lineage Differentiation and Dynamic Heterogeneity of Thymic Epithelial Cells During Thymus Organogenesis

Although much progress has been made recently in revealing the heterogeneity of the thymic stromal components, the molecular programs of cell lineage divergency and temporal dynamics of thymic epithelial cell (TEC) development are largely elusive. Here, we constructed a single-cell transcriptional landscape of non-hematopoietic cells from mouse thymus spanning embryonic to adult stages, producing transcriptomes of 30,959 TECs. We resolved the transcriptional heterogeneity of developing TECs and highlighted the molecular nature of early TEC lineage determination and cortico-medullary thymic epithelial cell lineage divergency. We further characterized the differentiation dynamics of TECs by clarification of molecularly distinct cell states in the thymus developing trajectory. We also identified a population of Bpifa1⁺ Plet1⁺ mTECs that was preserved during thymus organogenesis and highly expressed tissue-resident adult stem cell markers. Finally, we highlighted the expression of Aire-dependent tissue-restricted antigens mainly in Aire⁺ Csn2⁺ mTECs and Spink5⁺ Dmkn⁺ mTECs in postnatal thymus. Overall, our data provided a comprehensive characterization of cell lineage differentiation, maturation, and temporal dynamics of thymic epithelial cells during thymus organogenesis.

CD90 Marks a Mesenchymal Program in Human Thymic Epithelial Cells In Vitro and In Vivo

Thymic epithelium is critical for the structural integrity of the thymus and for T cell development. Within the fully formed thymus, large numbers of hematopoietic cells shape the thymic epithelium into a scaffold-like structure which bears little similarity to classical epithelial layers, such as those observed in the skin, intestine or pancreas. Here, we show that human thymic epithelial cells (TECs) possess an epithelial identity that also incorporates the expression of mesenchymal cell associated genes, whose expression levels vary between medullary and cortical TECs (m/cTECs). Using pluripotent stem cell (PSC) differentiation systems, we identified a unique population of cells that co-expressed the master TEC transcription factor FOXN1, as well as the epithelial associated marker EPCAM and the mesenchymal associated gene CD90. Using the same serum free culture conditions, we also observed co-expression of EPCAM and CD90 on cultured TECs derived from neonatal human thymus in vitro. Single cell RNA-sequencing revealed these cultured TECs possessed an immature mTEC phenotype and expressed epithelial and mesenchymal associated genes, such as EPCAM, CLDN4, CD90 and COL1A1. Importantly, flow cytometry and single cell RNA-sequencing analysis further confirmed the presence of an EPCAM+CD90+ population in the CD45- fraction of neonatal human thymic stromal cells in vivo. Using the human thymus cell atlas, we found that cTECs displayed more pronounced mesenchymal characteristics than mTECs during embryonic development. Collectively, these results suggest human TECs possess a hybrid gene expression program comprising both epithelial and mesenchymal elements, and provide a basis for the further exploration of thymus development from primary tissues and from the in vitro differentiation of PSCs.

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.

In vitro and in vivo functions of T cells produced in complemented thymi of chimeric mice generated by blastocyst complementation

Blastocyst complementation is an intriguing way of generating humanized animals for organ preparation in regenerative medicine and establishing novel models for drug development. Confirming that complemented organs and cells work normally in chimeric animals is critical to demonstrating the feasibility of blastocyst complementation. Here, we generated thymus-complemented chimeric mice, assessed the efficacy of anti-PD-L1 antibody in tumor-bearing chimeric mice, and then investigated T-cell function. Thymus-complemented chimeric mice were generated by injecting C57BL/6 (B6) embryonic stem cells into Foxn1^nu/nu morulae or blastocysts. Flow cytometry data showed that the chimeric mouse thymic epithelial cells (TECs) were derived from the B6 cells. T cells appeared outside the thymi. Single-cell RNA-sequencing analysis revealed that the TEC gene-expression profile was comparable to that in B6 mice. Splenic T cells of chimeric mice responded very well to anti-CD3 stimulation in vitro; CD4⁺ and CD8⁺ T cells proliferated and produced IFNγ, IL-2, and granzyme B, as in B6 mice. Anti-PD-L1 antibody treatment inhibited MC38 tumor growth in chimeric mice. Moreover, in the chimeras, anti-PD-L1 antibody restored T-cell activation by significantly decreasing PD-1 expression on T cells and increasing IFNγ-producing T cells in the draining lymph nodes and tumors. T cells produced by complemented thymi thus functioned normally in vitro and in vivo. To successfully generate humanized animals by blastocyst complementation, both verification of the function and gene expression profiling of complemented organs/cells in interspecific chimeras will be important in the near future.

Single-cell transcriptomics of LepR-positive skeletal cells reveals heterogeneous stress-dependent stem and progenitor pools

Leptin receptor (LepR)-positive cells are key components of the bone marrow hematopoietic microenvironment, and highly enrich skeletal stem and progenitor cells that maintain homeostasis of the adult skeleton. However, the heterogeneity and lineage hierarchy within this population has been elusive. Using genetic lineage tracing and single-cell RNA sequencing, we found that Lepr-Cre labels most bone marrow stromal cells and osteogenic lineage cells in adult long bones. Integrated analysis of Lepr-Cre-traced cells under homeostatic and stress conditions revealed dynamic changes of the adipogenic, osteogenic, and periosteal lineages. Importantly, we discovered a Notch3+ bone marrow sub-population that is slow-cycling and closely associated with the vasculatures, as well as key transcriptional networks promoting osteo-chondrogenic differentiation. We also identified a Sca-1+ periosteal sub-population with high clonogenic activity but limited osteo-chondrogenic potential. Together, we mapped the transcriptomic landscape of adult LepR+ stem and progenitor cells and uncovered cellular and molecular mechanisms underlying their maintenance and lineage specification.