BACKGROUND

Regulatory T cells (Tregs) were shown to be central in maintaining immunological homeostasis and preventing the development of autoimmune diseases. Several subsets of Tregs have been identified to date; however, the dynamics of the interactions between these subsets, and their implications on their regulatory functions are yet to be elucidated.

METHODOLOGY/PRINCIPAL FINDINGS: We employed a combination of mathematical modeling and frequent in vivo measurements of several T cell subsets. Healthy BALB/c mice received a single injection of either hCDR1--a tolerogenic peptide previously shown to induce Tregs, a control peptide or vehicle alone, and were monitored for 16 days. During this period, splenocytes from the treated mice were analyzed for the levels of CD4, CD25, CD8, CD28 and Foxp3. The collected data were then fitted to mathematical models, in order to test competing hypotheses regarding the interactions between the followed T cell subsets. In all 3 treatment groups, a significant, lasting, non-random perturbation of the immune system could be observed. Our analysis predicted the emergence of functional CD4 Tregs based on inverse oscillations of the latter and CD4(+)CD25(-) cells. Furthermore, CD4 Tregs seemed to require a sufficiently high level of CD8 Tregs in order to become functional, while conversion was unlikely to be their major source. Our results indicated in addition that Foxp3 is not a sufficient marker for regulatory activity.

CONCLUSIONS/SIGNIFICANCE: In this work, we unraveled the dynamics of the interplay between CD4, CD8 Tregs and effector T cells, using, for the first time, a mathematical-mechanistic perspective in the analysis of Treg kinetics. Furthermore, the results obtained from this interdisciplinary approach supported the notion that CD4 Tregs need to interact with CD8 Tregs in order to become functional. Finally, we generated predictions regarding the time-dependent function of Tregs, which can be further tested empirically in future work.