TCR beta allelic exclusion in dynamical models of V(D)J recombination based on allele independence

Allelic exclusion represents a major aspect of TCRbeta gene assembly by V(D)J recombination in developing T lymphocytes. Despite recent progress, its comprehension remains problematic when confronted with experimental data. Existing models fall short in terms of incorporating into a unique distribution all the cell subsets emerging from the TCRbeta assembly process. To revise this issue, we propose dynamical, continuous-time Markov chain-based modeling whereby essential steps in the biological procedure (D-J and V-DJ rearrangements and feedback inhibition) evolve independently on the two TCRbeta alleles in every single cell while displaying random modes of initiation and duration. By selecting parameters via fitting procedures, we demonstrate the capacity of the model to offer accurate fractions of all distinct TCRbeta genotypes observed in studies using developing and mature T cells from wild-type or mutant mice. Selected parameters in turn afford relative duration for each given step, hence updating TCRbeta recombination distinctive timings. Overall, our dynamical modeling integrating allele independence and noise in recombination and feedback-inhibition events illustrates how the combination of these ingredients alone may enforce allelic exclusion at the TCRbeta locus.