steven.kleinstein@yale.edu.

Analyses of somatic hypermutation (SHM) patterns in B cell immunoglobulin (Ig) sequences contribute to our basic understanding of adaptive immunity, and have broad applications not only for understanding the immune response to pathogens, but also to determining the role of SHM in autoimmunity and B cell cancers. Although stochastic, SHM displays intrinsic biases that can confound statistical analysis, especially when combined with the particular codon usage and base composition in Ig sequences. Analysis of B cell clonal expansion, diversification, and selection processes thus critically depends on an accurate background model for SHM micro-sequence targeting (i.e., hot/cold-spots) and nucleotide substitution. Existing models are based on small numbers of sequences/mutations, in part because they depend on data from non-coding regions or non-functional sequences to remove the confounding influences of selection. Here, we combine high-throughput Ig sequencing with new computational analysis methods to produce improved models of SHM targeting and substitution that are based only on synonymous mutations, and are thus independent of selection. The resulting “S5F” models are based on 806,860 Synonymous mutations in 5-mer motifs from 1,145,182 Functional sequences and account for dependencies on the adjacent four nucleotides (two bases upstream and downstream of the mutation). The estimated profiles can explain almost half of the variance in observed mutation patterns, and clearly show that both mutation targeting and substitution are significantly influenced by neighboring bases. While mutability and substitution profiles were highly conserved across individuals, the variability across motifs was found to be much larger than previously estimated. The model and method source code are made available at

Individual variation in germline and expressed B-cell immunoglobulin (Ig) repertoires has been associated with aging, disease susceptibility, and differential response to infection and vaccination. Repertoire properties can now be studied at large-scale through next-generation sequencing of rearranged Ig genes. Accurate analysis of these repertoire-sequencing (Rep-Seq) data requires identifying the germline variable (V), diversity (D), and joining (J) gene segments used by each Ig sequence. Current V(D)J assignment methods work by aligning sequences to a database of known germline V(D)J segment alleles. However, existing databases are likely to be incomplete and novel polymorphisms are hard to differentiate from the frequent occurrence of somatic hypermutations in Ig sequences. Here we develop a Tool for Ig Genotype Elucidation via Rep-Seq (TIgGER). TIgGER analyzes mutation patterns in Rep-Seq data to identify novel V segment alleles, and also constructs a personalized germline database containing the specific set of alleles carried by a subject. This information is then used to improve the initial V segment assignments from existing tools, like IMGT/HighV-QUEST. The application of TIgGER to Rep-Seq data from seven subjects identified 11 novel V segment alleles, including at least one in every subject examined. These novel alleles constituted 13% of the total number of unique alleles in these subjects, and impacted 3% of V(D)J segment assignments. These results reinforce the highly polymorphic nature of human Ig V genes, and suggest that many novel alleles remain to be discovered. The integration of TIgGER into Rep-Seq processing pipelines will increase the accuracy of V segment assignments, thus improving B-cell repertoire analyses.next-generation sequencing | B-cell repertoire | adaptive immunity | somatic hypermutation | variable gene segment T he production by B cells of immunoglobulin (Ig) proteins, which are expressed on the cell surface as B-cell receptors and secreted by subsets of B cells as antibodies, is a key component of the adaptive immune system in humans. Through their specific binding to an enormously diverse range of foreign bodies, Ig proteins are able to elicit further immunological response and provide protection. These proteins are assembled in B cells from two pairs of polypeptide chains, termed heavy and light. The antigen-binding portions of these genes are created through the somatic recombination of gene segments, termed variable (V), diversity (D), and joining (J). During the recombination process, one each of the ∼46 V, 23 D, and 6 J gene segments (1) recombine to make the antigen-binding region of the heavy chain; the light chain is created by a similar process, although involving one of two different loci (λ and κ) containing V and J genes only. Over three million different Ig sequences can be created through this V(D)J recombinatorial process alone (2). The potential diversity of these sequences is further expanded to the order of trillions (2) when combined with the random insert...

The adaptive immune system confers protection by generating a diverse repertoire of antibody receptors that are rapidly expanded and contracted in response to specific targets. Next-generation DNA sequencing now provides the opportunity to survey this complex and vast repertoire. In the present work, we describe a set of tools for the analysis of antibody repertoires and their application to elucidating the dynamics of the response to viral vaccination in human volunteers. By analyzing data from 38 separate blood samples across 2 y, we found that the use of the germ-line library of V and J segments is conserved between individuals over time. Surprisingly, there appeared to be no correlation between the use level of a particular VJ combination and degree of expansion. We found the antibody RNA repertoire in each volunteer to be highly dynamic, with each individual displaying qualitatively different response dynamics. By using combinatorial phage display, we screened selected VH genes paired with their corresponding VL library for affinity against the vaccine antigens. Altogether, this work presents an additional set of tools for profiling the human antibody repertoire and demonstrates characterization of the fast repertoire dynamics through time in multiple individuals responding to an immune challenge.next-generation sequencing | influenza | immunology T he immune system is able to rapidly sense and respond to a vast array of invading organisms. Its arsenal contains systems that are immediately effective against commonly seen patterns (innate immunity) and systems that are capable of responding to novel invaders (adaptive immunity). Given the acute nature and diversity of infections, the immune system must be capable of rapid recognition of a pathogen, amplification of the response, and subsequent contraction of the response after the resolution of the infection. Adaptive immune responses rely on the continuous selection and amplification of specific clones from an enormous library of immune receptors (antibodies and T cell receptors). Specifically, stimulation of B-cell immunity results in the synthesis of antibodies that are secreted into the blood stream or into the mucosa as well as the programming of B memory cells that play a crucial role in the generation of rapid protective responses upon reinfection.Currently, many immunology studies depend on characterizing lymphocyte subsets (e.g., assaying cell-surface receptors) and the ability to correlate them to encoded genetic information (1). Recent advances in next-generation sequencing (NGS) (2) have enabled any DNA-encodable assay to produce massive amounts of data. Indeed, NGS has enabled unprecedented views into the immune repertoire, as its immune receptor diversity is genetically encoded within a complex collection of lymphocytes (3-8).The present study set out to dissect the rapid dynamics of the complete human peripheral antibody response against a controlled immune challenge (vaccination), without the a priori notion of cell state markers or functions...

The adaptive immune system's capability to protect the body requires a highly diverse lymphocyte antigen receptor repertoire. However, the influence of individual genetic and epigenetic differences on these repertoires is not typically measured. By leveraging the unique characteristics of B, CD4+ T and CD8+ T-lymphocyte subsets from monozygotic twins, we quantify the impact of heritable factors on both the V(D)J recombination process and on thymic selection. We show that the resulting biases in both V(D)J usage and N/P addition lengths, which are found in naïve and antigen experienced cells, contribute to significant variation in the CDR3 region. Moreover, we show that the relative usage of V and J gene segments is chromosomally biased, with ∼1.5 times as many rearrangements originating from a single chromosome. These data refine our understanding of the heritable mechanisms affecting the repertoire, and show that biases are evident on a chromosome-wide level.