A mathematical model of the effects of aging on naive T-cell population and diversity

The human adaptive immune response is known to weaken in advanced age, resulting in increased severity of pathogen-born illness, poor vaccine efficacy, and a higher prevalence of cancer in the elderly. Age-related erosion of the T-cell compartment has been implicated as a likely cause, but the underlying mechanisms driving this immunosenescence have not been quantitatively modeled and systematically analyzed. T-cell receptor diversity, or the extent of pathogen-derived antigen responsiveness of the T-cell pool, is known to diminish with age, but inherent experimental difficulties preclude accurate analysis on the full organismal level. In this paper, we formulate a mechanistic mathematical model of T-cell population dynamics on the immunoclonal subpopulation level, which provides quantitative estimates of diversity. We define different estimates for diversity that depend on the individual number of cells in a specific immunoclone. We show that diversity decreases with age primarily due to diminished thymic output of new T-cells and the resulting overall loss of small immunoclones.

💡 Research Summary

The paper presents a mechanistic mathematical framework to quantify how aging impacts the naive T‑cell compartment and its receptor diversity, addressing a central aspect of immunosenescence. Beginning with a concise review of clinical observations—reduced vaccine efficacy, heightened infection severity, and increased cancer incidence in the elderly—the authors highlight that loss of naive T‑cells and the concomitant decline in T‑cell receptor (TCR) structural diversity are thought to underlie these phenomena. Experimental measurement of full‑organism TCR diversity is hampered by sampling limitations and sequencing errors, motivating a theoretical approach.

The core of the model is a birth‑death‑migration ordinary differential equation for the total naive T‑cell count N(t):  dN/dt = γ(t) + p N(t) – μ(N) N(t), where γ(t) = γ₀ e^{–a t} captures the exponentially decaying thymic output, p is the per‑cell peripheral proliferation rate, and μ(N) = μ₀ + μ₁ N²/(N²+K²) represents a density‑dependent death rate that incorporates IL‑7 mediated homeostatic regulation (μ₀ basal death, μ₁ maximal IL‑7‑induced increase, K a carrying‑capacity‑like constant). Six parameters (γ₀, a, p, μ₀, μ₁, K) are calibrated against published human and rodent data; non‑dimensionalisation reduces the dynamics to three key composite parameters: γ₀ a^{–1} K^{–1}, (p–μ₀) a^{–1}, and μ₁/(p–μ₀).

Numerical exploration reveals four qualitatively distinct regimes. When μ₁/(p–μ₀) < 1, proliferation outpaces death, leading to unbounded growth—a scenario deemed unrealistic for healthy aging. If (p–μ₀) a^{–1} < 0, death dominates and the population collapses. The biologically plausible region satisfies μ₁/(p–μ₀) ≥ 1 and (p–μ₀) a^{–1} > 0. Within this region, two sub‑regimes emerge: (i) a proliferation‑driven mode where γ₀ a^{–1} K^{–1} ≪ 1, causing N(t) to approach a steady‑state N_ss determined solely by the balance p = μ(N_ss); (ii) a thymus‑driven mode where γ₀ a^{–1} K^{–1} is larger, producing an early peak in N(t) followed by a gradual decline toward the same N_ss as thymic output wanes. Parameter choices consistent with human data (γ₀≈1.8×10¹⁰, a≈0.044 yr⁻¹, K≈10¹⁰–10¹¹, p≈0.018 day⁻¹, μ₀≈0.17 day⁻¹, μ₁≈0.04) reproduce the empirically observed ~50 % reduction in total naive T‑cell numbers between ages 20–30 and 70–80.

To address diversity, the authors introduce a clone‑level description. Each clone i carries n_i(t) cells; the distribution P(n,t) evolves under the same birth‑death‑migration processes. Diversity is quantified using Rényi entropy‑based indices D_q = (∑_i n_i^q)^{1/(1–q)} (q≠1), encompassing species richness (q→0), Shannon entropy (q→1), and Simpson’s index (q=2). Simulations show that aging disproportionately eliminates small clones (n_i = 1–2), which are the primary contributors to overall diversity. Consequently, D_0, D_1, and D_2 all decline to roughly 60‑80 % of their youthful values by age 70–80, even though the total cell count only halves. This finding underscores that loss of rare clones, rather than sheer cell number, drives the collapse of TCR repertoire breadth.

The paper also examines the impact of limited sampling. Real‑world studies typically analyze a minute blood fraction, leading to severe under‑estimation of true diversity. The authors compare their model‑based “ground truth” with common estimators (Chao1, ACE) and demonstrate that, because the clone‑size distribution is highly skewed, these estimators remain biased unless sampling depth is dramatically increased.

In the discussion, the authors argue that the model clarifies why elderly individuals exhibit reduced responsiveness to novel antigens: the dwindling pool of rare clones limits the repertoire’s ability to recognize unfamiliar epitopes. They suggest that therapeutic strategies aimed at boosting thymic output (e.g., IL‑7 supplementation, thymic rejuvenation) or preserving small clones could be evaluated within this quantitative framework. Moreover, the model can be extended to incorporate functional diversity (cross‑reactivity) and stochastic effects of clonal expansion during infections.

Overall, the study provides a rigorous, parameter‑grounded mathematical description of naive T‑cell dynamics across the human lifespan, linking mechanistic changes (thymic involution, homeostatic proliferation, density‑dependent death) to observable declines in both cell numbers and receptor diversity. It offers a valuable tool for immunologists and clinicians to predict the consequences of aging on adaptive immunity and to design interventions that may mitigate immunosenescence.