A Public Database of Memory and Naive B-Cell Receptor Sequences

Benchmarking and integrating human B-cell receptor genomic and antibody proteomic profiling

Immunoglobulins (Ig), which exist either as B-cell receptors (BCR) on the surface of B cells or as antibodies when secreted, play a key role in the recognition and response to antigenic threats. The capability to jointly characterize the BCR and antibody repertoire is crucial for understanding human adaptive immunity. From peripheral blood, bulk BCR sequencing (bulkBCR-seq) currently provides the highest sampling depth, single-cell BCR sequencing (scBCR-seq) allows for paired chain characterization, and antibody peptide sequencing by tandem mass spectrometry (Ab-seq) provides information on the composition of secreted antibodies in the serum. Yet, it has not been benchmarked to what extent the datasets generated by these three technologies overlap and complement each other. To address this question, we isolated peripheral blood B cells from healthy human donors and sequenced BCRs at bulk and single-cell levels, in addition to utilizing publicly available sequencing data. Integrated analysis was performed on these datasets, resolved by replicates and across individuals. Simultaneously, serum antibodies were isolated, digested with multiple proteases, and analyzed with Ab-seq. Systems immunology analysis showed high concordance in repertoire features between bulk and scBCR-seq within individuals, especially when replicates were utilized. In addition, Ab-seq identified clonotype-specific peptides using both bulk and scBCR-seq library references, demonstrating the feasibility of combining scBCR-seq and Ab-seq for reconstructing paired-chain Ig sequences from the serum antibody repertoire. Collectively, our work serves as a proof-of-principle for combining bulk sequencing, single-cell sequencing, and mass spectrometry as complementary methods towards capturing humoral immunity in its entirety.

Clonal associations between lymphocyte subsets and functional states in rheumatoid arthritis synovium

Rheumatoid arthritis (RA) is an autoimmune disease involving antigen-specific T and B cells. Here, we perform single-cell RNA and repertoire sequencing on paired synovial tissue and blood samples from 12 seropositive RA patients. We identify clonally expanded CD4 + T cells, including CCL5+ cells and T peripheral helper (Tph) cells, which show a prominent transcriptomic signature of recent activation and effector function. CD8 + T cells show higher oligoclonality than CD4 + T cells, with the largest synovial clones enriched in GZMK+ cells. CD8 + T cells with possibly virus-reactive TCRs are distributed across transcriptomic clusters. In the B cell compartment, NR4A1+ activated B cells, and plasma cells are enriched in the synovium and demonstrate substantial clonal expansion. We identify synovial plasma cells that share BCRs with synovial ABC, memory, and activated B cells. Receptor-ligand analysis predicted IFNG and TNFRSF members as mediators of synovial Tph-B cell interactions. Together, these results reveal clonal relationships between functionally distinct lymphocyte populations that infiltrate the synovium of patients with RA.

A VRC13-like bNAb response is associated with complex escape pathways in HIV-1 envelope

The rational design of HIV-1 immunogens to trigger the development of broadly neutralizing antibodies (bNAbs) requires understanding the viral evolutionary pathways influencing this process. An acute HIV-1-infected individual exhibiting >50% plasma neutralization breadth developed neutralizing antibody specificities against the CD4-binding site (CD4bs) and V1V2 regions of Env gp120. Comparison of pseudoviruses derived from early and late autologous env sequences demonstrated the development of >2 log resistance to VRC13 but not to other CD4bs-specific bNAbs. Mapping studies indicated that the V3 and CD4-binding loops of Env gp120 contributed significantly to developing resistance to the autologous neutralizing response and that the CD4-binding loop (CD4BL) specifically was responsible for the developing resistance to VRC13. Tracking viral evolution during the development of this cross-neutralizing CD4bs response identified amino acid substitutions arising at only 4 of 11 known VRC13 contact sites (K282, T283, K421, and V471). However, each of these mutations was external to the V3 and CD4BL regions conferring resistance to VRC13 and was transient in nature. Rather, complete resistance to VRC13 was achieved through the cooperative expression of a cluster of single amino acid changes within and immediately adjacent to the CD4BL, including a T359I substitution, exchange of a potential N-linked glycosylation (PNLG) site to residue S362 from N363, and a P369L substitution. Collectively, our data characterize complex HIV-1 env evolution in an individual developing resistance to a VRC13-like neutralizing antibody response and identify novel VRC13-associated escape mutations that may be important to inducing VRC13-like bNAbs for lineage-based immunogens.IMPORTANCEThe pursuit of eliciting broadly neutralizing antibodies (bNAbs) through vaccination and their use as therapeutics remains a significant focus in the effort to eradicate HIV-1. Key to our understanding of this approach is a more extensive understanding of bNAb contact sites and susceptible escape mutations in HIV-1 envelope (env). We identified a broad neutralizer exhibiting VRC13-like responses, a non-germline restricted class of CD4-binding site antibody distinct from the well-studied VRC01-class. Through longitudinal envelope sequencing and Env-pseudotyped neutralization assays, we characterized a complex escape pathway requiring the cooperative evolution of four amino acid changes to confer complete resistance to VRC13. This suggests that VRC13-class bNAbs may be refractory to rapid escape and attractive for therapeutic applications. Furthermore, the identification of longitudinal viral changes concomitant with the development of neutralization breadth may help identify the viral intermediates needed for the maturation of VRC13-like responses and the design of lineage-based immunogens.

Probabilities of developing HIV-1 bNAb sequence features in uninfected and chronically infected individuals

HIV-1 broadly neutralizing antibodies (bNAbs) are able to suppress viremia and prevent infection. Their induction by vaccination is therefore a major goal. However, in contrast to antibodies that neutralize other pathogens, HIV-1-specific bNAbs frequently carry uncommon molecular characteristics that might prevent their induction. Here, we perform unbiased sequence analyses of B cell receptor repertoires from 57 uninfected and 46 chronically HIV-1- or HCV-infected individuals and learn probabilistic models to predict the likelihood of bNAb development. We formally show that lower probabilities for bNAbs are predictive of higher HIV-1 neutralization activity. Moreover, ranking bNAbs by their probabilities allows to identify highly potent antibodies with superior generation probabilities as preferential targets for vaccination approaches. Importantly, we find equal bNAb probabilities across infected and uninfected individuals. This implies that chronic infection is not a prerequisite for the generation of bNAbs, fostering the hope that HIV-1 vaccines can induce bNAb development in uninfected people.

VH2+ Antigen-Experienced B Cells in the Cerebrospinal Fluid Are Expanded and Enriched in Pediatric Anti-NMDA Receptor Encephalitis

Pediatric and adult autoimmune encephalitis (AE) are often associated with Abs to the NR1 subunit of the N-methyl-d-aspartate (NMDA) receptor (NMDAR). Very little is known regarding the cerebrospinal fluid humoral immune profile and Ab genetics associated with pediatric anti-NMDAR-AE. Using a combination of cellular, molecular, and immunogenetics tools, we collected cerebrospinal fluid from pediatric subjects and generated 1) flow cytometry data to calculate the frequency of B cell subtypes in the cerebrospinal fluid of pediatric subjects with anti-NMDAR-AE and controls, 2) a panel of recombinant human Abs from a pediatric case of anti-NMDAR-AE that was refractory to treatment, and 3) a detailed analysis of the Ab genes that bound the NR1 subunit of the NMDAR. Ag-experienced B cells including memory cells, plasmablasts, and Ab-secreting cells were expanded in the pediatric anti-NMDAR-AE cohort, but not in the controls. These Ag-experienced B cells in the cerebrospinal fluid of a pediatric case of NMDAR-AE that was refractory to treatment had expanded use of variable H chain family 2 (VH2) genes with high somatic hypermutation that all bound to the NR1 subunit of the NMDAR. A CDR3 motif was identified in this refractory case that likely drove early stage activation and expansion of naive B cells to Ab-secreting cells, facilitating autoimmunity associated with pediatric anti-NMDAR-AE through the production of Abs that bind NR1. These features of humoral immune responses in the cerebrospinal fluid of pediatric anti-NMDAR-AE patients may be relevant for clinical diagnosis and treatment.

Mathematical Characterization of Private and Public Immune Receptor Sequences

Diverse T and B cell repertoires play an important role in mounting effective immune responses against a wide range of pathogens and malignant cells. The number of unique T and B cell clones is characterized by T and B cell receptors (TCRs and BCRs), respectively. Although receptor sequences are generated probabilistically by recombination processes, clinical studies found a high degree of sharing of TCRs and BCRs among different individuals. In this work, we use a general probabilistic model for T/B cell receptor clone abundances to define "publicness" or "privateness" and information-theoretic measures for comparing the frequency of sampled sequences observed across different individuals. We derive mathematical formulae to quantify the mean and the variances of clone richness and overlap. Our results can be used to evaluate the effect of different sampling protocols on abundances of clones within an individual as well as the commonality of clones across individuals. Using synthetic and empirical TCR amino acid sequence data, we perform simulations to study expected clonal commonalities across multiple individuals. Based on our formulae, we compare these simulated results with the analytically predicted mean and variances of the repertoire overlap. Complementing the results on simulated repertoires, we derive explicit expressions for the richness and its uncertainty for specific, single-parameter truncated power-law probability distributions. Finally, the information loss associated with grouping together certain receptor sequences, as is done in spectratyping, is also evaluated. Our approach can be, in principle, applied under more general and mechanistically realistic clone generation models.

Early activation of inflammatory pathways in UBA1-mutated hematopoietic stem and progenitor cells in VEXAS

VEXAS (vacuoles, E1 enzyme, X-linked, autoinflammatory, somatic) syndrome is a pleiotropic, severe autoinflammatory disease caused by somatic mutations in the ubiquitin-like modifier activating enzyme 1 (UBA1) gene. To elucidate VEXAS pathophysiology, we performed transcriptome sequencing of single bone marrow mononuclear cells and hematopoietic stem and progenitor cells (HSPCs) from VEXAS patients. HSPCs are biased toward myeloid (granulocytic) differentiation, and against lymphoid differentiation in VEXAS. Activation of multiple inflammatory pathways (interferons and tumor necrosis factor alpha) occurs ontogenically early in primitive hematopoietic cells and particularly in the myeloid lineage in VEXAS, and inflammation is prominent in UBA1-mutated cells. Dysregulation in protein degradation likely leads to higher stress response in VEXAS HSPCs, which positively correlates with inflammation. TCR usage is restricted and there are increased cytotoxicity and IFN-γ signaling in T cells. In VEXAS syndrome, both aberrant inflammation and myeloid predominance appear intrinsic to hematopoietic stem cells mutated in UBA1.

AIRRSHIP: simulating human B cell receptor repertoire sequences

SUMMARY

Adaptive Immune Receptor Repertoire Sequencing is a rapidly developing field that has advanced understanding of the role of the adaptive immune system in health and disease. Numerous tools have been developed to analyse the complex data produced by this technique but work to compare their accuracy and reliability has been limited. Thorough, systematic assessment of their performance is dependent on the ability to produce high quality simulated datasets with known ground truth. We have developed AIRRSHIP, a flexible and fast Python package that produces synthetic human B cell receptor sequences. AIRRSHIP uses a comprehensive set of reference data to replicate key mechanisms in the immunoglobulin recombination process, with a particular focus on junctional complexity. Repertoires generated by AIRRSHIP are highly similar to published data and all steps in the sequence generation process are recorded. These data can be used to not only determine the accuracy of repertoire analysis tools but can also, by tuning of the large number of user-controllable parameters, give insight into factors that contribute to inaccuracies in results.

AVAILABILITY AND IMPLEMENTATION

AIRRSHIP is implemented in Python. It is available via https://github.com/Cowanlab/airrship and on PyPI at https://pypi.org/project/airrship/. Documentation can be found at https://airrship.readthedocs.io/.

Clonal associations of lymphocyte subsets and functional states revealed by single cell antigen receptor profiling of T and B cells in rheumatoid arthritis synovium

Rheumatoid arthritis (RA) is an autoimmune disease initiated by antigen-specific T cells and B cells, which promote synovial inflammation through a complex set of interactions with innate immune and stromal cells. To better understand the phenotypes and clonal relationships of synovial T and B cells, we performed single-cell RNA and repertoire sequencing on paired synovial tissue and peripheral blood samples from 12 donors with seropositive RA ranging from early to chronic disease. Paired transcriptomic-repertoire analyses highlighted 3 clonally distinct CD4 T cells populations that were enriched in RA synovium: T peripheral helper (Tph) and T follicular helper (Tfh) cells, CCL5+ T cells, and T regulatory cells (Tregs). Among these cells, Tph cells showed a unique transcriptomic signature of recent T cell receptor (TCR) activation, and clonally expanded Tph cells expressed an elevated transcriptomic effector signature compared to non-expanded Tph cells. CD8 T cells showed higher oligoclonality than CD4 T cells, and the largest CD8 T cell clones in synovium were highly enriched in GZMK+ cells. TCR analyses revealed CD8 T cells with likely viral-reactive TCRs distributed across transcriptomic clusters and definitively identified MAIT cells in synovium, which showed transcriptomic features of TCR activation. Among B cells, non-naive B cells including age-associated B cells (ABC), NR4A1+ activated B cells, and plasma cells, were enriched in synovium and had higher somatic hypermutation rates compared to blood B cells. Synovial B cells demonstrated substantial clonal expansion, with ABC, memory, and activated B cells clonally linked to synovial plasma cells. Together, these results reveal clonal relationships between functionally distinct lymphocyte populations that infiltrate RA synovium.

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

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

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.

Unconstrained generation of synthetic antibody-antigen structures to guide machine learning methodology for antibody specificity prediction

Machine learning (ML) is a key technology for accurate prediction of antibody-antigen binding. Two orthogonal problems hinder the application of ML to antibody-specificity prediction and the benchmarking thereof: the lack of a unified ML formalization of immunological antibody-specificity prediction problems and the unavailability of large-scale synthetic datasets to benchmark real-world relevant ML methods and dataset design. Here we developed the Absolut! software suite that enables parameter-based unconstrained generation of synthetic lattice-based three-dimensional antibody-antigen-binding structures with ground-truth access to conformational paratope, epitope and affinity. We formalized common immunological antibody-specificity prediction problems as ML tasks and confirmed that for both sequence- and structure-based tasks, accuracy-based rankings of ML methods trained on experimental data hold for ML methods trained on Absolut!-generated data. The Absolut! framework has the potential to enable real-world relevant development and benchmarking of ML strategies for biotherapeutics design.

Humoral immunity to SARS-CoV-2 elicited by combination COVID-19 vaccination regimens

The SARS-CoV-2 pandemic prompted a global vaccination effort and the development of numerous COVID-19 vaccines at an unprecedented scale and pace. As a result, current COVID-19 vaccination regimens comprise diverse vaccine modalities, immunogen combinations, and dosing intervals. Here, we compare vaccine-specific antibody and memory B cell responses following two-dose mRNA, single-dose Ad26.COV.2S, and two-dose ChAdOx1, or combination ChAdOx1/mRNA vaccination. Plasma-neutralizing activity, as well as the magnitude, clonal composition, and antibody maturation of the RBD-specific memory B cell compartments, showed substantial differences between the vaccination regimens. While individual monoclonal antibodies derived from memory B cells exhibited similar binding affinities and neutralizing potency against Wuhan-Hu-1 SARS-CoV-2, there were significant differences in epitope specificity and neutralizing breadth against viral variants of concern. Although the ChAdOx1 vaccine was inferior to mRNA and Ad26.COV.2S in several respects, biochemical and structural analyses revealed enrichment in a subgroup of memory B cell neutralizing antibodies with distinct RBD-binding properties resulting in remarkable potency and breadth.

NoisET: Noise Learning and Expansion Detection of T-Cell Receptors

High-throughput sequencing of T- and B-cell receptors makes it possible to track immune repertoires across time, in different tissues, in acute and chronic diseases and in healthy individuals. However, quantitative comparison between repertoires is confounded by variability in the read count of each receptor clonotype due to sampling, library preparation, and expression noise. We review methods for accounting for both biological and experimental noise and present an easy-to-use python package NoisET that implements and generalizes a previously developed Bayesian method. It can be used to learn experimental noise models for repertoire sequencing from replicates, and to detect responding clones following a stimulus. We test the package on different repertoire sequencing technologies and data sets. We review how such approaches have been used to identify responding clonotypes in vaccination and disease data. Availability: NoisET is freely available to use with source code at github.com/statbiophys/NoisET.

Allelic polymorphism controls autoreactivity and vaccine elicitation of human broadly neutralizing antibodies against influenza virus

Human broadly neutralizing antibodies (bnAbs) targeting the hemagglutinin stalk of group 1 influenza A viruses (IAVs) are biased for IGHV1-69 alleles that use phenylalanine (F54) but not leucine (L54) within their CDRH2 loops. Despite this, we demonstrated that both alleles encode for human IAV bnAbs that employ structurally convergent modes of contact to the same epitope. To resolve differences in lineage expandability, we compared F54 versus L54 as substrate within humanized mice, where antibodies develop with human-like CDRH3 diversity but are restricted to single VH genes. While both alleles encoded for bnAb precursors, only F54 IGHV1-69 supported elicitation of heterosubtypic serum bnAbs following immunization with a stalk-only nanoparticle vaccine. L54 IGHV1-69 was unproductive, co-encoding for anergic B cells and autoreactive stalk antibodies that were cleared from B cell memory. Moreover, human stalk antibodies also demonstrated L54-dependent autoreactivity. Therefore, IGHV1-69 polymorphism, which is skewed ethnically, gates tolerance and vaccine expandability of influenza bnAbs.

Repertoire-scale measures of antigen binding

Antibodies and T cell receptors (TCRs) are the fundamental building blocks of adaptive immunity. Repertoire-scale functionality derives from their epitope-binding properties, just as macroscopic properties like temperature derive from microscopic molecular properties. However, most approaches to repertoire-scale measurement, including sequence diversity and entropy, are not based on antibody or TCR function in this way. Thus, they potentially overlook key features of immunological function. Here we present a framework that describes repertoires in terms of the epitope-binding properties of their constituent antibodies and TCRs, based on analysis of thousands of antibody-antigen and TCR-peptide-major-histocompatibility-complex binding interactions and over 400 high-throughput repertoires. We show that repertoires consist of loose overlapping classes of antibodies and TCRs with similar binding properties. We demonstrate the potential of this framework to distinguish specific responses vs. bystander activation in influenza vaccinees, stratify cytomegalovirus (CMV)-infected cohorts, and identify potential immunological "super-agers." Classes add a valuable dimension to the assessment of immune function.

Antibody evolution to SARS-CoV-2 after single-dose Ad26.COV2.S vaccine in humans

The single-dose Ad.26.COV.2 (Janssen) vaccine elicits lower levels of neutralizing antibodies and shows more limited efficacy in protection against infection than either of the two available mRNA vaccines. In addition, Ad.26.COV.2 has been less effective in protection against severe disease during the Omicron surge. Here, we examined the memory B cell response to single-dose Ad.26.COV.2 vaccination. Compared with mRNA vaccines, Ad.26.COV.2 recipients had significantly lower numbers of RBD-specific memory B cells 1.5 or 6 mo after vaccination. Despite the lower numbers, the overall quality of the memory B cell responses appears to be similar, such that memory antibodies elicited by both vaccine types show comparable neutralizing potency against SARS-CoV-2 Wuhan-Hu-1, Delta, and Omicron BA.1 variants. The data help explain why boosting Ad.26.COV.2 vaccine recipients with mRNA vaccines is effective and why the Ad26.COV2.S vaccine can maintain some protective efficacy against severe disease during the Omicron surge.

Benchmarking computational methods for B-cell receptor reconstruction from single-cell RNA-seq data

Multiple methods have recently been developed to reconstruct full-length B-cell receptors (BCRs) from single-cell RNA sequencing (scRNA-seq) data. This need emerged from the expansion of scRNA-seq techniques, the increasing interest in antibody-based drug development and the importance of BCR repertoire changes in cancer and autoimmune disease progression. However, a comprehensive assessment of performance-influencing factors such as the sequencing depth, read length or number of somatic hypermutations (SHMs) as well as guidance regarding the choice of methodology is still lacking. In this work, we evaluated the ability of six available methods to reconstruct full-length BCRs using one simulated and three experimental SMART-seq datasets. In addition, we validated that the BCRs assembled in silico recognize their intended targets when expressed as monoclonal antibodies. We observed that methods such as BALDR, BASIC and BRACER showed the best overall performance across the tested datasets and conditions, whereas only BASIC demonstrated acceptable results on very short read libraries. Furthermore, the de novo assembly-based methods BRACER and BALDR were the most accurate in reconstructing BCRs harboring different degrees of SHMs in the variable domain, while TRUST4, MiXCR and BASIC were the fastest. Finally, we propose guidelines to select the best method based on the given data characteristics.

Analysis of memory B cells identifies conserved neutralizing epitopes on the N-terminal domain of variant SARS-Cov-2 spike proteins

SARS-CoV-2 infection or vaccination produces neutralizing antibody responses that contribute to better clinical outcomes. The receptor-binding domain (RBD) and the N-terminal domain (NTD) of the spike trimer (S) constitute the two major neutralizing targets for antibodies. Here, we use NTD-specific probes to capture anti-NTD memory B cells in a longitudinal cohort of infected individuals, some of whom were vaccinated. We found 6 complementation groups of neutralizing antibodies. 58% targeted epitopes outside the NTD supersite, 58% neutralized either Gamma or Omicron, and 14% were broad neutralizers that also neutralized Omicron. Structural characterization revealed that broadly active antibodies targeted three epitopes outside the NTD supersite including a class that recognized both the NTD and SD2 domain. Rapid recruitment of memory B cells producing these antibodies into the plasma cell compartment upon re-infection likely contributes to the relatively benign course of subsequent infections with SARS-CoV-2 variants, including Omicron.

Learning the statistics and landscape of somatic mutation-induced insertions and deletions in antibodies

Affinity maturation is crucial for improving the binding affinity of antibodies to antigens. This process is mainly driven by point substitutions caused by somatic hypermutations of the immunoglobulin gene. It also includes deletions and insertions of genomic material known as indels. While the landscape of point substitutions has been extensively studied, a detailed statistical description of indels is still lacking. Here we present a probabilistic inference tool to learn the statistics of indels from repertoire sequencing data, which overcomes the pitfalls and biases of standard annotation methods. The model includes antibody-specific maturation ages to account for variable mutational loads in the repertoire. After validation on synthetic data, we applied our tool to a large dataset of human immunoglobulin heavy chains. The inferred model allows us to identify universal statistical features of indels in heavy chains. We report distinct insertion and deletion hotspots, and show that the distribution of lengths of indels follows a geometric distribution, which puts constraints on future mechanistic models of the hypermutation process.

Affinity maturation of B cell receptors is an important mechanism by which our body designs neutralizing antibodies to defend us against pathogens, including broadly neutralizing antibodies, which target a wide range of variants of the same pathogen. Such antibodies often contain key insertions and deletions to the germline gene, or “indels”, which are caused by somatic hypermutations. However, the mechanism, frequency and role of these indels are still elusive. We designed a computational method based on a probabilistic framework to infer the characteristics of this mutational process from high-throughput antibody sequencing experiments. Applied to human data, our approach provides a comprehensive quantitative description of insertions and deletions, opening avenues for better understanding the process of affinity maturation and the design of vaccines for eliciting a broad antibody response.

Current advances in biopharmaceutical informatics: guidelines, impact and challenges in the computational developability assessment of antibody therapeutics

Therapeutic monoclonal antibodies and their derivatives are key components of clinical pipelines in the global biopharmaceutical industry. The availability of large datasets of antibody sequences, structures, and biophysical properties is increasingly enabling the development of predictive models and computational tools for the "developability assessment" of antibody drug candidates. Here, we provide an overview of the antibody informatics tools applicable to the prediction of developability issues such as stability, aggregation, immunogenicity, and chemical degradation. We further evaluate the opportunities and challenges of using biopharmaceutical informatics for drug discovery and optimization. Finally, we discuss the potential of developability guidelines based on in silico metrics that can be used for the assessment of antibody stability and manufacturability.

Conserved Neutralizing Epitopes on the N-Terminal Domain of Variant SARS-CoV-2 Spike Proteins

SARS-CoV-2 infection or vaccination produces neutralizing antibody responses that contribute to better clinical outcomes. The receptor binding domain (RBD) and the N-terminal domain (NTD) of the spike trimer (S) constitute the two major neutralizing targets for the antibody system. Neutralizing antibodies targeting the RBD bind to several different sites on this domain. In contrast, most neutralizing antibodies to NTD characterized to date bind to a single supersite, however these antibodies were obtained by methods that were not NTD specific. Here we use NTD specific probes to focus on anti-NTD memory B cells in a cohort of pre-omicron infected individuals some of which were also vaccinated. Of 275 NTD binding antibodies tested 103 neutralized at least one of three tested strains: Wuhan-Hu-1, Gamma, or PMS20, a synthetic variant which is extensively mutated in the NTD supersite. Among the 43 neutralizing antibodies that were further characterized, we found 6 complementation groups based on competition binding experiments. 58% targeted epitopes outside the NTD supersite, and 58% neutralized either Gamma or Omicron, but only 14% were broad neutralizers. Three of the broad neutralizers were characterized structurally. C1520 and C1791 recognize epitopes on opposite faces of the NTD with a distinct binding pose relative to previously described antibodies allowing for greater potency and cross-reactivity with 7 different variants including Beta, Delta, Gamma and Omicron. Antibody C1717 represents a previously uncharacterized class of NTD-directed antibodies that recognizes the viral membrane proximal side of the NTD and SD2 domain, leading to cross-neutralization of Beta, Gamma and Omicron. We conclude SARS-CoV-2 infection and/or Wuhan-Hu-1 mRNA vaccination produces a diverse collection of memory B cells that produce anti-NTD antibodies some of which can neutralize variants of concern. Rapid recruitment of these cells into the antibody secreting plasma cell compartment upon re-infection likely contributes to the relatively benign course of subsequent infections with SARS-CoV-2 variants including omicron.

NGS Analysis of Clonality and Minimal Residual Disease in a Patient with Concurrent Richter's Transformation and CLL/SLL

B-cell lymphomas are neoplastic proliferations of clonal B lymphocytes. Clonality is generally determined by PCR amplification of VDJ rearrangements in the IgH heavy chain or VJ rearrangements in Igκ/Igλ light chain genes followed by capillary electrophoresis. More recently, next-generation sequencing (NGS) has been used to detect clonality in B-cell lymphomas because of the exponential amount of information that is obtained beyond just detecting a clonal population. The additional information obtained is useful for diagnostic confirmation, prognosis assessment, and response to therapy. In this study, we utilized NGS analysis to characterize two histologically distinct lymphomas (DLBCL and CLL/SLL) that were detected contemporaneously in an asymptomatic patient. NGS analysis showed that the same VDJ rearrangement was present in nodal (DLBCL) and marrow (CLL/SLL) biopsies confirming that the DLBCL resulted from Richter's transformation of a subclinical CLL/SLL. The V region of the rearrangement remained unmutated without somatic hypermutation. In silico analysis showed that the HCDR3 sequence was heterogeneous and not stereotypic. Minimal residual disease analysis by NGS showed that the tumor clone decreased by 2.84 logs in the bone marrow after R-CHOP therapy. However, a small number of tumor cells were still detected in the peripheral blood after R-CHOP therapy. Subsequent allogeneic transplantation was successful in eradicating the tumor clone and achieving deep molecular remission. We show that NGS analysis facilitated clinical management in our patient by helping to characterize the VDJ rearrangement in detail and by tracking minimal residual disease with high sensitivity and specificity.

The human antibody sequence space and structural design of the V, J regions, and CDRH3 with Rosetta

The human adaptive immune response enables the targeting of epitopes on pathogens with high specificity. Infection with a pathogen induces somatic hyper-mutation and B-cell selection processes that govern the shape and diversity of the antibody sequence landscape. To date, even the largest immunome repertoires of adaptive immune receptors acquired by next-generation sequencing cannot fully capture the vast antibody sequence space of a single individual, which is estimated to be at least 1012 potential sequences. Degeneracy of the genetic code means that the number of possible nucleotide triplets (64) is greater than the number of canonical amino acids (20), resulting in some amino acids being encoded by multiple triplets and different amino acids sharing the same nucleotide in 1 or 2 positions in the triplet. We hypothesize that the degeneracy of the genetic code can be used to statistically model an enlarged space of human antibody amino acid sequences, accommodating for the discrepancy between the observed and the hypothesized antibody sequence space. Facilitated by Bayesian statistics and immunome repertoire clustering, we calculated amino acid probabilities from single nucleotide frequencies to infer a human amino acid sequence space that is used to design human-like antibodies with Rosetta. We show that antibodies designed with our restraints are on average up to 16.6% more human-like in the V and J regions compared to the Rosetta designs produced without constraints. The human-likeness of the heavy-chain CDR3 region (CDRH3) could be increased for 8 of 27 antibodies compared to Rosetta designs with a similar number of mutations and could be successfully applied on Mus musculus antibodies to demonstrate humanization.

In silico proof of principle of machine learning-based antibody design at unconstrained scale

Generative machine learning (ML) has been postulated to become a major driver in the computational design of antigen-specific monoclonal antibodies (mAb). However, efforts to confirm this hypothesis have been hindered by the infeasibility of testing arbitrarily large numbers of antibody sequences for their most critical design parameters: paratope, epitope, affinity, and developability. To address this challenge, we leveraged a lattice-based antibody-antigen binding simulation framework, which incorporates a wide range of physiological antibody-binding parameters. The simulation framework enables the computation of synthetic antibody-antigen 3D-structures, and it functions as an oracle for unrestricted prospective evaluation and benchmarking of antibody design parameters of ML-generated antibody sequences. We found that a deep generative model, trained exclusively on antibody sequence (one dimensional: 1D) data can be used to design conformational (three dimensional: 3D) epitope-specific antibodies, matching, or exceeding the training dataset in affinity and developability parameter value variety. Furthermore, we established a lower threshold of sequence diversity necessary for high-accuracy generative antibody ML and demonstrated that this lower threshold also holds on experimental real-world data. Finally, we show that transfer learning enables the generation of high-affinity antibody sequences from low-N training data. Our work establishes a priori feasibility and the theoretical foundation of high-throughput ML-based mAb design.

Individualized VDJ recombination predisposes the available Ig sequence space

The process of recombination between variable (V), diversity (D), and joining (J) immunoglobulin (Ig) gene segments determines an individual's naive Ig repertoire and, consequently, (auto)antigen recognition. VDJ recombination follows probabilistic rules that can be modeled statistically. So far, it remains unknown whether VDJ recombination rules differ between individuals. If these rules differed, identical (auto)antigen-specific Ig sequences would be generated with individual-specific probabilities, signifying that the available Ig sequence space is individual specific. We devised a sensitivity-tested distance measure that enables inter-individual comparison of VDJ recombination models. We discovered, accounting for several sources of noise as well as allelic variation in Ig sequencing data, that not only unrelated individuals but also human monozygotic twins and even inbred mice possess statistically distinguishable immunoglobulin recombination models. This suggests that, in addition to genetic, there is also nongenetic modulation of VDJ recombination. We demonstrate that population-wide individualized VDJ recombination can result in orders of magnitude of difference in the probability to generate (auto)antigen-specific Ig sequences. Our findings have implications for immune receptor-based individualized medicine approaches relevant to vaccination, infection, and autoimmunity.

B Cell Receptor Repertoire Analysis in Malaria-Naive and Malaria-Experienced Individuals Reveals Unique Characteristics of Atypical Memory B Cells

Malaria, caused by parasites of the Plasmodium genus, is responsible for significant morbidity and mortality globally. Chronic Plasmodium falciparum exposure affects the B cell compartment, leading to the accumulation of atypical memory B cells (atMBCs). IgM-positive (IgM+) and IgG+ atMBCs have not been compared in-depth in the context of malaria, nor is it known if atMBCs in malaria-experienced individuals are different from phenotypically similar B cells in individuals with no known history of Plasmodium exposure. To address these questions, we characterized the B cell receptor (BCR) repertoire of naive B cells (NBCs), IgM+ and IgG+ classical MBCs (cMBCs), and IgM+ and IgG+ atMBCs from 13 malaria-naive American adults and 7 malaria-experienced Ugandan adults. Our results demonstrate that P. falciparum exposure mainly drives changes in atMBCs. In comparison to malaria-naive adults, the BCR repertoire of Plasmodium-exposed adults showed increased levels of somatic hypermutation in the heavy chain V region in IgM+ and IgG+ atMBCs, shorter heavy chain complementarity-determining region 3 (HCDR3) in IgG+ atMBCs, and increased usage of IGHV3-73 in IgG+ cMBCs and both IgM+ and IgG+ atMBCs. Irrespective of Plasmodium exposure, IgM+ atMBCs closely resembled NBCs, while IgG+ atMBCs resembled IgG+ cMBCs. Physicochemical properties of the HCDR3 seemed to be intrinsic to cell type and independent of malaria experience. The resemblance between atMBCs from Plasmodium-exposed and naive adults suggests similar differentiation pathways regardless of chronic antigen exposure. Moreover, these data demonstrate that IgM+ and IgG+ atMBCs are distinct populations that should be considered separately in future analyses. IMPORTANCE Malaria, caused by Plasmodium parasites, still contributes to a high global burden of disease, mainly in children under 5 years of age. Chronic and recurrent Plasmodium infections affect the development of B cell memory against the parasite and promote the accumulation of atypical memory B cells (atMBCs), which have an unclear function in the immune response. Understanding where these cells originate from and whether they are beneficial in the immune response to Plasmodium will help inform vaccination development efforts. We found differences in B cell receptor (BCR) properties of atMBCs between malaria-naive and malaria-experienced adults that are suggestive of divergent selection processes, resulting in more somatic hypermutation and differential immunoglobulin heavy chain V (IGHV) gene usage. Despite these differences, atMBCs from malaria-naive and malaria-experienced adults also showed many similarities in BCR characteristics, such as physicochemical properties of the HCDR3 region, suggesting that atMBCs undergo similar differentiation pathways in response to different pathogens. Our study provides new insights into the effects of malaria experience on the B cell compartment and the relationships between atMBCs and other B cell populations.

An Allogeneic Multiple Myeloma GM-CSF-Secreting Vaccine with Lenalidomide Induces Long-term Immunity and Durable Clinical Responses in Patients in Near Complete Remission

PURPOSE

This proof-of-principle clinical trial evaluated whether an allogeneic multiple myeloma GM-CSF-secreting vaccine (MM-GVAX) in combination with lenalidomide could deepen the clinical response in patients with multiple myeloma in sustained near complete remission (nCR).

PATIENTS AND METHODS

Fifteen patients on lenalidomide were treated with MM-GVAX and pneumococcal conjugate vaccine (PCV; Prevnar) at 1, 2, 3, and 6 months.

RESULTS

Eight patients (53.3%) achieved a true CR. With a median follow-up of 5 years, the median progression-free survival had not been reached, and the median overall survival was 7.8 years from enrollment. MM-GVAX induced clonal T-cell expansion and measurable cytokine responses that persisted up to 7 years in all patients. At baseline, a higher minimal residual disease was predictive of early relapse. After vaccination, a lack of both CD27^-DNAM1^-CD8⁺ T cells and antigen-presenting cells was associated with disease progression.

CONCLUSIONS

MM-GVAX, along with lenalidomide, effectively primed durable immunity and resulted in long-term disease control, as suggested by the reappearance of a detectable, fluctuating M-spike without meeting the criteria for clinical relapse. For patients in a nCR, MM-GVAX administration was safe and resulted in prolonged clinical responses.

Different B cell subpopulations show distinct patterns in their IgH repertoire metrics

Several human B cell subpopulations are recognised in the peripheral blood, which play distinct roles in the humoral immune response. These cells undergo developmental and maturational changes involving VDJ recombination, somatic hypermutation and class switch recombination, altogether shaping their immunoglobulin heavy chain (IgH) repertoire. Here, we sequenced the IgH repertoire of naïve, marginal zone, switched and plasma cells from 10 healthy adults along with matched unsorted and in silico separated CD19⁺ bulk B cells. Using advanced bioinformatic analysis and machine learning, we show that sorted B cell subpopulations are characterised by distinct repertoire characteristics on both the individual sequence and the repertoire level. Sorted subpopulations shared similar repertoire characteristics with their corresponding in silico separated subsets. Furthermore, certain IgH repertoire characteristics correlated with the position of the constant region on the IgH locus. Overall, this study provides unprecedented insight over mechanisms of B cell repertoire control in peripherally circulating B cell subpopulations.

TCR+/BCR+ dual-expressing cells and their associated public BCR clonotype are not enriched in type 1 diabetes

Ahmed and colleagues recently described a novel hybrid lymphocyte expressing both a B and T cell receptor, termed double expresser (DE) cells. DE cells in blood of type 1 diabetes (T1D) subjects were present at increased numbers and enriched for a public B cell clonotype. Here, we attempted to reproduce these findings. While we could identify DE cells by flow cytometry, we found no association between DE cell frequency and T1D status. We were unable to identify the reported public B cell clone, or any similar clone, in bulk B cells or sorted DE cells from T1D subjects or controls. We also did not observe increased usage of the public clone VH or DH genes in B cells or in sorted DE cells. Taken together, our findings suggest that DE cells and their alleged public clonotype are not enriched in T1D. This Matters Arising paper is in response to Ahmed et al. (2019), published in Cell. See also the response by Ahmed et al. (2021), published in this issue.

Characterization of ascites- and tumor-infiltrating γδ T cells reveals distinct repertoires and a beneficial role in ovarian cancer

The role of γδ T cells in antitumor immunity has been under investigation for the past two decades, but little is known about their contribution to clinical outcomes in patients. Here, we set out to define the clonotypic, phenotypic, and functional features of γδ T cells in peripheral blood, ascites, and metastatic tumor tissue from patients with advanced epithelial ovarian cancer. T cell receptor (TCR) sequencing of the γ chain revealed that tumor-infiltrating γδ T cells have a unique and skewed repertoire with high TCR diversity and low clonality. In contrast, ascites-derived γδ T cells presented a lower TCR diversity and higher clonality, suggesting a TCR-dependent clonal focusing at this site. Further investigation showed that tumor samples had abundant γδ T cells with a tissue-resident, activation-associated phenotype, less usage of Vγ9 and an impaired response to adaptive-associated stimuli, implying an innate-like activation pathway, rather than an adaptive TCR-engaging pathway, at these tumor sites. Furthermore, high γδ T cell cytokine responsiveness upon stimulation was associated with a favorable outcome for patients in terms of both overall survival and reduced residual tumor burden after primary surgery. Last, the functionality of γδ T cells and patient survival were negatively affected by the proportions of CD39-expressing T cells, highlighting the potential of CD39 as a target to improve γδ T cell responses and unleash their antitumor capabilities.

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

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

Human transitional and IgMlow mature naïve B cells preserve permissive B-cell receptors

The level of immunoglobulin M (IgM) displayed on the surface of peripheral blood B cells exhibits a broad dynamic range and has been associated with both development and selection. To determine whether IgM surface expression associates with distinct immunoglobulin heavy‐chain (IGH) repertoire properties, we performed deep IgM sequencing of peripheral blood transitional and mature naïve B cells in the upper and lower quartiles of surface IgM expression for 12 healthy donors. Mature naïve B cells within the lowest quartile for surface IgM expression displayed more diverse IGH features including increased complementarity‐determining region 3 length, IGHJ6 segment usage and aromatic amino acids compared with mature naïve B cells with high surface IgM. There were no differences between IGH repertoires for transitional B cells with high or low surface IgM. These findings suggest that a selection checkpoint during progression of transitional to mature naïve B cells reduces the breadth of the IGH repertoire among high surface IgM B cells but that diversity is preserved in B cells expressing low levels of surface IgM.

A Single Human VH-gene Allows for a Broad-Spectrum Antibody Response Targeting Bacterial Lipopolysaccharides in the Blood

B cell receptors (BCRs) display a combination of variable (V)-gene-encoded complementarity determining regions (CDRs) and adaptive/hypervariable CDR3 loops to engage antigens. It has long been proposed that the former tune for recognition of pathogens or groups of pathogens. To experimentally evaluate this within the human antibody repertoire, we perform immune challenges in transgenic mice that bear diverse human CDR3 and light chains but are constrained to different human V_H-genes. We find that, of six commonly deployed V_H sequences, only those CDRs encoded by IGHV1-2^∗02 enable polyclonal antibody responses against bacterial lipopolysaccharide (LPS) when introduced to the bloodstream. The LPS is from diverse strains of gram-negative bacteria, and the V_H-gene-dependent responses are directed against the non-variable and universal saccrolipid substructure of this antigen. This reveals a broad-spectrum anti-LPS response in which germline-encoded CDRs naturally hardwire the human antibody repertoire for recognition of a conserved microbial target.

Human-likeness of antibody biologics determined by back-translation and comparison with large antibody variable gene repertoires

The antibody (Ab) germline gene rearrangement of variable (V), diversity (D), and joining (J) gene segments, as well as somatic hypermutation, give rise to the human Ab variable gene sequence repertoire. It is common to characterize single nucleotide frequencies of the variable region by alignment to species-specific wildtype germline genes. The increasing application of next-generation sequencing to immune repertoire studies has led to the compilation of increasing large adaptive immunome receptor repertoire datasets. We have developed a method that maps the sequence of a target Ab onto an immunome dataset of 326 million human Ab sequences. For this purpose, we created a position- and gene-specific scoring matrix (PGSSM) and its corresponding antibody similarity score. We characterized our PGSSM score and found that it strongly correlated with the phylogenetic distance of 181,355 Ab sequences from GenBank across 20 species. The most likely human nucleotide back-translation was obtained given only PGSSMs and the amino acid sequence of an Ab achieving a nucleotide sequence recovery of 95.9% and 97.2% for human heavy and light chains, respectively. In conclusion, the scoring of our back-translation is a valuable estimate for the similarity of an Ab sequence to the natural human repertoire. As expected, Ab therapeutic molecules developed from a human source showed a higher similarity to the repertoire than engineered Abs. Thus, the PGSSM metric introduced here can be used to engineer human-like Ab therapeutics.

Broad and potent neutralizing human antibodies to tick-borne flaviviruses protect mice from disease

Tick-borne encephalitis virus (TBEV) is an emerging human pathogen that causes potentially fatal disease with no specific treatment. Mouse monoclonal antibodies are protective against TBEV, but little is known about the human antibody response to infection. Here, we report on the human neutralizing antibody response to TBEV in a cohort of infected and vaccinated individuals. Expanded clones of memory B cells expressed closely related anti-envelope domain III (EDIII) antibodies in both groups of volunteers. However, the most potent neutralizing antibodies, with IC50s below 1 ng/ml, were found only in individuals who recovered from natural infection. These antibodies also neutralized other tick-borne flaviviruses, including Langat, louping ill, Omsk hemorrhagic fever, Kyasanur forest disease, and Powassan viruses. Structural analysis revealed a conserved epitope near the lateral ridge of EDIII adjoining the EDI-EDIII hinge region. Prophylactic or early therapeutic antibody administration was effective at low doses in mice that were lethally infected with TBEV.

Landscape of Non-canonical Cysteines in Human VH Repertoire Revealed by Immunogenetic Analysis

Human antibody repertoire data captured through next-generation sequencing (NGS) has enabled deeper insights into B cell immunogenetics and paratope diversity. By analyzing large public NGS datasets, we map the landscape of non-canonical cysteines in human variable heavy-chain domains (V_Hs) at the repertoire level. We identify remarkable usage of non-canonical cysteines within the heavy-chain complementarity-determining region 3 (CDR-H3) and other CDRs and framework regions. Furthermore, our study reveals the diversity and location of non-canonical cysteines and their associated motifs in human V_Hs, which are reminiscent of and more complex than those found in other non-human species such as chicken, camel, llama, shark, and cow. These results explain how non-canonical cysteines strategically occur in the human antibodyome to expand its paratope space. This study will guide the design of human antibodies harboring disulfide-stabilized long CDR-H3s to access difficult-to-target epitopes and influence a paradigm shift in developability involving non-canonical cysteines.

•

NGS-based non-canonical cysteine landscape in human V_Hs

•

1 to 8 non-canonical cysteines and up to 30% in long CDR-H3s

•

An array of potential disulfide motifs adds paratope diversity

•

Non-canonical cysteines in human V_Hs are reminiscent of lower animals

Functional convergence of a germline-encoded neutralizing antibody response in rhesus macaques immunized with HCV envelope glycoproteins

Human IGHV1-69-encoded broadly neutralizing antibodies (bnAbs) that target the hepatitis C virus (HCV) envelope glycoprotein (Env) E2 are important for protection against HCV infection. An IGHV1-69 ortholog gene, VH1.36, is preferentially used for bnAbs isolated from HCV Env-immunized rhesus macaques (RMs). Here, we studied the genetic, structural, and functional properties of VH1.36-encoded bnAbs generated by vaccination, in comparison to IGHV1-69-encoded bnAbs from HCV patients. Global B cell repertoire analysis confirmed the expansion of VH1.36-derived B cells in immunized animals. Most E2-specific, VH1.36-encoded antibodies cross-neutralized HCV. Crystal structures of two RM bnAbs with E2 revealed that the RM bnAbs engaged conserved E2 epitopes using similar molecular features as human bnAbs but with a different binding mode. Longitudinal analyses of the RM antibody repertoire responses during immunization indicated rapid lineage development of VH1.36-encoded bnAbs with limited somatic hypermutation. Our findings suggest functional convergence of a germline-encoded bnAb response to HCV Env with implications for vaccination in humans.

mRNA vaccine-elicited antibodies to SARS-CoV-2 and circulating variants

Here we report on the antibody and memory B cell responses of a cohort of 20 volunteers who received the Moderna (mRNA-1273) or Pfizer-BioNTech (BNT162b2) vaccine against SARS-CoV-21-4. Eight weeks after the second injection of vaccine, volunteers showed high levels of IgM and IgG anti-SARS-CoV-2 spike protein (S) and receptor-binding-domain (RBD) binding titre. Moreover, the plasma neutralizing activity and relative numbers of RBD-specific memory B cells of vaccinated volunteers were equivalent to those of individuals who had recovered from natural infection5,6. However, activity against SARS-CoV-2 variants that encode E484K-, N501Y- or K417N/E484K/N501-mutant S was reduced by a small-but significant-margin. The monoclonal antibodies elicited by the vaccines potently neutralize SARS-CoV-2, and target a number of different RBD epitopes in common with monoclonal antibodies isolated from infected donors5-8. However, neutralization by 14 of the 17 most-potent monoclonal antibodies that we tested was reduced or abolished by the K417N, E484K or N501Y mutation. Notably, these mutations were selected when we cultured recombinant vesicular stomatitis virus expressing SARS-CoV-2 S in the presence of the monoclonal antibodies elicited by the vaccines. Together, these results suggest that the monoclonal antibodies in clinical use should be tested against newly arising variants, and that mRNA vaccines may need to be updated periodically to avoid a potential loss of clinical efficacy.

A single donor is sufficient to produce a highly functional in vitro antibody library

Antibody complementarity determining region diversity has been considered to be the most important metric for the production of a functional antibody library. Generally, the greater the antibody library diversity, the greater the probability of selecting a diverse array of high affinity leads. According to this paradigm, the primary means of elevating library diversity has been by increasing the number of donors. In the present study we explored the possibility of creating an in vitro antibody library from a single healthy individual, showing that the number of lymphocytes, rather than the number of donors, is the key criterion in the production of a diverse and functional antibody library. We describe the construction of a high-quality phage display library comprising 5 × 10⁹ human antibodies by applying an efficient B cell extraction protocol from a single donor and a targeted V-gene amplification strategy favoring specific antibody families for their improved developability profiles. Each step of the library generation process was followed and validated by next generation sequencing to monitor the library quality and diversity. The functionality of the library was tested using several therapeutically relevant targets for which a vast number of different antibodies with desired biophysical properties were obtained.

Antibody Responses in Hepatitis C Infection

Antibody responses in hepatitis C virus (HCV) have been a rather mysterious research topic for many investigators working in the field. Chronic HCV infection is often associated with dysregulation of immune functions particularly in B cells, leading to abnormal lymphoproliferation or the production of autoantibodies that exacerbate inflammation and extrahepatic diseases. When considering the antiviral function of antibody, it was difficult to endorse its role in HCV protection, whereas T-cell response has been shown unequivocally critical for natural recovery. Recent breakthroughs in the study of HCV and antigen-specific antibody responses provide important insights into viral vulnerability to antibodies and the immunogenetic and structural properties of the neutralizing antibodies. The new knowledge reinvigorates HCV vaccine research by illuminating a new path for the rational design of vaccine antigens to elicit broadly neutralizing antibodies for protection.

Evolution of antibody immunity to SARS-CoV-2

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has infected 78 million individuals and is responsible for over 1.7 million deaths to date. Infection is associated with the development of variable levels of antibodies with neutralizing activity, which can protect against infection in animal models^1,2. Antibody levels decrease with time, but, to our knowledge, the nature and quality of the memory B cells that would be required to produce antibodies upon reinfection has not been examined. Here we report on the humoral memory response in a cohort of 87 individuals assessed at 1.3 and 6.2 months after infection with SARS-CoV-2. We find that titres of IgM and IgG antibodies against the receptor-binding domain (RBD) of the spike protein of SARS-CoV-2 decrease significantly over this time period, with IgA being less affected. Concurrently, neutralizing activity in plasma decreases by fivefold in pseudotype virus assays. By contrast, the number of RBD-specific memory B cells remains unchanged at 6.2 months after infection. Memory B cells display clonal turnover after 6.2 months, and the antibodies that they express have greater somatic hypermutation, resistance to RBD mutations and increased potency, indicative of continued evolution of the humoral response. Immunofluorescence and PCR analyses of intestinal biopsies obtained from asymptomatic individuals at 4 months after the onset of coronavirus disease 2019 (COVID-19) revealed the persistence of SARS-CoV-2 nucleic acids and immunoreactivity in the small bowel of 7 out of 14 individuals. We conclude that the memory B cell response to SARS-CoV-2 evolves between 1.3 and 6.2 months after infection in a manner that is consistent with antigen persistence.

A reply to "TCR+/BCR+ dual-expressing cells and their associated public BCR clonotype are not enriched in type 1 diabetes"

We have recently identified a novel lymphocyte that is a dual expresser (DE) of TCRαβ and BCR. DEs in T1D patients are predominated by a public BCR clonotype (clone-x) that encodes a potent autoantigen that cross-activates insulin-reactive T cells. Betts and colleagues were able to detect DEs but alleged to not detect high DE frequency, clone-x, or similar clones in T1D patients. Unfortunately, the authors did not follow our methods and when they did, their flow cytometric data at two sites were conflicting. Moreover, contrary to their claim, we identified clones similar to clone-x in their data along with clones bearing the core motif (DTAMVYYFDYW). Additionally, their report of no increased usage of clone-x VH/DH genes by bulk B cells confirms rather than challenges our results. Finally, the authors failed to provide data verifying purity of their sorted DEs, making it difficult to draw reliable conclusion of their repertoire analysis. This Matters Arising Response paper addresses the Japp et al. (2021) Matters Arising paper, published concurrently in Cell.

mRNA vaccine-elicited antibodies to SARS-CoV-2 and circulating variants

To date severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has infected over 100 million individuals resulting in over two million deaths. Many vaccines are being deployed to prevent coronavirus disease 2019 (COVID-19) including two novel mRNA-based vaccines ^1,2 . These vaccines elicit neutralizing antibodies and appear to be safe and effective, but the precise nature of the elicited antibodies is not known ^3-6 . Here we report on the antibody and memory B cell responses in a cohort of 20 volunteers who received either the Moderna (mRNA-1273) or Pfizer-BioNTech (BNT162b2) vaccines. Consistent with prior reports, 8 weeks after the second vaccine injection volunteers showed high levels of IgM, and IgG anti-SARS-CoV-2 spike protein (S) and receptor binding domain (RBD) binding titers ^3,5,6 . Moreover, the plasma neutralizing activity, and the relative numbers of RBD-specific memory B cells were equivalent to individuals who recovered from natural infection ^7,8 . However, activity against SARS-CoV-2 variants encoding E484K or N501Y or the K417N:E484K:N501Y combination was reduced by a small but significant margin. Consistent with these findings, vaccine-elicited monoclonal antibodies (mAbs) potently neutralize SARS-CoV-2, targeting a number of different RBD epitopes in common with mAbs isolated from infected donors. Structural analyses of mAbs complexed with S trimer suggest that vaccine- and virus-encoded S adopts similar conformations to induce equivalent anti-RBD antibodies. However, neutralization by 14 of the 17 most potent mAbs tested was reduced or abolished by either K417N, or E484K, or N501Y mutations. Notably, the same mutations were selected when recombinant vesicular stomatitis virus (rVSV)/SARS-CoV-2 S was cultured in the presence of the vaccine elicited mAbs. Taken together the results suggest that the monoclonal antibodies in clinical use should be tested against newly arising variants, and that mRNA vaccines may need to be updated periodically to avoid potential loss of clinical efficacy.

Exploiting next-generation sequencing in antibody selections - a simple PCR method to recover binders

Antibody discovery using invitro display technologies such as phage and/or yeast display has become acornerstone in many research and development projects, including the creation of new drugs for clinical use. Traditionally, after the selection phase, random clones are isolated for binding validation and Sanger sequencing. More recently, next-generation sequencing (NGS) technology has allowed deeper insight into the antibody population after aselection campaign, enabling the identification of many more specific binders. However, this approach only provides the DNA sequences of potential binders, the properties of which need to be fully elucidated by obtaining corresponding clones and expressing them for further validation. Here we present arapid novel method to harvest potential clones identified by NGS that uses asimple PCR and yeast recombination approach. The protocol was tested in selections against three different targets and was able to recover clones at an abundance level that would be impractical to identify using traditional methods.

Evolution of Antibody Immunity to SARS-CoV-2

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has infected 78 million individuals and is responsible for over 1.7 million deaths to date. Infection is associated with development of variable levels of antibodies with neutralizing activity that can protect against infection in animal models. Antibody levels decrease with time, but the nature and quality of the memory B cells that would be called upon to produce antibodies upon re-infection has not been examined. Here we report on the humoral memory response in a cohort of 87 individuals assessed at 1.3 and 6.2 months after infection. We find that IgM, and IgG anti-SARS-CoV-2 spike protein receptor binding domain (RBD) antibody titers decrease significantly with IgA being less affected. Concurrently, neutralizing activity in plasma decreases by five-fold in pseudotype virus assays. In contrast, the number of RBD-specific memory B cells is unchanged. Memory B cells display clonal turnover after 6.2 months, and the antibodies they express have greater somatic hypermutation, increased potency and resistance to RBD mutations, indicative of continued evolution of the humoral response. Analysis of intestinal biopsies obtained from asymptomatic individuals 4 months after coronavirus disease-2019 (COVID-19) onset, using immunofluorescence, or polymerase chain reaction, revealed persistence of SARS-CoV-2 nucleic acids and immunoreactivity in the small bowel of 7 out of 14 volunteers. We conclude that the memory B cell response to SARS-CoV-2 evolves between 1.3 and 6.2 months after infection in a manner that is consistent with antigen persistence.

Sato A. Functional GLP-1R antibodies identified from a synthetic GPCR-focused library demonstrate potent blood glucose control

G protein-coupled receptors (GPCRs) are a group of seven-transmembrane receptor proteins that have proven to be successful drug targets. Antibodies are becoming an increasingly promising modality to target these receptors due to their unique properties, such as exquisite specificity, long half-life, and fewer side effects, and their improved pharmacokinetic and pharmacodynamic profiles compared to peptides and small molecules, which results from their more favorable biodistribution. To date, there are only two US Food and Drug Administration-approved GPCR antibody drugs, namely erenumab and mogamulizumab, and this highlights the challenges encountered in identifying functional antibodies against GPCRs. Utilizing Twist's precision DNA writing technologies, we have created a GPCR-focused phage display library with 1 × 1010 diversity. Specifically, we mined endogenous GPCR binding ligand and peptide sequences and incorporated these binding motifs into the heavy chain complementarity-determining region 3 in a synthetic antibody library. Glucagon-like peptide-1 receptor (GLP-1 R) is a class B GPCR that acts as the receptor for the incretin GLP-1, which is released to regulate insulin levels in response to food intake. GLP-1 R agonists have been widely used to increase insulin secretion to lower blood glucose levels for the treatment of type 1 and type 2 diabetes, whereas GLP-1 R antagonists have applications in the treatment of severe hypoglycemia associated with bariatric surgery and hyperinsulinomic hypoglycemia. Here we present the discovery and creation of both antagonistic and agonistic GLP-1 R antibodies by panning this GPCR-focused phage display library on a GLP-1 R-overexpressing Chinese hamster ovary cell line and demonstrate their in vitro and in vivo functional activity.

Antigen receptor repertoires of one of the smallest known vertebrates

The rules underlying the structure of antigen receptor repertoires are not yet fully defined, despite their enormous importance for the understanding of adaptive immunity. With current technology, the large antigen receptor repertoires of mice and humans cannot be comprehensively studied. To circumvent the problems associated with incomplete sampling, we have studied the immunogenetic features of one of the smallest known vertebrates, the cyprinid fish Paedocypris sp. "Singkep" ("minifish"). Despite its small size, minifish has the key genetic facilities characterizing the principal vertebrate lymphocyte lineages. As described for mammals, the frequency distributions of immunoglobulin and T cell receptor clonotypes exhibit the features of fractal systems, demonstrating that self-similarity is a fundamental property of antigen receptor repertoires of vertebrates, irrespective of body size. Hence, minifish achieve immunocompetence via a few thousand lymphocytes organized in robust scale-free networks, thereby ensuring immune reactivity even when cells are lost or clone sizes fluctuate during immune responses.

The Usage of Human IGHJ Genes Follows a Particular Non-random Selection: The Recombination Signal Sequence May Affect the Usage of Human IGHJ Genes

The formation of the B cell receptor (BCR) heavy chain variable region is derived from the germline V(D)J gene rearrangement according to the “12/23” rule and the “beyond 12/23” rule. The usage frequency of each V(D)J gene in the peripheral BCR repertoires is related to the initial recombination, self-tolerance selection, and the clonal proliferative response. However, their specific differences and possible mechanisms are still unknown. We analyzed in-frame and out-of-frame BCR-H repertoires from human samples with normal physiological and various pathological conditions by high-throughput sequencing. Our results showed that IGHJ gene frequency follows a similar pattern which is previously known, where IGHJ4 is used at high frequency (>40%), IGHJ6/IGHJ3/IGHJ5 is used at medium frequencies (10∼20%), and IGH2/IGHJ1 is used at low frequency (<4%) under whether normal physiological or various pathological conditions. However, our analysis of the recombination signal sequences suggested that the conserved non-amer and heptamer and certain 23 bp spacer length may affect the initial IGHD-IGHJ recombination, which results in different frequencies of IGHJ genes among the initial BCR-H repertoire. Based on this “initial repertoire,” we recommend that re-evaluation and further investigation are needed when analyzing the significance and mechanism of IGHJ gene frequency in self-tolerance selection and the clonal proliferative response.