Deep Sequencing of B Cell Receptor Repertoires From COVID-19 Patients Reveals Strong Convergent Immune Signatures

Computational detection of antigen-specific B cell receptors following immunization

B cell receptors (BCRs) play a crucial role in recognizing and fighting foreign antigens. High-throughput sequencing enables in-depth sampling of the BCRs repertoire after immunization. However, only a minor fraction of BCRs actively participate in any given infection. To what extent can we accurately identify antigen-specific sequences directly from BCRs repertoires? We present a computational method grounded on sequence similarity, aimed at identifying statistically significant responsive BCRs. This method leverages well-known characteristics of affinity maturation and expected diversity. We validate its effectiveness using longitudinally sampled human immune repertoire data following influenza vaccination and SARS-CoV-2 infections. We show that different lineages converge to the same responding Complementarity Determining Region 3, demonstrating convergent selection within an individual. The outcomes of this method hold promise for application in vaccine development, personalized medicine, and antibody-derived therapeutics.

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.

Computational mining of B cell receptor repertoires reveals antigen-specific and convergent responses to Ebola vaccination

Outbreaks of Ebolaviruses, such as Sudanvirus (SUDV) in Uganda in 2022, demonstrate that species other than the Zaire ebolavirus (EBOV), which is currently the sole virus represented in current licensed vaccines, remain a major threat to global health. There is a pressing need to develop effective pan-species vaccines and novel monoclonal antibody-based therapeutics for Ebolavirus disease. In response to recent outbreaks, the two dose, heterologous Ad26.ZEBOV/MVA-BN-Filo vaccine regimen was developed and was tested in a large phase II clinical trial (EBL2001) as part of the EBOVAC2 consortium. Here, we perform bulk sequencing of the variable heavy chain (VH) of B cell receptors (BCR) in forty participants from the EBL2001 trial in order to characterize the BCR repertoire in response to vaccination with Ad26.ZEBOV/MVA-BN-Filo. We develop a comprehensive database, EBOV-AbDab, of publicly available Ebolavirus-specific antibody sequences. We then use our database to predict the antigen-specific component of the vaccinee repertoires. Our results show striking convergence in VH germline gene usage across participants following the MVA-BN-Filo dose, and provide further evidence of the role of IGHV3-15 and IGHV3-13 antibodies in the B cell response to Ebolavirus glycoprotein. Furthermore, we found that previously described Ebola-specific mAb sequences present in EBOV-AbDab were sufficient to describe at least one of the ten most expanded BCR clonotypes in more than two thirds of our cohort of vaccinees following the boost, providing proof of principle for the utility of computational mining of immune repertoires.

An unconventional VH1-2 antibody tolerates escape mutations and shows an antigenic hotspot on SARS-CoV-2 spike

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike (S) protein continues to evolve antigenically, impacting antibody immunity. D1F6, an affinity-matured non-stereotypic VH1-2 antibody isolated from a patient infected with the SARS-CoV-2 ancestral strain, effectively neutralizes most Omicron variants tested, including XBB.1.5. We identify that D1F6 in the immunoglobulin G (IgG) form is able to overcome the effect of most Omicron mutations through its avidity-enhanced multivalent S-trimer binding. Cryo-electron microscopy (cryo-EM) and biochemical analyses show that three simultaneous epitope mutations are generally needed to substantially disrupt the multivalent S-trimer binding by D1F6 IgG. Antigenic mutations at spike positions 346, 444, and 445, which appeared in the latest variants, have little effect on D1F6 binding individually. However, these mutations are able to act synergistically with earlier Omicron mutations to impair neutralization by affecting the interaction between D1F6 IgG and the S-trimer. These results provide insight into the mechanism by which accumulated antigenic mutations facilitate evasion of affinity-matured antibodies.

Protocol for fast clonal family inference and analysis from large-scale B cell receptor repertoire sequencing data

The expeditious identification and comprehensive analysis of clonal families from extensive B cell receptor (BCR) repertoire sequencing data are imperative for elucidating the intricacies of B cell immune responses. Here, we introduce a computational pipeline designed to swiftly deduce clonal families from bulk BCR heavy-chain sequencing data, accompanied by a suite of functional modules tailored to streamline post-clustering analysis. The outlined methodology encompasses guidelines for software installation, meticulous data preparation, and the systematic inference and analysis of clonal families. For complete details on the use and execution of this protocol, please refer to Wang et al.1.

Predictability of B cell clonal persistence and immunosurveillance in breast cancer

B cells and T cells are important components of the adaptive immune system and mediate anticancer immunity. The T cell landscape in cancer is well characterized, but the contribution of B cells to anticancer immunosurveillance is less well explored. Here we show an integrative analysis of the B cell and T cell receptor repertoire from individuals with metastatic breast cancer and individuals with early breast cancer during neoadjuvant therapy. Using immune receptor, RNA and whole-exome sequencing, we show that both B cell and T cell responses seem to coevolve with the metastatic cancer genomes and mirror tumor mutational and neoantigen architecture. B cell clones associated with metastatic immunosurveillance and temporal persistence were more expanded and distinct from site-specific clones. B cell clonal immunosurveillance and temporal persistence are predictable from the clonal structure, with higher-centrality B cell antigen receptors more likely to be detected across multiple metastases or across time. This predictability was generalizable across other immune-mediated disorders. This work lays a foundation for prioritizing antibody sequences for therapeutic targeting in cancer.

Innovations and trends in antibody repertoire analysis

Monoclonal antibodies have revolutionized the treatment of human diseases, which has made them the fastest-growing class of therapeutics, with global sales expected to reach $346.6 billion USD by 2028. Advances in antibody engineering and development have led to the creation of increasingly sophisticated antibody-based therapeutics (e.g. bispecific antibodies and chimeric antigen receptor T cells). However, approaches for antibody discovery have remained comparatively grounded in conventional yet reliable in vitro assays. Breakthrough developments in high-throughput single B-cell sequencing and immunoglobulin proteomic serology, however, have enabled the identification of high-affinity antibodies directly from endogenous B cells or circulating immunoglobulin produced in vivo. Moreover, advances in artificial intelligence offer vast potential for antibody discovery and design with large-scale repertoire datasets positioned as the optimal source of training data for such applications. We highlight advances and recent trends in how these technologies are being applied to antibody repertoire analysis.

Benchmarking antibody clustering methods using sequence, structural, and machine learning similarity measures for antibody discovery applications

Antibodies are proteins produced by our immune system that have been harnessed as biotherapeutics. The discovery of antibody-based therapeutics relies on analyzing large volumes of diverse sequences coming from phage display or animal immunizations. Identification of suitable therapeutic candidates is achieved by grouping the sequences by their similarity and subsequent selection of a diverse set of antibodies for further tests. Such groupings are typically created using sequence-similarity measures alone. Maximizing diversity in selected candidates is crucial to reducing the number of tests of molecules with near-identical properties. With the advances in structural modeling and machine learning, antibodies can now be grouped across other diversity dimensions, such as predicted paratopes or three-dimensional structures. Here we benchmarked antibody grouping methods using clonotype, sequence, paratope prediction, structure prediction, and embedding information. The results were benchmarked on two tasks: binder detection and epitope mapping. We demonstrate that on binder detection no method appears to outperform the others, while on epitope mapping, clonotype, paratope, and embedding clusterings are top performers. Most importantly, all the methods propose orthogonal groupings, offering more diverse pools of candidates when using multiple methods than any single method alone. To facilitate exploring the diversity of antibodies using different methods, we have created an online tool-CLAP-available at (clap.naturalantibody.com) that allows users to group, contrast, and visualize antibodies using the different grouping methods.

Guidelines for reproducible analysis of adaptive immune receptor repertoire sequencing data

Enhancing the reproducibility and comprehension of adaptive immune receptor repertoire sequencing (AIRR-seq) data analysis is critical for scientific progress. This study presents guidelines for reproducible AIRR-seq data analysis, and a collection of ready-to-use pipelines with comprehensive documentation. To this end, ten common pipelines were implemented using ViaFoundry, a user-friendly interface for pipeline management and automation. This is accompanied by versioned containers, documentation and archiving capabilities. The automation of pre-processing analysis steps and the ability to modify pipeline parameters according to specific research needs are emphasized. AIRR-seq data analysis is highly sensitive to varying parameters and setups; using the guidelines presented here, the ability to reproduce previously published results is demonstrated. This work promotes transparency, reproducibility, and collaboration in AIRR-seq data analysis, serving as a model for handling and documenting bioinformatics pipelines in other research domains.

Systems biology of B cells in COVID-19

The integration of multi-'omic datasets into complex systems-wide assessments has become a mainstay in immunologic investigation. This focus on high-dimensional data collection and analysis was on full display in the investigation of COVID-19, the respiratory illness resulting from infection by the novel coronavirus SARS-CoV-2. Particularly in the area of B cell biology, tremendous efforts in both cellular and serologic investigation have resulted in an increasingly detailed mapping of the coordinated effector, memory, and antibody secreting cell responses that underpin the development of humoral immunity in response to primary viral infection. Further, the rapid development and deployment of effective vaccines has allowed for the assessment of developing memory responses across a wide variety of immune contexts, including in patients with compromised immune function. The result has been a period of rapid gains in the understanding of B cell biology unrestricted to the study of COVID-19. Here, we outline the systems-level technologies that have been routinely implemented in these investigations throughout the pandemic, and discuss how their use has led to clear and applicable gains in pursuance of the amelioration of human infectious disease and beyond.

The characteristics of BCR-CDR3 repertoire in COVID-19 patients and SARS-CoV-2 vaccinated volunteers

The global COVID-19 pandemic has caused more than 1 billion infections, and numerous SARS-CoV-2 vaccines developed rapidly have been administered over 10 billion doses. The world is continuously concerned about the cytokine storms induced by the interaction between SARS-CoV-2 and host, long COVID, breakthrough infections postvaccination, and the impact of SARS-CoV-2 variants. BCR-CDR3 repertoire serves as a molecular target for monitoring the antiviral response "trace" of B cells, evaluating the effects, mechanisms, and memory abilities of individual responses to B cells, and has been successfully applied in analyzing the infection mechanisms, vaccine improvement, and neutralizing antibodies preparation of influenza virus, HIV, MERS, and Ebola virus. Based on research on BCR-CDR3 repertoire of COVID-19 patients and volunteers who received different SARS-CoV-2 vaccines in multiple laboratories worldwide, we focus on analyzing the characteristics and changes of BCR-CDR3 repertoire, such as diversity, clonality, V&J genes usage and pairing, SHM, CSR, shared CDR3 clones, as well as the summary on BCR sequences targeting virus-specific epitopes in the preparation and application research of SARS-CoV-2 potential therapeutic monoclonal antibodies. This review provides comparative data and new research schemes for studying the possible mechanisms of differences in B cell response between SARS-CoV-2 infection or vaccination, and supplies a foundation for improving vaccines after SARS-CoV-2 mutations and potential antibody therapy for infected individuals.

Comparable CD8+ T-cell responses to SARS-CoV-2 vaccination in single-cell transcriptomics of recently allogeneic transplanted patients and healthy individuals

Despite extensive research on severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) vaccination responses in healthy individuals, there is comparatively little known beyond antibody titers and T-cell responses in the vulnerable cohort of patients after allogeneic hematopoietic stem cell transplantation (ASCT). In this study, we assessed the serological response and performed longitudinal multimodal analyses including T-cell functionality and single-cell RNA sequencing combined with T cell receptor (TCR)/B cell receptor (BCR) profiling in the context of BNT162b2 vaccination in ASCT patients. In addition, these data were compared to publicly available data sets of healthy vaccinees. Protective antibody titers were achieved in 40% of patients. We identified a distorted B- and T-cell distribution, a reduced TCR diversity, and increased levels of exhaustion marker expression as possible causes for the poorer vaccine response rates in ASCT patients. Immunoglobulin heavy chain gene rearrangement after vaccination proved to be highly variable in ASCT patients. Changes in TCRα and TCRβ gene rearrangement after vaccination differed from patterns observed in healthy vaccinees. Crucially, ASCT patients elicited comparable proportions of SARS-CoV-2 vaccine-induced (VI) CD8+ T-cells, characterized by a distinct gene expression pattern that is associated with SARS-CoV-2 specificity in healthy individuals. Our study underlines the impaired immune system and thus the lower vaccine response rates in ASCT patients. However, since protective vaccine responses and VI CD8+ T-cells can be induced in part of ASCT patients, our data advocate early posttransplant vaccination due to the high risk of infection in this vulnerable group.

Building Representation Learning Models for Antibody Comprehension

Antibodies are versatile proteins with both the capacity to bind a broad range of targets and a proven track record as some of the most successful therapeutics. However, the development of novel antibody therapeutics is a lengthy and costly process. It is challenging to predict the functional and biophysical properties of antibodies from their amino acid sequence alone, requiring numerous experiments for full characterization. Machine learning, specifically deep representation learning, has emerged as a family of methods that can complement wet lab approaches and accelerate the overall discovery and engineering process. Here, we review advances in antibody sequence representation learning, and how this has improved antibody structure prediction and facilitated antibody optimization. We discuss challenges in the development and implementation of such models, such as the lack of publicly available, well-curated antibody function data and highlight opportunities for improvement. These and future advances in machine learning for antibody sequences have the potential to increase the success rate in developing new therapeutics, resulting in broader access to transformative medicines and improved patient outcomes.

Characterizing adjuvants' effects at murine immunoglobulin repertoire level

Generating large-scale, high-fidelity sequencing data is challenging and, furthermore, not much has been done to characterize adjuvants' effects at the repertoire level. Thus, we introduced an IgSeq pipeline that standardized library prep protocols and data analysis functions for accurate repertoire profiling. We then studied systemically effects of CpG and Alum on the Ig heavy chain repertoire using the ovalbumin (OVA) murine model. Ig repertoires of different tissues (spleen and bone marrow) and isotypes (IgG and IgM) were examined and compared in IGHV mutation, gene usage, CDR3 length, clonal diversity, and clonal selection. We found Ig repertoires of different compartments exhibited distinguishable profiles at the non-immunized steady state, and distinctions became more pronounced upon adjuvanted immunizations. Notably, Alum and CpG effects exhibited different tissue- and isotype-preferences. The former led to increased diversity of abundant clones in bone marrow, and the latter promoted the selection of IgG clones in both tissues.

IgMAT: immunoglobulin sequence multi-species annotation tool for any species including those with incomplete antibody annotation or unusual characteristics

BACKGROUND

The advent and continual improvement of high-throughput sequencing technologies has made immunoglobulin repertoire sequencing accessible and informative regardless of study species. However, to fully map dynamic changes in polyclonal responses precise framework and complementarity determining region annotation of rearranging genes is pivotal. Most sequence annotation tools are designed primarily for use with human and mouse antibody sequences which use databases with fixed species lists, applying very specific assumptions which select against unique structural characteristics. For this reason, data agnostic tools able to learn from presented data can be very useful with new species or with novel datasets.

RESULTS

We have developed IgMAT, which utilises a reduced amino acid alphabet, that incorporates multiple HMM alignments into a single consensus to automatically annotate immunoglobulin sequences from most organisms. Additionally, the software allows the incorporation of user defined databases to better represent the species and/or antibody class of interest. To demonstrate the accuracy and utility of IgMAT, we present analysis of sequences extracted from structural data and immunoglobulin sequence datasets from several different species.

CONCLUSIONS

IgMAT is fully open-sourced and freely available on GitHub ( https://github.com/TPI-Immunogenetics/igmat ) for download under GPLv3 license. It can be used as a CLI application or as a python module to be integrated in custom scripts.

Interaction dynamics between innate and adaptive immune cells responding to SARS-CoV-2 vaccination in non-human primates

As SARS-CoV-2 variants continue evolving, testing updated vaccines in non-human primates remains important for guiding human clinical practice. To date, such studies have focused on antibody titers and antigen-specific B and T cell frequencies. Here, we extend our understanding by integrating innate and adaptive immune responses to mRNA-1273 vaccination in rhesus macaques. We sorted innate immune cells from a pre-vaccine time point, as well as innate immune cells and antigen-specific peripheral B and T cells two weeks after each of two vaccine doses and used single-cell sequencing to assess the transcriptomes and adaptive immune receptors of each cell. We show that a subset of S-specific T cells expresses cytokines critical for activating innate responses, with a concomitant increase in CCR5-expressing intermediate monocytes and a shift of natural killer cells to a more cytotoxic phenotype. The second vaccine dose, administered 4 weeks after the first, elicits an increase in circulating germinal center-like B cells 2 weeks later, which are more clonally expanded and enriched for epitopes in the receptor binding domain. Both doses stimulate inflammatory response genes associated with elevated antibody production. Overall, we provide a comprehensive picture of bidirectional signaling between innate and adaptive components of the immune system and suggest potential mechanisms for the enhanced response to secondary exposure.

Fast clonal family inference from large-scale B cell repertoire sequencing data

Advances in high-throughput sequencing technologies have facilitated the large-scale characterization of B cell receptor (BCR) repertoires. However, the vast amount and high diversity of the BCR sequences pose challenges for efficient and biologically meaningful analysis. Here, we introduce fastBCR, an efficient computational approach for inferring B cell clonal families from massive BCR heavy chain sequences. We demonstrate that fastBCR substantially reduces the running time while ensuring high accuracy on simulated datasets with diverse numbers of B cell lineages and varying mutation rates. We apply fastBCR to real BCR sequencing data from peripheral blood samples of COVID-19 patients, showing that the inferred clonal families display disease-associated features, as well as corresponding antigen-binding specificity and affinity. Overall, our results demonstrate the advantages of fastBCR for analyzing BCR repertoire data, which will facilitate the identification of disease-associated antibodies and improve our understanding of the B cell immune response.

Broad SARS-CoV-2 neutralization by monoclonal and bispecific antibodies derived from a Gamma-infected individual

The pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has remained a medical threat due to the evolution of multiple variants that acquire resistance to vaccines and prior infection. Therefore, it is imperative to discover monoclonal antibodies (mAbs) that neutralize a broad range of SARS-CoV-2 variants. A stabilized spike glycoprotein was used to enrich antigen-specific B cells from an individual with a primary Gamma variant infection. Five mAbs selected from those B cells showed considerable neutralizing potency against multiple variants, with COVA309-35 being the most potent against the autologous virus, as well as Omicron BA.1 and BA.2, and COVA309-22 having binding and neutralization activity against Omicron BA.4/5, BQ.1.1, and XBB.1. When combining the COVA309 mAbs as cocktails or bispecific antibodies, the breadth and potency were improved. In addition, the mechanism of cross-neutralization of the COVA309 mAbs was elucidated by structural analysis. Altogether these data indicate that a Gamma-infected individual can develop broadly neutralizing antibodies.

Post-hospitalization rehabilitation alleviates long-term immune repertoire alteration in COVID-19 convalescent patients

The global pandemic of Coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an once-in-a-lifetime public health crisis. Among hundreds of millions of people who have contracted with or are being infected with COVID-19, the question of whether COVID-19 infection may cause long-term health concern, even being completely recovered from the disease clinically, especially immune system damage, needs to be addressed. Here, we performed seven-chain adaptome immune repertoire analyses on convalescent COVID-19 patients who have been discharged from hospitals for at least 6 months. Surprisingly, we discovered lymphopenia, reduced number of unique CDR3s, and reduced diversity of the TCR/BCR immune repertoire in convalescent COVID-19 patients. In addition, the BCR repertoire appears to be activated, which is consistent with the protective antibody titres, but serological experiments reveal significantly lower IL-4 and IL-7 levels in convalescent patients compared to those in healthy controls. Finally, in comparison with convalescent patients who did not receive post-hospitalization rehabilitation, the convalescent patients who received post-hospitalization rehabilitation had attenuated immune repertoire abnormality, almost back to the level of healthy control, despite no detectable clinic demographic difference. Overall, we report the potential long-term immunological impairment for COVID-19 infection, and correction of this impairment via post-hospitalization rehabilitation may offer a new prospect for COVID-19 recovery strategy.

Characterization of B-cell receptor clonality and immunoglobulin gene usage at multiple time points during active SARS-CoV-2 infection

Although monoclonal antibodies to the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) are known, B-cell receptor repertoire and its change in patients during coronavirus disease-2019 (COVID-19) progression is underreported. We aimed to study this molecularly. We used immunoglobulin heavy chain (IGH) variable region (IGHV) spectratyping and next-generation sequencing of peripheral blood B-cell genomic DNA collected at multiple time points during disease evolution to study B-cell response to SARS-CoV-2 infection in 14 individuals with acute COVID-19. We found a broad distribution of responding B-cell clones. The IGH gene usage was not significantly skewed but frequencies of individual IGH genes changed repeatedly. We found predominant usage of unmutated and low mutation-loaded IGHV rearrangements characterizing naïve and extrafollicular B cells among the majority of expanded peripheral B-cell clonal lineages at most tested time points in most patients. IGH rearrangement usage showed no apparent relation to anti-SARS-CoV-2 antibody titers. Some patients demonstrated mono/oligoclonal populations carrying highly mutated IGHV rearrangements indicating antigen experience at some of the time points tested, including even before anti-SARS-CoV-2 antibodies were detected. We present evidence demonstrating that the B-cell response to SARS-CoV-2 is individual and includes different lineages of B cells at various time points during COVID-19 progression.

The immunogenicity of human-origin therapeutic antibodies are associated with V gene usage

Therapeutic antibodies can elicit unwanted immune responses in a subset of patients, which leads to the production of anti-drug antibodies (ADA). Some of these ADAs have been reported to effect the pharmacokinetics, efficacy and/or safety of the therapeutic antibodies. The sequence diversity of antibodies are generated by VDJ recombination and mutagenesis. While the antibody generation process can create a large candidate pool for identifying high-affinity antibodies, it also could produce sequences that are foreign to the human immune system. However, it is not clear how VDJ recombination and mutagenesis impact the clinical ADA rate of therapeutic antibodies. In this study, we identified a positive correlation between the clinical ADA rate and the number of introduced mutations in the antibody sequences. We also found that the use of rare V alleles in human-origin antibody therapeutics is associated with higher risk of immunogenicity. The results suggest that antibody engineering projects should start with frameworks that contain commonly used V alleles and prioritize antibody candidates with low number of mutations to reduce the risk of immunogenicity.

Phage display sequencing reveals that genetic, environmental, and intrinsic factors influence variation of human antibody epitope repertoire

Phage-displayed immunoprecipitation sequencing (PhIP-seq) has enabled high-throughput profiling of human antibody repertoires. However, a comprehensive overview of environmental and genetic determinants shaping human adaptive immunity is lacking. In this study, we investigated the effects of genetic, environmental, and intrinsic factors on the variation in human antibody repertoires. We characterized serological antibody repertoires against 344,000 peptides using PhIP-seq libraries from a wide range of microbial and environmental antigens in 1,443 participants from a population cohort. We detected individual-specificity, temporal consistency, and co-housing similarities in antibody repertoires. Genetic analyses showed the involvement of the HLA, IGHV, and FUT2 gene regions in antibody-bound peptide reactivity. Furthermore, we uncovered associations between phenotypic factors (including age, cell counts, sex, smoking behavior, and allergies, among others) and particular antibody-bound peptides. Our results indicate that human antibody epitope repertoires are shaped by both genetics and environmental exposures and highlight specific signatures of distinct phenotypes and genotypes.

The Longitudinal Analysis of Convergent Antibody VDJ Regions in SARS-CoV-2-Positive Patients Using RNA-Seq

Severe acute respiratory syndrome-related coronavirus-2 (SARS-CoV-2) is an ongoing pandemic that continues to evolve and reinfect individuals. To understand the convergent antibody responses that evolved over the course of the pandemic, we evaluated the immunoglobulin repertoire of individuals infected by different SARS-CoV-2 variants for similarity between patients. We utilized four public RNA-seq data sets collected between March 2020 and March 2022 from the Gene Expression Omnibus (GEO) in our longitudinal analysis. This covered individuals infected with Alpha and Omicron variants. In total, from 269 SARS-CoV-2-positive patients and 26 negative patients, 629,133 immunoglobulin heavy-chain variable region V(D)J sequences were reconstructed from sequencing data. We grouped samples based on the SARS-CoV-2 variant type and/or the time they were collected from patients. Our comparison of patients within each SARS-CoV-2-positive group found 1011 common V(D)Js (same V gene, J gene and CDR3 amino acid sequence) shared by more than one patient and no common V(D)Js in the noninfected group. Taking convergence into account, we clustered based on similar CDR3 sequence and identified 129 convergent clusters from the SARS-CoV-2-positive groups. Within the top 15 clusters, 4 contain known anti-SARS-CoV-2 immunoglobulin sequences with 1 cluster confirmed to cross-neutralize variants from Alpha to Omicron. In our analysis of longitudinal groups that include Alpha and Omicron variants, we find that 2.7% of the common CDR3s found within groups were also present in more than one group. Our analysis reveals common and convergent antibodies, which include anti-SARS-CoV-2 antibodies, in patient groups over various stages of the pandemic.

A Rapid Method for Detection of Antigen-Specific B Cells

The global SARS-CoV-2 pandemic has united the efforts of many scientists all over the world to develop wet-lab techniques and computational approaches aimed at the identification of antigen-specific T and B cells. The latter provide specific humoral immunity that is essential for the survival of COVID-19 patients, and vaccine development has essentially been based on these cells. Here, we implemented an approach that integrates the sorting of antigen-specific B cells and B-cell receptor mRNA sequencing (BCR-seq), followed by computational analysis. This rapid and cost-efficient method allowed us to identify antigen-specific B cells in the peripheral blood of patients with severe COVID-19 disease. Subsequently, specific BCRs were extracted, cloned, and produced as full antibodies. We confirmed their reactivity toward the spike RBD domain. Such an approach can be effective for the monitoring and identification of B cells participating in an individual immune response.

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

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

The evolutionary and functional significance of germline immunoglobulin gene variation

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

Altered somatic hypermutation patterns in COVID-19 patients classifies disease severity

Introduction

The success of the human body in fighting SARS-CoV2 infection relies on lymphocytes and their antigen receptors. Identifying and characterizing clinically relevant receptors is of utmost importance.

Methods

We report here the application of a machine learning approach, utilizing B cell receptor repertoire sequencing data from severely and mildly infected individuals with SARS-CoV2 compared with uninfected controls.

Results

In contrast to previous studies, our approach successfully stratifies non-infected from infected individuals, as well as disease level of severity. The features that drive this classification are based on somatic hypermutation patterns, and point to alterations in the somatic hypermutation process in COVID-19 patients.

Discussion

These features may be used to build and adapt therapeutic strategies to COVID-19, in particular to quantitatively assess potential diagnostic and therapeutic antibodies. These results constitute a proof of concept for future epidemiological challenges.

Immune repertoire profiling for disease pathobiology

Lymphocytes consist of highly heterogeneous populations, each expressing a specific cell surface receptor corresponding to a particular antigen. Lymphocytes are both the cause and regulator of various diseases, including autoimmune/allergic diseases, lifestyle diseases, neurodegenerative diseases, and cancers. Recently, immune repertoire sequencing has attracted much attention because it helps obtain global profiles of the immune receptor sequences of infiltrating T and B cells in specimens. Immune repertoire sequencing not only helps deepen our understanding of the molecular mechanisms of immune-related pathology but also assists in discovering novel therapeutic modalities for diseases, thereby shedding colorful light on otherwise tiny monotonous cells when observed under a microscope. In this review article, we introduce and detail the background and methodology of immune repertoire sequencing and summarize recent scientific achievements in association with human diseases. Future perspectives on this genetic technique in the field of histopathological research will also be discussed.

B cell receptor repertoire analysis from autopsy samples of COVID-19 patients

Neutralizing antibodies against the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are being developed world over. We investigated the possibility of producing artificial antibodies from the formalin fixation and paraffin-embedding (FFPE) lung lobes of a patient who died by coronavirus disease 2019 (COVID-19). The B-cell receptors repertoire in the lung tissue where SARS-CoV-2 was detected were considered to have highly sensitive virus-neutralizing activity, and artificial antibodies were produced by combining the most frequently detected heavy and light chains. Some neutralizing effects against the SARS-CoV-2 were observed, and mixing two different artificial antibodies had a higher tendency to suppress the virus. The neutralizing effects were similar to the immunoglobulin G obtained from healthy donors who had received a COVID-19 mRNA vaccine. Therefore, the use of FFPE lung tissue, which preserves the condition of direct virus sensitization, to generate artificial antibodies may be useful against future unknown infectious diseases.

simAIRR: simulation of adaptive immune repertoires with realistic receptor sequence sharing for benchmarking of immune state prediction methods

BACKGROUND

Machine learning (ML) has gained significant attention for classifying immune states in adaptive immune receptor repertoires (AIRRs) to support the advancement of immunodiagnostics and therapeutics. Simulated data are crucial for the rigorous benchmarking of AIRR-ML methods. Existing approaches to generating synthetic benchmarking datasets result in the generation of naive repertoires missing the key feature of many shared receptor sequences (selected for common antigens) found in antigen-experienced repertoires.

RESULTS

We demonstrate that a common approach to generating simulated AIRR benchmark datasets can introduce biases, which may be exploited for undesired shortcut learning by certain ML methods. To mitigate undesirable access to true signals in simulated AIRR datasets, we devised a simulation strategy (simAIRR) that constructs antigen-experienced-like repertoires with a realistic overlap of receptor sequences. simAIRR can be used for constructing AIRR-level benchmarks based on a range of assumptions (or experimental data sources) for what constitutes receptor-level immune signals. This includes the possibility of making or not making any prior assumptions regarding the similarity or commonality of immune state-associated sequences that will be used as true signals. We demonstrate the real-world realism of our proposed simulation approach by showing that basic ML strategies perform similarly on simAIRR-generated and real-world experimental AIRR datasets.

CONCLUSIONS

This study sheds light on the potential shortcut learning opportunities for ML methods that can arise with the state-of-the-art way of simulating AIRR datasets. simAIRR is available as a Python package: https://github.com/KanduriC/simAIRR.

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

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

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

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

Advances in antibody discovery from human BCR repertoires

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

Discovery and characterization of SARS-CoV-2 reactive and neutralizing antibodies from humanized CAMouseHG mice through rapid hybridoma screening and high-throughput single-cell V(D)J sequencing

The coronavirus disease 2019 pandemic has caused more than 532 million infections and 6.3 million deaths to date. The reactive and neutralizing fully human antibodies of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are effective detection tools and therapeutic measures. During SARS-CoV-2 infection, a large number of SARS-CoV-2 reactive and neutralizing antibodies will be produced. Most SARS-CoV-2 reactive and neutralizing fully human antibodies are isolated from human and frequently encoded by convergent heavy-chain variable genes. However, SARS-CoV-2 viruses can mutate rapidly during replication and the resistant variants of neutralizing antibodies easily survive and evade the immune response, especially in the face of such focused antibody responses in humans. Therefore, additional tools are needed to develop different kinds of fully human antibodies to compensate for current deficiency. In this study, we utilized antibody humanized CAMouse^HG mice to develop a rapid antibody discovery method and examine the antibody repertoire of SARS-CoV-2 RBD-reactive hybridoma cells derived from CAMouse^HG mice by using high-throughput single-cell V(D)J sequencing analysis. CAMouse^HG mice were immunized by 28-day rapid immunization method. After electrofusion and semi-solid medium screening on day 12 post-electrofusion, 171 hybridoma clones were generated based on the results of SARS-CoV-2 RBD binding activity assay. A rather obvious preferential usage of IGHV6-1 family was found in these hybridoma clones derived from CAMouse^HG mice, which was significantly different from the antibodies found in patients with COVID-19. After further virus neutralization screening and antibody competition assays, we generated a noncompeting two-antibody cocktail, which showed a potent prophylactic protective efficacy against SARS-CoV-2 in cynomolgus macaques. These results indicate that humanized CAMouse^HG mice not only provide a valuable platform to obtain fully human reactive and neutralizing antibodies but also have a different antibody repertoire from humans. Thus, humanized CAMouse^HG mice can be used as a good complementary tool in discovery of fully human therapeutic and diagnostic antibodies.

Effects of Prior Infection with SARS-CoV-2 on B Cell Receptor Repertoire Response during Vaccination

Understanding the B cell response to SARS-CoV-2 vaccines is a high priority. High-throughput sequencing of the B cell receptor (BCR) repertoire allows for dynamic characterization of B cell response. Here, we sequenced the BCR repertoire of individuals vaccinated by the Pfizer SARS-CoV-2 mRNA vaccine. We compared BCR repertoires of individuals with previous COVID-19 infection (seropositive) to individuals without previous infection (seronegative). We discovered that vaccine-induced expanded IgG clonotypes had shorter heavy-chain complementarity determining region 3 (HCDR3), and for seropositive individuals, these expanded clonotypes had higher somatic hypermutation (SHM) than seronegative individuals. We uncovered shared clonotypes present in multiple individuals, including 28 clonotypes present across all individuals. These 28 shared clonotypes had higher SHM and shorter HCDR3 lengths compared to the rest of the BCR repertoire. Shared clonotypes were present across both serotypes, indicating convergent evolution due to SARS-CoV-2 vaccination independent of prior viral exposure.

Global patterns of antigen receptor repertoire disruption across adaptive immune compartments in COVID-19

Whereas pathogen-specific T and B cells are a primary focus of interest during infectious disease, we have used COVID-19 to ask whether their emergence comes at a cost of broader B cell and T cell repertoire disruption. We applied a genomic DNA-based approach to concurrently study the immunoglobulin-heavy (IGH) and T cell receptor (TCR) β and δ chain loci of 95 individuals. Our approach detected anticipated repertoire focusing for the IGH repertoire, including expansions of clusters of related sequences temporally aligned with SARS-CoV-2-specific seroconversion, and enrichment of some shared SARS-CoV-2-associated sequences. No significant age-related or disease severity-related deficiencies were noted for the IGH repertoire. By contrast, whereas focusing occurred at the TCRβ and TCRδ loci, including some TCRβ sequence-sharing, disruptive repertoire narrowing was almost entirely limited to many patients aged older than 50 y. By temporarily reducing T cell diversity and by risking expansions of nonbeneficial T cells, these traits may constitute an age-related risk factor for COVID-19, including a vulnerability to new variants for which T cells may provide key protection.

COVID-19-specific transcriptomic signature detectable in blood across multiple cohorts

The coronavirus disease 2019 (COVID-19) caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is spreading across the world despite vast global vaccination efforts. Consequently, many studies have looked for potential human host factors and immune mechanisms associated with the disease. However, most studies have focused on comparing COVID-19 patients to healthy controls, while fewer have elucidated the specific host factors distinguishing COVID-19 from other infections. To discover genes specifically related to COVID-19, we reanalyzed transcriptome data from nine independent cohort studies, covering multiple infections, including COVID-19, influenza, seasonal coronaviruses, and bacterial pneumonia. The identified COVID-19-specific signature consisted of 149 genes, involving many signals previously associated with the disease, such as induction of a strong immunoglobulin response and hemostasis, as well as dysregulation of cell cycle-related processes. Additionally, potential new gene candidates related to COVID-19 were discovered. To facilitate exploration of the signature with respect to disease severity, disease progression, and different cell types, we also offer an online tool for easy visualization of the selected genes across multiple datasets at both bulk and single-cell levels.

B-Cell Receptor Repertoire Sequencing: Deeper Digging into the Mechanisms and Clinical Aspects of Immune-mediated Diseases

B cells play an essential role in adaptive immunity and are intimately correlated with pleiotropic immune-mediated diseases. Each B cell occupies a unique B cell receptor (BCR), and all BCRs throughout our body form “BCR repertoire.” With the development of sequencing technology and coupled bioinformatics, accumulating evidence indicates that BCR repertoire largely varies under physiological and pathological conditions. Therefore, comprehensive grasp of BCR repertoire will provide new insights into the pathogenesis of immune-mediated diseases and help exploit efficient diagnostic and treatment strategies. In this review, we start with an overview of BCR repertoire and related sequencing technologies and summarize their current applications in immune-mediated diseases. We also underscore the challenges of this emerging field and propose promising future directions in advancing BCR repertoire exploration.

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

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

Machine-designed biotherapeutics: opportunities, feasibility and advantages of deep learning in computational antibody discovery

Antibodies are versatile molecular binders with an established and growing role as therapeutics. Computational approaches to developing and designing these molecules are being increasingly used to complement traditional lab-based processes. Nowadays, in silico methods fill multiple elements of the discovery stage, such as characterizing antibody-antigen interactions and identifying developability liabilities. Recently, computational methods tackling such problems have begun to follow machine learning paradigms, in many cases deep learning specifically. This paradigm shift offers improvements in established areas such as structure or binding prediction and opens up new possibilities such as language-based modeling of antibody repertoires or machine-learning-based generation of novel sequences. In this review, we critically examine the recent developments in (deep) machine learning approaches to therapeutic antibody design with implications for fully computational antibody design.

Acute Surge of Atypical Memory and Plasma B-Cell Subsets Driven by an Extrafollicular Response in Severe COVID-19

Background

Despite the use of vaccines and therapeutics against the coronavirus disease 2019 (COVID-19) pandemic, this severe disease has been a critical burden on public health, whereas the pathogenic mechanism remains elusive. Recently, accumulating evidence underscores the potential role of the aberrant B-cell response and humoral immunity in disease progression, especially in high-risk groups.

Methods

Using single-cell RNA (scRNA) sequencing analysis, we investigated transcriptional features of B-cell population in peripheral blood from COVID-19 patients and compared them, according to clinical severity and disease course, against a public B-cell dataset.

Results

We confirmed that acute B cells differentiate into plasma cells, particularly in severe patients, potentially through enhanced extrafollicular (EF) differentiation. In severe groups, the elevated plasma B-cell response displayed increased B-cell receptor (BCR) diversity, as well as higher levels of anti–severe acute respiratory syndrome coronavirus 2 (anti–SARS-CoV-2) spike antibodies in plasma, than those in moderate cases, suggesting more robust and heterogeneous plasma cell response in severe COVID-19 patients. Trajectory analysis identified a differentiation pathway for the EF B-cell response from active naïve to atypical memory B cells (AM2), in addition to the emergence of an aberrant plasma cell subset (PC2), which was associated with COVID-19 progression and severity. The AM2 and PC2 subsets surged in the acute phase of the severe disease and presented multiple inflammatory features, including higher cytokine expression and humoral effector function, respectively. These features differ from other B-cell subsets, suggesting a pathogenic potential for disease progression.

Conclusion

The acute surge of AM2 and PC2 subsets with lower somatic hypermutation and higher inflammatory features may be driven by the EF B-cell response during the acute phase of severe COVID-19 and may represent one of the critical drivers in disease severity.

Development and application of a next-generation sequencing protocol and bioinformatics pipeline for the comprehensive analysis of the canine immunoglobulin repertoire

Profiling the adaptive immune repertoire using next generation sequencing (NGS) has become common in human medicine, showing promise in characterizing clonal expansion of B cell clones through analysis of B cell receptors (BCRs) in patients with lymphoid malignancies. In contrast, most work evaluating BCR repertoires in dogs has employed traditional PCR-based approaches analyzing the IGH locus only. The objectives of this study were to: (1) describe a novel NGS protocol to evaluate canine BCRs; (2) develop a bioinformatics pipeline for processing canine BCR sequencing data; and (3) apply these methods to derive insights into BCR repertoires of healthy dogs and dogs undergoing treatment for B-cell lymphoma. RNA from peripheral blood mononuclear cells of healthy dogs (n = 25) and dogs newly diagnosed with intermediate-to-large B-cell lymphoma (n = 18) with intent to pursue chemotherapy was isolated, converted into cDNA and sequenced by NGS. The BCR repertoires were identified and quantified using a novel analysis pipeline. The IGK repertoires of the healthy dogs were far less diverse compared to IGL which, as with IGH, was highly diverse. Strong biases at key positions within the CDR3 sequence were identified within the healthy dog BCR repertoire. For a subset of the dogs with B-cell lymphoma, clonal expansion of specific IGH sequences pre-treatment and reduction post-treatment was observed. The degree of expansion and reduction correlated with the clinical outcome in this subset. Future studies employing these techniques may improve disease monitoring, provide earlier recognition of disease progression, and ultimately lead to more targeted therapeutics.

Landscape of B Cell Receptor Repertoires in COVID-19 Patients Revealed Through CDR3 Sequencing of Immunoglobulin Heavy and Light Chains

The outbreak and persistence of coronavirus disease 2019 (COVID-19) threaten human health. B cells play a vital role in fighting the infections caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Despite many studies on the immune responses in COVID-19 patients, it is still unclear how B cell receptor (BCR) constituents, including immunoglobulin heavy (IGHs) and light chains (IGLs), respond to SARS-CoV-2 in patients with varying symptoms. In this study, we conducted complementarity-determining region 3 (CDR3) sequencing of BCR IGHs and IGLs from the peripheral blood of COVID-19 patients and healthy donors. The results showed significantly reduced clonal diversity, more expanded clones, and longer CDR3 lengths of IGH and IGL in COVID-19 patients than those in healthy individuals. The IGLs had a much higher percentage of VJ skew usage (47.83% IGLV and 42.86% IGLJ were significantly regulated) than the IGHs (12.09% IGHV and 0% IGHJ) between the healthy individuals and patients, which indicated the importance of BCR light chains. Furthermore, we found a largely expanded IGLV3-25 gene cluster mostly pairing with IGLJ1 and ILGJ2 in COVID-19 patients and a newly identified upregulated IGLJ1 gene and IGLJ2+IGLV13-21 recombination, both of which are potential sources of SARS-CoV-2-targeting antibodies. Our findings on specific immune B-cell signatures associated with COVID-19 have clinical implications for vaccine and biomarker development for disease diagnosis.

Potent human broadly SARS-CoV-2-neutralizing IgA and IgG antibodies effective against Omicron BA.1 and BA.2

Memory B-cell and antibody responses to the SARS-CoV-2 spike protein contribute to long-term immune protection against severe COVID-19, which can also be prevented by antibody-based interventions. Here, wide SARS-CoV-2 immunoprofiling in Wuhan COVID-19 convalescents combining serological, cellular, and monoclonal antibody explorations revealed humoral immunity coordination. Detailed characterization of a hundred SARS-CoV-2 spike memory B-cell monoclonal antibodies uncovered diversity in their repertoire and antiviral functions. The latter were influenced by the targeted spike region with strong Fc-dependent effectors to the S2 subunit and potent neutralizers to the receptor-binding domain. Amongst those, Cv2.1169 and Cv2.3194 antibodies cross-neutralized SARS-CoV-2 variants of concern, including Omicron BA.1 and BA.2. Cv2.1169, isolated from a mucosa-derived IgA memory B cell demonstrated potency boost as IgA dimers and therapeutic efficacy as IgG antibodies in animal models. Structural data provided mechanistic clues to Cv2.1169 potency and breadth. Thus, potent broadly neutralizing IgA antibodies elicited in mucosal tissues can stem SARS-CoV-2 infection, and Cv2.1169 and Cv2.3194 are prime candidates for COVID-19 prevention and treatment.

A Deep Learning Model for Accurate Diagnosis of Infection Using Antibody Repertoires

The adaptive immune receptor repertoire consists of the entire set of an individual's BCRs and TCRs and is believed to contain a record of prior immune responses and the potential for future immunity. Analyses of TCR repertoires via deep learning (DL) methods have successfully diagnosed cancers and infectious diseases, including coronavirus disease 2019. However, few studies have used DL to analyze BCR repertoires. In this study, we collected IgG H chain Ab repertoires from 276 healthy control subjects and 326 patients with various infections. We then extracted a comprehensive feature set consisting of 10 subsets of repertoire-level features and 160 sequence-level features and tested whether these features can distinguish between infected individuals and healthy control subjects. Finally, we developed an ensemble DL model, namely, DL method for infection diagnosis (https://github.com/chenyuan0510/DeepID), and used this model to differentiate between the infected and healthy individuals. Four subsets of repertoire-level features and four sequence-level features were selected because of their excellent predictive performance. The DL method for infection diagnosis outperformed traditional machine learning methods in distinguishing between healthy and infected samples (area under the curve = 0.9883) and achieved a multiclassification accuracy of 0.9104. We also observed differences between the healthy and infected groups in V genes usage, clonal expansion, the complexity of reads within clone, the physical properties in the α region, and the local flexibility of the CDR3 amino acid sequence. Our results suggest that the Ab repertoire is a promising biomarker for the diagnosis of various infections.

Antibody repertoire sequencing analysis

High-throughput sequencing for B cell receptor (BCR) repertoire provides useful insights for the adaptive immune system. With the continuous development of the BCR-seq technology, many efforts have been made to develop methods for analyzing the ever-increasing BCR repertoire data. In this review, we comprehensively outline different BCR repertoire library preparation protocols and summarize three major steps of BCR-seq data analysis, i. e., V(D)J sequence annotation, clonal phylogenetic inference, and BCR repertoire profiling and mining. Different from other reviews in this field, we emphasize background intuition and the statistical principle of each method to help biologists better understand it. Finally, we discuss data mining problems for BCR-seq data and with a highlight on recently emerging multiple-sample analysis.

Highly multiplexed immune repertoire sequencing links multiple lymphocyte classes with severity of response to COVID-19

BACKGROUND

Disease progression of subjects with coronavirus disease 2019 (COVID-19) varies dramatically. Understanding the various types of immune response to SARS-CoV-2 is critical for better clinical management of coronavirus outbreaks and to potentially improve future therapies. Disease dynamics can be characterized by deciphering the adaptive immune response.

METHODS

In this cross-sectional study we analyzed 117 peripheral blood immune repertoires from healthy controls and subjects with mild to severe COVID-19 disease to elucidate the interplay between B and T cells. We used an immune repertoire Primer Extension Target Enrichment method (immunoPETE) to sequence simultaneously human leukocyte antigen (HLA) restricted T cell receptor beta chain (TRB) and unrestricted T cell receptor delta chain (TRD) and immunoglobulin heavy chain (IgH) immune receptor repertoires. The distribution was analyzed of TRB, TRD and IgH clones between healthy and COVID-19 infected subjects. Using McFadden's Adjusted R2 variables were examined for a predictive model. The aim of this study is to analyze the influence of the adaptive immune repertoire on the severity of the disease (value on the World Health Organization Clinical Progression Scale) in COVID-19.

FINDINGS

Combining clinical metadata with clonotypes of three immune receptor heavy chains (TRB, TRD, and IgH), we found significant associations between COVID-19 disease severity groups and immune receptor sequences of B and T cell compartments. Logistic regression showed an increase in shared IgH clonal types and decrease of TRD in subjects with severe COVID-19. The probability of finding shared clones of TRD clonal types was highest in healthy subjects (controls). Some specific TRB clones seems to be present in severe COVID-19 (Figure S7b). The most informative models (McFadden´s Adjusted R2=0.141) linked disease severity with immune repertoire measures across all three cell types, as well as receptor-specific cell counts, highlighting the importance of multiple lymphocyte classes in disease progression.

INTERPRETATION

Adaptive immune receptor peripheral blood repertoire measures are associated with COVID-19 disease severity.

FUNDING

The study was funded with grants from the Berlin Institute of Health (BIH).

Rapid Hypermutation B Cell Trajectory Recruits Previously Primed B Cells Upon Third SARS-Cov-2 mRNA Vaccination

The COVID-19 pandemic shows that vaccination strategies building on an ancestral viral strain need to be optimized for the control of potentially emerging viral variants. Therefore, aiming at strong B cell somatic hypermutation to increase antibody affinity to the ancestral strain - not only at high antibody titers - is a priority when utilizing vaccines that are not targeted at individual variants since high affinity may offer some flexibility to compensate for strain-individual mutations. Here, we developed a next-generation sequencing based SARS-CoV-2 B cell tracking protocol to rapidly determine the level of immunoglobulin somatic hypermutation at distinct points during the immunization period. The percentage of somatically hypermutated B cells in the SARS-CoV-2 specific repertoire was low after the primary vaccination series, evolved further over months and increased steeply after boosting. The third vaccination mobilized not only naïve, but also antigen-experienced B cell clones into further rapid somatic hypermutation trajectories indicating increased affinity. Together, the strongly mutated post-booster repertoires and antibodies deriving from this may explain why the third, but not the primary vaccination series, offers some protection against immune-escape variants such as Omicron B.1.1.529.