OLGA: fast computation of generation probabilities of B- and T-cell receptor amino acid sequences and motifs

Abstract

Motivation: High-throughput sequencing of large immune repertoires has enabled the development of methods to predict the probability of generation by V(D)J recombination of T- and B-cell receptors of any specific nucleotide sequence. These generation probabilities are very non-homogeneous, ranging over 20 orders of magnitude in real repertoires. Since the function of a receptor really depends on its protein sequence, it is important to be able to predict this probability of generation at the amino acid level. However, brute-force summation over all the nucleotide sequences with the correct amino acid translation is computationally intractable. The purpose of this paper is to present a solution to this problem. Results: We use dynamic programming to construct an efficient and flexible algorithm, called OLGA (Optimized Likelihood estimate of immunoGlobulin Amino-acid sequences), for calculating the probability of generating a given CDR3 amino acid sequence or motif, with or without V/J restriction, as a result of V(D)J recombination in B or T cells. We apply it to databases of epitope-specific T-cell receptors to evaluate the probability that a typical human subject will possess T cells responsive to specific disease-associated epitopes. The model prediction shows an excellent agreement with published data. We suggest that OLGA may be a useful tool to guide vaccine design. Availability: Source code is available at https://github.com/zsethna/OLGA

References in corpus (10)

• Statistical inference of the generation probability of T-cell receptors from sequence repertoires
• IGoR: a tool for high-throughput immune repertoire analysis
• Inferring processes underlying B-cell repertoire diversity
• Replicate immunosequencing as a robust probe of B cell repertoire diversity
• Predicting the spectrum of TCR repertoire sharing with a data-driven model of recombination
• Precise tracking of vaccine-responding T-cell clones reveals convergent and personalized response in identical twins
• Persisting fetal clonotypes influence the structure and overlap of adult human T cell receptor repertoires
• Method for identification of condition-associated public antigen receptor sequences
• Mouse T cell repertoires as statistical ensembles: overall characterization and age dependence
• repgenHMM: a dynamic programming tool to infer the rules of immune receptor generation from sequence data

Cited by in corpus (13)

• Detecting T-cell receptors involved in immune responses from single repertoire snapshots
• Deep generative selection models of T and B cell receptor repertoires with soNNia
• How many different clonotypes do immune repertoires contain?
• Augmenting adaptive immunity: progress and challenges in the quantitative engineering and analysis of adaptive immune receptor repertoires
• Quantitative Immunology for Physicists
• Population variability in the generation and thymic selection of T-cell repertoires
• Early life imprints the hierarchy of T cell clone sizes
• Diversity in Biology: definitions, quantification, and models
• Limits on Inferring T-cell Specificity from Partial Information
• Immune Fingerprinting through Repertoire Similarity
• On generative models of T-cell receptor sequences
• SOS: Online probability estimation and generation of T and B cell receptors
• How heterogeneous thymic output and homeostatic proliferation shape naive T cell receptor clone abundance distributions