Sequence intrinsic somatic mutation mechanisms contribute to affinity maturation of VRC01-class HIV-1 broadly neutralizing antibodies

Mutation-guided vaccine design: A process for developing boosting immunogens for HIV broadly neutralizing antibody induction

A major goal of HIV-1 vaccine development is the induction of broadly neutralizing antibodies (bnAbs). Although success has been achieved in initiating bnAb B cell lineages, design of boosting immunogens that select for bnAb B cell receptors with improbable mutations required for bnAb affinity maturation remains difficult. Here, we demonstrate a process for designing boosting immunogens for a V3-glycan bnAb B cell lineage. The immunogens induced affinity-matured antibodies by selecting for functional improbable mutations in bnAb precursor knockin mice. Moreover, we show similar success in prime and boosting with nucleoside-modified mRNA-encoded HIV-1 envelope trimer immunogens, with improved selection by mRNA immunogens of improbable mutations required for bnAb binding to key envelope glycans. These results demonstrate the ability of both protein and mRNA prime-boost immunogens for selection of rare B cell lineage intermediates with neutralizing breadth after bnAb precursor expansion, a key proof of concept and milestone toward development of an HIV-1 vaccine.

DNA flexibility can shape the preferential hypermutation of antibody genes

Antibody-coding genes accumulate somatic mutations to achieve antibody affinity maturation. Genetic dissection using various mouse models has shown that intrinsic hypermutations occur preferentially and are predisposed in the DNA region encoding antigen-contacting residues. The molecular basis of nonrandom/preferential mutations is a long-sought question in the field. Here, we summarize recent findings on how single-strand (ss)DNA flexibility facilitates activation-induced cytidine deaminase (AID) activity and fine-tunes the mutation rates at a mesoscale within the antibody variable domain exon. We propose that antibody coding sequences are selected based on mutability during the evolution of adaptive immunity and that DNA mechanics play a noncoding role in the genome. The mechanics code may also determine other cellular DNA metabolism processes, which awaits future investigation.

Breakthrough infection elicits hypermutated IGHV3-53/3-66 public antibodies with broad and potent neutralizing activity against SARS-CoV-2 variants including the emerging EG.5 lineages

The rapid emergence of SARS-CoV-2 variants of concern (VOCs) calls for efforts to study broadly neutralizing antibodies elicited by infection or vaccination so as to inform the development of vaccines and antibody therapeutics with broad protection. Here, we identified two convalescents of breakthrough infection with relatively high neutralizing titers against all tested viruses. Among 50 spike-specific monoclonal antibodies (mAbs) cloned from their B cells, the top 6 neutralizing mAbs (KXD01-06) belong to previously defined IGHV3-53/3-66 public antibodies. Although most antibodies in this class are dramatically escaped by VOCs, KXD01-06 all exhibit broad neutralizing capacity, particularly KXD01-03, which neutralize SARS-CoV-2 from prototype to the emerging EG.5.1 and FL.1.5.1. Deep mutational scanning reveals that KXD01-06 can be escaped by current and prospective variants with mutations on D420, Y421, L455, F456, N460, A475 and N487. Genetic and functional analysis further indicates that the extent of somatic hypermutation is critical for the breadth of KXD01-06 and other IGHV3-53/3-66 public antibodies. Overall, the prevalence of broadly neutralizing IGHV3-53/3-66 public antibodies in these two convalescents provides rationale for novel vaccines based on this class of antibodies. Meanwhile, KXD01-06 can be developed as candidates of therapeutics against SARS-CoV-2 through further affinity maturation.

Strategies for HIV-1 vaccines that induce broadly neutralizing antibodies

After nearly four decades of research, a safe and effective HIV-1 vaccine remains elusive. There are many reasons why the development of a potent and durable HIV-1 vaccine is challenging, including the extraordinary genetic diversity of HIV-1 and its complex mechanisms of immune evasion. HIV-1 envelope glycoproteins are poorly recognized by the immune system, which means that potent broadly neutralizing antibodies (bnAbs) are only infrequently induced in the setting of HIV-1 infection or through vaccination. Thus, the biology of HIV-1-host interactions necessitates novel strategies for vaccine development to be designed to activate and expand rare bnAb-producing B cell lineages and to select for the acquisition of critical improbable bnAb mutations. Here we discuss strategies for the induction of potent and broad HIV-1 bnAbs and outline the steps that may be necessary for ultimate success.

Ancestral reconstruction reveals catalytic inactivation of activation-induced cytidine deaminase concomitant with cold water adaption in the Gadiformes bony fish

BACKGROUND

Antibody affinity maturation in vertebrates requires the enzyme activation-induced cytidine deaminase (AID) which initiates secondary antibody diversification by mutating the immunoglobulin loci. AID-driven antibody diversification is conserved across jawed vertebrates since bony and cartilaginous fish. Two exceptions have recently been reported, the Pipefish and Anglerfish, in which the AID-encoding aicda gene has been lost. Both cases are associated with unusual reproductive behavior, including male pregnancy and sexual parasitism. Several cold water fish in the Atlantic cod (Gadinae) family carry an aicda gene that encodes for a full-length enzyme but lack affinity-matured antibodies and rely on antibodies of broad antigenic specificity. Hence, we examined the functionality of their AID.

RESULTS

By combining genomics, transcriptomics, immune responsiveness, and functional enzymology of AID from 36 extant species, we demonstrate that AID of that Atlantic cod and related fish have extremely lethargic or no catalytic activity. Through ancestral reconstruction and functional enzymology of 71 AID enzymes, we show that this enzymatic inactivation likely took place relatively recently at the emergence of the true cod family (Gadidae) from their ancestral Gadiformes order. We show that this AID inactivation is not only concordant with the previously shown loss of key adaptive immune genes and expansion of innate and cell-based immune genes in the Gadiformes but is further reflected in the genomes of these fish in the form of loss of AID-favored sequence motifs in their immunoglobulin variable region genes.

CONCLUSIONS

Recent demonstrations of the loss of the aicda gene in two fish species challenge the paradigm that AID-driven secondary antibody diversification is absolutely conserved in jawed vertebrates. These species have unusual reproductive behaviors forming an evolutionary pressure for a certain loss of immunity to avoid tissue rejection. We report here an instance of catalytic inactivation and functional loss of AID rather than gene loss in a conventionally reproducing vertebrate. Our data suggest that an expanded innate immunity, in addition to lower pathogenic pressures in a cold environment relieved the pressure to maintain robust secondary antibody diversification. We suggest that in this unique scenario, the AID-mediated collateral genome-wide damage would form an evolutionary pressure to lose AID function.

A multi-objective based clustering for inferring BCR clonal lineages from high-throughput B cell repertoire data

The adaptive B cell response is driven by the expansion, somatic hypermutation, and selection of B cell clonal lineages. A high number of clonal lineages in a B cell population indicates a highly diverse repertoire, while clonal size distribution and sequence diversity reflect antigen selective pressure. Identifying clonal lineages is fundamental to many repertoire studies, including repertoire comparisons, clonal tracking, and statistical analysis. Several methods have been developed to group sequences from high-throughput B cell repertoire data. Current methods use clustering algorithms to group clonally-related sequences based on their similarities or distances. Such approaches create groups by optimizing a single objective that typically minimizes intra-clonal distances. However, optimizing several objective functions can be advantageous and boost the algorithm convergence rate. Here we propose MobiLLe, a new method based on multi-objective clustering. Our approach requires V(D)J annotations to obtain the initial groups and iteratively applies two objective functions that optimize cohesion and separation within clonal lineages simultaneously. We show that our method greatly improves clonal lineage grouping on simulated benchmarks with varied mutation rates compared to other tools. When applied to experimental repertoires generated from high-throughput sequencing, its clustering results are comparable to the most performing tools and can reproduce the results of previous publications. The method based on multi-objective clustering can accurately identify clonally-related antibody sequences and presents the lowest running time among state-of-art tools. All these features constitute an attractive option for repertoire analysis, particularly in the clinical context. MobiLLe can potentially help unravel the mechanisms involved in developing and evolving B cell malignancies.

High-throughput sequencing can produce a large set of sequences and has profoundly changed our ability to study immune repertoires, particularly B cell receptor sequences. An important application is the analysis of the clonal lineage composition of B cell populations; it is the starting point of many immune repertoire studies, for instance, to differentiate between healthy individuals and those with lymphoid malignancies or other diseases. Several computational methods have been developed to identify clonal lineages from a set of B cell receptor sequences. Most of them apply clustering algorithms and optimize a single objective function that typically minimizes intra-clonal distances. However, optimizing several objective functions in parallel can benefit and increase the clustering performance and efficiency. We propose MobiLLe, the first multi-objective clonal lineage grouping method, which simultaneously optimizes two objective functions for minimizing intra-clonal diversity and maximizing inter-clonal differences. Our approach greatly improved clonal grouping on simulated benchmarks and performed comparably to the most powerful and recent methods on experimental samples. MobiLLe is computationally more efficient than existing tools and does not require any training process or hyper-parameter optimization. It can easily manage large-scale experimental repertoires, providing useful plots to help researchers detect clonally-related sequences in high-throughput B cell repertoire data.

Contribution of rare mutational outcomes to broadly neutralizing antibodies

Antibodies are important immune molecules that are elicited by B cells to protect our bodies during viral infections or vaccinations. In humans, the antibody repertoire is diversified by programmed DNA lesion processes to ensure specific and high affinity binding to various antigens. Broadly neutralizing antibodies (bnAbs) are antibodies that have strong neutralizing activities against different variants of a virus. bnAbs such as anti-HIV bnAbs often have special characteristics including insertions and deletions, long complementarity determining region 3 (CDR3), and high frequencies of mutations, often at improbable sites of the variable regions. These unique features are rare mutational outcomes that are acquired during antibody diversification processes. In this review, we will discuss possible mechanisms that generate these rare antibody mutational outcomes. The understanding of the mechanisms that generate these rare mutational outcomes during antibody diversification will have implications in vaccine design strategies to elicit bnAbs.

Learning the statistics and landscape of somatic mutation-induced insertions and deletions in antibodies

Affinity maturation is crucial for improving the binding affinity of antibodies to antigens. This process is mainly driven by point substitutions caused by somatic hypermutations of the immunoglobulin gene. It also includes deletions and insertions of genomic material known as indels. While the landscape of point substitutions has been extensively studied, a detailed statistical description of indels is still lacking. Here we present a probabilistic inference tool to learn the statistics of indels from repertoire sequencing data, which overcomes the pitfalls and biases of standard annotation methods. The model includes antibody-specific maturation ages to account for variable mutational loads in the repertoire. After validation on synthetic data, we applied our tool to a large dataset of human immunoglobulin heavy chains. The inferred model allows us to identify universal statistical features of indels in heavy chains. We report distinct insertion and deletion hotspots, and show that the distribution of lengths of indels follows a geometric distribution, which puts constraints on future mechanistic models of the hypermutation process.

Affinity maturation of B cell receptors is an important mechanism by which our body designs neutralizing antibodies to defend us against pathogens, including broadly neutralizing antibodies, which target a wide range of variants of the same pathogen. Such antibodies often contain key insertions and deletions to the germline gene, or “indels”, which are caused by somatic hypermutations. However, the mechanism, frequency and role of these indels are still elusive. We designed a computational method based on a probabilistic framework to infer the characteristics of this mutational process from high-throughput antibody sequencing experiments. Applied to human data, our approach provides a comprehensive quantitative description of insertions and deletions, opening avenues for better understanding the process of affinity maturation and the design of vaccines for eliciting a broad antibody response.

Complementary Roles of Antibody Heavy and Light Chain Somatic Hypermutation in Conferring Breadth and Potency to the HIV-1-Specific CAP256-VRC26 bNAb Lineage

Some HIV-infected people develop broadly neutralizing antibodies (bNAbs) that block many diverse, unrelated strains of HIV from infecting target cells and, through passive immunization, protect animals and humans from infection. Therefore, understanding the development of bNAbs and their neutralization can inform the design of an HIV vaccine. Here, we extend our previous studies of the ontogeny of the CAP256-VRC26 V2-targeting bNAb lineage by defining the mutations that confer neutralization to the unmutated common ancestor (CAP256.UCA). Analysis of the sequence of the CAP256.UCA showed that many improbable mutations were located in the third complementarity-determining region of the heavy chain (CDRH3) and the heavy chain framework 3 (FR3). Transferring the CDRH3 from bNAb CAP256.25 (63% breadth and 0.003 μg/mL potency) into the CAP256.UCA introduced breadth and the ability to neutralize emerging viral variants. In addition, we showed that the framework and light chain contributed to potency and that the second CDR of the light chain forms part of the paratope of CAP256.25. Notably, a minimally mutated CAP256 antibody, with 41% of the mutations compared to bNAb CAP256.25, was broader (64% breadth) and more potent (0.39 μg/mL geometric potency) than many unrelated bNAbs. Together, we have identified key regions and mutations that confer breadth and potency in a V2-specific bNAb lineage. These data indicate that immunogens that target affinity maturation to key sites in CAP256-VRC26-like precursors, including the CDRHs and light chain, could rapidly elicit breadth through vaccination. IMPORTANCE A major focus in the search for an HIV vaccine is elucidating the ontogeny of broadly neutralizing antibodies (bNAbs), which prevent HIV infection in vitro and in vivo. The unmutated common ancestors (UCAs) of bNAbs are generally strain specific and acquire breadth through extensive, and sometimes redundant, somatic hypermutation during affinity maturation. We investigated which mutations in the CAP256-VRC26 bNAb lineage conferred neutralization capacity to the UCA. We found that mutations in the antibody heavy and light chains had complementary roles in neutralization breadth and potency, respectively. The heavy chain, particularly the third complementarity-determining region, was responsible for conferring breadth. In addition, previously uninvestigated mutations in the framework also contributed to breadth. Together, approximately half of the mutations in CAP256.25 were necessary for broader and more potent neutralization than many unrelated neutralizing antibodies. Vaccine approaches that promote affinity maturation at key sites could therefore more rapidly produce antibodies with neutralization breadth.

Improving antibody affinity through <i>in vitro</i> mutagenesis in complementarity determining regions

High-affinity antibodies are widely used in diagnostics and for the treatment of human diseases. However, most antibodies are isolated from semi-synthetic libraries by phage display and do not possess in vivo affinity maturation, which is triggered by antigen immunization. It is therefore necessary to engineer the affinity of these antibodies by way of in vitro assaying. In this study, we optimized the affinity of two human monoclonal antibodies which were isolated by phage display in a previous related study. For the 42A1 antibody, which targets the liver cancer antigen glypican-3, the variant T57H in the second complementarity-determining region of the heavy chain (CDR-H2) exhibited a 2.6-fold improvement in affinity, as well as enhanced cell-binding activity. For the I4A3 antibody to severe acute respiratory syndrome coronavirus 2, beneficial single mutations in CDR-H2 and CDR-H3 were randomly combined to select the best synergistic mutations. Among these, the mutation S53P-S98T improved binding affinity (about 3.7 fold) and the neutralizing activity (about 12 fold) compared to the parent antibody. Taken together, single mutations of key residues in antibody CDRs were enough to increase binding affinity with improved antibody functions. The mutagenic combination of key residues in different CDRs creates additive enhancements. Therefore, this study provides a safe and effective in vitro strategy for optimizing antibody affinity.

Rapid selection of HIV envelopes that bind to neutralizing antibody B cell lineage members with functional improbable mutations

Elicitation of broadly neutralizing antibodies (bnAbs) by an HIV vaccine will involve priming the immune system to activate antibody precursors, followed by boosting immunizations to select for antibodies with functional features required for neutralization breadth. The higher the number of acquired mutations necessary for function, the more convoluted are the antibody developmental pathways. HIV bnAbs acquire a large number of somatic mutations, but not all mutations are functionally important. In this study, we identify a minimal subset of mutations sufficient for the function of the naturally occurring V3-glycan bnAb DH270.6. Using antibody library screening, candidate envelope immunogens that interact with DH270.6-like antibodies containing this set of key mutations are identified and selected in vitro. Our results demonstrate that less complex B cell evolutionary pathways than those naturally observed exist for the induction of HIV bnAbs by vaccination, and they establish rational approaches to identify boosting candidate immunogens.

The unique biology of germinal center B cells

Germinal center (GC) B cells are the source of the high-affinity, class-switched antibodies required for protective immunity. The unique biology of GC B cells involves iterative rounds of antibody gene somatic hypermutation coupled to multiple selection and differentiation pathways. Recent advances in areas such as single cell and gene editing technologies have shed new light upon these complex and dynamic processes. We review these findings here and integrate them into the current understanding of GC B cell replication and death, the retention of high-affinity and class-switched B cells in the GC, and differentiation into plasma and memory cell effectors. We also discuss how the biology of GC responses relates to vaccine effectiveness and outline current and future challenges in the field.

In vitro evolution of antibody affinity via insertional scanning mutagenesis of an entire antibody variable region

We report a systematic combinatorial exploration of affinity enhancement of antibodies by insertions and deletions (InDels). Transposon-based introduction of InDels via the method TRIAD (transposition-based random insertion and deletion mutagenesis) was used to generate large libraries with random in-frame InDels across the entire single-chain variable fragment gene that were further recombined and screened by ribosome display. Knowledge of potential insertion points from TRIAD libraries formed the basis of exploration of length and sequence diversity of novel insertions by insertional-scanning mutagenesis (InScaM). An overall 256-fold affinity improvement of an anti-IL-13 antibody BAK1 as a result of InDel mutagenesis and combination with known point mutations validates this approach, and suggests that the results of this InDel mutagenesis and conventional exploration of point mutations can synergize to generate antibodies with higher affinity.

AID Overlapping and Polη Hotspots Are Key Features of Evolutionary Variation Within the Human Antibody Heavy Chain (IGHV) Genes

Somatic hypermutation (SHM) of the immunoglobulin variable (IgV) loci is a key process in antibody affinity maturation. The enzyme activation-induced deaminase (AID), initiates SHM by creating C → U mismatches on single-stranded DNA (ssDNA). AID has preferential hotspot motif targets in the context of WRC/GYW (W = A/T, R = A/G, Y = C/T) and particularly at WGCW overlapping hotspots where hotspots appear opposite each other on both strands. Subsequent recruitment of the low-fidelity DNA repair enzyme, Polymerase eta (Polη), during mismatch repair, creates additional mutations at WA/TW sites. Although there are more than 50 functional immunoglobulin heavy chain variable (IGHV) segments in humans, the fundamental differences between these genes and their ability to respond to all possible foreign antigens is still poorly understood. To better understand this, we generated profiles of WGCW hotspots in each of the human IGHV genes and found the expected high frequency in complementarity determining regions (CDRs) that encode the antigen binding sites but also an unexpectedly high frequency of WGCW in certain framework (FW) sub-regions. Principal Components Analysis (PCA) of these overlapping AID hotspot profiles revealed that one major difference between IGHV families is the presence or absence of WGCW in a sub-region of FW3 sometimes referred to as “CDR4.” Further differences between members of each family (e.g., IGHV1) are primarily determined by their WGCW densities in CDR1. We previously suggested that the co-localization of AID overlapping and Polη hotspots was associated with high mutability of certain IGHV sub-regions, such as the CDRs. To evaluate the importance of this feature, we extended the WGCW profiles, combining them with local densities of Polη (WA) hotspots, thus describing the co-localization of both types of hotspots across all IGHV genes. We also verified that co-localization is associated with higher mutability. PCA of the co-localization profiles showed CDR1 and CDR2 as being the main contributors to variance among IGHV genes, consistent with the importance of these sub-regions in antigen binding. Our results suggest that AID overlapping (WGCW) hotspots alone or in conjunction with Polη (WA/TW) hotspots are key features of evolutionary variation between IGHV genes.

VRC34-Antibody Lineage Development Reveals How a Required Rare Mutation Shapes the Maturation of a Broad HIV-Neutralizing Lineage

Rare mutations have been proposed to restrict the development of broadly neutralizing antibodies against HIV-1, but this has not been explicitly demonstrated. We hypothesized that such rare mutations might be identified by comparing broadly neutralizing and non-broadly neutralizing branches of an antibody-developmental tree. Because sequences of antibodies isolated from the fusion peptide (FP)-targeting VRC34-antibody lineage suggested it might be suitable for such rare mutation analysis, we carried out next-generation sequencing (NGS) on B cell transcripts from donor N123, the source of the VRC34 lineage, and functionally and structurally characterized inferred intermediates along broadly neutralizing and poorly neutralizing developmental branches. The broadly neutralizing VRC34.01 branch required the rare heavy-chain mutation Y33P to bind FP, whereas the early bifurcated VRC34.05 branch did not require this rare mutation and evolved less breadth. Our results demonstrate how a required rare mutation can restrict development and shape the maturation of a broad HIV-1-neutralizing antibody lineage.

Targeted selection of HIV-specific antibody mutations by engineering B cell maturation

INTRODUCTION

A major goal of HIV-1 vaccine development is the design of immunogens that induce broadly neutralizing antibodies (bnAbs). However, vaccination of humans has not resulted in the induction of affinity-matured and potent HIV-1 bnAbs. To devise effective vaccine strategies, we previously determined the maturation pathway of select HIV-1 bnAbs from acute infection through neutralizing antibody development. During their evolution, bnAbs acquire an abundance of improbable amino acid substitutions as a result of nucleotide mutations at variable region sequences rarely targeted by activation-induced cytidine deaminase, the enzyme responsible for antibody mutation. A subset of improbable mutations is essential for broad neutralization activity, and their acquisition represents a key roadblock to bnAb development.

RATIONALE

Current bnAb lineage-based vaccine strategies can initiate bnAb lineage development in animal models but have not specifically elicited the improbable mutations required for neutralization breadth. Induction of bnAbs requires vaccine strategies that specifically engage bnAb precursors and subsequently select for improbable mutations required for broadly neutralizing activity. We hypothesized that vaccination with immunogens that bind with moderate to high affinity to bnAb B cell precursors, and with higher affinity to precursors that have acquired improbable mutations, could initiate bnAb B cell lineages and select for key improbable mutations required for bnAb development.

RESULTS

We elicited serum neutralizing HIV-1 antibodies in human bnAb precursor knock-in mice and wild-type macaques vaccinated with immunogens designed to select for improbable mutations. We designed two HIV-1 envelope immunogens that bound precursor B cells of either a CD4 binding site or V3-glycan bnAb lineage. In vitro, these immunogens bound more strongly to bnAb precursors once the precursor acquired the desired improbable mutations. Vaccination of macaques with the CD4 binding site–targeting immunogen induced CD4 binding site serum neutralizing antibodies. Antibody sequences elicited in human bnAb precursor knock-in mice encoded functional improbable mutations critical for bnAb development. In bnAb precursor knock-in mice, we isolated a vaccine-elicited monoclonal antibody bearing functional improbable mutations that was capable of neutralizing multiple HIV-1 global isolates. Structures of a bnAb precursor, a bnAb, and the vaccine-elicited antibody revealed the precise roles that acquired improbable mutations played in recognizing the HIV-1 envelope. Thus, our immunogens elicited antibody responses in macaques and knock-in mice that exhibited the mutational patterns, structural characteristics, or neutralization profiles of nascent broadly neutralizing antibodies.

CONCLUSION

Our study represents a proof of concept for targeted selection of improbable mutations to guide antibody affinity maturation. Moreover, this study demonstrates a rational strategy for sequential immunogen design to circumvent the difficult roadblocks in HIV-1 bnAb induction by vaccination. We show that immunogens should exhibit differences in affinity across antibody maturation stages where improbable mutations are necessary for the desired antibody function. This strategy of selection of specific antibody nucleotides by immunogen design can be applied to B cell lineages targeting other pathogens where guided affinity maturation is needed for a protective antibody response.

Intrinsic Nucleotide Preference of Diversifying Base Editors Guides Antibody Ex Vivo Affinity Maturation

Base editors (BEs) are emerging tools used for precision correction or diversifying mutation. It provides a potential way to recreate somatic hypermutations (SHM) for generating high-affinity antibody, which is usually screened from antigen-challenged animal models or synthetic combinatorial libraries. By comparing somatic mutations in the same genomic context, we screened engineered deaminases and CRISPR-deaminase coupling approaches and updated diversifying base editors (DBEs) to generate SHM. The deaminase used in DBEs retains its intrinsic nucleotide preference and mutates cytidines at its preferred motifs. DBE with AID targets the same hotspots as physiological AID does in vivo, while DBE with other deaminases generates distinct mutation profiles from the same DNA substrate. Downstream DNA repair pathways further diversified the sequence, while Cas9-nickase restricted mutation spreading. Finally, application of DBE in an antibody display system achieved antibody affinity maturation ex vivo. Our findings provide insight of DBE working mechanism and an alternative antibody engineering approach.

A spectral clustering-based method for identifying clones from high-throughput B cell repertoire sequencing data

Motivation

B cells derive their antigen-specificity through the expression of Immunoglobulin (Ig) receptors on their surface. These receptors are initially generated stochastically by somatic re-arrangement of the DNA and further diversified following antigen-activation by a process of somatic hypermutation, which introduces mainly point substitutions into the receptor DNA at a high rate. Recent advances in next-generation sequencing have enabled large-scale profiling of the B cell Ig repertoire from blood and tissue samples. A key computational challenge in the analysis of these data is partitioning the sequences to identify descendants of a common B cell (i.e. a clone). Current methods group sequences using a fixed distance threshold, or a likelihood calculation that is computationally-intensive. Here, we propose a new method based on spectral clustering with an adaptive threshold to determine the local sequence neighborhood. Validation using simulated and experimental datasets demonstrates that this method has high sensitivity and specificity compared to a fixed threshold that is optimized for these measures. In addition, this method works on datasets where choosing an optimal fixed threshold is difficult and is more computationally efficient in all cases. The ability to quickly and accurately identify members of a clone from repertoire sequencing data will greatly improve downstream analyses. Clonally-related sequences cannot be treated independently in statistical models, and clonal partitions are used as the basis for the calculation of diversity metrics, lineage reconstruction and selection analysis. Thus, the spectral clustering-based method here represents an important contribution to repertoire analysis.

Availability and implementation

Source code for this method is freely available in the SCOPe (Spectral Clustering for clOne Partitioning) R package in the Immcantation framework: www.immcantation.org under the CC BY-SA 4.0 license.

Supplementary information

Supplementary data are available at Bioinformatics online.

cAb-Rep: A Database of Curated Antibody Repertoires for Exploring Antibody Diversity and Predicting Antibody Prevalence

The diversity of B cell receptors provides a basis for recognizing numerous pathogens. Antibody repertoire sequencing has revealed relationships between B cell receptor sequences, their diversity, and their function in infection, vaccination, and disease. However, many repertoire datasets have been deposited without annotation or quality control, limiting their utility. To accelerate investigations of B cell immunoglobulin sequence repertoires and to facilitate development of algorithms for their analysis, we constructed a comprehensive public database of curated human B cell immunoglobulin sequence repertoires, cAb-Rep (https://cab-rep.c2b2.columbia.edu), which currently includes 306 immunoglobulin repertoires from 121 human donors, who were healthy, vaccinated, or had autoimmune disease. The database contains a total of 267.9 million V(D)J heavy chain and 72.9 million VJ light chain transcripts. These transcripts are full-length or near full-length, have been annotated with gene origin, antibody isotype, somatic hypermutations, and other biological characteristics, and are stored in FASTA format to facilitate their direct use by most current repertoire-analysis programs. We describe a website to search cAb-Rep for similar antibodies along with methods for analysis of the prevalence of antibodies with specific genetic signatures, for estimation of reproducibility of somatic hypermutation patterns of interest, and for delineating frequencies of somatically introduced N-glycosylation. cAb-Rep should be useful for investigating attributes of B cell sequence repertoires, for understanding characteristics of affinity maturation, and for identifying potential barriers to the elicitation of effective neutralizing antibodies in infection or by vaccination.

Survival analysis of DNA mutation motifs with penalized proportional hazards

Antibodies, an essential part of our immune system, develop through an intricate process to bind a wide array of pathogens. This process involves randomly mutating DNA sequences encoding these antibodies to find variants with improved binding, though mutations are not distributed uniformly across sequence sites. Immunologists observe this nonuniformity to be consistent with "mutation motifs", which are short DNA subsequences that affect how likely a given site is to experience a mutation. Quantifying the effect of motifs on mutation rates is challenging: a large number of possible motifs makes this statistical problem high dimensional, while the unobserved history of the mutation process leads to a nontrivial missing data problem. We introduce an ℓ 1-penalized proportional hazards model to infer mutation motifs and their effects. In order to estimate model parameters, our method uses a Monte Carlo EM algorithm to marginalize over the unknown ordering of mutations. We show that our method performs better on simulated data compared to current methods and leads to more parsimonious models. The application of proportional hazards to mutation processes is, to our knowledge, novel and formalizes the current methods in a statistical framework that can be easily extended to analyze the effect of other biological features on mutation rates.

Selection and Neutral Mutations Drive Pervasive Mutability Losses in Long-Lived Anti-HIV B-Cell Lineages

High-affinity antibodies arise within weeks of infection from the evolution of B-cell receptors under selection to improve antigen recognition. This rapid adaptation is enabled by the distribution of highly mutable "hotspot" motifs in B-cell receptor genes. High mutability in antigen-binding regions (complementarity determining regions [CDRs]) creates variation in binding affinity, whereas low mutability in structurally important regions (framework regions [FRs]) may reduce the frequency of destabilizing mutations. During the response, loss of mutational hotspots and changes in their distribution across CDRs and FRs are predicted to compromise the adaptability of B-cell receptors, yet the contributions of different mechanisms to gains and losses of hotspots remain unclear. We reconstructed changes in anti-HIV B-cell receptor sequences and show that mutability losses were ∼56% more frequent than gains in both CDRs and FRs, with the higher relative mutability of CDRs maintained throughout the response. At least 21% of the total mutability loss was caused by synonymous mutations. However, nonsynonymous substitutions caused most (79%) of the mutability loss in CDRs. Because CDRs also show strong positive selection, this result suggests that selection for mutations that increase binding affinity contributed to loss of mutability in antigen-binding regions. Although recurrent adaptation to evolving viruses could indirectly select for high mutation rates, we found no evidence of indirect selection to increase or retain hotspots. Our results suggest mutability losses are intrinsic to both the neutral and adaptive evolution of B-cell populations and might constrain their adaptation to rapidly evolving pathogens such as HIV and influenza.

DNA Breaks in Ig V Regions Are Predominantly Single Stranded and Are Generated by UNG and MSH6 DNA Repair Pathways

Antibody diversity is initiated by activation-induced deaminase (AID), which deaminates cytosine to uracil in DNA. Uracils in the Ig gene loci can be recognized by uracil DNA glycosylase (UNG) or mutS homologs 2 and 6 (MSH2-MSH6) proteins, and then processed into DNA breaks. Breaks in switch regions of the H chain locus cause isotype switching and have been extensively characterized as staggered and blunt double-strand breaks. However, breaks in V regions that arise during somatic hypermutation are poorly understood. In this study, we characterize AID-dependent break formation in JH introns from mouse germinal center B cells. We used a ligation-mediated PCR assay to detect single-strand breaks and double-strand breaks that were either staggered or blunt. In contrast to switch regions, V regions contained predominantly single-strand breaks, which peaked 10 d after immunization. We then examined the pathways used to generate these breaks in UNG- and MSH6-deficient mice. Surprisingly, both DNA repair pathways contributed substantially to break formation, and in the absence of both UNG and MSH6, the frequency of breaks was severely reduced. When the breaks were sequenced and mapped, they were widely distributed over a 1000-bp intron region downstream of JH3 and JH4 exons and were unexpectedly located at all 4 nt. These data suggest that during DNA repair, nicks are generated at distal sites from the original deaminated cytosine, and these repair intermediates could generate both faithful and mutagenic repair. During mutagenesis, single-strand breaks would allow entry for low-fidelity DNA polymerases to generate somatic hypermutation.

Data-driven engineering of protein therapeutics

Protein therapeutics requires a series of properties beyond biochemical activity, including serum stability, low immunogenicity, and manufacturability. Mutations that improve one property often decrease one or more of the other essential requirements for therapeutic efficacy, making the protein engineering challenge difficult. The past decade has seen an explosion of new techniques centered around cheaply reading and writing DNA. This review highlights the recent use of such high throughput technologies for engineering protein therapeutics. Examples include the use of human antibody repertoire sequence data to pair antibody heavy and light chains, comprehensive mutational analysis for engineering antibody specificity, and the use of ancestral and inter-species sequence data to engineer simultaneous improvements in enzyme catalytic efficiency and stability. We conclude with a perspective on further ways to integrate mature protein engineering pipelines with the exponential increases in the volume of sequencing data expected in the forthcoming decade.

Cis- and trans-factors affecting AID targeting and mutagenic outcomes in antibody diversification

Antigen receptor diversification is a hallmark of adaptive immunity which allows specificity of the receptor to particular antigen. B cell receptor (BCR) or its secreted form, antibody, is diversified through antigen-independent and antigen-dependent mechanisms. During B cell development in bone marrow, BCR is diversified via V(D)J recombination mediated by RAG endonuclease. Upon stimulation by antigen, B cell undergo somatic hypermutation (SHM) to allow affinity maturation and class switch recombination (CSR) to change the effector function of the antibody. Both SHM and CSR are initiated by activation-induced cytidine deaminase (AID). Repair of AID-initiated lesions through different DNA repair pathways results in diverse mutagenic outcomes. Here, we focus on discussing cis- and trans-factors that target AID to its substrates and factors that affect different outcomes of AID-initiated lesions. The knowledge of mechanisms that govern AID targeting and outcomes could be harnessed to elicit rare functional antibodies and develop ex vivo antibody diversification approaches with diversifying base editors.

Signaling control of antibody isotype switching

Class switch recombination (CSR) generates isotype-switched antibodies with distinct effector functions essential for mediating effective humoral immunity. CSR is catalyzed by activation-induced deaminase (AID) that initiates DNA lesions in the evolutionarily conserved switch (S) regions at the immunoglobulin heavy chain (Igh) locus. AID-initiated DNA lesions are subsequently converted into DNA double stranded breaks (DSBs) in the S regions of Igh locus, repaired by non-homologous end-joining to effect CSR in mammalian B lymphocytes. While molecular mechanisms of CSR are well characterized, it remains less well understood how upstream signaling pathways regulate AID expression and CSR. B lymphocytes express multiple receptors including the B cell antigen receptor (BCR) and co-receptors (e.g., CD40). These receptors may share common signaling pathways or may use distinct signaling elements to regulate CSR. Here, we discuss how signals emanating from different receptors positively or negatively regulate AID expression and CSR.

Nomenclature of humanized mAbs: Early concepts, current challenges and future perspectives

Nomenclature of monoclonal antibodies traditionally followed a strict scheme indicating target and species information. Because of the rapid advances in this field, emphasized by approval of four humanized and six human antibodies in 2017, the International Nonproprietary Name of new antibodies was updated profoundly by removing the species substem completely. In this review we give an overview about what developments led to the preference of the scientific community towards human-like antibodies. We summarize the major updates in naming schemes that tried to classify antibodies according to their humanization technique or to the final primary sequence and how this led to the erroneous perception to indicate expected immunogenicity. Following the new 2017 nomenclature update, there will not be any information available about the species origin in the names of new antibodies, which emphasizes the need for providing additional supplemental information to the scientific community and develop tools to accurately estimate and control the safety of new monoclonal antibody molecules.

Beyond Hot Spots: Biases in Antibody Somatic Hypermutation and Implications for Vaccine Design

The evolution of antibodies in an individual during an immune response by somatic hypermutation (SHM) is essential for the ability of the immune system to recognize and remove the diverse spectrum of antigens that may be encountered. These mutations are not produced at random; nucleotide motifs that result in increased or decreased rates of mutation were first reported in 1992. Newer models that estimate the propensity for mutation for every possible 5- or 7-nucleotide motif have emphasized the complexity of SHM targeting and suggested possible new hot spot motifs. Even with these fine-grained approaches, however, non-local context matters, and the mutations observed at a specific nucleotide motif varies between species and even by locus, gene segment, and position along the gene segment within a single species. An alternative method has been provided to further abstract away the molecular mechanisms underpinning SHM, prompted by evidence that certain stereotypical amino acid substitutions are favored at each position of a particular V gene. These "substitution profiles," whether obtained from a single B cell lineage or an entire repertoire, offer a simplified approach to predict which substitutions will be well-tolerated and which will be disfavored, without the need to consider path-dependent effects from neighboring positions. However, this comes at the cost of merging the effects of two distinct biological processes, the generation of mutations, and the selection acting on those mutations. Since selection is contingent on the particular antigens an individual has been exposed to, this suggests that SHM may have evolved to prefer mutations that are most likely to be useful against pathogens that have co-evolved with us. Alternatively, the ability to select favorable mutations may be strongly limited by the biases of SHM targeting. In either scenario, the sequence space explored by SHM is significantly limited and this consequently has profound implications for the rational design of vaccine strategies.

The Bayesian optimist's guide to adaptive immune receptor repertoire analysis

Probabilistic modeling is fundamental to the statistical analysis of complex data. In addition to forming a coherent description of the data-generating process, probabilistic models enable parameter inference about given datasets. This procedure is well developed in the Bayesian perspective, in which one infers probability distributions describing to what extent various possible parameters agree with the data. In this paper, we motivate and review probabilistic modeling for adaptive immune receptor repertoire data then describe progress and prospects for future work, from germline haplotyping to adaptive immune system deployment across tissues. The relevant quantities in immune sequence analysis include not only continuous parameters such as gene use frequency but also discrete objects such as B-cell clusters and lineages. Throughout this review, we unravel the many opportunities for probabilistic modeling in adaptive immune receptor analysis, including settings for which the Bayesian approach holds substantial promise (especially if one is optimistic about new computational methods). From our perspective, the greatest prospects for progress in probabilistic modeling for repertoires concern ancestral sequence estimation for B-cell receptor lineages, including uncertainty from germline genotype, rearrangement, and lineage development.

Functional Relevance of Improbable Antibody Mutations for HIV Broadly Neutralizing Antibody Development

HIV-1 broadly neutralizing antibodies (bnAbs) require high levels of activation-induced cytidine deaminase (AID)-catalyzed somatic mutations for optimal neutralization potency. Probable mutations occur at sites of frequent AID activity, while improbable mutations occur where AID activity is infrequent. One bottleneck for induction of bnAbs is the evolution of viral envelopes (Envs) that can select bnAb B cell receptors (BCR) with improbable mutations. Here we define the probability of bnAb mutations and demonstrate the functional significance of key improbable mutations in three bnAb B cell lineages. We show that bnAbs are enriched for improbable mutations, which implies that their elicitation will be critical for successful vaccine induction of potent bnAb B cell lineages. We discuss a mutation-guided vaccine strategy for identification of Envs that can select B cells with BCRs that have key improbable mutations required for bnAb development.

•

HIV-1 broadly neutralizing antibodies are enriched with low-probability mutations

•

Improbable mutations can be functionally critical for bnAb neutralization breadth

•

Critical improbable mutations are high-value targets for selection with vaccines

Analyzing Immunoglobulin Repertoires

Somatic assembly of T cell receptor and B cell receptor (BCR) genes produces a vast diversity of lymphocyte antigen recognition capacity. The advent of efficient high-throughput sequencing of lymphocyte antigen receptor genes has recently generated unprecedented opportunities for exploration of adaptive immune responses. With these opportunities have come significant challenges in understanding the analysis techniques that most accurately reflect underlying biological phenomena. In this regard, sample preparation and sequence analysis techniques, which have largely been borrowed and adapted from other fields, continue to evolve. Here, we review current methods and challenges of library preparation, sequencing and statistical analysis of lymphocyte receptor repertoire studies. We discuss the general steps in the process of immune repertoire generation including sample preparation, platforms available for sequencing, processing of sequencing data, measurable features of the immune repertoire, and the statistical tools that can be used for analysis and interpretation of the data. Because BCR analysis harbors additional complexities, such as immunoglobulin (Ig) (i.e., antibody) gene somatic hypermutation and class switch recombination, the emphasis of this review is on Ig/BCR sequence analysis.