Abstract:The human body generates a diverse set of high affinity antibodies, the soluble form of B cell receptors (BCRs), that bind to and neutralize invading pathogens. The natural development of BCRs must be understood in order to design vaccines for highly mutable pathogens such as influenza and HIV. BCR diversity is induced by naturally occurring combinatorial "V(D)J" rearrangement, mutation, and selection processes. Most current methods for BCR sequence analysis focus on separately modeling the above processes. Statistical phylogenetic methods are often used to model the mutational dynamics of BCR sequence data, but these techniques do not consider all the complexities associated with B cell diversification such as the V(D)J rearrangement process. In particular, standard phylogenetic approaches assume the DNA bases of the progenitor (or "naive") sequence arise independently and according to the same distribution, ignoring the complexities of V(D)J rearrangement. In this paper, we introduce a novel approach to Bayesian phylogenetic inference for BCR sequences that is based on a phylogenetic hidden Markov model (phylo-HMM). This technique not only integrates a naive rearrangement model with a phylogenetic model for BCR sequence evolution but also naturally accounts for uncertainty in all unobserved variables, including the phylogenetic tree, via posterior distribution sampling.Rational vaccine design efforts depend on accurate inference of full evolutionary paths from a given naive sequence to the corresponding mature B cell receptor sequences in a germinal center. Before one can perform ancestral sequence inference for clonal sequences that result from the same naive rearrangement event, one must first obtain an estimate of the clonal phylogenetic tree. Currently, standard phylogenetic inference techniques are used to model the process of sequence evolution along the tree; however, these methods do not account for all the complexities associated with this evolutionary process. In this paper, we propose a Bayesian approach to phylogenetic inference for clonal sequences that is based on a phylogenetic hidden Markov model. Our phylo-HMM models both the naive rearrangement and somatic hypermutation processes and this Bayesian framework allows us to naturally account for uncertainty in all unobserved variables, including a phylogenetic tree, via posterior distribution sampling. We perform simulation-based experiments to show that naive sequence and phylogenetic inference performed jointly provides higher-quality estimates than those obtained by considering these inferences separately. Our application to real data reveals significant uncertainty in naive and ancestral sequences, confirming the importance of a Bayesian approach.

Combining mutation and recombination statistics to infer clonal families in antibody repertoires

Learning the heterogeneous hypermutation landscape of immunoglobulins from high-throughput repertoire data

Evaluating methods for B-cell clonal family assignment

Somatic hypermutation analysis for improved identification of B cell clonal families from next-generation sequencing data

Reconstructing a B-cell clonal lineage. I. Statistical inference of unobserved ancestors

Hierarchical Clustering Can Identify B Cell Clones with High Confidence in Ig Repertoire Sequencing Data

Inferring processes underlying B-cell repertoire diversity

Exploring the impact of clonal definition on B-cell diversity: implications for the analysis of immune repertoires

Clustering-based Identification of Clonally-Related Immunoglobulin Gene Sequence Sets

Systematic evaluation of B-cell clonal family inference approaches

Highly Sensitive And Unbiased Approach For Elucidating Antibody Repertoires

Using B cell receptor lineage structures to predict affinity

Applying phylogenetic methods for species delimitation to distinguish B-cell clonal families

A Bayesian phylogenetic hidden Markov model for B cell receptor sequence analysis

Statistical analysis of repertoire data demonstrates the influence of microhomology in V(D)J recombination

Consistency of VDJ Rearrangement and Substitution Parameters Enables Accurate B Cell Receptor Sequence Annotation

Predicting B cell receptor substitution profiles using public repertoire data

Inferring the immune response from repertoire sequencing

IGHV allele similarity clustering improves genotype inference from adaptive immune receptor repertoire sequencing data

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

Learning the statistics and landscape of somatic mutation-induced insertions and deletions in antibodies