Publications

Quantitative understanding of the kinetics of lymphocyte proliferation and death upon activation with an antigen is crucial for elucidating factors determining the magnitude, duration and efficiency of the immune response. Recent advances in quantitative experimental techniques, in particular intracellular labeling and multi-channel flow cytometry, allow one to measure the population structure of proliferating and dying lymphocytes for several generations with high precision. These new experimental techniques require novel quantitative methods of analysis. We review several recent mathematical approaches used to describe and analyze cell proliferation data. Using a rigorous mathematical framework, we show that two commonly used models that are based on the theories of age-structured cell populations and of branching processes, are mathematically identical. We provide several simple analytical solutions for a model in which the distribution of inter-division times follows a gamma distribution and show that this model can fit both simulated and experimental data. We also show that the estimates of some critical kinetic parameters, such as the average inter-division time, obtained by fitting models to data may depend on the assumed distribution of inter-division times, highlighting the challenges in quantitative understanding of cell kinetics.

We used ultra-deep sequencing to obtain tens of thousands of HIV-1 sequences from regions targeted by CD8+ T lymphocytes from longitudinal samples from three acutely infected subjects, and modeled viral evolution during the critical first weeks of infection. Previous studies suggested that a single virus established productive infection, but these conclusions were tempered because of limited sampling; now, we have greatly increased our confidence in this observation through modeling the observed earliest sample diversity based on vastly more extensive sampling. Conventional sequencing of HIV-1 from acute/early infection has shown different patterns of escape at different epitopes; we investigated the earliest escapes in exquisite detail. Over 3-6 weeks, ultradeep sequencing revealed that the virus explored an extraordinary array of potential escape routes in the process of evading the earliest CD8 T-lymphocyte responses–using 454 sequencing, we identified over 50 variant forms of each targeted epitope during early immune escape, while only 2-7 variants were detected in the same samples via conventional sequencing. In contrast to the diversity seen within epitopes, non-epitope regions, including the Envelope V3 region, which was sequenced as a control in each subject, displayed very low levels of variation. In early infection, in the regions sequenced, the consensus forms did not have a fitness advantage large enough to trigger reversion to consensus amino acids in the absence of immune pressure. In one subject, a genetic bottleneck was observed, with extensive diversity at the second time point narrowing to two dominant escape forms by the third time point, all within two months of infection. Traces of immune escape were observed in the earliest samples, suggesting that immune pressure is present and effective earlier than previously reported; quantifying the loss rate of the founder virus suggests a direct role for CD8 T-lymphocyte responses in viral containment after peak viremia. Dramatic shifts in the frequencies of epitope variants during the first weeks of infection revealed a complex interplay between viral fitness and immune escape.

Why some viruses establish chronic infections while others do not is poorly understood. One possibility is that the host's immune response is impaired during chronic infections and is unable to clear the virus from the host. In this report, we use a recently proposed framework to estimate the per capita killing efficacy of CD8(+) T cells, specific for the polyoma virus (PyV), which establishes a chronic infection in mice. Surprisingly, the estimated per cell killing efficacy of PyV-specific effector CD8(+) T cells during the acute phase of the infection was very similar to the efficacy of effector CD8(+) T cells specific to lymphocytic choriomeningitis virus (LCMV-Armstrong), which is cleared from the host. Our results suggest that persistence of PyV does not result from the generation of an inefficient PyV-specific CD8(+) T cell response, and that other host or viral factors are responsible for the ability of PyV to establish chronic infection.

BACKGROUND: Infection with Mycobacterium tuberculosis (Mtb) results in different clinical outcomes ranging from asymptomatic containment to rapidly progressing tuberculosis (TB). The mechanisms controlling TB progression in immunologically-competent hosts remain unclear.

METHODOLOGY/PRINCIPAL FINDINGS: To address these mechanisms, we analyzed TB progression in a panel of genetically heterogeneous (A/SnxI/St) F2 mice, originating from TB-highly-susceptible I/St and more resistant A/Sn mice. In F2 mice the rates of TB progression differed. In mice that did not reach terminal stage of infection, TB progression did not correlate with lung Mtb loads. Nor was TB progression correlated with lung expression of factors involved in antibacterial immunity, such as iNOS, IFN-gamma, or IL-12p40. The major characteristics of progressing TB was high lung expression of the inflammation-related factors IL-1beta, IL-6, IL-11 (p<0.0003); CCL3, CCL4, CXCL2 (p<0.002); MMP-8 (p<0.0001). The major predictors of TB progression were high expressions of IL-1beta and IL-11. TNF-alpha had both protective and harmful effects. Factors associated with TB progression were expressed mainly by macrophages (F4-80(+) cells) and granulocytes (Gr-1(hi)/Ly-6G(hi) cells). Macrophages and granulocytes from I/St and A/Sn parental strains exhibited intrinsic differences in the expression of inflammatory factors, suggesting that genetically determined peculiarities of phagocytes transcriptional response could account for the peculiarities of gene expression in the infected lungs. Another characteristic feature of progressing TB was the accumulation in the infected lungs of Gr-1(dim) cells that could contribute to TB progression.

CONCLUSIONS/SIGNIFICANCE: In a population of immunocompetent hosts, the outcome of TB depends on quantitatively- and genetically-controlled differences in the intensity of inflammatory responses, rather than being a direct consequence of mycobacterial colonization. Local accumulation of Gr-1(dim) cells is a newly identified feature of progressing TB. High expression of IL-1beta and IL-11 are potential risk factors for TB progression and possible targets for TB immunomodulation.

Estimation of division and death rates of lymphocytes in different conditions is vital for quantitative understanding of the immune system. Deuterium, in the form of deuterated glucose or heavy water, can be used to measure rates of proliferation and death of lymphocytes in vivo. Inferring these rates from labeling and delabeling curves has been subject to considerable debate with different groups suggesting different mathematical models for that purpose. We show that the three most common models, which are based on quite different biological assumptions, actually predict mathematically identical labeling curves with one parameter for the exponential up and down slope, and one parameter defining the maximum labeling level. By extending these previous models, we here propose a novel approach for the analysis of data from deuterium labeling experiments. We construct a model of "kinetic heterogeneity" in which the total cell population consists of many sub-populations with different rates of cell turnover. In this model, for a given distribution of the rates of turnover, the predicted fraction of labeled DNA accumulated and lost can be calculated. Our model reproduces several previously made experimental observations, such as a negative correlation between the length of the labeling period and the rate at which labeled DNA is lost after label cessation. We demonstrate the reliability of the new explicit kinetic heterogeneity model by applying it to artificially generated datasets, and illustrate its usefulness by fitting experimental data. In contrast to previous models, the explicit kinetic heterogeneity model 1) provides a novel way of interpreting labeling data; 2) allows for a non-exponential loss of labeled cells during delabeling, and 3) can be used to describe data with variable labeling length.

Quantitative understanding of immunology requires the development of experimental and mathematical techniques for estimation of rates of division and death of lymphocytes under different conditions. Here, we review the advantages and limitations of several labelling methods that are currently used to quantify turnover of lymphocytes in vivo. In addition to highlighting insights into lymphocyte kinetics which have recently been gained thanks to the development of novel techniques, we discuss important directions for future experimental and theoretical work in the field of lymphocyte turnover.

Identification of the transmitted/founder virus makes possible, for the first time, a genome-wide analysis of host immune responses against the infecting HIV-1 proteome. A complete dissection was made of the primary HIV-1-specific T cell response induced in three acutely infected patients. Cellular assays, together with new algorithms which identify sites of positive selection in the virus genome, showed that primary HIV-1-specific T cells rapidly select escape mutations concurrent with falling virus load in acute infection. Kinetic analysis and mathematical modeling of virus immune escape showed that the contribution of CD8 T cell-mediated killing of productively infected cells was earlier and much greater than previously recognized and that it contributed to the initial decline of plasma virus in acute infection. After virus escape, these first T cell responses often rapidly waned, leaving or being succeeded by T cell responses to epitopes which escaped more slowly or were invariant. These latter responses are likely to be important in maintaining the already established virus set point. In addition to mutations selected by T cells, there were other selected regions that accrued mutations more gradually but were not associated with a T cell response. These included clusters of mutations in envelope that were targeted by NAbs, a few isolated sites that reverted to the consensus sequence, and bystander mutations in linkage with T cell-driven escape.

Despite recent advances in immunology, several key parameters determining virus dynamics in infected hosts remain largely unknown. For example, the rate at which specific effector and memory CD8 T cells clear virus-infected cells in vivo is hardly known for any viral infection. We propose a framework to quantify T-cell-mediated killing of infected or peptide-pulsed target cells using the widely used in vivo cytotoxicity assay. We have reanalyzed recently published data on killing of peptide-pulsed splenocytes by cytotoxic T lymphocytes and memory CD8 T cells specific to NP396 and GP276 epitopes of lymphocytic choriomeningitis virus (LCMV) in the mouse spleen. Because there are so many effector CD8 T cells in spleens of mice at the peak of the immune response, NP396- and GP276-pulsed targets are estimated to have very short half-lives of 2 and 14 min, respectively. After the effector numbers have diminished, i.e., in LCMV-immune mice, the half-lives become 48 min and 2.8 h for NP396- and GP276-expressing targets, respectively. Analysis of several alternative models demonstrates that the estimates of half-life times of peptide-pulsed targets are not affected when changes are made in the model assumptions. Our report provides a unifying framework to compare killing efficacies of CD8 T-cell responses specific to different viral and bacterial infections in vivo, which may be used to compare efficacies of various cytotoxic-T-lymphocyte-based vaccines.

It is widely believed that the gut, and particularly the lamina propria (LP) of the gut, contains most of the lymphocytes in humans. The strong depletion of CD4(+) T cells from the gut LP of HIV-infected patients was, therefore, suggested to be such a large, irreversible insult that it could explain HIV disease progression. However, reviewing data from different mammalian species, we found that only 5%-20% of all lymphocytes reside in the gut, and that only 1%-9% of the total lymphocyte number is located in the gut LP. Our findings suggest that spleen and lymph nodes, rather than the gut, are the largest immune compartments in mammals.

Despite the rapid accumulation of quantitative data on the dynamics of CD8(+) T cell responses following acute viral or bacterial infections of mice, the pathways of differentiation of naive CD8(+) T cells into memory during an immune response remain controversial. Currently, three models have been proposed. In the "stem cell-associated differentiation" model, following activation, naive T cells differentiate into stem cell-like memory cells, which then convert into terminally differentiated short-lived effector cells. In the "linear differentiation" model, following activation, naive T cells first differentiate into effectors, and after Ag clearance, effectors convert into memory cells. Finally, in the "progressive differentiation" model, naive T cells differentiate into memory or effector cells depending on the amount of specific stimulation received, with weaker stimulation resulting in formation of memory cells. This study investigates whether the mathematical models formulated from these hypotheses are consistent with the data on the dynamics of the CD8(+) T cell response to lymphocytic choriomeningitis virus during acute infection of mice. Findings indicate that two models, the stem cell-associated differentiation model and the progressive differentiation model, in which differentiation of cells is strongly linked to the number of cell divisions, fail to describe the data at biologically reasonable parameter values. This work suggests additional experimental tests that may allow for further discrimination between different models of CD8(+) T cell differentiation in acute infections.