Abstract:Affinity maturation (AM) of B cells through somatic hypermutations (SHMs) enables the immune system to evolve to recognize diverse pathogens. The accumulation of SHMs leads to the formation of clonal lineages of antibody-secreting b cells that have evolved from a common naïve B cell. Advances in high-throughput sequencing have enabled deep scans of B cell receptor repertoires, paving the way for reconstructing clonal trees. However, it is not clear if clonal trees, which capture microevolutionary time scales, can be reconstructed using traditional phylogenetic reconstruction methods with adequate accuracy. In fact, several clonal tree reconstruction methods have been developed to fix supposed shortcomings of phylogenetic methods. Nevertheless, no consensus has been reached regarding the relative accuracy of these methods, partially because evaluation is challenging. Benchmarking the performance of existing methods and developing better methods would both benefit from realistic models of clonal lineage evolution specifically designed for emulating B cell evolution. In this paper, we propose a model for modeling B cell clonal lineage evolution and use this model to benchmark several existing clonal tree reconstruction methods. Our model, designed to be extensible, has several features: by evolving the clonal tree and sequences simultaneously, it allows modeling selective pressure due to changes in affinity binding; it enables scalable simulations of large numbers of cells; it enables several rounds of infection by an evolving pathogen; and, it models building of memory. In addition, we also suggest a set of metrics for comparing clonal trees and measuring their properties. Our results show that while maximum likelihood phylogenetic reconstruction methods can fail to capture key features of clonal tree expansion if applied naively, a simple post-processing of their results, where short branches are contracted, leads to inferences that are better than alternative methods.

Fluctuating fitness shapes the clone size distribution of immune repertoires

Population dynamics of immune repertoires

Affinity maturation for an optimal balance between long-term immune coverage and short-term resource constraints

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

Molecular tug of war reveals adaptive potential of an immune cell repertoire

The dynamic fitness landscape of ageing haematopoiesis through clonal competition

A dynamical model of the adaptive immune system: effects of cells promiscuity, antigens and B-B interactions

Identifying high power breakdowns in accelerating structures with acoustic sensors

Trait-space patterning and the role of feedback in antigen-immunity coevolution

Diversity against adversity: How adaptive immunity evolves potent antibodies

Regulation of T cell expansion by antigen presentation dynamics

Antibody Dynamics Simulation-A Mathematical Exploration of Clonal Deletion and Somatic Hypermutation.

Model-based assessment of the Role of Uneven Partitioning of Molecular Content on Heterogeneity and Regulation of Differentiation in CD8 T-cell Immune Responses

Assessing T Cell Clonal Size Distribution: A Non-Parametric Approach

A scalable model for simulating multi-round antibody evolution and benchmarking of clonal tree reconstruction methods

Emergence of clonal selection and affinity maturation in an ab initio microscopic model of immunity

Mathematical Characterization of Private and Public Immune Repertoire Sequences

Phenotypic determinism and stochasticity in antibody repertoires of clonally expanded plasma cells

Modeling Repertoire Behavior within the Energy-Based Framework: A Study on Cell Proliferation and Somatic Hypermutation during the Adaptive Immune Response

A branching stochastic evolutionary model of the B-cell repertoire

Persisting fetal clonotypes influence the structure and overlap of adult human T cell receptor repertoires