The concerted change in the distribution of cell cycle phases and zone composition in germinal centers is regulated by IL-21.

Humoral immune responses require germinal centres (GC) for antibody affinity maturation. Within GC, B cell proliferation and mutation are segregated from affinity-based positive selection in the dark zone (DZ) and light zone (LZ) substructures, respectively. While IL-21 is known to be important in affinity maturation and GC maintenance, here we show it is required for both establishing normal zone representation and preventing the accumulation of cells in the G1 cell cycle stage in the GC LZ. Cell cycle progression of DZ B cells is unaffected by IL-21 availability, as is the zone phenotype of the most highly proliferative GC B cells. Collectively, this study characterises the development of GC zones as a function of time and B cell proliferation and identifies IL-21 as an important regulator of these processes. These data help explain the requirement for IL-21 in normal antibody affinity maturation.

Cyton2: A Model of Immune Cell Population Dynamics That Includes Familial Instructional Inheritance.

Lymphocytes are the central actors in adaptive immune responses. When challenged with antigen, a small number of B and T cells have a cognate receptor capable of recognising and responding to the insult. These cells proliferate, building an exponentially growing, differentiating clone army to fight off the threat, before ceasing to divide and dying over a period of weeks, leaving in their wake memory cells that are primed to rapidly respond to any repeated infection. Due to the non-linearity of lymphocyte population dynamics, mathematical models are needed to interrogate data from experimental studies. Due to lack of evidence to the contrary and appealing to arguments based on Occam's Razor, in these models newly born progeny are typically assumed to behave independently of their predecessors. Recent experimental studies, however, challenge that assumption, making clear that there is substantial inheritance of timed fate changes from each cell by its offspring, calling for a revision to the existing mathematical modelling paradigms used for information extraction. By assessing long-term live-cell imaging of stimulated murine B and T cells in vitro, we distilled the key phenomena of these within-family inheritances and used them to develop a new mathematical model, Cyton2, that encapsulates them. We establish the model's consistency with these newly observed fine-grained features. Two natural concerns for any model that includes familial correlations would be that it is overparameterised or computationally inefficient in data fitting, but neither is the case for Cyton2. We demonstrate Cyton2's utility by challenging it with high-throughput flow cytometry data, which confirms the robustness of its parameter estimation as well as its ability to extract biological meaning from complex mixed stimulation experiments. Cyton2, therefore, offers an alternate mathematical model, one that is, more aligned to experimental observation, for drawing inferences on lymphocyte population dynamics.

Regulatory T Cells Suppress Effector T Cell Proliferation by Limiting Division Destiny.

Understanding how the strength of an effector T cell response is regulated is a fundamental problem in immunology with implications for immunity to pathogens, autoimmunity, and immunotherapy. The initial magnitude of the T cell response is determined by the sum of independent signals from antigen, co-stimulation and cytokines. By applying quantitative methods, the contribution of each signal to the number of divisions T cells undergo (division destiny) can be measured, and the resultant exponential increase in response magnitude accurately calculated. CD4+CD25+Foxp3+ regulatory T cells suppress self-reactive T cell responses and limit pathogen-directed immune responses before bystander damage occurs. Using a quantitative modeling framework to measure T cell signal integration and response, we show that Tregs modulate division destiny, rather than directly increasing the rate of death or delaying interdivision times. The quantitative effect of Tregs could be mimicked by modulating the availability of stimulatory co-stimuli and cytokines or through the addition of inhibitory signals. Thus, our analysis illustrates the primary effect of Tregs on the magnitude of effector T cell responses is mediated by modifying division destiny of responding cell populations.

DiSNE Movie Visualization and Assessment of Clonal Kinetics Reveal Multiple Trajectories of Dendritic Cell Development.

A thorough understanding of cellular development is incumbent on assessing the complexities of fate and kinetics of individual clones within a population. Here, we develop a system for robust periodical assessment of lineage outputs of thousands of transient clones and establishment of bona fide cellular trajectories. We appraise the development of dendritic cells (DCs) in fms-like tyrosine kinase 3 ligand culture from barcode-labeled hematopoietic stem and progenitor cells (HSPCs) by serially measuring barcode signatures and visualize these multidimensional data using developmental interpolated t-distributed stochastic neighborhood embedding (DiSNE) time-lapse movies. We identify multiple cellular trajectories of DC development that are characterized by distinct fate bias and expansion kinetics and determine that these are intrinsically programmed. We demonstrate that conventional DC and plasmacytoid DC trajectories are largely separated already at the HSPC stage. This framework allows systematic evaluation of clonal dynamics and can be applied to other steady-state or perturbed developmental systems.

Machine learning applications in cell image analysis.

Machine learning (ML) refers to a set of automatic pattern recognition methods that have been successfully applied across various problem domains, including biomedical image analysis. This review focuses on ML applications for image analysis in light microscopy experiments with typical tasks of segmenting and tracking individual cells, and modelling of reconstructed lineage trees. After describing a typical image analysis pipeline and highlighting challenges of automatic analysis (for example, variability in cell morphology, tracking in presence of clutters) this review gives a brief historical outlook of ML, followed by basic concepts and definitions required for understanding examples. This article then presents several example applications at various image processing stages, including the use of supervised learning methods for improving cell segmentation, and the application of active learning for tracking. The review concludes with remarks on parameter setting and future directions.

A Myc-dependent division timer complements a cell-death timer to regulate T cell and B cell responses.

T lymphocytes and B lymphocytes integrate activating signals to control the size of their proliferative response. Here we report that such control was achieved by timed changes in the production rate of cell-cycle-regulating proto-oncoprotein Myc, with division cessation occurring when Myc levels fell below a critical threshold. The changing pattern of the level of Myc was not affected by cell division, which identified the regulating mechanism as a cell-intrinsic, heritable temporal controller. Overexpression of Myc in stimulated T cells and B cells did not sustain cell proliferation indefinitely, as a separate 'time-to-die' mechanism, also heritable, was programmed after lymphocyte activation and led to eventual cell loss. Together the two competing cell-intrinsic timed fates created the canonical T cell and B cell immune-response pattern of rapid growth followed by loss of most cells. Furthermore, small changes in these timed processes by regulatory signals, or by oncogenic transformation, acted in synergy to greatly enhance cell numbers over time.

Stochastic Measurement Models for Quantifying Lymphocyte Responses Using Flow Cytometry.

Adaptive immune responses are complex dynamic processes whereby B and T cells undergo division and differentiation triggered by pathogenic stimuli. Deregulation of the response can lead to severe consequences for the host organism ranging from immune deficiencies to autoimmunity. Tracking cell division and differentiation by flow cytometry using fluorescent probes is a major method for measuring progression of lymphocyte responses, both in vitro and in vivo. In turn, mathematical modeling of cell numbers derived from such measurements has led to significant biological discoveries, and plays an increasingly important role in lymphocyte research. Fitting an appropriate parameterized model to such data is the goal of these studies but significant challenges are presented by the variability in measurements. This variation results from the sum of experimental noise and intrinsic probabilistic differences in cells and is difficult to characterize analytically. Current model fitting methods adopt different simplifying assumptions to describe the distribution of such measurements and these assumptions have not been tested directly. To help inform the choice and application of appropriate methods of model fitting to such data we studied the errors associated with flow cytometry measurements from a wide variety of experiments. We found that the mean and variance of the noise were related by a power law with an exponent between 1.3 and 1.8 for different datasets. This violated the assumptions inherent to commonly used least squares, linear variance scaling and log-transformation based methods. As a result of these findings we propose a new measurement model that we justify both theoretically, from the maximum entropy standpoint, and empirically using collected data. Our evaluation suggests that the new model can be reliably used for model fitting across a variety of conditions. Our work provides a foundation for modeling measurements in flow cytometry experiments thus facilitating progress in quantitative studies of lymphocyte responses.

T cell signaling. Antigen affinity, costimulation, and cytokine inputs sum linearly to amplify T cell expansion.

T cell responses are initiated by antigen and promoted by a range of costimulatory signals. Understanding how T cells integrate alternative signal combinations and make decisions affecting immune response strength or tolerance poses a considerable theoretical challenge. Here, we report that T cell receptor (TCR) and costimulatory signals imprint an early, cell-intrinsic, division fate, whereby cells effectively count through generations before returning automatically to a quiescent state. This autonomous program can be extended by cytokines. Signals from the TCR, costimulatory receptors, and cytokines add together using a linear division calculus, allowing the strength of a T cell response to be predicted from the sum of the underlying signal components. These data resolve a long-standing costimulation paradox and provide a quantitative paradigm for therapeutically manipulating immune response strength.

Mechanisms of cell division as regulators of acute immune response.

The acute adaptive immune response is complex, proceeding through phases of activation of quiescent lymphocytes, rapid expansion by cell division and cell differentiation, cessation of division and eventual death of greater than 95 % of the newly generated population. Control of the response is not central but appears to operate as a distributed process where global patterns reliably emerge as a result of collective behaviour of a large number of autonomous cells. In this review, we highlight evidence that competing intracellular timed processes underlie the distribution of individual fates and control cell proliferation, cessation and loss. These principles can be captured in a mathematical model to illustrate consistency with previously published experimentally observed data.

Stretched cell cycle model for proliferating lymphocytes.

Stochastic variation in cell cycle time is a consistent feature of otherwise similar cells within a growing population. Classic studies concluded that the bulk of the variation occurs in the G1 phase, and many mathematical models assume a constant time for traversing the S/G2/M phases. By direct observation of transgenic fluorescent fusion proteins that report the onset of S phase, we establish that dividing B and T lymphocytes spend a near-fixed proportion of total division time in S/G2/M phases, and this proportion is correlated between sibling cells. This result is inconsistent with models that assume independent times for consecutive phases. Instead, we propose a stretching model for dividing lymphocytes where all parts of the cell cycle are proportional to total division time. Data fitting based on a stretched cell cycle model can significantly improve estimates of cell cycle parameters drawn from DNA labeling data used to monitor immune cell dynamics.

Quantitative analysis of Plasmodium ookinete motion in three dimensions suggests a critical role for cell shape in the biomechanics of malaria parasite gliding motility.

Motility is a fundamental part of cellular life and survival, including for Plasmodium parasites--single-celled protozoan pathogens responsible for human malaria. The motile life cycle forms achieve motility, called gliding, via the activity of an internal actomyosin motor. Although gliding is based on the well-studied system of actin and myosin, its core biomechanics are not completely understood. Currently accepted models suggest it results from a specifically organized cellular motor that produces a rearward directional force. When linked to surface-bound adhesins, this force is passaged to the cell posterior, propelling the parasite forwards. Gliding motility is observed in all three life cycle stages of Plasmodium: sporozoites, merozoites and ookinetes. However, it is only the ookinetes--formed inside the midgut of infected mosquitoes--that display continuous gliding without the necessity of host cell entry. This makes them ideal candidates for invasion-free biomechanical analysis. Here we apply a plate-based imaging approach to study ookinete motion in three-dimensional (3D) space to understand Plasmodium cell motility and how movement facilitates midgut colonization. Using single-cell tracking and numerical analysis of parasite motion in 3D, our analysis demonstrates that ookinetes move with a conserved left-handed helical trajectory. Investigation of cell morphology suggests this trajectory may be based on the ookinete subpellicular cytoskeleton, with complementary whole and subcellular electron microscopy showing that, like their motion paths, ookinetes share a conserved left-handed corkscrew shape and underlying twisted microtubular architecture. Through comparisons of 3D movement between wild-type ookinetes and a cytoskeleton-knockout mutant we demonstrate that perturbation of cell shape changes motion from helical to broadly linear. Therefore, while the precise linkages between cellular architecture and actomyosin motor organization remain unknown, our analysis suggests that the molecular basis of cell shape may, in addition to motor force, be a key adaptive strategy for malaria parasite dissemination and, as such, transmission.