Multiple subpopulations in medullary thymic epithelial cells are regulated by Aire with distinctive modes of gene regulations [bulkRNA-seq]

Aire, a transcriptional regulator whose defect results in the development of autoimmunity, controls the transcriptome in the bulk of medullary thymic epithelial cells (mTECs) including the genes for self-antigens. Mechanisms for this process, however, remained incompletely understood, especially regarding the effects of Aire in each subpopulation in mTECs. Here, we combined single-cell multi-omics (scMulti-seq) technologies by which gene expression and open chromatin were simultaneously profiled from the same cell and deeper single-cell sequencing by RamDA-seq which enabled profiling of the minimum dropout. scMulti-seq effectively demonstrated the subpopulations among mTECs associated with diverse transcription factor activities which are affected by Aire and developmental stages, while RamDA-seq revealed the different modes of gene regulations by Aire in different subpopulations. Aire is not only essential for the development of a specific population in which Aire induces a set of genes in a stochastic manner but also promotes the development of multiple subpopulations in which Aire enhances the expressions of genes in a coordinated manner. The former group of Aire-induced genes was characterized by enhanced H3K27me3 chromatin modifications and Aire seems to unleash these silenced genes by deactivating polycomb repressive complex 2 (PRC2). We also found unique gene signatures for the neonatal stage such as high expression of chemokines and expression of cortical TEC (cTEC)-like genes in TECs, which were Aire-independent. Our approach has illuminated the multilayered transcriptional control in mTECs governing the multiple subpopulations organizing the thymic microenvironment. Overall design: Bulk RNA-seq analysis of FACS-sorted mouse thymic epithelial cells (TECs, CD45-EpCAM+Ly51-) from wild-type and knockout (KO) mice