Differential impact of self and environmental antigens on the ontogeny and maintenance of CD4+ T cell memory
Abstract
Laboratory mice develop populations of circulating memory CD4+ T cells in the absence of overt infection. We have previously shown that these populations are replenished from naive precursors at high levels throughout life (Gossel et al., 2017). However, the nature, relative importance and timing of the forces generating these cells remain unclear. Here, we tracked the generation of memory CD4+ T cell subsets in mice housed in facilities differing in their ‘dirtiness’. We found evidence for sequential naive to central memory to effector memory development, and confirmed that both memory subsets are heterogeneous in their rates of turnover. We also inferred that early exposure to self and environmental antigens establishes persistent memory populations at levels determined largely, although not exclusively, by the dirtiness of the environment. After the first few weeks of life, however, these populations are continuously supplemented by new memory cells at rates that are independent of environment.
Introduction
Conventional memory T cells are defined as lymphocytes which respond rapidly to previously encountered epitopes (Gourley et al., 2004; Kaech and Wherry, 2007). In mice, memory T cells exhibit considerable heterogeneity in their function, circulation patterns, response to re-challenge, and capacities for proliferative self-renewal and survival (Farber, 2000; Kaech and Wherry, 2007; Jameson and Masopust, 2009; Gossel et al., 2017). This phenotypic heterogeneity is reflected in differential expression of various cell-surface molecules. In uninfected naive mice, there are at least two distinct populations of recirculating cells distinguished by their expression of the lymph node homing receptor L-selectin (CD62L); CD44hi CD62L− effector memory (TEM) and CD44hi CD62L+ central memory (TCM) cells. During immune responses to active infection, there is an even more complex mix of effector and memory intermediates (Jameson and Masopust, 2018).
While it is clear that memory to infection resides amongst these CD44hi subsets, it is also evident that they are generated in naive mice in the absence of overt infection. The functional significance of these memory-phenotype (MP) CD4+ T cells is not fully understood, but there is evidence they can augment primary immune responses. They can facilitate rapid production of IFN- during the early inflammatory phase of the immune response to Toxoplasma gondii infection and enhance TH1-type CD4+ T cell responses later in infection (Kawabe et al., 2017). There is also evidence that MP cells are capable of making rapid cross-reactive responses during primary infections (Min and Paul, 2005). Given that MP cells represent the majority of the memory compartment in specific pathogen-free (SPF) mice (Kawabe et al., 2017), a better understanding of how these cells are generated and maintained is crucial for better understanding their function and impact upon conventional memory to defined challenges.
The precise nature of the forces driving the generation of MP cells remains unclear. Their development appears to require a TCR-mediated activation event; Cd28−/− mice have greatly reduced numbers of MP cells (Kotani et al., 2006), and mice lacking canonical NF-B signalling, an obligate pathway in T cell activation, are completely devoid of such cells (Webb et al., 2019). Whether the TCR stimuli derive from self or foreign recognition events is unknown. MP cell generation may reflect a stochastic process in which CD4+ naive T cells occasionally encounter homeostatic stimuli that are above an activation threshold (Sprent and Surh, 2011), and indeed MP cells are generated not only in lymphopenia but also constitutively under replete conditions throughout life (Gossel et al., 2017; Kawabe et al., 2017). There is a positive correlation between the affinity of naive CD4+ T cells to self-antigens and the potential for differentiating into MP cells (Kawabe et al., 2017), although self-reactivity is also positively correlated to reactivity to foreign antigens (Mandl et al., 2013). Stronger evidence supporting an autoreactive stimulus comes from the failure of broad-spectrum antibiotics to prevent conversion of naive T cells to MP cells following adoptive transfer in vivo, and the observation that mice raised in germ-free conditions contain similar numbers of MP cells as those in SPF conditions (Kawabe et al., 2017). However, there is also a role for foreign environmental antigen in generating T cell memory compartments, since mice raised in more antigenically diverse environments, but still in the absence of overt infection, exhibit larger peripheral memory CD8+ T cell pools than SPF mice (Beura et al., 2016). When exactly such environmental stimuli impact upon memory compartment development, or how foreign and self reactivity combine to form the memory compartments, is unknown.
The differentiation pathways of these MP cells are also poorly understood. They may derive either directly from naive T cells or through interconversion of other memory phenotypes. Amongst CD8+ lineage cells, there is evidence that some MP subpopulations are generated in the thymus (Lee et al., 2011). In the case of CD4+ lineage cells, it has been suggested that MP cells derive from peripheral naive phenotype cells in a thymus-independent fashion (Kawabe et al., 2017). While both CD4+ TCM and TEM are produced constitutively in adult mice (Gossel et al., 2017), it remains unclear how their differentiation patterns relate to those of classically antigen-stimulated naive T cells.
In this study, we aimed to characterise the development and maintenance of memory CD4+ T cell subsets in adult mice to identify the nature and timing of the signals driving these dynamics in the absence of overt infection. To do this, we quantified the homeostasis and ontogeny of memory CD4+ T cells in identical strains of mice raised in two different animal facilities with distinct caging environments; those housed in individually ventilated cages (IVCs) and fed irradiated water, and those fed untreated tap water and housed in open cages, who might consequently be exposed to a greater variety of environmental antigens. We made use of an established temporal fate-mapping model in both cohorts, which allowed us to track the development of T cells under replete conditions (Hogan et al., 2015; Gossel et al., 2017). We also used data from germ-free (GF) mice to dissect further the relative contributions of self and environmental antigens to the generation and maintenance of MP cells.
Results
Modelling the ontogeny and homeostasis of CD4+ MP T cell subsets
We employed a system described previously (Hogan et al., 2015; Gossel et al., 2017) to examine the flow of cells into memory subsets (Figure 1A). Briefly, CD45.1+ C57Bl6/SJL host mice of various ages were treated with optimised doses of the transplant conditioning drug busulfan to selectively deplete haematopoeitic stem cells (HSC). The HSC compartment was then reconstituted by bone marrow transplantation (BMT) with congenically labelled bone marrow from CD45.2+ C57Bl6/J donors. The progeny of donor HSC were then followed as they developed in the thymus and percolated into the peripheral T cell pools, initially replete with host-derived cells. Total numbers of CD4+ naive T cells and TCM and TEM cells in these busulfan chimeric mice were normal, in comparison to untreated wild-type (WT) controls (Figure 1—figure supplement 1), confirming that the busulfan treatment regime left peripheral compartments intact (Hogan et al., 2015; Gossel et al., 2017). The kinetics with which donor-derived cells infiltrate the peripheral compartments – first naive, and then memory subsets – are rich in information regarding developmental pathways, rates of turnover and differentiation of lymphocyte populations, and any heterogeneity in homeostatic dynamics within them (Hogan et al., 2015; Gossel et al., 2017; Rane et al., 2018).
We studied busulfan chimeric mice from two housing facilities that employed different levels of mouse containment. At the MRC National Institute for Medical Research (NIMR), mice were held in open cages and fed untreated tap water, while mice held at the UCL Comparative Biology Unit (UCL) were maintained in individually ventilated cages (IVCs) and fed irradiated water. Henceforth, we refer to UCL sourced mice as ‘clean’ and NIMR sourced mice as ‘dirty’, in reference to the presumed difference in health status of the mice. We use these terms for clarity, but emphasise that they are relative; mice co-housed with pet-store or feral mice would be expected to be substantially ‘dirtier’ (Beura et al., 2016), and those in turn are cleaner than truly feral mice. In both environments, the same C57Bl6/SJL strain was analysed by the same researcher and cells were enumerated using the same single CASY counter. From age 10 weeks onwards, the numbers of CD4+ naive T cells in mice from clean and dirty environments were broadly similar (Figure 1B, left panel). The total sizes (host+donor) of all circulating memory CD4+ T cell subsets remained relatively stable over the time frame of analysis, but were significantly larger in dirty mice (Figure 1B, right panels) at age 10 weeks. Following BMT, donor-derived memory T cells accumulated in similar numbers in the two environments (Figure 1C). Therefore, these two observations result in a lower proportional replacement of pre-existing memory cells with donor memory cells in dirty mice (Figure 1D).
To quantify the cellular processes underlying these kinetics, we first considered a simple mechanistic explanation shown schematically in Figure 2A. In this ‘homogeneous’ model, each memory population (CD4+ TCM or TEM) is fed at a constant per capita rate from a precursor population (source). We refer to this rate as the force of recruitment, . The total cellular flux into a memory subset per day is then multiplied by the size of the source population, which in principle could be CD4+ naive T cells, or the complementary memory subset. We assume that memory cells are then lost at a constant net per capita rate , which is the balance of loss (turnover) and proliferative self-renewal. In particular, the ‘clonal half-life’ is the average time taken for a population that undergoes any degree of self-renewal to halve in size, and may be much longer than the lifespan of any particular cell within it.
We also considered a ‘two phase’ model of memory dynamics (Figure 2B), which was motivated by three observations. First, newly generated donor CD4+ TCM and TEM in busulfan chimeras express Ki67, a marker of recent cell division, at higher levels than their established host-derived counterparts for some time after BMT (Gossel et al., 2017), although these levels eventually converge (data not shown). These observations suggest that memory CD4+ T cell populations become less proliferative, on average, with time since entry into the compartment. Second, we previously found evidence, using BrdU labelling in WT mice, that both CD4+ TEM and TCM appear to be kinetically heterogeneous, comprising at least two subpopulations turning over and dividing at different rates (Gossel et al., 2017). Third, and consistent with this picture, the increases in donor chimerism that we observed in both TCM and TEM with time post-BMT were suggestive of a biphasic kinetic, with a relatively rapid accumulation of donor cells followed by a slower increase (Figure 1D). As shown in Materials and methods, the rate of accumulation of new memory cells is dictated by both the dynamics of influx and the net loss rate of existing memory, . Therefore, all three observations are consistent with a mechanism in which cells newly recruited into memory comprise a subpopulation that both divides rapidly and has a high net loss rate . These cells then transition to a more quiescent state that divides more slowly and also has a lower net loss rate, (Figure 2B). The transient differences in Ki67 expression in donor and host memory cells could then be explained by an enrichment for ‘new’ memory (donor cells) in the fast phase in the weeks immediately following BMT, and not by any intrinsic differences in the behaviour of donor and host cells. Further, this transient difference in Ki67 expression implies a linear flow from fast to slow, rather than a branched process of establishment of the two populations separately; in the latter case, we would expect no differences in Ki67 expression between host and donor cells at any time, provided Ki67 expression within each subset derives entirely from self-renewal and is not inherited from the precursor population.
While the two-phase model is perhaps a minimal description of these observations, it seems plausible that any transition from active to quiescent memory might be more continuous. We previously found evidence for smooth changes in the rates of division and/or loss of naive T cells with their post-thymic age (Rane et al., 2018). We therefore also considered a model in which the net loss rate of a cohort of cells changes continuously with the time since their common ancestor entered memory, (the ‘age-dependent loss model’, Figure 2C). While the observations above are most consistent with , or a decreasing , when fitting the two-phase and age-dependent loss models we placed no constraints on their parameters and allowed the data to determine their values. When analysing the age-dependent loss model, we also explored a variety of forms for (see Materials and methods).
Finally, we considered an alternative form of heterogeneity in memory, in which subpopulations of CD4+ TCM and TEM generated early in life persist and are not replenished by newer cells (Gossel et al., 2017; Figure 2D). These ‘incumbent' or 'resistant memory’ populations, assumed to be stable in numbers and entirely host-derived (that is, established before 5 weeks of age, the earliest age at BMT in this study), could naturally explain the limited donor chimerism within memory subsets and, if they are less dynamic than memory generated later in life, might also be able to explain host/donor differences in Ki67 expression.
The kinetics of accumulation of CD4+ MP cells are consistent with a naive TCM TEM pathway, and both memory subsets are heterogeneous in their turnover
We compared the abilities of the four mechanisms to describe the replacement kinetics of memory subsets in the dirty and clean environments. The kinetic of donor chimerism in CD4+ TEM clearly lagged that of CD4+ TCM (Figure 1D), ruling out TEM as a direct predictor of TCM accumulation. We therefore considered only naive T cells as the source for TCM, but considered both naive and TCM cells as potential sources of TEM.
For each combination of source, environment (clean/dirty), and subset (TEM/TCM), we fitted each model simultaneously to the timecourses of the total cell numbers and the proportion of donor cells within the subset using a maximum likelihood approach. We then calculated the combined probabilities (joint likelihoods) that the replacement kinetics of a given subset in both clean and dirty environments derive from each combination of source and model, allowing for different parameters in clean and dirty mice. We then compared the support for each combination using the Akaike Information Criterion (Table 1, values in bold). Details of the model formulation, model fitting and inference procedures are given in Materials and methods.
We found clearly stronger support for TCM cells rather than naive T cells as a predictor of TEM production (Table 1). This conclusion contrasts with that of our earlier study (Gossel et al., 2017), which found evidence for a direct naive TEM transition; however, while these inferences may be model-dependent to some extent, the more detailed timecourses we studied here gave us greater power to discriminate between the two pathways.
We found almost no support for the homogeneous or resistant memory models. For TCM the age-dependent loss was strongly favoured statistically, while for TEM the two-phase model had the strongest support (fits shown in Figure 3). However, the two models gave visually very similar descriptions of each dataset (not shown). This similarity is perhaps unsurprising, as both describe a progressive increase in clonal persistence the longer cells or their progeny reside within memory. Therefore, we remain somewhat equivocal regarding the true nature of heterogeneity in each, and present parameter estimates below for both models (Table 2). A robust conclusion, however, is that we find a progressive lengthening of clonal lifetimes in both memory subsets and environments, with newly recruited memory being lost on timescales of days to weeks, and more established memory persisting for several months (Figure 4).
Constitutive generation of circulating CD4+ TEM and TCM cells in adult mice occurs at constant rates, irrespective of environment
The kinetics of accumulation of donor-derived memory cells were visually indistinguishable in the two environments from age 10 weeks onwards (Figure 1C), suggesting similar net rates of recruitment and loss. Consistent with this observation, we found no significant differences between clean and dirty mice in the forces of recruitment () into circulating CD4+ memory T cell subsets between the ages of 11–64 weeks, in either model, and no substantial differences in their rates of loss (Figure 4 and Table 2). Therefore, we infer that antigenic stimuli common to both environments drive the continuous generation of new CD4+ effector and central memory T cells in adult mice, at total rates that are proportional to the sizes of their precursor populations.
Quantifying the long-term dynamics of CD4+ MP T cell subsets
Our analyses give a quantitative picture of recruitment into memory and the cells’ subsequent life-histories, and allow us to identify features of their population dynamics that are common to both environments and model-independent. First, the donor chimerism in TCM reached substantially lower levels than that of their CD4+ naive T cell precursors (Figure 1D), suggesting that the rate of generation of new memory in both environments wanes with age, and/or that more established memory has a competitive advantage over recently recruited cells. We find evidence for both processes here. We show in Materials and methods that if influx declines faster than the average rate of turnover, a population will be unable to reach the same level of chimerism as its precursor – in effect, the flow from the source ‘dries up’ more quickly than the memory cells can be replaced by immigrants. We see signs of this effect; due to thymic involution, CD4+ naive T cell numbers decay exponentially in both clean and dirty adult mice (Figure 1B) with population half lives of 228 days (95% CI 227–231 days) and 143 (142-144) days, respectively. These timescales are comparable to or shorter than the estimated half lives of established TCM memory clones ( in the two-phase model; and for a > 30 weeks in the age-dependent loss model; Figure 4 and Table 2). In addition, both the two-phase and age-dependent models indicate that older memory clones have a fitness advantage over newer ones. Therefore, the limited replacement of host CD4+ TCM by donor cells derives from the decline in naive T cell numbers with age, slow average rates of turnover, and the increased persistence of more established memory cells. This slow rate of accumulation of new TCM in turn acts to limit the chimerism observed in TEM, which are also lost slowly.
One can also quantify the fates of populations after entering memory, although here our insights are more model-dependent. The two-phase model predicts that the establishment of memory is relatively inefficient, with ‘fast’ populations lost over timescales of days and only a small proportion of these surviving to become more persistent ‘slow’ memory (∼2–10% of TEM, and ∼3–5% of TCM; Table 2). Despite this inefficiency, the substantial constitutive influxes maintain the fast and slow populations at comparable sizes, consistent with our previous analysis of BrdU labelling of CD4+ TCM and TEM in WT mice (Gossel et al., 2017). In contrast, the age-dependent loss model makes lower estimates of the force of recruitment into memory (Figure 4, left panels) but predicts more efficient establishment, with newly generated memory having clonal half-lives of 20–40 days and a much greater proportion persisting longer term (Figure 4—figure supplement 1).
Larger memory populations in dirty mice derive from early antigen exposure
Given the similarity of the rates of generation of memory in clean and dirty adult mice, and of their rates of turnover, we infer that the larger, relatively stable TCM and TEM populations in dirty mice (Figure 3A) must derive from their establishment in greater numbers in the first few weeks of life. The differences in compartment sizes in the two environments are then sustained well into adulthood by the very slow loss of these early memory populations.
To explore this hypothesis, we used the parameters estimated in adult mice to predict the development of their CD4+ MP T cell populations early in life. To do this, we drew on measurements of naive and central memory CD4+ T cell numbers recovered from the spleen and lymph nodes of WT mice aged between 5 days and 14 weeks, kept in the clean facility (Figure 5A). We then used the naive T cell timecourse with the parameters estimated for the best-fitting (age-dependent loss) model of TCM development in clean adult mice, to predict their accumulation of TCM from age 5d onwards (Figure 5B, left panel), starting from the mean numbers of CD4+ TCM observed at age 5d. The model slightly underestimated TCM numbers in clean adult mice and failed to capture their rapid accumulation up to age 4 weeks. In contrast, using the observed timecourse of TCM in clean WT mice as a source (Figure 5A, right panel), the establishment of the TEM compartment was predicted remarkably well by the adult parameters from the favoured two-phase model (Figure 5B, right panel).
To predict the early kinetics of CD4+ MP T cell populations in dirty mice, we first assumed their accumulation of naive T cells closely approximated that in clean mice, given that naive T cell numbers were similar in young adults from the two facilities (Figure 1B). We then used this timecourse (Figure 5A, left panel) with the parameters from the best-fitting (age-dependent loss) model of TCM development in adult dirty mice, to predict their kinetics of accumulation. This prediction underestimated TCM numbers at age 14 weeks by a factor of approximately 4 (Figure 5C, left panel). In turn, using this trajectory as the source of TEM, and using the force of recruitment and loss rates derived from adults, led to a similarly substantial underestimate of their numbers (Figure 5C, right panel). Predictions using the alternative models in all cases were even poorer; and all predictions were insensitive to the presumed numbers of TCM or TEM at age 5d, which are small and rapidly outnumbered by the influx of new memory cells from their precursor population.
We conclude that to account for memory T cell numbers in adulthood, mice in the clean facility experience a slightly elevated force of recruitment into TCM early in life; and this force is much larger in dirtier mice, presumably deriving from greater levels of exposure to environmental antigens.
Analysis of germ free mice confirms roles for both autoreactive and commensal stimuli in the establishment of memory compartments
We found that the rate of constitutive recruitment into memory in adult mice was insensitive to variations in environmental commensals, and that these antigens must exert their biggest influence on the establishment of MP cells in neonates and young mice. However, earlier studies reported that numbers of memory cells in the spleens of SPF and germ-free (GF) mice are similar and argue that self-recognition is therefore the sole driver of MP cell generation in early life (Kawabe et al., 2017). To reconcile these apparent differences, and dissect the contributions of self and foreign antigens to the establishment of MP T cells in young mice, we compared the size and behaviour of memory CD4+ T cell subsets in C57Bl6/J and/or C57Bl6/SJL mice housed in a wider range of environments. In addition to the clean (UCL) and dirty (NIMR) mice analysed above, we enumerated cells from GF and SPF mice obtained from the Kennedy Institute (KI) in Oxford. Consistent with these earlier studies, substantial numbers of both TCM and TEM MP cells were recovered from GF mice aged between 40 and 200d, confirming that the generation of MP cells does not depend exclusively on commensal-derived foreign antigens (Figure 6A). However, the memory CD4+ T cell compartments of GF mice, enumerated from spleen and lymph nodes combined, were significantly smaller than in the SPF mice in all facilities. Clean mice from KI and UCL had similar-sized memory compartments, and in turn both were substantially smaller than those in mice from the dirty facility (NIMR) (Figure 6A). Indeed, the dirty mice played host to around five times the number of MP cells found in GF mice, indicating that antigens from commensal organisms are a substantial driver of MP cell expansion.
We also analysed the proliferative activity of MP cells in mice from the different environments, by measuring the expression of Ki67. Division of MP cells is dependent on TCR (Seddon et al., 2003) and CD28 costimulation signaling (Kawabe et al., 2017). A substantial fraction of MP cells in GF mice expressed Ki67 (Figure 6B), indicating cell cycle activity within the previous 3–4 days (Gossel et al., 2017). This proportion was broadly similar to that in mice from both clean and dirty environments, indicating that the level of proliferation of CD4+ MP T cells in adult mice was relatively insensitive to environmentally derived stimuli.
Quantifying the forces exerted by commensals on memory generation early in life
Finally, we estimated the rates of memory generation in young mice in the different environments. We began with the favoured age-dependent loss model of TCM dynamics. Using the parameters from clean adult mice (which were similar to those estimated for dirty mice, and for which no estimates were available for GF mice), and the empirical description of CD4+ naive T cell numbers in clean WT mice (Figure 5A), we then estimated the fold changes in the force of recruitment needed during the first 11 weeks of life needed to seed CD4+ TCM at the average numbers observed in the mice aged between 10 and 28 weeks (Figure 7, left panel). CD4+ TCM numbers were relatively stable in all facilities during this period. GF mice needed approximately 0.4 times the force of recruitment in clean adult mice, younger clean mice needed a force approximately 1.1 times greater, and dirty mice required a 2.7-fold increase. As before, uncertainty in memory cell numbers at age 5d had very little effect on the predicted levels of memory attained at week 11, or on these estimated correction factors.
We then used these environment-specific, corrected trajectories of TCM development to predict the accumulation of TEM by age 11 weeks, using the favoured two-phase model. Remarkably, after accounting for the different TCM population sizes, the force of recruitment from TCM to TEM estimated in clean adult mice was also sufficient to account for TEM numbers in all three environments (Figure 7, right panel).
In summary, this analysis showed that approximately 2- to 3-fold increases or decreases in the force of recruitment into CD4+ TCM observed in clean adult (UCL) mice were sufficient to explain their numbers in dirty mice from NIMR or GF mice. However the subsequent rate of development of TEM from TCM appeared to be independent of both mouse age and environment, and differences in the numbers of TEM could be explained simply by the differences in the size of the TCM precursor population. These results suggest that the rate of generation of CD4+ TCM from naive T cells in young mice reflects both self antigens and the level of exposure to environmental antigens, but that the rate of differentiation from TCM to TEM is largely insensitive to these forces.
Discussion
In this study, we compared mice housed in facilities with distinct antigenic burdens to investigate the nature, magnitude and timing of the forces that establish and maintain CD4+ MP T cell compartments. We examined (i) tonic recruitment into the TCM and TEM pools in adults, (ii) the kinetic substructure and maintenance of these compartments throughout life, and (iii) their generation/establishment early in life. Our analyses indicate that self recognition contributes to all these processes, but that the contribution of reactivity to commensal antigens is largely restricted to the neonatal period.
Our analysis of GF mice, which lack exposure to commensal organisms, confirmed earlier work showing that generation of MP cells in adults is not driven exclusively by foreign commensals (Kawabe et al., 2017) and suggesting that self-recognition is instead the key driver. GF mice are not entirely free of environmental antigens and it is possible that proteins in bedding material and diet could represent foreign antigenic stimuli. However, the argument for self-recognition is made through the correlation of the degree of MP cell conversion and steady state proliferation with affinity for self-MHC, as indicated by CD5 expression. We also found that memory CD4+ T cell division, as reported by Ki67, was substantial and largely independent of the level of commensals. We could not determine the extents to which Ki67 expression derived from the homeostatic proliferation of existing cells or the influx of newly generated (and recently divided) MP cells. However, the common levels of Ki67 across environments, together with our findings from the adult clean and dirty chimeric mice that the rates of memory generation and net loss were insensitive to environment, support our conclusion that self-recognition is the major driver of both recruitment and proliferative renewal in adult memory compartments.
Commensals did, however, have a substantial impact upon the sizes of the memory compartments generated early in life. CD4+ TCM and TEM numbers in mice raised in dirty environments were 3–5 times greater than those in either GF or cleaner IVC or SPF facilities, and these differences could be explained by differences in the forces of recruitment of CD4+ TCM during ontogeny. The estimated seven-fold difference in this force between GF and dirty mice prompts the simple interpretation that self-recognition only accounts for ∼14% of the memory compartment in dirty mice.
We ascribe differences in memory compartment sizes to different commensal burdens, but it is possible that their smaller sizes in GF mice also derive from their smaller lymph nodes. Bacterial stimulation of DCs is required for their migration into lymph nodes, and these DCs are required for their normal development (Wendland et al., 2011; Moussion and Girard, 2011; Zhang et al., 2016). Therefore, it is not straightforward to separate the indirect influence of commensals on lymphoid development from any direct influences upon memory generation. It is also possible that some of the additional force of recruitment in neonates derives not from commensals but from lymphopenia, which can drive naive T cells to acquire a memory phenotype (Min et al., 2003). However, this process was demonstrated by transferring naive cells from adults into very young mice; it is possible that such naive cells do not represent the activity of neonatal naive T cells, which are almost exclusively recent thymic emigrants.
We observed that environment not only impacted memory but also naive T cells. Their numbers naturally decline with age, but this decline was almost twice as fast in dirty mice than in genetically identical mice housed in cleaner facilities (Figure 1—figure supplement 2A; numbers halve every 228d (95% CI 227–231) in clean mice, and 143d (142-144) in dirty mice). This difference likely derives from a more rapid reduction in thymic output with age in the dirty environment, and not differences in lifespans of CD4+ naive T cells in the two environments, because we see similar rates of decline in the numbers of single positive thymocytes at the latest stage of thymic development (halving every every 282d (281-308) in clean mice, 151d (150-159) in dirty mice; data not shown). Therefore, it appears that another consequence of life in a more antigenically diverse environment is more rapid involution of the thymus. It is possible that this effect derives from the stress or inflammation associated with an increased commensal burden, or is somehow a consequence of expanded memory populations; but whatever the mechanism, our data clearly indicate that environmental factors can impact the maintenance of naive T cells.
Our models do not incorporate any homeostatic regulation in the sense of modulation of rates of division or loss through quorum sensing. Since memory cell numbers only vary by a factor of 2–3 between clean and dirty mice, we did not expect to detect any strong variation in net loss rate with pool size with these data, and indeed our estimates of loss rates were similar in the two environments. While we cannot rule it out at higher cell densities, there is arguably little evidence for homeostatic regulation of circulating memory T cells in SPF mice. We observed a range of relatively stable memory compartment sizes in the different containment facilities; these memory compartments do not appear to fill up rapidly with large clones in very young mice, which might occur if division or loss rates are sensitive to total pool size; and they appear to be expandable following multiple infections in older animals (Vezys et al., 2009).
The models yielded robust conclusions regarding the nature and magnitudes of the forces generating CD4+ TCM and TEM, and the existence of heterogeneity within both subsets, consistent with previous reports by ourselves and others of subpopulations of CD4+ MP cells with distinct rates of division and turnover (Younes et al., 2011; Gossel et al., 2017). However, these models are abstractions, and resolving the details of kinetic substructure in lymphocyte populations is challenging (Ganusov et al., 2010). For CD4+ TCM in adult mice we found evidence for slow and continuous changes in their net loss rates as cells age, and additional support for such a process comes from the slow rate at which Ki67 levels in donor and host memory cells converge in adult busulfan chimeras. In contrast, the data for TEM in adult mice, and the predictions of their accumulation in younger mice, more strongly supported a discrete two-phase model with a relatively rapid transition from fast to slow memory. Our previous study of CD4+ TCM and TEM homeostasis in adult mice used short-term BrdU labelling to identify populations in both subsets that divide and die rapidly (Gossel et al., 2017), although in that study we assumed the proliferative and more quiescent pools were maintained independently and so it is not straightforward to compare the rate estimates with those presented here. Overall, although, it seems likely that both MP cell subsets are more heterogeneous than any one of our models suggests. Indeed, there are other potential sources of heterogeneity. One possibility is that MP T cells are generated with a distribution of net loss rates, and those clones with greater intrinsic fitness (lower net loss rates) are simply selected for over time. Such a mechanism – a generalisation of the simple ‘resistant memory’ model we rejected here – could explain the under-representation of donor cells in the memory compartments of busulfan chimeric mice, and may be difficult to distinguish from our model of gradual changes in fitness with cell age; indeed the two mechanisms are not mutually exclusive. Characterising the homeostatic dynamics of CD4+ memory T cells, and ultimately how these dynamics relate to their functional capacity, requires further study.
Materials and methods
Generating busulfan chimeric mice
Request a detailed protocolMice were treated with optimised low doses of busulfan to deplete HSC but leave peripheral T cell subsets intact. HSC were reconstituted with congenically labelled, T-cell depleted bone marrow to generate stable chimeras (Figure 1A). Details of the protocols are given in Hogan et al. (2017b) and Hogan et al. (2017a).
Mice
Busulfan chimeric mice and wild-type control mice were housed in conventional animal facilities, either at the National Institute for Medical Research, London, UK (NIMR); or at the UCL Royal Free Campus, London, UK (UCL). At NIMR, mice were housed in open cages and drank tap water. At UCL, mice were housed in individually ventilated cages and drank irradiated water. Germ Free and SPF mice were housed at the Oxford Centre for Microbiome Studies, Oxford, UK.
Flow cytometry
Request a detailed protocolSingle cell suspensions were prepared from the thymus, spleen and lymph nodes of busulfan chimeric mice, wildtype control mice, or germ free mice. Cells were stained with the following monoclonal antibodies and cell dyes: CD45.1 FITC, CD45.2 FITC, CD45.2 AlexaFluor700, TCR- APC, CD4+ PerCP-eFluor710, CD44 APC-eFluor780, CD25 PE, CD25 eFluor450, CD25 PE-Cy7, CD62L eFluor450, NK1.1 PE-Cy7 (all eBioscience), CD45.1 BV650, CD45.2 PE-Dazzle, TCR- PerCP-Cy5.5 CD4+ BV711, CD44 BV785, CD25 BV650 (all Biolegend), CD62L BUV737 (BD Biosciences), LIVE/DEAD nearIR and LIVE/DEAD blue viability dyes. For Ki67 staining, cells were fixed using the eBioscience Foxp3/Transcription Factor Staining Buffer Set and stained with either anti-mouse Ki67 FITC or PE (both eBioscience). Cells were acquired on a BD LSR-II or a BD LSR-Fortessa flow cytometer and analysed with Flowjo software (Treestar). Conventional CD4+ cells were identified as live TCR-+ CD4+ CD25- NK1.1-, and then CD44 and CD62L were used to identify EM (CD44+CD62L-) and CM (CD44+CD62L+) subsets.
Modelling the fluxes between naive, central memory and effector memory CD4+ subsets
The homogeneous model
Request a detailed protocolOur simplest description of the kinetics of the generation and renewal of CD4+ TCM and TEM is illustrated in Figure 2A and was formulated as follows. We assume that cells flow into a memory subset of total size from a precursor population at total rate , where is the age of the animal. The rate constant is the ‘force of recruitment’, a compound parameter which is the product of the per capita rate of recruitment of cells from the source population per day multiplied by number representing any net expansion that occurs during recruitment. Memory is also lost at net per capita rate . This rate is the balance of loss through death and/or differentiation, and any compensatory cell production through division. It represents the rate of decline or growth of a population that self-renews to any extent, rather than the loss rate of individual cells. We place no constraints on this rate, and so may be positive or negative.
We assume host and donor cells each obey the same kinetics, so that
where the total population size is . Our strategy for parameter estimation was to fit this model simultaneously to the timecourses of total numbers of memory cells , and the donor chimerism within memory, , which is the fraction of cells in memory that are donor-derived. For reasons detailed below, we normalise this fraction to the proportion of early double-positive (DP1) thymocytes that are donor-derived, which is measured in the same mouse and denoted ;
Equations (2) and (3) then give
If observations are made sufficiently long after BMT (more than 3–4 weeks), chimerism among DP thymocytes can be assumed to have stabilised (Hogan et al., 2015) and so we can neglect the term in . Then
where we define
By normalising the chimerism of both the source and the memory populations to that in the thymus, we remove any variation in these quantities due to variation across individuals in the degree of chimerism achieved with busulfan treatment and BMT.
Fitting required initial conditions for the total numbers of memory cells and donor chimerism. We solved all of the models from host age , which was chosen to be the time at which donor chimerism in memory for the mouse with the youngest age at BMT could be assumed to be zero and donor chimerism in DP1 had stabilised. This was host age 66d for the clean mice and 84d for dirty mice (26d and 28d post-BMT, respectively). Our results were insensitive to changes of a few days in these baseline ages. We also required functional forms for the kinetics of the immediate precursor (source) population and . When considering CD4+ naive T cells as a source, their numbers from age onwards in both facilities were well described with an exponential decay curve, , though with different exponents (Figure 1—figure supplement 2A). With this form, we can solve Equation (1) for explicitly;
Using this expression for in Equation (6), the kinetics of normalised chimerism in memory are
The rate of increase in donor chimerism in memory then depends on the force of recruitment , the dynamics and chimerism of the source ( and ), the initial memory pool size , and the net loss rate . Note that Equation (10) predicts that given sufficiently long, and if the rate of decline of naive T cells is less than the rate of loss of memory (), the chimerism in memory will stabilise at the chimerism of CD4+ naive T cells.
When considering TCM as a source for TEM in clean and dirty mice, we described with sigmoid and exponential decay functions, respectively (Figure 1—figure supplement 2B). We described each source’s chimerism with the generalised logistic function (Figure 1—figure supplement 2C and D).
The homogeneous model is characterised by the three unknowns , and . To estimate them for a given subset, location and source population we solved Equations (1) and (6) numerically and fitted them simultaneously to the timecourses of total numbers and normalised chimerism of the memory subset, using a method detailed below. The clean mice underwent BMT at a range of ages, which were accounted for individually in the fitting; model predictions for a mouse which underwent BMT at age and was observed at age were generated by running the model from host age 26d (clean) or 28d (dirty) to time ; with calculated from using Equation (8), and the normalised chimerism at time assumed to be zero.
The two-phase model
Request a detailed protocolThe two-phase model (Figure 2B) describes the kinetics of CD4+ TCM and TEM assuming that both comprise two subpopulations with distinct rates of loss,
where is the mouse age. We assume that cells flow only into one subset from the precursor population at total rate , and transition to the next compartment at rate ;
Note that despite the nomenclature, when estimating the rates of loss of these subsets, we did not constrain them; rather, the model fits indicated that .
Assuming that host and donor cells obey the same kinetics, so that Equation (12) hold identically for both populations, then similar to the derivation of Equation (6) we obtain the following equations for the dynamics of donor chimerism in the fast and slow subsets, each normalised to the chimerism of DP1 thymocytes;
where
The normalised chimerism in the fast and slow populations combined is
We dealt with different ages at BMT using the same approach described for the simplest model. We determined the initial sizes of the subsets and by assuming that fast cells were in quasi-equilibrium with their source, because all T cell populations change slowly in adult mice (Figure 1B); and allowing to be free. Allowing both population sizes to be free yielded very similar results, at the cost of an additional parameter. The numbers of host-derived cells in the fast and slow memory subsets at each time were then generated from and by running the model forward from age using Equation (12). The two-phase model is then characterised by five unknowns; , and , the transition rate , and the force of recruitment . To estimate these parameters, we fitted the solutions of Equations (11), (12), and (15) simultaneously to the timecourses of total memory cell numbers and the normalised chimerism , using the empirical forms of and where is host age and is time post-BMT. To visualise the fits to data from the clean facility, we partitioned the mice into four groups based on age at BMT, and plotted the model predictions for the median age at BMT within each group (Figure 3B).
The age-dependent loss model
Request a detailed protocolIn this model (Figure 2C), the loss rate is a function of the time since entry of a cell or its ancestor into memory, which we denote its age . The time evolution of the population density of memory cells of age at host age is given by
where the population density of cells of age zero is the rate at which cells flow into memory from the source,
and we must specify the overall population density with respect to cell age at some initial mouse age, . As with the other models, we assumed all cells present at are host-derived; we model their age distribution as . The free parameter could be positive or negative, such that older cells can initially be over- or under-represented compared to younger cells. This definition ensures is the rate of influx of cells of age zero from the source at time , . We explored exponential () and sigmoid () forms for the dependence of the net loss rate on cell age, with = 1, 2, 3, 5 and 10. We found that consistently yielded the best fits, with the exponential performing the most poorly.
Solving this system using the method of characteristics allows us to track the fates of three cell populations – the host-derived population that was present at (), and host- and donor-derived cells that entered the population after ( and ). Total memory cell numbers at are then
The terms in this expression evolve according to
where and . These expressions give
The normalised donor chimerism is
We fitted Equations (20) and (21) to the timecourses of their observed counterparts from host age onwards. This model has four free parameters; and , which, together with the observed value of , specify the initial age distribution of host cells, ; and and , which specify the form of . The parameters and then determine the initial number of host-derived memory cells;
As described above, this model was fitted simultaneously to data from mice who underwent BMT at different ages, replacing in Equation (21) with the age at BMT plus 26d or 28d for clean and dirty mice, respectively.
Resistant memory model
Request a detailed protocolIn this model, proposed in Gossel et al. (2017), the CD4+ TCM and TEM populations are assumed to be heterogeneous, each consisting of a ‘displaceable’ subset turning over at rate and continuously supplemented from the source, and an ‘incumbent’ or 'resistant' population of host memory cells, ; these are assumed to be established early in life, not supplemented thereafter, and have a distinct net loss rate :
All donor-derived cells are assumed to be displaceable and obey the same kinetics as displaceable host-derived cells. We solved Equation (23) to obtain total memory cell numbers , and the normalised chimerism in memory,
For simplicity, we assumed that resistant memory cells were stable in number (). The resistant memory model then has four free parameters (, , , ). Multiple ages at BMT were handled as described for the homogeneous model.
Parameter estimation and model selection
Request a detailed protocolEach model (with its parameter set denoted ) was fitted simultaneously to the timecourses of total memory cell numbers and the normalised chimerism , , for a given source population (naive/TCM) and environment (dirty/clean). The cell counts and chimerism values were log- and arcsin-square root-transformed, respectively, such that these new variables (denoted and ) could be assumed to have normally distributed errors with constant variances and . We then maximised the joint likelihood of and with respect to the model parameters , and the unknown and . If and are the model predictions of the transformed observations and respectively, this likelihood is
where SSR denotes the sum of squared residuals, and both SSRx and SSRy are functions of the parameters . The log-likelihood is then
To reduce the number of unknowns to be estimated with a parameter search, we substituted the maximum likelihood estimates of the error variances, which can be calculated directly;
(and similarly for ), giving the following expression for the joint log-likelihood,
This quantity was then maximised with respect to the parameters using the scipy.optimize package in Python. We used the best-fitting model predictions to estimate 95% confidence intervals on parameters by bootstrapping residuals 1000 times, re-fitting and taking the 2.5% and 97.5% quantiles of the resulting distributions of parameter estimates.
For each model and source, we performed the above procedure separately for data from the clean and dirty environments, and calculated a combined, maximum log likelihood . We then used the corrected Akaike Information Criterion, AICc (Akaike, 1974; Burnham and Anderson, 2002) to assess the relative support for each model/source pairing, where
Here, is the total number of estimated parameters, which was double the number of parameters in each model (one set for clean and another for dirty); and is the total number of observations, which was 2 × the number of mice in the clean facility + 2 × the number in the dirty facility (each mouse yielded two measurements for each memory subset – total cell numbers, and donor chimerism).
Annoted code and data for performing all analyses are freely available from https://github.com/marianowicka/memory-CD4-and-dirt.git (Nowicka, 2019; copy archived at https://github.com/elifesciences-publications/memory-CD4-and-dirt).
Data availability
The code for mathematical models has been deposited on GitHub (https://github.com/marianowicka/memory-CD4-and-dirt; copy archived at https://github.com/elifesciences-publications/memory-CD4-and-dirt), and raw cell counts used for model fitting are provided as source data files, as featured in figures 1, 5 and 7.
References
-
A new look at the statistical model identificationIEEE Transactions on Automatic Control 19:716–723.https://doi.org/10.1109/TAC.1974.1100705
-
T cell memory: heterogeneity and mechanismsClinical Immunology 95:173–181.https://doi.org/10.1006/clim.2000.4858
-
Generation and maintenance of immunological memorySeminars in Immunology 16:323–333.https://doi.org/10.1016/j.smim.2004.08.013
-
Alternative memory in the CD8 T cell lineageTrends in Immunology 32:50–56.https://doi.org/10.1016/j.it.2010.12.004
-
Endogenous proliferation: burst-like CD4 T cell proliferation in lymphopenic settingsSeminars in Immunology 17:201–207.https://doi.org/10.1016/j.smim.2005.02.005
Article and author information
Author details
Funding
National Institutes of Health (R01 AI093870)
- Thea Hogan
- Maria Nowicka
- Daniel Cownden
- Andrew J Yates
Medical Research Council (MR/P011225/1)
- Benedict Seddon
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
This work was supported by the NIH (R01 AI093870) and the Medical Research Council (MC-PC-13055). We thank Fiona Powrie and the Oxford Centre for Microbiome Studies for germ-free mice.
Ethics
Animal experimentation: Animal experiments were performed according to the UCL Animal Welfare and Ethical Review Body and Home Office regulations under PPL 70-8310.
Copyright
© 2019, Hogan et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 1,831
- views
-
- 293
- downloads
-
- 15
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Immunology and Inflammation
The abundance and biological contribution of natural killer (NK) cells in cancer are controversial. Here, we aim to uncover clinical relevance and cellular roles of NK cells in colon cancer liver metastasis (CCLM). Here, we integrated single-cell RNA-sequencing, spatial transcriptomics (ST), and bulk RNA-sequencing datasets to investigate NK cells’ biological properties and functions in the microenvironment of primary and liver metastatic tumors. Results were validated through an in vitro co-culture experiment based on bioinformatics analysis. Useing single-cell RNA-sequencing and ST, we mapped the immune cellular landscape of colon cancer and well-matched liver metastatic cancer. We discovered that GZMK+ resting NK cells increased significantly in tumor tissues and were enriched in the tumor regions of both diseases. After combining bulk RNA and clinical data, we observed that these NK cell subsets contributed to a worse prognosis. Meanwhile, KIR2DL4+ activated NK cells exhibited the opposite position and relevance. Pseudotime cell trajectory analysis revealed the evolution of activated to resting NK cells. In vitro experiments further confirmed that tumor-cell-co-cultured NK cells exhibited a decidual-like status, as evidenced by remarkable increasing CD9 expression. Functional experiments finally revealed that NK cells exhibited tumor-activating characteristics by promoting the dissociation of SCF (stem cell factor) on the tumor cells membrane depending on cell-to-cell interaction, as the supernatant of the co-culture system enhanced tumor progression. In summary, our findings revealed resting NK cells exhibited a clinical relevance with CCLM, which may be exploited for novel strategies to improve therapeutic outcomes for patients with CCLM.
-
- Immunology and Inflammation
- Neuroscience
Acute retinal ischemia and ischemia-reperfusion injury are the primary causes of retinal neural cell death and vision loss in retinal artery occlusion (RAO). The absence of an accurate mouse model for simulating the retinal ischemic process has hindered progress in developing neuroprotective agents for RAO. We developed a unilateral pterygopalatine ophthalmic artery occlusion (UPOAO) mouse model using silicone wire embolization combined with carotid artery ligation. The survival of retinal ganglion cells and visual function were evaluated to determine the duration of ischemia. Immunofluorescence staining, optical coherence tomography, and haematoxylin and eosin staining were utilized to assess changes in major neural cell classes and retinal structure degeneration at two reperfusion durations. Transcriptomics was employed to investigate alterations in the pathological process of UPOAO following ischemia and reperfusion, highlighting transcriptomic differences between UPOAO and other retinal ischemia-reperfusion models. The UPOAO model successfully replicated the acute interruption of retinal blood supply observed in RAO. 60 min of Ischemia led to significant loss of major retinal neural cells and visual function impairment. Notable thinning of the inner retinal layer, especially the ganglion cell layer, was evident post-UPOAO. Temporal transcriptome analysis revealed various pathophysiological processes related to immune cell migration, oxidative stress, and immune inflammation during the non-reperfusion and reperfusion periods. A pronounced increase in microglia within the retina and peripheral leukocytes accessing the retina was observed during reperfusion periods. Comparison of differentially expressed genes (DEGs) between the UPOAO and high intraocular pressure models revealed specific enrichments in lipid and steroid metabolism-related genes in the UPOAO model. The UPOAO model emerges as a novel tool for screening pathogenic genes and promoting further therapeutic research in RAO.