  The human body is able to generate 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.
