Reconstructing a B-cell clonal lineage. I. Statistical inference of unobserved ancestors

Antigenic sin and multiple breakthrough infections drive converging evolution of COVID-19 neutralizing responses

Understanding the evolution of the B cell response to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants is fundamental to design the next generation of vaccines and therapeutics. We longitudinally analyze at the single-cell level almost 900 neutralizing human monoclonal antibodies (nAbs) isolated from vaccinated people and from individuals with hybrid and super hybrid immunity (SH), developed after three mRNA vaccine doses and two breakthrough infections. The most potent neutralization and Fc functions against highly mutated variants belong to the SH cohort. Repertoire analysis shows that the original Wuhan antigenic sin drives the convergent expansion of the same B cell germlines in vaccinated and SH cohorts. Only Omicron breakthrough infections expand previously unseen germ lines and generate broadly nAbs by restoring IGHV3-53/3-66 germ lines. Our analyses find that B cells initially expanded by the original antigenic sin continue to play a fundamental role in the evolution of the immune response toward an evolving virus.

Repeated Omicron exposures redirect SARS-CoV-2-specific memory B cell evolution toward the latest variants

Immunological imprinting by ancestral SARS-CoV-2 strains is thought to impede the robust induction of Omicron-specific humoral responses by Omicron-based booster vaccines. Here, we analyzed the specificity and neutralization activity of memory B (Bmem) cells after repeated BA.5 exposure in individuals previously imprinted by ancestral strain-based mRNA vaccines. After a second BA.5 exposure, Bmem cells with BA.5 spike protein-skewed reactivity were promptly elicited, correlating with preexisting antibody titers. Clonal lineage analysis identified BA.5-skewed Bmem cells that had redirected their specificity from the ancestral strain to BA.5 through somatic hypermutations. Moreover, Bmem cells with redirected BA.5 specificity exhibited accelerated development compared with de novo Bmem cells derived from naïve repertoires. This redirected BA.5 specificity demonstrated greater resilience to viral point mutation and adaptation to recent Omicron variants HK.3 and JN.1, months after the second BA.5 exposure, suggesting that existing Bmem cells elicited by older vaccines can redirect their specificity toward newly evolving variants.

Vaccine induction of heterologous HIV-1-neutralizing antibody B cell lineages in humans

A critical roadblock to HIV vaccine development is the inability to induce B cell lineages of broadly neutralizing antibodies (bnAbs) in humans. In people living with HIV-1, bnAbs take years to develop. The HVTN 133 clinical trial studied a peptide/liposome immunogen targeting B cell lineages of HIV-1 envelope (Env) membrane-proximal external region (MPER) bnAbs (NCT03934541). Here, we report MPER peptide-liposome induction of polyclonal HIV-1 B cell lineages of mature bnAbs and their precursors, the most potent of which neutralized 15% of global tier 2 HIV-1 strains and 35% of clade B strains with lineage initiation after the second immunization. Neutralization was enhanced by vaccine selection of improbable mutations that increased antibody binding to gp41 and lipids. This study demonstrates proof of concept for rapid vaccine induction of human B cell lineages with heterologous neutralizing activity and selection of antibody improbable mutations and outlines a path for successful HIV-1 vaccine development.

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.

Characterizing adjuvants' effects at murine immunoglobulin repertoire level

Generating large-scale, high-fidelity sequencing data is challenging and, furthermore, not much has been done to characterize adjuvants' effects at the repertoire level. Thus, we introduced an IgSeq pipeline that standardized library prep protocols and data analysis functions for accurate repertoire profiling. We then studied systemically effects of CpG and Alum on the Ig heavy chain repertoire using the ovalbumin (OVA) murine model. Ig repertoires of different tissues (spleen and bone marrow) and isotypes (IgG and IgM) were examined and compared in IGHV mutation, gene usage, CDR3 length, clonal diversity, and clonal selection. We found Ig repertoires of different compartments exhibited distinguishable profiles at the non-immunized steady state, and distinctions became more pronounced upon adjuvanted immunizations. Notably, Alum and CpG effects exhibited different tissue- and isotype-preferences. The former led to increased diversity of abundant clones in bone marrow, and the latter promoted the selection of IgG clones in both tissues.

Vaccine induction of CD4-mimicking HIV-1 broadly neutralizing antibody precursors in macaques

The CD4-binding site (CD4bs) is a conserved epitope on HIV-1 envelope (Env) that can be targeted by protective broadly neutralizing antibodies (bnAbs). HIV-1 vaccines have not elicited CD4bs bnAbs for many reasons, including the occlusion of CD4bs by glycans, expansion of appropriate naive B cells with immunogens, and selection of functional antibody mutations. Here, we demonstrate that immunization of macaques with a CD4bs-targeting immunogen elicits neutralizing bnAb precursors with structural and genetic features of CD4-mimicking bnAbs. Structures of the CD4bs nAb bound to HIV-1 Env demonstrated binding angles and heavy-chain interactions characteristic of all known human CD4-mimicking bnAbs. Macaque nAb were derived from variable and joining gene segments orthologous to the genes of human VH1-46-class bnAb. This vaccine study initiated in primates the B cells from which CD4bs bnAbs can derive, accomplishing the key first step in the development of an effective HIV-1 vaccine.

B cell phylogenetics in the single cell era

The widespread availability of single-cell RNA sequencing (scRNA-seq) has led to the development of new methods for understanding immune responses. Single-cell transcriptome data can now be paired with B cell receptor (BCR) sequences. However, RNA from BCRs cannot be analyzed like most other genes because BCRs are genetically diverse within individuals. In humans, BCRs are shaped through recombination followed by mutation and selection for antigen binding. As these processes co-occur with cell division, B cells can be studied using phylogenetic trees representing the mutations within a clone. B cell trees can link experimental timepoints, tissues, or cellular subtypes. Here, we review the current state and potential of how B cell phylogenetics can be combined with single-cell data to understand immune responses.

Engineering immunogens that select for specific mutations in HIV broadly neutralizing antibodies

Vaccine development targeting rapidly evolving pathogens such as HIV-1 requires induction of broadly neutralizing antibodies (bnAbs) with conserved paratopes and mutations, and, in some cases, the same Ig-heavy chains. The current trial-and-error search for immunogen modifications that improve selection for specific bnAb mutations is imprecise. To precisely engineer bnAb boosting immunogens, we used molecular dynamics simulations to examine encounter states that form when antibodies collide with the HIV-1 Envelope (Env). By mapping how bnAbs use encounter states to find their bound states, we identified Env mutations that were predicted to select for specific antibody mutations in two HIV-1 bnAb B cell lineages. The Env mutations encoded antibody affinity gains and selected for desired antibody mutations in vivo. These results demonstrate proof-of-concept that Env immunogens can be designed to directly select for specific antibody mutations at residue-level precision by vaccination, thus demonstrating the feasibility of sequential bnAb-inducing HIV-1 vaccine design.

A new class of antibodies that overcomes a steric barrier to cross-group neutralization of influenza viruses

Antibody titers that inhibit the influenza virus hemagglutinin (HA) from engaging its receptor are the accepted correlate of protection from infection. Many potent antibodies with broad, intra-subtype specificity bind HA at the receptor binding site (RBS). One barrier to broad H1-H3 cross-subtype neutralization is an insertion (133a) between positions 133 and 134 on the rim of the H1 HA RBS. We describe here a class of antibodies that overcomes this barrier. These genetically unrestricted antibodies are abundant in the human B cell memory compartment. Analysis of the affinities of selected members of this class for historical H1 and H3 isolates suggest that they were elicited by H3 exposure and broadened or diverted by later exposure(s) to H1 HA. RBS mutations in egg-adapted vaccine strains cause the new H1 specificity of these antibodies to depend on the egg adaptation. The results suggest that suitable immunogens might elicit 133a-independent, H1-H3 cross neutralization by RBS-directed antibodies.

Computing the inducibility of B cell lineages under a context-dependent model of affinity maturation: Applications to sequential vaccine design

A key challenge in B cell lineage-based vaccine design is understanding the inducibility of target neutralizing antibodies. We approach this problem through the use of detailed stochastic modeling of the somatic hypermutation process that occurs during affinity maturation. Under such a model, sequence mutation rates are context-dependent, rendering standard probability calculations for sequence evolution intractable. We develop an algorithmic approach to rapid, accurate approximation of key marginal sequence likelihoods required to inform modern sequential vaccine design strategies. These calculated probabilities are used to define an inducibility index for selecting among potential targets for immunogen design. We apply this approach to the problem of choosing targets for the design of boosting immunogens aimed at elicitation of the HIV broadly-neutralizing antibody DH270min11.

Engaging an HIV vaccine target through the acquisition of low B cell affinity

Low affinity is common for germline B cell receptors (BCR) seeding development of broadly neutralizing antibodies (bnAbs) that engage hypervariable viruses, including HIV. Antibody affinity selection is also non-homogenizing, insuring the survival of low affinity B cell clones. To explore whether this provides a natural window for expanding human B cell lineages against conserved vaccine targets, we deploy transgenic mice mimicking human antibody diversity and somatic hypermutation (SHM) and immunize with simple monomeric HIV glycoprotein envelope immunogens. We report an immunization regimen that focuses B cell memory upon the conserved CD4 binding site (CD4bs) through both conventional affinity maturation and reproducible expansion of low affinity BCR clones with public patterns in SHM. In the latter instance, SHM facilitates target acquisition by decreasing binding strength. This suggests that permissive B cell selection enables the discovery of antibody epitopes, in this case an HIV bnAb site.

Diversity of B Cell Populations and Ig Repertoire in Human Lungs

The human lung carries a unique microbiome adapted to the air-filled, mucous-lined environment, the presence of which requires an immune system capable of recognizing harmful populations while preventing reactions toward commensals. B cells in the lung play a key role in pulmonary immunity, generating Ag-specific Abs, as well as cytokine secretion for immune activation and regulation. In this study, we compared B cell subsets in human lungs versus circulating cells by analyzing patient-paired lung and blood samples. We found a significantly smaller pool of CD19+, CD20+ B cells in the lung relative to the blood. CD27+, IgD-, class-switched memory B cells (Bmems) composed a larger proportion of the pool of pulmonary B cells. The residency marker CD69 was also significantly higher in the lung. We also sequenced the Ig V region genes (IgVRGs) of class-switched Bmems that do, or do not, express CD69. We observed the IgVRGs of pulmonary Bmems to be as heavily mutated from the unmutated common ancestor as those in circulation. Furthermore, we found progenies within a quasi-clone can gain or lose CD69 expression, regardless of whether the parent clone expressed the residency marker. Overall, our results show that despite its vascularized nature, human lungs carry a unique proportion of B cell subsets. The IgVRGs of pulmonary Bmems are as diverse as those in blood, and progenies of Bmems retain the ability to gain or lose residency.

ARMADiLLO: a web server for analyzing antibody mutation probabilities

Antibodies are generated by B cells that evolve receptor specificity to pathogens through rounds of mutation and selection in a process called affinity maturation. Somatic hypermutation is mediated by an enzyme with DNA sequence context-dependent targeting and substitution resulting in variable probabilities of amino acid substitutions during affinity maturation. We have previously developed a program called Antigen Receptor Mutation Analyzer for the Detection of Low Likelihood Occurrences (ARMADiLLO) that performs simulations of the somatic hypermutation process to estimate the probabilities of observed antibody mutations. Here we describe the ARMADiLLO web server (https://armadillo.dhvi.duke.edu), an easy-to-use web interface that analyzes input antibody sequences and displays the probability estimates for all possible amino acid changes over the full length of an antibody sequence. The probability of antibody mutations can be used by immunologists studying B cell ontogenies and by vaccine designers that are pursuing strategies to elicit broadly neutralizing antibodies which are enriched with developmentally rate-limiting improbable mutations. The ARMADiLLO web server also contains precomputed results reporting the probability of amino acid substitutions in all human V gene segments and in a collection of HIV broadly neutralizing antibodies.

CD62L expression marks SARS-CoV-2 memory B cell subset with preference for neutralizing epitopes

Severe acute respiratory syndrome coronavirus 2-neutralizing antibodies primarily target the spike receptor binding domain (RBD). However, B cell antigen receptors (BCRs) on RBD-binding memory B (B_mem) cells have variation in the neutralizing activities. Here, by combining single B_mem cell profiling with antibody functional assessment, we dissected the phenotype of B_mem cell harboring the potently neutralizing antibodies in coronavirus disease 2019 (COVID-19)-convalescent individuals. The neutralizing subset was marked by an elevated CD62L expression and characterized by distinct epitope preference and usage of convergent V_H (variable region of immunoglobulin heavy chain) genes, accounting for the neutralizing activities. Concordantly, the correlation was observed between neutralizing antibody titers in blood and CD62L⁺ subset, despite the equivalent RBD binding of CD62L⁺ and CD62L^- subset. Furthermore, the kinetics of CD62L⁺ subset differed between the patients who recovered from different COVID-19 severities. Our B_mem cell profiling reveals the unique phenotype of B_mem cell subset that harbors potently neutralizing BCRs, advancing our understanding of humoral protection.

Structural basis for breadth development in the HIV-1 V3-glycan targeting DH270 antibody clonal lineage

Antibody affinity maturation enables adaptive immune responses to a wide range of pathogens. In some individuals broadly neutralizing antibodies develop to recognize rapidly mutating pathogens with extensive sequence diversity. Vaccine design for pathogens such as HIV-1 and influenza has therefore focused on recapitulating the natural affinity maturation process. Here, we determine structures of antibodies in complex with HIV-1 Envelope for all observed members and ancestral states of the broadly neutralizing HIV-1 V3-glycan targeting DH270 antibody clonal B cell lineage. These structures track the development of neutralization breadth from the unmutated common ancestor and define affinity maturation at high spatial resolution. By elucidating contacts mediated by key mutations at different stages of antibody development we identified sites on the epitope-paratope interface that are the focus of affinity optimization. Thus, our results identify bottlenecks on the path to natural affinity maturation and reveal solutions for these that will inform immunogen design aimed at eliciting a broadly neutralizing immune response by vaccination.

Exploring the impact of clonal definition on B-cell diversity: implications for the analysis of immune repertoires

The adaptive immune system has the extraordinary ability to produce a broad range of immunoglobulins that can bind a wide variety of antigens. During adaptive immune responses, activated B cells duplicate and undergo somatic hypermutation in their B-cell receptor (BCR) genes, resulting in clonal families of diversified B cells that can be related back to a common ancestor. Advances in high-throughput sequencing technologies have enabled the high-throughput characterization of B-cell repertoires, however, the accurate identification of clonally related BCR sequences remains a major challenge. In this study, we compare three different clone identification methods on both simulated and experimental data, and investigate their impact on the characterization of B-cell diversity. We observe that different methods lead to different clonal definitions, which affects the quantification of clonal diversity in repertoire data. Our analyses show that direct comparisons between clonal clusterings and clonal diversity of different repertoires should be avoided if different clone identification methods were used to define the clones. Despite this variability, the diversity indices inferred from the repertoires' clonal characterization across samples show similar patterns of variation regardless of the clonal identification method used. We find the Shannon entropy to be the most robust in terms of the variability of diversity rank across samples. Our analysis also suggests that the traditional germline gene alignment-based method for clonal identification remains the most accurate when the complete information about the sequence is known, but that alignment-free methods may be preferred for shorter sequencing read lengths. We make our implementation freely available as a Python library cdiversity.

mRNA vaccines and hybrid immunity use different B cell germlines against Omicron BA.4 and BA.5

Severe acute respiratory syndrome 2 Omicron BA.4 and BA.5 are characterized by high transmissibility and ability to escape natural and vaccine induced immunity. Here we test the neutralizing activity of 482 human monoclonal antibodies isolated from people who received two or three mRNA vaccine doses or from people vaccinated after infection. The BA.4 and BA.5 variants are neutralized only by approximately 15% of antibodies. Remarkably, the antibodies isolated after three vaccine doses target mainly the receptor binding domain Class 1/2, while antibodies isolated after infection recognize mostly the receptor binding domain Class 3 epitope region and the N-terminal domain. Different B cell germlines are used by the analyzed cohorts. The observation that mRNA vaccination and hybrid immunity elicit a different immunity against the same antigen is intriguing and its understanding may help to design the next generation of therapeutics and vaccines against coronavirus disease 2019.

Reconstructing B cell lineage trees with minimum spanning tree and genotype abundances

B cell receptor (BCR) genes exposed to an antigen undergo somatic hypermutations and Darwinian antigen selection, generating a large BCR-antibody diversity. This process, known as B cell affinity maturation, increases antibody affinity, forming a specific B cell lineage that includes the unmutated ancestor and mutated variants. In a B cell lineage, cells with a higher antigen affinity will undergo clonal expansion, while those with a lower affinity will not proliferate and probably be eliminated. Therefore, cellular (genotype) abundance provides a valuable perspective on the ongoing evolutionary process. Phylogenetic tree inference is often used to reconstruct B cell lineage trees and represents the evolutionary dynamic of BCR affinity maturation. However, such methods should process B-cell population data derived from experimental sampling that might contain different cellular abundances. There are a few phylogenetic methods for tracing the evolutionary events occurring in B cell lineages; best-performing solutions are time-demanding and restricted to analysing a reduced number of sequences, while time-efficient methods do not consider cellular abundances. We propose ClonalTree, a low-complexity and accurate approach to construct B-cell lineage trees that incorporates genotype abundances into minimum spanning tree (MST) algorithms. Using both simulated and experimental data, we demonstrate that ClonalTree outperforms MST-based algorithms and achieves a comparable performance to a method that explores tree-generating space exhaustively. Furthermore, ClonalTree has a lower running time, being more convenient for building B-cell lineage trees from high-throughput BCR sequencing data, mainly in biomedical applications, where a lower computational time is appreciable. It is hundreds to thousands of times faster than exhaustive approaches, enabling the analysis of a large set of sequences within minutes or seconds and without loss of accuracy. The source code is freely available at github.com/julibinho/ClonalTree.

A non-neutralizing glycoprotein B monoclonal antibody protects against herpes simplex virus disease in mice

There is an unmet need for monoclonal antibodies (mAbs) for prevention or as adjunctive treatment of herpes simplex virus (HSV) disease. Most vaccine and mAb efforts focus on neutralizing antibodies, but for HSV this strategy has proven ineffective. Preclinical studies with a candidate HSV vaccine strain, ΔgD-2, demonstrated that non-neutralizing antibodies that activate Fcγ receptors (FcγRs) to mediate antibody-dependent cellular cytotoxicity (ADCC) provide active and passive protection against HSV-1 and HSV-2. We hypothesized that this vaccine provides a tool to identify and characterize protective mAbs. We isolated HSV-specific mAbs from germinal center and memory B cells and bone marrow plasmacytes of ΔgD-2-vaccinated mice and evaluated these mAbs for binding, neutralizing, and FcγR-activating activity and for protective efficacy in mice. The most potent protective mAb, BMPC-23, was not neutralizing but activated murine FcγRIV, a biomarker of ADCC. The cryo-electron microscopic structure of the Fab-glycoprotein B (gB) assembly identified domain IV of gB as the epitope. A single dose of BMPC-23 administered 24 hours before or after viral challenge provided significant protection when configured as mouse IgG2c and protected mice expressing human FcγRIII when engineered as a human IgG1. These results highlight the importance of FcR-activating antibodies in protecting against HSV.

Rare, convergent antibodies targeting the stem helix broadly neutralize diverse betacoronaviruses

Humanity has faced three recent outbreaks of novel betacoronaviruses, emphasizing the need to develop approaches that broadly target coronaviruses. Here, we identify 55 monoclonal antibodies from COVID-19 convalescent donors that bind diverse betacoronavirus spike proteins. Most antibodies targeted an S2 epitope that included the K814 residue and were non-neutralizing. However, 11 antibodies targeting the stem helix neutralized betacoronaviruses from different lineages. Eight antibodies in this group, including the six broadest and most potent neutralizers, were encoded by IGHV1-46 and IGKV3-20. Crystal structures of three antibodies of this class at 1.5-1.75-Å resolution revealed a conserved mode of binding. COV89-22 neutralized SARS-CoV-2 variants of concern including Omicron BA.4/5 and limited disease in Syrian hamsters. Collectively, these findings identify a class of IGHV1-46/IGKV3-20 antibodies that broadly neutralize betacoronaviruses by targeting the stem helix but indicate these antibodies constitute a small fraction of the broadly reactive antibody response to betacoronaviruses after SARS-CoV-2 infection.

B cell analyses after SARS-CoV-2 mRNA third vaccination reveals a hybrid immunity like antibody response

The continuous evolution of SARS-CoV-2 generated highly mutated variants able to escape natural and vaccine-induced primary immunity. The administration of a third mRNA vaccine dose induces a secondary response with increased protection. Here we investigate the longitudinal evolution of the neutralizing antibody response in four donors after three mRNA doses at single-cell level. We sorted 4100 spike protein specific memory B cells identifying 350 neutralizing antibodies. The third dose increases the antibody neutralization potency and breadth against all SARS-CoV-2 variants as observed with hybrid immunity. However, the B cell repertoire generating this response is different. The increases of neutralizing antibody responses is largely due to the expansion of B cell germlines poorly represented after two doses, and the reduction of germlines predominant after primary immunization. Our data show that different immunization regimens induce specific molecular signatures which should be considered while designing new vaccines and immunization strategies.

Humanized V(D)J-rearranging and TdT-expressing mouse vaccine models with physiological HIV-1 broadly neutralizing antibody precursors

Antibody heavy chain (HC) and light chain (LC) variable region exons are assembled by V(D)J recombination. V(D)J junctional regions encode complementarity-determining-region 3 (CDR3), an antigen-contact region immensely diversified through nontemplated nucleotide additions ("N-regions") by terminal deoxynucleotidyl transferase (TdT). HIV-1 vaccine strategies seek to elicit human HIV-1 broadly neutralizing antibodies (bnAbs), such as the potent CD4-binding site VRC01-class bnAbs. Mice with primary B cells that express receptors (BCRs) representing bnAb precursors are used as vaccination models. VRC01-class bnAbs uniformly use human HC VH1-2 and commonly use human LCs Vκ3-20 or Vκ1-33 associated with an exceptionally short 5-amino-acid (5-aa) CDR3. Prior VRC01-class models had nonphysiological precursor levels and/or limited precursor diversity. Here, we describe VRC01-class rearranging mice that generate more physiological primary VRC01-class BCR repertoires via rearrangement of VH1-2, as well as Vκ1-33 and/or Vκ3-20 in association with diverse CDR3s. Human-like TdT expression in mouse precursor B cells increased LC CDR3 length and diversity and also promoted the generation of shorter LC CDR3s via N-region suppression of dominant microhomology-mediated Vκ-to-Jκ joins. Priming immunization with eOD-GT8 60mer, which strongly engages VRC01 precursors, induced robust VRC01-class germinal center B cell responses. Vκ3-20-based responses were enhanced by N-region addition, which generates Vκ3-20-to-Jκ junctional sequence combinations that encode VRC01-class 5-aa CDR3s with a critical E residue. VRC01-class-rearranging models should facilitate further evaluation of VRC01-class prime and boost immunogens. These new VRC01-class mouse models establish a prototype for the generation of vaccine-testing mouse models for other HIV-1 bnAb lineages that employ different HC or LC Vs.

A Zika virus-specific IgM elicited in pregnancy exhibits ultrapotent neutralization

Congenital Zika virus (ZIKV) infection results in neurodevelopmental deficits in up to 14% of infants born to ZIKV-infected mothers. Neutralizing antibodies are a critical component of protective immunity. Here, we demonstrate that plasma IgM contributes to ZIKV immunity in pregnancy, mediating neutralization up to 3 months post-symptoms. From a ZIKV-infected pregnant woman, we isolated a pentameric ZIKV-specific IgM (DH1017.IgM) that exhibited ultrapotent ZIKV neutralization dependent on the IgM isotype. DH1017.IgM targets an envelope dimer epitope within domain II. The epitope arrangement on the virion is compatible with concurrent engagement of all ten antigen-binding sites of DH1017.IgM, a solution not available to IgG. DH1017.IgM protected mice against viremia upon lethal ZIKV challenge more efficiently than when expressed as an IgG. Our findings identify a role for antibodies of the IgM isotype in protection against ZIKV and posit DH1017.IgM as a safe and effective candidate immunotherapeutic, particularly during pregnancy.

An epitope-enriched immunogen expands responses to a conserved viral site

Pathogens evade host humoral responses by accumulating mutations in surface antigens. While variable, there are conserved regions that cannot mutate without compromising fitness. Antibodies targeting these conserved epitopes are often broadly protective but remain minor components of the repertoire. Rational immunogen design leverages a structural understanding of viral antigens to modulate humoral responses to favor these responses. Here, we report an epitope-enriched immunogen presenting a higher copy number of the influenza hemagglutinin (HA) receptor-binding site (RBS) epitope relative to other B cell epitopes. Immunization in a partially humanized murine model imprinted with an H1 influenza shows H1-specific serum and >99% H1-specific B cells being RBS-directed. Single B cell analyses show a genetically restricted response that structural analysis defines as RBS-directed antibodies engaging the RBS with germline-encoded contacts. These data show how epitope enrichment expands B cell responses toward conserved epitopes and advances immunogen design approaches for next-generation viral vaccines.

Stabilized HIV-1 envelope immunization induces neutralizing antibodies to the CD4bs and protects macaques against mucosal infection

A successful HIV-1 vaccine will require induction of a polyclonal neutralizing antibody (nAb) response, yet vaccine-mediated induction of such a response in primates remains a challenge. We found that a stabilized HIV-1 CH505 envelope (Env) trimer formulated with a Toll-like receptor 7/8 agonist induced potent HIV-1 polyclonal nAbs that correlated with protection from homologous simian-human immunodeficiency virus (SHIV) infection. The serum dilution that neutralized 50% of virus replication (ID50 titer) required to protect 90% of macaques was 1:364 against the challenge virus grown in primary rhesus CD4+ T cells. Structural analyses of vaccine-induced nAbs demonstrated targeting of the Env CD4 binding site or the N156 glycan and the third variable loop base. Autologous nAb specificities similar to those elicited in macaques by vaccination were isolated from the human living with HIV from which the CH505 Env immunogen was derived. CH505 viral isolates were isolated that mutated the V1 to escape both the infection-induced and vaccine-induced antibodies. These results define the specificities of a vaccine-induced nAb response and the protective titers of HIV-1 vaccine-induced nAbs required to protect nonhuman primates from low-dose mucosal challenge by SHIVs bearing a primary transmitted/founder Env.

Rotavirus VP4 Epitope of a Broadly Neutralizing Human Antibody Defined by Its Structure Bound with an Attenuated-Strain Virion

Rotavirus live-attenuated vaccines, both mono- and pentavalent, generate broadly heterotypic protection. B-cells isolated from adults encode neutralizing antibodies, some with affinity for VP5*, that afford broad protection in mice. We have mapped the epitope of one such antibody by determining the high-resolution cryo-EM structure of its antigen-binding fragment (Fab) bound to the virion of a candidate vaccine strain, CDC-9. The Fab contacts both the distal end of a VP5* β-barrel domain and the two VP8* lectin-like domains at the tip of a projecting spike. Its interactions with VP8* do not impinge on the likely receptor-binding site, suggesting that the mechanism of neutralization is at a step subsequent to initial attachment. We also examined structures of CDC-9 virions from two different stages of serial passaging. Nearly all the VP4 (cleaved to VP8*/VP5*) spikes on particles from the earlier passage (wild-type isolate) had transitioned from the "upright" conformation present on fully infectious virions to the "reversed" conformation that is probably the end state of membrane insertion, unable to mediate penetration, consistent with the very low in vitro infectivity of the wild-type isolate. About half the VP4 spikes were upright on particles from the later passage, which had recovered substantial in vitro infectivity but had acquired an attenuated phenotype in neonatal rats. A mutation in VP4 that occurred during passaging appears to stabilize the interface at the apex of the spike and could account for the greater stability of the upright spikes on the late-passage, attenuated isolate. IMPORTANCE Rotavirus live-attenuated vaccines generate broadly heterotypic protection, and B-cells isolated from adults encode antibodies that are broadly protective in mice. Determining the structural and mechanistic basis of broad protection can contribute to understanding the current limitations of vaccine efficacy in developing countries. The structure of an attenuated human rotavirus isolate (CDC-9) bound with the Fab fragment of a broadly heterotypic protective antibody shows that protection is probably due to inhibition of the conformational transition in the viral spike protein (VP4) critical for viral penetration, rather than to inhibition of receptor binding. A comparison of structures of CDC-9 virus particles at two stages of serial passaging supports a proposed mechanism for initial steps in rotavirus membrane penetration.

Anatomy of Omicron BA.1 and BA.2 neutralizing antibodies in COVID-19 mRNA vaccinees

SARS-CoV-2 vaccines, administered to billions of people worldwide, mitigate the effects of the COVID-19 pandemic, however little is known about the molecular basis of antibody cross-protection to emerging variants, such as Omicron BA.1, its sublineage BA.2, and other coronaviruses. To answer this question, 276 neutralizing monoclonal antibodies (nAbs), previously isolated from seronegative and seropositive donors vaccinated with BNT162b2 mRNA vaccine, were tested for neutralization against the Omicron BA.1 and BA.2 variants, and SARS-CoV-1 virus. Only 14.2, 19.9 and 4.0% of tested antibodies neutralize BA.1, BA.2, and SARS-CoV-1 respectively. These nAbs recognize mainly the SARS-CoV-2 receptor binding domain (RBD) and target Class 3 and Class 4 epitope regions on the SARS-CoV-2 spike protein. Interestingly, around 50% of BA.2 nAbs did not neutralize BA.1 and among these, several targeted the NTD. Cross-protective antibodies derive from a variety of germlines, the most frequents of which were the IGHV1-58;IGHJ3-1, IGHV2-5;IGHJ4-1 and IGHV1-69;IGHV4-1. Only 15.6, 20.3 and 7.8% of predominant gene-derived nAbs elicited against the original Wuhan virus cross-neutralize Omicron BA.1, BA.2 and SARS-CoV-1 respectively. Our data provide evidence, at molecular level, of the presence of cross-neutralizing antibodies induced by vaccination and map conserved epitopes on the S protein that can inform vaccine design.

Antibody repertoire sequencing analysis

High-throughput sequencing for B cell receptor (BCR) repertoire provides useful insights for the adaptive immune system. With the continuous development of the BCR-seq technology, many efforts have been made to develop methods for analyzing the ever-increasing BCR repertoire data. In this review, we comprehensively outline different BCR repertoire library preparation protocols and summarize three major steps of BCR-seq data analysis, i. e., V(D)J sequence annotation, clonal phylogenetic inference, and BCR repertoire profiling and mining. Different from other reviews in this field, we emphasize background intuition and the statistical principle of each method to help biologists better understand it. Finally, we discuss data mining problems for BCR-seq data and with a highlight on recently emerging multiple-sample analysis.

Recall of B cell memory depends on relative locations of prime and boost immunization

Immunization or microbial infection can establish long-term B cell memory not only systemically but also locally. Evidence has suggested that local B cell memory contributes to early local plasmacytic responses after secondary challenge. However, it is unclear whether locality of immunization plays any role in memory B cell participation in recall germinal centers (GCs), which is essential for updating their B cell antigen receptors (BCRs). Using single B cell culture and fate mapping, we have characterized BCR repertoires in recall GCs after boost immunizations at sites local or distal to the priming. Local boosts with homologous antigen recruit the progeny of primary GC B cells to recall GCs more efficiently than do distal boosts. Recall GCs elicited by local boosts contain significantly more B cells with elevated levels of immunoglobulin (Ig) mutation and higher avidity BCRs. This local preference is unaffected by blocking CD40:CD154 interaction to terminate active, GC responses. Local boosts with heterologous antigens elicit secondary GCs with B cell populations enriched for cross-reactivity to the prime and boost antigens; in contrast, cross-reactive GC B cells are rare after distal boosts. Our results suggest that local B cell memory is retained in the form of memory B cells, GC B cells, and GC phenotype B cells that are independent of organized GC structures and that these persistent "primed B cells" contribute to recall GC responses at local sites. Our findings indicate the importance of locality in humoral immunity and inform serial vaccination strategies for evolving viruses.

Development of Neutralization Breadth against Diverse HIV-1 by Increasing Ab-Ag Interface on V2

Understanding maturation pathways of broadly neutralizing antibodies (bnAbs) against HIV-1 can be highly informative for HIV-1 vaccine development. A lineage of J038 bnAbs is now obtained from a long-term SHIV-infected macaque. J038 neutralizes 54% of global circulating HIV-1 strains. Its binding induces a unique "up" conformation for one of the V2 loops in the trimeric envelope glycoprotein and is heavily dependent on glycan, which provides nearly half of the binding surface. Their unmutated common ancestor neutralizes the autologous virus. Continuous maturation enhances neutralization potency and breadth of J038 lineage antibodies via expanding antibody-Env contact areas surrounding the core region contacted by germline-encoded residues. Developmental details and recognition features of J038 lineage antibodies revealed here provide a new pathway for elicitation and maturation of V2-targeting bnAbs.

Affinity maturation for an optimal balance between long-term immune coverage and short-term resource constraints

Humoral immunity relies on the mutation and selection of B cells to better recognize pathogens. This affinity maturation process produces cells with diverse recognition capabilities. Examining optimal immune strategies that maximize the long-term immune coverage at a minimal metabolic cost, we show when the immune system should mount a de novo response rather than rely on existing memory cells. Our theory recapitulates known modes of the B cell response, predicts the empirical form of the distribution of clone sizes, and rationalizes as a trade-off between metabolic and immune costs the antigenic imprinting effects that limit the efficacy of vaccines (original antigenic sin). Our predictions provide a framework to interpret experimental results that could be used to inform vaccination strategies.

In order to target threatening pathogens, the adaptive immune system performs a continuous reorganization of its lymphocyte repertoire. Following an immune challenge, the B cell repertoire can evolve cells of increased specificity for the encountered strain. This process of affinity maturation generates a memory pool whose diversity and size remain difficult to predict. We assume that the immune system follows a strategy that maximizes the long-term immune coverage and minimizes the short-term metabolic costs associated with affinity maturation. This strategy is defined as an optimal decision process on a finite dimensional phenotypic space, where a preexisting population of cells is sequentially challenged with a neutrally evolving strain. We show that the low specificity and high diversity of memory B cells—a key experimental result—can be explained as a strategy to protect against pathogens that evolve fast enough to escape highly potent but narrow memory. This plasticity of the repertoire drives the emergence of distinct regimes for the size and diversity of the memory pool, depending on the density of de novo responding cells and on the mutation rate of the strain. The model predicts power-law distributions of clonotype sizes observed in data and rationalizes antigenic imprinting as a strategy to minimize metabolic costs while keeping good immune protection against future strains.

Different adjuvanted pediatric HIV envelope vaccines induced distinct plasma antibody responses despite similar B cell receptor repertoires in infant rhesus macaques

Different HIV vaccine regimens elicit distinct plasma antibody responses in both human and nonhuman primate models. Previous studies in human and non-human primate infants showed that adjuvants influenced the quality of plasma antibody responses induced by pediatric HIV envelope vaccine regimens. We recently reported that use of the 3M052-SE adjuvant and longer intervals between vaccinations are associated with higher magnitude of antibody responses in infant rhesus macaques. However, the impact of different adjuvants in HIV vaccine regimens on the developing infant B cell receptor (BCR) repertoire has not been studied. This study evaluated whether pediatric HIV envelope vaccine regimens with different adjuvants induced distinct antigen-specific memory B cell repertoires and whether specific immunoglobulin (Ig) immunogenetic characteristics are associated with higher magnitude of plasma antibody responses in vaccinated infant rhesus macaques. We utilized archived preclinical pediatric HIV vaccine studies PBMCs and tissue samples from 19 infant rhesus macaques immunized either with (i) HIV Env protein with a squalene adjuvant, (ii) MVA-HIV and Env protein co-administered using a 3-week interval, (iii) MVA-HIV prime/ protein boost with an extended 6-week interval between immunizations, or (iv) with HIV Env administered with 3M-052-SE adjuvant. Frequencies of vaccine-elicited HIV Env-specific memory B cells from PBMCs and tissues were similar across vaccination groups (frequency range of 0.06–1.72%). There was no association between vaccine-elicited antigen-specific memory B cell frequencies and plasma antibody titer or avidity. Moreover, the epitope specificity and Ig immunogenetic features of vaccine-elicited monoclonal antibodies did not differ between the different vaccine regimens. These data suggest that pediatric HIV envelope vaccine candidates with different adjuvants that previously induced higher magnitude and quality of plasma antibody responses in infant rhesus macaques were not driven by distinct antigen-specific memory BCR repertoires.

Bulk gDNA Sequencing of Antibody Heavy-Chain Gene Rearrangements for Detection and Analysis of B-Cell Clone Distribution: A Method by the AIRR Community

In this method we illustrate how to amplify, sequence, and analyze antibody/immunoglobulin (IG) heavy-chain gene rearrangements from genomic DNA that is derived from bulk populations of cells by next-generation sequencing (NGS). We focus on human source material and illustrate how bulk gDNA-based sequencing can be used to examine clonal architecture and networks in different samples that are sequenced from the same individual. Although bulk gDNA-based sequencing can be performed on both IG heavy (IGH) or kappa/lambda light (IGK/IGL) chains, we focus here on IGH gene rearrangements because IG heavy chains are more diverse, tend to harbor higher levels of somatic hypermutations (SHM), and are more reliable for clone identification and tracking. We also provide a procedure, including code, and detailed instructions for processing and annotation of the NGS data. From these data we show how to identify expanded clones, visualize the overall clonal landscape, and track clonal lineages in different samples from the same individual. This method has a broad range of applications, including the identification and monitoring of expanded clones, the analysis of blood and tissue-based clonal networks, and the study of immune responses including clonal evolution.

Adaptive Immune Receptor Repertoire (AIRR) Community Guide to Repertoire Analysis

Adaptive immune receptor repertoires (AIRRs) are rich with information that can be mined for insights into the workings of the immune system. Gene usage, CDR3 properties, clonal lineage structure, and sequence diversity are all capable of revealing the dynamic immune response to perturbation by disease, vaccination, or other interventions. Here we focus on a conceptual introduction to the many aspects of repertoire analysis and orient the reader toward the uses and advantages of each. Along the way, we note some of the many software tools that have been developed for these investigations and link the ideas discussed to chapters on methods provided elsewhere in this volume.

Hybrid immunity improves B cells and antibodies against SARS-CoV-2 variants

The emergence of SARS-CoV-2 variants is jeopardizing the effectiveness of current vaccines and limiting the application of monoclonal antibody-based therapy for COVID-19 (refs. ^1,2). Here we analysed the memory B cells of five naive and five convalescent people vaccinated with the BNT162b2 mRNA vaccine to investigate the nature of the B cell and antibody response at the single-cell level. Almost 6,000 cells were sorted, over 3,000 cells produced monoclonal antibodies against the spike protein and more than 400 cells neutralized the original SARS-CoV-2 virus first identified in Wuhan, China. The B.1.351 (Beta) and B.1.1.248 (Gamma) variants escaped almost 70% of these antibodies, while a much smaller portion was impacted by the B.1.1.7 (Alpha) and B.1.617.2 (Delta) variants. The overall loss of neutralization was always significantly higher in the antibodies from naive people. In part, this was due to the IGHV2-5;IGHJ4-1 germline, which was found only in people who were convalescent and generated potent and broadly neutralizing antibodies. Our data suggest that people who are seropositive following infection or primary vaccination will produce antibodies with increased potency and breadth and will be able to better control emerging SARS-CoV-2 variants.

Single-cell-level analysis of memory B cells and their response to vaccination against all SARS-CoV-2 variants of concern in individuals who either had or had not been previously exposed to the virus.

The respiratory syncytial virus (RSV) prefusion F-protein functional antibody repertoire in adult healthy donors

Respiratory syncytial virus (RSV) is the leading cause of death from lower respiratory tract infection in infants and children, and is responsible for considerable morbidity and mortality in older adults. Vaccines for pregnant women and elderly which are in phase III clinical studies target people with pre-existing natural immunity against RSV. To investigate the background immunity which will be impacted by vaccination, we single cell-sorted human memory B cells and dissected functional and genetic features of neutralizing antibodies (nAbs) induced by natural infection. Most nAbs recognized both the prefusion and postfusion conformations of the RSV F-protein (cross-binders) while a smaller fraction bound exclusively to the prefusion conformation. Cross-binder nAbs used a wide array of gene rearrangements, while preF-binder nAbs derived mostly from the expansion of B-cell clonotypes from the IGHV1 germline. This latter class of nAbs recognizes an epitope located between Site Ø, Site II, and Site V on the F-protein, identifying an important site of pathogen vulnerability.

Structural and genetic convergence of HIV-1 neutralizing antibodies in vaccinated non-human primates

A primary goal of HIV-1 vaccine development is the consistent elicitation of protective, neutralizing antibodies. While highly similar neutralizing antibodies (nAbs) have been isolated from multiple HIV-infected individuals, it is unclear whether vaccination can consistently elicit highly similar nAbs in genetically diverse primates. Here, we show in three outbred rhesus macaques that immunization with Env elicits a genotypically and phenotypically conserved nAb response. From these vaccinated macaques, we isolated four antibody lineages that had commonalities in immunoglobulin variable, diversity, and joining gene segment usage. Atomic-level structures of the antigen binding fragments of the two most similar antibodies showed nearly identical paratopes. The Env binding modes of each of the four vaccine-induced nAbs were distinct from previously known monoclonal HIV-1 neutralizing antibodies, but were nearly identical to each other. The similarities of these antibodies show that the immune system in outbred primates can respond to HIV-1 Env vaccination with a similar structural and genotypic solution for recognizing a particular neutralizing epitope. These results support rational vaccine design for HIV-1 that aims to reproducibly elicit, in genetically diverse primates, nAbs with specific paratope structures capable of binding conserved epitopes.

A Prevalent Focused Human Antibody Response to the Influenza Virus Hemagglutinin Head Interface

Novel animal influenza viruses emerge, initiate pandemics, and become endemic seasonal variants that have evolved to escape from prevalent herd immunity. These processes often outpace vaccine-elicited protection. Focusing immune responses on conserved epitopes may impart durable immunity. We describe a focused, protective antibody response, abundant in memory and serum repertoires, to a conserved region at the influenza virus hemagglutinin (HA) head interface. Structures of 11 examples, 8 reported here, from seven human donors demonstrate the convergence of responses on a single epitope. The 11 are genetically diverse, with one class having a common, IGκV1-39, light chain. All of the antibodies bind HAs from multiple serotypes. The lack of apparent genetic restriction and potential for elicitation by more than one serotype may explain their abundance. We define the head interface as a major target of broadly protective antibodies with the potential to influence the outcomes of influenza virus infection. IMPORTANCE The rapid appearance of mutations in circulating human influenza viruses and selection for escape from herd immunity require prediction of likely variants for an annual updating of influenza vaccines. The identification of human antibodies that recognize conserved surfaces on the influenza virus hemagglutinin (HA) has prompted efforts to design immunogens that might selectively elicit such antibodies. The recent discovery of a widely prevalent antibody response to the conserved interface between two HA "heads" (the globular, receptor-binding domains at the apex of the spike-like trimer) has added a new target for these efforts. We report structures of eight such antibodies, bound with HA heads, and compare them with each other and with three others previously described. Although genetically diverse, they all converge on a common binding site. The analysis here can guide immunogen design for preclinical trials.

Fab-dimerized glycan-reactive antibodies are a structural category of natural antibodies

Natural antibodies (Abs) can target host glycans on the surface of pathogens. We studied the evolution of glycan-reactive B cells of rhesus macaques and humans using glycosylated HIV-1 envelope (Env) as a model antigen. 2G12 is a broadly neutralizing Ab (bnAb) that targets a conserved glycan patch on Env of geographically diverse HIV-1 strains using a unique heavy-chain (VH) domain-swapped architecture that results in fragment antigen-binding (Fab) dimerization. Here, we describe HIV-1 Env Fab-dimerized glycan (FDG)-reactive bnAbs without VH-swapped domains from simian-human immunodeficiency virus (SHIV)-infected macaques. FDG Abs also recognized cell-surface glycans on diverse pathogens, including yeast and severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike. FDG precursors were expanded by glycan-bearing immunogens in macaques and were abundant in HIV-1-naive humans. Moreover, FDG precursors were predominately mutated IgM+IgD+CD27+, thus suggesting that they originated from a pool of antigen-experienced IgM+ or marginal zone B cells.

Functional Homology for Antibody-Dependent Phagocytosis Across Humans and Rhesus Macaques

Analyses of human clinical HIV-1 vaccine trials and preclinical vaccine studies performed in rhesus macaque (RM) models have identified associations between non-neutralizing Fc Receptor (FcR)-dependent antibody effector functions and reduced risk of infection. Specifically, antibody-dependent phagocytosis (ADP) has emerged as a common correlate of reduced infection risk in multiple RM studies and the human HVTN505 trial. This recurrent finding suggests that antibody responses with the capability to mediate ADP are most likely a desirable component of vaccine responses aimed at protecting against HIV-1 acquisition. As use of RM models is essential for development of the next generation of candidate HIV-1 vaccines, there is a need to determine how effectively ADP activity observed in RMs translates to activity in humans. In this study we compared ADP activity of human and RM monocytes and polymorphonuclear leukocytes (PMN) to bridge this gap in knowledge. We observed considerable variability in the magnitude of monocyte and PMN ADP activity across individual humans and RM that was not dependent on FcR alleles, and only modestly impacted by cell-surface levels of FcRs. Importantly, we found that for both human and RM phagocytes, ADP activity of antibodies targeting the CD4 binding site was greatest when mediated by human IgG3, followed by RM and human IgG1. These results demonstrate that there is functional homology between antibody and FcRs from these two species for ADP. We also used novel RM IgG1 monoclonal antibodies engineered with elongated hinge regions to show that hinge elongation augments RM ADP activity. The RM IgGs with engineered hinge regions can achieve ADP activity comparable to that observed with human IgG3. These novel modified antibodies will have utility in passive immunization studies aimed at defining the role of IgG3 and ADP in protection from virus challenge or control of disease in RM models. Our results contribute to a better translation of human and macaque antibody and FcR biology, and may help to improve testing accuracy and evaluations of future active and passive prevention strategies.

Mining HIV controllers for broad and functional antibodies to recognize and eliminate HIV-infected cells

HIV monoclonal antibodies for viral reservoir eradication strategies will likely need to recognize reactivated infected cells and potently drive Fc-mediated innate effector cell activity. We systematically characterize a library of 185 HIV-envelope-specific antibodies derived from 15 spontaneous HIV controllers (HCs) that selectively exhibit robust serum Fc functionality and compared them to broadly neutralizing antibodies (bNAbs) in clinical development. Within the 10 antibodies with the broadest cell-recognition capability, seven originated from HCs and three were bNAbs. V3-loop-targeting antibodies are enriched among the top cell binders, suggesting the V3-loop may be selectively exposed and accessible on the cell surface. Fc functionality is more variable across antibodies, which is likely influenced by distinct binding topology and corresponding Fc accessibility, highlighting not only the importance of target-cell recognition but also the need to optimize for Fc-mediated elimination. Ultimately, our results demonstrate that this comprehensive selection process can identify monoclonal antibodies poised to eliminate infected cells.

Rossignol et al. characterize 185 HIV-envelope-specific antibodies derived from spontaneous HIV controllers, downselecting antibodies based on their ability to broadly recognize infected cells and potently drive Fc-mediated innate effector cell activity. This comprehensive selection process can identify monoclonal antibodies poised to eliminate infected cells for viral reservoir eradication strategies.

Alignment free identification of clones in B cell receptor repertoires

Following antigenic challenge, activated B cells rapidly expand and undergo somatic hypermutation, yielding groups of clonally related B cells with diversified immunoglobulin receptors. Inference of clonal relationships based on the receptor sequence is an essential step in many adaptive immune receptor repertoire sequencing studies. These relationships are typically identified by a multi-step process that involves: (i) grouping sequences based on shared V and J gene assignments, and junction lengths and (ii) clustering these sequences using a junction-based distance. However, this approach is sensitive to the initial gene assignments, which are error-prone, and fails to identify clonal relatives whose junction length has changed through accumulation of indels. Through defining a translation-invariant feature space in which we cluster the sequences, we develop an alignment free clonal identification method that does not require gene assignments and is not restricted to a fixed junction length. This alignment free approach has higher sensitivity compared to a typical junction-based distance method without loss of specificity and PPV. While the alignment free procedure identifies clones that are broadly consistent with the junction-based distance method, it also identifies clones with characteristics (multiple V or J gene assignments or junction lengths) that are not detectable with the junction-based distance method.

Correlations in Somatic Hypermutation Between Sites in IGHV Genes Can Be Explained by Interactions Between AID and/or Polη Hotspots

The somatic hypermutation (SHM) of Immunoglobulin (Ig) genes is a key process during antibody affinity maturation in B cells. The mutagenic enzyme activation induced deaminase (AID) is required for SHM and has a preference for WRC hotspots in DNA. Error-prone repair mechanisms acting downstream of AID introduce further mutations, including DNA polymerase eta (Polη), part of the non-canonical mismatch repair pathway (ncMMR), which preferentially generates mutations at WA hotspots. Previously proposed mechanistic models lead to a variety of predictions concerning interactions between hotspots, for example, how mutations in one hotspot will affect another hotspot. Using a large, high-quality, Ig repertoire sequencing dataset, we evaluated pairwise correlations between mutations site-by-site using an unbiased measure similar to mutual information which we termed “mutational association” (MA). Interactions are dominated by relatively strong correlations between nearby sites (short-range MAs), which can be almost entirely explained by interactions between overlapping hotspots for AID and/or Polη. We also found relatively weak dependencies between almost all sites throughout each gene (longer-range MAs), although these arise mostly as a statistical consequence of high pairwise mutation frequencies. The dominant short-range interactions are also highest within the most highly mutating IGHV sub-regions, such as the complementarity determining regions (CDRs), where there is a high hotspot density. Our results suggest that the hotspot preferences for AID and Polη have themselves evolved to allow for greater interactions between AID and/or Polη induced mutations.

The respiratory syncytial virus fusion protein-specific B cell receptor repertoire reshaped by post-fusion subunit vaccination

Respiratory syncytial virus (RSV) is the major cause of acute lower respiratory illness in children of less than 5 years of age which usually results in hospitalization or even in death. Vaccine development is hampered in consequence of a failed vaccine trial with fatalities in the 1960s. Even though research has been more focused on the RSV fusion protein in its pre-fusion conformation, maternal vaccination with post-fusion protein (post F) was considered as a promising vaccine strategy for passive immunization of babies, because post F preserves very potent neutralizing epitopes. We extensively analyzed post F-binding B cell receptor (BCR) repertoires of three vaccinees who received a post F-subunit vaccine in the context of a first-in-human, Phase 1, randomized, observer-blind, placebo-controlled clinical trial (ClinicalTrials.gov Identifier: NCT02298179). In order to compare the vaccine-induced BCR repertoires with BCR repertoires induced by natural infection, we also analyzed pre F- and post F-binding BCRs isolated from a healthy blood donor with relatively high F-binding memory B cell (MBC) frequencies. Analysis of the vaccine-induced repertoires revealed that preferentially V_H4-encoded BCRs were expanded in response to vaccination. Estimation of antigen-driven selection further demonstrated that expanded BCRs accumulated positively selected replacement mutations which substantiated the hypothesis that post F-vaccination induces diversification of V_H4-encoded BCRs in germinal centers. Comparison of the vaccine-induced BCR repertoires with clonally related pre and post F-binding BCRs of the healthy blood donor suggested that the vaccine expanded pre/post F cross-reactive MBCs. Interestingly, several vaccine-induced BCRs shared stereotypic VDJ gene junctions with known neutralizing Abs. Once expressed for functional characterization, the selected monoclonal Abs demonstrated the predicted neutralization activities in plaque reduction neutralization assays indicating that the post F-vaccine induced expansion of neutralizing BCRs.

Dataset of antibody variable region sequence features inferred from a respiratory syncytial virus fusion protein-specific B cell receptor repertoire induced by natural infection of a healthy adult

Respiratory syncytial virus (RSV) is the primary cause for acute lower respiratory syndrome in children younger than 5 years. Research on B cell repertoires and antibodies binding the RSV fusion protein (RSV F) is of major interest in the development of potential vaccine candidates and therapies. B cell receptors (BCRs) which have higher affinities for a specific antigen are preferentially selected for B cell clonal expansion in germinal center reactions. Consequently, antigen-specific BCR repertoires share common features, as for instance preferential variable gene usage, variable region mutation levels or lengths of the heavy chain complementarity-determining region 3. Since RSV repeatedly infects every person throughout life, memory B cells (MBC) expressing RSV F-binding BCRs circulate in the blood of healthy adults. This dataset of BCR variable region sequence features was derived from single cell-sorted RSV F-directed MBCs of a healthy adult blood donor [1]. The dataset was produced with publicly available data analysis software programs and scripts, which facilitates integration or comparison with antibody sequence repertoire data of different individuals derived with the same or comparable data analysis approaches and tools.

A Bayesian phylogenetic hidden Markov model for B cell receptor sequence analysis

The human body generates a diverse set of high affinity antibodies, the soluble form of B cell receptors (BCRs), that bind to and neutralize invading pathogens. The natural development of BCRs must be understood in order to design vaccines for highly mutable pathogens such as influenza and HIV. BCR diversity is induced by naturally occurring combinatorial "V(D)J" rearrangement, mutation, and selection processes. Most current methods for BCR sequence analysis focus on separately modeling the above processes. Statistical phylogenetic methods are often used to model the mutational dynamics of BCR sequence data, but these techniques do not consider all the complexities associated with B cell diversification such as the V(D)J rearrangement process. In particular, standard phylogenetic approaches assume the DNA bases of the progenitor (or "naive") sequence arise independently and according to the same distribution, ignoring the complexities of V(D)J rearrangement. In this paper, we introduce a novel approach to Bayesian phylogenetic inference for BCR sequences that is based on a phylogenetic hidden Markov model (phylo-HMM). This technique not only integrates a naive rearrangement model with a phylogenetic model for BCR sequence evolution but also naturally accounts for uncertainty in all unobserved variables, including the phylogenetic tree, via posterior distribution sampling.

Antibody Feedback Limits the Expansion of B Cell Responses to Malaria Vaccination but Drives Diversification of the Humoral Response

Generating sufficient antibody to block infection is a key challenge for vaccines against malaria. Here, we show that antibody titers to a key target, the repeat region of the Plasmodium falciparum circumsporozoite protein (PfCSP), plateaued after two immunizations in a clinical trial of the radiation-attenuated sporozoite vaccine. To understand the mechanisms limiting vaccine responsiveness, we developed immunoglobulin (Ig)-knockin mice with elevated numbers of PfCSP-binding B cells. We determined that recall responses were inhibited by antibody feedback, potentially via epitope masking of the immunodominant PfCSP repeat region. Importantly, the amount of antibody that prevents boosting is below the amount of antibody required for protection. Finally, while antibody feedback limited responses to the PfCSP repeat region in vaccinated volunteers, potentially protective subdominant responses to PfCSP C-terminal regions expanded with subsequent boosts. These data suggest that antibody feedback drives the diversification of immune responses and that vaccination for malaria will require targeting multiple antigens.

AncesTree: An interactive immunoglobulin lineage tree visualizer

High-throughput sequencing of human immunoglobulin genes allows analysis of antibody repertoires and the reconstruction of clonal lineage evolution. The study of antibodies (Abs) affinity maturation is of specific interest to understand the generation of Abs with high affinity or broadly neutralizing activities. Moreover, phylogenic analysis enables the identification of the key somatic mutations required to achieve optimal antigen binding. The Immcantation framework provides a start-to-finish set of analytical methods for high-throughput adaptive immune receptor repertoire sequencing (AIRR-Seq; Rep-Seq) data. Furthermore, Immcantation's Change-O package has developed IgPhyML, an algorithm designed to build specifically immunoglobulin (Ig) phylogenic trees. Meanwhile Phylip, an algorithm that has been originally developed for applications in ecology and macroevolution, can also be used for the phylogenic reconstruction of antibodies maturation pathway. To complement Ig lineages made by IgPhyML or Dnaml (Phylip), we developed AncesTree, a graphic user interface (GUI) that aims to give researchers the opportunity to interactively explore antibodies clonal evolution. AncesTree displays interactive immunoglobulins phylogenic tree, Ig related mutations and sequence alignments using additional information coming from specialized antibody tools. The GUI is a Java standalone application allowing interaction with Ig tree that can run under Windows, Linux and Mac OS.

Somatic hypermutation analysis for improved identification of B cell clonal families from next-generation sequencing data

Adaptive immune receptor repertoire sequencing (AIRR-Seq) offers the possibility of identifying and tracking B cell clonal expansions during adaptive immune responses. Members of a B cell clone are descended from a common ancestor and share the same initial V(D)J rearrangement, but their B cell receptor (BCR) sequence may differ due to the accumulation of somatic hypermutations (SHMs). Clonal relationships are learned from AIRR-seq data by analyzing the BCR sequence, with the most common methods focused on the highly diverse junction region. However, clonally related cells often share SHMs which have been accumulated during affinity maturation. Here, we investigate whether shared SHMs in the V and J segments of the BCR can be leveraged along with the junction sequence to improve the ability to identify clonally related sequences. We develop independent distance functions that capture junction similarity and shared mutations, and combine these in a spectral clustering framework to infer the BCR clonal relationships. Using both simulated and experimental data, we show that this model improves both the sensitivity and specificity for identifying B cell clones. Source code for this method is freely available in the SCOPer (Spectral Clustering for clOne Partitioning) R package (version 0.2 or newer) in the Immcantation framework: www.immcantation.org under the AGPLv3 license.

Neonatal Rhesus Macaques Have Distinct Immune Cell Transcriptional Profiles following HIV Envelope Immunization

HIV-1-infected infants develop broadly neutralizing antibodies (bnAbs) more rapidly than adults, suggesting differences in the neonatal versus adult responses to the HIV-1 envelope (Env). Here, trimeric forms of HIV-1 Env immunogens elicit increased gp120- and gp41-specific antibodies more rapidly in neonatal macaques than adult macaques. Transcriptome analyses of neonatal versus adult immune cells after Env vaccination reveal that neonatal macaques have higher levels of the apoptosis regulator BCL2 in T cells and lower levels of the immunosuppressive interleukin-10 (IL-10) receptor alpha (IL10RA) mRNA transcripts in T cells, B cells, natural killer (NK) cells, and monocytes. In addition, immunized neonatal macaques exhibit increased frequencies of activated blood T follicular helper-like (Tfh) cells compared to adults. Thus, neonatal macaques have transcriptome signatures of decreased immunosuppression and apoptosis compared with adult macaques, providing an immune landscape conducive to early-life immunization prior to sexual debut.