Combined transient ablation and single-cell RNA-sequencing reveals the development of medullary thymic epithelial cells

Abstract

Medullary thymic epithelial cells (mTECs) play a critical role in central immune tolerance by mediating negative selection of autoreactive T cells through the collective expression of the peripheral self-antigen compartment, including tissue-specific antigens (TSAs). Recent work has shown that gene-expression patterns within the mTEC compartment are heterogenous and include multiple differentiated cell states. To further define mTEC development and medullary epithelial lineage relationships, we combined lineage tracing and recovery from transient in vivo mTEC ablation with single-cell RNA-sequencing in <i>Mus musculus</i>. The combination of bioinformatic and experimental approaches revealed a non-stem transit-amplifying population of cycling mTECs that preceded <i>Aire</i> expression. We propose a branching model of mTEC development wherein a heterogeneous pool of transit-amplifying cells gives rise to <i>Aire</i>- and <i>Ccl21a</i>-expressing mTEC subsets. We further use experimental techniques to show that within the <i>Aire-</i>expressing developmental branch, TSA expression peaked as <i>Aire</i> expression decreased, implying <i>Aire</i> expression must be established before TSA expression can occur. Collectively, these data provide a roadmap of mTEC development and demonstrate the power of combinatorial approaches leveraging both in vivo models and high-dimensional datasets.