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.

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.

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.

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.

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.

A stochastic epigenetic switch controls the dynamics of T-cell lineage commitment.

Cell fate decisions occur through the switch-like, irreversible activation of fate-specifying genes. These activation events are often assumed to be tightly coupled to changes in upstream transcription factors, but could also be constrained by cis-epigenetic mechanisms at individual gene loci. Here, we studied the activation of Bcl11b, which controls T-cell fate commitment. To disentangle cis and trans effects, we generated mice where two Bcl11b copies are tagged with distinguishable fluorescent proteins. Quantitative live microscopy of progenitors from these mice revealed that Bcl11b turned on after a stochastic delay averaging multiple days, which varied not only between cells but also between Bcl11b alleles within the same cell. Genetic perturbations, together with mathematical modeling, showed that a distal enhancer controls the rate of epigenetic activation, while a parallel Notch-dependent trans-acting step stimulates expression from activated loci. These results show that developmental fate transitions can be controlled by stochastic cis-acting events on individual loci.

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.

A simplified Bcl-2 network model reveals quantitative determinants of cell-to-cell variation in sensitivity to anti-mitotic chemotherapeutics.

Anti-mitotic drugs constitute a major class of cytotoxic chemotherapeutics used in the clinic, killing cancer cells by inducing prolonged mitotic arrest that activates intrinsic apoptosis. Anti-mitotics-induced apoptosis is known to involve degradation of anti-apoptotic Bcl-2 proteins during mitotic arrest; however, it remains unclear how this mechanism accounts for significant heterogeneity observed in the cell death responses both within and between cancer cell types. To unravel quantitative determinants underlying variability in anti-mitotic drug response, we constructed a single-cell dynamical Bcl-2 network model describing cell death control during mitotic arrest, and constrained the model using experimental data from four representative cancer cell lines. The modeling analysis revealed that, given a variable, slowly accumulating pro-apoptotic signal arising from anti-apoptotic protein degradation, generation of a switch-like apoptotic response requires formation of pro-apoptotic Bak complexes with hundreds of subunits, suggesting a crucial role for high-order cooperativity. Moreover, we found that cell-type variation in susceptibility to drug-induced mitotic death arises primarily from differential expression of the anti-apoptotic proteins Bcl-xL and Mcl-1 relative to Bak. The dependence of anti-mitotic drug response on Bcl-xL and Mcl-1 that we derived from the modeling analysis provides a quantitative measure to predict sensitivity of distinct cancer cells to anti-mitotic drug treatment.

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.

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.

Transcriptional establishment of cell-type identity: dynamics and causal mechanisms of T-cell lineage commitment.

Precursor cell entry into the T-cell developmental pathway can be divided into two phases by the closure of T-lineage commitment. As cells decide against the last alternative options to the T-cell fate, they turn on the transcription factor Bcl11b and silence expression of a group of multipotent progenitor regulatory factors that include hematopoietic transcription factor PU.1. Functional perturbation tests show that Bcl11b is needed for commitment while PU.1 actively participates in keeping open access to alternative fates, until it is silenced; however, PU.1 and Bcl11b both contribute positively to T-cell development. Our recent work reviewed here sheds light on the transcriptional regulatory network that determines the timing and irreversibility of Bcl11b activation, the ways that Notch signaling from the thymic microenvironment restricts the action of PU.1 to prevent it from diverting cells to non-T fates, and the target genes that PU.1 still regulates under the influence of Notch signaling to contribute to T-cell generation. We argue that T-cell development depends on the sequential operation of two interlaced, but mutually antagonistic, gene regulatory networks, one initially supporting expansion before commitment and the other imposing a "terminal" differentiation process on committed cells.

Positive feedback between PU.1 and the cell cycle controls myeloid differentiation.

Regulatory gene circuits with positive-feedback loops control stem cell differentiation, but several mechanisms can contribute to positive feedback. Here, we dissect feedback mechanisms through which the transcription factor PU.1 controls lymphoid and myeloid differentiation. Quantitative live-cell imaging revealed that developing B cells decrease PU.1 levels by reducing PU.1 transcription, whereas developing macrophages increase PU.1 levels by lengthening their cell cycles, which causes stable PU.1 accumulation. Exogenous PU.1 expression in progenitors increases endogenous PU.1 levels by inducing cell cycle lengthening, implying positive feedback between a regulatory factor and the cell cycle. Mathematical modeling showed that this cell cycle-coupled feedback architecture effectively stabilizes a slow-dividing differentiated state. These results show that cell cycle duration functions as an integral part of a positive autoregulatory circuit to control cell fate.

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.

Structural plasticity in actin and tubulin polymer dynamics.

Actin filaments and microtubules polymerize and depolymerize by adding and removing subunits at polymer ends, and these dynamics drive cytoplasmic organization, cell division, and cell motility. Since Wegner proposed the treadmilling theory for actin in 1976, it has largely been assumed that the chemical state of the bound nucleotide determines the rates of subunit addition and removal. This chemical kinetics view is difficult to reconcile with observations revealing multiple structural states of the polymer that influence polymerization dynamics but that are not strictly coupled to the bound nucleotide state. We refer to these phenomena as "structural plasticity" and discuss emerging evidence that they play a central role in polymer dynamics and function.

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.

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.