Estimating Residence Times of Lymphocytes in Ovine Lymph Nodes

The ability of lymphocytes to recirculate between blood and secondary lymphoid tissues such as lymph nodes (LNs) and spleen is well established. Sheep have been used as an experimental system to study lymphocyte recirculation for decades and multiple studies document accumulation and loss of intravenously (i.v.) transferred lymphocytes in efferent lymph of various ovine LNs. Yet, surprisingly little work has been done to accurately quantify the dynamics of lymphocyte exit from the LNs and to estimate the average residence times of lymphocytes in ovine LNs. In this work we developed a series of mathematical models based on fundamental principles of lymphocyte recirculation in the body under non-inflammatory (resting) conditions. Our analysis suggested that in sheep, recirculating lymphocytes spend on average 3 h in the spleen and 20 h in skin or gut-draining LNs with a distribution of residence times in LNs following a skewed gamma (lognormal-like) distribution. Our mathematical models also suggested an explanation for a puzzling observation of the long-term persistence of i.v. transferred lymphocytes in the efferent lymph of the prescapular LN (pLN); the model predicted that this is a natural consequence of long-term persistence of the transferred lymphocytes in circulation. We also found that lymphocytes isolated from the skin-draining pLN have a 2-fold increased entry rate into the pLN as opposed to the mesenteric (gut-draining) LN (mLN). Likewise, lymphocytes from mLN had a 3-fold increased entry rate into the mLN as opposed to entry rate into pLN. In contrast, these cannulation data could not be explained by preferential retention of cells in LNs of their origin. Taken together, our work illustrates the power of mathematical modeling in describing the kinetics of lymphocyte migration in sheep and provides quantitative estimates of lymphocyte residence times in ovine LNs

Mathematical Modeling Suggests Cooperation of Plant-Infecting Viruses

Viruses are major pathogens of agricultural crops. Viral infections often start after the virus enters the outer layer of a tissue, and many successful viruses, after local replication in the infected tissue, are able to spread systemically. Quantitative details of virus dynamics in plants, however, are poorly understood, in part, because of the lack of experimental methods which allow the accurate measurement of the degree of infection in individual plant tissues. Recently, a group of researchers followed the kinetics of infection of individual cells in leaves of Nicotiana tabacum plants using Tobacco etch virus (TEV) expressing either Venus or blue fluorescent protein (BFP). Assuming that viral spread occurs from lower to upper leaves, the authors fitted a simple mathematical model to the frequency of cellular infection by the two viral variants found using flow cytometry. While the original model could accurately describe the kinetics of viral spread locally and systemically, we found that many alternative versions of the model, for example, if viral spread starts at upper leaves and progresses to lower leaves or when virus dissemination is stopped due to an immune response, fit the data with reasonable quality, and yet with different parameter estimates. These results strongly suggest that experimental measurements of the virus infection in individual leaves may not be sufficient to identify the pathways of viral dissemination between different leaves and reasons for viral control. We propose experiments that may allow discrimination between the alternatives. By analyzing the kinetics of coinfection of individual cells by Venus and BFP strains of TEV we found a strong deviation from the random infection model, suggesting cooperation between the two strains when infecting plant cells. Importantly, we showed that many mathematical models on the kinetics of coinfection of cells with two strains could not adequately describe the data, and the best fit model needed to assume (i) different susceptibility of uninfected cells to infection by two viruses locally in the leaf vs. systemically from other leaves, and (ii) decrease in the infection rate depending on the fraction of uninfected cells which could be due to a systemic immune response. Our results thus demonstrate the difficulty in reaching definite conclusions from extensive and yet limited experimental data and provide evidence of potential cooperation between different viral variants infecting individual cells in plants

Spatial heterogeneity and peptide availability determine CTL killing efficiency in vivo

The rate at which a cytotoxic T lymphocyte (CTL) can survey for infected cells is a key ingredient of models of vertebrate immune responses to intracellular pathogens. Estimates have been obtained using in vivo cytotoxicity assays in which peptide-pulsed splenocytes are killed by CTL in the spleens of immunised mice. However the spleen is a heterogeneous environment and splenocytes comprise multiple cell types. Are some cell types intrinsically more susceptible to lysis than others? Quantitatively, what impacts are made by the spatial distribution of targets and effectors, and the level of peptide-MHC on the target cell surface? To address these questions we revisited the splenocyte killing assay, using CTL specific for an epitope of influenza virus. We found that at the cell population level T cell targets were killed more rapidly than B cells. Using modeling, quantitative imaging and in vitro killing assays we conclude that this difference in vivo likely reflects different migratory patterns of targets within the spleen and a heterogeneous distribution of CTL, with no detectable difference in the intrinsic susceptibilities of the two populations to lysis. Modeling of the stages involved in the detection and killing of peptide-pulsed targets in vitro revealed that peptide dose influenced the ability of CTL to form conjugates with targets but had no detectable effect on the probability that conjugation resulted in lysis, and that T cell targets took longer to lyse than B cells. We also infer that incomplete killing in vivo of cells pulsed with low doses of peptide may be due to a combination of heterogeneity in peptide uptake and the dissociation, but not internalisation, of peptide-MHC complexes. Our analyses demonstrate how population-averaged parameters in models of immune responses can be dissected to account for both spatial and cellular heterogeneity

Evaluating contribution of the cellular and humoral immune responses to the control of shedding of \u3cem\u3eMycobacterium avium\u3c/em\u3e spp. \u3cem\u3eparatuberculosis\u3c/em\u3e in cattle

Mycobacterium avium spp. paratuberculosis (MAP) causes a persistent infection and chronic inflammation of the gut in ruminants leading to bacterial shedding in feces in many infected animals. Although there are often strong MAP-specific immune responses in infected animals, immunological correlates of protection against progression to disease remain poorly defined. Analysis of cross-sectional data has suggested that the cellular immune response observed early in infection is effective at containing bacterial growth and shedding, in contrast to humoral immune responses. In this study, 20 MAP-infected calves were followed for nearly 5 years during which MAP shedding, antigen-specific cellular (LPT) and humoral (ELISA) immune responses were measured. We found that MAP-specific cellular immune response developed slowly, with the peak of the immune response occurring one year post infection. MAP-specific humoral immunity expanded only in a subset of animals. Only in a subset of animals there was a statistically significant negative correlation between the amount of MAP shedding and magnitude of the MAP-specific cellular immune response. Direct fitting of simple mechanistic mathematical models to the shedding data suggested that MAP-specific immune responses contributed significantly to the kinetics of MAP shedding in most animals. However, whereas the MAP-specific cellular immune response was predicted to suppress shedding in some animals, in other animals it was predicted to increase shedding. In contrast, MAP-specific humoral response was always predicted to increase shedding. Our results illustrate the use of mathematical methods to understand relationships between mycobacteria and immunity in vivo but also highlight problems with establishing cause-effect links from observational data

Clustering of activated CD8 T cells around Malaria-infected hepatocytes is rapid and is driven by antigen-specific cells

Malaria, a disease caused by parasites of the Plasmodium genus, begins when Plasmodium-infected mosquitoes inject malaria sporozoites while searching for blood. Sporozoites migrate from the skin via blood to the liver, infect hepatocytes, and form liver stages which in mice 48 h later escape into blood and cause clinical malaria. Vaccine-induced activated or memory CD8 T cells are capable of locating and eliminating all liver stages in 48 h, thus preventing the blood-stage disease. However, the rules of how CD8 T cells are able to locate all liver stages within a relatively short time period remains poorly understood. We recently reported formation of clusters consisting of variable numbers of activated CD8 T cells around Plasmodium yoelii (Py)-infected hepatocytes. Using a combination of experimental data and mathematical models we now provide additional insights into mechanisms of formation of these clusters. First, we show that a model in which cluster formation is driven exclusively by T-cell-extrinsic factors, such as variability in “attractiveness” of different liver stages, cannot explain distribution of cluster sizes in different experimental conditions. In contrast, the model in which cluster formation is driven by the positive feedback loop (i.e., larger clusters attract more CD8 T cells) can accurately explain the available data. Second, while both Py-specific CD8 T cells and T cells of irrelevant specificity (non-specific CD8 T cells) are attracted to the clusters, we found no evidence that non-specific CD8 T cells play a role in cluster formation. Third and finally, mathematical modeling suggested that formation of clusters occurs rapidly, within few hours after adoptive transfer of CD8 T cells, thus illustrating high efficiency of CD8 T cells in locating their targets in complex peripheral organs, such as the liver. Taken together, our analysis provides novel insights into and attempts to discriminate between alternative mechanisms driving the formation of clusters of antigen-specific CD8 T cells in the liver

Cellular and population plasticity of helper CD4(+) T cell responses

Vertebrates are constantly exposed to pathogens, and the adaptive immunity has most likely evolved to control and clear such infectious agents. CD4(+) T cells are the major players in the adaptive immune response to pathogens. Following recognition of pathogen-derived antigens naïve CD4(+) T cells differentiate into effectors which then control pathogen replication either directly by killing pathogen-infected cells or by assisting with generation of cytotoxic T lymphocytes (CTLs) or pathogen-specific antibodies. Pathogen-specific effector CD4(+) T cells are highly heterogeneous in terms of cytokines they produce. Three major subtypes of effector CD4(+) T cells have been identified: T-helper 1 (Th1) cells producing IFN-γ and TNF-α, Th2 cells producing IL-4 and IL-10, and Th17 cells producing IL-17. How this heterogeneity is maintained and what regulates changes in effector T cell composition during chronic infections remains poorly understood. In this review we discuss recent advances in our understanding of CD4(+) T cell differentiation in response to microbial infections. We propose that a change in the phenotype of pathogen-specific effector CD4(+) T cells during chronic infections, for example, from Th1 to Th2 response as observed in Mycobactrium avium ssp. paratuberculosis (MAP) infection of ruminants, can be achieved by conversion of T cells from one effector subset to another (cellular plasticity) or due to differences in kinetics (differentiation, proliferation, death) of different effector T cell subsets (population plasticity). We also shortly review mathematical models aimed at describing CD4(+) T cell differentiation and outline areas for future experimental and theoretical research. doi: 10.3389/fphys.2013.0020

Estimating Costs and Benefits of CTL Escape Mutations in SIV/HIV Infection

Mutations that allow SIV/HIV to avoid the cytotoxic T lymphocyte (CTL) response are well documented. Recently, there have been a few attempts at estimating the costs of CTL escape mutations in terms of the reduction in viral fitness and the killing rate at which the CTL response specific to one viral epitope clears virus-infected cells. Using a mathematical model we show that estimation of both parameters depends critically on the underlying changes in the replication rate of the virus and the changes in the killing rate over time (which in previous studies were assumed to be constant). We provide a theoretical basis for estimation of these parameters using in vivo data. In particular, we show that 1) by assuming unlimited virus growth one can obtain a minimal estimate of the fitness cost of the escape mutation, and 2) by assuming no virus growth during the escape, one can obtain a minimal estimate of the average killing rate. We also discuss the conditions under which better estimates of the average killing rate can be obtained

How Does Cross-Reactive Stimulation Affect the Longevity of CD8+ T Cell Memory?

Immunological memory—the ability to “remember” previously encountered pathogens and respond faster upon re-exposure is a central feature of the immune response in vertebrates. The cross-reactive stimulation hypothesis for the maintenance of memory proposes that memory cells specific for a given pathogen are maintained by cross-reactive stimulation following infections with other (unrelated) pathogens. We use mathematical models to examine the cross-reactive stimulation hypothesis. We find that: (i) the direct boosting of cross-reactive lineages only provides a very small increase in the average longevity of immunological memory; (ii) the expansion of cross-reactive lineages can indirectly increase the longevity of memory by reducing the magnitude of expansion of new naive lineages which occupy space in the memory compartment and are responsible for the decline in memory; (iii) cross-reactive stimulation results in variation in the rates of decline of different lineages of memory cells and enrichment of memory cell population for cells that are cross-reactive for the pathogens to which the individual has been exposed