Multi-Donor Longitudinal Antibody Repertoire Sequencing Reveals the Existence of Public Antibody Clonotypes in HIV-1 Infection

Novel polymorphic and copy number diversity in the antibody IGH locus of South African individuals

The heavy chain of an antibody is crucial for mediating antigen binding. IGHV genes, which partially encode the heavy chain of antibodies, exhibit vast genetic diversity largely through polymorphism and copy number variation (CNV). These genetic variations impact population-level expression levels. In this study, we analyzed expressed antibody transcriptomes and matched germline IGHV genes from donors from KwaZulu-Natal, South Africa. Amplicon NGS targeting germline IGHV sequences was performed on genomic DNA from 70 participants, eight of whom had matched datasets of expressed antibody transcriptomes. Germline IGHV sequencing identified 161 unique IGHV alleles, of which 32 were novel. A further 21 novel IGHV alleles were detected in the expressed transcriptomes of these donors. We also examined the datasets for CNV, uncovering gene duplications of 10 IGHV genes from germline sequencing and 33 genes in the expressed transcriptomes. Many of the IGHV gene duplications have not been described in other populations. This study expands our understanding of genetic differences in distinct populations and suggests the potential impact of genetic diversity on immune responses.

Deep immunoglobulin repertoire sequencing depicts a comprehensive atlas of spike-specific antibody lineages shared among COVID-19 convalescents

Neutralizing antibodies are a key component in protective humoral immunity against SARS-CoV-2. Currently, available technologies cannot track epitope-specific antibodies in global antibody repertoires. Thus, the comprehensive repertoire of spike-specific neutralizing antibodies elicited by SARS-CoV-2 infection is not fully understood. We therefore combined high-throughput immunoglobulin heavy chain (IgH) repertoire sequencing, and structural and bioinformatics analysis to establish an antibodyomics pipeline, which enables tracking spike-specific antibody lineages that target certain neutralizing epitopes. We mapped the neutralizing epitopes on the spike and determined the epitope-preferential antibody lineages. This analysis also revealed numerous overlaps between immunodominant neutralizing antibody-binding sites and mutation hotspots on spikes as observed so far in SARS-CoV-2 variants. By clustering 2677 spike-specific antibodies with 360 million IgH sequences that we sequenced, a total of 329 shared spike-specific antibody clonotypes were identified from 33 COVID-19 convalescents and 24 SARS-CoV-2-naïve individuals. Epitope mapping showed that the shared antibody responses target not only neutralizing epitopes on RBD and NTD but also non-neutralizing epitopes on S2. The immunodominance of neutralizing antibody response is determined by the occurrence of specific precursors in human naïve B-cell repertoires. We identified that only 28 out of the 329 shared spike-specific antibody clonotypes persisted for at least 12 months. Among them, long-lived IGHV3-53 antibodies are likely to evolve cross-reactivity to Omicron variants through accumulating somatic hypermutations. Altogether, we created a comprehensive atlas of spike-targeting antibody lineages in COVID-19 convalescents and antibody precursors in human naïve B cell repertoires, providing a valuable reference for future vaccine design and evaluation.

Antibody sequence determinants of viral antigen specificity

Throughout life, humans experience repeated exposure to viral antigens through infection and vaccination, resulting in the generation of diverse, antigen-specific antibody repertoires. A paramount feature of antibodies that enables their critical contributions in counteracting recurrent and novel pathogens, and consequently fostering their utility as valuable targets for therapeutic and vaccine development, is the exquisite specificity displayed against their target antigens. Yet, there is still limited understanding of the determinants of antibody-antigen specificity, particularly as a function of antibody sequence. In recent years, experimental characterization of antibody repertoires has led to novel insights into fundamental properties of antibody sequences but has been largely decoupled from at-scale antigen specificity analysis. Here, using the LIBRA-seq technology, we generated a large data set mapping antibody sequence to antigen specificity for thousands of B cells, by screening the repertoires of a set of healthy individuals against 20 viral antigens representing diverse pathogens of biomedical significance. Analysis uncovered virus-specific patterns in variable gene usage, gene pairing, somatic hypermutation, as well as the presence of convergent antiviral signatures across multiple individuals, including the presence of public antibody clonotypes. Notably, our results showed that, for B-cell receptors originating from different individuals but leveraging an identical combination of heavy and light chain variable genes, there is a specific CDRH3 identity threshold above which B cells appear to exclusively share the same antigen specificity. This finding provides a quantifiable measure of the relationship between antibody sequence and antigen specificity and further defines experimentally grounded criteria for defining public antibody clonality.IMPORTANCEThe B-cell compartment of the humoral immune system plays a critical role in the generation of antibodies upon new and repeated pathogen exposure. This study provides an unprecedented level of detail on the molecular characteristics of antibody repertoires that are specific to each of the different target pathogens studied here and provides empirical evidence in support of a 70% CDRH3 amino acid identity threshold in pairs of B cells encoded by identical IGHV:IGL(K)V genes, as a means of defining public clonality and therefore predicting B-cell antigen specificity in different individuals. This is of exceptional importance when leveraging public clonality as a method to annotate B-cell receptor data otherwise lacking antigen specificity information. Understanding the fundamental rules of antibody-antigen interactions can lead to transformative new approaches for the development of antibody therapeutics and vaccines against current and emerging viruses.

Ultrapotent class I neutralizing antibodies post Omicron breakthrough infection overcome broad SARS-CoV-2 escape variants

BACKGROUND

The spread of emerging SARS-CoV-2 immune escape sublineages, especially JN.1 and KP.2, has resulted in new waves of COVID-19 globally. The evolving memory B cell responses elicited by the parental Omicron variants to subvariants with substantial antigenic drift remain incompletely investigated.

METHODS

Using the single B cell antibody cloning technology, we isolated single memory B cells, delineated the B cell receptor repertoire and conducted the pseudovirus-based assay for recovered neutralizing antibodies (NAb) screening. We analyzed the cryo-EM structures of top broadly NAbs (bnAbs) and evaluated their in vivo efficacy (golden Syrian hamster model).

FINDINGS

By investigating the evolution of human B cell immunity, we discovered a new panel of bnAbs arising from vaccinees after Omicron BA.2/BA.5 breakthrough infections. Two lead bnAbs neutralized major Omicron subvariants including JN.1 and KP.2 with IC50 values less than 10 ng/mL, representing ultrapotent receptor binding domain (RBD)-specific class I bnAbs. They belonged to the IGHV3-53/3-66 clonotypes instead of evolving from the pre-existing vaccine-induced IGHV1-58/IGKV3-20 bnAb ZCB11. Despite sequence diversity, they targeted previously unrecognized, highly conserved conformational epitopes in the receptor binding motif (RBM) for ultrapotent ACE2 blockade. The lead bnAb ZCP3B4 not only protected the lungs of hamsters intranasally challenged with BA.5.2, BQ.1.1 and XBB.1.5 but also prevented their contact transmission.

INTERPRETATION

Our findings demonstrated that class I bnAbs have evolved an ultrapotent mode of action protecting against highly transmissible and broad Omicron escape variants, and their epitopes are potential targets for novel bnAbs and vaccine development.

FUNDING

A full list of funding bodies that contributed to this study can be found in the Acknowledgements section.

Rapid affinity optimization of an anti-TREM2 clinical lead antibody by cross-lineage immune repertoire mining

We describe a process for rapid antibody affinity optimization by repertoire mining to identify clones across B cell clonal lineages based on convergent immune responses where antigen-specific clones with the same heavy (VH) and light chain germline segment pairs, or parallel lineages, bind a single epitope on the antigen. We use this convergence framework to mine unique and distinct VH lineages from rat anti-triggering receptor on myeloid cells 2 (TREM2) antibody repertoire datasets with high diversity in the third complementarity-determining loop region (CDR H3) to further affinity-optimize a high-affinity agonistic anti-TREM2 antibody while retaining critical functional properties. Structural analyses confirm a nearly identical binding mode of anti-TREM2 variants with subtle but significant structural differences in the binding interface. Parallel lineage repertoire mining is uniquely tailored to rationally explore the large CDR H3 sequence space in antibody repertoires and can be easily and generally applied to antibodies discovered in vivo.

High throughput long-read sequencing of circulating lymphocytes of the evolutionarily distant sea lamprey reveals diversity and common elements of the variable lymphocyte receptor B (VLRB) repertoire

The leucine-rich repeat-based variable lymphocyte receptor B (VLRB) antibody system of jawless vertebrates is capable of generating an antibody repertoire equal to or exceeding the diversity of antibody repertoires of jawed vertebrates. Unlike immunoglobulin-based immune repertoires, the VLRB repertoire diversity is characterized by variable lengths of VLRB encoding transcripts, rendering conventional immunoreceptor repertoire sequencing approaches unsuitable for VLRB repertoire sequencing. Here we demonstrate that long-read single-molecule real-time (SMRT) sequencing (PacBio) approaches permit the efficient large-scale assessment of the VLRB repertoire. We present a computational pipeline for sequence data processing and provide the first repertoire-based analysis of VLRB protein characteristics including properties of its subunits and regions of diversity within each structural leucine-rich repeat subunit. Our study provides a template to explore changes in the VLRB repertoire during immune responses and to establish large scale VLRB repertoire databases for computational approaches aimed at isolating monoclonal VLRB reagents for biomedical research and clinical applications.

Potent HPIV3-neutralizing IGHV5-51 Antibodies Identified from Multiple Individuals Show L Chain and CDRH3 Promiscuity

Human parainfluenza virus 3 (HPIV3) is a widespread pathogen causing severe and lethal respiratory illness in at-risk populations. Effective countermeasures are in various stages of development; however, licensed therapeutic and prophylactic options are not available. The fusion glycoprotein (HPIV3 F), responsible for facilitating viral entry into host cells, is a major target of neutralizing Abs that inhibit infection. Although several neutralizing Abs against a small number of HPIV3 F epitopes have been identified to date, relatively little is known about the Ab response to HPIV3 compared with other pathogens, such as influenza virus and SARS-CoV-2. In this study, we aimed to characterize a set of HPIV3-specific Abs identified in multiple individuals for genetic signatures, epitope specificity, neutralization potential, and publicness. We identified 12 potently neutralizing Abs targeting three nonoverlapping epitopes on HPIV3 F. Among these, six Abs identified from two different individuals used Ig heavy variable gene IGHV 5-51, with five of the six Abs targeting the same epitope. However, despite the use of the same H chain variable (VH) gene, these Abs used multiple different L chain variable genes (VL) and diverse H chain CDR 3 (CDRH3) sequences. Together, these results provide further information about the genetic and functional characteristics of HPIV3-neutralizing Abs and suggest the existence of a reproducible VH-dependent Ab response associated with VL and CDRH3 promiscuity. Understanding sites of HPIV3 F vulnerability and the genetic and molecular characteristics of Abs targeting these sites will help guide efforts for effective vaccine and therapeutic development.

Evidence of antigenic drift in the fusion machinery core of SARS-CoV-2 spike

Antigenic drift of SARS-CoV-2 is typically defined by mutations in the N-terminal domain and receptor binding domain of spike protein. In contrast, whether antigenic drift occurs in the S2 domain remains largely elusive. Here, we perform a deep mutational scanning experiment to identify S2 mutations that affect binding of SARS-CoV-2 spike to three S2 apex public antibodies. Our results indicate that spatially diverse mutations, including D950N and Q954H, which are observed in Delta and Omicron variants, respectively, weaken the binding of spike to these antibodies. Although S2 apex antibodies are known to be nonneutralizing, we show that they confer protection in vivo through Fc-mediated effector functions. Overall, this study indicates that the S2 domain of SARS-CoV-2 spike can undergo antigenic drift, which represents a potential challenge for the development of more universal coronavirus vaccines.

Systematic evaluation of B-cell clonal family inference approaches

The reconstruction of clonal families (CFs) in B-cell receptor (BCR) repertoire analysis is a crucial step to understand the adaptive immune system and how it responds to antigens. The BCR repertoire of an individual is formed throughout life and is diverse due to several factors such as gene recombination and somatic hypermutation. The use of Adaptive Immune Receptor Repertoire sequencing (AIRR-seq) using next generation sequencing enabled the generation of full BCR repertoires that also include rare CFs. The reconstruction of CFs from AIRR-seq data is challenging and several approaches have been developed to solve this problem. Currently, most methods use the heavy chain (HC) only, as it is more variable than the light chain (LC). CF reconstruction options include the definition of appropriate sequence similarity measures, the use of shared mutations among sequences, and the possibility of reconstruction without preliminary clustering based on V- and J-gene annotation. In this study, we aimed to systematically evaluate different approaches for CF reconstruction and to determine their impact on various outcome measures such as the number of CFs derived, the size of the CFs, and the accuracy of the reconstruction. The methods were compared to each other and to a method that groups sequences based on identical junction sequences and another method that only determines subclones. We found that after accounting for data set variability, in particular sequencing depth and mutation load, the reconstruction approach has an impact on part of the outcome measures, including the number of CFs. Simulations indicate that unique junctions and subclones should not be used as substitutes for CF and that more complex methods do not outperform simpler methods. Also, we conclude that different approaches differ in their ability to correctly reconstruct CFs when not considering the LC and to identify shared CFs. The results showed the effect of different approaches on the reconstruction of CFs and highlighted the importance of choosing an appropriate method.

Computational Analysis of B-Cell Receptor (BCR) Immune Repertoires with Abalign

The widespread application of high-throughput sequencing technology has generated massive sequences of B-cell receptor (BCR) immune repertoires. Computational analysis of these data has gained significant attention due to the increasing importance of immunotherapy and precision medicine. It not only reveals the diversity and dynamic changes in immune responses, contributing to the study of associated diseases, but also provides valuable information for immunodiagnostics and drug development. Recently, we introduced a BCR-specific multiple sequence alignment (MSA) method along with a comprehensive platform software called Abalign, which stands out as an excellent choice for analyzing BCR immune repertoires due to its unique high-throughput processing capability. It offers ultra-fast MSA functionality and a wide range of analytical features, including BCR/antibody extraction, clonal grouping, lineage tree construction, mutation profiling, diversity statistics, VJ gene assignment, antibody humanization, and more. Importantly, users can perform these analyses using the graphical user interface without any programming skills or scripts. In this article, we present a series of protocols that integrate Abalign's analysis modules into a cohesive workflow. This step-by-step workflow provides detailed instructions for software installation, data preparation, and comprehensive analysis of BCR immune repertoires. This workflow facilitates the efficient acquisition of comprehensive results in profiling BCR immune repertoires, offering insights into the impacts of infectious diseases, allergies, autoimmune disorders, tumor immunology, and antibody drugs. Abalign is freely available at http://cao.labshare.cn/abalign/. © 2024 Wiley Periodicals LLC. Basic Protocol 1: Resource preparation Basic Protocol 2: Analyzing BCR immune repertoires Support Protocol 1: Aiding antibody humanization Support Protocol 2: Constructing B-cell lineage trees Alternate Protocol: Running with Linux command line Basic Protocol 3: Comparing BCR immune repertoires.

SARS-CoV-2 antibodies from children exhibit broad neutralization and belong to adult public clonotypes

From the beginning of the COVID-19 pandemic, children have exhibited different susceptibility to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, reinfection, and disease compared with adults. Motivated by the established significance of SARS-CoV-2-neutralizing antibodies in adults, here we characterize SARS-CoV-2-specific antibody repertoires in a young cohort of individuals aged from 5 months to 18 years old. Our results show that neutralizing antibodies in children possess similar genetic features compared to antibodies identified in adults, with multiple antibodies from children belonging to previously established public antibody clonotypes in adults. Notably, antibodies from children show potent neutralization of circulating SARS-CoV-2 variants that have cumulatively resulted in resistance to virtually all approved monoclonal antibody therapeutics. Our results show that children can rely on similar SARS-CoV-2 antibody neutralization mechanisms compared to adults and are an underutilized source for the discovery of effective antibody therapeutics to counteract the ever-evolving pandemic.

Recognition determinants of improved HIV-1 neutralization by a heavy chain matured pediatric antibody

The structural and characteristic features of HIV-1 broadly neutralizing antibodies (bnAbs) from chronically infected pediatric donors are currently unknown. Herein, we characterized a heavy chain matured HIV-1 bnAb 44m, identified from a pediatric elite-neutralizer. Interestingly, in comparison to its wild-type AIIMS-P01 bnAb, 44m exhibited moderately higher level of somatic hypermutations of 15.2%. The 44m neutralized 79% of HIV-1 heterologous viruses (n = 58) tested, with a geometric mean IC50 titer of 0.36 μg/mL. The cryo-EM structure of 44m Fab in complex with fully cleaved glycosylated native-like BG505.SOSIP.664.T332N gp140 envelope trimer at 4.4 Å resolution revealed that 44m targets the V3-glycan N332-supersite and GDIR motif to neutralize HIV-1 with improved potency and breadth, plausibly attributed by a matured heavy chain as compared to that of wild-type AIIMS-P01. This study further improves our understanding on pediatric HIV-1 bnAbs and structural basis of broad HIV-1 neutralization by 44m may be useful blueprint for vaccine design in future.

Molecular basis of SARS-CoV-2 Omicron variant evasion from shared neutralizing antibody response

Understanding the molecular features of neutralizing epitopes is important for developing vaccines/therapeutics against emerging SARS-CoV-2 variants. We describe three monoclonal antibodies (mAbs) generated from COVID-19 recovered individuals during the first wave of the pandemic in India. These mAbs had publicly shared near germline gene usage and potently neutralized Alpha and Delta, poorly neutralized Beta, and failed to neutralize Omicron BA.1 SARS-CoV-2 variants. Structural analysis of these mAbs in complex with trimeric spike protein showed that all three mAbs bivalently bind spike with two mAbs targeting class 1 and one targeting a class 4 receptor binding domain epitope. The immunogenetic makeup, structure, and function of these mAbs revealed specific molecular interactions associated with the potent multi-variant binding/neutralization efficacy. This knowledge shows how mutational combinations can affect the binding or neutralization of an antibody, which in turn relates to the efficacy of immune responses to emerging SARS-CoV-2 escape variants.

Mixed Origins: HIV gp120-Specific Memory Develops from Pre-Existing Memory and Naive B Cells Following Vaccination in Humans

The most potent and broad HIV envelope (Env)-specific antibodies often when reverted to their inferred germline versions representing the naive B cell receptor, fail to bind Env, suggesting that the initial responding B cell population not only exclusively comprises a naive population, but also a pre-existing cross-reactive antigen-experienced B cell pool that expands following Env exposure. Previously we isolated gp120-reactive monoclonal antibodies (mAbs) from participants in HVTN 105, an HIV vaccine trial. Using deep sequencing, focused on immunoglobulin G (IgG), IgA, and IgM, VH-lineage tracking, we identified four of these mAb lineages in pre-immune peripheral blood. We also looked through the ∼7 month postvaccination bone marrow, and interestingly, several of these lineages that were found in prevaccination blood were still persistent in the postvaccination bone marrow, including the CD138+ long-lived plasma cell compartment. The majority of the pre-immune lineage members included IgM, however, IgG and IgA members were also prevalent and exhibited somatic hypermutation. These results suggest that vaccine-induced gp120-specific antibody lineages originate from both naive and cross-reactive memory B cells. ClinicalTrials.gov NCT02207920.

Identification of CDRH3 loops in the B cell receptor repertoire that can be engaged by candidate immunogens

A major goal for the development of vaccines against rapidly mutating viruses, such as influenza or HIV, is to elicit antibodies with broad neutralization capacity. However, B cell precursors capable of maturing into broadly neutralizing antibodies (bnAbs) can be rare in the immune repertoire. Due to the stochastic nature of B cell receptor (BCR) rearrangement, a limited number of third heavy chain complementary determining region (CDRH3) sequences are identical between different individuals. Thus, in order to successfully engage broadly neutralizing antibody precursors that rely on their CDRH3 loop for antigen recognition, immunogens must be able to tolerate sequence diversity in the B cell receptor repertoire across an entire vaccinated population. Here, we present a combined experimental and computational approach to identify BCRs in the human repertoire with CDRH3 loops predicted to be engaged by a target immunogen. For a given antibody/antigen pair, deep mutational scanning was first used to measure the effect of CDRH3 loop substitution on binding. BCR sequences, isolated experimentally or generated in silico, were subsequently evaluated to identify CDRH3 loops expected to be bound by the candidate immunogen. We applied this method to characterize two HIV-1 germline-targeting immunogens and found differences in the frequencies with which they are expected to engage target B cells, thus illustrating how this approach can be used to evaluate candidate immunogens towards B cell precursors engagement and to inform immunogen optimization strategies for more effective vaccine design.

Systematic analysis of human antibody response to ebolavirus glycoprotein shows high prevalence of neutralizing public clonotypes

Understanding the human antibody response to emerging viral pathogens is key to epidemic preparedness. As the size of the B cell response to a pathogenic-virus-protective antigen is poorly defined, we perform deep paired heavy- and light-chain sequencing in Ebola virus glycoprotein (EBOV-GP)-specific memory B cells, allowing analysis of the ebolavirus-specific antibody repertoire both genetically and functionally. This approach facilitates investigation of the molecular and genetic basis for the evolution of cross-reactive antibodies by elucidating germline-encoded properties of antibodies to EBOV and identification of the overlap between antibodies in the memory B cell and serum repertoire. We identify 73 public clonotypes of EBOV, 20% of which encode antibodies with neutralization activity and capacity to protect mice in vivo. This comprehensive analysis of the public and private antibody repertoire provides insight into the molecular basis of the humoral immune response to EBOV GP, which informs the design of vaccines and improved therapeutics.

Modeling and predicting the overlap of B- and T-cell receptor repertoires in healthy and SARS-CoV-2 infected individuals

Adaptive immunity's success relies on the extraordinary diversity of protein receptors on B and T cell membranes. Despite this diversity, the existence of public receptors shared by many individuals gives hope for developing population-wide vaccines and therapeutics. Using probabilistic modeling, we show many of these public receptors are shared by chance in healthy individuals. This predictable overlap is driven not only by biases in the random generation process of receptors, as previously reported, but also by their common functional selection. However, the model underestimates sharing between repertoires of individuals infected with SARS-CoV-2, suggesting strong specific antigen-driven convergent selection. We exploit this discrepancy to identify COVID-associated receptors, which we validate against datasets of receptors with known viral specificity. We study their properties in terms of sequence features and network organization, and use them to design an accurate diagnostic tool for predicting SARS-CoV-2 status from repertoire data.

The evolutionary and functional significance of germline immunoglobulin gene variation

The recombination between immunoglobulin (IG) gene segments determines an individual's naïve antibody repertoire and, consequently, (auto)antigen recognition. Emerging evidence suggests that mammalian IG germline variation impacts humoral immune responses associated with vaccination, infection, and autoimmunity - from the molecular level of epitope specificity, up to profound changes in the architecture of antibody repertoires. These links between IG germline variants and immunophenotype raise the question on the evolutionary causes and consequences of diversity within IG loci. We discuss why the extreme diversity in IG loci remains a mystery, why resolving this is important for the design of more effective vaccines and therapeutics, and how recent evidence from multiple lines of inquiry may help us do so.

Comparative global B cell receptor repertoire difference induced by SARS-CoV-2 infection or vaccination via single-cell V(D)J sequencing

Dynamic changes of the paired heavy and light chain B cell receptor (BCR) repertoire provide an essential insight into understanding the humoral immune response post-SARS-CoV-2 infection and vaccination. However, differences between the endogenous paired BCR repertoire kinetics in SARS-CoV-2 infection and previously recovered/naïve subjects treated with the inactivated vaccine remain largely unknown. We performed single-cell V(D)J sequencing of B cells from six healthy donors with three shots of inactivated SARS-CoV-2 vaccine (BBIBP-CorV), five people who received the BBIBP-CorV vaccine after having recovered from COVID-19, five unvaccinated COVID-19 recovered patients and then integrated with public data of B cells from four SARS-CoV-2-infected subjects. We discovered that BCR variable (V) genes were more prominently used in the SARS-CoV-2 exposed groups (both in the group with active infection and in the group that had recovered) than in the vaccinated groups. The VH gene that expanded the most after SARS-CoV-2 infection was IGHV3-33, while IGHV3-23 in the vaccinated groups. SARS-CoV-2-infected group enhanced more BCR clonal expansion and somatic hypermutation than the vaccinated healthy group. A small proportion of public clonotypes were shared between the SARS-CoV-2 infected, vaccinated healthy, and recovered groups. Moreover, several public antibodies had been identified against SARS-CoV-2 spike protein. We comprehensively characterize the paired heavy and light chain BCR repertoire from SARS-CoV-2 infection to vaccination, providing further guidance for the development of the next-generation precision vaccine.

Functional antibodies exhibit light chain coherence

The vertebrate adaptive immune system modifies the genome of individual B cells to encode antibodies that bind particular antigens1. In most mammals, antibodies are composed of heavy and light chains that are generated sequentially by recombination of V, D (for heavy chains), J and C gene segments. Each chain contains three complementarity-determining regions (CDR1-CDR3), which contribute to antigen specificity. Certain heavy and light chains are preferred for particular antigens2-22. Here we consider pairs of B cells that share the same heavy chain V gene and CDRH3 amino acid sequence and were isolated from different donors, also known as public clonotypes23,24. We show that for naive antibodies (those not yet adapted to antigens), the probability that they use the same light chain V gene is around 10%, whereas for memory (functional) antibodies, it is around 80%, even if only one cell per clonotype is used. This property of functional antibodies is a phenomenon that we call light chain coherence. We also observe this phenomenon when similar heavy chains recur within a donor. Thus, although naive antibodies seem to recur by chance, the recurrence of functional antibodies reveals surprising constraint and determinism in the processes of V(D)J recombination and immune selection. For most functional antibodies, the heavy chain determines the light chain.

Molecular basis of SARS-CoV-2 Omicron variant evasion from shared neutralizing antibody response

A detailed understanding of the molecular features of the neutralizing epitopes developed by viral escape mutants is important for predicting and developing vaccines or therapeutic antibodies against continuously emerging SARS-CoV-2 variants. Here, we report three human monoclonal antibodies (mAbs) generated from COVID-19 recovered individuals during first wave of pandemic in India. These mAbs had publicly shared near germline gene usage and potently neutralized Alpha and Delta, but poorly neutralized Beta and completely failed to neutralize Omicron BA.1 SARS-CoV-2 variants. Structural analysis of these three mAbs in complex with trimeric spike protein showed that all three mAbs are involved in bivalent spike binding with two mAbs targeting class-1 and one targeting class-4 Receptor Binding Domain (RBD) epitope. Comparison of immunogenetic makeup, structure, and function of these three mAbs with our recently reported class-3 RBD binding mAb that potently neutralized all SARS-CoV-2 variants revealed precise antibody footprint, specific molecular interactions associated with the most potent multi-variant binding / neutralization efficacy. This knowledge has timely significance for understanding how a combination of certain mutations affect the binding or neutralization of an antibody and thus have implications for predicting structural features of emerging SARS-CoV-2 escape variants and to develop vaccines or therapeutic antibodies against these.

Molecular basis of SARS-CoV-2 Omicron variant evasion from shared neutralizing antibody response

A detailed understanding of the molecular features of the neutralizing epitopes developed by viral escape mutants is important for predicting and developing vaccines or therapeutic antibodies against continuously emerging SARS-CoV-2 variants. Here, we report three human monoclonal antibodies (mAbs) generated from COVID-19 recovered individuals during first wave of pandemic in India. These mAbs had publicly shared near germline gene usage and potently neutralized Alpha and Delta, but poorly neutralized Beta and completely failed to neutralize Omicron BA.1 SARS-CoV-2 variants. Structural analysis of these three mAbs in complex with trimeric spike protein showed that all three mAbs are involved in bivalent spike binding with two mAbs targeting class-1 and one targeting class-4 Receptor Binding Domain (RBD) epitope. Comparison of immunogenetic makeup, structure, and function of these three mAbs with our recently reported class-3 RBD binding mAb that potently neutralized all SARS-CoV-2 variants revealed precise antibody footprint, specific molecular interactions associated with the most potent multi-variant binding / neutralization efficacy. This knowledge has timely significance for understanding how a combination of certain mutations affect the binding or neutralization of an antibody and thus have implications for predicting structural features of emerging SARS-CoV-2 escape variants and to develop vaccines or therapeutic antibodies against these.

Advances in antibody discovery from human BCR repertoires

Antibodies make up an important and growing class of compounds used for the diagnosis or treatment of disease. While traditional antibody discovery utilized immunization of animals to generate lead compounds, technological innovations have made it possible to search for antibodies targeting a given antigen within the repertoires of B cells in humans. Here we group these innovations into four broad categories: cell sorting allows the collection of cells enriched in specificity to one or more antigens; BCR sequencing can be performed on bulk mRNA, genomic DNA or on paired (heavy-light) mRNA; BCR repertoire analysis generally involves clustering BCRs into specificity groups or more in-depth modeling of antibody-antigen interactions, such as antibody-specific epitope predictions; validation of antibody-antigen interactions requires expression of antibodies, followed by antigen binding assays or epitope mapping. Together with innovations in Deep learning these technologies will contribute to the future discovery of diagnostic and therapeutic antibodies directly from humans.

AIRRscape: An interactive tool for exploring B-cell receptor repertoires and antibody responses

The sequencing of antibody repertoires of B-cells at increasing coverage and depth has led to the identification of vast numbers of immunoglobulin heavy and light chains. However, the size and complexity of these Adaptive Immune Receptor Repertoire sequencing (AIRR-seq) datasets makes it difficult to perform exploratory analyses. To aid in data exploration, we have developed AIRRscape, an R Shiny-based interactive web browser application that enables B-cell receptor (BCR) and antibody feature discovery through comparisons among multiple repertoires. Using AIRR-seq data as input, AIRRscape starts by aggregating and sorting repertoires into interactive and explorable bins of germline V-gene, germline J-gene, and CDR3 length, providing a high-level view of the entire repertoire. Interesting subsets of repertoires can be quickly identified and selected, and then network topologies of CDR3 motifs can be generated for further exploration. Here we demonstrate AIRRscape using patient BCR repertoires and sequences of published monoclonal antibodies to investigate patterns of humoral immunity to three viral pathogens: SARS-CoV-2, HIV-1, and DENV (dengue virus). AIRRscape reveals convergent antibody sequences among datasets for all three pathogens, although HIV-1 antibody datasets display limited convergence and idiosyncratic responses. We have made AIRRscape available as a web-based Shiny application, along with code on GitHub to encourage its open development and use by immuno-informaticians, virologists, immunologists, vaccine developers, and other scientists that are interested in exploring and comparing multiple immune receptor repertoires.

Association of B cell profile and receptor repertoire with the progression of Alzheimer's disease

Alzheimer's disease (AD) is the most prevalent type of dementia. Reports have revealed that the peripheral immune system is linked to neuropathology; however, little is known about the contribution of B lymphocytes in AD. For this longitudinal study, 133 participants are included at baseline and second-year follow-up. Also, we analyze B cell receptor (BCR) repertoire data generated from a public dataset of three normal and 10 AD samples and perform BCR repertoire profiling and pairwise sharing analysis. As a result, longitudinal increase in B lymphocytes is associated with increased cerebral amyloid deposition and hyperactivates induced pluripotent stem cell-derived microglia with loss-of-function for beta-amyloid clearance. Patients with AD share similar class-switched BCR sequences with identical isotypes, despite the high somatic hypermutation rate. Thus, BCR repertoire profiling can lead to the development of individualized immune-based therapeutics and treatment. We provide evidence of both quantitative and qualitative changes in B lymphocytes during AD pathogenesis.

Optimization of peripheral blood volume for in silico reconstitution of the human B cell receptor repertoire

B cells recognize antigens via membrane‐expressed B‐cell receptors (BCR) and antibodies. Similar human BCR sequences are frequently found at a significantly higher frequency than that theoretically calculated. Patients infected with SARS‐CoV2 and HIV or with autoimmune diseases share very similar BCRs. Therefore, in silico reconstitution of BCR repertoires and identification of stereotypical BCR sequences related to human pathology have diagnostic potential. Furthermore, monitoring changes of clinically significant BCR sequences and isotype conversion has prognostic potential. For BCR repertoire analysis, peripheral blood (PB) is the most convenient source. However, the optimal human PB volume for in silico reconstitution of the BCR repertoire has not been studied in detail. Here, we sampled 5, 10, and 20 mL PB from the left arm and 40 mL PB from the right arm of two volunteers, reconstituted in silico PB BCR repertoires, and compared their composition. In both volunteers, PB sampling over 20 mL resulted in slight increases in functional unique sequences (FUSs) or almost no increase in repertoire diversity. All FUSs with a frequency above 0.08% or 0.03% in the 40 mL PB BCR repertoire were detected even in the 5 mL PB BCR repertoire from each volunteer. FUSs with a higher frequency were more likely to be found in BCR repertoires from reduced PB volume, and those coexisting in two repertoires showed a statistically significant correlation in frequency irrespective of sampled anatomical site. The correlation was more significant in higher‐frequency FUSs. These observations support the potential of BCR repertoire analysis for diagnosis.

A large-scale systematic survey reveals recurring molecular features of public antibody responses to SARS-CoV-2

Global research to combat the COVID-19 pandemic has led to the isolation and characterization of thousands of human antibodies to the SARS-CoV-2 spike protein, providing an unprecedented opportunity to study the antibody response to a single antigen. Using the information derived from 88 research publications and 13 patents, we assembled a dataset of ∼8,000 human antibodies to the SARS-CoV-2 spike protein from >200 donors. By analyzing immunoglobulin V and D gene usages, complementarity-determining region H3 sequences, and somatic hypermutations, we demonstrated that the common (public) responses to different domains of the spike protein were quite different. We further used these sequences to train a deep-learning model to accurately distinguish between the human antibodies to SARS-CoV-2 spike protein and those to influenza hemagglutinin protein. Overall, this study provides an informative resource for antibody research and enhances our molecular understanding of public antibody responses.

Clonal expansion and markers of directed mutation of IGHV4-34 B cells in plasmablasts during Kawasaki disease

Kawasaki disease (KD) is the leading cause of acquired heart disease in children. The cause remains unknown; however, epidemiologic and demographic data support a single preceding infectious agent may lead to KD. A variety of pathophysiologic responses have been proposed, including direct invasion of the coronary arteries, a superantigen response, and a post-infectious autoimmune phenomenon. A role for B cell responses during KD are supported by numerous findings including B cell specific markers identified in genome wide association studies. We have recently published data showing children with KD have similar plasmablast (PB) responses to children with infections. Since during other infections, cells expressing antibodies against the preceding infection are enriched in PBs, we sought to explore the specific antibodies encoded by PBs during KD. In one child we see a massive expansion in IGHV4-34 utilizing antibodies, which has been associated with autoimmunity in the past. We further explored this expansion of IGHV4-34 utilization during the peripheral PB rise with next generation sequencing (NGS) analysis and utilizing newer techniques of chromium chip single cell separation (10x Genomics®). We also utilized peptide array screening to attempt to identify an antigen to the most prolific clones.

Characterization of human IgM and IgG repertoires in individuals with chronic HIV-1 infection

Advancements in high-throughput sequencing (HTS) of antibody repertoires (Ig-Seq) have unprecedentedly improved our ability to characterize the antibody repertoires on a large scale. However, currently, only a few studies explored the influence of chronic HIV-1 infection on human antibody repertoires and many of them reached contradictory conclusions, possibly limited by inadequate sequencing depth and throughput. To better understand how HIV-1 infection would impact humoral immune system, in this study, we systematically analyzed the differences between the IgM (HIV-IgM) and IgG (HIV-IgG) heavy chain repertoires of HIV-1 infected patients, as well as between antibody repertoires of HIV-1 patients and healthy donors (HH). Notably, the public unique clones accounted for only a negligible proportion between the HIV-IgM and HIV-IgG repertoires libraries, and the diversity of unique clones in HIV-IgG remarkably reduced. In aspect of somatic mutation rates of CDR1 and CDR2, the HIV-IgG repertoire was higher than HIV-IgM. Besides, the average length of CDR3 region in HIV-IgM was significant longer than that in the HH repertoire, presumably caused by the great number of novel VDJ rearrangement patterns, especially a massive use of IGHJ6. Moreover, some of the B cell clonotypes had numerous clones, and somatic variants were detected within the clonotype lineage in HIV-IgG, indicating HIV-1 neutralizing activities. The in-depth characterization of HIV-IgG and HIV-IgM repertoires enriches our knowledge in the profound effect of HIV-1 infection on human antibody repertoires and may have practical value for the discovery of therapeutic antibodies.

•

Ultra-deep sequencing of both IgM and IgG repertoires in chronic HIV-1 infection.

•

VDJ gene rearrangement patterns can be dramatically changed by HIV-1 infection.

•

Multiple mechanisms cause the high complexity of HIV-1-experienced antibodies.

•

Discovery of promising neutralizing HIV-1 antibodies from antibody repertoires.

Sequence and functional characterization of a public HIV-specific antibody clonotype

Public antibody clonotypes shared among multiple individuals have been identified for several pathogens. However, little is known about the determinants of antibody "publicness". Here, we characterize the sequence and functional properties of antibodies from a public clonotype targeting the CD4 binding site on HIV-1 Env. Our results showed that HIV-1 specificity for the public antibodies studied here, comprising sequences from three individuals, was modulated by the VH, but not VL, germline gene. Non-native pairing of public heavy and light chains from different individuals suggested functional complementation of sequences within this public antibody clonotype. The strength of antigen recognition appeared to be dependent on the specific antibody light chain used, but not on other sequence features such as native-antibody or germline sequence identity. Understanding the determinants of antibody clonotype "publicness" can provide insights into the fundamental rules of host-pathogen interactions at the population level, with implications for clonotype-specific vaccine development.

Epitope profiling using computational structural modelling demonstrated on coronavirus-binding antibodies

Identifying the epitope of an antibody is a key step in understanding its function and its potential as a therapeutic. Sequence-based clonal clustering can identify antibodies with similar epitope complementarity, however, antibodies from markedly different lineages but with similar structures can engage the same epitope. We describe a novel computational method for epitope profiling based on structural modelling and clustering. Using the method, we demonstrate that sequence dissimilar but functionally similar antibodies can be found across the Coronavirus Antibody Database, with high accuracy (92% of antibodies in multiple-occupancy structural clusters bind to consistent domains). Our approach functionally links antibodies with distinct genetic lineages, species origins, and coronavirus specificities. This indicates greater convergence exists in the immune responses to coronaviruses than is suggested by sequence-based approaches. Our results show that applying structural analytics to large class-specific antibody databases will enable high confidence structure-function relationships to be drawn, yielding new opportunities to identify functional convergence hitherto missed by sequence-only analysis.

Antibodies are a key component of the immune system that combat pathogens by binding to a defined region of their molecular surface (known as an ‘epitope’). The ability to map which antibodies target the same epitopes is crucial when designing non-competing antibody therapeutics or predicting the influence of pathogen mutation on population immunity. While one can use laboratory experiments to deduce when pairs of antibodies engage the same epitope, such experiments are very expensive and time consuming if used to compare on the order of thousands of antibodies. In this work, we report a new computational algorithm (SPACE) that clusters antibodies that target the same epitope based on their predicted 3D structure, as binding site structure is a property often conserved between binders complementary to the same epitope. Unlike existing antibody epitope profiling tools which assume two antibodies must share a high sequence identity/similar genetic basis to engage the same region, our orthogonal method can detect broader patterns of convergent evolution across binders to different pathogen strains, and between antibodies with different genetic and even species origins.

Analysis of B cell receptor repertoires reveals key signatures of systemic B cell response after SARS-CoV-2 infection

A comprehensive study of the B cell response against SARS-CoV-2 could be significant for understanding the immune response and developing therapeutical antibodies and vaccines. To define the dynamics and characteristics of the antibody repertoire following SARS-CoV-2 infection, we analyzed the mRNA transcripts of immunoglobulin heavy chain (IgH) repertoires of 24 peripheral blood samples collected between 3 and 111 days after symptom onset from 10 COVID-19 patients. Massive clonal expansion of naive B cells with limited somatic hypermutation (SHM) was observed in the second week after symptom onset. The proportion of low-SHM IgG clones strongly correlated with spike-specific IgG antibody titers, highlighting the significant activation of naive B cells in response to a novel virus infection. The antibody isotype switching landscape showed a transient IgA surge in the first week after symptom onset, followed by a sustained IgG elevation that lasted for at least 3 months. SARS-CoV-2 infection elicited poly-germ line reactive antibody responses. Interestingly, 17 different IGHV germ line genes recombined with IGHJ6 showed significant clonal expansion. By comparing the IgH repertoires that we sequenced with the 774 reported SARS-CoV-2–reactive monoclonal antibodies (MAbs), 13 shared spike-specific IgH clusters were found. These shared spike-specific IgH clusters are derived from the same lineage of several recently published neutralizing MAbs, including CC12.1, CC12.3, C102, REGN10977, and 4A8. Furthermore, identical spike-specific IgH sequences were found in different COVID-19 patients, suggesting a highly convergent antibody response to SARS-CoV-2. Our analysis based on sequencing antibody repertoires from different individuals revealed key signatures of the systemic B cell response induced by SARS-CoV-2 infection.

IMPORTANCE Although the canonical delineation of serum antibody responses following SARS-CoV-2 infection has been well established, the dynamics of antibody repertoire at the mRNA transcriptional level has not been well understood, especially the correlation between serum antibody titers and the antibody mRNA transcripts. In this study, we analyzed the IgH transcripts and characterized the B cell clonal expansion and differentiation, isotype switching, and somatic hypermutation in COVID-19 patients. This study provided insights at the repertoire level for the B cell response after SARS-CoV-2 infection.

A large-scale systematic survey of SARS-CoV-2 antibodies reveals recurring molecular features

In the past two years, the global research in combating COVID-19 pandemic has led to isolation and characterization of numerous human antibodies to the SARS-CoV-2 spike. This enormous collection of antibodies provides an unprecedented opportunity to study the antibody response to a single antigen. From mining information derived from 88 research publications and 13 patents, we have assembled a dataset of ∼8,000 human antibodies to the SARS-CoV-2 spike from >200 donors. Analysis of antibody targeting of different domains of the spike protein reveals a number of common (public) responses to SARS-CoV-2, exemplified via recurring IGHV/IGK(L)V pairs, CDR H3 sequences, IGHD usage, and somatic hypermutation. We further present a proof-of-concept for prediction of antigen specificity using deep learning to differentiate sequences of antibodies to SARS-CoV-2 spike and to influenza hemagglutinin. Overall, this study not only provides an informative resource for antibody and vaccine research, but fundamentally advances our molecular understanding of public antibody responses to a viral pathogen.

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.

RAPID: A Rep-Seq Dataset Analysis Platform With an Integrated Antibody Database

The antibody repertoire is a critical component of the adaptive immune system and is believed to reflect an individual’s immune history and current immune status. Delineating the antibody repertoire has advanced our understanding of humoral immunity, facilitated antibody discovery, and showed great potential for improving the diagnosis and treatment of disease. However, no tool to date has effectively integrated big Rep-seq data and prior knowledge of functional antibodies to elucidate the remarkably diverse antibody repertoire. We developed a Rep-seq dataset Analysis Platform with an Integrated antibody Database (RAPID; https://rapid.zzhlab.org/), a free and web-based tool that allows researchers to process and analyse Rep-seq datasets. RAPID consolidates 521 WHO-recognized therapeutic antibodies, 88,059 antigen- or disease-specific antibodies, and 306 million clones extracted from 2,449 human IGH Rep-seq datasets generated from individuals with 29 different health conditions. RAPID also integrates a standardized Rep-seq dataset analysis pipeline to enable users to upload and analyse their datasets. In the process, users can also select set of existing repertoires for comparison. RAPID automatically annotates clones based on integrated therapeutic and known antibodies, and users can easily query antibodies or repertoires based on sequence or optional keywords. With its powerful analysis functions and rich set of antibody and antibody repertoire information, RAPID will benefit researchers in adaptive immune studies.

Convergent antibody responses to the SARS-CoV-2 spike protein in convalescent and vaccinated individuals

Unrelated individuals can produce genetically similar clones of antibodies, known as public clonotypes, which have been seen in responses to different infectious diseases, as well as healthy individuals. Here we identify 37 public clonotypes in memory B cells from convalescent survivors of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection or in plasmablasts from an individual after vaccination with mRNA-encoded spike protein. We identify 29 public clonotypes, including clones recognizing the receptor-binding domain (RBD) in the spike protein S1 subunit (including a neutralizing, angiotensin-converting enzyme 2 [ACE2]-blocking clone that protects in vivo) and others recognizing non-RBD epitopes that bind the S2 domain. Germline-revertant forms of some public clonotypes bind efficiently to spike protein, suggesting these common germline-encoded antibodies are preconfigured for avid recognition. Identification of large numbers of public clonotypes provides insight into the molecular basis of efficacy of SARS-CoV-2 vaccines and sheds light on the immune pressures driving the selection of common viral escape mutants.

The methods and advances of adaptive immune receptors repertoire sequencing

The adaptive immune response is a powerful tool, capable of recognizing, binding to, and neutralizing a vast number of internal and external threats via T or B lymphatic receptors with widespread sets of antigen specificities. The emergence of high-throughput sequencing technology and bioinformatics provides opportunities for research in the fields of life sciences and medicine. The analysis and annotation for immune repertoire data can reveal biologically meaningful information, including immune prediction, target antigens, and effective evaluation. Continuous improvements of the immunological repertoire sequencing methods and analysis tools will help to minimize the experimental and calculation errors and realize the immunological information to meet the clinical requirements. That said, the clinical application of adaptive immune repertoire sequencing requires appropriate experimental methods and standard analytical tools. At the population cell level, we can acquire the overview of cell groups, but the information about a single cell is not obtained accurately. The information that is ignored may be crucial for understanding the heterogeneity of each cell, gene expression and drug response. The combination of high-throughput sequencing and single-cell technology allows us to obtain single-cell information with low-cost and high-throughput. In this review, we summarized the current methods and progress in this area.

Potent and protective IGHV3-53/3-66 public antibodies and their shared escape mutant on the spike of SARS-CoV-2

Neutralizing antibodies (nAbs) to SARS-CoV-2 hold powerful potentials for clinical interventions against COVID-19 disease. However, their common genetic and biologic features remain elusive. Here we interrogate a total of 165 antibodies from eight COVID-19 patients, and find that potent nAbs from different patients have disproportionally high representation of IGHV3-53/3-66 usage, and therefore termed as public antibodies. Crystal structural comparison of these antibodies reveals they share similar angle of approach to RBD, overlap in buried surface and binding residues on RBD, and have substantial spatial clash with receptor angiotensin-converting enzyme-2 (ACE2) in binding to RBD. Site-directed mutagenesis confirms these common binding features although some minor differences are found. One representative antibody, P5A-3C8, demonstrates extraordinarily protective efficacy in a golden Syrian hamster model against SARS-CoV-2 infection. However, virus escape analysis identifies a single natural mutation in RBD, namely K417N found in B.1.351 variant from South Africa, abolished the neutralizing activity of these public antibodies. The discovery of public antibodies and shared escape mutation highlight the intricate relationship between antibody response and SARS-CoV-2, and provide critical reference for the development of antibody and vaccine strategies to overcome the antigenic variation of SARS-CoV-2.

Sequence signatures of two public antibody clonotypes that bind SARS-CoV-2 receptor binding domain

Since the COVID-19 pandemic onset, the antibody response to SARS-CoV-2 has been extensively characterized. Antibodies to the receptor binding domain (RBD) on the spike protein are frequently encoded by IGHV3-53/3-66 with a short complementarity-determining region (CDR) H3. Germline-encoded sequence motifs in heavy chain CDRs H1 and H2 have a major function, but whether any common motifs are present in CDR H3, which is often critical for binding specificity, is not clear. Here, we identify two public clonotypes of IGHV3-53/3-66 RBD antibodies with a 9-residue CDR H3 that pair with different light chains. Distinct sequence motifs on CDR H3 are present in the two public clonotypes that seem to be related to differential light chain pairing. Additionally, we show that Y58F is a common somatic hypermutation that results in increased binding affinity of IGHV3-53/3-66 RBD antibodies with a short CDR H3. These results advance understanding of the antibody response to SARS-CoV-2.

Germline IGHV3-53-encoded RBD-targeting neutralizing antibodies are commonly present in the antibody repertoires of COVID-19 patients

Monoclonal antibodies (mAbs) encoded by IGHV3-53 (VH3-53) targeting the spike receptor-binding domain (RBD) have been isolated from different COVID-19 patients. However, the existence and prevalence of shared VH3-53-encoded antibodies in the antibody repertoires is not clear. Using antibody repertoire sequencing, we found that the usage of VH3-53 increased after SARS-CoV-2 infection. A highly shared VH3-53-J6 clonotype was identified in 9 out of 13 COVID-19 patients. This clonotype was derived from convergent gene rearrangements with few somatic hypermutations and was evolutionary conserved. We synthesized 34 repertoire-deduced novel VH3-53-J6 heavy chains and paired with a common IGKV1-9 light chain to produce recombinant mAbs. Most of these recombinant mAbs (23/34) possess RBD binding and virus-neutralizing activities, and recognize ACE2 binding site via the same molecular interface. Our computational analysis, validated by laboratory experiments, revealed that VH3-53 antibodies targeting RBD are commonly present in COVID-19 patients’ antibody repertoires, indicating many people have germline-like precursor sequences to rapidly generate SARS-CoV-2 neutralizing antibodies. Moreover, antigen-specific mAbs can be digitally obtained through antibody repertoire sequencing and computational analysis.

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.

Rapid isolation and immune profiling of SARS-CoV-2 specific memory B cell in convalescent COVID-19 patients via LIBRA-seq

B cell response plays a critical role against SARS-CoV-2 infection. However, little is known about the diversity and frequency of the paired SARS-CoV-2 antigen-specific BCR repertoire after SARS-CoV-2 infection. Here, we performed single-cell RNA sequencing and VDJ sequencing using the memory and plasma B cells isolated from five convalescent COVID-19 patients, and analyzed the spectrum and transcriptional heterogeneity of antibody immune responses. Via linking BCR to antigen specificity through sequencing (LIBRA-seq), we identified a distinct activated memory B cell subgroup (CD11c^high CD95^high) had a higher proportion of SARS-CoV-2 antigen-labeled cells compared with memory B cells. Our results revealed the diversity of paired BCR repertoire and the non-stochastic pairing of SARS-CoV-2 antigen-specific immunoglobulin heavy and light chains after SARS-CoV-2 infection. The public antibody clonotypes were shared by distinct convalescent individuals. Moreover, several antibodies isolated by LIBRA-seq showed high binding affinity against SARS-CoV-2 receptor-binding domain (RBD) or nucleoprotein (NP) via ELISA assay. Two RBD-reactive antibodies C14646P3S and C2767P3S isolated by LIBRA-seq exhibited high neutralizing activities against both pseudotyped and authentic SARS-CoV-2 viruses in vitro. Our study provides fundamental insights into B cell response following SARS-CoV-2 infection at the single-cell level.

Diverse immunoglobulin gene usage and convergent epitope targeting in neutralizing antibody responses to SARS-CoV-2

It is unclear whether individuals with enormous diversity in B cell receptor repertoires are consistently able to mount effective antibody responses against SARS-CoV-2. We analyzed antibody responses in a cohort of 55 convalescent patients and isolated 54 potent neutralizing monoclonal antibodies (mAbs). While most of the mAbs target the angiotensin-converting enzyme 2 (ACE2) binding surface on the receptor binding domain (RBD) of SARS-CoV-2 spike protein, mAb 47D1 binds only to one side of the receptor binding surface on the RBD. Neutralization by 47D1 is achieved independent of interfering RBD-ACE2 binding. A crystal structure of the mAb-RBD complex shows that the IF motif at the tip of 47D1 CDR H2 interacts with a hydrophobic pocket in the RBD. Diverse immunoglobulin gene usage and convergent epitope targeting characterize neutralizing antibody responses to SARS-CoV-2, suggesting that vaccines that effectively present the receptor binding site on the RBD will likely elicit neutralizing antibody responses in a large fraction of the population.

From Structural Studies to HCV Vaccine Design

Hepatitis C virus (HCV) is a serious and growing public health problem despite recent developments of antiviral therapeutics. To achieve global elimination of HCV, an effective cross-genotype vaccine is needed. The failure of previous vaccination trials to elicit an effective cross-reactive immune response demands better vaccine antigens to induce a potent cross-neutralizing response to improve vaccine efficacy. HCV E1 and E2 envelope (Env) glycoproteins are the main targets for neutralizing antibodies (nAbs), which aid in HCV clearance and protection. Therefore, a molecular-level understanding of the nAb responses against HCV is imperative for the rational design of cross-genotype vaccine antigens. Here we summarize the recent advances in structural studies of HCV Env and Env-nAb complexes and how they improve our understanding of immune recognition of HCV. We review the structural data defining HCV neutralization epitopes and conformational plasticity of the Env proteins, and the knowledge applicable to rational vaccine design.

Convergent antibody responses to the SARS-CoV-2 spike protein in convalescent and vaccinated individuals

Unrelated individuals can produce genetically similar clones of antibodies, known as public clonotypes, which have been seen in responses to different infectious diseases as well as healthy individuals. Here we identify 37 public clonotypes in memory B cells from convalescent survivors of SARS-CoV-2 infection or in plasmablasts from an individual after vaccination with mRNA-encoded spike protein. We identified 29 public clonotypes, including clones recognizing the receptor-binding domain (RBD) in the spike protein S1 subunit (including a neutralizing, ACE2-blocking clone that protects in vivo ), and others recognizing non-RBD epitopes that bound the heptad repeat 1 region of the S2 domain. Germline-revertant forms of some public clonotypes bound efficiently to spike protein, suggesting these common germline-encoded antibodies are preconfigured for avid recognition. Identification of large numbers of public clonotypes provides insight into the molecular basis of efficacy of SARS-CoV-2 vaccines and sheds light on the immune pressures driving the selection of common viral escape mutants.

Large-scale analysis of 2,152 Ig-seq datasets reveals key features of B cell biology and the antibody repertoire

Antibody repertoire sequencing enables researchers to acquire millions of B cell receptors and investigate these molecules at the single-nucleotide level. This power and resolution in studying humoral responses have led to its wide applications. However, most of these studies were conducted with a limited number of samples. Given the extraordinary diversity, assessment of these key features with a large sample set is demanded. Thus, we collect and systematically analyze 2,152 high-quality heavy-chain antibody repertoires. Our study reveals that 52 core variable genes universally contribute to more than 99% of each individual's repertoire; a distal interspersed preferences characterize V gene recombination; the number of public clones between two repertoires follows a linear model, and the positive selection dominates at RGYW motif in somatic hypermutations. Thus, this population-level analysis resolves some critical features of the antibody repertoire and may have significant value to the large cadre of scientists.

Overcoming Immunological Challenges Limiting Capsid-Mediated Gene Therapy With Machine Learning

A key hurdle to making adeno-associated virus (AAV) capsid mediated gene therapy broadly beneficial to all patients is overcoming pre-existing and therapy-induced immune responses to these vectors. Recent advances in high-throughput DNA synthesis, multiplexing and sequencing technologies have accelerated engineering of improved capsid properties such as production yield, packaging efficiency, biodistribution and transduction efficiency. Here we outline how machine learning, advances in viral immunology, and high-throughput measurements can enable engineering of a new generation of de-immunized capsids beyond the antigenic landscape of natural AAVs, towards expanding the therapeutic reach of gene therapy.

Applications of Machine and Deep Learning in Adaptive Immunity

Adaptive immunity is mediated by lymphocyte B and T cells, which respectively express a vast and diverse repertoire of B cell and T cell receptors and, in conjunction with peptide antigen presentation through major histocompatibility complexes (MHCs), can recognize and respond to pathogens and diseased cells. In recent years, advances in deep sequencing have led to a massive increase in the amount of adaptive immune receptor repertoire data; additionally, proteomics techniques have led to a wealth of data on peptide-MHC presentation. These large-scale data sets are now making it possible to train machine and deep learning models, which can be used to identify complex and high-dimensional patterns in immune repertoires. This article introduces adaptive immune repertoires and machine and deep learning related to biological sequence data and then summarizes the many applications in this field, which span from predicting the immunological status of a host to the antigen specificity of individual receptors and the engineering of immunotherapeutics.

Shared B cell memory to coronaviruses and other pathogens varies in human age groups and tissues

It remains unclear whether B cell repertoires against coronaviruses and other pathogens differ between adults and children and how important these distinctions are. Yang et al. analyzed blood samples from young children and adults, as well as tissues from deceased organ donors, characterizing the B cell receptor (BCR) repertoires specific to six common pathogens and two viruses that they had not seen before: Ebola virus and severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Children had higher frequencies of B cells with convergent BCR heavy chains against previously encountered pathogens and higher frequencies of class-switched convergent B cell clones against SARS-CoV-2 and related coronaviruses. These findings suggest that encounters with coronaviruses in early life may produce cross-reactive memory B cell populations that contribute to divergent COVID-19 susceptibilities.

Science, this issue p. 738

Vaccination and infection promote the formation, tissue distribution, and clonal evolution of B cells, which encode humoral immune memory. We evaluated pediatric and adult blood and deceased adult organ donor tissues to determine convergent antigen-specific antibody genes of similar sequences shared between individuals. B cell memory varied for different pathogens. Polysaccharide antigenspecific clones were not exclusive to the spleen. Adults had higher clone frequencies and greater class switching in lymphoid tissues than blood, while pediatric blood had abundant class-switched convergent clones. Consistent with reported serology, prepandemic children had class-switched convergent clones to severe acute respiratory syndrome coronavirus 2 with weak cross-reactivity to other coronaviruses, while adult blood or tissues showed few such clones. These results highlight the prominence of early childhood B cell clonal expansions and cross-reactivity for future responses to novel pathogens.