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.
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.
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.
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).
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.
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).
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.
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.
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.
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.
Mice 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).
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.
Single 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.
Our 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 (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;
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).
In 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.
In 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.
Each 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).
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.
A new look at the statistical model identificationIEEE Transactions on Automatic Control 19:716–723.https://doi.org/10.1109/TAC.1974.1100705
BookModel Selection and Multimodel Inference (Second edition)Springer.
Generation and maintenance of immunological memorySeminars in Immunology 16:323–333.https://doi.org/10.1016/j.smim.2004.08.013
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
Softwarememory-CD4-and-dirt, version 115e88dGitHub.
Satyajit RathSenior Editor; Indian Institute of Science Education and Research (IISER), India
Rob J de BoerReviewing Editor; Utrecht University, Netherlands
Rob J de BoerReviewer; Utrecht University, Netherlands
Ruy RibeiroReviewer; Universidade de Lisboa, United States
Jose BorghansReviewer; UMCU, Netherlands
In the interests of transparency, eLife publishes the most substantive revision requests and the accompanying author responses.
In this paper Hogan et al. analyze the development of CD4+ T cell memory populations (central and effector) in mice housed under different conditions. They aim to understand the main driver of the establishment of memory population in the absence of infections, and provide insight into relative roles of antigens coming from commensals, self-antigens and pathogens in the establishment and maintenance of the memory T-cell pool throughout life of a mouse. A major finding is that after an early period where commensal / environmental antigens are important in driving memory, the replenishment of new memory has similar dynamics independent of the environment.
Decision letter after peer review:
Thank you for submitting your article "Differential impact of self and environmental antigens on the ontogeny and maintenance of CD4+ T cell memory" for consideration by eLife. Your article has been reviewed by three peer reviewers, including Rob de Boer as the Guest Editor and Reviewer #1, and the evaluation has been overseen by Satyajit Rath as the Senior Editor. The following individuals involved in review of your submission have agreed to reveal their identity: Ruy Ribeiro (Reviewer #2); Jose Borghans (Reviewer #3).
The reviewers have discussed the reviews with one another and the Guest Editor has drafted this decision to help you prepare a revised submission.
In this paper Hogan et al. analyze the development of CD4+ T cell memory populations (central and effector) in mice housed under different conditions. They aim to understand the main driver of the establishment of memory population in the absence of infections. This is an important and insightful paper that perfectly fits the scope of eLife. It is highly original work, very clearly written, and reports important insights into the relative roles of antigens coming from commensals, self-antigens and pathogens in the establishment and maintenance of the memory T-cell pool throughout life (of a mouse). It is also a very timely paper given the recent discussions in the immunological field on the role of exposure to "dirt". The main finding is that after an early period where commensal / environmental antigens are important in driving memory, after a few weeks continuous replenishment of new memory has similar dynamics independent of the environment where the mice are housed. The paper is well-written and clear.
The reviewers were particularly impressed with the novel insights into the maintenance and formation of CD4+ memory T cells in mice housed under various conditions.
1) The main question is regarding the modeling approach. We appreciate that the authors have tested three different models and chosen the best one. But all the models are linear, and so to achieve the balance observed between source, fast and slow cells requires fine tuning of the dynamical parameters. The question is whether it would make more biological sense to have some kind of feedback or density-dependent mechanism. Although we don't think you have to "invent a more complex model" when the simpler one works well, this could at least be a point of discussion.
2) Another question is that the authors model central and effector memory cells independently. Could these two cell populations be coupled, for example, through interconversion? How would this affect the modeling approach? And more importantly, could it affect the conclusions?https://doi.org/10.7554/eLife.48901.sa1
1) The main question is regarding the modeling approach. We appreciate that the authors have tested three different models and chosen the best one. But all the models are linear, and so to achieve the balance observed between source, fast and slow cells requires fine tuning of the dynamical parameters. The question is whether it would make more biological sense to have some kind of feedback or density-dependent mechanism. Although we don't think you have to "invent a more complex model" when the simpler one works well, this could at least be a point of discussion.
This is a nice point. We considered it but concluded that a treatment of feedback regulation of cell numbers will not be informative here, for three reasons;
1) The memory cell numbers we observe from age 10 weeks onwards only vary by factor of 2 or 3; the average pool sizes in clean and dirty mice differ by a similar factor (Figure 1B); and we don’t find significant differences in the average net loss rates of each subset in the two environments (Figure 4B and Table 2). Without larger perturbations to cell numbers, then, we do not expect to be able to identify any density dependence of the net loss rates with our data. For the two-compartment model, we are also already at the limit of the number of parameters we are comfortable attempting to estimate.
2) None of the models require fine-tuning to achieve stability. Numbers are simply the balance of slowly declining influx from the naive pool (constant force of recruitment * naive numbers, which decline naturally); and first order decay from one or more sequential compartments. So in these models all memory (sub)compartments can achieve a stable steady state without feedback regulation.
3) Arguably there is little experimental evidence for homeostatic regulation of memory numbers, at least in SPF mice — we see three levels of stable memory in the three facilities, memory pools appears not to fill up rapidly with large clones in very young mice; and memory cell numbers are expandable following multiple infections in older animals (Vezys, Yates et al. Nature 2009).
We have added a discussion of these points.
2) Another question is that the authors model central and effector memory cells independently. Could these two cell populations be coupled, for example, through interconversion? How would this affect the modeling approach? And more importantly, could it affect the conclusions?
In response to other comments we revisited the data with some alternative models (thank you for suggesting this) and in the process improved our descriptor functions for the sources and our strategy for parameter estimation; detailed in Materials and methods. As a result we now clearly favour a naive – CM – EM pathway.
We reached this conclusion by fitting naive -> CM and CM -> EM separately. We were not able to fit all three simultaneously – the number of parameters was too large, and given that we find heterogeneity within both pools there is some uncertainty in how to connect the two models. Even with a constant rate of conversion, solving two coupled models where survival changes continuously with cell age is technically beyond us. Also – since there may be an expansion factor between CM and EM, the two are only coupled weakly as regards fitting.
We did model a Naive-EM-CM pathway but the fits are extremely poor (see Author response image 1). Indeed we argue that such a pathway is inconsistent with the timecourses of chimerism in the two pools. So we choose not to include these fits.https://doi.org/10.7554/eLife.48901.sa2
- Thea Hogan
- Maria Nowicka
- Daniel Cownden
- Andrew J Yates
- Benedict Seddon
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
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.
Animal experimentation: Animal experiments were performed according to the UCL Animal Welfare and Ethical Review Body and Home Office regulations under PPL 70-8310.
- Satyajit Rath, Indian Institute of Science Education and Research (IISER), India
- Rob J de Boer, Utrecht University, Netherlands
- Rob J de Boer, Utrecht University, Netherlands
- Ruy Ribeiro, Universidade de Lisboa, United States
- Jose Borghans, UMCU, Netherlands
© 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.
Partial/complete pathologic response following neoadjuvant chemotherapy (NAC) in pancreatic cancer (PDAC) patients undergoing pancreatectomy is associated with improved survival. We sought to determine whether neutrophil-to-lymphocyte ratio (NLR) dynamics predict pathologic response following chemotherapy in PDAC, and if manipulating NLR impacts chemosensitivity in preclinical models and uncovers potential mechanistic underpinnings underlying these effects.
Pathologic response in PDAC patients (n=94) undergoing NAC and pancreatectomy (7/2015-12/2019) was dichotomized as partial/complete or poor/absent. Bootstrap-validated multivariable models assessed associations between pre-chemotherapy NLR (%neutrophils÷%lymphocytes) or NLR dynamics during chemotherapy (ΔNLR = pre-surgery—pre-chemotherapy NLR) and pathologic response, disease-free survival (DFS), and overall survival (OS). To preclinically model effects of NLR attenuation on chemosensitivity, Ptf1aCre/+; KrasLSL-G12D/+;Tgfbr2flox/flox (PKT) mice and C57BL/6 mice orthotopically injected with KrasLSL-G12D/+;Trp53LSL-R172H/+;Pdx1Cre(KPC) cells were randomized to vehicle, gemcitabine/paclitaxel alone, and NLR-attenuating anti-Ly6G with/without gemcitabine/paclitaxel treatment.
In 94 PDAC patients undergoing NAC (median:4 months), pre-chemotherapy NLR (p<0.001) and ΔNLR attenuation during NAC (p=0.002) were independently associated with partial/complete pathologic response. An NLR score = pre-chemotherapy NLR+ΔNLR correlated with DFS (p=0.006) and OS (p=0.002). Upon preclinical modeling, combining NLR-attenuating anti-Ly6G treatment with gemcitabine/paclitaxel—compared with gemcitabine/paclitaxel or anti-Ly6G alone—not only significantly reduced tumor burden and metastatic outgrowth, but also augmented tumor-infiltrating CD107a+-degranulating CD8+ T-cells (p<0.01) while dampening inflammatory cancer-associated fibroblast (CAF) polarization (p=0.006) and chemoresistant IL-6/STAT-3 signaling in vivo. Neutrophil-derived IL-1β emerged as a novel mediator of stromal inflammation, inducing inflammatory CAF polarization and CAF-tumor cell IL-6/STAT-3 signaling in ex vivo co-cultures.
Therapeutic strategies to mitigate neutrophil-CAF-tumor cell IL-1β/IL-6/STAT-3 signaling during NAC may improve pathologic responses and/or survival in PDAC.
Supported by KL2 career development grant by Miami CTSI under NIH Award UL1TR002736, Stanley Glaser Foundation, American College of Surgeons Franklin Martin Career Development Award, and Association for Academic Surgery Joel J. Roslyn Faculty Award (to J. Datta); NIH R01 CA161976 (to N.B. Merchant); and NCI/NIH Award P30CA240139 (to J. Datta and N.B. Merchant).
During obesity and high fat-diet (HFD) feeding in mice, sustained low-grade inflammation includes not only increased pro-inflammatory macrophages in the expanding adipose tissue, but also bone marrow (BM) production of invasive Ly6Chigh monocytes. As BM adiposity also accrues with HFD, we explored the relationship between the gains in BM white adipocytes and invasive Ly6Chigh monocytes by in vivo and ex vivo paradigms. We find a temporal and causal link between BM adipocyte whitening and the Ly6Chigh monocyte surge, preceding the adipose tissue macrophage rise during HFD in mice. Phenocopying this, ex vivo treatment of BM cells with conditioned media from BM adipocytes or bona fide white adipocytes favoured Ly6Chigh monocyte preponderance. Notably, Ly6Chigh skewing was preceded by monocyte metabolic reprogramming towards glycolysis, reduced oxidative potential and increased mitochondrial fission. In sum, short-term HFD changes BM cellularity, resulting in local adipocyte whitening driving a gradual increase and activation of invasive Ly6Chigh monocytes.