Overlapping gene coexpression patterns in human medullary thymic epithelial cells generate self-antigen diversity

Promiscuous expression of numerous tissue-restricted self-antigens (TRAs) in medullary thymic epithelial cells (mTECs) is essential to safeguard self-tolerance. A distinct feature of promiscuous gene expression is its mosaic pattern (i.e., at a given time, each self-antigen is expressed only in 1-3% of mTECs). How this mosaic pattern is generated at the single-cell level is currently not understood. Here, we show that subsets of human mTECs expressing a particular TRA coexpress distinct sets of genes. We identified three coexpression groups comprising overlapping and complementary gene sets, which preferentially mapped to certain chromosomes and intrachromosomal gene clusters. Coexpressed gene loci tended to colocalize to the same nuclear subdomain. The TRA subsets aligned along progressive differentiation stages within the mature mTEC subset and, in vitro, interconverted along this sequence. Our data suggest that single mTECs shift through distinct gene pools, thus scanning a sizeable fraction of the overall repertoire of promiscuously expressed self-antigens. These findings have implications for the temporal and spatial (re)presentation of self-antigens in the medulla in the context of tolerance induction.