Language model-based B cell receptor sequence embeddings can effectively encode receptor specificity.

High throughput sequencing of B cell receptors (BCRs) is increasingly applied to study the immense diversity of antibodies. Learning biologically meaningful embeddings of BCR sequences is beneficial for predictive modeling. Several embedding methods have been developed for BCRs, but no direct performance benchmarking exists. Moreover, the impact of the input sequence length and paired-chain information on the prediction remains to be explored. We evaluated the performance of multiple embedding models to predict BCR sequence properties and receptor specificity. Despite the differences in model architectures, most embeddings effectively capture BCR sequence properties and specificity. BCR-specific embeddings slightly outperform general protein language models in predicting specificity. In addition, incorporating full-length heavy chains and paired light chain sequences improves the prediction performance of all embeddings. This study provides insights into the properties of BCR embeddings to improve downstream prediction applications for antibody analysis and discovery.

Immunoglobulin M perception by FcµR.

Immunoglobulin M (IgM) is the first antibody to emerge during embryonic development and the humoral immune response1. IgM can exist in several distinct forms, including monomeric, membrane-bound IgM within the B cell receptor (BCR) complex, pentameric and hexameric IgM in serum and secretory IgM on the mucosal surface. FcµR, the only IgM-specific receptor in mammals, recognizes different forms of IgM to regulate diverse immune responses2-5. However, the underlying molecular mechanisms remain unknown. Here we delineate the structural basis of the FcµR-IgM interaction by crystallography and cryo-electron microscopy. We show that two FcµR molecules interact with a Fcµ-Cµ4 dimer, suggesting that FcµR can bind to membrane-bound IgM with a 2:1 stoichiometry. Further analyses reveal that FcµR-binding sites are accessible in the context of IgM BCR. By contrast, pentameric IgM can recruit four FcµR molecules to bind on the same side and thereby facilitate the formation of an FcµR oligomer. One of these FcµR molecules occupies the binding site of the secretory component. Nevertheless, four FcµR molecules bind to the other side of secretory component-containing secretory IgM, consistent with the function of FcµR in the retrotransport of secretory IgM. These results reveal intricate mechanisms of IgM perception by FcµR.

Structural principles of B cell antigen receptor assembly.

The B cell antigen receptor (BCR) is composed of a membrane-bound class M, D, G, E or A immunoglobulin for antigen recognition1-3 and a disulfide-linked Igα (also known as CD79A) and Igß (also known as CD79B) heterodimer (Igα/ß) that functions as the signalling entity through intracellular immunoreceptor tyrosine-based activation motifs (ITAMs)4,5. The organizing principle of the BCR remains unknown. Here we report cryo-electron microscopy structures of mouse full-length IgM BCR and its Fab-deleted form. At the ectodomain (ECD), the Igα/ß heterodimer mainly uses Igα to associate with Cµ3 and Cµ4 domains of one heavy chain (µHC) while leaving the other heavy chain (µHC') unbound. The transmembrane domain (TMD) helices of µHC and µHC' interact with those of the Igα/ß heterodimer to form a tight four-helix bundle. The asymmetry at the TMD prevents the recruitment of two Igα/ß heterodimers. Notably, the connecting peptide between the ECD and TMD of µHC intervenes in between those of Igα and Igß to guide TMD assembly through charge complementarity. Weaker but distinct density for the Igß ITAM nestles next to the TMD, suggesting potential autoinhibition of ITAM phosphorylation. Interfacial analyses suggest that all BCR classes utilize a general organizational architecture. Our studies provide a structural platform for understanding B cell signalling and designing rational therapies against BCR-mediated diseases.

Characterization of interactions within the Igα/Igß transmembrane domains of the human B-cell receptor provides insights into receptor assembly.

The B-cell receptor (BCR), a complex comprised of a membrane-associated immunoglobulin and the Igα/ß heterodimer, is one of the most important immune receptors in humans and controls B-cell development, activity, selection, and death. BCR signaling plays key roles in autoimmune diseases and lymphoproliferative disorders, yet, despite the clinical significance of this protein complex, key regions (i.e., the transmembrane domains) have yet to be structurally characterized. The mechanism for BCR signaling also remains unclear and has been variously described by the mutually exclusive cross-linking and dissociation activation models. Common to these models is the significance of local plasma membrane composition, which implies that interactions between BCR transmembrane domains (TMDs) play a role in receptor functionality. Here we used an in vivo assay of TMD oligomerization called GALLEX alongside spectroscopic and computational methods to characterize the structures and interactions of human Igα and Igß TMDs in detergent micelles and natural membranes. We observed weak self-association of the Igß TMD and strong self-association of the Igα TMD, which scanning mutagenesis revealed was entirely stabilized by an E-X10-P motif. We also demonstrated strong heterotypic interactions between the Igα and Igß TMDs both in vitro and in vivo, which scanning mutagenesis and computational models suggest is multiconfigurational but can accommodate distinct interaction sites for self-interactions and heterotypic interactions of the Igα TMD. Taken together, these results demonstrate that the TMDs of the human BCR are sites of strong protein-protein interactions that may direct BCR assembly, endoplasmic reticulum retention, and immune signaling.

Detection, sorting, and immortalization of dual expresser lymphocytes from human peripheral blood samples.

This protocol describes how to identify Dual Expressers (DEs), a rare type of lymphocyte that co-expresses B-cell receptors and T-cell receptors, by flow cytometry using a cocktail of four antibodies. It also shows the subsequent gating strategy for detecting and sorting DEs and the generation of EBV-immortalized DE lymphoblastoid cell lines and clones for antibody production and cloning antigen receptors. Use of this protocol maximizes detection of DEs and minimizes inclusion of doublets. For complete details on the use and execution of this protocol, please refer to Ahmed et al. (2019).

Potent neutralization of SARS-CoV-2 variants of concern by an antibody with an uncommon genetic signature and structural mode of spike recognition.

The emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) lineages that are more transmissible and resistant to currently approved antibody therapies poses a considerable challenge to the clinical treatment of coronavirus disease (COVID-19). Therefore, the need for ongoing discovery efforts to identify broadly reactive monoclonal antibodies to SARS-CoV-2 is of utmost importance. Here, we report a panel of SARS-CoV-2 antibodies isolated using the linking B cell receptor to antigen specificity through sequencing (LIBRA-seq) technology from an individual who recovered from COVID-19. Of these antibodies, 54042-4 shows potent neutralization against authentic SARS-CoV-2 viruses, including variants of concern (VOCs). A cryoelectron microscopy (cryo-EM) structure of 54042-4 in complex with the SARS-CoV-2 spike reveals an epitope composed of residues that are highly conserved in currently circulating SARS-CoV-2 lineages. Further, 54042-4 possesses uncommon genetic and structural characteristics that distinguish it from other potently neutralizing SARS-CoV-2 antibodies. Together, these findings provide motivation for the development of 54042-4 as a lead candidate to counteract current and future SARS-CoV-2 VOCs.

The molecular characterization of antibody binding to a superantigen-like protein from a commensal microbe.

Microorganisms have coevolved diverse mechanisms to impair host defenses. A major one, superantigens, can result in devastating effects on the immune system. While all known superantigens induce vast immune cell proliferation and come from opportunistic pathogens, recently, proteins with similar broad specificity to antibody variable (V) domain families were identified in a commensal microbiota. These proteins, identified in the human commensal Ruminococcus gnavus, are called immunoglobulin-binding protein (Ibp) A and B and have been shown to activate B cells in vitro expressing either human VH3 or murine VH5/6/7. Here, we provide molecular and functional studies revealing the basis of this Ibp/immunoglobulin (Ig) interaction. The crystal structure and biochemical assays of a truncated IbpA construct in complex with mouse VH5 antigen-binding fragment (Fab) shows a binding of Ig heavy chain framework residues to the Ibp Domain D and the C-terminal heavy chain binding domain (HCBD). We used targeted mutagenesis of contact residues and affinity measurements and performed studies of the Fab-IbpA complex to determine the stoichiometry between Ibp and VH domains, suggesting Ibp may serve to cluster full-length IgA antibodies in vivo. Furthermore, in vitro stimulation experiments indicate that binding of the Ibp HCBD alone is sufficient to activate responsive murine B cell receptors. The presence of these proteins in a commensal microbe suggest that binding a broad repertoire of immunoglobulins, particularly in the gut/microbiome environment, may provide an important function in the maintenance of host/microbiome homeostasis contrasting with the pathogenic role of structurally homologous superantigens expressed by pathogens.

Self-assembling synthetic nanoadjuvant scaffolds cross-link B cell receptors and represent new platform technology for therapeutic antibody production.

Host antibody responses are pivotal for providing protection against infectious agents. We have pioneered a new class of self-assembling micelles based on pentablock copolymers that enhance antibody responses while providing a low inflammatory environment compared to traditional adjuvants. This type of "just-right" immune response is critical in the rational design of vaccines for older adults. Here, we report on the mechanism of enhancement of antibody responses by pentablock copolymer micelles, which act as scaffolds for antigen presentation to B cells and cross-link B cell receptors, unlike other micelle-forming synthetic block copolymers. We exploited this unique mechanism and developed these scaffolds as a platform technology to produce antibodies in vitro. We show that this novel approach can be used to generate laboratory-scale quantities of therapeutic antibodies against multiple antigens, including those associated with SARS-CoV-2 and Yersinia pestis, further expanding the value of these nanomaterials to rapidly develop countermeasures against infectious diseases.

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

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

Cryo-EM structure of the B cell co-receptor CD19 bound to the tetraspanin CD81.

Signaling through the CD19-CD81 co-receptor complex, in combination with the B cell receptor, is a critical determinant of B cell development and activation. It is unknown how CD81 engages CD19 to enable co-receptor function. Here, we report a 3.8-angstrom structure of the CD19-CD81 complex bound to a therapeutic antigen-binding fragment, determined by cryo-electron microscopy (cryo-EM). The structure includes both the extracellular domains and the transmembrane helices of the complex, revealing a contact interface between the ectodomains that drives complex formation. Upon binding to CD19, CD81 opens its ectodomain to expose a hydrophobic CD19-binding surface and reorganizes its transmembrane helices to occlude a cholesterol binding pocket present in the apoprotein. Our data reveal the structural basis for CD19-CD81 complex assembly, providing a foundation for rational design of therapies for B cell dysfunction.

Use of machine learning to identify a T cell response to SARS-CoV-2.

The identification of SARS-CoV-2-specific T cell receptor (TCR) sequences is critical for understanding T cell responses to SARS-CoV-2. Accordingly, we reanalyze publicly available data from SARS-CoV-2-recovered patients who had low-severity disease (n = 17) and SARS-CoV-2 infection-naive (control) individuals (n = 39). Applying a machine learning approach to TCR beta (TRB) repertoire data, we can classify patient/control samples with a training sensitivity, specificity, and accuracy of 88.2%, 100%, and 96.4% and a testing sensitivity, specificity, and accuracy of 82.4%, 97.4%, and 92.9%, respectively. Interestingly, the same machine learning approach cannot separate SARS-CoV-2 recovered from SARS-CoV-2 infection-naive individual samples on the basis of B cell receptor (immunoglobulin heavy chain; IGH) repertoire data, suggesting that the T cell response to SARS-CoV-2 may be more stereotyped and longer lived. Following validation in larger cohorts, our method may be useful in detecting protective immunity acquired through natural infection or in determining the longevity of vaccine-induced immunity.

Higher-order immunoglobulin repertoire restrictions in CLL: the illustrative case of stereotyped subsets 2 and 169.

Chronic lymphocytic leukemia (CLL) major stereotyped subset 2 (IGHV3-21/IGLV3-21, ∼2.5% of all cases of CLL) is an aggressive disease variant, irrespective of the somatic hypermutation (SHM) status of the clonotypic IGHV gene. Minor stereotyped subset 169 (IGHV3-48/IGLV3-21, ∼0.2% of all cases of CLL) is related to subset 2, as it displays a highly similar variable antigen-binding site. We further explored this relationship through next-generation sequencing and crystallographic analysis of the clonotypic B-cell receptor immunoglobulin. Branching evolution of the predominant clonotype through intraclonal diversification in the context of ongoing SHM was evident in both heavy and light chain genes of both subsets. Molecular similarities between the 2 subsets were highlighted by the finding of shared SHMs within both the heavy and light chain genes in all analyzed cases at either the clonal or subclonal level. Particularly noteworthy in this respect was a ubiquitous SHM at the linker region between the variable and the constant domain of the IGLV3-21 light chains, previously reported as critical for immunoglobulin homotypic interactions underlying cell-autonomous signaling capacity. Notably, crystallographic analysis revealed that the IGLV3-21-bearing CLL subset 169 immunoglobulin retains the same geometry and contact residues for the homotypic intermolecular interaction observed in subset 2, including the SHM at the linker region, and, from a molecular standpoint, belong to a common structural mode of autologous recognition. Collectively, our findings document that stereotyped subsets 2 and 169 are very closely related, displaying shared immunoglobulin features that can be explained only in the context of shared functional selection.

Reading the B-cell receptor immunome in chronic lymphocytic leukemia: revelations and applications.

B-Cell receptor (BCR) sequencing has been the force driving many recent advances in chronic lymphocytic leukemia (CLL) research. Here, we discuss the general principles, revelations, and applications of reading the BCR immunome in the context of CLL. First, IGHV mutational status, obtained by measuring the mutational imprint on the IGHV gene of the CLL clonotype, is the cornerstone of CLL risk stratification. Furthermore, the discovery of "BCR-stereotyped" groups of unrelated patients that share not only a highly similar BCR on their leukemic clone, but also certain clinical characteristics has provided insights key to understanding disease ontogeny. Additionally, whereas the BCR repertoire of most CLL patients is characterized by a single dominant rearrangement, next-generation sequencing (NGS) has revealed a rich subclonal landscape in a larger than previously expected proportion of CLL patients. We review the mechanisms underlying these "multiple dominant" cases, including V(D)J-recombination errors, failure of allelic exclusion, intraclonal diversification, and "true" bi- or oligoclonality, and their implications, in detail. Finally, BCR repertoire sequencing can be used for sensitive quantification of minimal residual disease to potentially unprecedented depth. To surmount pitfalls inherent to this approach and develop internationally harmonized protocols, the EuroClonality-NGS Working Group has been established.

Using B cell receptor lineage structures to predict affinity.

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

TRIP - T cell receptor/immunoglobulin profiler.

BACKGROUND: Antigen receptors are characterized by an extreme diversity of specificities, which poses major computational and analytical challenges, particularly in the era of high-throughput immunoprofiling by next generation sequencing (NGS). The T cell Receptor/Immunoglobulin Profiler (TRIP) tool offers the opportunity for an in-depth analysis based on the processing of the output files of the IMGT/HighV-Quest tool, a standard in NGS immunoprofiling, through a number of interoperable modules. These provide detailed information about antigen receptor gene rearrangements, including variable (V), diversity (D) and joining (J) gene usage, CDR3 amino acid and nucleotide composition and clonality of both T cell receptors (TR) and B cell receptor immunoglobulins (BcR IG), and characteristics of the somatic hypermutation within the BcR IG genes. TRIP is a web application implemented in R shiny. RESULTS: Two sets of experiments have been performed in order to evaluate the efficiency and performance of the TRIP tool. The first used a number of synthetic datasets, ranging from 250k to 1M sequences, and established the linear response time of the tool (about 6 h for 1M sequences processed through the entire BcR IG data pipeline). The reproducibility of the tool was tested comparing the results produced by the main TRIP workflow with the results from a previous pipeline used on the Galaxy platform. As expected, no significant differences were noted between the two tools; although the preselection process seems to be stricter within the TRIP pipeline, about 0.1% more rearrangements were filtered out, with no impact on the final results. CONCLUSIONS: TRIP is a software framework that provides analytical services on antigen receptor gene sequence data. It is accurate and contains functions for data wrangling, cleaning, analysis and visualization, enabling the user to build a pipeline tailored to their needs. TRIP is publicly available at https://bio.tools/TRIP_-_T-cell_Receptor_Immunoglobulin_Profiler .

Hydrophobicity and the Physico-Chemical Basis of Immunotolerance.

This study analyzes the primary electrostatic interaction between the membrane-bound B-cell receptor and antigen peptide determinants, and identifies peptide frequency and hydrophobicity as the main factors that govern and shape immunotolerance versus immunoreactivity.

Isolation of Tailor-Made Antibody Fragments from Yeast-Displayed B-Cell Receptor Repertoires by Multiparameter Fluorescence-Activated Cell Sorting.

In the past decades, monoclonal antibodies have made an unprecedented transformation from research tools to a rapidly growing class of therapeutics. Advancements in the yeast surface display platform enable the selection of favorable mouse or human antibody variants from large B-cell receptor (BCR) gene repertoires that are derived from immunized normal or transgenic animals. Application of high-throughput fluorescence-activated cell sorting (FACS) screening along with well-chosen selection settings can be utilized to identify variants with desired affinities and predefined epitope binding properties. In the following chapter, we describe in detail a multiparameter screening protocol for the selection of antibody variants from yeast libraries generated from BCR gene repertoires from immunized transgenic rats. The procedure provides guidance for the selection of antigen-specific, high-affinity binding, and species cross-reactive human antibodies with a broad epitope coverage. Essentially, this can accelerate target-specific antibody characterization as multiple desirable antibody features can be easily integrated into the selection procedure. In addition, we provide information on how to validate binding behavior of selected candidates after expression as soluble, full-length IgG molecules.

Isolation of Anti-Hapten Antibodies by Fluorescence-Activated Cell Sorting of Yeast-Displayed B-Cell Receptor Gene Repertoires.

Anti-hapten antibodies are widely used as specific immunochemical detection tools in a variety of clinical and environmental analyses. The sensitivity, however, is limited due to the resulting antibody affinities to the haptens which, in turn, leads to a high demand for specific affinity reagents. A well-established path for the generation of high-affinity antibodies is the immunization of animals with the target antigen. However, the generation of anti-hapten antibodies via immunization remains challenging as small molecule haptens usually possess low immunogenicity and, therefore, must be coupled to an immunogenic and high molecular weight carrier to provoke an immune response.Consequently, antibodies are primarily raised against the carrier molecule or structural features of the hapten-linker fused to the carrier protein. This turns the generation of antibodies which bind exclusively to the hapten structure into a search for the needle in a haystack. In the following chapter, we describe how yeast surface display and high-throughput fluorescence-activated cell sorting can be used to isolate anti-hapten antibodies from a large, yeast-displayed B-cell receptor gene library derived from immunized animals. For this, we describe in detail the preparation of protein-hapten conjugates, the immunization procedure, and the subsequent screening process. Moreover, we provide a simple flow cytometry protocol that allows for a rapid analysis of the enriched clones toward free hapten binding.

High-Throughput Mapping of B Cell Receptor Sequences to Antigen Specificity.

B cell receptor (BCR) sequencing is a powerful tool for interrogating immune responses to infection and vaccination, but it provides limited information about the antigen specificity of the sequenced BCRs. Here, we present LIBRA-seq (linking B cell receptor to antigen specificity through sequencing), a technology for high-throughput mapping of paired heavy- and light-chain BCR sequences to their cognate antigen specificities. B cells are mixed with a panel of DNA-barcoded antigens so that both the antigen barcode(s) and BCR sequence are recovered via single-cell next-generation sequencing. Using LIBRA-seq, we mapped the antigen specificity of thousands of B cells from two HIV-infected subjects. The predicted specificities were confirmed for a number of HIV- and influenza-specific antibodies, including known and novel broadly neutralizing antibodies. LIBRA-seq will be an integral tool for antibody discovery and vaccine development efforts against a wide range of antigen targets.

Ionic protein-lipid interactions at the plasma membrane regulate the structure and function of immunoreceptors.

Adaptive lymphocytes express a panel of immunoreceptors on the cell surface. Phospholipids are the major components of cell membranes, but they have functional roles beyond forming lipid bilayers. In particular, acidic phospholipids forming microdomains in the plasma membrane can ionically interact with proteins via polybasic sequences, which can have functional consequences for the protein. We have shown that negatively charged acidic phospholipids can interact with positively charged juxtamembrane polybasic regions of immunoreceptors, such as TCR-CD3, CD28 and IgG-BCR, to regulate protein structure and function. Furthermore, we pay our attention to protein transmembrane domains. We show that a membrane-snorkeling Lys residue in integrin αLß2 regulates transmembrane heterodimer formation and integrin adhesion through ionic interplay with acidic phospholipids and calcium ions (Ca2+) in T cells, thus providing a new mechanism of integrin activation. Here, we review our recent progress showcasing the importance of both juxtamembrane and intramembrane ionic protein-lipid interactions.