Converging evolution leads to near maximal junction diversity through parallel mechanisms in B and T cell receptors

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.

Counting is almost all you need

The immune memory repertoire encodes the history of present and past infections and immunological attributes of the individual. As such, multiple methods were proposed to use T-cell receptor (TCR) repertoires to detect disease history. We here show that the counting method outperforms two leading algorithms. We then show that the counting can be further improved using a novel attention model to weigh the different TCRs. The attention model is based on the projection of TCRs using a Variational AutoEncoder (VAE). Both counting and attention algorithms predict better than current leading algorithms whether the host had CMV and its HLA alleles. As an intermediate solution between the complex attention model and the very simple counting model, we propose a new Graph Convolutional Network approach that obtains the accuracy of the attention model and the simplicity of the counting model. The code for the models used in the paper is provided at: https://github.com/louzounlab/CountingIsAlmostAllYouNeed.

Two Step Selection for Bias in β Chain V-J Pairing

The β chain rearrangement in T cells is a two-step process where first Dβ and Jβ bind, and only then Vβ is joined to the complex. We here show that the frequency of human and mouse Vβ Jβ combinations deviates from the one expected based on each gene usage frequency. This bias is observed mainly in functional (F) rearrangements, but also slightly in non-functional (NF) rearrangements. Preferred Vβ Jβ combinations in F clones are shared between donors and samples, suggesting a common structural mechanism for these biases in addition to any host-specific antigen-induced peripheral selection. The sharing holds even in clones with Jβ1 that share the same Dβ1 gene. Vβ Jβ usage is correlated with the Molecular Weight and Isoelectric Point in F clones. The pairing is also observed in the Double Positive cells in mice thymocytes, suggesting that the selection leading to such a pairing occurs before thymic selection. These results suggest an additional structural checkpoint in the beta chain development prior to thymic selection during the T cell receptor expression. Understanding this structural selection is important for the distinction between normal and aberrant T cell development, and crucial for the design of engineered TCRs.

An Analysis of the Effects of Spaceflight and Vaccination on Antibody Repertoire Diversity

Ab repertoire diversity plays a critical role in the host's ability to fight pathogens. CDR3 is partially responsible for Ab-Ag binding and is a significant source of diversity in the repertoire. CDR3 diversity is generated during VDJ rearrangement because of gene segment selection, gene segment trimming and splicing, and the addition of nucleotides. We analyzed the Ab repertoire diversity across multiple experiments examining the effects of spaceflight on the Ab repertoire after vaccination. Five datasets from four experiments were analyzed using rank-abundance curves and Shannon indices as measures of diversity. We discovered a trend toward lower diversity as a result of spaceflight but did not find the same decrease in our physiological model of microgravity in either the spleen or bone marrow. However, the bone marrow repertoire showed a reduction in diversity after vaccination. We also detected differences in Shannon indices between experiments and tissues. We did not detect a pattern of CDR3 usage across the experiments. Overall, we were able to find differences in the Ab repertoire diversity across experimental groups and tissues.

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.

Using High-Throughput Sequencing to Characterize the Development of the Antibody Repertoire During Infections: A Case Study of HIV-1

High throughput sequencing (HTS) approaches have only recently been applied to describing the antibody/B-cell repertoire in fine detail, but these data sets have already become critical to the design of vaccines and therapeutics, and monitoring of cancer immunotherapy. As a case study, we describe the potential and present limitations of HTS studies of the Ab repertoire during infection with HIV-1. Most of the present studies restrict their analyses to lineages of specific bnAbs. We discuss future initiatives to expand this type of analysis to more complete repertoires and to improve comparing and sharing of these Ab repertoire data across studies and institutions.

Evidence for Shaping of Light Chain Repertoire by Structural Selection

The naïve immunoglobulin (IG) repertoire in the blood differs from the direct output of the rearrangement process. These differences stem from selection that affects the germline gene usage and the junctional nucleotides. A major complication obscuring the details of the selection mechanism in the heavy chain is the failure to properly identify the D germline and determine the nucleotide addition and deletion in the junction region. The selection affecting junctional diversity can, however, be studied in the light chain that has no D gene. We use probabilistic and deterministic models to infer and disentangle generation and selection of the light chain, using large samples of light chains sequenced from healthy donors and transgenic mice. We have previously used similar models for the beta chain of T-cell receptors and the heavy chain of IGs. Selection is observed mainly in the CDR3. The CDR3 length and mass distributions are narrower after selection than before, indicating stabilizing selection for mid-range values. Within the CDR3, proline and cysteine undergo negative selection, while glycine undergoes positive selection. The results presented here suggest structural selection maintaining the size of the CDR3 within a limited range, and preventing turns in the CDR3 region.