Reversible, tunable epigenetic silencing of TCF1 generates flexibility in the T cell memory decision.

The immune system encodes information about the severity of a pathogenic threat in the quantity and type of memory cells it forms. This encoding emerges from lymphocyte decisions to maintain or lose self-renewal and memory potential during a challenge. By tracking CD8+ T cells at the single-cell and clonal lineage level using time-resolved transcriptomics, quantitative live imaging, and an acute infection model, we find that T cells will maintain or lose memory potential early after antigen recognition. However, following pathogen clearance, T cells may regain memory potential if initially lost. Mechanistically, this flexibility is implemented by a stochastic cis-epigenetic switch that tunably and reversibly silences the memory regulator, TCF1, in response to stimulation. Mathematical modeling shows how this flexibility allows memory T cell numbers to scale robustly with pathogen virulence and immune response magnitudes. We propose that flexibility and stochasticity in cellular decisions ensure optimal immune responses against diverse threats.

Asynchronous combinatorial action of four regulatory factors activates Bcl11b for T cell commitment.

During T cell development, multipotent progenitors relinquish competence for other fates and commit to the T cell lineage by turning on Bcl11b, which encodes a transcription factor. To clarify lineage commitment mechanisms, we followed developing T cells at the single-cell level using Bcl11b knock-in fluorescent reporter mice. Notch signaling and Notch-activated transcription factors collaborate to activate Bcl11b expression irrespectively of Notch-dependent proliferation. These inputs work via three distinct, asynchronous mechanisms: an early locus 'poising' function dependent on TCF-1 and GATA-3, a stochastic-permissivity function dependent on Notch signaling, and a separate amplitude-control function dependent on Runx1, a factor already present in multipotent progenitors. Despite their necessity for Bcl11b expression, these inputs act in a stage-specific manner, providing a multitiered mechanism for developmental gene regulation.

Antigen perception in T cells by long-term Erk and NFAT signaling dynamics.

Immune system threat detection hinges on T cells' ability to perceive varying peptide-major histocompatibility complex (pMHC) antigens. As the Erk and NFAT pathways link T cell receptor engagement to gene regulation, their signaling dynamics may convey information about pMHC inputs. To test this idea, we developed a dual reporter mouse strain and a quantitative imaging assay that, together, enable simultaneous monitoring of Erk and NFAT dynamics in live T cells over day-long timescales as they respond to varying pMHC inputs. Both pathways initially activate uniformly across various pMHC inputs but diverge only over longer (9+ h) timescales, enabling independent encoding of pMHC affinity and dose. These late signaling dynamics are decoded via multiple temporal and combinatorial mechanisms to generate pMHC-specific transcriptional responses. Our findings underscore the importance of long timescale signaling dynamics in antigen perception and establish a framework for understanding T cell responses under diverse contexts.

In search of lost time: Enhancers as modulators of timing in lymphocyte development and differentiation.

Proper timing of gene expression is central to lymphocyte development and differentiation. Lymphocytes often delay gene activation for hours to days after the onset of signaling components, which act on the order of seconds to minutes. Such delays play a prominent role during the intricate choreography of developmental events and during the execution of an effector response. Though a number of mechanisms are sufficient to explain timing at short timescales, it is not known how timing delays are implemented over long timescales that may span several cell generations. Based on the literature, we propose that a class of cis-regulatory elements, termed "timing enhancers," may explain how timing delays are controlled over these long timescales. By considering chromatin as a kinetic barrier to state switching, the timing enhancer model explains experimentally observed dynamics of gene expression where other models fall short. In this review, we elaborate on features of the timing enhancer model and discuss the evidence for its generality throughout development and differentiation. We then discuss potential molecular mechanisms underlying timing enhancer function. Finally, we explore recent evidence drawing connections between timing enhancers and genetic risk for immunopathology. We argue that the timing enhancer model is a useful framework for understanding how cis-regulatory elements control the central dimension of timing in lymphocyte biology.

Multiplexed single-cell profiling of chromatin states at genomic loci by expansion microscopy.

Proper regulation of genome architecture and activity is essential for the development and function of multicellular organisms. Histone modifications, acting in combination, specify these activity states at individual genomic loci. However, the methods used to study these modifications often require either a large number of cells or are limited to targeting one histone mark at a time. Here, we developed a new method called Single Cell Evaluation of Post-TRanslational Epigenetic Encoding (SCEPTRE) that uses Expansion Microscopy (ExM) to visualize and quantify multiple histone modifications at non-repetitive genomic regions in single cells at a spatial resolution of ∼75 nm. Using SCEPTRE, we distinguished multiple histone modifications at a single housekeeping gene, quantified histone modification levels at multiple developmentally-regulated genes in individual cells, and evaluated the relationship between histone modifications and RNA polymerase II loading at individual loci. We find extensive variability in epigenetic states between individual gene loci hidden from current population-averaged measurements. These findings establish SCEPTRE as a new technique for multiplexed detection of combinatorial chromatin states at single genomic loci in single cells.

Unsupervised discovery of dynamic cell phenotypic states from transmitted light movies.

Identification of cell phenotypic states within heterogeneous populations, along with elucidation of their switching dynamics, is a central challenge in modern biology. Conventional single-cell analysis methods typically provide only indirect, static phenotypic readouts. Transmitted light images, on the other hand, provide direct morphological readouts and can be acquired over time to provide a rich data source for dynamic cell phenotypic state identification. Here, we describe an end-to-end deep learning platform, UPSIDE (Unsupervised Phenotypic State IDEntification), for discovering cell states and their dynamics from transmitted light movies. UPSIDE uses the variational auto-encoder architecture to learn latent cell representations, which are then clustered for state identification, decoded for feature interpretation, and linked across movie frames for transition rate inference. Using UPSIDE, we identified distinct blood cell types in a heterogeneous dataset. We then analyzed movies of patient-derived acute myeloid leukemia cells, from which we identified stem-cell associated morphological states as well as the transition rates to and from these states. UPSIDE opens up the use of transmitted light movies for systematic exploration of cell state heterogeneity and dynamics in biology and medicine.

Bcl11b and combinatorial resolution of cell fate in the T-cell gene regulatory network.

T-cell development from hematopoietic progenitors depends on multiple transcription factors, mobilized and modulated by intrathymic Notch signaling. Key aspects of T-cell specification network architecture have been illuminated through recent reports defining roles of transcription factors PU.1, GATA-3, and E2A, their interactions with Notch signaling, and roles of Runx1, TCF-1, and Hes1, providing bases for a comprehensively updated model of the T-cell specification gene regulatory network presented herein. However, the role of lineage commitment factor Bcl11b has been unclear. We use self-organizing maps on 63 RNA-seq datasets from normal and perturbed T-cell development to identify functional targets of Bcl11b during commitment and relate them to other regulomes. We show that both activation and repression target genes can be bound by Bcl11b in vivo, and that Bcl11b effects overlap with E2A-dependent effects. The newly clarified role of Bcl11b distinguishes discrete components of commitment, resolving how innate lymphoid, myeloid, and dendritic, and B-cell fate alternatives are excluded by different mechanisms.

Hematopoiesis and T-cell specification as a model developmental system.

The pathway to generate T cells from hematopoietic stem cells guides progenitors through a succession of fate choices while balancing differentiation progression against proliferation, stage to stage. Many elements of the regulatory system that controls this process are known, but the requirement for multiple, functionally distinct transcription factors needs clarification in terms of gene network architecture. Here, we compare the features of the T-cell specification system with the rule sets underlying two other influential types of gene network models: first, the combinatorial, hierarchical regulatory systems that generate the orderly, synchronized increases in complexity in most invertebrate embryos; second, the dueling 'master regulator' systems that are commonly used to explain bistability in microbial systems and in many fate choices in terminal differentiation. The T-cell specification process shares certain features with each of these prevalent models but differs from both of them in central respects. The T-cell system is highly combinatorial but also highly dose-sensitive in its use of crucial regulatory factors. The roles of these factors are not always T-lineage-specific, but they balance and modulate each other's activities long before any mutually exclusive silencing occurs. T-cell specification may provide a new hybrid model for gene networks in vertebrate developmental systems.

A far downstream enhancer for murine Bcl11b controls its T-cell specific expression.

Bcl11b is a T-cell specific gene in hematopoiesis that begins expression during T-lineage commitment and is required for this process. Aberrant expression of BCL11B or proto-oncogene translocation to the vicinity of BCL11B can be a contributing factor in human T-ALL. To identify the mechanism that controls its distinctive T-lineage expression, we corrected the identified Bcl11b transcription start site and mapped a cell-type-specific differentially methylated region bracketing the Bcl11b promoter. We identified a 1.9-kb region 850 kb downstream of Bcl11b, "Major Peak," distinguished by its dynamic histone marking pattern in development that mirrors the pattern at the Bcl11b promoter. Looping interactions between promoter-proximal elements including the differentially methylated region and downstream elements in the Major Peak are required to recapitulate the T-cell specific expression of Bcl11b in stable reporter assays. Functional dissection of the Major Peak sequence showed distinct subregions, in which TCF-1 sites and a conserved element were required for T-lineage-specific activation and silencing in non-T cells. A bacterial artificial chromosome encompassing the full Bcl11b gene still required the addition of the Major Peak to exhibit T-cell specific expression. Thus, promoter-proximal and Major Peak sequences are cis-regulatory elements that interact over 850 kb to control expression of Bcl11b in hematopoietic cells.

Early generation of a precursor CD8 T cell that can adapt to acute or chronic viral infection.

Virus specific PD-1+ TCF-1+ TOX+ stem-like CD8+ T cells are essential for maintaining T cell responses during chronic infection and are also critical for PD-1 directed immunotherapy. In this study we have used the mouse model of chronic LCMV infection to examine when these virus specific stem-like CD8+ T cells are generated during the course of chronic infection and what is the role of antigen in maintaining the stem-like program. We found that these stem-like CD8+ T cells are generated early (day 5) during chronic infection and that antigen is essential for maintaining their stem-like program. This early generation of stem-like CD8+ T cells suggested that the fate commitment to this cell population was agnostic to the eventual outcome of infection and the immune system prepares a priori for a potential chronic infection. Indeed, we found that an identical virus specific stem-cell like CD8+ T cell population was also generated during acute LCMV infection but these cells were lost once the virus was cleared. To determine the fate of these early PD-1+TCF-1+TOX+ stem-like CD8+ T cells that are generated during both acute and chronic LCMV infection we set up two reciprocal adoptive transfer experiments. In the first experiment we transferred day 5 stem-like CD8+ T cells from chronically infected into acutely infected mice and examined their differentiation after viral clearance. We found that these early stem-like CD8+ T cells downregulated canonical markers of the chronic stem-like CD8+ T cells and expressed markers (CD127 and CD62L) associated with central memory CD8+ T cells. In the second experiment, we transferred day 5 stem-like cells from acutely infected mice into chronically infected mice and found that these CD8+ T cells could function like resource cells after transfer into a chronic environment by generating effector CD8+ T cells in both lymphoid and non-lymphoid tissues while also maintaining the number of stem-like CD8+ T cells. These findings provide insight into the generation and maintenance of virus specific stem-like CD8+ T cells that play a critical role in chronic viral infection. In particular, our study highlights the early generation of stem-like CD8+ T cells and their ability to adapt to either an acute or chronic infection. These findings are of broad significance since these novel stem-like CD8+ T cells play an important role in not only viral infections but also in cancer and autoimmunity.

Maintenance of mitochondrial oxygen homeostasis by cosubstrate compensation.

Mitochondria maintain a constant rate of aerobic respiration over a wide range of oxygen levels. However, the control strategies underlying oxygen homeostasis are still unclear. Using mathematical modeling, we found that the mitochondrial electron transport chain (ETC) responds to oxygen level changes by undergoing compensatory changes in reduced electron carrier levels. This emergent behavior, which we named cosubstrate compensation (CSC), enables the ETC to maintain homeostasis over a wide of oxygen levels. When performing CSC, our ETC models recapitulated a classic scaling relationship discovered by Chance [Chance B (1965) J. Gen. Physiol. 49:163-165] relating the extent of oxygen homeostasis to the kinetics of mitochondrial electron transport. Analysis of an in silico mitochondrial respiratory system further showed evidence that CSC constitutes the dominant control strategy for mitochondrial oxygen homeostasis during active respiration. Our findings indicate that CSC constitutes a robust control strategy for homeostasis and adaptation in cellular biochemical networks.

Actin behavior in bulk cytoplasm is cell cycle regulated in early vertebrate embryos.

The mechanical properties of cells change as they proceed through the cell cycle, primarily owing to regulation of actin and myosin II. Most models for cell mechanics focus on actomyosin in the cortex and ignore possible roles in bulk cytoplasm. We explored cell cycle regulation of bulk cytoplasmic actomyosin in Xenopus egg extracts, which is almost undiluted cytoplasm from unfertilized eggs. We observed dramatic gelation-contraction of actomyosin in mitotic (M phase) extract where Cdk1 activity is high, but not in interphase (I-phase) extract. In spread droplets, M-phase extract exhibited regular, periodic pulses of gelation-contraction a few minutes apart that continued for many minutes. Comparing actin nucleation, disassembly and myosin II activity between M-phase and I-phase extracts, we conclude that regulation of nucleation is likely to be the most important for cell cycle regulation. We then imaged F-actin in early zebrafish blastomeres using a GFP-Utrophin probe. Polymerization in bulk cytoplasm around vesicles increased dramatically during mitosis, consistent with enhanced nucleation. We conclude that F-actin polymerization in bulk cytoplasm is cell cycle regulated in early vertebrate embryos and discuss possible biological functions of this regulation.

Antigen perception in T cells by long-term Erk and NFAT signaling dynamics.

Immune system threat detection hinges on T cells' ability to perceive varying peptide major-histocompatibility complex (pMHC) antigens. As the Erk and NFAT pathways link T cell receptor engagement to gene regulation, their signaling dynamics may convey information about pMHC inputs. To test this idea, we developed a dual reporter mouse strain and a quantitative imaging assay that, together, enable simultaneous monitoring of Erk and NFAT dynamics in live T cells over day-long timescales as they respond to varying pMHC inputs. Both pathways initially activate uniformly across various pMHC inputs, but diverge only over longer (9+ hrs) timescales, enabling independent encoding of pMHC affinity and dose. These late signaling dynamics are decoded via multiple temporal and combinatorial mechanisms to generate pMHC-specific transcriptional responses. Our findings underscore the importance of long timescale signaling dynamics in antigen perception, and establish a framework for understanding T cell responses under diverse contexts. SIGNIFICANCE STATEMENT: To counter diverse pathogens, T cells mount distinct responses to varying peptide-major histocompatibility complex ligands (pMHCs). They perceive the affinity of pMHCs for the T cell receptor (TCR), which reflects its foreignness, as well as pMHC abundance. By tracking signaling responses in single living cells to different pMHCs, we find that T cells can independently perceive pMHC affinity vs dose, and encode this information through the dynamics of Erk and NFAT signaling pathways downstream of the TCR. These dynamics are jointly decoded by gene regulatory mechanisms to produce pMHC-specific activation responses. Our work reveals how T cells can elicit tailored functional responses to diverse threats and how dysregulation of these responses may lead to immune pathologies.

Spatial patterning of metabolism by mitochondria, oxygen, and energy sinks in a model cytoplasm.

Metabolite gradients might guide mitochondrial localization in cells and angiogenesis in tissues. It is unclear whether they can exist in single cells, because the length scale of most cells is small compared to the expected diffusion times of metabolites. For investigation of metabolic gradients, we need experimental systems in which spatial patterns of metabolism can be systematically measured and manipulated. We used concentrated cytoplasmic extracts from Xenopus eggs as a model cytoplasm, and visualized metabolic gradients formed in response to spatial stimuli. Restriction of oxygen supply to the edge of a drop mimicked distance to the surface of a single cell, or distance from a blood vessel in tissue. We imaged a step-like increase of Nicotinamide adenine dinucleotide (NAD) reduction approximately 600 microm distant from the oxygen source. This oxic-anoxic switch was preceded on the oxic side by a gradual rise of mitochondrial transmembrane potential (Deltapsi) and reactive oxygen species (ROS) production, extending over approximately 600 microm and approximately 300 microm, respectively. Addition of Adenosine triphosphate (ATP)-consuming beads mimicked local energy sinks in the cell. We imaged Deltapsi gradients with a decay length of approximately 50-300 microm around these beads, in the first visualization of an energy demand signaling gradient. Our study demonstrates that mitochondria can pattern the cytoplasm over length scales that are suited to convey morphogenetic information in large cells and tissues and provides a versatile model system for probing of the formation and function of metabolic gradients.

Rapid actin monomer-insensitive depolymerization of Listeria actin comet tails by cofilin, coronin, and Aip1.

Actin filaments in cells depolymerize rapidly despite the presence of high concentrations of polymerizable G actin. Cofilin is recognized as a key regulator that promotes actin depolymerization. In this study, we show that although pure cofilin can disassemble Listeria monocytogenes actin comet tails, it cannot efficiently disassemble comet tails in the presence of polymerizable actin. Thymus extracts also rapidly disassemble comet tails, and this reaction is more efficient than pure cofilin when normalized to cofilin concentration. By biochemical fractionation, we identify Aip1 and coronin as two proteins present in thymus extract that facilitate the cofilin-mediated disassembly of Listeria comet tails. Together, coronin and Aip1 lower the amount of cofilin required to disassemble the comet tail and permit even low concentrations of cofilin to depolymerize actin in the presence of polymerizable G actin. The cooperative activities of cofilin, coronin, and Aip1 should provide a biochemical basis for understanding how actin filaments can grow in some places in the cell while shrinking in others.

Dynamic stabilization of actin filaments.

We report here that actin filaments in vitro exist in two populations with significantly different shrinkage rates. Newly polymerized filaments shrink rapidly, primarily from barbed ends, at 1.8/s, but as they age they switch to a stable state that shrinks slowly, primarily from pointed ends, at approximately 0.1/s. This dynamic filament stabilization runs opposite to the classical prediction that actin filaments become more unstable with age as they hydrolyze their bound ATP and release phosphate. Upon cofilin treatment, aged filaments revert to a dynamic state that shows accelerated shrinkage from both ends at a combined rate of 5.9/s. In light of recent electron microscopy studies [Orlova A, et al. (2004) Actin-destabilizing factors disrupt filaments by means of a time reversal of polymerization. Proc Natl Acad Sci USA 101:17664-17668], we propose that dynamic stabilization arises from rearrangement of the filament structure from a relatively disordered state immediately after polymerization to the canonical Holmes helix, a change that is reversed by cofilin binding. Our results suggest that plasticity in the internal structure of the actin filament may play a fundamental role in regulating actin dynamics and may help cells build actin assemblies with vastly different turnover rates.

Scalable control of developmental timetables by epigenetic switching networks.

During development, progenitor cells follow timetables for differentiation that span many cell generations. These developmental timetables are robustly encoded by the embryo, yet scalably adjustable by evolution, facilitating variation in organism size and form. Epigenetic switches, involving rate-limiting activation steps at regulatory gene loci, control gene activation timing in diverse contexts, and could profoundly impact the dynamics of gene regulatory networks controlling developmental lineage specification. Here, we develop a mathematical framework to model regulatory networks with genes controlled by epigenetic switches. Using this framework, we show that such epigenetic switching networks uphold developmental timetables that robustly span many cell generations, and enable the generation of differentiated cells in precisely defined numbers and fractions. Changes to epigenetic switching networks can readily alter the timing of developmental events within a timetable, or alter the overall speed at which timetables unfold, enabling scalable control over differentiated population sizes. With their robust, yet flexibly adjustable nature, epigenetic switching networks could represent central targets on which evolution acts to manufacture diversity in organism size and form.

Tunable, division-independent control of gene activation timing by a polycomb switch.

During development, progenitors often differentiate many cell generations after receiving signals. These delays must be robust yet tunable for precise population size control. Polycomb repressive mechanisms, involving histone H3 lysine-27 trimethylation (H3K27me3), restrain the expression of lineage-specifying genes in progenitors and may delay their activation and ensuing differentiation. Here, we elucidate an epigenetic switch controlling the T cell commitment gene Bcl11b that holds its locus in a heritable inactive state for multiple cell generations before activation. Integrating experiments and modeling, we identify a mechanism where H3K27me3 levels at Bcl11b, regulated by methyltransferase and demethylase activities, set the time delay at which the locus switches from a compacted, silent state to an extended, active state. This activation delay robustly spans many cell generations, is tunable by chromatin modifiers and transcription factors, and is independent of cell division. With their regulatory flexibility, such timed epigenetic switches may broadly control timing in development.

Quantitative analysis of actin turnover in Listeria comet tails: evidence for catastrophic filament turnover.

Rapid assembly and disassembly (turnover) of actin filaments in cytoplasm drives cell motility and shape remodeling. While many biochemical processes that facilitate filament turnover are understood in isolation, it remains unclear how they work together to promote filament turnover in cells. Here, we studied cellular mechanisms of actin filament turnover by combining quantitative microscopy with mathematical modeling. Using live cell imaging, we found that actin polymer mass decay in Listeria comet tails is very well fit by a simple exponential. By analyzing candidate filament turnover pathways using stochastic modeling, we found that exponential polymer mass decay is consistent with either slow treadmilling, slow Arp2/3-dissociation, or catastrophic bursts of disassembly, but is inconsistent with acceleration of filament turnover by severing. Imaging of single filaments in Xenopus egg extract provided evidence that disassembly by bursting dominates isolated filament turnover in a cytoplasmic context. Taken together, our results point to a pathway where filaments grow transiently from barbed ends, rapidly terminate growth to enter a long-lived stable state, and then undergo a catastrophic burst of disassembly. By keeping filament lengths largely constant over time, such catastrophic filament turnover may enable cellular actin assemblies to maintain their mechanical integrity as they are turning over.

Resolving Cell Cycle Speed in One Snapshot with a Live-Cell Fluorescent Reporter.

Cell proliferation changes concomitantly with fate transitions during reprogramming, differentiation, regeneration, and oncogenesis. Methods to resolve cell cycle length heterogeneity in real time are currently lacking. Here, we describe a genetically encoded fluorescent reporter that captures live-cell cycle speed using a single measurement. This reporter is based on the color-changing fluorescent timer (FT) protein, which emits blue fluorescence when newly synthesized before maturing into a red fluorescent protein. We generated a mouse strain expressing an H2B-FT fusion reporter from a universally active locus and demonstrate that faster cycling cells can be distinguished from slower cycling ones on the basis of the intracellular fluorescence ratio between the FT's blue and red states. Using this reporter, we reveal the native cell cycle speed distributions of fresh hematopoietic cells and demonstrate its utility in analyzing cell proliferation in solid tissues. This system is broadly applicable for dissecting functional heterogeneity associated with cell cycle dynamics in complex tissues.