Evidence for Shaping of Light Chain Repertoire by Structural Selection

Shared bias in H chain V-J pairing in naive and memory B cells

Introduction

H chain rearrangement in B cells is a two-step process where first DH binds JH, and only then VH is joined to the complex. As such, there is no direct rearrangement between VH and JH.

Results

Nevertheless, we here show that the VHJH combinations frequency in humans deviates from the one expected based on each gene usage frequency. This bias is observed mainly in functional rearrangements, and much less in out-of-frame rearrangements. The bias cannot be explained by preferred binding for DH genes or a preferred reading frame. Preferred VH JH combinations are shared between donors.

Discussion

These results suggest a common structural mechanism for these biases. Through development, thepreferred VH JH combinations evolve during peripheral selection to become stronger, but less shared. We propose that peripheral Heavy chain VH JH usage is initially shaped by a structural selection before the naive B cellstate, followed by pathogen-induced selection for host specific VH-JH pairs.

Autoencoder based local T cell repertoire density can be used to classify samples and T cell receptors

Recent advances in T cell repertoire (TCR) sequencing allow for the characterization of repertoire properties, as well as the frequency and sharing of specific TCR. However, there is no efficient measure for the local density of a given TCR. TCRs are often described either through their Complementary Determining region 3 (CDR3) sequences, or theirV/J usage, or their clone size. We here show that the local repertoire density can be estimated using a combined representation of these components through distance conserving autoencoders and Kernel Density Estimates (KDE). We present ELATE-an Encoder-based LocAl Tcr dEnsity and show that the resulting density of a sample can be used as a novel measure to study repertoire properties. The cross-density between two samples can be used as a similarity matrix to fully characterize samples from the same host. Finally, the same projection in combination with machine learning algorithms can be used to predict TCR-peptide binding through the local density of known TCRs binding a specific target.

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

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 AND IMPLEMENTATION

Source code is available at https://github.com/zsethna/OLGA.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genesis of the αβ T-cell receptor

The T-cell (TCR) repertoire relies on the diversity of receptors composed of two chains, called α and β, to recognize pathogens. Using results of high throughput sequencing and computational chain-pairing experiments of human TCR repertoires, we quantitively characterize the αβ generation process. We estimate the probabilities of a rescue recombination of the β chain on the second chromosome upon failure or success on the first chromosome. Unlike β chains, α chains recombine simultaneously on both chromosomes, resulting in correlated statistics of the two genes which we predict using a mechanistic model. We find that ∼35% of cells express both α chains. Altogether, our statistical analysis gives a complete quantitative mechanistic picture that results in the observed correlations in the generative process. We learn that the probability to generate any TCRαβ is lower than 10(-12) and estimate the generation diversity and sharing properties of the αβ TCR repertoire.