B cells within germinal centers migrate preferentially from dark to light zone

A multiscale spatial modeling framework for the germinal center response

The germinal center response or reaction (GCR) is a hallmark event of adaptive humoral immunity. Unfolding in the B cell follicles of the secondary lymphoid organs, a GC culminates in the production of high-affinity antibody-secreting plasma cells along with memory B cells. By interacting with follicular dendritic cells (FDC) and T follicular helper (Tfh) cells, GC B cells exhibit complex spatiotemporal dynamics. Driving the B cell dynamics are the intracellular signal transduction and gene regulatory network that responds to cell surface signaling molecules, cytokines, and chemokines. As our knowledge of the GC continues to expand in depth and in scope, mathematical modeling has become an important tool to help disentangle the intricacy of the GCR and inform novel mechanistic and clinical insights. While the GC has been modeled at different granularities, a multiscale spatial simulation framework - integrating molecular, cellular, and tissue-level responses - is still rare. Here, we report our recent progress toward this end with a hybrid stochastic GC framework developed on the Cellular Potts Model-based CompuCell3D platform. Tellurium is used to simulate the B cell intracellular molecular network comprising NF-κB, FOXO1, MYC, AP4, CXCR4, and BLIMP1 that responds to B cell receptor (BCR) and CD40-mediated signaling. The molecular outputs of the network drive the spatiotemporal behaviors of B cells, including cyclic migration between the dark zone (DZ) and light zone (LZ) via chemotaxis; clonal proliferative bursts, somatic hypermutation, and DNA damage-induced apoptosis in the DZ; and positive selection, apoptosis via a death timer, and emergence of plasma cells in the LZ. Our simulations are able to recapitulate key molecular, cellular, and morphological GC events, including B cell population growth, affinity maturation, and clonal dominance. This novel modeling framework provides an open-source, customizable, and multiscale virtual GC simulation platform that enables qualitative and quantitative in silico investigations of a range of mechanistic and applied research questions on the adaptive humoral immune response in the future.

A Multiscale Spatial Modeling Framework for the Germinal Center Response

The germinal center response or reaction (GCR) is a hallmark event of adaptive humoral immunity. Unfolding in the B cell follicles of the secondary lymph organs, a GC culminates in the production of high-affinity antibody-secreting plasma cells along with memory B cells. By interacting with follicular dendritic cells (FDC) and T follicular helper (Tfh) cells, GC B cells exhibit complex spatiotemporal dynamics. Driving the B cell dynamics are the intracellular signal transduction and gene regulatory network that responds to cell surface signaling molecules, cytokines, and chemokines. As our knowledge of the GC continues to expand in depth and in scope, mathematical modeling has become an important tool to help disentangle the intricacy of the GCR and inform novel mechanistic and clinical insights. While the GC has been modeled at different granularities, a multiscale spatial simulation framework - integrating molecular, cellular, and tissue-level responses - is still rare. Here, we report our recent progress toward this end with a hybrid stochastic GC framework developed on the Cellular Potts Model-based CompuCell3D platform. Tellurium is used to simulate the B cell intracellular molecular network comprising NF-κB, FOXO1, MYC, AP4, CXCR4, and BLIMP1 that responds to B cell receptor (BCR) and CD40-mediated signaling. The molecular outputs of the network drive the spatiotemporal behaviors of B cells, including cyclic migration between the dark zone (DZ) and light zone (LZ) via chemotaxis; clonal proliferative bursts, somatic hypermutation, and DNA damage-induced apoptosis in the DZ; and positive selection, apoptosis via a death timer, and emergence of plasma cells in the LZ. Our simulations are able to recapitulate key molecular, cellular, and morphological GC events including B cell population growth, affinity maturation, and clonal dominance. This novel modeling framework provides an open-source, customizable, and multiscale virtual GC simulation platform that enables qualitative and quantitative in silico investigations of a range of mechanic and applied research questions in future.

Single-cell spatial metabolomics with cell-type specific protein profiling for tissue systems biology

Metabolic reprogramming in cancer and immune cells occurs to support their increasing energy needs in biological tissues. Here we propose Single Cell Spatially resolved Metabolic (scSpaMet) framework for joint protein-metabolite profiling of single immune and cancer cells in male human tissues by incorporating untargeted spatial metabolomics and targeted multiplexed protein imaging in a single pipeline. We utilized the scSpaMet to profile cell types and spatial metabolomic maps of 19507, 31156, and 8215 single cells in human lung cancer, tonsil, and endometrium tissues, respectively. The scSpaMet analysis revealed cell type-dependent metabolite profiles and local metabolite competition of neighboring single cells in human tissues. Deep learning-based joint embedding revealed unique metabolite states within cell types. Trajectory inference showed metabolic patterns along cell differentiation paths. Here we show scSpaMet's ability to quantify and visualize the cell-type specific and spatially resolved metabolic-protein mapping as an emerging tool for systems-level understanding of tissue biology.

A New Method Based on the von Mises-Fisher Distribution Shows that a Minority of Liver-Localized CD8 T Cells Display Hard-To-Detect Attraction to Plasmodium-Infected Hepatocytes

Malaria is a disease caused by Plasmodium parasites, resulting in over 200 million infections and 400,000 deaths every year. A critical step of malaria infection is when sporozoites, injected by mosquitoes, travel to the liver and form liver stages. Malaria vaccine candidates which induce large numbers of malaria-specific CD8 T cells in mice are able to eliminate all liver stages, preventing fulminant malaria. However, how CD8 T cells find all parasites in 48 h of the liver stage lifespan is not well understood. Using intravital microscopy of murine livers, we generated unique data on T cell search for malaria liver stages within a few hours after infection. To detect attraction of T cells to an infection site, we used the von Mises-Fisher distribution in 3D, similar to the 2D von Mises distribution previously used in ecology. Our results suggest that the vast majority (70-95%) of malaria-specific and non-specific liver-localized CD8 T cells did not display attraction towards the infection site, suggesting that the search for malaria liver stages occurs randomly. However, a small fraction (15-20%) displayed weak but detectable attraction towards parasites which already had been surrounded by several T cells. We found that speeds and turning angles correlated with attraction, suggesting that understanding mechanisms that determine the speed of T cell movement in the liver may improve the efficacy of future T cell-based vaccines. Stochastic simulations suggest that a small movement bias towards the parasite dramatically reduces the number of CD8 T cells needed to eliminate all malaria liver stages, but to detect such attraction by individual cells requires data from long imaging experiments which are not currently feasible. Importantly, as far as we know this is the first demonstration of how activated/memory CD8 T cells might search for the pathogen in nonlymphoid tissues a few hours after infection. We have also established a framework for how attraction of individual T cells towards a location in 3D can be rigorously evaluated.

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.

Compartments and Connections Within the Germinal Center

Protective high affinity antibody responses emerge through an orchestrated developmental process that occurs in germinal centers (GCs). While GCs have been appreciated since 1930, a wealth of recent progress provides new insights into the molecular and cellular dynamics governing humoral immunity. In this review, we highlight advances that demonstrate that fundamental GC B cell function, selection, proliferation and SHM occur within distinct cell states. The resulting new model provides new opportunities to understand the evolution of immunity in infectious, autoimmune and neoplastic diseases.

Capturing change in clonal composition amongst single mouse germinal centers

Understanding cellular processes occurring in vivo on time scales of days to weeks requires repeatedly interrogating the same tissue without perturbing homeostasis. We describe a novel setup for longitudinal intravital imaging of murine peripheral lymph nodes (LNs). The formation and evolution of single germinal centers (GCs) was visualized over days to weeks. Naïve B cells encounter antigen and form primary foci, which subsequently seed GCs. These experience widely varying rates of homogenizing selection, even within closely confined spatial proximity. The fluidity of GCs is greater than previously observed with large shifts in clonality over short time scales; and loss of GCs is a rare, observable event. The observation of contemporaneous, congruent shifts in clonal composition between GCs within the same animal suggests inter-GC trafficking of memory B cells. This tool refines approaches to resolving immune dynamics in peripheral LNs with high temporospatial resolution and minimal perturbation of homeostasis.

Cell Cycle-Associated CXCR4 Expression in Germinal Center B Cells and Its Implications on Affinity Maturation

Adaptation of antibody-mediated immunity occurs in germinal centers (GC). It is where affinity maturation, class switching, memory and plasma cell differentiation synergize to generate specific high-affinity antibodies that aid both to clear and protect against reinfection of invading pathogens. Within GCs, light and dark zone are two compartments instrumental in regulating this process, by segregating T cell-dependent selection and differentiation from generation of GC B cells bearing hypermutated antigen receptors. Spatial segregation of GC B cells into the two zones relies on the chemokine receptor CXCR4, with textbooks attributing high and low expression to a dark and light zone phenotype. Interestingly, this bipolarity is not reflected in the CXCR4 expression profile of GC B cells, which is highly variable and unimodal, indicating a continuum of intermediate CXCR4 levels rather than a binary dark or light zone phenotype. Here, analysis of published BrdU pulse-chase data reveals that throughout cell cycle, average CXCR4 expression in GC B cells steadily increases close to twofold, scaling with cell surface area. CXCR4 expression in recently divided GC B cells in G0/G1 or early S phase shows intermediate levels compared to cells in G2M phase, consistent with their smaller size. The lowest number of CXCR4 receptors are displayed by relatively aged GC B cells in G0/G1 or early S phase. The latter, upon progressing through S phase, however, ramp up relative CXCR4 expression twice as much as recently divided cells. Twelve hours after the BrdU pulse, labeled GC B cells, while initially in S phase, are desynchronized in terms of cell cycle and match the CXCR4 profile of unlabeled cells. A model is discussed in which CXCR4 expression in GC B cell increases with cell cycle and cell surface area, with highest levels in G2 and M phase, coinciding with GC B cell receptor signaling in G2 and immediately preceding activation-induced cytidine deaminase (AID) activity in early G1. In the model, GC B cells compete for CXCL12 expression on the basis of their CXCR4 expression, gaining a relative advantage as they progress in cell cycle, but loosing the advantage at the moment they divide.

Third-Kind Encounters in Biomedicine: Immunology Meets Mathematics and Informatics to Become Quantitative and Predictive

The understanding of the immune response is right now at the center of biomedical research. There are growing expectations that immune-based interventions will in the midterm provide new, personalized, and targeted therapeutic options for many severe and highly prevalent diseases, from aggressive cancers to infectious and autoimmune diseases. To this end, immunology should surpass its current descriptive and phenomenological nature, and become quantitative, and thereby predictive.Immunology is an ideal field for deploying the tools, methodologies, and philosophy of systems biology, an approach that combines quantitative experimental data, computational biology, and mathematical modeling. This is because, from an organism-wide perspective, the immunity is a biological system of systems, a paradigmatic instance of a multi-scale system. At the molecular scale, the critical phenotypic responses of immune cells are governed by large biochemical networks, enriched in nested regulatory motifs such as feedback and feedforward loops. This network complexity confers them the ability of highly nonlinear behavior, including remarkable examples of homeostasis, ultra-sensitivity, hysteresis, and bistability. Moving from the cellular level, different immune cell populations communicate with each other by direct physical contact or receiving and secreting signaling molecules such as cytokines. Moreover, the interaction of the immune system with its potential targets (e.g., pathogens or tumor cells) is far from simple, as it involves a number of attack and counterattack mechanisms that ultimately constitute a tightly regulated multi-feedback loop system. From a more practical perspective, this leads to the consequence that today's immunologists are facing an ever-increasing challenge of integrating massive quantities from multi-platforms.In this chapter, we support the idea that the analysis of the immune system demands the use of systems-level approaches to ensure the success in the search for more effective and personalized immune-based therapies.

Phenotypic and Morphological Properties of Germinal Center Dark Zone Cxcl12-Expressing Reticular Cells

The germinal center (GC) is divided into a dark zone (DZ) and a light zone (LZ). GC B cells must cycle between these zones to achieve efficient Ab affinity maturation. Follicular dendritic cells (FDCs) are well characterized for their role in supporting B cell Ag encounter in primary follicles and in the GC LZ. However, the properties of stromal cells supporting B cells in the DZ are relatively unexplored. Recent work identified a novel stromal population of Cxcl12-expressing reticular cells (CRCs) in murine GC DZs. In this article, we report that CRCs have diverse morphologies, appearing in open and closed networks, with variable distribution in lymphoid tissue GCs. CRCs are also present in splenic and peripheral lymph node primary follicles. Real-time two-photon microscopy of Peyer's patch GCs demonstrates B cells moving in close association with CRC processes. CRCs are gp38(+) with low to undetectable expression of FDC markers, but CRC-like cells in the DZ are lineage marked, along with FDCs and fibroblastic reticular cells, by CD21-Cre- and Ccl19-Cre-directed fluorescent reporters. In contrast to FDCs, CRCs do not demonstrate dependence on lymphotoxin or TNF for chemokine expression or network morphology. CRC distribution in the DZ does require CXCR4 signaling, which is necessary for GC B cells to access the DZ and likely to interact with CRC processes. Our findings establish CRCs as a major stromal cell type in the GC DZ and suggest that CRCs support critical activities of GC B cells in the DZ niche through Cxcl12 expression and direct cell-cell interactions.

IL-22 regulates lymphoid chemokine production and assembly of tertiary lymphoid organs

The series of events leading to tertiary lymphoid organ (TLO) formation in mucosal organs following tissue damage remain unclear. Using a virus-induced model of autoantibody formation in the salivary glands of adult mice, we demonstrate that IL-22 provides a mechanistic link between mucosal infection, B-cell recruitment, and humoral autoimmunity. IL-22 receptor engagement is necessary and sufficient to promote differential expression of chemokine (C-X-C motif) ligand 12 and chemokine (C-X-C motif) ligand 13 in epithelial and fibroblastic stromal cells that, in turn, is pivotal for B-cell recruitment and organization of the TLOs. Accordingly, genetic and therapeutic blockade of IL-22 impairs and reverses TLO formation and autoantibody production. Our work highlights a critical role for IL-22 in TLO-induced pathology and provides a rationale for the use of IL-22-blocking agents in B-cell-mediated autoimmune conditions.

Stabilizing 3D in vivo intravital microscopy images with an iteratively refined soft-tissue model for immunology experiments

We describe a set of new algorithms and a software tool, StabiTissue, for stabilizing in vivo intravital microscopy images that suffer from soft-tissue background movement. Because these images lack predetermined anchors and are dominated by noise, we use a pixel weighted image alignment together with a correction for nonlinear tissue deformations. We call this correction a poor man׳s diffeomorphic map since it ascertains the nonlinear regions of the image without resorting to a full integral equation method. To determine the quality of the image stabilization, we developed an ensemble sampling method that quantifies the coincidence between image pairs from randomly distributed image regions. We obtain global stabilization alignment through an iterative constrained simulated annealing optimization procedure. To show the accuracy of our algorithm with existing software, we measured the misalignment error rate in datasets taken from two different organs and compared the results to a similar and popular open-source solution. Present open-source stabilization software tools perform poorly because they do not treat the specific needs of the IV-2pM datasets with soft-tissue deformation, speckle noise, full 5D inter- and intra-stack motion error correction, and undefined anchors. In contrast, the results of our tests demonstrate that our method is more immune to noise and provides better performance for datasets' possessing nonlinear tissue deformations. As a practical application of our software, we show how our stabilization improves cell tracking, where the presence of background movement would degrade track information. We also provide a qualitative comparison of our software with other open-source libraries/applications. Our software is freely available at the open source repository http://sourceforge.net/projects/stabitissue/.

Dendritic cells: key players in human herpesvirus 8 infection and pathogenesis

Human herpesvirus 8 (HHV-8; Kaposi's sarcoma-associated herpesvirus) is an oncogenic gammaherpesvirus that primarily infects cells of the immune and vascular systems. HHV-8 interacts with and targets professional antigen presenting cells and influences their function. Infection alters the maturation, antigen presentation, and immune activation capabilities of certain dendritic cells (DC) despite non-robust lytic replication in these cells. DC sustains a low level of antiviral functionality during HHV-8 infection in vitro. This may explain the ability of healthy individuals to effectively control this virus without disease. Following an immune compromising event, such as organ transplantation or human immunodeficiency virus type 1 infection, a reduced cellular antiviral response against HHV-8 compounded with skewed DC cytokine production and antigen presentation likely contributes to the development of HHV-8 associated diseases, i.e., Kaposi's sarcoma and certain B cell lymphomas. In this review we focus on the role of DC in the establishment of HHV-8 primary and latent infection, the functional state of DC during HHV-8 infection, and the current understanding of the factors influencing virus-DC interactions in the context of HHV-8-associated disease.

Automated characterization and parameter-free classification of cell tracks based on local migration behavior

Cell migration is the driving force behind the dynamics of many diverse biological processes. Even though microscopy experiments are routinely performed today by which populations of cells are visualized in space and time, valuable information contained in image data is often disregarded because statistical analyses are performed at the level of cell populations rather than at the single-cell level. Image-based systems biology is a modern approach that aims at quantitatively analyzing and modeling biological processes by developing novel strategies and tools for the interpretation of image data. In this study, we take first steps towards a fully automated characterization and parameter-free classification of cell track data that can be generally applied to tracked objects as obtained from image data. The requirements to achieve this aim include: (i) combination of different measures for single cell tracks, such as the confinement ratio and the asphericity of the track volume, and (ii) computation of these measures in a staggered fashion to retrieve local information from all possible combinations of track segments. We demonstrate for a population of synthetic cell tracks as well as for in vitro neutrophil tracks obtained from microscopy experiment that the information contained in the track data is fully exploited in this way and does not require any prior knowledge, which keeps the analysis unbiased and general. The identification of cells that show the same type of migration behavior within the population of all cells is achieved via agglomerative hierarchical clustering of cell tracks in the parameter space of the staggered measures. The recognition of characteristic patterns is highly desired to advance our knowledge about the dynamics of biological processes.

IL-17RA is essential for optimal localization of follicular Th cells in the germinal center light zone to promote autoantibody-producing B cells

Germinal centers (GCs) provide a microenvironment that promotes and regulates the interactions of B cells with follicular Th (TFH) cells. In this study, we show that there are significantly higher frequencies of CXCR5(+)ICOS(+) TFH cells in autoimmune BXD2 mice, and these cells express both IL-21R and IL-17RA. Although IL-17 and IL-21 are both important for the formation of spontaneous GCs and development of pathogenic autoantibodies, IL-21, but not IL-17, is required for the proper development of TFH cells in BXD2 mice. The total numbers of TFH cells and their ability to induce B cell responses in vitro were not affected by a deficiency of IL-17RA in BXD2-Il17ra(-/-) mice, the majority of CXCR5(+) TFH cells from BXD2-Il17ra(-/-) mice were, however, not localized in the GC light zone (LZ). Interruption of IL-17 signaling, either acutely by AdIL-17R:Fc or chronically by Il17ra(-/-), disrupted TFH-B interactions and abrogated the generation of autoantibody-forming B cells in BXD2 mice. IL-17 upregulated the expression of regulator of G-protein signaling 16 (RGS16) to promote the ability of TFH to form conjugates with B cells, which was abolished in TFH cells from BXD2-Rgs16(-/-) mice. The results suggests that IL-17 is an extrinsic stop signal that it acts on postdifferentiated IL-17RA(+) TFH to enable its interaction with responder B cells in the LZ niche. These data suggest a novel concept that TFH differentiation and its stabilization in the LZ are two separate checkpoints and that IL-21 and IL-17 act at each checkpoint to enable pathogenic GC development.

B-lymphocyte lineage cells and the respiratory system

Adaptive humoral immune responses in the airways are mediated by B cells and plasma cells that express highly evolved and specific receptors and produce immunoglobulins of most isotypes. In some cases, such as autoimmune diseases or inflammatory diseases caused by excessive exposure to foreign antigens, these same immune cells can cause disease by virtue of overly vigorous responses. This review discusses the generation, differentiation, signaling, activation, and recruitment pathways of B cells and plasma cells, with special emphasis on unique characteristics of subsets of these cells functioning within the respiratory system. The primary sensitization events that generate B cells responsible for effector responses throughout the airways usually occur in the upper airways, tonsils, and adenoid structures that make up the Waldeyer ring. On secondary exposure to antigen in the airways, antigen-processing dendritic cells migrate into secondary lymphoid organs, such as lymph nodes, that drain the upper and lower airways, and further B-cell expansion takes place at those sites. Antigen exposure in the upper or lower airways can also drive expansion of B-lineage cells in the airway mucosal tissue and lead to the formation of inducible lymphoid follicles or aggregates that can mediate local immunity or disease.

A one-shot germinal center model under protein structural stability constraints

The germinal center reaction is the process by which low-affinity B cells evolve into potent, immunoglobulin-secreting plasma and memory B cells. Since the recycling hypothesis was created, experimental studies have both tracked movement of a small population of B cells from the light zone into the dark zone, supporting the recycling model, and parallel to the light zone-dark zone interface, indicating a one-way trajectory. We present a novel, sequence-based ab initio model of protein stability and protein interactions. Our model contains a dark zone region of clonal expansion and somatic hypermutation and a light zone site of antigenic selection. We show not only that a one-shot model is sufficient to achieve biologically-realistic rates of affinity growth, population dynamics, and silent:non-silent mutation ratios in the complementary determining region and framework region of antibodies, but also that a stochastic recycling program with or without realistic constraints on the structural stabilities of GC antibodies cannot produce biologically-observed affinity growth, population dynamics or silent:non-silent mutation profiles. The effect of recycling erases affinity gains made by potent antibodies cycling back from the light zone and causes B cells to pool in the dark zone under high replication rates.

A global "imaging'' view on systems approaches in immunology

The immune system exhibits an enormous complexity. High throughput methods such as the "-omic'' technologies generate vast amounts of data that facilitate dissection of immunological processes at ever finer resolution. Using high-resolution data-driven systems analysis, causal relationships between complex molecular processes and particular immunological phenotypes can be constructed. However, processes in tissues, organs, and the organism itself (so-called higher level processes) also control and regulate the molecular (lower level) processes. Reverse systems engineering approaches, which focus on the examination of the structure, dynamics and control of the immune system, can help to understand the construction principles of the immune system. Such integrative mechanistic models can properly describe, explain, and predict the behavior of the immune system in health and disease by combining both higher and lower level processes. Moving from molecular and cellular levels to a multiscale systems understanding requires the development of methodologies that integrate data from different biological levels into multiscale mechanistic models. In particular, 3D imaging techniques and 4D modeling of the spatiotemporal dynamics of immune processes within lymphoid tissues are central for such integrative approaches. Both dynamic and global organ imaging technologies will be instrumental in facilitating comprehensive multiscale systems immunology analyses as discussed in this review.

A model for migratory B cell oscillations from receptor down-regulation induced by external chemokine fields

A long-standing paradigm in B cell immunology is that effective somatic hypermutation and affinity maturation require cycling between the dark zone and light zone of the germinal center. The cyclic re-entry hypothesis was first proposed based on considerations of the efficiency of affinity maturation using an ordinary differential equations model for B cell population dynamics. More recently, two-photon microscopy studies of B cell motility within lymph nodes in situ have revealed the complex migration patterns of B lymphocytes both in the preactivation follicle and post-activation germinal center. There is strong evidence that chemokines secreted by stromal cells and the regulation of cognate G-protein coupled receptors by these chemokines are necessary for the observed spatial cell distributions. For example, the distribution of B cells within the light and dark zones of the germinal center appears to be determined by the reciprocal interaction between the level of the CXCR4 and CXCR5 receptors and the spatial distribution of their respective chemokines CXCL12 and CXCL13. Computer simulations of individual-based models have been used to study the complex biophysical and mechanistic processes at the individual cell level, but such simulations can be challenging to parameterize and analyze. In contrast, ordinary differential equations are more tractable, but traditional compartment model formalizations ignore the spatial chemokine distribution that drives B cell redistribution. Motivated by the desire to understand the motility patterns observed in an individual-based simulation of B cell migration in the lymph node, we propose and analyze the dynamics of an ordinary differential equation model incorporating explicit chemokine spatial distributions. While there is experimental evidence that B cell migration patterns in the germinal center are driven by extrinsically regulated differentiation programs, the model shows, perhaps surprisingly, that feedback from receptor down-regulation induced by external chemokine fields can give rise to spontaneous interzonal and intrazonal oscillations in the absence of any extrinsic regulation. While the extent to which such simple feedback mechanisms contributes to B cell migration patterns in the germinal center is unknown, the model provides an alternative hypothesis for how complex B cell migration patterns might arise from very simple mechanisms.

Expression of N-cadherin in salivary gland tumors

OBJECTIVE

To detect immunohistochemically the N-cadherin expression in different types of benign and malignant salivary gland tumors in an attempt to note any possible correlation to their development, stage and invasive properties.

MATERIALS AND METHODS

N-cadherin expression was examined in tissue specimens from 49 salivary gland tumors including: pleomorphic adenomas (4), Warthin's tumors (10), and myoepitheliomas (4) (benign tumors), as well as adenoid cystic carcinomas (14), mucoepidermoid carcinomas (4), polymorphous low-grade adenocarcinomas (6), and adenocarcinomas not otherwise specified (5) (malignant tumors). Twelve specimens of normal salivary glands were used as control. The perineural invasion and stage of malignant salivary gland tumors were evaluated. Immunohistochemical procedure was performed automatically using the Bond Polymer Refine Detection Kit.

RESULTS

N-cadherin expression was not found in normal salivary glands. In benign salivary gland tumors, N-cadherin along membranes of neoplastic cells as well as in centrocytes of lymphoid germinal centers was seen in 1 and 4 cases of Warthin's tumors, respectively. Varied degree of N-cadherin expression was found in 13 (45%) cases of malignant salivary gland tumors. N-cadherin expression was significantly correlated with perineural invasion (χ(2) = 11.7, p < 0.0001), but not with stage of malignant salivary gland tumors.

CONCLUSION

N-cadherin expression was observed in malignant salivary gland tumors and could be an indicator of potentially aggressive behavior. N-cadherin expression by tumor cells could be attributed to perineural invasion.

Understanding natural killer cell regulation by mathematical approaches

The activity of natural killer (NK) cells is regulated by various processes including education/licensing, priming, integration of positive and negative signals through an array of activating and inhibitory receptors, and the development of memory-like functionality. These processes are often very complex due to the large number of different receptors and signaling pathways involved. Understanding these complex mechanisms is therefore a challenge, but is critical for understanding NK cell regulation. Mathematical approaches can facilitate the analysis and understanding of complex systems. Therefore, they may be instrumental for studies in NK cell biology. Here we provide a review of the different mathematical approaches to the analysis of NK cell signal integration, activation, proliferation, and the acquisition of inhibitory receptors. These studies show how mathematical methods can aid the analysis of NK cell regulation.

Chemotactic migration of T cells towards dendritic cells promotes the detection of rare antigens

In many immunological processes chemoattraction is thought to play a role in guiding cells to their sites of action. However, based on in vivo two-photon microscopy experiments in the absence of cognate antigen, T cell migration in lymph nodes (LNs) has been roughly described as a random walk. Although it has been shown that dendritic cells (DCs) carrying cognate antigen in some circumstances attract T cells chemotactically, it is currently still unclear whether chemoattraction of T cells towards DCs helps or hampers scanning. Chemoattraction towards DCs could on the one hand help T cells to rapidly find DCs. On the other hand, it could be deleterious if DCs become shielded by a multitude of attracted yet non-specific T cells. Results from a recent simulation study suggested that the deleterious effect dominates. We re-addressed the question whether T cell chemoattraction towards DCs is expected to promote or hamper the detection of rare antigens using the Cellular Potts Model, a formalism that allows for dynamic, flexible cellular shapes and cell migration. Our simulations show that chemoattraction of T cells enhances the DC scanning efficiency, leading to an increased probability that rare antigen-specific T cells find DCs carrying cognate antigen. Desensitization of T cells after contact with a DC further improves the scanning efficiency, yielding an almost threefold enhancement compared to random migration. Moreover, the chemotaxis-driven migration still roughly appears as a random walk, hence fine-tuned analysis of cell tracks will be required to detect chemotaxis within microscopy data.

Models and methods for analysis of lymphocyte repertoire generation, development, selection and evolution

T and B cell receptor repertoires are diversified by variable region gene rearrangement and selected based on functionality and lack of self-reactivity. Repertoires can also be defined based on phenotype and function rather than receptor specificity - such as the diversity of T helper cell subsets. Natural killer (NK) cell repertoires, in which each cell expresses a randomly chosen subset of its inhibitory receptor genes, and is educated based on self-MHC recognition by yet unknown mechanisms, are also phenotypic repertoires. Studying the generation, development and selection of lymphocyte repertoires, and their functions during immune responses, is essential for understanding the function of the immune system in healthy individuals and in immune deficient, autoimmune or cancer patients. The study of lymphocyte repertoires will enable clinical immunologists to develop better therapeutic monoclonal antibodies, vaccines, transplantation donor-recipient matching protocols, and other immune intervention strategies. The recent development of high-throughput methods for repertoire data collection - from multicolor flow cytometry through single-cell imaging to deep sequencing - presents us now, for the first time, with the ability to analyze and compare large samples of lymphocyte repertoires in health, aging and disease. The exponential growth of these datasets, however, challenges the theoretical immunology community to develop methods for data organization and analysis. Furthermore, the need to test hypotheses regarding immune function, and generate predictions regarding the outcomes of medical interventions, necessitates the development of complex mathematical and computational models, covering processes on multiple scales, from the genetic and molecular to the cellular and system scales.

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.

S1PR2 links germinal center confinement and growth regulation

Germinal centers (GCs) are sites of rapid B-cell proliferation and somatic mutation. These ovoid structures develop within the center of follicles and grow to a stereotypic size. The cell migration and interaction dynamics underlying GC B-cell selection events are currently under intense scrutiny. In recent study, we identified a role for a migration inhibitory receptor, S1PR2, in promoting GC B-cell confinement to GCs. S1PR2 also dampens Akt activation and deficiency in S1PR2 or components of its signaling pathway result in a loss of growth control in chronically stimulated mucosal GCs. Herein, we detail present understanding of S1PR2 and S1P biology as it pertains to GC B cells and place this information in the context of a current model of GC function.

CX3CL1/fractalkine is a novel regulator of normal and malignant human B cell function

CX(3)CL1, or fractalkine, the unique member of the CX(3)C chemokine family, exists as a transmembrane glycoprotein, as well as in soluble form, each mediating different biological activities, and is constitutively expressed in many hematopoietic and nonhematopoietic tissues. CX(3)CR1, the CX(3)CL1 exclusive receptor, is a classical GPCR, expressed on NK cells, CD14(+) monocytes, and some subpopulation of T cells, B cells, and mast cells. A recent paper by our group has demonstrated for the first time that highly purified human B cells from tonsil and peripheral blood expressed CX(3)CR1 at mRNA and protein levels. In particular, tonsil naïve, GC, and memory B cells expressed CX(3)CR1, but only GC centrocytes were attracted by soluble CX(3)CL1, which with its receptor, are also involved in the pathogenesis of several inflammatory disorders, as well as of cancer. Previous studies have shown that CX(3)CR1 is up-regulated in different types of B cell lymphoma, as well as in B-CLL. Recently, we have demonstrated that the CX(3)CL1/CX(3)CR1 axis is involved in the interaction of B-CLL cells with their microenvironment. Taken together, our data delineate a novel role for the CX(3)CL1/CX(3)CR1 complex in the biology of normal B cells and B-CLL cells. These topics are the subject of this review article.

Molecular programming of B cell memory

The development of high-affinity B cell memory is regulated through three separable phases, each involving antigen recognition by specific B cells and cognate T helper cells. Initially, antigen-primed B cells require cognate T cell help to gain entry into the germinal centre pathway to memory. Once in the germinal centre, B cells with variant B cell receptors must access antigens and present them to germinal centre T helper cells to enter long-lived memory B cell compartments. Following antigen recall, memory B cells require T cell help to proliferate and differentiate into plasma cells. A recent surge of information - resulting from dynamic B cell imaging in vivo and the elucidation of T follicular helper cell programmes - has reshaped the conceptual landscape surrounding the generation of memory B cells. In this Review, we integrate this new information about each phase of antigen-specific B cell development to describe the newly unravelled molecular dynamics of memory B cell programming.

Defining the quantitative limits of intravital two-photon lymphocyte tracking

Two-photon microscopy has substantially advanced our understanding of cellular dynamics in the immune system. Cell migration can now be imaged in real time in the living animal. Strikingly, the migration of naive lymphocytes in secondary lymphoid tissue appears predominantly random. It is unclear, however, whether directed migration may escape detection in this random background. Using a combination of mathematical modeling and experimental data, we investigate the extent to which modern two-photon imaging can rule out biologically relevant directed migration. For naive T cells migrating in uninfected lymph nodes (LNs) at average 3D speeds of around 18 μm/min, we rule out uniform directed migration of more than 1.7 μm/min at the 95% confidence level, confirming that T cell migration is indeed mostly random on a timescale of minutes. To investigate whether this finding still holds for longer timescales, we use a 3D simulation of the naive T cell LN transit. A pure random walk predicts a transit time of around 16 h, which is in good agreement with experimental results. A directional bias of only 0.5 μm/min-less than 3% of the cell speed-would already accelerate the transit twofold. These results jointly strengthen the random walk analogy for naive T cell migration in LNs, but they also emphasize that very small deviations from random migration can still be important. Our methods are applicable to cells of any type and can be used to reanalyze existing datasets.