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.

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.

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 ~10¹³ 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.

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.

Time to Viral Suppression in Perinatally HIV-Infected Infants Depends on the Viral Load and CD4 T-Cell Percentage at the Start of Treatment

BACKGROUND

Interventions aiming for an HIV cure would benefit from rapid elimination of virus after the onset of antiretroviral therapy (ART), by keeping the latent HIV reservoir small.

SETTING

We investigated HIV suppression in 312 perinatally infected infants starting ART within 6 months after birth from the EPPICC (European Pregnancy and Paediatric HIV Cohort Collaboration).

METHODS

To better understand kinetic differences in HIV suppression among infants, we investigated their individual viral load (VL) decay dynamics. We identified VL decay patterns and determined times to viral suppression (TTS). For infants with strictly declining VLs (n = 188), we used parameter fitting methods to estimate baseline VLs, decay rates, and TTS. We subsequently identified the parameters determining TTS by linear modeling.

RESULTS

The majority of infants suppress HIV VL after the onset of ART. Some children experienced a long TTS due to an "erratic" VL decay pattern. We cannot exclude that this is partly due to treatment complications and subsequent treatment changes, but these children were characterized by significantly lower CD4 percentages (CD4%) at start of treatment compared with those with a "clean" VL decline. Focusing on this "clean" subset, the TTS could be predicted by mathematical modeling, and we identified baseline VL and CD4% as the major factors determining the TTS.

CONCLUSIONS

As VL steeply increases and CD4% constantly decreases in untreated HIV-infected infants, the progression of an HIV infection is largely determined by these 2 factors. To prevent a further disease progression, treatment should be initiated early after contracting HIV, which consequently shortens TTS.

CXCL4 Links Inflammation and Fibrosis by Reprogramming Monocyte-Derived Dendritic Cells in vitro

Fibrosis is a condition shared by numerous inflammatory diseases. Our incomplete understanding of the molecular mechanisms underlying fibrosis has severely hampered effective drug development. CXCL4 is associated with the onset and extent of fibrosis development in multiple inflammatory and fibrotic diseases. Here, we used monocyte-derived cells as a model system to study the effects of CXCL4 exposure on dendritic cell development by integrating 65 longitudinal and paired whole genome transcriptional and methylation profiles. Using data-driven gene regulatory network analyses, we demonstrate that CXCL4 dramatically alters the trajectory of monocyte differentiation, inducing a novel pro-inflammatory and pro-fibrotic phenotype mediated via key transcriptional regulators including CIITA. Importantly, these pro-inflammatory cells directly trigger a fibrotic cascade by producing extracellular matrix molecules and inducing myofibroblast differentiation. Inhibition of CIITA mimicked CXCL4 in inducing a pro-inflammatory and pro-fibrotic phenotype, validating the relevance of the gene regulatory network. Our study unveils that CXCL4 acts as a key secreted factor driving innate immune training and forming the long-sought link between inflammation and fibrosis.

The naive T-cell receptor repertoire has an extremely broad distribution of clone sizes

The clone size distribution of the human naive T-cell receptor (TCR) repertoire is an important determinant of adaptive immunity. We estimated the abundance of TCR sequences in samples of naive T cells from blood using an accurate quantitative sequencing protocol. We observe most TCR sequences only once, consistent with the enormous diversity of the repertoire. However, a substantial number of sequences were observed multiple times. We detect abundant TCR sequences even after exclusion of methodological confounders such as sort contamination, and multiple mRNA sampling from the same cell. By combining experimental data with predictions from models we describe two mechanisms contributing to TCR sequence abundance. TCRα abundant sequences can be primarily attributed to many identical recombination events in different cells, while abundant TCRβ sequences are primarily derived from large clones, which make up a small percentage of the naive repertoire, and could be established early in the development of the T-cell repertoire.

Toxin production spontaneously becomes regulated by local cell density in evolving bacterial populations

The production of anticompetitor toxins is widespread among bacteria. Because production of such toxins is costly, it is typically regulated. In particular, many toxins are produced only when the local cell density is high. It is unclear which selection pressures shaped the evolution of density-dependent regulation of toxin production. Here, we study the evolution of toxin production, resistance and the response to a cell-density cue in a model of an evolving bacterial population with spatial structure. We present results for two growth regimes: (i) an undisturbed, fixed habitat in which only small fluctuations of cell density occur, and (ii) a serial-transfer regime with large fluctuations in cell density. We find that density-dependent toxin production can evolve under both regimes. However, the selection pressures driving the evolution of regulation differ. In the fixed habitat, regulation evolves because it allows cells to produce toxin only when opportunities for reproduction are highly limited (because of a high local cell density), and the effective fitness costs of toxin production are hence low. Under serial transfers, regulation evolves because it allows cells to switch from a fast-growing non-toxic phenotype when colonising a new habitat, to a slower-growing competitive toxic phenotype when the cell density increases. Colonies of such regulating cells rapidly expand into unoccupied space because their edges consist of fast-growing, non-toxin-producing cells, but are also combative because cells at the interfaces with competing colonies do produce toxin. Because under the two growth regimes different types of regulation evolve, our results underscore the importance of growth conditions in the evolution of social behaviour in bacteria.

Bacteria live in microbial communities, in which they compete with many other microbes for nutrients and space. In this competitive environment, almost all known bacterial strains produce toxins that impair or kill other bacteria. This chemical warfare is thought to be one of the major factors shaping microbial diversity. Many toxins are produced only if the local density of bacteria is high. To achieve this, bacteria respond to cell-density cues: signalling molecules or other indicators of the presence of other cells. Here, we use a computational model to study the evolution of density-based regulation of toxin production in bacterial populations. We show that such regulation can arise under various growth conditions, and analyse the selection pressures driving its evolution. In particular, we find that if bacteria regularly need to colonise a new habitat, density-based regulation allows them to express a fast-growing, non-toxic phenotype when expanding into uncolonised territory, and a slower-growing, toxin-producing phenotype when competing with other strains. Colonies of regulating cells show a typical structure, with cells of the fast-growing, sensitive phenotype at their expanding edges, and toxin-producing cells in the colony interior and at interfaces between colonies.

Current best estimates for the average lifespans of mouse and human leukocytes: reviewing two decades of deuterium-labeling experiments

Deuterium is a non-toxic, stable isotope that can safely be administered to humans and mice to study their cellular turnover rates in vivo. It is incorporated into newly synthesized DNA strands during cell division, without interference with the kinetics of cells, and the accumulation and loss of deuterium in the DNA of sorted (sub-)populations of leukocytes can be used to estimate their cellular production rates and lifespans. In the past two decades, this powerful technology has been used to estimate the turnover rates of various types of leukocytes. Although it is the most reliable technique currently available to study leukocyte turnover, there are remarkable differences between the cellular turnover rates estimated by some of these studies. We have recently established that part of this variation is due to (a) difficulties in estimating deuterium availability in some deuterium-labeling studies, and (b) assumptions made by the mathematical models employed to fit the data. Being aware of these two problems, we here aim to approach a consensus on the life expectancies of different types of T cells, B cells, monocytes, and neutrophils in mice and men. We address remaining outstanding problems whenever appropriate and discuss for which immune subpopulations we currently have too little information to draw firm conclusions about their turnover.

Local Attachment Explains Small World-like Properties of Fibroblastic Reticular Cell Networks in Lymph Nodes

Fibroblastic reticular cells (FRCs) form a cellular network that serves as the structural backbone of lymph nodes and facilitates lymphocyte migration. In mice, this FRC network has been found to have small-world properties. Using a model based on geographical preferential attachment, we simulated the formation of a variety of cellular networks and show that similar small-world properties robustly emerge under such natural conditions. By estimating the parameters of this model, we generated FRC network representations with realistic topological properties. We found that the topological properties change markedly when the network is expanded from a thin slice to a three-dimensional cube. Typical small-world properties were found to persist as network size was increased. The simulated networks were very similar to two-dimensional and three-dimensional lattice networks. According to the used metrics, these lattice networks also have small-world properties, indicating that lattice likeness is sufficient to become classified as a small-world network. Our results explain why FRC networks have small-world properties and provide a framework for simulating realistic FRC networks.

Short Lifespans of Memory T-cells in Bone Marrow, Blood, and Lymph Nodes Suggest That T-cell Memory Is Maintained by Continuous Self-Renewal of Recirculating Cells

Memory T-cells are essential to maintain long-term immunological memory. It is widely thought that the bone marrow (BM) plays an important role in the long-term maintenance of memory T-cells. There is controversy however on the longevity and recirculating kinetics of BM memory T-cells. While some have proposed that the BM is a reservoir for long-lived, non-circulating memory T-cells, it has also been suggested to be the preferential site for memory T-cell self-renewal. In this study, we used in vivo deuterium labeling in goats to simultaneously quantify the average turnover rates—and thereby expected lifespans—of memory T-cells from BM, blood and lymph nodes (LN). While the fraction of Ki-67 positive cells, a snapshot marker for recent cell division, was higher in memory T-cells from blood compared to BM and LN, in vivo deuterium labeling revealed no substantial differences in the expected lifespans of memory T-cells between these compartments. Our results support the view that the majority of memory T-cells in the BM are self-renewing as fast as those in the periphery, and are continuously recirculating between the blood, BM, and LN.

Dynamics of Recent Thymic Emigrants in Young Adult Mice

The peripheral naive T-cell pool is generally thought to consist of a subpopulation of recent thymic emigrants (RTEs) and a subpopulation of mature naive (MN) T cells with different dynamics. Thymus transplantation and adoptive transfer studies in mice have provided contradicting results, with some studies suggesting that RTEs are relatively short-lived cells, while another study suggested that RTEs have a survival advantage. We here estimate the death rates of RTE and MN T cells by performing both thymus transplantations and deuterium labeling experiments in mice of at least 12 weeks old, an age at which the size of the T-cell pool has stabilized. For CD4⁺ T cells, we found the total loss rate from the RTE compartment (by death and maturation) to be fourfold faster than that of MN T cells. We estimate the death rate of CD4⁺ RTE to be 0.046 per day, which is threefold faster than the total loss rate from the MN T-cell compartment. For CD8⁺ T cells, we found no evidence for kinetic differences between RTE and MN T cells. Thus, our data support the notion that in young adult mice, CD4⁺ RTE are relatively short-lived cells within the naive CD4⁺ T-cell pool.

An evolutionary perspective on the systems of adaptive immunity

We propose an evolutionary perspective to classify and characterize the diverse systems of adaptive immunity that have been discovered across all major domains of life. We put forward a new function-based classification according to the way information is acquired by the immune systems: Darwinian immunity (currently known from, but not necessarily limited to, vertebrates) relies on the Darwinian process of clonal selection to 'learn' by cumulative trial-and-error feedback; Lamarckian immunity uses templated targeting (guided adaptation) to internalize heritable information on potential threats; finally, shotgun immunity operates through somatic mechanisms of variable targeting without feedback. We argue that the origin of Darwinian (but not Lamarckian or shotgun) immunity represents a radical innovation in the evolution of individuality and complexity, and propose to add it to the list of major evolutionary transitions. While transitions to higher-level units entail the suppression of selection at lower levels, Darwinian immunity re-opens cell-level selection within the multicellular organism, under the control of mechanisms that direct, rather than suppress, cell-level evolution for the benefit of the individual. From a conceptual point of view, the origin of Darwinian immunity can be regarded as the most radical transition in the history of life, in which evolution by natural selection has literally re-invented itself. Furthermore, the combination of clonal selection and somatic receptor diversity enabled a transition from limited to practically unlimited capacity to store information about the antigenic environment. The origin of Darwinian immunity therefore comprises both a transition in individuality and the emergence of a new information system - the two hallmarks of major evolutionary transitions. Finally, we present an evolutionary scenario for the origin of Darwinian immunity in vertebrates. We propose a revival of the concept of the 'Big Bang' of vertebrate immunity, arguing that its origin involved a 'difficult' (i.e. low-probability) evolutionary transition that might have occurred only once, in a common ancestor of all vertebrates. In contrast to the original concept, we argue that the limiting innovation was not the generation of somatic diversity, but the regulatory circuitry needed for the safe operation of amplifiable immune responses with somatically acquired targeting. Regulatory complexity increased abruptly by genomic duplications at the root of the vertebrate lineage, creating a rare opportunity to establish such circuitry. We discuss the selection forces that might have acted at the origin of the transition, and in the subsequent stepwise evolution leading to the modern immune systems of extant vertebrates.

Specificity of inhibitory KIRs enables NK cells to detect changes in an altered peptide environment

The activity of natural killer (NK) cells is tightly regulated by inhibitory and activating receptors. Inhibitory killer immunoglobulin-like receptors (iKIRs) survey the surface of target cells by monitoring the expression of human leukocyte antigen (HLA) class I. The binding of iKIRs has been shown to be sensitive to the peptides presented by HLA class I, implying that iKIRs have the ability to detect the changes in the repertoire of peptide-HLA class I complexes (pHLA), a process occurring during viral infection and in tumor cells. To study how the pHLA repertoire changes upon infection, and whether an iKIR is able to detect these changes, we study peptides eluted from cells prior and after infection with measles virus (MV). Remarkably, most changes in the repertoire of potential iKIR ligands are predicted to be caused by the altered expression of self-peptides. We show that an iKIR can detect these changes in the presented peptides only if it is sufficiently specific, e.g., if iKIRs can distinguish between different amino acids in the contact residues (e.g., position 7 and 8). Our analysis further indicates that one single iKIR per host is not sufficient to detect changes in the peptide repertoire, suggesting that a multigene family encoding for different iKIRs is required for successful peptide recognition.

A Sigmoid Functional Response Emerges When Cytotoxic T Lymphocytes Start Killing Fresh Target Cells

Cytotoxic T lymphocyte (CTL)-mediated killing involves the formation of a synapse with a target cell, followed by delivery of perforin and granzymes. Previously, we derived a general functional response for CTL killing while considering that CTLs form stable synapses (i.e., single-stage) and that the number of conjugates remains at steady state. However, the killing of target cells sometimes requires multiple engagements (i.e., multistage). To study how multistage killing and a lack of steady state influence the functional response, we here analyze a set of differential equations as well as simulations employing the cellular Potts model, in both cases describing CTLs that kill target cells. We find that at steady state the total killing rate (i.e., the number of target cells killed by all CTLs) is well described by the previously derived double saturation function. Compared to single-stage killing, the total killing rate during multistage killing saturates at higher CTL and target cell densities. Importantly, when the killing is measured before the steady state is approached, a qualitatively different functional response emerges for two reasons: First, the killing signal of each CTL gets diluted over several targets and because this dilution effect is strongest at high target cell densities; this can result in a peak in the dependence of the total killing rate on the target cell density. Second, the total killing rate exhibits a sigmoid dependence on the CTL density when killing is a multistage process, because it takes typically more than one CTL to kill a target. In conclusion, a sigmoid dependence of the killing rate on the CTLs during initial phases of killing may be indicative of a multistage killing process. Observation of a sigmoid functional response may thus arise from a dilution effect and is not necessarily due to cooperative behavior of the CTLs.

Tissue Dimensionality Influences the Functional Response of Cytotoxic T Lymphocyte-Mediated Killing of Targets

Cytotoxic T lymphocyte (CTL)-mediated killing of virus infections and tumors occurs over a wide range of conditions. The spatial environments in which CTLs encounter target cells vary from narrow vessels, to two-dimensional epithelial tissues, to densely populated 3-dimensional (3D) T cell areas within lymphoid tissues. How such spatial environments alter the functional response of CTL-mediated killing, i.e., how the killing efficiency depends on cell densities, is unclear. In this study, we perform cellular Potts model simulations in different spatial configurations to investigate how the dimensionality of the space affects the functional response of CTL-mediated killing. Irrespective of the spatial configuration, the function with separate saturation constants for CTL and for target cell densities that we previously proposed can in all cases describe the response, demonstrating its generality. However, the tissue dimensionality determines at which cell densities the killing rate starts to saturate. We show that saturation in a fully 3D environment is stronger than in a "flat" 3D environment, which is largely due to accompanying differences in the CTL-target encounter rates.

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.

RTCR: a pipeline for complete and accurate recovery of T cell repertoires from high throughput sequencing data

Motivation: High Throughput Sequencing (HTS) has enabled researchers to probe the human T cell receptor (TCR) repertoire, which consists of many rare sequences. Distinguishing between true but rare TCR sequences and variants generated by polymerase chain reaction (PCR) and sequencing errors remains a formidable challenge. The conventional approach to handle errors is to remove low quality reads, and/or rare TCR sequences. Such filtering discards a large number of true and often rare TCR sequences. However, accurate identification and quantification of rare TCR sequences is essential for repertoire diversity estimation.

Results: We devised a pipeline, called Recover TCR (RTCR), that accurately recovers TCR sequences, including rare TCR sequences, from HTS data (including barcoded data) even at low coverage. RTCR employs a data-driven statistical model to rectify PCR and sequencing errors in an adaptive manner. Using simulations, we demonstrate that RTCR can easily adapt to the error profiles of different types of sequencers and exhibits consistently high recall and high precision even at low coverages where other pipelines perform poorly. Using published real data, we show that RTCR accurately resolves sequencing errors and outperforms all other pipelines.

Availability and Implementation: The RTCR pipeline is implemented in Python (v2.7) and C and is freely available at http://uubram.github.io/RTCR/along with documentation and examples of typical usage.

Contact:b.gerritsen@uu.nl

The Branching Point in Erythro-Myeloid Differentiation

Development of mature blood cell progenies from hematopoietic stem cells involves the transition through lineage-restricted progenitors. The first branching point along this developmental process is thought to separate the erythro-myeloid and lymphoid lineage fate by yielding two intermediate progenitors, the common myeloid and the common lymphoid progenitors (CMPs and CLPs). Here, we use single-cell lineage tracing to demonstrate that so-called CMPs are highly heterogeneous with respect to cellular output, with most individual CMPs yielding either only erythrocytes or only myeloid cells after transplantation. Furthermore, based on the labeling of earlier progenitors, we show that the divergence between the myeloid and erythroid lineage develops within multipotent progenitors (MPP). These data provide evidence for a model of hematopoietic branching in which multiple distinct lineage commitments occur in parallel within the MPP pool.

Quantifying the effect of Vpu on the promotion of HIV-1 replication in the humanized mouse model

BACKGROUND

Tetherin is an intrinsic anti-viral factor impairing the release of nascent HIV-1 particles from infected cells. Vpu, an HIV-1 accessory protein, antagonizes the anti-viral action of tetherin. Although previous studies using in vitro cell culture systems have revealed the molecular mechanisms of the anti-viral action of tetherin and the antagonizing action of Vpu against tetherin, it still remains unclear how Vpu affects the kinetics of HIV-1 replication in vivo.

RESULTS

To quantitatively assess the role of Vpu in viral replication in vivo, we analyzed time courses of experimental data with viral load and target cell levels in the peripheral blood of humanized mice infected with wild-type and vpu-deficient HIV-1. Our recently developed mathematical model describes the acute phase of this infection reasonably, and allowed us to estimate several parameters characterizing HIV-1 infection in mice. Using a technique of Bayesian parameter estimation, we estimate distributions of the basic reproduction number of wild-type and vpu-deficient HIV-1. This reveals that Vpu markedly increases the rate of viral replication in vivo.

CONCLUSIONS

Combining experiments with mathematical modeling, we provide an estimate for the contribution of Vpu to viral replication in humanized mice.

Reconciling Longitudinal Naive T-Cell and TREC Dynamics during HIV-1 Infection

Naive T cells in untreated HIV-1 infected individuals have a reduced T-cell receptor excision circle (TREC) content. Previous mathematical models have suggested that this is due to increased naive T-cell division. It remains unclear, however, how reduced naive TREC contents can be reconciled with a gradual loss of naive T cells in HIV-1 infection. We performed longitudinal analyses in humans before and after HIV-1 seroconversion, and used a mathematical model to investigate which processes could explain the observed changes in naive T-cell numbers and TRECs during untreated HIV-1 disease progression. Both CD4⁺ and CD8⁺ naive T-cell TREC contents declined biphasically, with a rapid loss during the first year and a much slower loss during the chronic phase of infection. While naive CD8⁺ T-cell numbers hardly changed during follow-up, naive CD4⁺ T-cell counts continually declined. We show that a fine balance between increased T-cell division and loss in the peripheral naive T-cell pool can explain the observed short- and long-term changes in TRECs and naive T-cell numbers, especially if T-cell turnover during the acute phase is more increased than during the chronic phase of infection. Loss of thymic output, on the other hand, does not help to explain the biphasic loss of TRECs in HIV infection. The observed longitudinal changes in TRECs and naive T-cell numbers in HIV-infected individuals are most likely explained by a tight balance between increased T-cell division and death, suggesting that these changes are intrinsically linked in HIV infection.

HIV-1 CCR5 gene therapy will fail unless it is combined with a suicide gene

Highly active antiretroviral therapy (ART) has successfully turned Human immunodeficiency virus type 1 (HIV-1) from a deadly pathogen into a manageable chronic infection. ART is a lifelong therapy which is both expensive and toxic, and HIV can become resistant to it. An alternative to lifelong ART is gene therapy that targets the CCR5 co-receptor and creates a population of genetically modified host cells that are less susceptible to viral infection. With generic mathematical models we show that gene therapy that only targets the CCR5 co-receptor fails to suppress HIV-1 (which is in agreement with current data). We predict that the same gene therapy can be markedly improved if it is combined with a suicide gene that is only expressed upon HIV-1 infection.

Subtle CXCR3-Dependent Chemotaxis of CTLs within Infected Tissue Allows Efficient Target Localization

It is well established how effector T cells exit the vasculature to enter the peripheral tissues in which an infection is ongoing. However, less is known regarding how CTLs migrate toward infected cells after entry into peripheral organs. Recently, it was shown that the chemokine receptor CXCR3 on T cells has an important role in their ability to localize infected cells and to control vaccinia virus infection. However, the search strategy of T cells for virus-infected targets has not been investigated in detail and could involve chemotaxis toward infected cells, chemokinesis (i.e., increased motility) combined with CTL arrest when targets are detected, or both. In this study, we describe and analyze the migration of CTLs within HSV-1-infected epidermis in vivo. We demonstrate that activated T cells display a subtle distance-dependent chemotaxis toward clusters of infected cells and confirm that this is mediated by CXCR3 and its ligands. Although the chemotactic migration is weak, computer simulations based on short-term experimental data, combined with subsequent long-term imaging indicate that this behavior is crucial for efficient target localization and T cell accumulation at effector sites. Thus, chemotactic migration of effector T cells within peripheral tissue forms an important factor in the speed with which T cells are able to arrive at sites of infection.

Crawling and Gliding: A Computational Model for Shape-Driven Cell Migration

Cell migration is a complex process involving many intracellular and extracellular factors, with different cell types adopting sometimes strikingly different morphologies. Modeling realistically behaving cells in tissues is computationally challenging because it implies dealing with multiple levels of complexity. We extend the Cellular Potts Model with an actin-inspired feedback mechanism that allows small stochastic cell rufflings to expand to cell protrusions. This simple phenomenological model produces realistically crawling and deforming amoeboid cells, and gliding half-moon shaped keratocyte-like cells. Both cell types can migrate randomly or follow directional cues. They can squeeze in between other cells in densely populated environments or migrate collectively. The model is computationally light, which allows the study of large, dense and heterogeneous tissues containing cells with realistic shapes and migratory properties.

Cell migration is involved in vital processes like morphogenesis, regeneration and immune system responses, but can also play a central role in pathological processes like metastasization. Computational models have been successfully employed to explain how single cells migrate, and to study how diverse cell-cell interactions contribute to tissue level behavior. However, there are few models that implement realistic cell shapes in multicellular simulations. The method we present here is able to reproduce two different types of motile cells—amoeboid and keratocyte-like cells. Amoeboid cells are highly motile and deform frequently; many cells can act amoeboid in certain circumstances e.g., immune system cells, epithelial cells, individually migrating cancer cells. Keratocytes are (fish) epithelial cells which are famous for their ability to preserve their shape and direction when migrating individually; during wound healing, keratocytes migrate collectively, in sheets, to the site needing reepithelialization. Our method is computationally simple, improves the realism of multicellular simulations and can help assess the tissue level impact of specific cell shapes. For example, it can be employed to study the tissue scanning strategies of leukocytes, the circumstances in which cancer cells adopt amoeboid migration strategies, or the collective migration of keratocytes.

Reconciling Estimates of Cell Proliferation from Stable Isotope Labeling Experiments

Stable isotope labeling is the state of the art technique for in vivo quantification of lymphocyte kinetics in humans. It has been central to a number of seminal studies, particularly in the context of HIV-1 and leukemia. However, there is a significant discrepancy between lymphocyte proliferation rates estimated in different studies. Notably, deuterated ²H₂-glucose (D₂-glucose) labeling studies consistently yield higher estimates of proliferation than deuterated water (D₂O) labeling studies. This hampers our understanding of immune function and undermines our confidence in this important technique. Whether these differences are caused by fundamental biochemical differences between the two compounds and/or by methodological differences in the studies is unknown. D₂-glucose and D₂O labeling experiments have never been performed by the same group under the same experimental conditions; consequently a direct comparison of these two techniques has not been possible. We sought to address this problem. We performed both in vitro and murine in vivo labeling experiments using identical protocols with both D₂-glucose and D₂O. This showed that intrinsic differences between the two compounds do not cause differences in the proliferation rate estimates, but that estimates made using D₂-glucose in vivo were susceptible to difficulties in normalization due to highly variable blood glucose enrichment. Analysis of three published human studies made using D₂-glucose and D₂O confirmed this problem, particularly in the case of short term D₂-glucose labeling. Correcting for these inaccuracies in normalization decreased proliferation rate estimates made using D₂-glucose and slightly increased estimates made using D₂O; thus bringing the estimates from the two methods significantly closer and highlighting the importance of reliable normalization when using this technique.

Stable isotope labeling is used to quantify the rate at which living cells proliferate and die in humans. It has been central to a number of seminal studies, particularly in viral infections such as HIV-1, and leukemia. However, different labels (deuterated water or deuterated glucose) yield different estimates for the rate of cell proliferation and loss; this hampers our understanding and weakens our confidence in this important technique. We performed in vitro and in vivo experiments as well as a new analysis of existing data to directly compare the two labels. This reveals that a major source of the discrepancy lies in the difficulty of assessing deuterated glucose availability. We reconcile published studies and provide recommendations to avoid this problem in the future.

What do mathematical models tell us about killing rates during HIV-1 infection?

Over the past few decades the extent to which cytotoxic T lymphocytes (CTLs) control human immunodeficiency virus (HIV) replication has been studied extensively, yet their role and mode of action remain controversial. In some studies, CTLs were found to kill a large fraction of the productively infected cells relative to the viral cytopathicity, whereas in others CTLs were suggested to kill only a small fraction of infected cells. In this review, we compile published estimates of CTL-mediated death rates, and examine whether these studies permit determining the rate at which CTLs kill HIV-1 infected cells. We highlight potential misinterpretations of the CTL-killing rates from the escape rates of mutants, and from perturbations of the steady state viral load during chronic infection. Our major conclusion is that CTL-mediated killing rates remain unknown. But contrary to current consensus, we argue that killing rates higher than one per day are perfectly consistent with the experimental data, which would imply that the majority of the productively infected cells could still die from CTL-mediated killing rather than from viral cytopathicity.

Can Selective MHC Downregulation Explain the Specificity and Genetic Diversity of NK Cell Receptors?

Natural killer (NK) cells express inhibiting receptors (iNKRs), which specifically bind MHC-I molecules on the surface of healthy cells. When the expression of MHC-I on the cell surface decreases, which might occur during certain viral infections and cancer, iNKRs lose inhibiting signals and the infected cells become target for NK cell activation (missing-self detection). Although the detection of MHC-I deficient cells can be achieved by conserved receptor-ligand interactions, several iNKRs are encoded by gene families with a remarkable genetic diversity, containing many haplotypes varying in gene content and allelic polymorphism. So far, the biological function of this expansion within the NKR cluster has remained poorly understood. Here, we investigate whether the evolution of diverse iNKRs genes can be driven by a specific viral immunoevasive mechanism: selective MHC downregulation. Several viruses, including EBV, CMV, and HIV, decrease the expression of MHC-I to escape from T cell responses. This downregulation does not always affect all MHC loci in the same way, as viruses target particular MHC molecules. To study the selection pressure of selective MHC downregulation on iNKRs, we have developed an agent-based model simulating an evolutionary scenario of hosts infected with herpes-like viruses, which are able to selectively downregulate the expression of MHC-I molecules on the cell surface. We show that iNKRs evolve specificity and, depending on the similarity of MHC alleles within each locus and the differences between the loci, they can specialize to a particular MHC-I locus. The easier it is to classify an MHC allele to its locus, the lower the required diversity of the NKRs. Thus, the diversification of the iNKR cluster depends on the locus specific MHC structure.

A Coevolutionary Arms Race between Hosts and Viruses Drives Polymorphism and Polygenicity of NK Cell Receptors

Natural killer cell receptors (NKRs) monitor the expression of major histocompatibility class I (MHC-I) and stress molecules to detect unhealthy tissue, such as infected or tumor cells. The NKR gene family shows a remarkable genetic diversity, containing several genes encoding receptors with activating and inhibiting signaling, and varying in gene content and allelic polymorphism. The expansion of the NKR genes is species-specific, with different species evolving alternative expanded NKR genes, which encode structurally different proteins, yet perform comparable functions. So far, the biological function of this expansion within the NKR cluster has remained poorly understood. To study the evolution of NKRs, we have developed an agent-based model implementing a coevolutionary scenario between hosts and herpes-like viruses that are able to evade the immune response by downregulating the expression of MHC-I on the cell surface. We show that hosts evolve specific inhibitory NKRs, specialized to particular MHC-I alleles in the population. Viruses in our simulations readily evolve proteins mimicking the MHC molecules of their host, even in the absence of MHC-I downregulation. As a result, the NKR locus becomes polygenic and polymorphic, encoding both specific inhibiting and activating receptors to optimally protect the hosts from coevolving viruses.

Lymphocyte maintenance during healthy aging requires no substantial alterations in cellular turnover

In healthy humans, lymphocyte populations are maintained at a relatively constant size throughout life, reflecting a balance between lymphocyte production and loss. Given the profound immunological changes that occur during healthy aging, including a significant decline in T-cell production by the thymus, lymphocyte maintenance in the elderly is generally thought to require homeostatic alterations in lymphocyte dynamics. Surprisingly, using in vivo²H₂O labeling, we find similar dynamics of most lymphocyte subsets between young adult and elderly healthy individuals. As the contribution of thymic output to T-cell production is only minor from young adulthood onward, compensatory increases in peripheral T-cell division rates are not required to maintain the T-cell pool, despite a tenfold decline in thymic output. These fundamental insights will aid the interpretation of further research into aging and clinical conditions related to disturbed lymphocyte dynamics.

Immuno-epidemiological modeling of HIV-1 predicts high heritability of the set-point virus load, while selection for CTL escape dominates virulence evolution

It has been suggested that HIV-1 has evolved its set-point virus load to be optimized for transmission. Previous epidemiological models and studies into the heritability of set-point virus load confirm that this mode of adaptation within the human population is feasible. However, during the many cycles of replication between infection of a host and transmission to the next host, HIV-1 is under selection for escape from immune responses, and not transmission. Here we investigate with computational and mathematical models how these two levels of selection, within-host and between-host, are intertwined. We find that when the rate of immune escape is comparable to what has been observed in patients, immune selection within hosts is dominant over selection for transmission. Surprisingly, we do find high values for set-point virus load heritability, and argue that high heritability estimates can be caused by the 'footprints' left by differing hosts' immune systems on the virus.

Random migration and signal integration promote rapid and robust T cell recruitment

To fight infections, rare T cells must quickly home to appropriate lymph nodes (LNs), and reliably localize the antigen (Ag) within them. The first challenge calls for rapid trafficking between LNs, whereas the second may require extensive search within each LN. Here we combine simulations and experimental data to investigate which features of random T cell migration within and between LNs allow meeting these two conflicting demands. Our model indicates that integrating signals from multiple random encounters with Ag-presenting cells permits reliable detection of even low-dose Ag, and predicts a kinetic feature of cognate T cell arrest in LNs that we confirm using intravital two-photon data. Furthermore, we obtain the most reliable retention if T cells transit through LNs stochastically, which may explain the long and widely distributed LN dwell times observed in vivo. Finally, we demonstrate that random migration, both between and within LNs, allows recruiting the majority of cognate precursors within a few days for various realistic infection scenarios. Thus, the combination of two-scale stochastic migration and signal integration is an efficient and robust strategy for T cell immune surveillance.

Reliable reconstruction of HIV-1 whole genome haplotypes reveals clonal interference and genetic hitchhiking among immune escape variants

BACKGROUND

Following transmission, HIV-1 evolves into a diverse population, and next generation sequencing enables us to detect variants occurring at low frequencies. Studying viral evolution at the level of whole genomes was hitherto not possible because next generation sequencing delivers relatively short reads.

RESULTS

We here provide a proof of principle that whole HIV-1 genomes can be reliably reconstructed from short reads, and use this to study the selection of immune escape mutations at the level of whole genome haplotypes. Using realistically simulated HIV-1 populations, we demonstrate that reconstruction of complete genome haplotypes is feasible with high fidelity. We do not reconstruct all genetically distinct genomes, but each reconstructed haplotype represents one or more of the quasispecies in the HIV-1 population. We then reconstruct 30 whole genome haplotypes from published short sequence reads sampled longitudinally from a single HIV-1 infected patient. We confirm the reliability of the reconstruction by validating our predicted haplotype genes with single genome amplification sequences, and by comparing haplotype frequencies with observed epitope escape frequencies.

CONCLUSIONS

Phylogenetic analysis shows that the HIV-1 population undergoes selection driven evolution, with successive replacement of the viral population by novel dominant strains. We demonstrate that immune escape mutants evolve in a dependent manner with various mutations hitchhiking along with others. As a consequence of this clonal interference, selection coefficients have to be estimated for complete haplotypes and not for individual immune escapes.

Quantifying the Protection of Activating and Inhibiting NK Cell Receptors during Infection with a CMV-Like Virus

The responsiveness of natural killer (NK) cells is controlled by balancing signals from activating and inhibitory receptors. The most important ligands of inhibitory NK cell receptors are the highly polymorphic major histocompatibility complex (MHC) class I molecules, which allow NK cells to screen the cellular health of target cells. Although these inhibitory receptor–ligand interactions have been well characterized, the ligands for most activating receptors are still unknown. The mouse cytomegalovirus (MCMV) represents a helpful model to study NK cell-driven immune responses. Many studies have demonstrated that CMV infection can be controlled by NK cells via their activating receptors, but the exact contribution of the different signaling potential (i.e., activating vs. inhibiting) remains puzzling. In this study, we have developed a probabilistic model, which predicts the optimal specificity of inhibitory and activating NK cell receptors needed to offer the best protection against a CMV-like virus. We confirm our analytical predictions with an agent-based model of an evolving host population. Our analysis quantifies the degree of protection of each receptor type, revealing that mixed haplotypes (i.e., haplotypes composed of activating and inhibiting receptors) are most protective against CMV-like viruses, and that the protective effect depends on the number of MHC loci per individual.

Virus encoded MHC-like decoys diversify the inhibitory KIR repertoire

Natural killer (NK) cells are circulating lymphocytes that play an important role in the control of viral infections and tumors. Their functions are regulated by several activating and inhibitory receptors. A subset of these receptors in human NK cells are the killer immunoglobulin-like receptors (KIRs), which interact with the highly polymorphic MHC class I molecules. One important function of NK cells is to detect cells that have down-regulated MHC expression (missing-self). Because MHC molecules have non polymorphic regions, their expression could have been monitored with a limited set of monomorphic receptors. Surprisingly, the KIR family has a remarkable genetic diversity, the function of which remains poorly understood. The mouse cytomegalovirus (MCMV) is able to evade NK cell responses by coding "decoy" molecules that mimic MHC class I. This interaction was suggested to have driven the evolution of novel NK cell receptors. Inspired by the MCMV system, we develop an agent-based model of a host population infected with viruses that are able to evolve MHC down-regulation and decoy molecules. Our simulations show that specific recognition of MHC class I molecules by inhibitory KIRs provides excellent protection against viruses evolving decoys, and that the diversity of inhibitory KIRs will subsequently evolve as a result of the required discrimination between host MHC molecules and decoy molecules.

Immune activation and collateral damage in AIDS pathogenesis

In the past decade, evidence has accumulated that human immunodeficiency virus (HIV)-induced chronic immune activation drives progression to AIDS. Studies among different monkey species have shown that the difference between pathological and non-pathological infection is determined by the response of the immune system to the virus, rather than its cytopathicity. Here we review the current understanding of the various mechanisms driving chronic immune activation in HIV infection, the cell types involved, its effects on HIV-specific immunity, and how persistent inflammation may cause AIDS and the wide spectrum of non-AIDS related pathology. We argue that therapeutic relief of inflammation may be beneficial to delay HIV-disease progression and to reduce non-AIDS related pathological side effects of HIV-induced chronic immune stimulation.

The rate of immune escape vanishes when multiple immune responses control an HIV infection

During the first months of HIV infection, the virus typically evolves several immune escape mutations. These mutations are found in epitopes in viral proteins and reduce the impact of the CD8⁺ T cells specific for these epitopes. Recent data show that only a subset of the epitopes escapes, that most of these escapes evolve early, and that the rate of immune escape slows down considerably. To investigate why the evolution of immune escape slows down over the time of infection, we have extended a consensus mathematical model to allow several immune responses to control the virus together. In the extended model, most escapes also occur early, and the immune escape rate becomes small later, and typically only a minority of the epitopes escape. We show that escaping one of the many immune responses provides little advantage after viral setpoint has been approached because the total killing rate hardly depends on the breadth of the immune response. If the breadth of the immune response slowly wanes during disease progression, the model predicts an increase in the rate of immune escape at late stages of infection. Overall, the most striking prediction of the model is that HIV evolves a small number of immune escapes, in both relative and absolute terms, when the CTL immune response is broad.

Polar auxin transport: models and mechanisms

Spatial patterns of the hormone auxin are important drivers of plant development. The observed feedback between the active, directed transport that generates auxin patterns and the auxin distribution that influences transport orientation has rendered this a popular subject for modelling studies. Here we propose a new mathematical framework for the analysis of polar auxin transport and present a detailed mathematical analysis of published models. We show that most models allow for self-organised patterning for similar biological assumptions, and find that the pattern generated is typically unidirectional, unless additional assumptions or mechanisms are incorporated. Our analysis thus suggests that current models cannot explain the bidirectional fountain-type patterns found in plant meristems in a fully self-organised manner, and we discuss future research directions to address the gaps in our understanding of auxin transport mechanisms.