Dynamics of B cell repertoires and emergence of cross-reactive responses in patients with different severities of COVID-19

NK- and T-cell repertoire is established early after allogeneic HSCT and is imprinted by CMV reactivation

ABSTRACT

Besides genetic influences, nongenetic factors such as graft-versus-host disease and viral infections have been shown to be important shapers of the immune reconstitution and diversification processes after hematopoietic stem cell transplantation (HSCT). However, differential susceptibility to immune modulation by nongenetic factors is not fully understood. We determined to follow the reconstitution of the T-cell receptor (TCR) repertoire through immune sequencing of natural killer (NK) cells using a 35-marker spectral flow cytometry panel and in relation to clinical events. A longitudinal investigation was performed on samples derived from 54 HSCT recipients during the first year after HSCT. We confirmed a significant contraction in TCR repertoire diversity, with remarkable stability over time. Cytomegalovirus (CMV) reactivation had the ability to significantly change TCR repertoire clonality and composition, with a long-lasting imprint. Our data further revealed skewing of NK-cell reconstitution in CMV reactivated recipients, with an increased frequency of KIR2DL2L3S2+ adaptive, cytolytic, and functional CD107a+ NK cells, concomitant with a reduced pool of NKG2A+ NK cells. We provided support that CMV might act as an important driver of peripheral homeostatic proliferation of circulating specific T and NK cells, which can be viewed as a compensatory mechanism to establish a new peripheral repertoire.

Diverse BCR usage and T cell activation induced by different COVID-19 sequential vaccinations

Limited knowledge is available on the differences in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) specific antibody breadth and T cell differentiation among different COVID-19 sequential vaccination strategies. In this study, we compared the immunogenicity of the third different dose of COVID-19 vaccines, such as mRNA (I-I-M), adenoviral vector (I-I-A), and recombinant protein (I-I-R) vaccines, in terms of the magnitude and breadth of antibody response and differentiation of SARS-CoV-2-specific T and B cells. These studies were performed in the same clinical trial, and the samples were assessed in the same laboratory. IGHV1-69, IGHV3-9, and IGHV4-34 were the dominant B cell receptor (BCR) usages of the I-I-M, I-I-A, and I-I-R groups, respectively; the RBD+ B cell activation capacities were comparable. Additionally, the I-I-R group was characterized by higher numbers of regulatory T cells, circulating T follicular helper cells (cTFH) - cTFH1 (CXRC3+CCR6-), cTFH1-17 (CXRC3+CCR6+), cTFH17 (CXRC3-CCR6+), and cTFH-CM (CD45RA-CCR7+), and lower SMNE+ T cell proliferative capacity than the other two groups, whereas I-I-A showed a higher proportion and number of virus-specific CD4+ T cells than I-I-R, as determined in ex vivo experiments. Our data confirmed different SARS-CoV-2-specific antibody profiles among the three different vaccination strategies and also provided insights regarding BCR usage and T/B cell activation and differentiation, which will guide a better selection of vaccination strategies in the future.

IMPORTANCE

Using the same laboratory test to avoid unnecessary interference due to cohort ethnicity, and experimental and statistical errors, we have compared the T/B cell immune response in the same cohort sequential vaccinated by different types of COVID-19 vaccine. We found that different sequential vaccinations can induce different dominant BCR usage with no significant neutralizing titers and RBD+ B-cell phenotype. Recombinant protein vaccine can induce higher numbers of regulatory T cells, circulating TFH (CTFH)1, CTFH17, and CTFH-CM, and lower SMNE+ T-cell proliferative capacity than the other two groups, whereas I-I-A showed higher proportion and number of virus-specific CD4+ T cells than I-I-R. Overall, our study provides a deep insight about the source of differences in immune protection of different types of COVID-19 vaccines, which further improves our understanding of the mechanisms underlying the immune response to SARS-CoV-2.

Antibody sequence determinants of viral antigen specificity

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

bcRflow: a Nextflow pipeline for characterizing B cell receptor repertoires from non-targeted transcriptomic data

B cells play a critical role in the adaptive recognition of foreign antigens through diverse receptor generation. While targeted immune sequencing methods are commonly used to profile B cell receptors (BCRs), they have limitations in cost and tissue availability. Analyzing B cell receptor profiling from non-targeted transcriptomics data is a promising alternative, but a systematic pipeline integrating tools for accurate immune repertoire extraction is lacking. Here, we present bcRflow, a Nextflow pipeline designed to characterize BCR repertoires from non-targeted transcriptomics data, with functional modules for alignment, processing, and visualization. bcRflow is a comprehensive, reproducible, and scalable pipeline that can run on high-performance computing clusters, cloud-based computing resources like Amazon Web Services (AWS), the Open OnDemand framework, or even local desktops. bcRflow utilizes institutional configurations provided by nf-core to ensure maximum portability and accessibility. To demonstrate the functionality of the bcRflow pipeline, we analyzed a public dataset of bulk transcriptomic samples from COVID-19 patients and healthy controls. We have shown that bcRflow streamlines the analysis of BCR repertoires from non-targeted transcriptomics data, providing valuable insights into the B cell immune response for biological and clinical research. bcRflow is available at https://github.com/Bioinformatics-Core-at-Childrens/bcRflow.

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.

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.

Current advancements in B-cell receptor sequencing fast-track the development of synthetic antibodies

Synthetic antibodies (Abs) are a class of engineered proteins designed to mimic the functions of natural Abs. These are produced entirely in vitro, eliminating the need for an immune response. As such, synthetic Abs have transformed the traditional methods of raising Abs. Likewise, deep sequencing technologies have revolutionized genomics and molecular biology. These enable the rapid and cost-effective sequencing of DNA and RNA molecules. They have allowed for accurate and inexpensive analysis of entire genomes and transcriptomes. Notably, via deep sequencing it is now possible to sequence a person's entire B-cell receptor immune repertoire, termed BCR sequencing. This procedure allows for big data explorations of natural Abs associated with an immune response. Importantly, the identified sequences have the ability to improve the design and engineering of synthetic Abs by offering an initial sequence framework for downstream optimizations. Additionally, machine learning algorithms can be introduced to leverage the vast amount of BCR sequencing datasets to rapidly identify patterns hidden in big data to effectively make in silico predictions of antigen selective synthetic Abs. Thus, the convergence of BCR sequencing, machine learning, and synthetic Ab development has effectively promoted a new era in Ab therapeutics. The combination of these technologies is driving rapid advances in precision medicine, diagnostics, and personalized treatments.

Steering and controlling evolution - from bioengineering to fighting pathogens

Control interventions steer the evolution of molecules, viruses, microorganisms or other cells towards a desired outcome. Applications range from engineering biomolecules and synthetic organisms to drug, therapy and vaccine design against pathogens and cancer. In all these instances, a control system alters the eco-evolutionary trajectory of a target system, inducing new functions or suppressing escape evolution. Here, we synthesize the objectives, mechanisms and dynamics of eco-evolutionary control in different biological systems. We discuss how the control system learns and processes information about the target system by sensing or measuring, through adaptive evolution or computational prediction of future trajectories. This information flow distinguishes pre-emptive control strategies by humans from feedback control in biotic systems. We establish a cost-benefit calculus to gauge and optimize control protocols, highlighting the fundamental link between predictability of evolution and efficacy of pre-emptive control.

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.

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.

Potent neutralizing broad-spectrum antibody against SARS-CoV-2 generated from dual-antigen-specific B cells from convalescents

Several antibody therapeutics have been developed against SARS-CoV-2; however, they have attenuated neutralizing ability against variants. In this study, we generated multiple broadly neutralizing antibodies from B cells of convalescents, by using two types of receptor-binding domains, Wuhan strain and the Gamma variant as bait. From 172 antibodies generated, six antibodies neutralized all strains prior to the Omicron variant, and the five antibodies were able to neutralize some of the Omicron sub-strains. Structural analysis showed that these antibodies have a variety of characteristic binding modes, such as ACE2 mimicry. We subjected a representative antibody to the hamster infection model after introduction of the N297A modification, and observed a dose-dependent reduction of the lung viral titer, even at a dose of 2 mg/kg. These results demonstrated that our antibodies have certain antiviral activity as therapeutics, and highlighted the importance of initial cell-screening strategy for the efficient development of therapeutic antibodies.

Approaches for single-cell RNA sequencing across tissues and cell types

Single-cell sequencing of RNA (scRNA-seq) has advanced our understanding of cellular heterogeneity and signaling in developmental biology and disease. A large number of complementary assays have been developed to profile transcriptomes of individual cells, also in combination with other readouts, such as chromatin accessibility or antibody-based analysis of protein surface markers. As scRNA-seq technologies are advancing fast, it is challenging to establish robust workflows and up-to-date protocols that are best suited to address the large range of research questions. Here, we review scRNA-seq techniques from mRNA end-counting to total RNA in relation to their specific features and outline the necessary sample preparation steps and quality control measures. Based on our experience in dealing with the continuously growing portfolio from the perspective of a central single-cell facility, we aim to provide guidance on how workflows can be best automatized and share our experience in coping with the continuous expansion of scRNA-seq techniques.

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.

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.

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.

IgTreeZ, A Toolkit for Immunoglobulin Gene Lineage Tree-Based Analysis, Reveals CDR3s Are Crucial for Selection Analysis

Somatic hypermutation (SHM) is an important diversification mechanism that plays a part in the creation of immune memory. Immunoglobulin (Ig) variable region gene lineage trees were used over the last four decades to model SHM and the selection mechanisms operating on B cell clones. We hereby present IgTreeZ (Immunoglobulin Tree analyZer), a python-based tool that analyses many aspects of Ig gene lineage trees and their repertoires. Using simulations, we show that IgTreeZ can be reliably used for mutation and selection analyses. We used IgTreeZ on empirical data, found evidence for different mutation patterns in different B cell subpopulations, and gained insights into antigen-driven selection in corona virus disease 19 (COVID-19) patients. Most importantly, we show that including the CDR3 regions in selection analyses - which is only possible if these analyses are lineage tree-based - is crucial for obtaining correct results. Overall, we present a comprehensive lineage tree analysis tool that can reveal new biological insights into B cell repertoire dynamics.

Delineating COVID-19 immunological features using single-cell RNA sequencing

Understanding the molecular mechanisms of coronavirus disease 2019 (COVID-19) pathogenesis and immune response is vital for developing therapies. Single-cell RNA sequencing has been applied to delineate the cellular heterogeneity of the host response toward COVID-19 in multiple tissues and organs. Here, we review the applications and findings from over 80 original COVID-19 single-cell RNA sequencing studies as well as many secondary analysis studies. We describe that single-cell RNA sequencing reveals multiple features of COVID-19 patients with different severity, including cell populations with proportional alteration, COVID-19-induced genes and pathways, severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) infection in single cells, and adaptation of immune repertoire. We also collect published single-cell RNA sequencing datasets from original studies. Finally, we discuss the limitations in current studies and perspectives for future advance.

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.

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.

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.

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).

High-resolution analysis of long-term serum antibodies in humans following convalescence of SARS-CoV-2 infection

Long-term antibody responses to SARS-CoV-2 have focused on responses to full-length spike protein, specific domains within spike, or nucleoprotein. In this study, we used high-density peptide microarrays representing the complete proteome of SARS-CoV-2 to identify binding sites (epitopes) targeted by antibodies present in the blood of COVID-19 resolved cases at 5 months post-diagnosis. Compared to previous studies that evaluated epitope-specific responses early post-diagnosis (< 60 days), we found that epitope-specific responses to nucleoprotein and spike protein have contracted, and that responses to membrane protein have expanded. Although antibody titers to full-length spike and nucleoprotein remain steady over months, taken together our data suggest that the population of epitope-specific antibodies that contribute to this reactivity is dynamic and evolves over time. Further, the spike epitopes bound by polyclonal antibodies in COVID-19 convalescent serum samples aligned with known target sites that can neutralize viral activity suggesting that the maintenance of these antibodies might provide rapid serological immunity. Finally, the most dominant epitopes for membrane protein and spike showed high diagnostic accuracy providing novel biomarkers to refine blood-based antibody tests. This study provides new insights into the specific regions of SARS-CoV-2 targeted by serum antibodies long after infection.

MetaLogo: a heterogeneity-aware sequence logo generator and aligner

Sequence logos are used to visually display conservations and variations in short sequences. They can indicate the fixed patterns or conserved motifs in a batch of DNA or protein sequences. However, most of the popular sequence logo generators are based on the assumption that all the input sequences are from the same homologous group, which will lead to an overlook of the heterogeneity among the sequences during the sequence logo making process. Heterogeneous groups of sequences may represent clades of different evolutionary origins, or genes families with different functions. Therefore, it is essential to divide the sequences into different phylogenetic or functional groups to reveal their specific sequence motifs and conservation patterns. To solve these problems, we developed MetaLogo, which can automatically cluster the input sequences after multiple sequence alignment and phylogenetic tree construction, and then output sequence logos for multiple groups and aligned them in one figure. User-defined grouping is also supported by MetaLogo to allow users to investigate functional motifs in a more delicate and dynamic perspective. MetaLogo can highlight both the homologous and nonhomologous sites among sequences. MetaLogo can also be used to annotate the evolutionary positions and gene functions of unknown sequences, together with their local sequence characteristics. We provide users a public MetaLogo web server (http://metalogo.omicsnet.org), a standalone Python package (https://github.com/labomics/MetaLogo), and also a built-in web server available for local deployment. Using MetaLogo, users can draw informative, customized and publishable sequence logos without any programming experience to present and investigate new knowledge on specific sequence sets.

SARS-CoV-2 Host Immunogenetic Biomarkers

SARS-CoV-2 causes generally mild symptoms, with approximately 10-20% of cases progressing to severe disease. The pathophysiologic mechanisms by which SARS-CoV-2 causes severe disease are largely unknown. Data have indicated the involvement of different immunogenetic markers such as HLA, T, and B cells, to be associated with disease outcome. This has led to interest in these genes as potential biomarkers of SARS-CoV-2 susceptibility and for predicting prognosis and response to vaccines and other therapeutic strategies. In this chapter, we discussed outline protocols for characterizing these potential biomarkers and methods for identifying SARS-CoV-2 biomarkers using the Luminex^® 100/200 technology and next-generation sequencing.

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

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

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

The immunology of asymptomatic SARS-CoV-2 infection: what are the key questions?

An important challenge during the COVID-19 pandemic has been to understand asymptomatic disease and the extent to which this may be a source of transmission. As asymptomatic disease is by definition hard to screen for, there is a lack of clarity about this aspect of the COVID-19 spectrum. Studies have considered whether the prevalence of asymptomatic disease is determined by differences in age, demographics, viral load, duration of shedding, and magnitude or durability of immunity. It is clear that adaptive immunity is strongly activated during asymptomatic infection, but some features of the T cell and antibody response may differ from those in symptomatic disease. Areas that need greater clarity include the extent to which asymptomatic disease leads to persistent symptoms (long COVID), and the quality, quantity and durability of immune priming required to confer subsequent protection.

Longitudinal Dynamics of Human B-Cell Response at the Single-Cell Level in Response to Tdap Vaccination

To mount an adequate immune response against pathogens, stepwise mutation and selection processes are crucial functions of the adaptive immune system. To better characterize a successful vaccination response, we performed longitudinal (days 0, 5, 7, 10, and 14 after Boostrix vaccination) analysis of the single-cell transcriptome as well as the B-cell receptor (BCR) repertoire (scBCR-rep) in plasma cells of an immunized donor and compared it with baseline B-cell characteristics as well as flow cytometry findings. Based on the flow cytometry knowledge and literature findings, we discriminated individual B-cell subsets in the transcriptomics data and traced over-time maturation of plasmablasts/plasma cells (PB/PCs) and identified the pathways associated with the plasma cell maturation. We observed that the repertoire in PB/PCs differed from the baseline B-cell repertoire e.g., regarding expansion of unique clones in post-vaccination visits, high usage of IGHG1 in expanded clones, increased class-switching events post-vaccination represented by clonotypes spanning multiple IGHC classes and positive selection of CDR3 sequences over time. Importantly, the Variable gene family-based clustering of BCRs represented a similar measure as the gene-based clustering, but certainly improved the clustering of BCRs, as BCRs from duplicated Variable gene families could be clustered together. Finally, we developed a query tool to dissect the immune response to the components of the Boostrix vaccine. Using this tool, we could identify the BCRs related to anti-tetanus and anti-pertussis toxoid BCRs. Collectively, we developed a bioinformatic workflow which allows description of the key features of an ongoing (longitudinal) immune response, such as activation of PB/PCs, Ig class switching, somatic hypermutation, and clonal expansion, all of which are hallmarks of antigen exposure, followed by mutation & selection processes.

Single-cell immune profiling reveals distinct immune response in asymptomatic COVID-19 patients

While some individuals infected by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) present mild-to-severe disease, many SARS-CoV-2-infected individuals are asymptomatic. We sought to identify the distinction of immune response between asymptomatic and moderate patients. We performed single-cell transcriptome and T-cell/B-cell receptor (TCR/BCR) sequencing in 37 longitudinal collected peripheral blood mononuclear cell samples from asymptomatic, moderate, and severe patients with healthy controls. Asymptomatic patients displayed increased CD56^briCD16⁻ natural killer (NK) cells and upregulation of interferon-gamma in effector CD4⁺ and CD8⁺ T cells and NK cells. They showed more robust TCR clonal expansion, especially in effector CD4⁺ T cells, but lack strong BCR clonal expansion compared to moderate patients. Moreover, asymptomatic patients have lower interferon-stimulated genes (ISGs) expression in general but large interpatient variability, whereas moderate patients showed various magnitude and temporal dynamics of the ISGs expression across multiple cell populations but lower than a patient with severe disease. Our data provide evidence of different immune signatures to SARS-CoV-2 in asymptomatic infections.

Advances in understanding the formation and fate of B-cell memory in response to immunization or infection

Immunological memory has the potential to provide lifelong protection against recurrent infections. As such, it has been crucial to the success of vaccines. Yet, the recent pandemic has illuminated key gaps in our knowledge related to the factors influencing effective memory formation and the inability to predict the longevity of immune protection. In recent decades, researchers have acquired a number of novel and powerful tools with which to study the factors underpinning humoral memory. These tools have been used to study the B-cell fate decisions that occur within the germinal centre (GC), a site where responding B cells undergo affinity maturation and are one of the major routes for memory B cell and high-affinity long-lived plasma cell formation. The advent of single-cell sequencing technology has provided an enhanced resolution for studying fate decisions within the GC and cutting-edge techniques have enabled researchers to model this reaction with more accuracy both in vitro and in silico. Moreover, modern approaches to studying memory B cells have allowed us to gain a better appreciation for the heterogeneity and adaptability of this vital class of B cells. Together, these studies have facilitated important breakthroughs in our understanding of how these systems operate to ensure a successful immune response. In this review, we describe recent advances in the field of GC and memory B-cell biology in order to provide insight into how humoral memory is formed, as well as the potential for generating lasting immunity to novel pathogens such as severe acute respiratory syndrome coronavirus 2.

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

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

Massive surge of mRNA expression of clonal B-cell receptor in patients with COVID-19

Background

Antibody production is one of the primary mechanisms for recovery from coronavirus disease 2019 (COVID-19). It is speculated that massive clonal expansion of B cells, which can produce clinically meaningful neutralizing antibodies, occurs in patients who recover on the timing of acquiring adaptive immunity.

Methods

To evaluate fluctuations in clonal B cells and the size of the clones, we chronologically assessed the B-cell receptor (BCR) repertoire in three patients with COVID-19 who recovered around 10 days after symptom onset.

Results

We focused on the three dominant clonotypes (top 3) in each individual. The percentage frequencies of the top 3 clonotypes increased rapidly and accounted for 27.8 % on day 9 in patient 1, 10.4 % on day 12 in patient 2, and 10.8 % on day 11 in patient 3, respectively. The frequencies of these top 3 clonotypes rapidly decreased as the patients’ clinical symptoms improved. Furthermore, BCR network analysis revealed that accumulation of clusters composed of similar complementarity-determining region 3 (CDR3) sequences were rapidly formed, grew, and reached their maximum size around 10 days after symptom onset.

Conclusions

BCR repertoire analysis revealed that a massive surge of some unique BCRs occurs during the acquisition of adaptive immunity and recovery. The peaks were more prominent than expected. These results provide insight into the important role of BCRs in the recovery from COVID-19 and raise the possibility of developing neutralizing antibodies as COVID-19 immunotherapy.

[Dengue amidst COVID-19: challenges & control measures for the double burden]

Dengue, an arbovirus, is a public health treat in the tropics and sub-tropical climates worldwide. The disease incidence has grown dramatically worldwide, with an estimated 390 million dengue virus infection per year. Dengue has distinct epidemiological patterns which are associated with the four virus serotypes. All four serotypes can co-circulate within a region, in which a number of regions are hyperendemic for all 4 serotypes. Currently, there are no specific treatment or vaccine for the disease. Dengue prevention depends on vector control measures and early interventions. The COVID-19 pandemic has placed immense pressure on health care and management systems worldwide. During the COVID-19 pandemic, the situation was aggravated by the simultaneous dengue outbreaks, that has led to a double burden which has further impacted the healthcare sector, particularly in resource limited settings. This review article will focus on dengue epidemics during the COVID-19 pandemic and discuss on recent findings on immunological cascades between dengue and COVID-19 and, the impact on vaccine development.