How heterogeneous thymic output and homeostatic proliferation shape naive T cell receptor clone abundance distributions.

The set of T cells that express the same T cell receptor (TCR) sequence represent a T cell clone. The number of different naive T cell clones in an organism reflects the number of different T cell receptors (TCRs) arising from recombination of the V(D)J gene segments during T cell development in the thymus. TCR diversity and more specifically, the clone abundance distribution is an important factor in immune function. Specific recombination patterns occur more frequently than others while subsequent interactions between TCRs and self-antigens are known to trigger proliferation and sustain naive T cell survival. These processes are TCR-dependent, leading to clone-dependent thymic export and naive T cell proliferation rates. Using a mean-field approximation to the solution of a regulated birth-death-immigration model, we systematically quantify how TCR-dependent heterogeneities in immigration and proliferation rates affect the shape of clone abundance distributions (the number of different clones that are represented by a specific number of cells). By comparing predicted clone abundances derived from our heterogeneous birth-death-immigration model with experimentally sampled clone abundances, we quantify the heterogeneity necessary to generate the observed abundances. Our findings indicate that heterogeneity in proliferation rates is more likely the mechanism underlying the observed clone abundance distributions than heterogeneity in immigration rates.