Synergy and antagonism in the integration of BCR and CD40 signals that control B-cell proliferation

In response to infection or vaccination, a successful antibody response must enrich high-affinity antigen-reactive B-cells through positive selection, but eliminate auto-reactive B-cells through negative selection. B-cells receive signals from the B-cell receptor (BCR) which binds the antigen, and the CD40 receptor which is stimulated by neighboring T-cells that also recognize the antigen. How BCR and CD40 signaling are integrated quantitatively to jointly determine B-cell fate decision and proliferation remains unclear. To investigate this, we developed a differential-equations-based model of the BCR and CD40 signaling networks activating NFκB. Our model accurately recapitulates the NFκB dynamics of B-cells stimulated through their BCR and CD40 receptors, correctly predicting that costimulation induces more NFκB activity. However, when linking it to established cell fate decision models of cell survival and cell cycle control, it predicted potentiated population expansion that was not observed experimentally. We found that this discrepancy was due to a time-dependent functional antagonism exacerbated by BCR-induced caspase activity that can trigger apoptosis in founder cells, unless NFκB-induced survival gene expression protects B-cells in time. Guided by model predictions, sequential co-stimulation experiments revealed how the temporal dynamics of BCR and CD40 signaling control the fate decision between negative and positive selection of B-cell clonal expansion. Our quantitative findings highlight a complex non-monotonic integration of BCR and CD40 signals that is controlled by a balance between NFκB and cell-death pathways, and suggest a mechanism for regulating the stringency of B-cell selection during an antibody response.

Programmed death receptor 1 (PD-1) ligand Fc fusion proteins reduce T-cell proliferation in vitro independently of PD-1

Programmed death receptor 1 (PD-1) is an inhibitory receptor on T cells shown to restrain T-cell proliferation. PD-1 immune checkpoint blockade has emerged as a highly promising approach in cancer treatment. Much of our understanding of the function of PD-1 is derived from in vitro T-cell activation assays. Here we set out to further investigate how T cells integrate inhibitory signals such as PD-1 in vitro using the PD-1 agonist, PD-1 ligand 1 (PD-L1) fusion protein (PD-L1.Fc), coimmobilized alongside anti-CD3 agonist monoclonal antibody (mAb) on plates to deliver PD-1 signals to wild-type and PD-1-/- CD8+ T cells. Surprisingly, we found that the PD-L1.Fc fusion protein inhibited T-cell proliferation independently of PD-1. This PD-L1.Fc inhibition was observed in the presence and absence of CD28 and interleukin-2 signaling. Binding of PD-L1.Fc was restricted to PD-1-expressing T cells and thus inhibition was not mediated by the interaction of PD-L1.Fc with CD80 or other yet unknown binding partners. Furthermore, a similar PD-1-independent reduction of T-cell proliferation was observed with plate-bound PD-L2.Fc. Hence, our results suggest that the coimmobilization of PD-1 ligand fusion proteins with anti-CD3 mAb leads to a reduction of T-cell engagement with plate-bound anti-CD3 mAb. This study demonstrates a nonspecific mechanism of T-cell inhibition when PD-L1.Fc or PD-L2.Fc fusion proteins are delivered in a plate-bound coimmobilization assay and highlights the importance of careful optimization of assay systems and reagents when interpreting their influence on T-cell proliferation.

Quantifying Human Naïve B Cell Proliferation Kinetics and Differentiation in Controlled In Vitro Cell Culture

Division tracking dyes like Cell Trace Violet (CTV) enable the quantification of cell proliferation, division, and survival kinetics of human naïve B cell responses in vitro. Human naïve B cells exhibit distinct responses to different stimuli, with CpG and anti-Ig inducing a T cell-independent (TI) response, while CD40L and IL-21 promote a T cell-dependent (TD) response that induces isotype switching and differentiation into antibody-secreting cells (ASCs). Both stimulation methods yield valuable insights into the intrinsic programming of B cell health within individuals, making them useful for clinical investigations. For instance, quantitative analysis from these B cell populations could reveal biologically meaningful measurements such as the average number of division rounds and the time to cells' fate. Here, we describe a novel in vitro culture setup for CTV-labelled human naïve B cells and a method for obtaining precise time-based data on proliferation, division-linked isotype switching, and differentiation.

Epigenetic modulators of B cell fate identified through coupled phenotype-transcriptome analysis

High-throughput methodologies are the cornerstone of screening approaches to identify novel compounds that regulate immune cell function. To identify novel targeted therapeutics to treat immune disorders and haematological malignancies, there is a need to integrate functional cellular information with the molecular mechanisms that regulate changes in immune cell phenotype. We facilitate this goal by combining quantitative methods for dissecting complex simultaneous cell phenotypic effects with genomic analysis. This combination strategy we term Multiplexed Analysis of Cells sequencing (MAC-seq), a modified version of Digital RNA with perturbation of Genes (DRUGseq). We applied MAC-seq to screen compounds that target the epigenetic machinery of B cells and assess altered humoral immunity by measuring changes in proliferation, survival, differentiation and transcription. This approach revealed that polycomb repressive complex 2 (PRC2) inhibitors promote antibody secreting cell (ASC) differentiation in both murine and human B cells in vitro. This is further validated using T cell-dependent immunization in mice. Functional dissection of downstream effectors of PRC2 using arrayed CRISPR screening uncovered novel regulators of B cell differentiation, including Mybl1, Myof, Gas7 and Atoh8. Together, our findings demonstrate that integrated phenotype-transcriptome analyses can be effectively combined with drug screening approaches to uncover the molecular circuitry that drives lymphocyte fate decisions.

Multimericity Amplifies the Synergy of BCR and TLR4 for B Cell Activation and Antibody Class Switching

Sustained signaling through the B cell antigen receptor (BCR) is thought to occur only when antigen(s) crosslink or disperse multiple BCR units, such as by multimeric antigens found on the surfaces of viruses or bacteria. B cell-intrinsic Toll-like receptor (TLR) signaling synergizes with the BCR to induce and shape antibody production, hallmarked by immunoglobulin (Ig) class switch recombination (CSR) of constant heavy chains from IgM/IgD to IgG, IgA or IgE isotypes, and somatic hypermutation (SHM) of variable heavy and light chains. Full B cell differentiation is essential for protective immunity, where class switched high affinity antibodies neutralize present pathogens, memory B cells are held in reserve for future encounters, and activated B cells also serve as semi-professional APCs for T cells. But the rules that fine-tune B cell differentiation remain partially understood, despite their being essential for naturally acquired immunity and for guiding vaccine development. To address this in part, we have developed a cell culture system using splenic B cells from naive mice stimulated with several biotinylated ligands and antibodies crosslinked by streptavidin reagents. In particular, biotinylated lipopolysaccharide (LPS), a Toll-like receptor 4 (TLR4) agonist, and biotinylated anti-IgM were pre-assembled (multimerized) using streptavidin, or immobilized on nanoparticles coated with streptavidin, and used to active B cells in this precisely controlled, high throughput assay. Using B cell proliferation and Ig class switching as metrics for successful B cell activation, we show that the stimuli are both synergistic and dose-dependent. Crucially, the multimerized immunoconjugates are most active over a narrow concentration range. These data suggest that multimericity is an essential requirement for B cell BCR/TLRs ligands, and clarify basic rules for B cell activation. Such studies highlight the importance in determining the choice of single vs multimeric formats of antigen and PAMP agonists during vaccine design and development.

Postmitotic G1 phase survivin drives mitogen-independent cell division of B lymphocytes

The prevailing dogma is that renewed mitogenic signaling is essential to traverse G1 phase of the cell cycle after each division. B lymphocytes undergo multiple mitotic divisions, termed clonal expansion, to expand antigen-specific cells that mediate effective immunity. Here we demonstrate that B cells that have undergone one cell division continue to proliferate even in absence of further mitogenic signals. This mitogen-independent proliferation is accompanied by an altered G1 phase marked by transcriptomic and proteomic features of G2/M. Survivin, a G2/M-specific oncogene, is required in G1 to achieve mitogen-independent proliferation.

B and T lymphocytes of the adaptive immune system undergo proliferative bursts to generate pools of antigen-specific cells for effective immunity. Here we show that in contrast to the canonical view that G1 progression signals are essential after mitosis to reenter S phase, B lymphocytes sustain several rounds of mitogen-independent cell division following the first mitosis. Such division appears to be driven by unique characteristics of the postmitotic G1 phase that has features of S and G2/M phases. Birc5 (survivin), a protein associated with chromosome segregation in G2/M, is expressed in the G1 phase of divided B cells and is necessary for mitogen-independent divisions. The partially active G1 phase and propensity for apoptosis inherited after each division may underlie rapid proliferation and cell death, which are hallmarks of B cell proliferative responses.

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.

Single cell resolution of Plasma Cell fate programming in health and disease

Long considered a homogeneous population dedicated to antibody secretion, plasma cell phenotypic and functional heterogeneity is increasingly recognised. Plasma cells were first segregated based on their maturation level, but the complexity of this subset might well be underestimated by this simple dichotomy. Indeed, in the last decade new functions have been attributed to plasma cells including but not limited to cytokine secretion. However, a proper characterization of plasma cell heterogeneity has remained elusive partly due to technical issues and cellular features that are specific to this cell type. Cell intrinsic and cell extrinsic signals could be at the origin of this heterogeneity. Recent advances in technologies like single cell RNA-seq, ATAC-seq or ChIP-seq on low cell numbers helped to elucidate the fate decision in other cell lineages and similar approaches could be implemented to evaluate the heterogeneous fate of activated B cells in health and disease. Here, we summarized published work shedding some lights on the stimuli and genetic program shaping B cell terminal differentiation at the single cell level in mice and men. We also discuss the fate and heterogeneity of plasma cells during immune responses, vaccination and in the frame of human plasma cell disorders. This article is protected by copyright. All rights reserved.

An IRF4-MYC-mTORC1 Integrated Pathway Controls Cell Growth and the Proliferative Capacity of Activated B Cells during B Cell Differentiation In Vivo

Cell division is an essential component of B cell differentiation to Ab-secreting plasma cells, with critical reprogramming occurring during the initial stages of B cell activation. However, a complete understanding of the factors that coordinate early reprogramming events in vivo remain to be determined. In this study, we examined the initial reprogramming by IRF4 in activated B cells using an adoptive transfer system and mice with a B cell-specific deletion of IRF4. IRF4-deficient B cells responding to influenza, 4-hydroxy-3-nitrophenylacetyl-Ficoll, and LPS divided but stalled during the proliferative response. Gene expression profiling of IRF4-deficient B cells at discrete divisions revealed IRF4 was critical for inducing MYC target genes, oxidative phosphorylation, and glycolysis. Moreover, IRF4-deficient B cells maintained an inflammatory gene expression signature. Complementary chromatin accessibility analyses established a hierarchy of IRF4 activity and identified networks of dysregulated transcription factor families in IRF4-deficient B cells, including E-box binding bHLH family members. Indeed, B cells lacking IRF4 failed to fully induce Myc after stimulation and displayed aberrant cell cycle distribution. Furthermore, IRF4-deficient B cells showed reduced mTORC1 activity and failed to initiate the B cell activation unfolded protein response and grow in cell size. Myc overexpression in IRF4-deficient cells was sufficient to overcome the cell growth defect. Together, these data reveal an IRF4-MYC-mTORC1 relationship critical for controlling cell growth and the proliferative response during B cell differentiation.

Pre-mitotic genome re-organisation bookends the B cell differentiation process

During cellular differentiation chromosome conformation is intricately remodelled to support the lineage-specific transcriptional programs required for initiating and maintaining lineage identity. When these changes occur in relation to cell cycle, division and time in response to cellular activation and differentiation signals has yet to be explored, although it has been proposed to occur during DNA synthesis or after mitosis. Here, we elucidate the chromosome conformational changes in B lymphocytes as they differentiate and expand from a naive, quiescent state into antibody secreting plasma cells. We find gene-regulatory chromosome reorganization in late G1 phase before the first division, and that this configuration is remarkably stable as the cells massively and rapidly clonally expand. A second wave of conformational change occurs as cells terminally differentiate into plasma cells, coincident with increased time in G1 phase. These results provide further explanation for how lymphocyte fate is imprinted prior to the first division. They also suggest that chromosome reconfiguration occurs prior to DNA replication and mitosis, and is linked to a gene expression program that controls the differentiation process required for the generation of immunity.

Three-dimensional genome rewiring during the development of antibody-secreting cells

The development of B lymphocytes into antibody-secreting plasma cells is central to the adaptive immune system in that it confers protective and specific antibody response against invading pathogen. This developmental process involves extensive morphological and functional alterations that begin early after antigenic stimulation. These include chromatin restructuring that is critical in regulating gene expression, DNA rearrangement and other cellular processes. Here we outline the recent understanding of the three-dimensional architecture of the genome, specifically focused on its contribution to the process of B cell activation and terminal differentiation into antibody-secreting cells.

Temporal Analysis of Brd4 Displacement in the Control of B Cell Survival, Proliferation, and Differentiation

JQ1 is a BET-bromodomain inhibitor that has immunomodulatory effects. However, the precise molecular mechanism that JQ1 targets to elicit changes in antibody production is not understood. Our results show that JQ1 induces apoptosis, reduces cell proliferation, and as a consequence, inhibits antibody-secreting cell differentiation. ChIP-sequencing reveals a selective displacement of Brd4 in response to acute JQ1 treatment (<2 h), resulting in specific transcriptional repression. After 8 h, subsequent alterations in gene expression arise as a result of the global loss of Brd4 occupancy. We demonstrate that apoptosis induced by JQ1 is solely attributed to the pro-apoptotic protein Bim (Bcl2l11). Conversely, cell-cycle regulation by JQ1 is associated with multiple Myc-associated gene targets. Our results demonstrate that JQ1 drives temporal changes in Brd4 displacement that results in a specific transcriptional profile that directly affects B cell survival and proliferation to modulate the humoral immune response.

Speculations on the evolution of humoral adaptive immunity

The protection of a multicellular organism from infection, at both cell and humoral levels, has been a tremendous driver of gene selection and cellular response strategies. Here we focus on a critical event in the development of humoral immunity: The transition from principally innate responses to a system of adaptive cell selection, with all the attendant mechanical problems that must be solved in order for it to work effectively. Here we review recent advances, but our major goal is to highlight that the development of adaptive immunity resulted from the adoption, reuse and repurposing of an ancient, autonomous cellular program that combines and exploits three titratable cellular fate timers. We illustrate how this common cell machinery recurs and appears throughout biology, and has been essential for the evolution of complex organisms, at many levels of scale.

Quantifying NK cell growth and survival changes in response to cytokines and regulatory checkpoint blockade helps identify optimal culture and expansion conditions

NK cells are innate lymphocytes critical for immune surveillance, particularly in eradication of metastatic cancer cells and acute antiviral responses. In contrast to T cells, NK cell-mediated immunity is rapid, with spontaneous cytotoxicity and cytokine/chemokine production upon pathogen detection. The renaissance in cancer immunology has cast NK cell biology back into the spotlight with an urgent need for deeper understanding of the regulatory networks that govern NK cell antitumor activity. To this end, we have adapted and refined a series of quantitative cellular calculus methods, previously applied to T and B lymphocytes, to dissect the biologic outcomes of NK cells following stimulation with cytokines (IL-15, IL-12, IL-18) or deletion of genes that regulate NK cell proliferation (Cish), survival (Bcl2l11), and activation-induced-cell-death (AICD; Fas). Our methodology is well suited to delineate effects on division rate, intrinsic apoptosis, and AICD, permitting variables such as population half-life, rate of cell division, and their combined influence on population numbers in response to stimuli to be accurately measured and modelled. Changes in these variables that result from gene deletion, concentration of stimuli, time, and cell density give insight into the dynamics of NK cell responses and serve as a platform to dissect the mechanism of action of putative checkpoints in NK cell activation and novel NK cell immunotherapy agents.

Sample path properties of the average generation of a Bellman-Harris process

Motivated by a recently proposed design for a DNA coded randomised algorithm that enables inference of the average generation of a collection of cells descendent from a common progenitor, here we establish strong convergence properties for the average generation of a super-critical Bellman-Harris process. We further extend those results to a two-type Bellman-Harris process where one type can give rise to the other, but not vice versa. These results further affirm the estimation method's potential utility by establishing its long run accuracy on individual sample-paths, and significantly expanding its remit to encompass cellular development that gives rise to differentiated offspring with distinct population dynamics.

TLR9 Signaling Suppresses the Canonical Plasma Cell Differentiation Program in Follicular B Cells

The relative potency and quality of mouse B cell response to Toll-like receptors (TLRs) signaling varies significantly depending on the B cell subset and on the TLR member being engaged. Although it has been shown that marginal zone cells respond faster than follicular (FO) splenic B cells to TLR4 stimulus, FO B cells retain full capacity to proliferate and generate plasmablasts and plasma cells (PBs/PCs) with 2–3 days delayed kinetics. It is not clear whether this scenario could be extended to other members of the TLR family. Here, using quantitative cell culture conditions optimized for B cell growth and differentiation, we show that TLR9 signaling by CpG, while promoting vigorous proliferation, completely fails to induce differentiation of FO B cells into PBs/PCs. Little or absent Ig secretion following TLR9 stimulus was accompanied by lack of expression of cell surface markers and canonical transcription factors involved in PB/PC differentiation. Moreover, not only TLR9 did not induce plasmocyte differentiation, but it also strongly inhibited the massive PB/PC differentiation of FO B cells triggered by LPS/TLR4. Our study reveals unexpected opposite roles for TLR4 and TLR9 in the control of plasma cell differentiation program and disagrees with previous conclusions obtained in high-density cultures conditions on the generation of plasmocytes by TRL9 signaling. The potential implications of these findings on the role of TLR9 in controlling self-tolerance, clonal sizes and regulation of humoral responses are discussed.

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.

Site-1 protease function is essential for the generation of antibody secreting cells and reprogramming for secretory activity

The unfolded protein response (UPR) and activation of XBP1 is necessary for high secretory efficiency and functional differentiation of antibody secreting cells (ASCs). The UPR additionally includes a branch in which membrane-bound transcription factors, exemplified by ATF6, undergo intramembrane-proteolysis by the sequential action of site-1 (MBTPS1/S1P) and site-2 proteases (MBTPS2/S2P) and release of the cytoplasmic domain as an active transcription factor. Such regulation is shared with a family of CREB3-related transcription factors and sterol regulatory element-binding proteins (SREBPs). Of these, we identify that the CREB3 family member CREB3L2 is strongly induced and activated during the transition from B-cell to plasma cell state. Inhibition of site-1 protease leads to a profound reduction in plasmablast number linked to induction of autophagy. Plasmablasts generated in the presence of site-1 protease inhibitor segregated into CD38high and CD38low populations, the latter characterized by a marked reduction in the capacity to secrete IgG. Site-1 protease inhibition is accompanied by a distinctive change in gene expression associated with amino acid, steroid and fatty acid synthesis pathways. These results demonstrate that transcriptional control of metabolic programs necessary for secretory activity can be targeted via site-1 protease inhibition during ASC differentiation.

Converse Smith-Martin cell cycle kinetics by transformed B lymphocytes

Recent studies using direct live cell imaging have reported that individual B lymphocytes have correlated transit times between their G1 and S/G2/M phases. This finding is in contradiction with the influential model of Smith and Martin that assumed the bulk of the total cell cycle time variation arises in the G1 phase of the cell cycle with little contributed by the S/G2/M phase. Here we extend these studies to examine the relation between cell cycle phase lengths in two B lymphoma cell lines. We report that transformed B lymphoma cells undergo a short G1 period that displays little correlation with the time taken for the subsequent S/G2/M phase. Consequently, the bulk of the variation noted for total division times within a population is found in the S/G2/M phases and not the G1 phase. Models that reverse the expected source of variation and assume a single deterministic time in G1 followed by a lag + exponential distribution for S/G2/M fit the data well. These models can be improved further by adopting two sequential distributions or by using the stretched lognormal model developed for primary lymphocytes. We propose that shortening of G1 transit times and uncoupling from other cell cycle phases may be a hallmark of lymphocyte transformation that could serve as an observable phenotypic marker of cancer evolution.

Modifying clonal selection theory with a probabilistic cell

Problem-solving strategies in immunology currently utilize a series of ad hoc, qualitative variations on a foundation of Burnet's formulation of clonal selection theory. These modifications, including versions of two-signal theory, describe how signals regulate lymphocytes to make important decisions governing self-tolerance and changes to their effector and memory states. These theories are useful but are proving inadequate to explain the observable genesis and control of heterogeneity in cell types, the nonlinear passage of cell fate trajectories and how the input from multiple environmental signals can be integrated at different times and strengths. Here, I argue for a paradigm change to place immune theory on a firmer philosophical and quantitative foundation to resolve these difficulties. This change rejects the notion of identical cell subsets and substitutes the concept of a cell as comprised of autonomous functional mechanical components subject to stochastic variations in construction and operation. The theory aims to explain immunity in terms of cell population dynamics, dictated by the operation of cell machinery, such as randomizing elements, division counters, and fate timers. The effect of communicating signals alone and in combination within this system is determined with a cellular calculus. A series of models developed with these principles can resolve logical cell fate and signaling paradoxes and offer a reinterpretation for how self-non-self discrimination and immune response class are controlled.

The regulation of lymphocyte activation and proliferation

Activation induced proliferation and clonal expansion of antigen specific lymphocytes is a hallmark of the adaptive immune response to pathogens. Recent studies identify two distinct control phases. In the first T and B lymphocytes integrate antigen and additional costimuli to motivate a programmed proliferative burst that ceases with a return to cell quiescence and eventual death. This proliferative burst is autonomously timed, ensuring an appropriate response magnitude whilst preventing uncontrolled expansion. This initial response is subject to further modification and extension by a range of signals that modify, expand and direct the emergence of a rich array of new cell types. Thus, both robust clonal expansion of a small number of antigen specific T cells, and the concurrent emergence of extensive cellular diversity, confers immunity to a vast array of different pathogens. The in vivo response to a given pathogen is made up by the sum of all responding clones and is reproducible and pathogen specific. Thus, a precise description of the regulatory principles governing lymphocyte proliferation, differentiation and survival is essential to a unified understanding of the immune system.

T-cell stimuli independently sum to regulate an inherited clonal division fate

In the presence of antigen and costimulation, T cells undergo a characteristic response of expansion, cessation and contraction. Previous studies have revealed that population-level reproducibility is a consequence of multiple clones exhibiting considerable disparity in burst size, highlighting the requirement for single-cell information in understanding T-cell fate regulation. Here we show that individual T-cell clones resulting from controlled stimulation in vitro are strongly lineage imprinted with highly correlated expansion fates. Progeny from clonal families cease dividing in the same or adjacent generations, with inter-clonal variation producing burst-size diversity. The effects of costimulatory signals on individual clones sum together with stochastic independence; therefore, the net effect across multiple clones produces consistent, but heterogeneous population responses. These data demonstrate that substantial clonal heterogeneity arises through differences in experience of clonal progenitors, either through stochastic antigen interaction or by differences in initial receptor sensitivities.

Why do populations of highly similar T cells have heterogeneous division destinies in response to antigenic stimulus? Here the authors develop a multiplex-dye assay and a mathematical framework to test clonal heterogeneity and show distinction in division destiny is a result of inter-clonal variability as lineage imprinting ensures clones share similar proliferation fates.

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.

Mathematics in modern immunology

Mathematical and statistical methods enable multidisciplinary approaches that catalyse discovery. Together with experimental methods, they identify key hypotheses, define measurable observables and reconcile disparate results. We collect a representative sample of studies in T-cell biology that illustrate the benefits of modelling–experimental collaborations and that have proven valuable or even groundbreaking. We conclude that it is possible to find excellent examples of synergy between mathematical modelling and experiment in immunology, which have brought significant insight that would not be available without these collaborations, but that much remains to be discovered.

Manipulation of B-cell responses with histone deacetylase inhibitors

Histone deacetylase inhibitors (HDACi) are approved for treating certain haematological malignancies, however, recent evidence also illustrates they are modulators of the immune system. In experimental models, HDACi are particularly potent against malignancies originating from the B-lymphocyte lineage. Here we examine the ability of this class of compounds to modify both protective and autoimmune antibody responses. In vitro, HDACi affect B-cell proliferation, survival and differentiation in an HDAC-class-dependent manner. Strikingly, treatment of lupus-prone Mrl/lpr mice with the HDACi panobinostat significantly reduces autoreactive plasma-cell numbers, autoantibodies and nephritis, while other immune parameters remain largely unaffected. Immunized control mice treated with panobinostat or the clinically approved HDACi vorinostat have significantly impaired primary antibody responses, but these treatments surprisingly spare circulating memory B cells. These studies indicate that panobinostat is a potential therapy for B-cell-driven autoimmune conditions and HDACi do not induce major long-term detrimental effects on B-cell memory.

The generation of antibody-secreting plasma cells

The regulation of antibody production is linked to the generation and maintenance of plasmablasts and plasma cells from their B cell precursors. Plasmablasts are the rapidly produced and short-lived effector cells of the early antibody response, whereas plasma cells are the long-lived mediators of lasting humoral immunity. An extraordinary number of control mechanisms, at both the cellular and molecular levels, underlie the regulation of this essential arm of the immune response. Despite this complexity, the terminal differentiation of B cells can be described as a simple probabilistic process that is governed by a central gene-regulatory network and modified by environmental stimuli.

Division time-based amplifiers for stochastic gene expression

While cell-to-cell variability is a phenotypic consequence of gene expression noise, sources of this noise may be complex – apart from intrinsic sources such as the random birth/death of mRNA and stochastic switching between promoter states, there are also extrinsic sources of noise such as cell division where division times are either constant or random.

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.

Cyclooxygenase inhibitors inhibit antibody response through interference with MAPK/ERK pathways and BLIMP-1 inhibition

Fever is a common symptom of illness in children, and although not harmful in itself, fever and its associated symptoms are often treated with antipyretic drugs. A number of national and other guidelines now recommend against their routine use; a conclusion that was initially supported by a study showing that the prophylactic use of paracetamol might reduce antibody response to some vaccine antigens, although data from booster vaccinations are more equivocal. Although in vivo data on the cause of this inhibition are scarce, in vitro data suggests that the cause may be due to inhibition of the mitogen activated protein kinase/extracellular regulated protein kinase pathways, and a subsequent reduction in the process of plasma cell differentiation at the beginning of the antibody response. This suggests that in high-risk patients these drugs could be avoided in the early part of an infection when plasma-cell differentiation is occurring. More data are needed to define this period; until then existing data support the recommendation against the routine use of these drugs.

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.

Quantal and graded stimulation of B lymphocytes as alternative strategies for regulating adaptive immune responses

Lymphocytes undergo a typical response pattern following stimulation in vivo: they proliferate, differentiate to effector cells, cease dividing and predominantly die, leaving a small proportion of long-lived memory and effector cells. This pattern results from cell-intrinsic processes following activation and the influence of external regulation. Here we apply quantitative methods to study B-cell responses in vitro. Our results reveal that B cells stimulated through two Toll-like receptors (TLRs) require minimal external direction to undergo the basic pattern typical of immunity. Altering the stimulus strength regulates the outcome in a quantal manner by varying the number of cells that participate in the response. In contrast, the T-cell-dependent CD40 activation signal induces a response where division times and differentiation rates vary in relation to stimulus strength. These studies offer insight into how the adaptive antibody response may have evolved from simple autonomous response patterns to the highly regulable state that is now observed in mammals.

Quantifying T lymphocyte turnover

Peripheral T cell populations are maintained by production of naive T cells in the thymus, clonal expansion of activated cells, cellular self-renewal (or homeostatic proliferation), and density dependent cell life spans. A variety of experimental techniques have been employed to quantify the relative contributions of these processes. In modern studies lymphocytes are typically labeled with 5-bromo-2'-deoxyuridine (BrdU), deuterium, or the fluorescent dye carboxy-fluorescein diacetate succinimidyl ester (CFSE), their division history has been studied by monitoring telomere shortening and the dilution of T cell receptor excision circles (TRECs) or the dye CFSE, and clonal expansion has been documented by recording changes in the population densities of antigen specific cells. Proper interpretation of such data in terms of the underlying rates of T cell production, division, and death has proven to be notoriously difficult and involves mathematical modeling. We review the various models that have been developed for each of these techniques, discuss which models seem most appropriate for what type of data, reveal open problems that require better models, and pinpoint how the assumptions underlying a mathematical model may influence the interpretation of data. Elaborating various successful cases where modeling has delivered new insights in T cell population dynamics, this review provides quantitative estimates of several processes involved in the maintenance of naive and memory, CD4(+) and CD8(+) T cell pools in mice and men.

A novel statistical analysis and interpretation of flow cytometry data

A recently developed class of models incorporating the cyton model of population generation structure into a conservation-based model of intracellular label dynamics is reviewed. Statistical aspects of the data collection process are quantified and incorporated into a parameter estimation scheme. This scheme is then applied to experimental data for PHA-stimulated CD4+T and CD8+T cells collected from two healthy donors. This novel mathematical and statistical framework is shown to form the basis for accurate, meaningful analysis of cellular behaviour for a population of cells labelled with the dye carboxyfluorescein succinimidyl ester and stimulated to divide.

A division-dependent compartmental model for computing cell numbers in CFSE-based lymphocyte proliferation assays

Some key features of a mathematical description of an immune response are an estimate of the number of responding cells and the manner in which those cells divide, differentiate, and die. The intracellular dye CFSE is a powerful experimental tool for the analysis of a population of dividing cells, and numerous mathematical treatments have been aimed at using CFSE data to describe an immune response [30,31,32,37,38,42,48,49]. Recently, partial differential equation structured population models, with measured CFSE fluorescence intensity as the structure variable, have been shown to accurately fit histogram data obtained from CFSE flow cytometry experiments [18,19,52,54]. In this report, the population of cells is mathematically organized into compartments, with all cells in a single compartment having undergone the same number of divisions. A system of structured partial differential equations is derived which can be fit directly to CFSE histogram data. From such a model, cell counts (in terms of the number of divisions undergone) can be directly computed and thus key biological parameters such as population doubling time and precursor viability can be determined. Mathematical aspects of this compartmental model are discussed, and the model is fit to a data set. As in [18,19], we find temporal and division dependence in the rates of proliferation and death to be essential features of a structured population model for CFSE data. Variability in cellular autofluorescence is found to play a significant role in the data, as well. Finally, the compartmental model is compared to previous work, and statistical aspects of the experimental data are discussed.

Dynamics of the splenic innate-like CD19⁺CD45Rlo cell population from adult mice in homeostatic and activated conditions

In the adult spleen, CD19⁺CD45R(-/lo) (19⁺45R(lo)) lymphocytes of embryonic origin exist as a distinct population to that of the conventional B cell lineage. These cells display a plasmablast phenotype, and they spontaneously secrete IgG1 and IgA, whereas the bone marrow population of 19⁺45R(lo) cells contains B1 progenitors. In this study, we show that 19⁺45R(lo) cells are also present in Peyer's patches and in the spleen throughout the life span of wild-type mice, beginning at postnatal day 7. Although this population is heterogeneous, the surface phenotype of most of these cells distinguishes them from follicular, transitional, marginal zone, and B1 cells. In CBA/CaHN mice, few 19⁺45R(lo) cells were detected at postnatal day 7, and none was observed in the adult spleen. Splenic 19⁺45R(lo) cells exhibited homeostatic BrdU uptake in vivo and actively transcribed cell cycle genes. When transferred to immunodeficient RAG2⁻/⁻γchain⁻/⁻ recipient mice, 19⁺45R(lo) cells survived and differentiated into IgG1- and IgA-plasma cells. Moreover, in vitro stimulation of splenic 19⁺45R(lo) cells with LPS, CpG, BAFF/IL4, and CD40/IL4 induced cell proliferation, IgG1/IgA secretion and the release of IL-10, suggesting a potential immunoregulatory role for this subset of innate-like B cells.

A theory of germinal center B cell selection, division, and exit

High-affinity antibodies are generated in germinal centers in a process involving mutation and selection of B cells. Information processing in germinal center reactions has been investigated in a number of recent experiments. These have revealed cell migration patterns, asymmetric cell divisions, and cell-cell interaction characteristics, used here to develop a theory of germinal center B cell selection, division, and exit (the LEDA model). According to this model, B cells selected by T follicular helper cells on the basis of successful antigen processing always return to the dark zone for asymmetric division, and acquired antigen is inherited by one daughter cell only. Antigen-retaining B cells differentiate to plasma cells and leave the germinal center through the dark zone. This theory has implications for the functioning of germinal centers because compared to previous models, high-affinity antibodies appear one day earlier and the amount of derived plasma cells is considerably larger.

Activated mouse B cells lack expression of granzyme B

Recently, it has been reported that human B cells express and secrete the cytotoxic protease granzyme B (GrB) after stimulation with IL-21 and BCR cross-linking. To date, there are few clues on the function of GrB in B cell biology. As experimental transgenic murine systems should provide insights into these issues, we assayed for GrB in C57BL/6 B cells using an extensive array of physiologically relevant stimuli but were unable to detect either GrB expression or its proteolytic activity, even when Ag-specific transgenic BCRs were engaged. Similar results were also obtained with B cells from DBA/2, CBA, or BALB/c mice. In vivo, infection with either influenza virus or murine γ-herpesvirus induced the expected expression of GrB in CTLs, but not in B cell populations. We also investigated a possible role of GrB on the humoral immune response to the model Ag 4-hydroxy-3-nitrophenylacetyl-keyhole limpet hemocyanin, but GrB-deficient mice produced normal amounts of Ab with typical affinity maturation and a heightened secondary response, demonstrating conclusively the redundancy of GrB for Ab responses. Our results highlight the complex evolutionary differences that have shaped the immune systems of mice and humans. The physiological consequences of GrB expression in human B cells remain unclear, and the current study suggests that experimental mouse models will not be helpful in addressing this issue.

Fc receptor-like 5 promotes B cell proliferation and drives the development of cells displaying switched isotypes

The biological roles of B cell membrane proteins in the FCRL family are enigmatic. FCRL proteins, including FCRL5, were shown to modulate early BCR signaling, although the subsequent, functional consequences of receptor engagement are poorly understood. We found that FCRL5 surface protein itself was induced temporarily upon BCR stimulation of human, naive B cells, indicating precise control over timing of FCRL5 engagement. Cross-linking of FCRL5 on cells induced to express FCRL5 enhanced B cell proliferation significantly. This enhancement required costimulation of the BCR and TLR9, two signals required for optimal proliferation of naive B cells, whereas T cell help in the form of anti-CD40 and IL-2 was dispensable. In addition, we found that FCRL5 stimulation generated a high proportion of cells displaying surface IgG and IgA. Optimal development of cells expressing switched isotypes required T cell help, in addition to stimuli found necessary for enhanced proliferation. Surprisingly, cells that developed upon FCRL5 stimulation simultaneously displayed surface IgM, IgG, and IgA. Cells expressing multiple Ig isotypes were described in hairy cell leukemia, a disease in which FCRL5 is overexpressed. Enhanced proliferation and downstream isotype expression upon FCRL5 stimulation could reflect a physiological role for FCRL5 in the expansion and development of antigen-primed B cells. In addition, FCRL5 may promote growth of malignant cells in hairy cell leukemia and other FCRL5-expressing tumors.

The genetic network controlling plasma cell differentiation

Upon activation by antigen, mature B cells undergo immunoglobulin class switch recombination and differentiate into antibody-secreting plasma cells, the endpoint of the B cell developmental lineage. Careful quantitation of these processes, which are stochastic, independent and strongly linked to the division history of the cell, has revealed that populations of B cells behave in a highly predictable manner. Considerable progress has also been made in the last few years in understanding the gene regulatory network that controls the B cell to plasma cell transition. The mutually exclusive transcriptomes of B cells and plasma cells are maintained by the antagonistic influences of two groups of transcription factors, those that maintain the B cell program, including Pax5, Bach2 and Bcl6, and those that promote and facilitate plasma cell differentiation, notably Irf4, Blimp1 and Xbp1. In this review, we discuss progress in the definition of both the transcriptional and cellular events occurring during late B cell differentiation, as integrating these two approaches is crucial to defining a regulatory network that faithfully reflects the stochastic features and complexity of the humoral immune response.

A new model for the estimation of cell proliferation dynamics using CFSE data

CFSE analysis of a proliferating cell population is a popular tool for the study of cell division and divisionlinked changes in cell behavior. Recently Banks et al. (2011), Luzyanina et al. (2009), Luzyanina et al. (2007), a partial differential equation (PDE) model to describe lymphocyte dynamics in a CFSE proliferation assay was proposed. We present a significant revision of this model which improves the physiological understanding of several parameters. Namely, the parameter used previously as a heuristic explanation for the dilution of CFSE dye by cell division is replaced with a more physical component, cellular autofluorescence. The rate at which label decays is also quantified using a Gompertz decay process. We then demonstrate a revised method of fitting the model to the commonly used histogram representation of the data. It is shown that these improvements result in a model with a strong physiological basis which is fully capable of replicating the behavior observed in the data.

B cells do not take up bacterial DNA: an essential role for antigen in exposure of DNA to toll-like receptor-9

Murine dendritic cells (DC) and macrophages respond to bacterial CpG DNA through toll-like receptor 9 (TLR9). Although it is frequently assumed that bacterial DNA is a direct stimulus for B cells, published work does not reliably show responses of purified B cells. Here we show that purified splenic B cells did not respond to Escherichia coli DNA with induction of CD86, despite readily responding to single-stranded (ss) phosphodiester CpG oligodeoxynucleotides (ODN). This was due to a combination of weak responses to both long and double-stranded (ds) DNA. B-cell DNA uptake was greatly reduced with increasing DNA length. This contrasts with macrophages where DNA uptake and subsequent responses were enhanced with increasing DNA length. However, when DNA was physically linked to hen egg lysozyme (HEL), HEL-specific B cells showed efficient uptake of DNA, and limited proliferation in response to the HEL-DNA complex. We propose that, in the absence of other signals, B cells have poor uptake and responses to long dsDNA to prevent polyclonal activation. Conversely, when DNA is physically linked to a B-cell receptor (BCR) ligand, its uptake is increased, allowing TLR9-dependent B-cell activation in an antigen-specific manner. We could not generate fragments of E. coli DNA by limited DNaseI digestion that could mimic the stimulatory effect of ss CpG ODN on naïve B cells. We suggest that the frequently studied polyclonal B-cell responses to CpG ODN are relevant to therapeutic applications of phosphorothioate-modified CpG-containing ODN, but not to natural responses to foreign or host dsDNA.

Control systems and decision making for antibody production

This paper synthesizes recent progress toward understanding the integrated control systems and fail-safes that guide the quality and quantity of antibody produced by B cells. We focus on four key decisions: (1) the choice between proliferation or death in perifollicular B cells in the first 3 days after antigen encounter; (2) differentiation of proliferating perifollicular B cells into extrafollicular plasma cells or germinal center B cells; (3) positive selection of B cell antigen receptor (BCR) affinity for foreign antigen versus negative selection of BCR affinity for self antigen in germinal center B cells; and (4) survival versus death of antibody-secreting plasma cells. Understanding the engineering of these control systems represents a challenging future step for treating disorders of antibody production in autoimmunity, allergy and immunodeficiency.

Elevated Mcl-1 perturbs lymphopoiesis, promotes transformation of hematopoietic stem/progenitor cells, and enhances drug resistance

Diverse human cancers with poor prognosis, including many lymphoid and myeloid malignancies, exhibit high levels of Mcl-1. To explore the impact of Mcl-1 overexpression on the hematopoietic compartment, we have generated vavP-Mcl-1 transgenic mice. Their lymphoid and myeloid cells displayed increased resistance to a variety of cytotoxic agents. Myelopoiesis was relatively normal, but lymphopoiesis was clearly perturbed, with excess mature B and T cells accumulating. Rather than the follicular lymphomas typical of vavP-BCL-2 mice, aging vavP-Mcl-1 mice were primarily susceptible to lymphomas having the phenotype of a stem/progenitor cell (11 of 30 tumors) or pre-B cell (12 of 30 tumors). Mcl-1 overexpression dramatically accelerated Myc-driven lymphomagenesis. Most vavP-Mcl-1/ Eμ-Myc mice died around birth, and transplantation of blood from bitransgenic E18 embryos into unirradiated mice resulted in stem/progenitor cell tumors. Furthermore, lethally irradiated mice transplanted with E13 fetal liver cells from Mcl-1/Myc bitransgenic mice uniformly died of stem/progenitor cell tumors. When treated in vivo with cyclophosphamide, tumors coexpressing Mcl-1 and Myc transgenes were significantly more resistant than conventional Eμ-Myc lymphomas. Collectively, these results demonstrate that Mcl-1 overexpression renders hematopoietic cells refractory to many cytotoxic insults, perturbs lymphopoiesis and promotes malignant transformation of hematopoietic stem and progenitor cells.

Innate pathways to B-cell activation and tolerance

B cells represent an important link between the adaptive and innate immune systems as they express both antigen-specific B-cell receptors (BCRs) as well as various Toll-like receptors (TLRs). Several checkpoints in B-cell development ensure that self-specific cells are eliminated from the mature B-cell repertoire to avoid harmful autoreactive responses. These checkpoints are controlled by BCR-mediated events but are also influenced by TLR-dependent signals from the innate immune system. Additionally, B-cell-intrinsic and extrinsic TLR signaling are critical for inflammatory events required for the clearance of microbial infections. Factors secreted by TLR-activated macrophages or dendritic cells directly influence the fate of protective and autoreactive B cells. Additionally, naive and memory B cells respond differentially to TLR ligands, as do different B-cell subsets. We review here recent literature describing intrinsic and extrinsic effects of TLR stimulation on the fate of B cells, with particular attention to autoimmune diseases.