Stochastic modeling of T cell receptor gamma gene rearrangement

The mechanisms controlling the recombination process of the gamma genes that encode the gamma chain of the antigen receptor of the gammadelta T lymphocytes are unclear. Based on experimental data on the recombination status of the two major TCR gamma genes expressed in V(gamma)4+ and V(gamma)1+ thymocytes, we tested the plausibility of three possible rearrangement mechanisms: (1) a time window mechanism according to which the two chromosomes are accessible to the recombination machinery during a defined period of time; (2) a feedback mechanism in which recombination stops shortly after the first in-frame rearrangement event anywhere in both chromosomes; and (3) a feedback mechanism with asynchronous chromosome accessibility, in which there is a first period when only one chromosome is accessible for recombination, followed by a second period when both chromosomes are accessible; shortly after the first in-frame rearrangement event, during any of these two periods, recombination will definitely stop. We model the time window mechanism using a pure probabilistic approach and the two feedback mechanisms using a continuous-time Markov chain formalism. We used maximum likelihood methodology to infer the goodness-of-fit of the models showing evidence for the last model, which best fits the data. Further analysis of this model suggests an evolutionary tradeoff between allelic and isotypic exclusion and the probability that a precursor differentiates into a mature gammadelta T lymphocyte.