T-cell movement on the reticular network

Quantitative analyses of T cell motion in tissue reveals factors driving T cell search in tissues

T cells are required to clear infection, and T cell motion plays a role in how quickly a T cell finds its target, from initial naive T cell activation by a dendritic cell to interaction with target cells in infected tissue. To better understand how different tissue environments affect T cell motility, we compared multiple features of T cell motion including speed, persistence, turning angle, directionality, and confinement of T cells moving in multiple murine tissues using microscopy. We quantitatively analyzed naive T cell motility within the lymph node and compared motility parameters with activated CD8 T cells moving within the villi of small intestine and lung under different activation conditions. Our motility analysis found that while the speeds and the overall displacement of T cells vary within all tissues analyzed, T cells in all tissues tended to persist at the same speed. Interestingly, we found that T cells in the lung show a marked population of T cells turning at close to 180o, while T cells in lymph nodes and villi do not exhibit this "reversing" movement. T cells in the lung also showed significantly decreased meandering ratios and increased confinement compared to T cells in lymph nodes and villi. These differences in motility patterns led to a decrease in the total volume scanned by T cells in lung compared to T cells in lymph node and villi. These results suggest that the tissue environment in which T cells move can impact the type of motility and ultimately, the efficiency of T cell search for target cells within specialized tissues such as the lung.

Assessing the impact of HIV treatment interruptions using stochastic cellular Automata

Chronic HIV infection causes a progressive decrease in the ability to maintain homeostasis resulting, after some time, in eventual break down of immune functions. Recent clinical research has shed light on a significant contribution of the lymphatic tissues, where HIV causes accumulation of collagen, (fibrosis). Specifically, where tissue is populated by certain types of functional stromal cells designated Fibroblastic Reticular Cells (FRCs), these have been found to play a crucial role in balancing out apoptosis and regeneration of naïve T-cells through 2-way cellular signaling. Tissue fibrosis not only impedes this signaling, effectively reducing T-cell levels through increased apoptosis of cells of both T- and FRC type but has been found to be irreversible by current HIV standard treatment (cART). While the therapy aims to block the viral lifecycle, cART-associated increase of T-cell levels in blood appears to conceal existing FRC impairment through fibrosis. This hidden impairment can lead to adverse consequences if treatment is interrupted, e.g. due to poor adherence (missing doses) or through periods recovering from drug toxicities. Formal clinical studies on treatment interruption have indicated possible adverse effects, but quantification of those effects in relation to interruption protocol and patient predisposition remains unclear. Accordingly, the impact of treatment interruption on lymphatic tissue structure and T-cell levels is explored here by means of computer simulation. A novel Stochastic Cellular Automata model is proposed, which utilizes all sources of clinical detail available to us (though sparse in part) for model parametrization. Sources are explicitly referenced and conflicting evidence from previous studies explored. The main focus is on (i) spatial aspects of collagen build up, together with (ii) collagen increase after repeated treatment interruptions to explore the dynamics of HIV-induced fibrosis and T-cell loss.

High-resolution 3D imaging and topological mapping of the lymph node conduit system

The conduit network is a hallmark of lymph node microanatomy, but lack of suitable imaging technology has prevented comprehensive investigation of its topology. We employed an extended-volume imaging system to capture the conduit network of an entire murine lymph node (comprising over 280,000 segments). The extensive 3D images provide a comprehensive overview of the regions supplied by conduits, including perivascular sleeves and distinctive “follicular reservoirs” within B cell follicles, surrounding follicular dendritic cells. A 3D topology map of conduits within the T-cell zone showed homogeneous branching, but conduit density was significantly higher in the superficial T-cell zone compared with the deep zone, where distances between segments are sufficient for T cells to lose contact with fibroblastic reticular cells. This topological mapping of the conduit anatomy can now aid modeling of its roles in lymph node function, as we demonstrate by simulating T-cell motility in the different T-cell zones.

Extended-volume confocal imaging allowed 3D visualisation of the fine network of conduits within lymph nodes; the resulting map of conduit topology underscores structural differences between the deep and superficial T cell zone and identifies "follicular reservoirs" within B cell follicles that concentrate lymphoid fluid around follicular dendritic cells.

Fate of a Naive T Cell: A Stochastic Journey

The homeostasis of T cell populations depends on migration, division and death of individual cells (1). T cells migrate between spatial compartments (spleen, lymph nodes, lung, liver, etc.), where they may divide or differentiate, and eventually die (2). The kinetics of recirculation influences the speed at which local infections are detected and controlled (3). New experimental techniques have been developed to measure the lifespan of cells, and their migration dynamics; for example, fluorescence-activated cell sorting (4), in vitro time-lapse microscopy (5), or in vivo stable isotope labeling (e.g., deuterium) (6). When combined with mathematical and computational models, they allow estimation of rates of migration, division, differentiation and death (6, 7). In this work, we develop a stochastic model of a single cell migrating between spatial compartments, dividing and eventually dying. We calculate the number of division events during a T cell's journey, its lifespan, the probability of dying in each compartment and the number of progeny cells. A fast-migration approximation allows us to compute these quantities when migration rates are larger than division and death rates. Making use of published rates: (i) we analyse how perturbations in a given spatial compartment impact the dynamics of a T cell, (ii) we study the accuracy of the fast-migration approximation, and (iii) we quantify the role played by direct migration (not via the blood) between some compartments.

Integrative Computational Modeling of the Lymph Node Stromal Cell Landscape

Adaptive immune responses develop in secondary lymphoid organs such as lymph nodes (LNs) in a well-coordinated series of interactions between migrating immune cells and resident stromal cells. Although many processes that occur in LNs are well understood from an immunological point of view, our understanding of the fundamental organization and mechanisms that drive these processes is still incomplete. The aim of systems biology approaches is to unravel the complexity of biological systems and describe emergent properties that arise from interactions between individual constituents of the system. The immune system is greater than the sum of its parts, as is the case with any sufficiently complex system. Here, we review recent work and developments of computational LN models with focus on the structure and organization of the stromal cells. We explore various mathematical studies of intranodal T cell motility and migration, their interactions with the LN-resident stromal cells, and computational models of functional chemokine gradient fields and lymph flow dynamics. Lastly, we discuss briefly the importance of hybrid and multi-scale modeling approaches in immunology and the technical challenges involved.

The Reticular Cell Network: A Robust Backbone for Immune Responses

Lymph nodes are meeting points for circulating immune cells. A network of reticular cells that ensheathe a mesh of collagen fibers crisscrosses the tissue in each lymph node. This reticular cell network distributes key molecules and provides a structure for immune cells to move around on. During infections, the network can suffer damage. A new study has now investigated the network's structure in detail, using methods from graph theory. The study showed that the network is remarkably robust to damage: it can still support immune responses even when half of the reticular cells are destroyed. This is a further important example of how network connectivity achieves tolerance to failure, a property shared with other important biological and nonbiological networks.

Leukocyte Motility Models Assessed through Simulation and Multi-objective Optimization-Based Model Selection

The advent of two-photon microscopy now reveals unprecedented, detailed spatio-temporal data on cellular motility and interactions in vivo. Understanding cellular motility patterns is key to gaining insight into the development and possible manipulation of the immune response. Computational simulation has become an established technique for understanding immune processes and evaluating hypotheses in the context of experimental data, and there is clear scope to integrate microscopy-informed motility dynamics. However, determining which motility model best reflects in vivo motility is non-trivial: 3D motility is an intricate process requiring several metrics to characterize. This complicates model selection and parameterization, which must be performed against several metrics simultaneously. Here we evaluate Brownian motion, Lévy walk and several correlated random walks (CRWs) against the motility dynamics of neutrophils and lymph node T cells under inflammatory conditions by simultaneously considering cellular translational and turn speeds, and meandering indices. Heterogeneous cells exhibiting a continuum of inherent translational speeds and directionalities comprise both datasets, a feature significantly improving capture of in vivo motility when simulated as a CRW. Furthermore, translational and turn speeds are inversely correlated, and the corresponding CRW simulation again improves capture of our in vivo data, albeit to a lesser extent. In contrast, Brownian motion poorly reflects our data. Lévy walk is competitive in capturing some aspects of neutrophil motility, but T cell directional persistence only, therein highlighting the importance of evaluating models against several motility metrics simultaneously. This we achieve through novel application of multi-objective optimization, wherein each model is independently implemented and then parameterized to identify optimal trade-offs in performance against each metric. The resultant Pareto fronts of optimal solutions are directly contrasted to identify models best capturing in vivo dynamics, a technique that can aid model selection more generally. Our technique robustly determines our cell populations’ motility strategies, and paves the way for simulations that incorporate accurate immune cell motility dynamics.

Advances in imaging technology allow investigators to monitor the movements and interactions of immune cells in a live animal, processes essential to understanding and manipulating how an immune response is generated. T cells in the brains of Toxoplasma gondii-infected mice have previously been described as performing a Lévy walk, an optimal strategy for locating sparsely, randomly distributed targets. Determining which motility model best characterizes a population of cells is problematic; multiple metrics are required to specify as intricate and nuanced a process as 3D motility, and the tools to evaluate model-parameter combinations have been lacking. We have developed a novel framework to perform this model evaluation through simulation, a popular tool for exploring complex immune system phenomena. We find that Lévy walk offers an inferior capture of our data to another class of motility model, the correlated random walk, and this determination was possible because we are able to explicitly evaluate several motility metrics simultaneously. Further, we find evidence that leukocytes differ in their inherent translational and rotational speeds. These findings facilitate more accurate immune system simulations aimed at unravelling the processes underpinning this critical biological function.

Immune-inspired search strategies for robot swarms

Detection of targets distributed randomly in space is a task common to both robotic and biological systems. Lévy search has previously been used to characterize T cell search in the immune system. We use a robot swarm to evaluate the effectiveness of a Lévy search strategy and map the relationship between search parameters and target configurations. We show that the fractal dimension of the Lévy search which optimizes search efficiency depends strongly on the distribution of targets but only weakly on the number of agents involved in search. Lévy search can therefore be tuned to the target configuration while also being scalable. Implementing search behaviors observed in T cells in a robot swarm provides an effective, adaptable, and scalable swarm robotic search strategy. Additionally, the adaptability and scalability of Lévy search may explain why Lévy-like movement has been observed in T cells in multiple immunological contexts.

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.

Persistence and Adaptation in Immunity: T Cells Balance the Extent and Thoroughness of Search

Effective search strategies have evolved in many biological systems, including the immune system. T cells are key effectors of the immune response, required for clearance of pathogenic infection. T cell activation requires that T cells encounter antigen-bearing dendritic cells within lymph nodes, thus, T cell search patterns within lymph nodes may be a crucial determinant of how quickly a T cell immune response can be initiated. Previous work suggests that T cell motion in the lymph node is similar to a Brownian random walk, however, no detailed analysis has definitively shown whether T cell movement is consistent with Brownian motion. Here, we provide a precise description of T cell motility in lymph nodes and a computational model that demonstrates how motility impacts T cell search efficiency. We find that both Brownian and Lévy walks fail to capture the complexity of T cell motion. Instead, T cell movement is better described as a correlated random walk with a heavy-tailed distribution of step lengths. Using computer simulations, we identify three distinct factors that contribute to increasing T cell search efficiency: 1) a lognormal distribution of step lengths, 2) motion that is directionally persistent over short time scales, and 3) heterogeneity in movement patterns. Furthermore, we show that T cells move differently in specific frequently visited locations that we call "hotspots" within lymph nodes, suggesting that T cells change their movement in response to the lymph node environment. Our results show that like foraging animals, T cells adapt to environmental cues, suggesting that adaption is a fundamental feature of biological search.

T cell and reticular network co-dependence in HIV infection

Fibroblastic reticular cells (FRC) are arranged on a network in the T cell zone of lymph nodes, forming a scaffold for T cell migration, and providing survival factors, especially interleukin-7 (IL-7). Conversely, CD4(+) T cells are the major producers of lymphotoxin-β (LT-β), necessary for the construction and maintenance of the FRC network. This interdependence creates the possibility of a vicious cycle, perpetuating loss of both FRC and T cells. Furthermore, evidence that HIV infection is responsible for collagenation of the network suggests that long term loss of network function might be responsible for the attenuated recovery in T cell count seen in HIV patients undergoing antiretroviral therapy (ART). We present computational and mathematical models of this interaction mechanism and subsequent naive CD4(+) T-cell depletion in which (1) collagen deposition impedes access of naive T cells to IL-7 on the FRC and loss of IL-7 production by loss of FRC network itself, leading to the depletion of naive T cells through increased apoptosis; and (2) depletion of naive T cells as the source of LT-β on which the FRC depend for survival leads to loss of the network, thereby amplifying and perpetuating the cycle of depletion of both naive T cells and stromal cells. Our computational model explicitly includes an FRC network and its cytokine exchange with a heterogeneous T-cell population. We also derive lumped models, in terms of partial differential equations and reduced to ordinary differential equations, that provide additional insight into the mechanisms at work. The central conclusions are that (1) damage to the reticular network, caused by HIV infection is a plausible mechanism for attenuated recovery post-ART; (2) within this, the production of T cell survival factors by FRCs may be the key rate-limiting step; and (3) the methods of model reduction and analysis presented are useful for both immunological studies and other contexts in which agent-based models are severely limited by computational cost.

Analytical results on the Beauchemin model of lymphocyte migration

The Beauchemin model is a simple particle-based description of stochastic lymphocyte migration in tissue, which has been successfully applied to studying immunological questions. In addition to being easy to implement, the model is also to a large extent mathematically tractable. This article provides a comprehensive overview of both existing and new analytical results on the Beauchemin model within a common mathematical framework. Specifically, we derive the motility coefficient, the mean square displacement, and the confinement ratio, and discuss four different methods for simulating biased migration of pre-defined speed. The results provide new insight into published studies and a reference point for future research based on this simple and popular lymphocyte migration model.

Influence of the fibroblastic reticular network on cell-cell interactions in lymphoid organs

Secondary lymphoid organs (SLO), such as lymph nodes and the spleen, display a complex micro-architecture. In the T cell zone the micro-architecture is provided by a network of fibroblastic reticular cells (FRC) and their filaments. The FRC network is thought to enhance the interaction between immune cells and their cognate antigen. However, the effect of the FRC network on cell interaction cannot be quantified to date because of limitations in immunological methodology. We use computational models to study the influence of different densities of FRC networks on the probability that two cells meet. We developed a 3D cellular automaton model to simulate cell movements and interactions along the FRC network inside lymphatic tissue. We show that the FRC network density has only a small effect on the probability of a cell to come into contact with a static or motile target. However, damage caused by a disruption of the FRC network is greatest at FRC densities corresponding to densities observed in the spleen of naïve mice. Our analysis suggests that the FRC network as a guiding structure for moving T cells has only a minor effect on the probability to find a corresponding dendritic cell. We propose alternative hypotheses by which the FRC network might influence the functionality of immune responses in a more significant way.