Is the exquisite specificity of lymphocytes generated by thymic selection or due to evolution?

We have previously argued that the antigen receptors of T and B lymphocytes evolved to be sufficiently specific to avoid massive deletion of clonotypes by negative selection. Their optimal 'specificity' level, i.e., probability of binding any particular epitope, was shown to be inversely related to the number of self-antigens that the cells have to be tolerant to. Experiments have demonstrated that T lymphocytes also become more specific during negative selection in the thymus, because cells expressing the most crossreactive receptors have the highest likelihood of binding a self-antigen, and hence to be tolerized (i.e., deleted, anergized, or diverted into a regulatory T cell phenotype). Thus, there are two -not mutually exclusive- explanations for the exquisite specificity of T cells, one involving evolution and the other thymic selection. To better understand the impact of both, we extend a previously developed mathematical model by allowing for T cells with very different binding probabilities in the pre-selection repertoire. We confirm that negative selection tends to tolerize the most crossreactive clonotypes. As a result, the average level of specificity in the functional post-selection repertoire depends on the number of self-antigens, even if there is no evolutionary optimization of binding probabilities. However, the evolutionary optimal range of binding probabilities in the pre-selection repertoire also depends on the number of self-antigens. Species with more self antigens need more specific pre-selection repertoires to avoid excessive loss of T cells during thymic selection, and hence mount protective immune responses. We conclude that both evolution and negative selection are responsible for the high level of specificity of lymphocytes.

Better safe than sorry: Naive T-cell dynamics in healthy ageing.

It is well-known that the functioning of the immune system gradually deteriorates with age, and we are increasingly confronted with its consequences as the life expectancy of the human population increases. Changes in the T-cell pool are among the most prominent features of the changing immune system during healthy ageing, and changes in the naive T-cell pool in particular are generally held responsible for its gradual deterioration. These changes in the naive T-cell pool are thought to be due to involution of the thymus. It is commonly believed that the gradual loss of thymic output induces compensatory mechanisms to maintain the number of naive T cells at a relatively constant level, and induces a loss of diversity in the T-cell repertoire. Here we review the studies that support or challenge this widely-held view of immune ageing and discuss the implications for vaccination strategies.

On the feasibility of using TCR sequencing to follow a vaccination response - lessons learned.

T cells recognize pathogens by their highly specific T-cell receptor (TCR), which can bind small fragments of an antigen presented on the Major Histocompatibility Complex (MHC). Antigens that are provided through vaccination cause specific T cells to respond by expanding and forming specific memory to combat a future infection. Quantification of this T-cell response could improve vaccine monitoring or identify individuals with a reduced ability to respond to a vaccination. In this proof-of-concept study we use longitudinal sequencing of the TCRß repertoire to quantify the response in the CD4+ memory T-cell pool upon pneumococcal conjugate vaccination. This comes with several challenges owing to the enormous size and diversity of the T-cell pool, the limited frequency of vaccine-specific TCRs in the total repertoire, and the variation in sample size and quality. We defined quantitative requirements to classify T-cell expansions and identified critical parameters that aid in reliable analysis of the data. In the context of pneumococcal conjugate vaccination, we were able to detect robust T-cell expansions in a minority of the donors, which suggests that the T-cell response against the conjugate in the pneumococcal vaccine is small and/or very broad. These results indicate that there is still a long way to go before TCR sequencing can be reliably used as a personal biomarker for vaccine-induced protection. Nevertheless, this study highlights the importance of having multiple samples containing sufficient T-cell numbers, which will support future studies that characterize T-cell responses using longitudinal TCR sequencing.

Challenges in cybersecurity: Lessons from biological defense systems.

Defending against novel, repeated, or unpredictable attacks, while avoiding attacks on the 'self', are the central problems of both mammalian immune systems and computer systems. Both systems have been studied in great detail, but with little exchange of information across the different disciplines. Here, we present a conceptual framework for structured comparisons across the fields of biological immunity and cybersecurity, by framing the context of defense, considering different (combinations of) defensive strategies, and evaluating defensive performance. Throughout this paper, we pose open questions for further exploration. We hope to spark the interdisciplinary discovery of general principles of optimal defense, which can be understood and applied in biological immunity, cybersecurity, and other defensive realms.

What explains the poor contraction of the viral load during paediatric HIV infection?

An acute HIV infection in young children differs markedly from that in adults: Children have higher viral loads (VL), and a poor contraction to a setpoint VL that is not much lower than the peak VL. As a result, children progress faster towards AIDS in the absence of treatment. We used a classical ordinary differential equation model for viral infection dynamics to study why children have a lower viral contraction ratio than adults. We performed parameter sweeps to identify factors explaining the observed difference between children and adults. We grouped parameters associated with the host, the infection, or the immune response. Based on paediatric data available from datasets within the EPIICAL project (https://www.epiical.org/), we refuted that viral replication rates differ between young children and adults, and therefore these cannot be responsible for the low VL contraction ratios seen in children. The major differences in lowering VL contraction ratio resulted from sweeping the parameters linked to the immune response. Thus, we postulate that an "ineffective" (late and/or weak) immune response is the most parsimonious explanation for the higher setpoint VL in young children, and hence the reason for their fast disease progression.

Modeling T Cell Fate.

Many of the pathways that underlie the diversification of naive T cells into effector and memory subsets, and the maintenance of these populations, remain controversial. In recent years a variety of experimental tools have been developed that allow us to follow the fates of cells and their descendants. In this review we describe how mathematical models provide a natural language for describing the growth, loss, and differentiation of cell populations. By encoding mechanistic descriptions of cell behavior, models can help us interpret these new datasets and reveal the rules underpinning T cell fate decisions, both at steady state and during immune responses.

Age-Dependent Normalization Functions for T Lymphocytes in Healthy Individuals.

Lymphocyte numbers naturally change through age. Normalization functions to account for this are sparse and mostly disregard measurements from children in which these changes are most prominent. In this study, we analyze cross-sectional numbers of mainly T lymphocytes (CD3+, CD3+CD4+, and CD3+CD8+) and their subpopulations (naive and memory) from 673 healthy Dutch individuals ranging from infancy to adulthood (0-62 y). We fitted the data by a delayed exponential function and estimated parameters for each lymphocyte subset. Our modeling approach follows general laboratory measurement procedures in which absolute cell counts of T lymphocyte subsets are calculated from observed percentages within a reference population that is truly counted (typically the total lymphocyte count). Consequently, we obtain one set of parameter estimates per T cell subset representing both the trajectories of their counts and percentages. We allow for an initial time delay of half a year before the total lymphocyte counts per microliter of blood start to change exponentially, and we find that T lymphocyte trajectories tend to increase during the first half a year of life. Thus, our study provides functions describing the general trajectories of T lymphocyte counts and percentages of the Dutch population. These functions provide important references to study T lymphocyte dynamics in disease, and they allow one to quantify losses and gains in longitudinal data, such as the CD4+ T cell decline in HIV-infected children and/or the rate of T cell recovery after the onset of treatment.

Compartmentalization and persistence of dominant (regulatory) T cell clones indicates antigen skewing in juvenile idiopathic arthritis.

Autoimmune inflammation is characterized by tissue infiltration and expansion of antigen-specific T cells. Although this inflammation is often limited to specific target tissues, it remains yet to be explored whether distinct affected sites are infiltrated with the same, persistent T cell clones. Here, we performed CyTOF analysis and T cell receptor (TCR) sequencing to study immune cell composition and (hyper-)expansion of circulating and joint-derived Tregs and non-Tregs in juvenile idiopathic arthritis (JIA). We studied different joints affected at the same time, as well as over the course of relapsing-remitting disease. We found that the composition and functional characteristics of immune infiltrates are strikingly similar between joints within one patient, and observed a strong overlap between dominant T cell clones, especially Treg, of which some could also be detected in circulation and persisted over the course of relapsing-remitting disease. Moreover, these T cell clones were characterized by a high degree of sequence similarity, indicating the presence of TCR clusters responding to the same antigens. These data suggest that in localized autoimmune disease, there is autoantigen-driven expansion of both Teffector and Treg clones that are highly persistent and are (re)circulating. These dominant clones might represent interesting therapeutic targets.

Effect of cellular aging on memory T-cell homeostasis.

The fact that T-cell numbers remain relatively stable throughout life, and that T-cell proliferation rates increase during lymphopenia, has led to the consensus that T-cell numbers are regulated in a density-dependent manner. Competition for resources among memory T cells has been proposed to underlie this 'homeostatic' regulation. We first review how two classic models of resource competition affect the T-cell receptor (TCR) diversity of the memory T-cell pool. First, 'global' competition for cytokines leads to a skewed repertoire that tends to be dominated by the very first immune response. Second, additional 'cognate' competition for specific antigens results in a very diverse and stable memory T-cell pool, allowing every antigen to be remembered, which we therefore define as the 'gold-standard'. Because there is limited evidence that memory T cells of the same specificity compete more strongly with each other than with memory T cells of different specificities, i.e., for 'cognate' competition, we investigate whether cellular aging could account for a similar level of TCR diversity. We define cellular aging as a declining cellular fitness due to reduced proliferation. We find that the gradual erosion of previous T-cell memories due to cellular aging allows for better establishment of novel memories and for a much higher level of TCR diversity compared to global competition. A small continual source (either from stem-cell-like memory T-cells or from naive T-cells due to repeated antigen exposure) improves the diversity of the memory T-cell pool, but remarkably, only in the cellular aging model. We further show that the presence of a source keeps the inflation of chronic memory responses in check by maintaining the immune memories to non-chronic antigens. We conclude that cellular aging along with a small source provides a novel and immunologically realistic mechanism to achieve and maintain the 'gold-standard' level of TCR diversity in the memory T-cell pool.

Replicative history marks transcriptional and functional disparity in the CD8+ T cell memory pool.

Clonal expansion is a core aspect of T cell immunity. However, little is known with respect to the relationship between replicative history and the formation of distinct CD8+ memory T cell subgroups. To address this issue, we developed a genetic-tracing approach, termed the DivisionRecorder, that reports the extent of past proliferation of cell pools in vivo. Using this system to genetically 'record' the replicative history of different CD8+ T cell populations throughout a pathogen-specific immune response, we demonstrate that the central memory T (TCM) cell pool is marked by a higher number of prior divisions than the effector memory T cell pool, owing to the combination of strong proliferative activity during the acute immune response and selective proliferative activity after pathogen clearance. Furthermore, by combining DivisionRecorder analysis with single-cell transcriptomics and functional experiments, we show that replicative history identifies distinct cell pools within the TCM compartment. Specifically, we demonstrate that lowly divided TCM cells display enriched expression of stem-cell-associated genes, exist in a relatively quiescent state, and are superior in eliciting a proliferative recall response upon activation. These data provide the first evidence that a stem-cell-like memory T cell pool that reconstitutes the CD8+ T cell effector pool upon reinfection is marked by prior quiescence.

Quantification of CD4 Recovery in Early-Treated Infants Living With HIV.

BACKGROUND: Perinatally HIV-acquired infants benefit from an early antiretroviral treatment initiation. Thanks to a short viral exposure time, their immune system can be maintained or reconstituted, allowing a "normal" immune development. METHODS: In this study, we mathematically modeled and quantified individual CD4+ T-cell reconstitution of a subset of 276 children who started treatment within 6 months of age and achieved sustained viral suppression. Considering natural age differences in CD4+ T-cell dynamics, we fitted distances to age-matched healthy reference values with a linear model approaching an asymptote. RESULTS: Depleted CD4+ percentages (CD4%) and CD4+ counts (CD4ct) restored healthy levels during treatment. CD4ct recovered with a median rate of 4 cells/µL/d, and individual recovery rates were correlated negatively with their initial CD4ct. CD4 values at onset of treatment decrease with age, whereas recovery times and levels seem to be age-independent. CD4 recovery correlates positively with viral suppression, and the stabilization of CD4 levels usually occurs after viral suppression. CD4 levels stabilize within 3-13 months after treatment initiation. The recovery dynamics of the CD4% is comparable with those of the CD4ct. CONCLUSIONS: In early-treated children with successful viral suppression, the CD4 depletion is typically mild and CD4+ T cells tend to "fully" recover in numbers.

Turnover of Murine Cytomegalovirus-Expanded CD8+ T Cells Is Similar to That of Memory Phenotype T Cells and Independent of the Magnitude of the Response.

The potential of memory T cells to provide protection against reinfection is beyond question. Yet, it remains debated whether long-term T cell memory is due to long-lived memory cells. There is ample evidence that blood-derived memory phenotype CD8+ T cells maintain themselves through cell division, rather than through longevity of individual cells. It has recently been proposed, however, that there may be heterogeneity in the lifespans of memory T cells, depending on factors such as exposure to cognate Ag. CMV infection induces not only conventional, contracting T cell responses, but also inflationary CD8+ T cell responses, which are maintained at unusually high numbers, and are even thought to continue to expand over time. It has been proposed that such inflating T cell responses result from the accumulation of relatively long-lived CMV-specific memory CD8+ T cells. Using in vivo deuterium labeling and mathematical modeling, we found that the average production rates and expected lifespans of mouse CMV-specific CD8+ T cells are very similar to those of bulk memory-phenotype CD8+ T cells. Even CMV-specific inflationary CD8+ T cell responses that differ 3-fold in size were found to turn over at similar rates.

Quantification of T-cell dynamics during latent cytomegalovirus infection in humans.

Cytomegalovirus (CMV) infection has a major impact on the T-cell pool, which is thought to be associated with ageing of the immune system. The effect on the T-cell pool has been interpreted as an effect of CMV on non-CMV specific T-cells. However, it remains unclear whether the effect of CMV could simply be explained by the presence of large, immunodominant, CMV-specific memory CD8+ T-cell populations. These have been suggested to establish through gradual accumulation of long-lived cells. However, little is known about their maintenance. We investigated the effect of CMV infection on T-cell dynamics in healthy older adults, and aimed to unravel the mechanisms of maintenance of large numbers of CMV-specific CD8+ T-cells. We studied the expression of senescence, proliferation, and apoptosis markers and quantified the in vivo dynamics of CMV-specific and other memory T-cell populations using in vivo deuterium labelling. Increased expression of late-stage differentiation markers by CD8+ T-cells of CMV+ versus CMV- individuals was not solely explained by the presence of large, immunodominant CMV-specific CD8+ T-cell populations. The lifespans of circulating CMV-specific CD8+ T-cells did not differ significantly from those of bulk memory CD8+ T-cells, and the lifespans of bulk memory CD8+ T-cells did not differ significantly between CMV- and CMV+ individuals. Memory CD4+ T-cells of CMV+ individuals showed increased expression of late-stage differentiation markers and decreased Ki-67 expression. Overall, the expression of senescence markers on T-cell populations correlated positively with their expected in vivo lifespan. Together, this work suggests that i) large, immunodominant CMV-specific CD8+ T-cell populations do not explain the phenotypical differences between CMV+ and CMV- individuals, ii) CMV infection hardly affects the dynamics of the T-cell pool, and iii) large numbers of CMV-specific CD8+ T-cells are not due to longer lifespans of these cells.

Hematopoiesis in numbers.

Hematopoiesis is a dynamic process in which stem and progenitor cells give rise to the ~1013 blood and immune cells distributed throughout the human body. We argue that a quantitative description of hematopoiesis can help consolidate existing data, identify knowledge gaps, and generate new hypotheses. Here, we review known numbers in murine and, where possible, human hematopoiesis, and consolidate murine numbers into a set of reference values. We present estimates of cell numbers, division and differentiation rates, cell size, and macromolecular composition for each hematopoietic cell type. We also propose guidelines to improve the reporting of measurements and highlight areas in which quantitative data are lacking. Overall, we show how quantitative approaches can be used to understand key properties of hematopoiesis.

The limitations, dangers, and benefits of simple methods for testing identifiability.

In their Commentary paper, Villaverde and Massonis (On testing structural identifiability by a simple scaling method: relying on scaling symmetries can be misleading) have commented on our paper in which we proposed a simple scaling method to test structural identifiability. Our scaling invariance method (SIM) tests for scaling symmetries only, and Villaverde and Massonis correctly show the SIM may fail to detect identifiability problems when a model has other types of symmetries. We agree with the limitations raised by these authors but, also, we emphasize that the method is still valuable for its applicability to a wide variety of models, its simplicity, and even as a tool to introduce the problem of identifiability to investigators with little training in mathematics.

TCRß rearrangements without a D segment are common, abundant, and public.

T cells play an important role in adaptive immunity. An enormous clonal diversity of T cells with a different specificity, encoded by the T cell receptor (TCR), protect the body against infection. Most TCRß chains are generated from a V, D, and J segment during recombination in the thymus. Although complete absence of the D segment is not easily detectable from sequencing data, we find convincing evidence for a substantial proportion of TCRß rearrangements lacking a D segment. Additionally, sequences without a D segment are more likely to be abundant within individuals and/or shared between individuals. Our analysis indicates that such sequences are preferentially generated during fetal development and persist within the elderly. Summarizing, TCRß rearrangements without a D segment are not uncommon, and tend to allow for TCRß chains with a high abundance in the naive repertoire.

Cytotoxic T cells are able to efficiently eliminate cancer cells by additive cytotoxicity.

Lethal hit delivery by cytotoxic T lymphocytes (CTL) towards B lymphoma cells occurs as a binary, "yes/no" process. In non-hematologic solid tumors, however, CTL often fail to kill target cells during 1:1 conjugation. Here we describe a mechanism of "additive cytotoxicity" by which time-dependent integration of sublethal damage events, delivered by multiple CTL transiting between individual tumor cells, mediates effective elimination. Reversible sublethal damage includes perforin-dependent membrane pore formation, nuclear envelope rupture and DNA damage. Statistical modeling reveals that 3 serial hits delivered with decay intervals below 50 min discriminate between tumor cell death or survival after recovery. In live melanoma lesions in vivo, sublethal multi-hit delivery is most effective in interstitial tissue where high CTL densities and swarming support frequent serial CTL-tumor cell encounters. This identifies CTL-mediated cytotoxicity by multi-hit delivery as an incremental and tunable process, whereby accelerating damage magnitude and frequency may improve immune efficacy.

Local actin dynamics couple speed and persistence in a cellular Potts model of cell migration.

Cell migration is astoundingly diverse. Molecular signatures, cell-cell interactions, and environmental structures each play their part in shaping cell motion, yielding numerous morphologies and migration modes. Nevertheless, in recent years, a simple unifying law was found to describe cell migration across many different cell types and contexts: faster cells turn less frequently. This universal coupling between speed and persistence (UCSP) was explained by retrograde actin flow from front to back, but it remains unclear how this mechanism generalizes to cells with complex shapes and cells migrating in structured environments, which may not have a well-defined front-to-back orientation. Here, we present an in-depth characterization of an existing cellular Potts model, in which cells polarize dynamically from a combination of local actin dynamics (stimulating protrusions) and global membrane tension along the perimeter (inhibiting protrusions). We first show that the UCSP emerges spontaneously in this model through a cross talk of intracellular mechanisms, cell shape, and environmental constraints, resembling the dynamic nature of cell migration in vivo. Importantly, we find that local protrusion dynamics suffice to reproduce the UCSP-even in cases in which no clear global, front-to-back polarity exists. We then harness the spatial nature of the cellular Potts model to show how cell shape dynamics limit both the speed and persistence a cell can reach and how a rigid environment such as the skin can restrict cell motility even further. Our results broaden the range of potential mechanisms underlying the speed-persistence coupling that has emerged as a fundamental property of migrating cells.

Cell-density independent increased lymphocyte production and loss rates post-autologous HSCT.

Lymphocyte numbers need to be quite tightly regulated. It is generally assumed that lymphocyte production and lifespan increase homeostatically when lymphocyte numbers are low and, vice versa, return to normal once cell numbers have normalized. This widely accepted concept is largely based on experiments in mice, but is hardly investigated in vivo in humans. Here we quantified lymphocyte production and loss rates in vivo in patients 0.5-1 year after their autologous hematopoietic stem cell transplantation (autoHSCT). We indeed found that the production rates of most T- and B-cell subsets in autoHSCT-patients were two to eight times higher than in healthy controls, but went hand in hand with a threefold to ninefold increase in cell loss rates. Both rates also did not normalize when cell numbers did. This shows that increased lymphocyte production and loss rates occur even long after autoHSCT and can persist in the face of apparently normal cell numbers.