Biologically indeterminate yet ordered promiscuous gene expression in single medullary thymic epithelial cells

During thymic negative selection, medullary thymic epithelial cells (mTEC) collectively express most protein coding genes, a process termed promiscuous gene expression (PGE). Although PGE is crucial for inducing central T-cell tolerance, this process has not been established definitively as being either stochastic or coordinated. To resolve this question, we sequenced the transcriptomes of 6,894 single mTEC, including 1,795 rare cells expressing either of two tissue-restricted antigens, TSPAN8 or GP2. Transcriptional heterogeneity allowed partitioning of mTEC into 15 robustly-defined subpopulations representing distinct maturational stages and subtypes. Although 50 gene co-expression groups were robustly identified, few could be explained by chromosomal location, biological pathway, or tissue specificity. Further, GP2+ mTEC were randomly dispersed spatially within medullary islands. Thus although PGE exhibits ordered co-expression, biologically it is indeterminate. This likely enhances the presentation of diverse antigens to passing thymocytes during their medullary residency, while simultaneously maintaining mTEC identity throughout PGE.