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

Maturation of germinal center B cells after influenza virus vaccination in humans

Germinal centers (GC) are microanatomical lymphoid structures where affinity-matured memory B cells and long-lived bone marrow plasma cells are primarily generated. It is unclear how the maturation of B cells within the GC impacts the breadth and durability of B cell responses to influenza vaccination in humans. We used fine needle aspiration of draining lymph nodes to longitudinally track antigen-specific GC B cell responses to seasonal influenza vaccination. Antigen-specific GC B cells persisted for at least 13 wk after vaccination in two out of seven individuals. Monoclonal antibodies (mAbs) derived from persisting GC B cell clones exhibit enhanced binding affinity and breadth to influenza hemagglutinin (HA) antigens compared with related GC clonotypes isolated earlier in the response. Structural studies of early and late GC-derived mAbs from one clonal lineage in complex with H1 and H5 HAs revealed an altered binding footprint. Our study shows that inducing sustained GC reactions after influenza vaccination in humans supports the maturation of responding B cells.

Broadly neutralizing antibodies for HIV prevention: a comprehensive review and future perspectives

SUMMARYThe human immunodeficiency virus (HIV) epidemic remains a formidable global health concern, with 39 million people living with the virus and 1.3 million new infections reported in 2022. Despite anti-retroviral therapy's effectiveness in pre-exposure prophylaxis, its global adoption is limited. Broadly neutralizing antibodies (bNAbs) offer an alternative strategy for HIV prevention through passive immunization. Historically, passive immunization has been efficacious in the treatment of various diseases ranging from oncology to infectious diseases. Early clinical trials suggest bNAbs are safe, tolerable, and capable of reducing HIV RNA levels. Although challenges such as bNAb resistance have been noted in phase I trials, ongoing research aims to assess the additive or synergistic benefits of combining multiple bNAbs. Researchers are exploring bispecific and trispecific antibodies, and fragment crystallizable region modifications to augment antibody efficacy and half-life. Moreover, the potential of other antibody isotypes like IgG3 and IgA is under investigation. While promising, the application of bNAbs faces economic and logistical barriers. High manufacturing costs, particularly in resource-limited settings, and logistical challenges like cold-chain requirements pose obstacles. Preliminary studies suggest cost-effectiveness, although this is contingent on various factors like efficacy and distribution. Technological advancements and strategic partnerships may mitigate some challenges, but issues like molecular aggregation remain. The World Health Organization has provided preferred product characteristics for bNAbs, focusing on optimizing their efficacy, safety, and accessibility. The integration of bNAbs in HIV prophylaxis necessitates a multi-faceted approach, considering economic, logistical, and scientific variables. This review comprehensively covers the historical context, current advancements, and future avenues of bNAbs in HIV prevention.

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.

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.

Inference of B cell clonal families using heavy/light chain pairing information

Next generation sequencing of B cell receptor (BCR) repertoires has become a ubiquitous tool for understanding the antibody-mediated immune response: it is now common to have large volumes of sequence data coding for both the heavy and light chain subunits of the BCR. However, until the recent development of high throughput methods of preserving heavy/light chain pairing information, these samples contained no explicit information on which heavy chain sequence pairs with which light chain sequence. One of the first steps in analyzing such BCR repertoire samples is grouping sequences into clonally related families, where each stems from a single rearrangement event. Many methods of accomplishing this have been developed, however, none so far has taken full advantage of the newly-available pairing information. This information can dramatically improve clustering performance, especially for the light chain. The light chain has traditionally been challenging for clonal family inference because of its low diversity and consequent abundance of non-clonal families with indistinguishable naive rearrangements. Here we present a method of incorporating this pairing information into the clustering process in order to arrive at a more accurate partition of the data into clonally related families. We also demonstrate two methods of fixing imperfect pairing information, which may allow for simplified sample preparation and increased sequencing depth. Finally, we describe several other improvements to the partis software package.

Phylogenetic analysis of migration, differentiation, and class switching in B cells

B cells undergo rapid mutation and selection for antibody binding affinity when producing antibodies capable of neutralizing pathogens. This evolutionary process can be intermixed with migration between tissues, differentiation between cellular subsets, and switching between functional isotypes. B cell receptor (BCR) sequence data has the potential to elucidate important information about these processes. However, there is currently no robust, generalizable framework for making such inferences from BCR sequence data. To address this, we develop three parsimony-based summary statistics to characterize migration, differentiation, and isotype switching along B cell phylogenetic trees. We use simulations to demonstrate the effectiveness of this approach. We then use this framework to infer patterns of cellular differentiation and isotype switching from high throughput BCR sequence datasets obtained from patients in a study of HIV infection and a study of food allergy. These methods are implemented in the R package dowser, available at https://dowser.readthedocs.io.

B cells produce high affinity antibodies through an evolutionary process of mutation and selection during adaptive immune responses. Migration between tissues, differentiation to cellular subtypes, and switching between different antibody isotypes can be important factors in shaping the role B cells play in response to infection, autoimmune disease, and allergies. B cell receptor (BCR) sequence data has the potential to elucidate important information about these processes. However, there is currently no robust, generalizable framework for making such inferences from BCR sequence data. Here, we develop three parsimony-based summary statistics to characterize migration, differentiation, and isotype switching along B cell phylogenetic trees. We confirm the effectiveness of our approach using simulations and further apply our method to data from patients with HIV and allergy. Our methods are released in the R package dowser: https://dowser.readthedocs.io.

Ecological Processes Shaping Microbiomes of Extremely Low Birthweight Infants

The human microbiome has been implicated in affecting health outcomes in premature infants, but the ecological processes governing early life microbiome assembly remain poorly understood. Here, we investigated microbial community assembly and dynamics in extremely low birth weight infants (ELBWI) over the first 2 weeks of life. We profiled the gut, oral cavity and skin microbiomes over time using 16S rRNA gene amplicon sequencing and evaluated the ecological forces shaping these microbiomes. Though microbiomes at all three body sites were characterized by compositional instability over time and had low body-site specificity (PERMANOVA, r 2 = 0.09, p = 0.001), they could nonetheless be clustered into four discrete community states. Despite the volatility of these communities, deterministic assembly processes were detectable in this period of initial microbial colonization. To further explore these deterministic dynamics, we developed a probabilistic approach in which we modeled microbiome state transitions in each ELBWI as a Markov process, or a "memoryless" shift, from one community state to another. This analysis revealed that microbiomes from different body sites had distinctive dynamics as well as characteristic equilibrium frequencies. Time-resolved microbiome sampling of premature infants may help to refine and inform clinical practices. Additionally, this work provides an analysis framework for microbial community dynamics based on Markov modeling that can facilitate new insights, not only into neonatal microbiomes but also other human-associated or environmental microbiomes.

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.

Immunizations with diverse sarbecovirus receptor-binding domains elicit SARS-CoV-2 neutralizing antibodies against a conserved site of vulnerability

Viral mutations are an emerging concern in reducing SARS-CoV-2 vaccination efficacy. Second-generation vaccines will need to elicit neutralizing antibodies against sites that are evolutionarily conserved across the sarbecovirus subgenus. Here, we immunized mice containing a human antibody repertoire with diverse sarbecovirus receptor-binding domains (RBDs) to identify antibodies targeting conserved sites of vulnerability. Antibodies with broad reactivity against diverse clade B RBDs targeting the conserved class 4 epitope, with recurring IGHV/IGKV pairs, were readily elicited but were non-neutralizing. However, rare class 4 antibodies binding this conserved RBD supersite showed potent neutralization of SARS-CoV-2 and all variants of concern. Structural analysis revealed that the neutralizing ability of cross-reactive antibodies was reserved only for those with an elongated CDRH3 that extends the antiparallel beta-sheet RBD core and orients the antibody light chain to obstruct ACE2-RBD interactions. These results identify a structurally defined pathway for vaccine strategies eliciting escape-resistant SARS-CoV-2 neutralizing antibodies.

Functional development of a V3/glycan-specific broadly neutralizing antibody isolated from a case of HIV superinfection

Stimulating broadly neutralizing antibodies (bnAbs) directly from germline remains a barrier for HIV vaccines. HIV superinfection elicits bnAbs more frequently than single infection, providing clues of how to elicit such responses. We used longitudinal antibody sequencing and structural studies to characterize bnAb development from a superinfection case. BnAb QA013.2 bound initial and superinfecting viral Env, despite its probable naive progenitor only recognizing the superinfecting strain, suggesting both viruses influenced this lineage. A 4.15 Å cryo-EM structure of QA013.2 bound to native-like trimer showed recognition of V3 signatures (N301/N332 and GDIR). QA013.2 relies less on CDRH3 and more on framework and CDRH1 for affinity and breadth compared to other V3/glycan-specific bnAbs. Antigenic profiling revealed that viral escape was achieved by changes in the structurally-defined epitope and by mutations in V1. These results highlight shared and novel properties of QA013.2 relative to other V3/glycan-specific bnAbs in the setting of sequential, diverse antigens.

Rationalizing Random Walks: Replicating Protective Antibody Trajectories

'Reverse vaccinology 2.0' aims to rationally reproduce template antibody responses, such as broadly neutralizing antibodies against human immunodeficiency virus-1. While observations of antibody convergence across individuals support the assumption that responses may be replicated, the diversity of humoral immunity and the process of antibody selection are rooted in stochasticity. Drawing from experience with in vitro antibody engineering by directed evolution, we consider how antibody selection may be driven, as in germline-targeting vaccine approaches to elicit broadly neutralizing antibodies and illustrate the potential consequences of over-defining a template antibody response. We posit that the prospective definition of template antibody responses and the odds of replicating them must be considered within the randomness of humoral immunity.

Development of antibody-dependent cell cytotoxicity function in HIV-1 antibodies

A prerequisite for the design of an HIV vaccine that elicits protective antibodies is understanding the developmental pathways that result in desirable antibody features. The development of antibodies that mediate antibody-dependent cellular cytotoxicity (ADCC) is particularly relevant because such antibodies have been associated with HIV protection in humans. We reconstructed the developmental pathways of six human HIV-specific ADCC antibodies using longitudinal antibody sequencing data. Most of the inferred naive antibodies did not mediate detectable ADCC. Gain of antigen binding and ADCC function typically required mutations in complementarity determining regions of one or both chains. Enhancement of ADCC potency often required additional mutations in framework regions. Antigen binding affinity and ADCC activity were correlated, but affinity alone was not sufficient to predict ADCC potency. Thus, elicitation of broadly active ADCC antibodies may require mutations that enable high-affinity antigen recognition along with mutations that optimize factors contributing to functional ADCC activity.

Learning the heterogeneous hypermutation landscape of immunoglobulins from high-throughput repertoire data

Somatic hypermutations of immunoglobulin (Ig) genes occurring during affinity maturation drive B-cell receptors’ ability to evolve strong binding to their antigenic targets. The landscape of these mutations is highly heterogeneous, with certain regions of the Ig gene being preferentially targeted. However, a rigorous quantification of this bias has been difficult because of phylogenetic correlations between sequences and the interference of selective forces. Here, we present an approach that corrects for these issues, and use it to learn a model of hypermutation preferences from a recently published large IgH repertoire dataset. The obtained model predicts mutation profiles accurately and in a reproducible way, including in the previously uncharacterized Complementarity Determining Region 3, revealing that both the sequence context of the mutation and its absolute position along the gene are important. In addition, we show that hypermutations occurring concomittantly along B-cell lineages tend to co-localize, suggesting a possible mechanism for accelerating affinity maturation.

Using B cell receptor lineage structures to predict affinity

We are frequently faced with a large collection of antibodies, and want to select those with highest affinity for their cognate antigen. When developing a first-line therapeutic for a novel pathogen, for instance, we might look for such antibodies in patients that have recovered. There exist effective experimental methods of accomplishing this, such as cell sorting and baiting; however they are time consuming and expensive. Next generation sequencing of B cell receptor (BCR) repertoires offers an additional source of sequences that could be tapped if we had a reliable method of selecting those coding for the best antibodies. In this paper we introduce a method that uses evolutionary information from the family of related sequences that share a naive ancestor to predict the affinity of each resulting antibody for its antigen. When combined with information on the identity of the antigen, this method should provide a source of effective new antibodies. We also introduce a method for a related task: given an antibody of interest and its inferred ancestral lineage, which branches in the tree are likely to harbor key affinity-increasing mutations? We evaluate the performance of these methods on a wide variety of simulated samples, as well as two real data samples. These methods are implemented as part of continuing development of the partis BCR inference package, available at https://github.com/psathyrella/partis. Comments Please post comments or questions on this paper as new issues at https://git.io/Jvxkn.