Impact of new sequencing technologies on studies of the human B cell repertoire

AntiRef: reference clusters of human antibody sequences

Motivation

Genetic biases in the human antibody repertoire result in publicly available antibody sequence datasets that contain many duplicate or highly similar sequences. Available datasets are further skewed by the predominance of studies focused on specific disease states, primarily cancer, autoimmunity, and a small number of infectious diseases that includes HIV, influenza, and SARS-CoV-2. These biases and redundancies are a barrier to rapid similarity searches and reduce the efficiency with which these datasets can be used to train statistical or machine-learning models. Identity-based clustering provides a solution; however, the extremely large size of available antibody sequence datasets makes such clustering operations computationally intensive and potentially out of reach for many scientists and researchers who would benefit from such data.

Results

Antibody Reference Clusters (AntiRef), which is modeled after UniRef, provides clustered datasets of filtered human antibody sequences. Due to the modular nature of recombined antibody genes, the clustering thresholds used by UniRef for general protein sequences are suboptimal for antibody clustering. Starting with an input dataset of ∼451M full-length, productive human antibody sequences, AntiRef provides reference datasets clustered at a range of antibody-optimized identity thresholds. AntiRef90 is one-third the size of the input dataset and less than half the size of the non-redundant AntiRef100.

Availability and implementation

AntiRef datasets are available on Zenodo (zenodo.org/record/7474336). All code used to generate AntiRef is available on GitHub (github.com/briney/antiref). The AntiRef versioning scheme (current version: v2022.12.14) refers to the date on which sequences were retrieved from OAS.

Normal IgH Repertoire Diversity in an Infant with ADA Deficiency After Gene Therapy

PURPOSE

Adenosine deaminase (ADA) deficiency causes severe combined immunodeficiency (SCID) through an accumulation of toxic metabolites within lymphocytes. Recently, ADA deficiency has been successfully treated using lentiviral-transduced autologous CD34+ cells carrying the ADA gene. T and B cell function appears to be fully restored, but in many patients' B cell numbers remain low, and assessments of the immunoglobulin heavy (IgHV) repertoire following gene therapy are lacking.

METHODS

We performed deep sequencing of IgHV repertoire in peripheral blood lymphocytes from a child following lentivirus-based gene therapy for ADA deficiency and compared to the IgHV repertoire in healthy infants and adults.

RESULTS

After gene therapy, Ig diversity increased over time as evidenced by V, D, and J gene usage, N-additions, CDR3 length, extent of somatic hypermutation, and Ig class switching. There was the emergence of predominant IgHM, IgHG, and IgHA CDR3 lengths after gene therapy indicating successful oligoclonal expansion in response to antigens. This provides proof of concept for the feasibility and utility of molecular monitoring in following B cell reconstitution following gene therapy for ADA deficiency.

CONCLUSION

Based on deep sequencing, gene therapy resulted in an IgHV repertoire with molecular diversity similar to healthy infants.

Immune repertoire fingerprinting by principal component analysis reveals shared features in subject groups with common exposures

BACKGROUND

Advances in next-generation sequencing (NGS) of antibody repertoires have led to an explosion in B cell receptor sequence data from donors with many different disease states. These data have the potential to detect patterns of immune response across populations. However, to this point it has been difficult to interpret such patterns of immune response between disease states in the absence of functional data. There is a need for a robust method that can be used to distinguish general patterns of immune responses at the antibody repertoire level.

RESULTS

We developed a method for reducing the complexity of antibody repertoire datasets using principal component analysis (PCA) and refer to our method as "repertoire fingerprinting." We reduce the high dimensional space of an antibody repertoire to just two principal components that explain the majority of variation in those repertoires. We show that repertoires from individuals with a common experience or disease state can be clustered by their repertoire fingerprints to identify common antibody responses.

CONCLUSIONS

Our repertoire fingerprinting method for distinguishing immune repertoires has implications for characterizing an individual disease state. Methods to distinguish disease states based on pattern recognition in the adaptive immune response could be used to develop biomarkers with diagnostic or prognostic utility in patient care. Extending our analysis to larger cohorts of patients in the future should permit us to define more precisely those characteristics of the immune response that result from natural infection or autoimmunity.

Structural Basis for the Broad, Antibody-Mediated Neutralization of H5N1 Influenza Virus

Human infection with highly pathogenic avian influenza A viruses causes severe disease and fatalities. We previously identified a potent and broadly neutralizing antibody (bnAb), 13D4, against the H5N1 virus. Here, we report the co-crystal structure of 13D4 in complex with the hemagglutinin (HA) of A/Vietnam/1194/2004 (H5N1). We show that heavy-chain complementarity-determining region 3 (HCDR3) of 13D4 confers broad yet specific neutralization against H5N1, undergoing conformational rearrangement to bind to the receptor binding site (RBS). Further, we show that mutating four critical residues within the RBS-Trp153, Lys156, Lys193, and Leu194-disrupts the binding between 13D4 and HA. Viruses bearing Asn193 instead of Lys/Arg can evade 13D4 neutralization, indicating that Lys193 polymorphism might be, at least in part, involved in the antigenicity of recent H5 genotypes (such as H5N6 and H5N8) as distinguished from H5N1. BnAb 13D4 may offers a template for therapeutic RBS inhibitor design and serve as an indicator of antigenic change for current H5 viruses.IMPORTANCE Infection by highly pathogenic avian influenza A virus remains a threat to public health. Our broadly neutralizing antibody, 13D4, is capable of neutralizing all representative H5N1 viruses and protecting mice against lethal challenge. Structural analysis revealed that 13D4 uses heavy-chain complementarity-determining region 3 (HCDR3) to fit the receptor binding site (RBS) via conformational rearrangement. Four conserved residues within the RBS are critical for the broad potency of 13D4. Importantly, polymorphism of Lys193 on the RBS may be associated with the antigenicity shift from H5N1 to other newly emerging viruses, such as H5N6 and H5N8. Our findings may pave the way for highly pathogenic avian influenza virus vaccine development and therapeutic RBS inhibitor design.

Biosensor-based epitope mapping of antibodies targeting the hemagglutinin and neuraminidase of influenza A virus

Characterization of the epitopes on antigen recognized by monoclonal antibodies (mAb) is useful for the development of therapeutic antibodies, diagnostic tools, and vaccines. Epitope mapping also provides functional information for sequence-based repertoire analysis of antibody response to pathogen infection and/or vaccination. However, development of mapping strategies has lagged behind mAb discovery. We have developed a site-directed mutagenesis approach that can be used in conjunction with bio-layer interferometry (BLI) biosensors to map mAb epitopes. By generating a panel of single point mutants in the recombinant hemagglutinin (HA) and neuraminidase (NA) proteins of influenza A viruses, we have characterized the epitopes of hundreds of mAbs targeting the H1 and H3 subtypes of HA and the N9 subtype of NA.

Quantitative Analysis of Repertoire-Scale Immunoglobulin Properties in Vaccine-Induced B-Cell Responses

Recent advances in the next-generation sequencing of B-cell receptors (BCRs) enable the characterization of humoral responses at a repertoire-wide scale and provide the capability for identifying unique features of immune repertoires in response to disease, vaccination, or infection. Immunosequencing now readily generates 10³–10⁵ sequences per sample; however, statistical analysis of these repertoires is challenging because of the high genetic diversity of BCRs and the elaborate clonal relationships among them. To date, most immunosequencing analyses have focused on reporting qualitative trends in immunoglobulin (Ig) properties, such as usage or somatic hypermutation (SHM) percentage of the Ig heavy chain variable (IGHV) gene segment family, and on reducing complex Ig property distributions to simple summary statistics. However, because Ig properties are typically not normally distributed, any approach that fails to assess the distribution as a whole may be inadequate in (1) properly assessing the statistical significance of repertoire differences, (2) identifying how two repertoires differ, and (3) determining appropriate confidence intervals for assessing the size of the differences and their potential biological relevance. To address these issues, we have developed a technique that uses Wilcox’ robust statistics toolbox to identify statistically significant vaccine-specific differences between Ig repertoire properties. The advantage of this technique is that it can determine not only whether but also where the distributions differ, even when the Ig repertoire properties are non-normally distributed. We used this technique to characterize murine germinal center (GC) B-cell repertoires in response to a complex Ebola virus-like particle (eVLP) vaccine candidate with known protective efficacy. The eVLP-mediated GC B-cell responses were highly diverse, consisting of thousands of clonotypes. Despite this staggering diversity, we identified statistically significant differences between non-immunized, vaccine only, and vaccine-plus-adjuvant groups in terms of Ig properties, including IGHV-family usage, SHM percentage, and characteristics of the BCR complementarity-determining region. Most notably, our analyses identified a robust eVLP-specific feature—enhanced IGHV8-family usage in B-cell repertoires. These findings demonstrate the utility of our technique in identifying statistically significant BCR repertoire differences following vaccination. More generally, our approach is potentially applicable to a wide range of studies in infection, vaccination, auto-immunity, and cancer.

ADAMTS13 autoantibodies cloned from patients with acquired thrombotic thrombocytopenic purpura: 1. Structural and functional characterization in vitro

BACKGROUND

Acquired thrombotic thrombocytopenia purpura (TTP) is a life-threatening illness caused by autoantibodies that decrease the activity of ADAMTS13, the von Willebrand factor-cleaving protease. Despite efficacy of plasma exchange, mortality remains high and relapse is common. Improved therapies may come from understanding the diversity of pathogenic autoantibodies on a molecular or genetic level. Cloning comprehensive repertoires of patient autoantibodies can provide the necessary tools for studying immunobiology of disease and developing animal models.

STUDY DESIGN AND METHODS

Anti-ADAMTS13 antibodies were cloned from four patients with acquired TTP using phage display and characterized with respect to genetic origin, inhibition of ADAMTS13 proteolytic activity, and epitope specificity. Anti-idiotypic antisera raised to a subset of autoantibodies enabled comparison of their relatedness to each other and to polyclonal immunoglobulin (Ig)G in patient plasma.

RESULTS

Fifty-one unique antibodies were isolated comprising epitope specificities resembling the diversity found in circulating patient IgG. Antibodies directed both to the amino terminal domains and to those requiring the ADAMTS13 cysteine-rich/spacer region for binding inhibited proteolytic activity, while those solely targeting carboxy-terminal domains were noninhibitory. Anti-idiotypic antisera raised to a subset of antibody clones crossreacted with and reduced the inhibitory activity of polyclonal IgG from a set of unrelated patients.

CONCLUSIONS

Anti-ADAMTS13 autoantibodies isolated by repertoire cloning display the diversity of epitope specificities found in patient plasma and provide tools for developing animal models of acquired TTP. Shared idiotypes of inhibitory clones with circulating IgG from multiple patients suggest common features of pathogenic autoantibodies that could be exploited for developing more targeted therapies.

Improving Loop Modeling of the Antibody Complementarity-Determining Region 3 Using Knowledge-Based Restraints

Structural restrictions are present even in the most sequence diverse portions of antibodies, the complementary determining region (CDR) loops. Previous studies identified robust rules that define canonical structures for five of the six CDR loops, however the heavy chain CDR 3 (HCDR3) defies standard classification attempts. The HCDR3 loop can be subdivided into two domains referred to as the “torso” and the “head” domains and two major families of canonical torso structures have been identified; the more prevalent “bulged” and less frequent “non-bulged” torsos. In the present study, we found that Rosetta loop modeling of 28 benchmark bulged HCDR3 loops is improved with knowledge-based structural restraints developed from available antibody crystal structures in the PDB. These restraints restrict the sampling space Rosetta searches in the torso domain, limiting the φ and ψ angles of these residues to conformations that have been experimentally observed. The application of these restraints in Rosetta result in more native-like structure sampling and improved score-based differentiation of native-like HCDR3 models, significantly improving our ability to model antibody HCDR3 loops.

Clonify: unseeded antibody lineage assignment from next-generation sequencing data

Defining the dynamics and maturation processes of antibody clonal lineages is crucial to understanding the humoral response to infection and immunization. Although individual antibody lineages have been previously analyzed in isolation, these studies provide only a narrow view of the total antibody response. Comprehensive study of antibody lineages has been limited by the lack of an accurate clonal lineage assignment algorithm capable of operating on next-generation sequencing datasets. To address this shortcoming, we developed Clonify, which is able to perform unseeded lineage assignment on very large sets of antibody sequences. Application of Clonify to IgG+ memory repertoires from healthy individuals revealed a surprising lack of influence of large extended lineages on the overall repertoire composition, indicating that this composition is driven less by the order and frequency of pathogen encounters than previously thought. Clonify is freely available at www.github.com/briney/clonify-python.

Long antibody HCDR3s from HIV-naïve donors presented on a PG9 neutralizing antibody background mediate HIV neutralization

Development of broadly neutralizing antibodies (bnAbs) against HIV-1 usually requires prolonged infection and induction of Abs with unusual features, such as long heavy-chain complementarity-determining region 3 (HCDR3) loops. Here we sought to determine whether the repertoires of HIV-1-naïve individuals contain Abs with long HCDR3 loops that could mediate HIV-1 neutralization. We interrogated at massive scale the structural properties of long Ab HCDR3 loops in HIV-1-naïve donors, searching for structured HCDR3s similar to those of the HIV-1 bnAb PG9. We determined the nucleotide sequences encoding 2.3 × 10(7)unique HCDR3 amino acid regions from 70 different HIV-1-naïve donors. Of the 26,917 HCDR3 loops with 30-amino acid length identified, we tested 30 for further study that were predicted to have PG9-like structure when chimerized onto PG9. Three of these 30 PG9 chimeras bound to the HIV-1 gp120 monomer, and two were neutralizing. In addition, we found 14 naturally occurring HCDR3 sequences that acquired the ability to bind to the HIV-1 gp120 monomer when adding 2- to 7-amino acid mutations via computational design. Of those 14 designed Abs, 8 neutralized HIV-1, with IC50values ranging from 0.7 to 98 µg/mL. These data suggest that the repertoire of HIV-1-naïve individuals contains rare B cells that encode HCDR3 loops that bind or neutralize HIV-1 when presented on a PG9 background with relatively few or no additional mutations. Long HCDR3 sequences are present in the HIV-naïve B-cell repertoire, suggesting that this class of bnAbs is a favorable target for rationally designed preventative vaccine efforts.

Reproducibility of Illumina platform deep sequencing errors allows accurate determination of DNA barcodes in cells

Background

Next generation sequencing (NGS) of amplified DNA is a powerful tool to describe genetic heterogeneity within cell populations that can both be used to investigate the clonal structure of cell populations and to perform genetic lineage tracing. For applications in which both abundant and rare sequences are biologically relevant, the relatively high error rate of NGS techniques complicates data analysis, as it is difficult to distinguish rare true sequences from spurious sequences that are generated by PCR or sequencing errors. This issue, for instance, applies to cellular barcoding strategies that aim to follow the amount and type of offspring of single cells, by supplying these with unique heritable DNA tags.

Results

Here, we use genetic barcoding data from the Illumina HiSeq platform to show that straightforward read threshold-based filtering of data is typically insufficient to filter out spurious barcodes. Importantly, we demonstrate that specific sequencing errors occur at an approximately constant rate across different samples that are sequenced in parallel. We exploit this observation by developing a novel approach to filter out spurious sequences.

Conclusions

Application of our new method demonstrates its value in the identification of true sequences amongst spurious sequences in biological data sets.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-0999-4) contains supplementary material, which is available to authorized users.

Deciphering the human immunome

Technological advances in next-generation DNA sequencing offer great potential to probe the human immune system. On 3 April 2015, leading immunologists and bioinformatics scientists met to consider how best to harness these advances for decoding the human immunome.

A Simple Flow-Cytometric Method Measuring B Cell Surface Immunoglobulin Avidity Enables Characterization of Affinity Maturation to Influenza A Virus

Antibody (Ab) affinity maturation enables an individual to maintain immunity to an increasing number of pathogens within the limits of a total Ig production threshold. A better understanding of this process is critical for designing vaccines that generate optimal Ab responses to pathogens. Our study describes a simple flow-cytometric method that enumerates virus-specific germinal center (GC) B cells as well as their AC₅₀, a measure of Ab avidity, defined as the antigen concentration required to detect 50% of specific B cells. Using a model of mouse Ab responses to the influenza A virus hemagglutinin (IAV HA), we obtained data indicating that AC₅₀ decreases with time postinfection in an affinity maturation-dependent process. As proof of principle of the utility of the method, our data clearly show that relative to intranasal IAV infection, intramuscular immunization against inactivated IAV in adjuvant results in a diminished GC HA B cell response, with increased AC₅₀ correlating with an increased serum Ab off-rate. Enabling simultaneous interrogation of both GC HA B cell quantity and quality, this technique should facilitate study of affinity maturation and rational vaccine design.

Though it was first described 50 years ago, little is known about how antibody affinity maturation contributes to immunity. This question is particularly relevant to developing more effective vaccines for influenza A virus (IAV) and other viruses that are difficult vaccine targets. Limitations in methods for characterizing antigen-specific B cells have impeded progress in characterizing the quality of immune responses to vaccine and natural immunogens. In this work, we describe a simple flow cytometry-based approach that measures both the number and affinity of IAV-binding germinal center B cells specific for the IAV HA, the major target of IAV-neutralizing antibodies. Using this method, we showed that the route and form of immunization significantly impacts the quality and quantity of B cell antibody responses. This method provides a relatively simple yet powerful tool for better understanding the contribution of affinity maturation to viral immunity.

Association of VH4-59 Antibody Variable Gene Usage with Recognition of an Immunodominant Epitope on the HIV-1 Gag Protein

The human antibody response against HIV-1 infection recognizes diverse antigenic subunits of the virion, and includes a high level of antibodies to the Gag protein. We report here the isolation and characterization of a subset of Gag-specific human monoclonal antibodies (mAbs) that were prevalent in the antibody repertoire of an HIV-infected individual. Several lineages of Gag-specifc mAbs were encoded by a single antibody heavy chain variable region, V_H4-59, and a representative antibody from this group designated mAb 3E4 recognized a linear epitope on the globular head of the p17 subunit of Gag. We found no evidence that mAb 3E4 exhibited any function in laboratory studies aimed at elucidating the immunologic activity, including assays for neutralization, Ab-dependent cell-mediated virus inhibition, or enhanced T cell reactivity caused by Gag-3E4 complexes. The findings suggest this immunodominant epitope in Gag protein, which is associated with V_H4-59 germline gene usage, may induce a high level of B cells that encode binding but non-functional antibodies that occupy significant repertoire space following HIV infection. The studies define an additional specific molecular mechanism in the immune distraction activity of the HIV virion.

Predominant structural configuration of natural antibody repertoires enables potent antibody responses against protein antigens

Humoral immunity against diverse pathogens is rapidly elicited from natural antibody repertoires of limited complexity. But the organizing principles underlying the antibody repertoires that facilitate this immunity are not well-understood. We used HER2 as a model immunogen and reverse-engineered murine antibody response through constructing an artificial antibody library encoded with rudimentary sequence and structural characteristics learned from high throughput sequencing of antibody variable domains. Antibodies selected in vitro from the phage-displayed synthetic antibody library bound to the model immunogen with high affinity and specificities, which reproduced the specificities of natural antibody responses. We conclude that natural antibody structural repertoires are shaped to allow functional antibodies to be encoded efficiently, within the complexity limit of an individual antibody repertoire, to bind to diverse protein antigens with high specificity and affinity. Phage-displayed synthetic antibody libraries, in conjunction with high-throughput sequencing, can thus be designed to replicate natural antibody responses and to generate novel antibodies against diverse antigens.

Sequencing the functional antibody repertoire--diagnostic and therapeutic discovery

The development of high-throughput DNA sequencing technologies has enabled large-scale characterization of functional antibody repertoires, a new method of understanding protective and pathogenic immune responses. Important parameters to consider when sequencing antibody repertoires include the methodology, the B-cell population and clinical characteristics of the individuals analysed, and the bioinformatic analysis. Although focused sequencing of immunoglobulin heavy chains or complement determining regions can be utilized to monitor particular immune responses and B-cell malignancies, high-fidelity analysis of the full-length paired heavy and light chains expressed by individual B cells is critical for characterizing functional antibody repertoires. Bioinformatic identification of clonal antibody families and recombinant expression of representative members produces recombinant antibodies that can be used to identify the antigen targets of functional immune responses and to investigate the mechanisms of their protective or pathogenic functions. Integrated analysis of coexpressed functional genes provides the potential to further pinpoint the most important antibodies and clonal families generated during an immune response. Sequencing antibody repertoires is transforming our understanding of immune responses to autoimmunity, vaccination, infection and cancer. We anticipate that antibody repertoire sequencing will provide next-generation biomarkers, diagnostic tools and therapeutic antibodies for a spectrum of diseases, including rheumatic diseases.

Immunology of bovine respiratory syncytial virus in calves

Bovine respiratory syncytial virus (BRSV) is an important cause of respiratory disease in young calves. The virus is genetically and antigenically closely related to human (H)RSV, which is a major cause of respiratory disease in young infants. As a natural pathogen of calves, BRSV infection recapitulates the pathogenesis of respiratory disease in man more faithfully than semi-permissive, animal models of HRSV infection. With the increasing availability of immunological reagents, the calf can be used to dissect the pathogenesis of and mechanisms of immunity to RSV infection, to analyse the ways in which the virus proteins interact with components of the innate response, and to evaluate RSV vaccine strategies. Passively transferred, neutralising bovine monoclonal antibodies, which recognise the same epitopes in the HRSV and BRSV fusion (F) protein, can protect calves against BRSV infection, and depletion of different T cells subsets in calves has highlighted the importance of CD8(+) T cells in viral clearance. Calves can be used to model maternal-antibody mediated suppression of RSV vaccine efficacy, and to increase understanding of the mechanisms responsible for RSV vaccine-enhanced respiratory disease.

In-depth determination and analysis of the human paired heavy- and light-chain antibody repertoire

High-throughput immune repertoire sequencing has emerged as a critical step in the understanding of adaptive responses following infection or vaccination or in autoimmunity. However, determination of native antibody variable heavy-light pairs (VH-VL pairs) remains a major challenge, and no technologies exist to adequately interrogate the >1 × 10(6) B cells in typical specimens. We developed a low-cost, single-cell, emulsion-based technology for sequencing antibody VH-VL repertoires from >2 × 10(6) B cells per experiment with demonstrated pairing precision >97%. A simple flow-focusing apparatus was used to sequester single B cells into emulsion droplets containing lysis buffer and magnetic beads for mRNA capture; subsequent emulsion RT-PCR generated VH-VL amplicons for next-generation sequencing. Massive VH-VL repertoire analyses of three human donors provided new immunological insights including (i) the identity, frequency and pairing propensity of shared, or 'public', VL genes, (ii) the detection of allelic inclusion (an implicated autoimmune mechanism) in healthy individuals and (iii) the occurrence of antibodies with features, in terms of gene usage and CDR3 length, associated with broadly neutralizing antibodies to rapidly evolving viruses such as HIV-1 and influenza.

A high-throughput shotgun mutagenesis approach to mapping B-cell antibody epitopes

Characterizing the binding sites of monoclonal antibodies (mAbs) on protein targets, their 'epitopes', can aid in the discovery and development of new therapeutics, diagnostics and vaccines. However, the speed of epitope mapping techniques has not kept pace with the increasingly large numbers of mAbs being isolated. Obtaining detailed epitope maps for functionally relevant antibodies can be challenging, particularly for conformational epitopes on structurally complex proteins. To enable rapid epitope mapping, we developed a high-throughput strategy, shotgun mutagenesis, that enables the identification of both linear and conformational epitopes in a fraction of the time required by conventional approaches. Shotgun mutagenesis epitope mapping is based on large-scale mutagenesis and rapid cellular testing of natively folded proteins. Hundreds of mutant plasmids are individually cloned, arrayed in 384-well microplates, expressed within human cells, and tested for mAb reactivity. Residues are identified as a component of a mAb epitope if their mutation (e.g. to alanine) does not support candidate mAb binding but does support that of other conformational mAbs or allows full protein function. Shotgun mutagenesis is particularly suited for studying structurally complex proteins because targets are expressed in their native form directly within human cells. Shotgun mutagenesis has been used to delineate hundreds of epitopes on a variety of proteins, including G protein-coupled receptor and viral envelope proteins. The epitopes mapped on dengue virus prM/E represent one of the largest collections of epitope information for any viral protein, and results are being used to design better vaccines and drugs.

Molecular signatures of hemagglutinin stem-directed heterosubtypic human neutralizing antibodies against influenza A viruses

Recent studies have shown high usage of the IGHV1-69 germline immunoglobulin gene for influenza hemagglutinin stem-directed broadly-neutralizing antibodies (HV1-69-sBnAbs). Here we show that a major structural solution for these HV1-69-sBnAbs is achieved through a critical triad comprising two CDR-H2 loop anchor residues (a hydrophobic residue at position 53 (Ile or Met) and Phe54), and CDR-H3-Tyr at positions 98±1; together with distinctive V-segment CDR amino acid substitutions that occur in positions sparse in AID/polymerase-η recognition motifs. A semi-synthetic IGHV1-69 phage-display library screen designed to investigate AID/polη restrictions resulted in the isolation of HV1-69-sBnAbs that featured a distinctive Ile52Ser mutation in the CDR-H2 loop, a universal CDR-H3 Tyr at position 98 or 99, and required as little as two additional substitutions for heterosubtypic neutralizing activity. The functional importance of the Ile52Ser mutation was confirmed by mutagenesis and by BCR studies. Structural modeling suggests that substitution of a small amino acid at position 52 (or 52a) facilitates the insertion of CDR-H2 Phe54 and CDR-H3-Tyr into adjacent pockets on the stem. These results support the concept that activation and expansion of a defined subset of IGHV1-69-encoded B cells to produce potent HV1-69-sBnAbs does not necessarily require a heavily diversified V-segment acquired through recycling/reentry into the germinal center; rather, the incorporation of distinctive amino acid substitutions by Phase 2 long-patch error-prone repair of AID-induced mutations or by random non-AID SHM events may be sufficient. We propose that these routes of B cell maturation should be further investigated and exploited as a pathway for HV1-69-sBnAb elicitation by vaccination.

The quest for universal influenza vaccine has gained wide interest with the discovery of human neutralizing antibodies that are able to variably cross neutralize and protect against different influenza strains, subtypes, groups and lineages. These antibodies, which bind to a highly conserved epitope in the hemagglutinin stem, are often encoded by rearranged IGHV1-69 germline genes that alone make contact with HA and prevent virus entry and emergence of escape mutants. Our study was undertaken to gain an understanding of what structural requirements enable a rearranged IGHV1-69 Ab to become a potent cross-neutralizing antibody. We found that in addition to a critical amino acid triad consisting of a pair of anchor residues in CDR-H2 and a properly positioned CDR-H3 Tyr, distinctive V-segment substitutions that arise in positions that are distinct from phase I AID somatic hypermutation (SHM) hotspot motifs are often required. As few as two V-segment SHM can fulfill this role which appears to facilitate the optimal binding of CDR-H2 Phe54 and CHR-H3-Tyr into adjacent hydrophobic pockets in the HA stem. These studies provide new information on the SHM requirements for IGHV1-69-encoded B cells to produce HV1-69-sBnAbs and suggest that there may exist alternative routes to their elicitation by vaccination.