Comprehensive assessment and mathematical modeling of T cell population dynamics and homeostasis

Our current view of T cell differentiation and population dynamics is assembled from pieces of data obtained from separate experimental systems and is thus patchy. We reassessed homeostasis and dynamics of T cells 1) by generating a mathematical model describing the spatiotemporal features of T cell differentiation, and 2) by fitting this model to experimental data generated by disturbing T cell differentiation through transient depletion of dividing T cells in mice. This specific depletion was obtained by administration of ganciclovir to mice expressing the conditional thymidine kinase suicide gene in T cells. With this experimental approach, we could derive quantitative parameters describing the cell fluxes, residence times, and rates of import, export, proliferation, and death across cell compartments for thymocytes and recent thymic emigrants (RTEs). Among other parameters, we show that 93% of thymocytes produced before single-positive stages are eliminated through the selection process. Then, a postselection peripheral expansion of naive T cells contributes three times more to naive T cell production than the thymus, with half of the naive T cells consisting of dividing RTEs. Altogether, this work provides a quantitative population dynamical framework of thymocyte development, RTEs, and naive T cells.