A Probabilistic Model of the Germinal Center Reaction

In vivo analysis of CRISPR-edited germinal center murine B cells

The germinal center (GC) reaction is crucial for somatic hypermutation, affinity maturation, and the selection of high-affinity B cells, all of which are hallmarks of the humoral immune response. Understanding the distinct roles of various B cell genes is essential for elucidating the selection mechanisms within the GC reaction. Traditionally, studying B cell gene function in the GC reaction involved generating knock-out mice, a highly time-consuming method that necessitates complex vectors. The advent of Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) technology has simplified the creation of knock-out mice. However, even with CRISPR, the generation of knock-out mice still faces challenges, including being time-consuming, costly, having low knock-out efficiency, and raising ethical concerns regarding animal use. To address these challenges, we developed an alternative method to traditional knock-out mouse generation. Our approach entails the ex vivo CRISPR editing of B cells from transgenic donor mice with different B cell receptor affinities followed by their adoptive transfer into recipient mice. We present a cost-effective, rapid, versatile, and adaptable CRISPR-Cas9 method for in vivo loss-of-function studies of individual murine B cell genes within the context of the GC reaction. This method provides a valuable tool for investigating the complex roles of different B cell genes in the GC selection process. As proof of concept, we validated our approach by examining the role of the pro-apoptotic gene Fas in the GC selection process. We adoptively transferred a mix of Fas knock-out (FasKO) low-affinity B cells, Fas wild-type (FasWT) low-affinity B cells, and FasWT high-affinity B cells into recipient mice. From our results, FasKO low-affinity B cells were still outcompeted by the FasWT high-affinity B cells for selection in the GC. An important observation was the accumulation of FasKO low-affinity GC B cells when compared to the FasWT low-affinity B cells, which suggested a role of Fas in the GC selection process.

B cells present a double-sided effect in digestive system tumors: a review for tumor microenvironment

Over the past few years, there has been an increasing interest in investigating tumor-infiltrating lymphocytes. B lymphocytes (B cells) are extensively distributed within tertiary lymphoid structure (TLS) as multifaceted subgroups and are intimately linked to the anti-tumor properties of TLS, as well as the survival and prognostication of individuals. While the investigation of T lymphocytes in the TLS has advanced to the level of clinical practice, the study of B cells remains limited. The principal impediment to the utilization of B cells in immunotherapy is their notable dual impact on tumors. Compared with tumors in other parts and systems, the function of B cells in the microenvironment of digestive system tumors to promote tumors proliferation, differentiation and migration cannot be ignored. Therefore, this review collects the studies of B cell subsets in tumor microenvironments, particularly related single cell sequencing research. The multifaceted role and function of B cells are investigated in esophageal, liver, colorectal, gastric and pancreatic cancers. And through the identification of B cell subsets and specific markers, this review attempts to explain the reasons why B cells produce different tumor-promoting effects in those tumors. The insights gleaned from this review may provide potential help and support the development of B cell-based immunotherapies.

Lymphatic distribution considerations for subunit vaccine design and development

Subunit vaccines are an important platform for controlling current and emerging infectious diseases. The lymph nodes are the primary site generating the humoral response and delivery of antigens to these sites is critical to effective immunization. Indeed, the duration of antigen exposure within the lymph node is correlated with the antibody response. While current licensed vaccines are typically given through the intramuscular route, injecting vaccines subcutaneously allows for direct access to lymphatic vessels and therefore can enhance the transfer of antigen to the lymph nodes. However, protein subunit antigen uptake into the lymph nodes is inefficient, and subunit vaccines require adjuvants to stimulate the initial immune response. Therefore, formulation strategies have been developed to enhance the exposure of subunit proteins and adjuvants to the lymph nodes by increasing lymphatic uptake or prolonging the retention at the injection site. Given that lymph node exposure is a crucial consideration in vaccine design, in depth analyses of the pharmacokinetics of antigens and adjuvants should be the focus of future preclinical and clinical studies. This review will provide an overview of formulation strategies for targeting the lymphatics and prolonging antigen exposure and will discuss pharmacokinetic evaluations which can be applied toward vaccine development.

B cell memory: from generation to reactivation: a multipronged defense wall against pathogens

Development of B cell memory is a conundrum that scientists are still exploring. Studies have been conducted in vitro and using advanced animal models to elucidate the mechanism underlying the generation of memory B cells (MBCs), the precise roles of MBCs against pathogens, and their protective functions against repeated infections throughout life. Lifelong immunity against invading diseases is mainly the result of overcoming a single infection. This protection is largely mediated by the two main components of B cell memory-MBCs and long-lived plasma cells (PCs). The chemical and cellular mechanisms that encourage fat selection for MBCs or long-lived PCs are an area of active research. Despite the fact that nearly all available vaccinations rely on the capacity to elicit B-cell memory, we have yet to develop successful vaccines that can induce broad-scale protective MBCs against some of the deadliest diseases, including malaria and AIDS. A deeper understanding of the specific cellular and molecular pathways that govern the generation, function, and reactivation of MBCs is critical for overcoming the challenges associated with vaccine development. Here, we reviewed literature on the development of MBCs and their reactivation, interaction with other cell types, strategies against invading pathogens, and function throughout life and discussed the recent advances regarding the key signals and transcription factors which regulate B cell memory and their relevance to the quest for vaccine development.

Efficient evolution of human antibodies from general protein language models

Natural evolution must explore a vast landscape of possible sequences for desirable yet rare mutations, suggesting that learning from natural evolutionary strategies could guide artificial evolution. Here we report that general protein language models can efficiently evolve human antibodies by suggesting mutations that are evolutionarily plausible, despite providing the model with no information about the target antigen, binding specificity or protein structure. We performed language-model-guided affinity maturation of seven antibodies, screening 20 or fewer variants of each antibody across only two rounds of laboratory evolution, and improved the binding affinities of four clinically relevant, highly mature antibodies up to sevenfold and three unmatured antibodies up to 160-fold, with many designs also demonstrating favorable thermostability and viral neutralization activity against Ebola and severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pseudoviruses. The same models that improve antibody binding also guide efficient evolution across diverse protein families and selection pressures, including antibiotic resistance and enzyme activity, suggesting that these results generalize to many settings.

Convergent evolution and B-cell recirculation in germinal centers in a human lymph node

Germinal centers (GCs) play a central role in generating an effective immune response against infectious pathogens, and failures in their regulating mechanisms can lead to the development of autoimmune diseases and cancer. Although previous works study experimental systems of the immune response with mouse models that are immunized with specific antigens, our study focused on a real-life situation, with an ongoing GC response in a human lymph node (LN) involving multiple asynchronized GCs reacting simultaneously to unknown antigens. We combined laser capture microdissection of individual GCs from human LN with next-generation repertoire sequencing to characterize individual GCs as distinct evolutionary spaces. In line with well-characterized GC responses in mice, elicited by immunization with model antigens, we observe a heterogeneous clonal diversity across individual GCs from the same human LN. Still, we identify shared clones in several individual GCs, and phylogenetic tree analysis combined with paratope modeling suggest the re-engagement and rediversification of B-cell clones across GCs and expanded clones exhibiting shared antigen responses across distinct GCs, indicating convergent evolution of the GCs.

New tools for immunologists: models of lymph node function from cells to tissues

The lymph node is a highly structured organ that mediates the body's adaptive immune response to antigens and other foreign particles. Central to its function is the distinct spatial assortment of lymphocytes and stromal cells, as well as chemokines that drive the signaling cascades which underpin immune responses. Investigations of lymph node biology were historically explored in vivo in animal models, using technologies that were breakthroughs in their time such as immunofluorescence with monoclonal antibodies, genetic reporters, in vivo two-photon imaging, and, more recently spatial biology techniques. However, new approaches are needed to enable tests of cell behavior and spatiotemporal dynamics under well controlled experimental perturbation, particularly for human immunity. This review presents a suite of technologies, comprising in vitro, ex vivo and in silico models, developed to study the lymph node or its components. We discuss the use of these tools to model cell behaviors in increasing order of complexity, from cell motility, to cell-cell interactions, to organ-level functions such as vaccination. Next, we identify current challenges regarding cell sourcing and culture, real time measurements of lymph node behavior in vivo and tool development for analysis and control of engineered cultures. Finally, we propose new research directions and offer our perspective on the future of this rapidly growing field. We anticipate that this review will be especially beneficial to immunologists looking to expand their toolkit for probing lymph node structure and function.

Notch Signaling in Acute Inflammation and Sepsis

Notch signaling, a highly conserved pathway in mammals, is crucial for differentiation and homeostasis of immune cells. Besides, this pathway is also directly involved in the transmission of immune signals. Notch signaling per se does not have a clear pro- or anti-inflammatory effect, but rather its impact is highly dependent on the immune cell type and the cellular environment, modulating several inflammatory conditions including sepsis, and therefore significantly impacts the course of disease. In this review, we will discuss the contribution of Notch signaling on the clinical picture of systemic inflammatory diseases, especially sepsis. Specifically, we will review its role during immune cell development and its contribution to the modulation of organ-specific immune responses. Finally, we will evaluate to what extent manipulation of the Notch signaling pathway could be a future therapeutic strategy.

Holistic View on the Structure of Immune Response: Petri Net Model

The simulation of immune response is a challenging task because quantitative data are scarce. Quantitative theoretical models either focus on specific cell-cell interactions or have to make assumptions about parameters. The broad variation of, e.g., the dimensions and abundance between lymph nodes as well as between individual patients hampers conclusive quantitative modeling. No theoretical model has been established representing a consensus on the set of major cellular processes involved in the immune response. In this paper, we apply the Petri net formalism to construct a semi-quantitative mathematical model of the lymph nodes. The model covers the major cellular processes of immune response and fulfills the formal requirements of Petri net models. The intention is to develop a model taking into account the viewpoints of experienced pathologists and computer scientists in the field of systems biology. In order to verify formal requirements, we discuss invariant properties and apply the asynchronous firing rule of a place/transition net. Twenty-five transition invariants cover the model, and each is assigned to a functional mode of the immune response. In simulations, the Petri net model describes the dynamic modes of the immune response, its adaption to antigens, and its loss of memory.

The role of geometric features in a germinal center

The germinal center (GC) is a self-organizing structure produced in the lymphoid follicle during the T-dependent immune response and is an important component of the humoral immune system. However, the impact of the special structure of GC on antibody production is not clear. According to the latest biological experiments, we establish a spatiotemporal stochastic model to simulate the whole self-organization process of the GC including the appearance of two specific zones: the dark zone (DZ) and the light zone (LZ), the development of which serves to maintain an effective competition among different cells and promote affinity maturation. A phase transition is discovered in this process, which determines the critical GC volume for a successful growth in both the stochastic and the deterministic model. Further increase of the volume does not make much improvement on the performance. It is found that the critical volume is determined by the distance between the activated B cell receptor (BCR) and the target epitope of the antigen in the shape space. The observation is confirmed in both 2D and 3D simulations and explains partly the variability of the observed GC size.

From Antibody Repertoires to Cell-Cell Interactions to Molecular Networks: Bridging Scales in the Germinal Center

Antibody-mediated adaptive immunity must provide effective long-term protection with minimal adverse effects, against rapidly mutating pathogens, in a human population with diverse ages, genetics, and immune histories. In order to grasp and leverage the complexities of the antibody response, we advocate for a mechanistic understanding of the multiscale germinal center (GC) reaction - the process by which precursor B-cells evolve high-affinity antigen-specific antibodies, forming an effector repertoire of plasma and memory cells for decades-long protection. The regulatory dynamics of B-cells within the GC are complex, and unfold across multiple interacting spatial and temporal scales. At the organism scale, over weeks to years, the antibody sequence repertoire formed by various B-cell clonal lineages modulates antibody quantity and quality over time. At the tissue and cellular scale, over hours to weeks, B-cells undergo selection via spatially distributed interactions with local stroma, antigen, and helper T-cells. At the molecular scale, over seconds to days, intracellular signaling, transcriptional, and epigenetic networks modulate B-cell fates and shape their clonal lineages. We summarize our current understanding within each of these scales, and identify missing links in connecting them. We suggest that quantitative multi-scale mathematical models of B-cell and GC reaction dynamics provide predictive frameworks that can apply basic immunological knowledge to practical challenges such as rational vaccine design.

Mathematical Modeling of Proliferative Immune Response Initiated by Interactions Between Classical Antigen-Presenting Cells Under Joint Antagonistic IL-2 and IL-4 Signaling

During an adaptive immune response from pathogen invasion, multiple cytokines are produced by various immune cells interacting jointly at the cellular level to mediate several processes. For example, studies have shown that regulation of interleukin-4 (IL-4) correlates with interleukin-2 (IL-2) induced lymphocyte proliferation. This motivates the need to better understand and model the mechanisms driving the dynamic interplay of proliferation of lymphocytes with the complex interaction effects of cytokines during an immune response. To address this challenge, we adopt a hybrid computational approach comprising of continuous, discrete and stochastic non-linear model formulations to predict a system-level immune response as a function of multiple dependent signals and interacting agents including cytokines and targeted immune cells. We propose a hybrid ordinary differential equation-based (ODE) multicellular model system with a stochastic component of antigen microscopic states denoted as Multiscale Multicellular Quantitative Evaluator (MMQE) implemented using MATLAB. MMQE combines well-defined immune response network-based rules and ODE models to capture the complex dynamic interactions between the proliferation levels of different types of communicating lymphocyte agents mediated by joint regulation of IL-2 and IL-4 to predict the emergent global behavior of the system during an immune response. We model the activation of the immune system in terms of different activation protocols of helper T cells by the interplay of independent biological agents of classic antigen-presenting cells (APCs) and their joint activation which is confounded by the exposure time to external pathogens. MMQE quantifies the dynamics of lymphocyte proliferation during pathogen invasion as bivariate distributions of IL-2 and IL-4 concentration levels. Specifically, by varying activation agents such as dendritic cells (DC), B cells and their joint mechanism of activation, we quantify how lymphocyte activation and differentiation protocols boost the immune response against pathogen invasion mediated by a joint downregulation of IL-4 and upregulation of IL-2. We further compare our in-silico results to in-vivo and in-vitro experimental studies for validation. In general, MMQE combines intracellular and extracellular effects from multiple interacting systems into simpler dynamic behaviors for better interpretability. It can be used to aid engineering of anti-infection drugs or optimizing drug combination therapies against several diseases.

Advances in understanding the formation and fate of B-cell memory in response to immunization or infection

Immunological memory has the potential to provide lifelong protection against recurrent infections. As such, it has been crucial to the success of vaccines. Yet, the recent pandemic has illuminated key gaps in our knowledge related to the factors influencing effective memory formation and the inability to predict the longevity of immune protection. In recent decades, researchers have acquired a number of novel and powerful tools with which to study the factors underpinning humoral memory. These tools have been used to study the B-cell fate decisions that occur within the germinal centre (GC), a site where responding B cells undergo affinity maturation and are one of the major routes for memory B cell and high-affinity long-lived plasma cell formation. The advent of single-cell sequencing technology has provided an enhanced resolution for studying fate decisions within the GC and cutting-edge techniques have enabled researchers to model this reaction with more accuracy both in vitro and in silico. Moreover, modern approaches to studying memory B cells have allowed us to gain a better appreciation for the heterogeneity and adaptability of this vital class of B cells. Together, these studies have facilitated important breakthroughs in our understanding of how these systems operate to ensure a successful immune response. In this review, we describe recent advances in the field of GC and memory B-cell biology in order to provide insight into how humoral memory is formed, as well as the potential for generating lasting immunity to novel pathogens such as severe acute respiratory syndrome coronavirus 2.

The Multiple Dimensions of Networks in Cancer: A Perspective

This perspective article gathers the latest developments in mathematical and computational oncology tools that exploit network approaches for the mathematical modelling, analysis, and simulation of cancer development and therapy design. It instigates the community to explore new paths and synergies under the umbrella of the Special Issue “Networks in Cancer: From Symmetry Breaking to Targeted Therapy”. The focus of the perspective is to demonstrate how networks can model the physics, analyse the interactions, and predict the evolution of the multiple processes behind tumour-host encounters across multiple scales. From agent-based modelling and mechano-biology to machine learning and predictive modelling, the perspective motivates a methodology well suited to mathematical and computational oncology and suggests approaches that mark a viable path towards adoption in the clinic.

B Cell Speed and B-FDC Contacts in Germinal Centers Determine Plasma Cell Output via Swiprosin-1/EFhd2

Plasma cells secreting affinity-matured antibodies develop in germinal centers (GCs), where B cells migrate persistently and directionally over defined periods of time. How modes of GC B cell migration influence plasma cell development remained unclear. Through genetic deletion of the F-actin bundling protein Swiprosin-1/EF-hand domain family member 2 (EFhd2) and by two-photon microscopy, we show that EFhd2 restrains B cell speed in GCs and hapten-specific plasma cell output. Modeling the GC reaction reveals that increasing GC B cell speed promotes plasma cell generation. Lack of EFhd2 also reduces contacts of GC B cells with follicular dendritic cells in vivo. Computational modeling uncovers that both GC output and antibody affinity depend quantitatively on contacts of GC B cells with follicular dendritic cells when B cells migrate more persistently. Collectively, our data explain how GC B cells integrate speed and persistence of cell migration with B cell receptor affinity.

Transcriptional regulation of memory B cell differentiation

Memory B cells (MBCs) are critical for the rapid development of protective immunity following re-infection. MBCs capable of neutralizing distinct subclasses of pathogens, such as influenza and HIV, have been identified in humans. However, efforts to develop vaccines that induce broadly protective MBCs to rapidly mutating pathogens have not yet been successful. Better understanding of the signals regulating MBC development and function are essential to overcome current challenges hindering successful vaccine development. Here, we discuss recent advancements regarding the signals and transcription factors regulating germinal centre-derived MBC development and function.

Multiscale Modeling of Germinal Center Recapitulates the Temporal Transition From Memory B Cells to Plasma Cells Differentiation as Regulated by Antigen Affinity-Based Tfh Cell Help

Germinal centers play a key role in the adaptive immune system since they are able to produce memory B cells and plasma cells that produce high affinity antibodies for an effective immune protection. The mechanisms underlying cell-fate decisions are not well understood but asymmetric division of antigen, B-cell receptor affinity, interactions between B-cells and T follicular helper cells (triggering CD40 signaling), and regulatory interactions of transcription factors have all been proposed to play a role. In addition, a temporal switch from memory B-cell to plasma cell differentiation during the germinal center reaction has been shown. To investigate if antigen affinity-based Tfh cell help recapitulates the temporal switch we implemented a multiscale model that integrates cellular interactions with a core gene regulatory network comprising BCL6, IRF4, and BLIMP1. Using this model we show that affinity-based CD40 signaling in combination with asymmetric division of B-cells result in switch from memory B-cell to plasma cell generation during the course of the germinal center reaction. We also show that cell fate division is unlikely to be (solely) based on asymmetric division of Ag but that BLIMP1 is a more important factor. Altogether, our model enables to test the influence of molecular modulations of the CD40 signaling pathway on the production of germinal center output cells.

Cytoskeleton protein 4.1R regulates B-cell fate by modulating the canonical NF-κB pathway

During the immune response, B cells can enter the memory pathway, which is characterized by class switch recombination (CSR), or they may undergo plasma cell differentiation (PCD) to secrete immunoglobulin. Both of these processes occur in activated B cells, which are reported to relate to membrane-association proteins and adaptors. Protein 4.1R acts as an adaptor, linking membrane proteins to the cytoskeleton, and is involved in many cell events such as cell activation and differentiation, and cytokine secretion. However, the effect of 4.1R on regulating B-cell fate is unclear. Here, we show an important association between B-cell fate and 4.1R. In vitro, primary B cells were stimulated with lipopolysaccharide combined with interleukin-4; results showed that 4.1R-deficient (4.1R^-/- ) cells compared with wild-type (4.1R^+/+ ) B cells augmented expression of activation-induced cytidine deaminase and germline, resulting in increased IgG1⁺ B cells, whereas the secretion of IgG1 and IgM was reduced, and CD138⁺ B cells were also decreased. Throughout the process, 4.1R regulated canonical nuclear factor (NF-κB) rather than non-canonical NF-κB to promote the expression of CSR complex components, leading to up-regulation of B-cell CSR. In contrast, 4.1R-deficient B cells showed reduced expression of Blimp-1, which caused B cells to down-regulate PCD. Furthermore, over-activation of canonical NF-κB may induce apoptosis signaling to cause PCD apoptosis to reduce PCD number. In summary, our results suggest that 4.1R acts as a B-cell fate regulator by inhibiting the canonical NF-κB signaling pathway.

Computational Model Reveals a Stochastic Mechanism behind Germinal Center Clonal Bursts

Germinal centers (GCs) are specialized compartments within the secondary lymphoid organs where B cells proliferate, differentiate, and mutate their antibody genes in response to the presence of foreign antigens. Through the GC lifespan, interclonal competition between B cells leads to increased affinity of the B cell receptors for antigens accompanied by a loss of clonal diversity, although the mechanisms underlying clonal dynamics are not completely understood. We present here a multi-scale quantitative model of the GC reaction that integrates an intracellular component, accounting for the genetic events that shape B cell differentiation, and an extracellular stochastic component, which accounts for the random cellular interactions within the GC. In addition, B cell receptors are represented as sequences of nucleotides that mature and diversify through somatic hypermutations. We exploit extensive experimental characterizations of the GC dynamics to parameterize our model, and visualize affinity maturation by means of evolutionary phylogenetic trees. Our explicit modeling of B cell maturation enables us to characterise the evolutionary processes and competition at the heart of the GC dynamics, and explains the emergence of clonal dominance as a result of initially small stochastic advantages in the affinity to antigen. Interestingly, a subset of the GC undergoes massive expansion of higher-affinity B cell variants (clonal bursts), leading to a loss of clonal diversity at a significantly faster rate than in GCs that do not exhibit clonal dominance. Our work contributes towards an in silico vaccine design, and has implications for the better understanding of the mechanisms underlying autoimmune disease and GC-derived lymphomas.

Maintenance of Germinal Center B Cells by Caspase-9 through Promotion of Apoptosis and Inhibition of Necroptosis

In response to T cell-dependent Ag encounter, naive B cells develop into germinal center (GC) B cells, which can further differentiate into Ab-secreting plasma cells or memory B cells. GC B cells are short lived and are prone to caspase-mediated apoptosis. However, how apoptotic caspases regulate GC B cell fate has not been fully characterized. In this study, we show that mice with B cell-specific knockout of caspase-9 had decreases in GC B cells and Ab production after immunization. Caspase-9-deficient B cells displayed defects in caspase-dependent apoptosis but increases in necroptosis signaling. Additional deletion of Ripk3 restored GC B cells and Ab production in mice with B cell-specific knockout of caspase-9. Our results indicate that caspase-9 plays an important role in the maintenance of Ab responses by promoting apoptosis and inhibiting necroptosis in B cells.

What Will B Will B: Identifying Molecular Determinants of Diverse B-Cell Fate Decisions Through Systems Biology

B-cells are the poster child for cellular diversity and heterogeneity. The diverse repertoire of B lymphocytes, each expressing unique antigen receptors, provides broad protection against pathogens. However, B-cell diversity goes beyond unique antigen receptors. Side-stepping B-cell receptor (BCR) diversity through BCR-independent stimuli or engineered organisms with monoclonal BCRs still results in seemingly identical B-cells reaching a wide variety of fates in response to the same challenge. Identifying to what extent the molecular state of a B-cell determines its fate is key to gaining a predictive understanding of B-cells and consequently the ability to control them with targeted therapies. Signals received by B-cells through transmembrane receptors converge on intracellular molecular signaling networks, which control whether each B-cell divides, dies, or differentiates into a number of antibody-secreting distinct B-cell subtypes. The signaling networks that interpret these signals are well known to be susceptible to molecular variability and noise, providing a potential source of diversity in cell fate decisions. Iterative mathematical modeling and experimental studies have provided quantitative insight into how B-cells achieve distinct fates in response to pathogenic stimuli. Here, we review how systems biology modeling of B-cells, and the molecular signaling networks controlling their fates, is revealing the key determinants of cell-to-cell variability in B-cell destiny.

Synchronous Germinal Center Onset Impacts the Efficiency of Antibody Responses

The germinal center reaction is an important target for modulating antibody responses. Antibody production from germinal centers is regulated by a negative feedback mechanism termed antibody feedback. By imposing antibody feedback, germinal centers can interact and regulate the output of other germinal centers. Using an agent-based model of the germinal center reaction, we studied the impact of antibody feedback on kinetics and efficiency of a germinal center. Our simulations predict that high feedback of antibodies from germinal centers reduces the production of plasma cells and subsequently the efficiency of the germinal center reaction by promoting earlier termination. Affinity maturation is only weakly improved by increased antibody feedback and ultimately interrupted because of premature termination of the reaction. The model predicts that the asynchronous onset and changes in number of germinal centers could alter the efficiency of antibody response due to changes in feedback by soluble antibodies. Consequently, late initialized germinal centers have a compromised output due to higher antibody feedback from the germinal centers formed earlier. The results demonstrate potential effects of germinal center intercommunication and highlight the importance of understanding germinal center interactions for optimizing the antibody response, in particular, in the elderly and in the context of vaccination.