Maturation of germinal center B cells after influenza virus vaccination in humans.

Germinal centers (GC) are microanatomical lymphoid structures where affinity-matured memory B cells and long-lived bone marrow plasma cells are primarily generated. It is unclear how the maturation of B cells within the GC impacts the breadth and durability of B cell responses to influenza vaccination in humans. We used fine needle aspiration of draining lymph nodes to longitudinally track antigen-specific GC B cell responses to seasonal influenza vaccination. Antigen-specific GC B cells persisted for at least 13 wk after vaccination in two out of seven individuals. Monoclonal antibodies (mAbs) derived from persisting GC B cell clones exhibit enhanced binding affinity and breadth to influenza hemagglutinin (HA) antigens compared with related GC clonotypes isolated earlier in the response. Structural studies of early and late GC-derived mAbs from one clonal lineage in complex with H1 and H5 HAs revealed an altered binding footprint. Our study shows that inducing sustained GC reactions after influenza vaccination in humans supports the maturation of responding B cells.

Inferring B Cell Phylogenies from Paired H and L Chain BCR Sequences with Dowser.

Abs are vital to human immune responses and are composed of genetically variable H and L chains. These structures are initially expressed as BCRs. BCR diversity is shaped through somatic hypermutation and selection during immune responses. This evolutionary process produces B cell clones, cells that descend from a common ancestor but differ by mutations. Phylogenetic trees inferred from BCR sequences can reconstruct the history of mutations within a clone. Until recently, BCR sequencing technologies separated H and L chains, but advancements in single-cell sequencing now pair H and L chains from individual cells. However, it is unclear how these separate genes should be combined to infer B cell phylogenies. In this study, we investigated strategies for using paired H and L chain sequences to build phylogenetic trees. We found that incorporating L chains significantly improved tree accuracy and reproducibility across all methods tested. This improvement was greater than the difference between tree-building methods and persisted even when mixing bulk and single-cell sequencing data. However, we also found that many phylogenetic methods estimated significantly biased branch lengths when some L chains were missing, such as when mixing single-cell and bulk BCR data. This bias was eliminated using maximum likelihood methods with separate branch lengths for H and L chain gene partitions. Thus, we recommend using maximum likelihood methods with separate H and L chain partitions, especially when mixing data types. We implemented these methods in the R package Dowser: https://dowser.readthedocs.io.

CD23+IgG1+ memory B cells are poised to switch to pathogenic IgE production in food allergy.

Food allergy is caused by allergen-specific immunoglobulin E (IgE) antibodies, but little is known about the B cell memory of persistent IgE responses. Here, we describe, in human pediatric peanut allergy, a population of CD23+IgG1+ memory B cells arising in type 2 immune responses that contain high-affinity peanut-specific clones and generate IgE-producing cells upon activation. The frequency of CD23+IgG1+ memory B cells correlated with circulating concentrations of IgE in children with peanut allergy. A corresponding population of "type 2-marked" IgG1+ memory B cells was identified in single-cell RNA sequencing experiments. These cells differentially expressed interleukin-4 (IL-4)- and IL-13-regulated genes, such as FCER2/CD23+, IL4R, and germline IGHE, and carried highly mutated B cell receptors (BCRs). In children with high concentrations of serum peanut-specific IgE, high-affinity B cells that bind the main peanut allergen Ara h 2 mapped to the population of "type 2-marked" IgG1+ memory B cells and included clones with convergent BCRs across different individuals. Our findings indicate that CD23+IgG1+ memory B cells transcribing germline IGHE are a unique memory population containing precursors of high-affinity pathogenic IgE-producing cells that are likely to be involved in the long-term persistence of peanut allergy.

B cell phylogenetics in the single cell era.

The widespread availability of single-cell RNA sequencing (scRNA-seq) has led to the development of new methods for understanding immune responses. Single-cell transcriptome data can now be paired with B cell receptor (BCR) sequences. However, RNA from BCRs cannot be analyzed like most other genes because BCRs are genetically diverse within individuals. In humans, BCRs are shaped through recombination followed by mutation and selection for antigen binding. As these processes co-occur with cell division, B cells can be studied using phylogenetic trees representing the mutations within a clone. B cell trees can link experimental timepoints, tissues, or cellular subtypes. Here, we review the current state and potential of how B cell phylogenetics can be combined with single-cell data to understand immune responses.

Inferring B cell phylogenies from paired heavy and light chain BCR sequences with Dowser.

Antibodies are vital to human immune responses and are composed of genetically variable heavy and light chains. These structures are initially expressed as B cell receptors (BCRs). BCR diversity is shaped through somatic hypermutation and selection during immune responses. This evolutionary process produces B cell clones, cells that descend from a common ancestor but differ by mutations. Phylogenetic trees inferred from BCR sequences can reconstruct the history of mutations within a clone. Until recently, BCR sequencing technologies separated heavy and light chains, but advancements in single cell sequencing now pair heavy and light chains from individual cells. However, it is unclear how these separate genes should be combined to infer B cell phylogenies. In this study, we investigated strategies for using paired heavy and light chain sequences to build phylogenetic trees. We found incorporating light chains significantly improved tree accuracy and reproducibility across all methods tested. This improvement was greater than the difference between tree building methods and persisted even when mixing bulk and single cell sequencing data. However, we also found that many phylogenetic methods estimated significantly biased branch lengths when some light chains were missing, such as when mixing single cell and bulk BCR data. This bias was eliminated using maximum likelihood methods with separate branch lengths for heavy and light chain gene partitions. Thus, we recommend using maximum likelihood methods with separate heavy and light chain partitions, especially when mixing data types. We implemented these methods in the R package Dowser: https://dowser.readthedocs.io.

Tracking B cell responses to the SARS-CoV-2 mRNA-1273 vaccine.

Protective immunity following vaccination is sustained by long-lived antibody-secreting cells and resting memory B cells (MBCs). Responses to two-dose SARS-CoV-2 mRNA-1273 vaccination are evaluated longitudinally by multimodal single-cell analysis in three infection-naïve individuals. Integrated surface protein, transcriptomics, and B cell receptor (BCR) repertoire analysis of sorted plasmablasts and spike+ (S-2P+) and S-2P- B cells reveal clonal expansion and accumulating mutations among S-2P+ cells. These cells are enriched in a cluster of immunoglobulin G-expressing MBCs and evolve along a bifurcated trajectory rooted in CXCR3+ MBCs. One branch leads to CD11c+ atypical MBCs while the other develops from CD71+ activated precursors to resting MBCs, the dominant population at month 6. Among 12 evolving S-2P+ clones, several are populated with plasmablasts at early timepoints as well as CD71+ activated and resting MBCs at later timepoints, and display intra- and/or inter-cohort BCR convergence. These relationships suggest a coordinated and predictable evolution of SARS-CoV-2 vaccine-generated MBCs.

The memory of pathogenic IgE is contained within CD23 + IgG1 + memory B cells poised to switch to IgE in food allergy.

Food allergy is caused by allergen-specific IgE antibodies but little is known about the B cell memory of persistent IgE responses. Here we describe in human pediatric peanut allergy CD23 + IgG1 + memory B cells arising in type 2 responses that contain peanut specific clones and generate IgE cells on activation. These 'type2-marked' IgG1 + memory B cells differentially express IL-4/IL-13 regulated genes FCER2 / CD23, IL4R , and germline IGHE and carry highly mutated B cell receptors (BCRs). Further, high affinity memory B cells specific for the main peanut allergen Ara h 2 mapped to the population of 'type2-marked' IgG1 + memory B cells and included convergent BCRs across different individuals. Our findings indicate that CD23 + IgG1 + memory B cells transcribing germline IGHE are a unique memory population containing precursors of pathogenic IgE. One-Sentence Summary: We describe a unique population of IgG + memory B cells poised to switch to IgE that contains high affinity allergen-specific clones in peanut allergy.

Age-associated B cells are heterogeneous and dynamic drivers of autoimmunity in mice.

Age-associated B cells (ABCs) are formed under inflammatory conditions and are considered a type of memory B cell (MBC) expressing the transcription factor T-bet. In SLE, ABC frequency is correlated with disease, and they are thought to be the source of autoantibody-secreting cells. However, in inflammatory conditions, whether autoreactive B cells can become resting MBCs is uncertain. Further, the phenotypic identity of ABCs and their relationship to other B cell subsets, such as plasmablasts, is unclear. Whether ABCs directly promote disease is untested. Here we report, in the MRL/lpr SLE model, unexpected heterogeneity among ABC-like cells for expression of the integrins CD11b and CD11c, T-bet, and memory or plasmablast markers. Transfer and labeling studies demonstrated that ABCs are dynamic, rapidly turning over. scRNA-seq identified B cell clones present in multiple subsets, revealing that ABCs can be plasmablast precursors or undergo cycles of reactivation. Deletion of CD11c-expressing B cells revealed a direct role for ABC-like B cells in lupus pathogenesis.

PD-1highCXCR5-CD4+ peripheral helper T cells promote CXCR3+ plasmablasts in human acute viral infection.

T cell-B cell interaction is the key immune response to protect the host from severe viral infection. However, how T cells support B cells to exert protective humoral immunity in humans is not well understood. Here, we use COVID-19 as a model of acute viral infections and analyze CD4+ T cell subsets associated with plasmablast expansion and clinical outcome. Peripheral helper T cells (Tph cells; denoted as PD-1highCXCR5-CD4+ T cells) are significantly increased, as are plasmablasts. Tph cells exhibit "B cell help" signatures and induce plasmablast differentiation in vitro. Interestingly, expanded plasmablasts show increased CXCR3 expression, which is positively correlated with higher frequency of activated Tph cells and better clinical outcome. Mechanistically, Tph cells help B cell differentiation and produce more interferon γ (IFNγ), which induces CXCR3 expression on plasmablasts. These results elucidate a role for Tph cells in regulating protective B cell response during acute viral infection.

BTG1 mutation yields supercompetitive B cells primed for malignant transformation.

Multicellular life requires altruistic cooperation between cells. The adaptive immune system is a notable exception, wherein germinal center B cells compete vigorously for limiting positive selection signals. Studying primary human lymphomas and developing new mouse models, we found that mutations affecting BTG1 disrupt a critical immune gatekeeper mechanism that strictly limits B cell fitness during antibody affinity maturation. This mechanism converted germinal center B cells into supercompetitors that rapidly outstrip their normal counterparts. This effect was conferred by a small shift in MYC protein induction kinetics but resulted in aggressive invasive lymphomas, which in humans are linked to dire clinical outcomes. Our findings reveal a delicate evolutionary trade-off between natural selection of B cells to provide immunity and potentially dangerous features that recall the more competitive nature of unicellular organisms.

Adaptive immune responses to SARS-CoV-2 persist in the pharyngeal lymphoid tissue of children.

Most studies of adaptive immunity to SARS-CoV-2 infection focus on peripheral blood, which may not fully reflect immune responses at the site of infection. Using samples from 110 children undergoing tonsillectomy and adenoidectomy during the COVID-19 pandemic, we identified 24 samples with evidence of previous SARS-CoV-2 infection, including neutralizing antibodies in serum and SARS-CoV-2-specific germinal center and memory B cells in the tonsils and adenoids. Single-cell B cell receptor (BCR) sequencing indicated virus-specific BCRs were class-switched and somatically hypermutated, with overlapping clones in the two tissues. Expanded T cell clonotypes were found in tonsils, adenoids and blood post-COVID-19, some with CDR3 sequences identical to previously reported SARS-CoV-2-reactive T cell receptors (TCRs). Pharyngeal tissues from COVID-19-convalescent children showed persistent expansion of germinal center and antiviral lymphocyte populations associated with interferon (IFN)-γ-type responses, particularly in the adenoids, and viral RNA in both tissues. Our results provide evidence for persistent tissue-specific immunity to SARS-CoV-2 in the upper respiratory tract of children after infection.

Reemergence of pathogenic, autoantibody-producing B cell clones in myasthenia gravis following B cell depletion therapy.

Myasthenia gravis (MG) is an autoantibody-mediated autoimmune disorder of the neuromuscular junction. A small subset of patients (<10%) with MG, have autoantibodies targeting muscle-specific tyrosine kinase (MuSK). MuSK MG patients respond well to CD20-mediated B cell depletion therapy (BCDT); most achieve complete stable remission. However, relapse often occurs. To further understand the immunomechanisms underlying relapse, we studied autoantibody-producing B cells over the course of BCDT. We developed a fluorescently labeled antigen to enrich for MuSK-specific B cells, which was validated with a novel Nalm6 cell line engineered to express a human MuSK-specific B cell receptor. B cells (≅ 2.6 million) from 12 different samples collected from nine MuSK MG patients were screened for MuSK specificity. We successfully isolated two MuSK-specific IgG4 subclass-expressing plasmablasts from two of these patients, who were experiencing a relapse after a BCDT-induced remission. Human recombinant MuSK mAbs were then generated to validate binding specificity and characterize their molecular properties. Both mAbs were strong MuSK binders, they recognized the Ig1-like domain of MuSK, and showed pathogenic capacity when tested in an acetylcholine receptor (AChR) clustering assay. The presence of persistent clonal relatives of these MuSK-specific B cell clones was investigated through B cell receptor repertoire tracing of 63,977 unique clones derived from longitudinal samples collected from these two patients. Clonal variants were detected at multiple timepoints spanning more than five years and reemerged after BCDT-mediated remission, predating disease relapse by several months. These findings demonstrate that a reservoir of rare pathogenic MuSK autoantibody-expressing B cell clones survive BCDT and reemerge into circulation prior to manifestation of clinical relapse. Overall, this study provides both a mechanistic understanding of MuSK MG relapse and a valuable candidate biomarker for relapse prediction.

Adaptive Immune Receptor Repertoire (AIRR) Community Guide to Repertoire Analysis.

Adaptive immune receptor repertoires (AIRRs) are rich with information that can be mined for insights into the workings of the immune system. Gene usage, CDR3 properties, clonal lineage structure, and sequence diversity are all capable of revealing the dynamic immune response to perturbation by disease, vaccination, or other interventions. Here we focus on a conceptual introduction to the many aspects of repertoire analysis and orient the reader toward the uses and advantages of each. Along the way, we note some of the many software tools that have been developed for these investigations and link the ideas discussed to chapters on methods provided elsewhere in this volume.

Phylogenetic analysis of migration, differentiation, and class switching in B cells.

B cells undergo rapid mutation and selection for antibody binding affinity when producing antibodies capable of neutralizing pathogens. This evolutionary process can be intermixed with migration between tissues, differentiation between cellular subsets, and switching between functional isotypes. B cell receptor (BCR) sequence data has the potential to elucidate important information about these processes. However, there is currently no robust, generalizable framework for making such inferences from BCR sequence data. To address this, we develop three parsimony-based summary statistics to characterize migration, differentiation, and isotype switching along B cell phylogenetic trees. We use simulations to demonstrate the effectiveness of this approach. We then use this framework to infer patterns of cellular differentiation and isotype switching from high throughput BCR sequence datasets obtained from patients in a study of HIV infection and a study of food allergy. These methods are implemented in the R package dowser, available at https://dowser.readthedocs.io.

Robust, persistent adaptive immune responses to SARS-CoV-2 in the oropharyngeal lymphoid tissue of children.

SARS-CoV-2 infection triggers adaptive immune responses from both T and B cells. However, most studies focus on peripheral blood, which may not fully reflect immune responses in lymphoid tissues at the site of infection. To evaluate both local and systemic adaptive immune responses to SARS-CoV-2, we collected peripheral blood, tonsils, and adenoids from 110 children undergoing tonsillectomy/adenoidectomy during the COVID-19 pandemic and found 24 with evidence of prior SARS-CoV-2 infection, including detectable neutralizing antibodies against multiple viral variants. We identified SARS-CoV-2-specific germinal center (GC) and memory B cells; single cell BCR sequencing showed that these virus-specific B cells were class-switched and somatically hypermutated, with overlapping clones in the adenoids and tonsils. Oropharyngeal tissues from COVID-19-convalescent children showed persistent expansion of GC and anti-viral lymphocyte populations associated with an IFN-γ-type response, with particularly prominent changes in the adenoids, as well as evidence of persistent viral RNA in both tonsil and adenoid tissues of many participants. Our results show robust, tissue-specific adaptive immune responses to SARS-CoV-2 in the upper respiratory tract of children weeks to months after acute infection, providing evidence of persistent localized immunity to this respiratory virus.

Single-cell multi-omics reveals dyssynchrony of the innate and adaptive immune system in progressive COVID-19.

Dysregulated immune responses against the SARS-CoV-2 virus are instrumental in severe COVID-19. However, the immune signatures associated with immunopathology are poorly understood. Here we use multi-omics single-cell analysis to probe the dynamic immune responses in hospitalized patients with stable or progressive course of COVID-19, explore V(D)J repertoires, and assess the cellular effects of tocilizumab. Coordinated profiling of gene expression and cell lineage protein markers shows that S100Ahi/HLA-DRlo classical monocytes and activated LAG-3hi T cells are hallmarks of progressive disease and highlights the abnormal MHC-II/LAG-3 interaction on myeloid and T cells, respectively. We also find skewed T cell receptor repertories in expanded effector CD8+ clones, unmutated IGHG+ B cell clones, and mutated B cell clones with stable somatic hypermutation frequency over time. In conclusion, our in-depth immune profiling reveals dyssynchrony of the innate and adaptive immune interaction in progressive COVID-19.

B Cell Mobilization, Dissemination, Fine Tuning of Local Antigen Specificity and Isotype Selection in Asthma.

In order to better understand how the immune system interacts with environmental triggers to produce organ-specific disease, we here address the hypothesis that B and plasma cells are free to migrate through the mucosal surfaces of the upper and lower respiratory tracts, and that their total antibody repertoire is modified in a common respiratory tract disease, in this case atopic asthma. Using Adaptive Immune Receptor Repertoire sequencing (AIRR-seq) we have catalogued the antibody repertoires of B cell clones retrieved near contemporaneously from multiple sites in the upper and lower respiratory tract mucosa of adult volunteers with atopic asthma and non-atopic controls and traced their migration. We show that the lower and upper respiratory tracts are immunologically connected, with trafficking of B cells directionally biased from the upper to the lower respiratory tract and points of selection when migrating from the nasal mucosa and into the bronchial mucosa. The repertoires are characterized by both IgD-only B cells and others undergoing class switch recombination, with restriction of the antibody repertoire distinct in asthmatics compared with controls. We conclude that B cells and plasma cells migrate freely throughout the respiratory tract and exhibit distinct antibody repertoires in health and disease.

Human B cell lineages associated with germinal centers following influenza vaccination are measurably evolving.

The poor efficacy of seasonal influenza virus vaccines is often attributed to pre-existing immunity interfering with the persistence and maturation of vaccine-induced B cell responses. We previously showed that a subset of vaccine-induced B cell lineages are recruited into germinal centers (GCs) following vaccination, suggesting that affinity maturation of these lineages against vaccine antigens can occur. However, it remains to be determined whether seasonal influenza vaccination stimulates additional evolution of vaccine-specific lineages, and previous work has found no significant increase in somatic hypermutation among influenza-binding lineages sampled from the blood following seasonal vaccination in humans. Here, we investigate this issue using a phylogenetic test of measurable immunoglobulin sequence evolution. We first validate this test through simulations and survey measurable evolution across multiple conditions. We find significant heterogeneity in measurable B cell evolution across conditions, with enrichment in primary response conditions such as HIV infection and early childhood development. We then show that measurable evolution following influenza vaccination is highly compartmentalized: while lineages in the blood are rarely measurably evolving following influenza vaccination, lineages containing GC B cells are frequently measurably evolving. Many of these lineages appear to derive from memory B cells. We conclude from these findings that seasonal influenza virus vaccination can stimulate additional evolution of responding B cell lineages, and imply that the poor efficacy of seasonal influenza vaccination is not due to a complete inhibition of vaccine-specific B cell evolution.

When the immune system encounters a disease-causing pathogen, it releases antibodies that can bind to specific regions of the bacterium or virus and help to clear the infection. These proteins are generated by B cells which, upon detecting the pathogen, can begin to mutate and alter the structure of the antibody they produce: the better the antibody is at binding to the pathogen, the more likely the B cell is to survive. This process of evolution produces B cells that make more effective antibodies. After the infection, some of these cells become 'memory B cells' which can be stimulated in to action when the pathogen invades again. Many vaccines also depend on this process to trigger the production of memory B cells that can fight off a specific disease-causing agent. However, it is unclear to what extent memory B cells that already exist are able to continue to evolve and modify their antibodies. This is particularly important for the flu vaccine, as the virus that causes influenza rapidly mutates. To provide high levels of protection, the memory B cells formed following the vaccine may therefore need to evolve to make different antibodies that recognize mutated forms of the virus. It is thought that the low effectiveness of the flu vaccine is partially because the response it triggers does not stimulate additional evolution of memory B cells. To test this theory, Hoehn et al. developed a computational method that can detect the evolution of B cells over time. The tool was applied to samples collected from the blood and lymph nodes (organ where immune cells reside) of people who recently received the flu vaccine. The results were then compared to B cells taken from people after different infections, vaccinations, and other conditions. Hoehn et al. found the degree to which B cells evolve varies significantly between conditions. For example, B cells produced during chronic HIV infections frequently evolved over time, while such evolution was rarely observed during the autoimmune disease myasthenia gravis. The analysis also showed that memory B cells produced by the flu vaccine were able to evolve if recruited to the lymph nodes, but this was rarely detected in B cells in the blood. These findings suggest the low efficacy of the flu vaccine is not due to a complete lack of B cell evolution, but likely due to other factors. For instance, it is possible the evolutionary process it stimulates is not as robust as in other conditions, or is less likely to produce long-lived B cells that release antibodies. More research is needed to explore these ideas and could lead to the development of more effective flu vaccines.

Elevated N-Linked Glycosylation of IgG V Regions in Myasthenia Gravis Disease Subtypes.

Elevated N-linked glycosylation of IgG V regions (IgG-VN-Glyc) is an emerging molecular phenotype associated with autoimmune disorders. To test the broader specificity of elevated IgG-VN-Glyc, we studied patients with distinct subtypes of myasthenia gravis (MG), a B cell-mediated autoimmune disease. Our experimental design focused on examining the B cell repertoire and total IgG. It specifically included adaptive immune receptor repertoire sequencing to quantify and characterize N-linked glycosylation sites in the circulating BCR repertoire, proteomics to examine glycosylation patterns of the total circulating IgG, and an exploration of human-derived recombinant autoantibodies, which were studied with mass spectrometry and Ag binding assays to respectively confirm occupation of glycosylation sites and determine whether they alter binding. We found that the frequency of IgG-VN-Glyc motifs was increased in the total BCR repertoire of patients with MG when compared with healthy donors. The elevated frequency was attributed to both biased V gene segment usage and somatic hypermutation. IgG-VN-Glyc could be observed in the total circulating IgG in a subset of patients with MG. Autoantigen binding, by four patient-derived MG autoantigen-specific mAbs with experimentally confirmed presence of IgG-VN-Glyc, was not altered by the glycosylation. Our findings extend prior work on patterns of Ig V region N-linked glycosylation in autoimmunity to MG subtypes.

Single-cell immunophenotyping of the skin lesion erythema migrans identifies IgM memory B cells.

The skin lesion erythema migrans (EM) is an initial sign of the Ixodes tick-transmitted Borreliella spirochetal infection known as Lyme disease. T cells and innate immune cells have previously been shown to predominate the EM lesion and promote the reaction. Despite the established importance of B cells and antibodies in preventing infection, the role of B cells in the skin immune response to Borreliella is unknown. Here, we used single-cell RNA-Seq in conjunction with B cell receptor (BCR) sequencing to immunophenotype EM lesions and their associated B cells and BCR repertoires. We found that B cells were more abundant in EM in comparison with autologous uninvolved skin; many were clonally expanded and had circulating relatives. EM-associated B cells upregulated the expression of MHC class II genes and exhibited preferential IgM isotype usage. A subset also exhibited low levels of somatic hypermutation despite a gene expression profile consistent with memory B cells. Our study demonstrates that single-cell gene expression with paired BCR sequencing can be used to interrogate the sparse B cell populations in human skin and reveals that B cells in the skin infection site in early Lyme disease expressed a phenotype consistent with local antigen presentation and antibody production.