Effects of Darwinian Selection and Mutability on Rate of Broadly Neutralizing Antibody Evolution during HIV-1 Infection

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.

Immunoglobulin VDJ repertoires reveal hallmarks of germinal centers in unique cell clusters isolated from zebrafish (Danio rerio) lymphoid tissues

DNA mutagenesis during antibody affinity maturation has potentially oncogenic or autoimmune outcomes if not tightly controlled as it is in mammalian germinal centers. Cold blooded vertebrates lack germinal centers, yet have a functional Ig gene mutator enzyme, Aicda. In fish there are clusters of Aicda⁺ cells encircled by pigmented 'melano-macrophages' and we test the hypothesis that these clusters are functionally analogous to germinal centers. Sequenced IgH VDJ repertoire libraries from individual isolated clusters showed evidence of B-cell clonal expansion and VDJ somatic hypermutation. Construction of Ig clonal lineage trees revealed that unlike surrounding lymphoid tissue, each cluster is dominated by a few B-cell VDJ clonotypes having hundreds of mutated variants. Recruitment of B-cells to the clusters appears to be ongoing, as there are additional Ig clones having smaller lineages. Finally, we show evidence for positive selection for replacement mutations in regions encoding the antigen contact loops, but not in the framework regions, consistent with functional antibody modification. Melano-macrophages appear to trap the Ag used for post-mutation B-cell selection, performing a role analogous to the follicular dendritic cells of mammalian germinal centers. These findings provide insights into the evolution of the affinity maturation process, the improvement of fish vaccines and possibly also the workings of atypical ectopic germinal centers generated in several human diseases.

Human B cell lineages associated with germinal centers following influenza vaccination are measurably evolving

The poor efficacy of seasonal influenza virus vaccines is often attributed to pre-existing immunity interfering with the persistence and maturation of vaccine-induced B cell responses. We previously showed that a subset of vaccine-induced B cell lineages are recruited into germinal centers (GCs) following vaccination, suggesting that affinity maturation of these lineages against vaccine antigens can occur. However, it remains to be determined whether seasonal influenza vaccination stimulates additional evolution of vaccine-specific lineages, and previous work has found no significant increase in somatic hypermutation among influenza-binding lineages sampled from the blood following seasonal vaccination in humans. Here, we investigate this issue using a phylogenetic test of measurable immunoglobulin sequence evolution. We first validate this test through simulations and survey measurable evolution across multiple conditions. We find significant heterogeneity in measurable B cell evolution across conditions, with enrichment in primary response conditions such as HIV infection and early childhood development. We then show that measurable evolution following influenza vaccination is highly compartmentalized: while lineages in the blood are rarely measurably evolving following influenza vaccination, lineages containing GC B cells are frequently measurably evolving. Many of these lineages appear to derive from memory B cells. We conclude from these findings that seasonal influenza virus vaccination can stimulate additional evolution of responding B cell lineages, and imply that the poor efficacy of seasonal influenza vaccination is not due to a complete inhibition of vaccine-specific B cell evolution.

Genome edited B cells: a new frontier in immune cell therapies

Cell therapies based on reprogrammed adaptive immune cells have great potential as "living drugs." As first demonstrated clinically for engineered chimeric antigen receptor (CAR) T cells, the ability of such cells to undergo clonal expansion in response to an antigen promotes both self-renewal and self-regulation in vivo. B cells also have the potential to be developed as immune cell therapies, but engineering their specificity and functionality is more challenging than for T cells. In part, this is due to the complexity of the immunoglobulin (Ig) locus, as well as the requirement for regulated expression of both cell surface B cell receptor and secreted antibody isoforms, in order to fully recapitulate the features of natural antibody production. Recent advances in genome editing are now allowing reprogramming of B cells by site-specific engineering of the Ig locus with preformed antibodies. In this review, we discuss the potential of engineered B cells as a cell therapy, the challenges involved in editing the Ig locus and the advances that are making this possible, and envision future directions for this emerging field of immune cell engineering.

Coevolutionary transitions emerging from flexible molecular recognition and eco-evolutionary feedback

Highly mutable viruses evolve to evade host immunity that exerts selective pressure and adapts to viral dynamics. Here, we provide a framework for identifying key determinants of the mode and fate of viral-immune coevolution by linking molecular recognition and eco-evolutionary dynamics. We find that conservation level and initial diversity of antigen jointly determine the timing and efficacy of narrow and broad antibody responses, which in turn control the transition between viral persistence, clearance, and rebound. In particular, clearance of structurally complex antigens relies on antibody evolution in a larger antigenic space than where selection directly acts; viral rebound manifests binding-mediated feedback between ecology and rapid evolution. Finally, immune compartmentalization can slow viral escape but also delay clearance. This work suggests that flexible molecular binding allows a plastic phenotype that exploits potentiating neutral variations outside direct contact, opening new and shorter paths toward highly adaptable states.

•

A scale-crossing framework identifies key determinants of viral-immune coevolution

•

Fast specific response influences slow broad response by shaping antigen dynamics

•

Antibody footprint shift enables breadth acquisition and viral clearance

•

Model explains divergent kinetics and outcomes of HCV infection in humans

Prospective Artificial Intelligence to Dissect the Dengue Immune Response and Discover Therapeutics

Dengue virus (DENV) poses a serious threat to global health as the causative agent of dengue fever. The virus is endemic in more than 128 countries resulting in approximately 390 million infection cases each year. Currently, there is no approved therapeutic for treatment nor a fully efficacious vaccine. The development of therapeutics is confounded and hampered by the complexity of the immune response to DENV, in particular to sequential infection with different DENV serotypes (DENV1-5). Researchers have shown that the DENV envelope (E) antigen is primarily responsible for the interaction and subsequent invasion of host cells for all serotypes and can elicit neutralizing antibodies in humans. The advent of high-throughput sequencing and the rapid advancements in computational analysis of complex data, has provided tools for the deconvolution of the DENV immune response. Several types of complex statistical analyses, machine learning models and complex visualizations can be applied to begin answering questions about the B- and T-cell immune responses to multiple infections, antibody-dependent enhancement, identification of novel therapeutics and advance vaccine research.

Distinct clonal evolution of B-cells in HIV controllers with neutralizing antibody breadth

A minor subset of individuals infected with HIV-1 develop antibody neutralization breadth during the natural course of the infection, often linked to chronic, high-level viremia. Despite significant efforts, vaccination strategies have been unable to induce similar neutralization breadth and the mechanisms underlying neutralizing antibody induction remain largely elusive. Broadly neutralizing antibody responses can also be found in individuals who control HIV to low and even undetectable plasma levels in the absence of antiretroviral therapy, suggesting that high antigen exposure is not a strict requirement for neutralization breadth. We therefore performed an analysis of paired heavy and light chain B-cell receptor (BCR) repertoires in 12,591 HIV-1 envelope-specific single memory B-cells to determine alterations in the BCR immunoglobulin gene repertoire and B-cell clonal expansions that associate with neutralizing antibody breadth in 22 HIV controllers. We found that the frequency of genomic mutations in IGHV and IGLV was directly correlated with serum neutralization breadth. The repertoire of the most mutated antibodies was dominated by a small number of large clones with evolutionary signatures suggesting that these clones had reached peak affinity maturation. These data demonstrate that even in the setting of low plasma HIV antigenemia, similar to what a vaccine can potentially achieve, BCR selection for extended somatic hypermutation and clonal evolution can occur in some individuals suggesting that host-specific factors might be involved that could be targeted with future vaccine strategies.

Vaccine targeting SIVmac251 protease cleavage sites protects macaques against vaginal infection

After over 3 decades of research, an effective anti-HIV vaccine remains elusive. The recently halted HVTN702 clinical trial not only further stresses the challenge to develop an effective HIV vaccine but also emphasizes that unconventional and novel vaccine strategies are urgently needed. Here, we report that a vaccine focusing the immune response on the sequences surrounding the 12 viral protease cleavage sites (PCSs) provided greater than 80% protection to Mauritian cynomolgus macaques against repeated intravaginal SIVmac251 challenges. The PCS-specific T cell responses correlated with vaccine efficacy. The PCS vaccine did not induce immune activation or inflammation known to be associated with increased susceptibility to HIV infection. Machine learning analyses revealed that the immune microenvironment generated by the PCS vaccine was predictive of vaccine efficacy. Our study demonstrates, for the first time to our knowledge, that a vaccine which targets only viral maturation, but lacks full-length Env and Gag immunogens, can prevent intravaginal infection in a stringent macaque/SIV challenge model. Targeting HIV maturation thus offers a potentially novel approach to developing an effective HIV vaccine.

Shaping Polyclonal Responses via Antigen-Mediated Antibody Interference

Broadly neutralizing antibodies (bnAbs) recognize conserved features of rapidly mutating pathogens and confer universal protection, but they emerge rarely in natural infection. Increasing evidence indicates that seemingly passive antibodies may interfere with natural selection of B cells. Yet, how such interference modulates polyclonal responses is unknown. Here we provide a framework for understanding the role of antibody interference—mediated by multi-epitope antigens—in shaping B cell clonal makeup and the fate of bnAb lineages. We find that, under heterogeneous interference, clones with different intrinsic fitness can collectively persist. Furthermore, antagonism among fit clones (specific for variable epitopes) promotes expansion of unfit clones (targeting conserved epitopes), at the cost of repertoire potency. This trade-off, however, can be alleviated by synergy toward the unfit. Our results provide a physical basis for antigen-mediated clonal interactions, stress system-level impacts of molecular synergy and antagonism, and offer principles to amplify naturally rare clones.

•

Multi-epitope antigens mediate antibody interference that couples B cell lineages

•

Trade-off exists between repertoire potency and persistence of broad lineages

•

Antigen-mediated synergy toward intrinsically unfit clones alleviates the trade-off

•

Amplifying rare clones by leveraging molecular interference structure

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

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

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

INTRODUCTION

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

RATIONALE

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

RESULTS

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

CONCLUSION

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

Repertoire-wide phylogenetic models of B cell molecular evolution reveal evolutionary signatures of aging and vaccination

In order to produce effective antibodies, B cells undergo rapid somatic hypermutation (SHM) and selection for binding affinity to antigen via a process called affinity maturation. The similarities between this process and evolution by natural selection have led many groups to use phylogenetic methods to characterize the development of immunological memory, vaccination, and other processes that depend on affinity maturation. However, these applications are limited by the fact that most phylogenetic models are designed to be applied to individual lineages comprising genetically diverse sequences, while B cell repertoires often consist of hundreds to thousands of separate low-diversity lineages. Further, several features of affinity maturation violate important assumptions in standard phylogenetic models. Here, we introduce a hierarchical phylogenetic framework that integrates information from all lineages in a repertoire to more precisely estimate model parameters while simultaneously incorporating the unique features of SHM. We demonstrate the power of this repertoire-wide approach by characterizing previously undescribed phenomena in affinity maturation. First, we find evidence consistent with age-related changes in SHM hot-spot targeting. Second, we identify a consistent relationship between increased tree length and signs of increased negative selection, apparent in the repertoires of recently vaccinated subjects and those without any known recent infections or vaccinations. This suggests that B cell lineages shift toward negative selection over time as a general feature of affinity maturation. Our study provides a framework for undertaking repertoire-wide phylogenetic testing of SHM hypotheses and provides a means of characterizing dynamics of mutation and selection during affinity maturation.

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

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

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

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

Predicting B cell receptor substitution profiles using public repertoire data

B cells develop high affinity receptors during the course of affinity maturation, a cyclic process of mutation and selection. At the end of affinity maturation, a number of cells sharing the same ancestor (i.e. in the same “clonal family”) are released from the germinal center; their amino acid frequency profile reflects the allowed and disallowed substitutions at each position. These clonal-family-specific frequency profiles, called “substitution profiles”, are useful for studying the course of affinity maturation as well as for antibody engineering purposes. However, most often only a single sequence is recovered from each clonal family in a sequencing experiment, making it impossible to construct a clonal-family-specific substitution profile. Given the public release of many high-quality large B cell receptor datasets, one may ask whether it is possible to use such data in a prediction model for clonal-family-specific substitution profiles. In this paper, we present the method “Substitution Profiles Using Related Families” (SPURF), a penalized tensor regression framework that integrates information from a rich assemblage of datasets to predict the clonal-family-specific substitution profile for any single input sequence. Using this framework, we show that substitution profiles from similar clonal families can be leveraged together with simulated substitution profiles and germline gene sequence information to improve prediction. We fit this model on a large public dataset and validate the robustness of our approach on two external datasets. Furthermore, we provide a command-line tool in an open-source software package (https://github.com/krdav/SPURF) implementing these ideas and providing easy prediction using our pre-fit models.

Antibody engineering can be greatly informed by knowledge about the underlying affinity maturation process. As such this can be probed by sequencing, but unfortunately, in practice often only one member of the clonal family is sequenced, making it difficult to determine a set of possible amino acid mutations that would retain the original antibody antigen binding affinity. We overcome this data sparsity by developing a statistical learning approach that leverages vast information about amino acid preferences available in public immune system repertoire data. We use a penalized regression approach to devise a flexible statistical model that integrates multiple sources of information into a coherent prediction framework and validate our prediction algorithm using subsampling and held out data.

The expanding array of HIV broadly neutralizing antibodies

A large array of broadly neutralizing antibodies (bnAbs) against HIV have been isolated and described, particularly in the last decade. This continually expanding array of bnAbs has crucially led to the identification of novel epitopes on the HIV envelope protein via which antibodies can block a broad range of HIV strains. Moreover, these studies have produced high-resolution understanding of these sites of vulnerability on the envelope protein. They have also clarified the mechanisms of action of bnAbs and provided detailed descriptions of B cell ontogenies from which they arise. However, it is still not possible to predict which HIV-infected individuals will go onto develop breath nor is it possible to induce neutralization breadth by immunization in humans. This review aims to discuss the major insights gained so far and also to evaluate the requirement to continue isolating and characterizing new bnAbs. While new epitopes may remain to be uncovered, a clearer probable benefit of further bnAb characterization is a greater understanding of key decision points in bnAb development within the anti-HIV immune response. This in turn may lead to new insights into how to trigger bnAbs by immunization and more clearly define the challenges to using bnAbs as therapeutic agents.

Sequencing HIV-neutralizing antibody exons and introns reveals detailed aspects of lineage maturation

The developmental pathways of broadly neutralizing antibodies (bNAbs) against HIV are of great importance for the design of immunogens that can elicit protective responses. Here we show the maturation features of the HIV-neutralizing anti-V1V2 VRC26 lineage by simultaneously sequencing the exon together with the downstream intron of VRC26 members. Using the mutational landscapes of both segments and the selection-free nature of the intron region, we identify multiple events of amino acid mutational convergence in the complementarity-determining region 3 (CDR3) of VRC26 members, and determine potential intermediates with diverse CDR3s to a late stage bNAb from 2 years prior to its isolation. Moreover, we functionally characterize the earliest neutralizing intermediates with critical CDR3 mutations, with some emerging only 14 weeks after initial lineage detection and containing only ~6% V gene mutations. Our results thus underscore the utility of analyzing exons and introns simultaneously for studying antibody maturation and repertoire selection.

Knowledge on how antibody responses have evolved is critical for the induction of protective immunity. Here the authors analyse, using high-throughput sequencing of both exon and intron regions, the mutation and lineage development of an HIV-neutralizing antibody to find an unexpected early emergence of broadly neutralizing species.

Mucosal antibody responses to vaccines targeting SIV protease cleavage sites or full-length Gag and Env proteins in Mauritian cynomolgus macaques

HIV mutates rapidly and infects CD4⁺ T cells, especially when they are activated. A vaccine targeting conserved, essential viral elements while limiting CD4⁺ T cell activation could be effective. Learning from natural immunity observed in a group of highly HIV-1 exposed seronegative Kenyan female sex workers, we are testing a novel candidate HIV vaccine targeting the 12 viral protease cleavage sites (PCSs) (the PCS vaccine), in comparison with a vaccine targeting full-length Gag and Env (the Gag/Env vaccine) in a Mauritian cynomolgus macaque/SIV model. In this study we evaluated these vaccines for induction of mucosal antibodies to SIV immunogens at the female genital tract. Bio-Plex and Western blot analyses of cervicovaginal lavage samples showed that both the PCS and Gag/Env vaccines can elicit mucosal IgG antibody responses to SIV immunogens. Significantly higher increase of anti-PCS antibodies was induced by the PCS vaccine than by the Gag/Env vaccine (p<0.0001). The effect of the mucosal antibody responses in protection from repeated low dose pathogenic SIVmac251 challenges is being evaluated.

Computational Strategies for Dissecting the High-Dimensional Complexity of Adaptive Immune Repertoires

The adaptive immune system recognizes antigens via an immense array of antigen-binding antibodies and T-cell receptors, the immune repertoire. The interrogation of immune repertoires is of high relevance for understanding the adaptive immune response in disease and infection (e.g., autoimmunity, cancer, HIV). Adaptive immune receptor repertoire sequencing (AIRR-seq) has driven the quantitative and molecular-level profiling of immune repertoires, thereby revealing the high-dimensional complexity of the immune receptor sequence landscape. Several methods for the computational and statistical analysis of large-scale AIRR-seq data have been developed to resolve immune repertoire complexity and to understand the dynamics of adaptive immunity. Here, we review the current research on (i) diversity, (ii) clustering and network, (iii) phylogenetic, and (iv) machine learning methods applied to dissect, quantify, and compare the architecture, evolution, and specificity of immune repertoires. We summarize outstanding questions in computational immunology and propose future directions for systems immunology toward coupling AIRR-seq with the computational discovery of immunotherapeutics, vaccines, and immunodiagnostics.

HIV Envelope Glycoform Heterogeneity and Localized Diversity Govern the Initiation and Maturation of a V2 Apex Broadly Neutralizing Antibody Lineage

Understanding how broadly neutralizing antibodies (bnAbs) to HIV envelope (Env) develop during natural infection can help guide the rational design of an HIV vaccine. Here, we described a bnAb lineage targeting the Env V2 apex and the Ab-Env co-evolution that led to development of neutralization breadth. The lineage Abs bore an anionic heavy chain complementarity-determining region 3 (CDRH3) of 25 amino acids, among the shortest known for this class of Abs, and achieved breadth with only 10% nucleotide somatic hypermutation and no insertions or deletions. The data suggested a role for Env glycoform heterogeneity in the activation of the lineage germline B cell. Finally, we showed that localized diversity at key V2 epitope residues drove bnAb maturation toward breadth, mirroring the Env evolution pattern described for another donor who developed V2-apex targeting bnAbs. Overall, these findings suggest potential strategies for vaccine approaches based on germline-targeting and serial immunogen design.

Approaches to the induction of HIV broadly neutralizing antibodies

PURPOSE OF REVIEW

A vaccine that elicits antibody responses that can neutralize the diversity of HIV clades has not yet been achieved, and is a major focus of HIV vaccine research. Here, we provide an update on the barriers to eliciting such antibodies, and how advances in immunogen design may circumvent these roadblocks, focusing on data published in the last year.

RECENT FINDINGS

Studies of how broadly neutralizing antibodies (bNAbs) develop in HIV-infected donors continue to produce key insights, suggesting that for some viral targets there are common pathways to developing breadth. Germline-targeting strategies, that aim to recruit rare precursors of bNAbs, have shown promise in immunogenicity studies, and structural biology has led to advances in immunogen design. Mapping of strain-specific tier 2 vaccine responses has highlighted the challenges that remain in driving antibodies toward breadth.

SUMMARY

Elucidation of the HIV envelope structure, together with an understanding of how bNAbs emerge in vivo has guided the design of new immunogens and vaccine strategies that show promise for eliciting protective antibodies.

Ontogeny-based immunogens for the induction of V2-directed HIV broadly neutralizing antibodies

The development of a preventative HIV vaccine able to elicit broadly neutralizing antibodies (bNAbs) remains a major challenge. Antibodies that recognize the V2 region at the apex of the HIV envelope trimer are among the most common bNAb specificities during chronic infection and many exhibit remarkable breadth and potency. Understanding the developmental pathway of these antibodies has provided insights into their precursors, and the viral strains that engage them, as well as defined how such antibodies mature to acquire breadth. V2‐apex bNAbs are derived from rare precursors with long anionic CDR H3s that are often deleted in the B cell repertoire. However, longitudinal studies suggest that once engaged, these precursors contain many of the structural elements required for neutralization, and can rapidly acquire breadth through moderate levels of somatic hypermutation in response to emerging viral variants. These commonalities in the precursors and mechanism of neutralization have enabled the identification of viral strains that show enhanced reactivity for V2 precursors from multiple donors, and may form the basis of germline targeting approaches. In parallel, new structural insights into the HIV trimer, the target of these quaternary antibodies, has created invaluable new opportunities for ontogeny‐based immunogens designed to select for rare V2‐bNAb precursors, and drive them toward breadth.

Gene-Specific Substitution Profiles Describe the Types and Frequencies of Amino Acid Changes during Antibody Somatic Hypermutation

Somatic hypermutation (SHM) plays a critical role in the maturation of antibodies, optimizing recognition initiated by recombination of V(D)J genes. Previous studies have shown that the propensity to mutate is modulated by the context of surrounding nucleotides and that SHM machinery generates biased substitutions. To investigate the intrinsic mutation frequency and substitution bias of SHMs at the amino acid level, we analyzed functional human antibody repertoires and developed mGSSP (method for gene-specific substitution profile), a method to construct amino acid substitution profiles from next-generation sequencing-determined B cell transcripts. We demonstrated that these gene-specific substitution profiles (GSSPs) are unique to each V gene and highly consistent between donors. We also showed that the GSSPs constructed from functional antibody repertoires are highly similar to those constructed from antibody sequences amplified from non-productively rearranged passenger alleles, which do not undergo functional selection. This suggests the types and frequencies, or mutational space, of a majority of amino acid changes sampled by the SHM machinery to be well captured by GSSPs. We further observed the rates of mutational exchange between some amino acids to be both asymmetric and context dependent and to correlate weakly with their biochemical properties. GSSPs provide an improved, position-dependent alternative to standard substitution matrices, and can be utilized to developing software for accurately modeling the SHM process. GSSPs can also be used for predicting the amino acid mutational space available for antigen-driven selection and for understanding factors modulating the maturation pathways of antibody lineages in a gene-specific context. The mGSSP method can be used to build, compare, and plot GSSPs¹; we report the GSSPs constructed for 69 common human V genes (DOI: 10.6084/m9.figshare.3511083) and provide high-resolution logo plots for each (DOI: 10.6084/m9.figshare.3511085).

Prospecting for an HIV vaccine

Human immunodeficiency virus (HIV) sets several challenges for the development of a preventative HIV vaccine. Predictable, protective natural immunity against HIV does not occur and so unlike most other diseases for which vaccines exist, there are few guideposts from natural infection. Nonetheless, six vaccine efficacy trials have occurred. One in particular, the Thai trial called RV144, showed partial protective efficacy and potential ways ahead to a better vaccine approach. This coupled with other lessons from studies of acute infections as well as an increasingly complex knowledge of HIV-related vaccine immunology bring hope that a vaccine solution might be reached for this pervasive and deadly pandemic.

A Phylogenetic Codon Substitution Model for Antibody Lineages

Phylogenetic methods have shown promise in understanding the development of broadly neutralizing antibody lineages (bNAbs). However, the mutational process that generates these lineages, somatic hypermutation, is biased by hotspot motifs which violates important assumptions in most phylogenetic substitution models. Here, we develop a modified GY94-type substitution model that partially accounts for this context dependency while preserving independence of sites during calculation. This model shows a substantially better fit to three well-characterized bNAb lineages than the standard GY94 model. We also demonstrate how our model can be used to test hypotheses concerning the roles of different hotspot and coldspot motifs in the evolution of B-cell lineages. Further, we explore the consequences of the idea that the number of hotspot motifs, and perhaps the mutation rate in general, is expected to decay over time in individual bNAb lineages.

Antibodyomics: bioinformatics technologies for understanding B-cell immunity to HIV-1

Numerous antibodies have been identified from HIV-1-infected donors that neutralize diverse strains of HIV-1. These antibodies may provide the basis for a B cell-mediated HIV-1 vaccine. However, it has been unclear how to elicit similar antibodies by vaccination. To address this issue, we have undertaken an informatics-based approach to understand the genetic and immunologic processes controlling the development of HIV-1-neutralizing antibodies. As DNA sequencing comprises the fastest growing database of biological information, we focused on incorporating next-generation sequencing of B-cell transcripts to determine the origin, maturation pathway, and prevalence of broadly neutralizing antibody lineages (Antibodyomics1, 2, 4, and 6). We also incorporated large-scale robotic analyses of serum neutralization to identify and quantify neutralizing antibodies in donor cohorts (Antibodyomics3). Statistical analyses furnish another layer of insight (Antibodyomics5), with physical characteristics of antibodies and their targets through molecular dynamics simulations (Antibodyomics7) and free energy perturbation analyses (Antibodyomics8) providing information-rich output. Functional interrogation of individual antibodies (Antibodyomics9) and synthetic antibody libraries (Antibodyomics10) also yields multi-dimensional data by which to understand and improve antibodies. Antibodyomics, described here, thus comprise resolution-enhancing tools, which collectively embody an information-driven discovery engine aimed toward the development of effective B cell-based vaccines.

SONAR: A High-Throughput Pipeline for Inferring Antibody Ontogenies from Longitudinal Sequencing of B Cell Transcripts

The rapid advance of massively parallel or next-generation sequencing technologies has made possible the characterization of B cell receptor repertoires in ever greater detail, and these developments have triggered a proliferation of software tools for processing and annotating these data. Of especial interest, however, is the capability to track the development of specific antibody lineages across time, which remains beyond the scope of most current programs. We have previously reported on the use of techniques such as inter- and intradonor analysis and CDR3 tracing to identify transcripts related to an antibody of interest. Here, we present Software for the Ontogenic aNalysis of Antibody Repertoires (SONAR), capable of automating both general repertoire analysis and specialized techniques for investigating specific lineages. SONAR annotates next-generation sequencing data, identifies transcripts in a lineage of interest, and tracks lineage development across multiple time points. SONAR also generates figures, such as identity-divergence plots and longitudinal phylogenetic "birthday" trees, and provides interfaces to other programs such as DNAML and BEAST. SONAR can be downloaded as a ready-to-run Docker image or manually installed on a local machine. In the latter case, it can also be configured to take advantage of a high-performance computing cluster for the most computationally intensive steps, if available. In summary, this software provides a useful new tool for the processing of large next-generation sequencing datasets and the ontogenic analysis of neutralizing antibody lineages. SONAR can be found at https://github.com/scharch/SONAR, and the Docker image can be obtained from https://hub.docker.com/r/scharch/sonar/.