A branching stochastic evolutionary model of the B-cell repertoire

Abstract

We propose a stochastic framework to describe the evolution of the B-cell repertoire during germinal center (GC) reactions. Our model is formulated as a multitype age-dependent branching process with time-varying immigration. The immigration process captures the mechanism by which founder B cells initiate clones by gradually seeding GC over time, while the branching process describes the temporal evolution of the composition of these clones. The model assigns a type to each cell to represent attributes of interest. Examples of attributes include the binding affinity class of the B cells, their clonal family, or the nucleotide sequence of the heavy and light chains of their receptors. The process is generally non-Markovian. We present its properties, including as t→∞\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$t\\rightarrow \\infty $$\\end{document} when the process is supercritical, the most relevant case to study expansion of GC B cells. We introduce temporal alpha and beta diversity indices for multitype branching processes. We focus on the dynamics of clonal dominance, highlighting its non-stationarity, and the accumulation of somatic hypermutations in the context of sequential immunization. We evaluate the impact of the ongoing seeding of GC by founder B cells on the dynamics of the B-cell repertoire, and quantify the effect of precursor frequency and antigen availability on the timing of GC entry. An application of the model illustrates how it may help with interpretation of BCR sequencing data.