Estimates and impact of lymphocyte division parameters from CFSE data using mathematical modelling

Counting generations in birth and death processes with competing Erlang and exponential waiting times

Lymphocyte populations, stimulated in vitro or in vivo, grow as cells divide. Stochastic models are appropriate because some cells undergo multiple rounds of division, some die, and others of the same type in the same conditions do not divide at all. If individual cells behave independently, then each cell can be imagined as sampling from a probability density of times to division and death. The exponential density is the most mathematically and computationally convenient choice. It has the advantage of satisfying the memoryless property, consistent with a Markov process, but it overestimates the probability of short division times. With the aim of preserving the advantages of a Markovian framework while improving the representation of experimentally-observed division times, we consider a multi-stage model of cellular division and death. We use Erlang-distributed (or, more generally, phase-type distributed) times to division, and exponentially distributed times to death. We classify cells into generations, using the rule that the daughters of cells in generation n are in generation \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n+1$$\end{document}n+1. In some circumstances, our representation is equivalent to established models of lymphocyte dynamics. We find the growth rate of the cell population by calculating the proportions of cells by stage and generation. The exponent describing the late-time cell population growth, and the criterion for extinction of the population, differs from what would be expected if N steps with rate \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda$$\end{document}λ were equivalent to a single step of rate \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda /N$$\end{document}λ/N. We link with a published experimental dataset, where cell counts were reported after T cells were transferred to lymphopenic mice, using Approximate Bayesian Computation. In the comparison, the death rate is assumed to be proportional to the generation and the Erlang time to division for generation 0 is allowed to differ from that of subsequent generations. The multi-stage representation is preferred to a simple exponential in posterior distributions, and the mean time to first division is estimated to be longer than the mean time to subsequent divisions.

Cyton2: A Model of Immune Cell Population Dynamics That Includes Familial Instructional Inheritance

Lymphocytes are the central actors in adaptive immune responses. When challenged with antigen, a small number of B and T cells have a cognate receptor capable of recognising and responding to the insult. These cells proliferate, building an exponentially growing, differentiating clone army to fight off the threat, before ceasing to divide and dying over a period of weeks, leaving in their wake memory cells that are primed to rapidly respond to any repeated infection. Due to the non-linearity of lymphocyte population dynamics, mathematical models are needed to interrogate data from experimental studies. Due to lack of evidence to the contrary and appealing to arguments based on Occam’s Razor, in these models newly born progeny are typically assumed to behave independently of their predecessors. Recent experimental studies, however, challenge that assumption, making clear that there is substantial inheritance of timed fate changes from each cell by its offspring, calling for a revision to the existing mathematical modelling paradigms used for information extraction. By assessing long-term live-cell imaging of stimulated murine B and T cells in vitro, we distilled the key phenomena of these within-family inheritances and used them to develop a new mathematical model, Cyton2, that encapsulates them. We establish the model’s consistency with these newly observed fine-grained features. Two natural concerns for any model that includes familial correlations would be that it is overparameterised or computationally inefficient in data fitting, but neither is the case for Cyton2. We demonstrate Cyton2’s utility by challenging it with high-throughput flow cytometry data, which confirms the robustness of its parameter estimation as well as its ability to extract biological meaning from complex mixed stimulation experiments. Cyton2, therefore, offers an alternate mathematical model, one that is, more aligned to experimental observation, for drawing inferences on lymphocyte population dynamics.

Prediction of PD-L1 inhibition effects for HIV-infected individuals

The novel therapies with immune checkpoint inhibitors hold great promises for patients with chronic virus infections and cancers. This is based mainly on the partial reversal of the exhausted phenotype of antigen-specific cytotoxic CD8 T cells (CTL). Recently, we have shown that the restoration of HIV-specific T cell function depends on the HIV infection stage of an infected individual. Here we aimed to answer two fundamental questions: (i) Can one estimate growth parameters for the HIV-specific proliferative responsiveness upon PD-L1 blockade ex vivo? (ii) Can one use these parameter estimates to predict clinical benefit for HIV-infected individuals displaying diverse infection phenotypes? To answer these questions, we first analyzed HIV-1 Gag-specific CD8 T cell proliferation by time-resolved CFSE assays and estimated the effect of PD-L1 blockade on division and death rates, and specific precursor frequencies. These values were then incorporated into a model for CTL-mediated HIV control and the effects on CTL frequencies, viral loads and CD4 T cell counts were predicted for different infection phenotypes. The biggest absolute increase in CD4 T cell counts was in the group of slow progressors while the strongest reduction in virus loads was observed in progressor patients. These results suggest a significant clinical benefit only for a subgroup of HIV-infected individuals. However, as PD1 is a marker of lymphocyte activation and expressed on several lymphocyte subsets including also CD4 T cells and B cells, we subsequently examined the multiple effects of anti-PD-L1 blockade beyond those on CD8 T cells. This extended model then predicts that the net effect on HIV load and CD4 T cell number depends on the interplay between positive and negative effects of lymphocyte subset activation. For a physiologically relevant range of affected model parameters, PD-L1 blockade is likely to be overall beneficial for HIV-infected individuals.