Differentiation of naive T lymphocytes into type I T helper (Th1) cells requires interferon-gamma and interleukin-12. It is puzzling that interferon-gamma induces the Th1 transcription factor T-bet, whereas interleukin-12 mediates Th1 cell lineage differentiation. We use mathematical modeling to analyze the expression kinetics of T-bet, interferon-gamma, and the IL-12 receptor beta2 chain (IL-12Rbeta2) during Th1 cell differentiation, in the presence or absence of interleukin-12 or interferon-gamma signaling. We show that interferon-gamma induced initial T-bet expression, whereas IL-12Rbeta2 was repressed by T cell receptor (TCR) signaling. The termination of TCR signaling permitted upregulation of IL-12Rbeta2 by T-bet and interleukin-12 signaling that maintained T-bet expression. This late expression of T-bet, accompanied by the upregulation of the transcription factors Runx3 and Hlx, was required to imprint the Th cell for interferon-gamma re-expression. Thus initial polarization and subsequent imprinting of Th1 cells are mediated by interlinked, sequentially acting positive feedback loops of TCR-interferon-gamma-Stat1-T-bet and interleukin-12-Stat4-T-bet signaling.

Understanding how regulatory sequences interact in the context of chromosomal architecture is a central challenge in biology. Chromosome conformation capture revealed that mammalian chromosomes possess a rich hierarchy of structural layers, from multi-megabase compartments to sub-megabase topologically associating domains (TADs) and sub-TAD contact domains. TADs appear to act as regulatory microenvironments by constraining and segregating regulatory interactions across discrete chromosomal regions. However, it is unclear whether other (or all) folding layers share similar properties, or rather TADs constitute a privileged folding scale with maximal impact on the organization of regulatory interactions. Here, we present a novel algorithm named CaTCH that identifies hierarchical trees of chromosomal domains in Hi-C maps, stratified through their reciprocal physical insulation, which is a single and biologically relevant parameter. By applying CaTCH to published Hi-C data sets, we show that previously reported folding layers appear at different insulation levels. We demonstrate that although no structurally privileged folding level exists, TADs emerge as a functionally privileged scale defined by maximal boundary enrichment in CTCF and maximal cell-type conservation. By measuring transcriptional output in embryonic stem cells and neural precursor cells, we show that the likelihood that genes in a domain are coregulated during differentiation is also maximized at the scale of TADs. Finally, we observe that regulatory sequences occur at genomic locations corresponding to optimized mutual interactions at the same scale. Our analysis suggests that the architectural functionality of TADs arises from the interplay between their ability to partition interactions and the specific genomic position of regulatory sequences.