Consistency of VDJ rearrangement and substitution parameters enables accurate B cell receptor sequence annotation

Abstract

VDJ rearrangement and somatic hypermutation work together to produce antibody-coding B cell receptor (BCR) sequences for a remarkable diversity of antigens. It is now possible to sequence these BCRs in high throughput; analysis of these sequences is bringing new insight into how antibodies develop, in particular for broadly-neutralizing antibodies against HIV and influenza. A fundamental step in such sequence analysis is to annotate each base as coming from a specific one of the V, D, or J genes, or from an N-addition (a.k.a. non-templated insertion). Previous work has used simple parametric distributions to model transitions from state to state in a hidden Markov model (HMM) of VDJ recombination, and assumed that mutations occur via the same process across sites. However, codon frame and other effects have been observed to violate these parametric assumptions for such coding sequences, suggesting that a non-parametric approach to modeling the recombination process could be useful. In our paper, we find that indeed large modern data sets suggest a model using parameter-rich per-allele categorical distributions for HMM transition probabilities and per-allele-per-position mutation probabilities, and that using such a model for inference leads to significantly improved results. We present an accurate and efficient BCR sequence annotation software package using a novel HMM "factorization" strategy. This package, called partis (https://github.com/psathyrella/partis/), is built on a new general-purpose HMM compiler that can perform efficient inference given a simple text description of an HMM.

References in corpus (3)

• An introduction to Docker for reproducible research, with examples from the R environment
• Statistical inference of the generation probability of T-cell receptors from sequence repertoires
• Inferring processes underlying B-cell repertoire diversity

Cited by in corpus (16)

• IGoR: a tool for high-throughput immune repertoire analysis
• Computational strategies for dissecting the high-dimensional complexity of adaptive immune repertoires
• Predicting the spectrum of TCR repertoire sharing with a data-driven model of recombination
• Likelihood-based inference of B-cell clonal families
• Deep generative selection models of T and B cell receptor repertoires with soNNia
• Per-sample immunoglobulin germline inference from B cell receptor deep sequencing data
• Using genotype abundance to improve phylogenetic inference
• Quantitative Immunology for Physicists
• repgenHMM: a dynamic programming tool to infer the rules of immune receptor generation from sequence data
• Inference of B cell clonal families using heavy/light chain pairing information
• A Bayesian Phylogenetic Hidden Markov Model for B Cell Receptor Sequence Analysis
• Using B cell receptor lineage structures to predict affinity
• Predicting B Cell Receptor Substitution Profiles Using Public Repertoire Data
• SOS: Online probability estimation and generation of T and B cell receptors
• Survival analysis of DNA mutation motifs with penalized proportional hazards
• Fierce selection and interference in B-cell repertoire response to chronic HIV-1