Determining the expected variability of immune responses using the Cyton Model

During an adaptive immune response, lymphocytes proliferate for ve to twenty cell divisions, then stop and die over a period of weeks. The cyton model for regulation of lymphocyte proliferation and survival was introduced by Hawkins et al. [17] to provide a framework for understanding this response and its regulation. The model assumes stochastic values for division and survival times for each cell in a responding population. Experimental evidence indicates that the choice of times is drawn from a skewed distribution such as the lognormal, with the fate of individual cells being potentially highly variable. For this reason we calculate the higher moments of the model so that the expected variability can be determined. To do this we formulate a new analytic framework for the cyton model by introducing a generalization to the Bellman-Harris branching process

A minimum of two distinct heritable factors are required to explain correlation structures in proliferating lymphocytes

During the adaptive immune response, lymphocyte populations undergo a characteristic three-phase process: expansion through a series of cell divisions; cessation of expansion; and, finally, most of the accumulated lymphocytes die by apoptosis. The data used, thus far, to inform understanding of these processes, both in vitro and in vivo, are taken from flow cytometry experiments. One significant drawback of flow cytometry is that individual cells cannot be tracked, so that it is not possible to investigate interdependencies in the fate of cells within a family tree. This deficit in experimental information has recently been overcome by Hawkins et al. (Hawkins et al. 2009 Proc. Natl Acad. Sci. USA 106, 13 457–13 462 (doi:10.1073/pnas.0905629106)), who reported on time-lapse microscopy experiments in which B-cells were stimulated through the TLR-9 receptor. Cells stimulated in this way do not aggregate, so that data regarding family trees can be recorded. In this article, we further investigate the Hawkins et al. data. Our conclusions are striking: in order to explain the familial correlation structure in division times, death times and propensity to divide, a minimum of two distinct heritable factors are necessary. As the data show that two distinct factors are necessary, we develop a stochastic model that has two heritable factors and demonstrate that it can reproduce the key features of the data. This model shows that two heritable factors are sufficient. These deductions have a clear impact upon biological understanding of the adaptive immune response. They also necessitate changes to the fundamental premises behind the tools developed by statisticians to draw deductions from flow cytometry data. Finally, they affect the mathematical modelling paradigms that are used to study these systems, as these are widely developed based on assumptions of cellular independence that are not accurate

Early appearance of germinal center–derived memory B cells and plasma cells in blood after primary immunization

Immunization with a T cell–dependent antigen elicits production of specific memory B cells and antibody-secreting cells (ASCs). The kinetic and developmental relationships between these populations and the phenotypic forms they and their precursors may take remain unclear. Therefore, we examined the early stages of a primary immune response, focusing on the appearance of antigen-specific B cells in blood. Within 1 wk, antigen-specific B cells appear in the blood with either a memory phenotype or as immunoglobulin (Ig)G1 ASCs expressing blimp-1. The memory cells have mutated VH genes; respond to the chemokine CXCL13 but not CXCL12, suggesting recirculation to secondary lymphoid organs; uniformly express B220; show limited differentiation potential unless stimulated by antigen; and develop independently of blimp-1 expression. The antigen-specific IgG1 ASCs in blood show affinity maturation paralleling that of bone marrow ASCs, raising the possibility that this compartment is established directly by blood-borne ASCs. We find no evidence for a blimp-1–expressing preplasma memory compartment, suggesting germinal center output is restricted to ASCs and B220+ memory B cells, and this is sufficient to account for the process of affinity maturation

Stochastically Timed Competition Between Division and Differentiation Fates Regulates the Transition From B Lymphoblast to Plasma Cell

In response to external stimuli, naïve B cells proliferate and take on a range of fates important for immunity. How their fate is determined is a topic of much recent research, with candidates including asymmetric cell division, lineage priming, stochastic assignment, and microenvironment instruction. Here we manipulate the generation of plasmablasts from B lymphocytes in vitro by varying CD40 stimulation strength to determine its influence on potential sources of fate control. Using long-term live cell imaging, we directly measure times to differentiate, divide, and die of hundreds of pairs of sibling cells. These data reveal that while the allocation of fates is significantly altered by signal strength, the proportion of siblings identified with asymmetric fates is unchanged. In contrast, we find that plasmablast generation is enhanced by slowing times to divide, which is consistent with a hypothesis of competing timed stochastic fate outcomes. We conclude that this mechanistically simple source of alternative fate regulation is important, and that useful quantitative models of signal integration can be developed based on its principles

Plasma Cell Ontogeny Defined by Quantitative Changes in Blimp-1 Expression

Plasma cells comprise a population of terminally differentiated B cells that are dependent on the transcriptional regulator B lymphocyte–induced maturation protein 1 (Blimp-1) for their development. We have introduced a gfp reporter into the Blimp-1 locus and shown that heterozygous mice express the green fluorescent protein in all antibody-secreting cells (ASCs) in vivo and in vitro. In vitro, these cells display considerable heterogeneity in surface phenotype, immunoglobulin secretion rate, and Blimp-1 expression levels. Importantly, analysis of in vivo ASCs induced by immunization reveals a developmental pathway in which increasing levels of Blimp-1 expression define developmental stages of plasma cell differentiation that have many phenotypic and molecular correlates. Thus, maturation from transient plasmablast to long-lived ASCs in bone marrow is predicated on quantitative increases in Blimp-1 expression

B Cell Receptor–independent Stimuli Trigger Immunoglobulin (Ig) Class Switch Recombination and Production of IgG Autoantibodies by Anergic Self-Reactive B Cells

In both humans and animals, immunoglobulin (Ig)G autoantibodies are less frequent but more pathogenic than IgM autoantibodies, suggesting that controls over Ig isotype switching are required to reinforce B cell self-tolerance. We have used gene targeting to produce mice in which hen egg lysozyme (HEL)-specific B cells can switch to all Ig isotypes (SWHEL mice). When crossed with soluble HEL transgenic (Tg) mice, self-reactive SWHEL B cells became anergic. However, in contrast to anergic B cells from the original nonswitching anti-HEL × soluble HEL double Tg model, self-reactive SWHEL B cells also displayed an immature phenotype, reduced lifespan, and exclusion from the splenic follicle. These differences were not related to their ability to Ig class switch, but instead to competition with non-HEL–binding B cells generated by VH gene replacement in SWHEL mice. When activated in vitro with B cell receptor (BCR)-independent stimuli such as anti-CD40 monoclonal antibody plus interleukin 4 or lipopolysaccharide (LPS), anergic SWHEL double Tg B cells proliferated and produced IgG anti-HEL antibodies as efficiently as naive HEL-binding B cells from SWHEL Ig Tg mice. These results demonstrate that no intrinsic constraints to isotype switching exist in anergic self-reactive B cells. Instead, production of IgG autoantibodies is prevented by separate controls that reduce the likelihood of anergic B cells encountering BCR-independent stimuli. That bacteria-derived LPS could circumvent these controls may explain the well-known association between autoantibody-mediated diseases and episodes of systemic infection

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