Change-O: a toolkit for analyzing large-scale B cell immunoglobulin repertoire sequencing data

IL-27 elicits a cytotoxic CD8+ T cell program to enforce tumour control

Although cytotoxic CD8+ T lymphocytes (CTLs) are essential for anti-tumour immunity, they are frequently dysfunctional in tumours1. Cytokines that sustain CTL activity are attractive for cancer immunotherapy, but avoiding inflammatory toxicity remains a challenge for their clinical use2. Here we show that expression of a CTL signature is strongly associated with IL27 expression in human and mouse tumours. In mice, IL-27 acts directly on tumour-specific CTLs to promote their persistence and effector function in the tumour microenvironment. Moreover, treatment with inducible IL-27 overexpression or a half-life-extended IL-27 protein in vivo is well tolerated, induces regression of established tumours, drives an enhanced cytotoxic program in anti-tumour CTLs and synergizes with PD-L1 blockade. In patients with cancer who were treated with anti-PD-1/PD-L1 therapy, high expression of IL-27 correlates with a favourable clinical response, and IL-27 supports human CTL function during chronic antigen stimulation ex vivo. Our data demonstrate that endogenous IL-27 is essential for anti-tumour immunity and that IL-27 receptor agonism can safely improve anti-tumour T cell responses alone or in combination with PD-L1 blockade.

Robust mucosal SARS-CoV-2-specific T cells effectively combat COVID-19 and establish polyfunctional resident memory in patient lungs

Mucosal antigen-specific T cells are pivotal for pathogen clearance and immune modulation in respiratory infections. Dysregulated T cell responses exacerbate coronavirus disease 2019 severity, marked by cytokine storms and respiratory failure. Despite extensive description in peripheral blood, the characteristics of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)-specific T cells in the lungs remain elusive. Here we conducted integrated single-cell profiling of SARS-CoV-2-specific T cells in 122 bronchoalveolar lavage fluid (BALF) and 280 blood samples from 159 patients, including 27 paired BALF and blood samples from 24 patients. SARS-CoV-2-specific T cells were robustly elicited in BALF irrespective of prior vaccination, correlating with diminished viral loads, lessened systemic inflammation and improved respiratory function. SARS-CoV-2-specific T cells in BALF exhibited profound activation, along with proliferative and multi-cytokine-producing capabilities and a glycolysis-driven metabolic signature, which were distinct from those observed in peripheral blood mononuclear cells. After viral clearance, these specific T cells maintained a polyfunctional tissue-resident memory phenotype, highlighting their critical roles in infection control and long-term protection.

Vaccination with mRNA-encoded membrane-bound HIV Envelope trimer induces neutralizing antibodies in animal models

A protective vaccine against HIV will likely need to induce broadly neutralizing antibodies (bnAbs) that engage relatively conserved epitopes on the HIV envelope glycoprotein (Env) trimer. Nearly all vaccine strategies to induce bnAbs require the use of relatively complex immunization regimens involving a series of different immunogens, most of which are Env trimers. Producing protein-based clinical material to evaluate such relatively complex regimens in humans presents major challenges in cost and time. Furthermore, immunization with HIV trimers as soluble proteins induces strong non-neutralizing responses to the trimer base, which is normally occluded on the virion. These base responses could potentially detract from the induction of nAbs and the eventual induction of bnAbs. mRNA vaccine platforms offer potential advantages over protein delivery for HIV vaccine development, including increased production speed, reduced cost, and the ability to deliver membrane-bound trimers that might facilitate improved immuno-focusing to non-base epitopes. We report the design of mRNA-delivered soluble and membrane-bound forms of a stabilized native-like Env trimer (BG505 MD39.3), initial immunogenicity evaluation in rabbits that triggered clinical evaluation, and more comprehensive evaluation of B cell, T cell, and antibody responses in non-human primates. mRNA-encoded membrane-bound Env immunization elicited reduced off-target base-directed Env responses and stronger neutralizing antibody responses, compared with mRNA-encoded soluble Env. Overall, mRNA delivery of membrane-bound Env appears promising for enhancing B cell responses to subdominant epitopes and facilitating rapid translation to clinical testing, which should assist HIV vaccine development.

One Sentence Summary

HIV envelope trimer mRNA enables membrane-bound expression and represents a functional immunogen in pre-clinical mammalian models.

Single-cell analysis of cerebrospinal fluid reveals common features of neuroinflammation

Neuroinflammation is often characterized by immune cell infiltrates in the cerebrospinal fluid (CSF). Here, we apply single-cell RNA sequencing to explore the functional characteristics of these cells in patients with various inflammatory, infectious, and non-inflammatory neurological disorders. We show that CSF is distinct from the peripheral blood in terms of both cellular composition and gene expression. We report that the cellular and transcriptional landscape of CSF is altered in neuroinflammation but is strikingly similar across different neuroinflammatory disorders. We find clonal expansion of CSF lymphocytes in all disorders but most pronounced in inflammatory diseases, and we functionally characterize the transcriptional features of these cells. Finally, we explore the genetic control of gene expression in CSF lymphocytes. Our results highlight the common features of immune cells in the CSF compartment across diverse neurological diseases and may help to identify new targets for drug development or repurposing in multiple sclerosis (MS).

Systematic evaluation of intratumoral and peripheral BCR repertoires in three cancers

The current understanding of humoral immune response in cancer patients suggests that tumors may be infiltrated with diffuse B cells of extra-tumoral origin or may develop organized lymphoid structures, where somatic hypermutation and antigen-driven selection occur locally. These processes are believed to be significantly influenced by the tumor microenvironment through secretory factors and biased cell-cell interactions. To explore the manifestation of this influence, we used deep unbiased immunoglobulin profiling and systematically characterized the relationships between B cells in circulation, draining lymph nodes (draining LNs), and tumors in 14 patients with three human cancers. We demonstrated that draining LNs are differentially involved in the interaction with the tumor site, and that significant heterogeneity exists even between different parts of a single lymph node (LN). Next, we confirmed and elaborated upon previous observations regarding intratumoral immunoglobulin heterogeneity. We identified B cell receptor (BCR) clonotypes that were expanded in tumors relative to draining LNs and blood and observed that these tumor-expanded clonotypes were less hypermutated than non-expanded (ubiquitous) clonotypes. Furthermore, we observed a shift in the properties of complementarity-determining region 3 of the BCR heavy chain (CDR-H3) towards less mature and less specific BCR repertoire in tumor-infiltrating B-cells compared to circulating B-cells, which may indicate less stringent control for antibody-producing B cell development in tumor microenvironment (TME). In addition, we found repertoire-level evidence that B-cells may be selected according to their CDR-H3 physicochemical properties before they activate somatic hypermutation (SHM). Altogether, our work outlines a broad picture of the differences in the tumor BCR repertoire relative to non-tumor tissues and points to the unexpected features of the SHM process.

Profiling the broad antibody diversity of lymphatic filariasis immune antibody repertoire by deep sequencing

Lymphatic filariasis is caused by infections of thread-like filarial worms, namely Wuchereria bancrofti, Brugia Malayi and Brugia timori. However, in-depth analysis of the antibody repertoire against Lymphatic filariasis is lacking. Using high-throughput sequencing of antibody repertoires, immunome analysis of IgG (LG) and IgM (LM) repertoires were studied. Despite significant differences between LG and LM in V(D)J gene usage, IGHV4-34, IGHV6-1, IGHD3-10 and IGHJ4 were preferred in both repertoires. The CDR3 in the LG repertoire showed a longer length than LM. Higher SHM level were observed in LG sequences and presence of oligoclonal sequences indicates the extent of clonal expansion. The prevalence of rare clonotypes in LM repertoire depicts the high clonal diversity when compared to LG repertoire. Monoclonal antibodies against closely related parasitic infections were present within the LG repertoire suggesting that immune repertoires may not be as exclusive and biased against the target infection as initially thought. The characterization of the immune repertoire can provide critical insight into the antibody response patterns in disease state, antibody generation process during infections and future antibody designs.

Mutability and hypermutation antagonize immunoglobulin codon optimality

The efficacy of antibody responses is inherently linked to paratope diversity, as generated through V(D)J recombination and somatic hypermutation. Despite this, it is unclear how genetic diversification mechanisms evolved alongside codon optimality and affect antibody expression. Here, we analyze germline immunoglobulin (IG) genes, natural V(D)J repertoires, serum IgG, and monoclonal antibody (mAb) expression through the lens of codon optimality. Germline variable genes (IGVs) exhibit diverse optimality that is inversely related to mutability. Hypermutation deoptimizes heavy-chain (IGH) VDJ repertoires within human tonsils, bone marrow, lymph nodes (including SARS-CoV-2-specific clones), blood (HIV-1-specific clones), mice, and zebrafish. Analyses of mutation-affected codons show that targeting to complementarity-determining regions constrains deoptimization. Germline IGHV optimality correlates with serum variable fragment (VH) usage after influenza vaccination, while synonymous deoptimization attenuated mAb yield. These findings provide unanticipated insights into an antagonistic relationship between diversification mechanisms and codon optimality. Ultimately, the need for diversity takes precedence over that for the most optimal codon usage.

IL-33-activated ILC2s induce tertiary lymphoid structures in pancreatic cancer

Tertiary lymphoid structures (TLSs) are de novo ectopic lymphoid aggregates that regulate immunity in chronically inflamed tissues, including tumours. Although TLSs form due to inflammation-triggered activation of the lymphotoxin (LT)-LTβ receptor (LTβR) pathway1, the inflammatory signals and cells that induce TLSs remain incompletely identified. Here we show that interleukin-33 (IL-33), the alarmin released by inflamed tissues2, induces TLSs. In mice, Il33 deficiency severely attenuates inflammation- and LTβR-activation-induced TLSs in models of colitis and pancreatic ductal adenocarcinoma (PDAC). In PDAC, the alarmin domain of IL-33 activates group 2 innate lymphoid cells (ILC2s) expressing LT that engage putative LTβR+ myeloid organizer cells to initiate tertiary lymphoneogenesis. Notably, lymphoneogenic ILC2s migrate to PDACs from the gut, can be mobilized to PDACs in different tissues and are modulated by gut microbiota. Furthermore, we detect putative lymphoneogenic ILC2s and IL-33-expressing cells within TLSs in human PDAC that correlate with improved prognosis. To harness this lymphoneogenic pathway for immunotherapy, we engineer a recombinant human IL-33 protein that expands intratumoural lymphoneogenic ILC2s and TLSs and demonstrates enhanced anti-tumour activity in PDAC mice. In summary, we identify the molecules and cells of a druggable pathway that induces inflammation-triggered TLSs. More broadly, we reveal a lymphoneogenic function for alarmins and ILC2s.

Hypomorphic RAG2 Deficiency Promotes Selection of Self-Reactive B Cells

Reduced function or hypomorphic variants in recombination-activating genes (RAG) 1 or 2 result in a broad clinical phenotype including common variable immunodeficiency (CVID) and even adult-onset disease. Milder RAG variants are less characterized. Here we describe the longitudinal course of a milder combined RAG deficiency in 3 of 7 siblings sharing the same RAG2 mutations over a 50-year study. Whole-genome and repertoire sequencing, bacteriophage immunizations, and deep immunophenotyping were used to compare affected and unaffected family members. The clinical phenotype of three affected siblings with hypomorphic RAG deficiency ranged from combined immunodeficiency and early mortality to a late-onset CID with hyper-IgM phenotype. T cells were remarkably similar across affected siblings, yet CDR3 skewing and regulatory T cell defects were not observed. B cell analysis showed elevated unswitched CD27+ and CD21low cells as well as features of an autoreactive antibody repertoire and presence of secreted autoantibodies, yet no clinical autoimmunity was present. Most striking was an expanded polyclonal marginal zone-like B cell population (IgM+IgD+CD27+) utilizing the self-reactive unmutated VH4-34 receptor demonstrating that hypomorphic RAG deficiency can promote expansion of self-reactive B cells. This process, however, was not sufficient to trigger clinical autoimmunity. Utilizing multiple approaches, we functionally measured the specific RAG2 variant effects and assessed how selection and secondary triggers altered the BCR repertoire and immunophenotype overtime. Overall, we demonstrate a broad disease spectrum in siblings with identical hypomorphic RAG deficiency, highlighting that phenotypic divergence can result from expansion of IgM + memory B cells.

A Broad Survey and Functional Analysis of Immunoglobulin Loci Variation in Rhesus Macaques

Rhesus macaques (RMs) are vital models for studying human disease, and are invaluable to pre-clinical pipelines for vaccine discovery and testing. Particularly in this regard, they are often used to study infection and vaccine-associated broadly neutralizing antibody responses. This has resulted in an increasing demand for improved genetic resources for the immunoglobulin (IG) loci, which harbor antibody-encoding genes. However, the highly polymorphic and structurally variable nature of these loci have them historically challenging to sequence and characterize at the level of both the genome and expressed repertoire. To address these challenges, we have developed a novel integrated analysis workflow for conducting the combined processing of B cell receptor repertoire sequencing data with matched whole-genome and targeted long-read genomic sequencing data. Using this novel approach, we have assembled the largest collection of IG germline alleles reported to date, amassed from 106 Indian origin RMs. Using a conservative annotation approach, requiring sample-level internal validation from both genomic and expressed datasets, we created a comprehensive resource that captures extensive diversity of IG heavy and light chain variable (V), diversity (D), and joining (J) alleles, as well as leader, intronic, and recombination signal sequences (RSSs). This publicly available, continually updated database will advance vaccine research for infectious disease, and provide a robust foundation for immunogenomics and future translational research.

Single-Cell Transcriptomic Analysis of Normal and Malignant B Cells

In the past decade, single-cell (sc) transcriptomics has overcome the limitations of bulk analysis by measuring gene expression in individual cells, not just a population average. This can identify diverse cell types and states within a sample with high resolution, even without prior purification. Various technologies exist, each with its own capture, barcoding, and library preparation methods. This chapter focuses on the analysis of normal and malignant mature B cells using the 10× Genomics 5' sc-gene expression in parallel with B cell immune repertoire profiling. By integrating the gene expression data from similar cells, the complete transcriptome for each population can be reconstructed, while the identification of the expressed immunoglobulin genes allows investigating clonotype evolution and the detection of tumor clones that share the same clonally rearranged B cell receptor sequence. Researchers are guided through both the experimental protocols and data analysis with a comprehensive, step-by-step walkthrough of how to use some of the more popular single-cell software tools.

Vaccination with different group 2 influenza subtypes alters epitope targeting and breadth of hemagglutinin stem-specific human B cells

The conserved influenza hemagglutinin stem, which is a target of cross-neutralizing antibodies, is now used in vaccine strategies focused on protecting against influenza pandemics. Antibody responses to group 1 stem have been extensively characterized, but little is known about group 2. Here, we characterized the stem-specific repertoire of individuals vaccinated with one of three group 2 influenza subtypes (H3, H7, or H10). Epitope mapping revealed two epitope supersites on the group 2 stem. Antibodies targeting the central epitope were broadly cross-reactive, whereas antibodies targeting the lower epitope had narrower breadth but higher potency against H3 subtypes. The ratio of B cells targeting each of the supersites varied with the vaccine subtype, leading to differences in the cross-reactivity of the B cell response. Our findings suggest that vaccine strategies targeting both group 2 stem epitopes would be complementary, eliciting broader and more potent protection against both seasonal and pandemic influenza strains.

Evolving antibody response to SARS-CoV-2 antigenic shift from XBB to JN.1

The continuous evolution of SARS-CoV-2, particularly the emergence of the BA.2.86/JN.1 lineage replacing XBB, necessitates re-evaluation of vaccine compositions1-3. Here, we provide a comprehensive analysis of the humoral immune response to XBB and JN.1 human exposure. We demonstrate the antigenic distinctiveness of XBB and JN.1 lineages in SARS-CoV-2-naive individuals and show that infection with JN.1 elicits superior plasma neutralization against its subvariants. We highlight the strong immune evasion and receptor-binding capability of KP.3, supporting its foreseeable prevalence. Extensive analysis of the B cell receptor repertoire, in which we isolate approximately 2,000 receptor-binding-domain-specific antibodies, with targeting epitopes characterized by deep mutational scanning, underscores the superiority of JN.1-elicited memory B cells4,5. Class 1 IGHV3-53/3-66-derived neutralizing antibodies (NAbs) are important contributors to the wild-type reactivity of NAbs against JN.1. However, KP.2 and KP.3 evade a substantial subset of these antibodies, even those induced by JN.1, supporting a need for booster updates. JN.1-induced Omicron-specific antibodies also demonstrate high potency across Omicron. Escape hotspots for these NAbs have already been mutated, resulting in a higher immune barrier to escape and indicating probable recovery of escaped NAbs. In addition, the prevalence of IGHV3-53/3-66-derived antibodies and their ability to compete with all Omicron-specific NAbs suggests that they have an inhibitory effect on the activation of Omicron-specific naive B cells, potentially explaining the heavy immune imprinting in mRNA-vaccinated individuals6-8. These findings delineate the evolving antibody response to the antigenic shift of Omicron from XBB to JN.1 and highlight the importance of developing the JN.1 lineage, especially KP.2- and KP.3-based vaccine boosters.

Dynamic O-GlcNAcylation governs long-range chromatin interactions in V(D)J recombination during early B-cell development

V(D)J recombination secures the production of functional immunoglobulin (Ig) genes and antibody diversity during the early stages of B-cell development through long-distance interactions mediated by cis-regulatory elements and trans-acting factors. O-GlcNAcylation is a dynamic and reversible posttranslational modification of nuclear and cytoplasmic proteins that regulates various protein functions, including DNA-binding affinity and protein-protein interactions. However, the effects of O-GlcNAcylation on proteins involved in V(D)J recombination remain largely unknown. To elucidate this relationship, we downregulated O-GlcNAcylation in a mouse model by administering an O-GlcNAc inhibitor or restricting the consumption of a regular diet. Interestingly, the inhibition of O-GlcNAcylation in mice severely impaired Ig heavy-chain (IgH) gene rearrangement. We identified several factors crucial for V(D)J recombination, including YY1, CTCF, SMC1, and SMC3, as direct targets of O-GlcNAc modification. Importantly, O-GlcNAcylation regulates the physical interaction between SMC1 and SMC3 and the DNA-binding patterns of YY1 at the IgH gene locus. Moreover, O-GlcNAc inhibition downregulated DDX5 protein expression, affecting the functional association of CTCF with its DNA-binding sites at the IgH locus. Our results showed that locus contraction and long-range interactions throughout the IgH locus are disrupted in a manner dependent on the cellular O-GlcNAc level. In this study, we established that V(D)J recombination relies on the O-GlcNAc status of stage-specific proteins during early B-cell development and identified O-GlcNAc-dependent mechanisms as new regulatory components for the development of a diverse antibody repertoire.

The Type 1 Diabetes T Cell Receptor and B Cell Receptor Repository in the AIRR Data Commons: a practical guide for access, use and contributions through the Type 1 Diabetes AIRR Consortium

Human molecular genetics has brought incredible insights into the variants that confer risk for the development of tissue-specific autoimmune diseases, including type 1 diabetes. The hallmark cell-mediated immune destruction that is characteristic of type 1 diabetes is closely linked with risk conferred by the HLA class II gene locus, in combination with a broad array of additional candidate genes influencing islet-resident beta cells within the pancreas, as well as function, phenotype and trafficking of immune cells to tissues. In addition to the well-studied germline SNP variants, there are critical contributions conferred by T cell receptor (TCR) and B cell receptor (BCR) genes that undergo somatic recombination to yield the Adaptive Immune Receptor Repertoire (AIRR) responsible for autoimmunity in type 1 diabetes. We therefore created the T1D TCR/BCR Repository (The Type 1 Diabetes T Cell Receptor and B Cell Receptor Repository) to study these highly variable and dynamic gene rearrangements. In addition to processed TCR and BCR sequences, the T1D TCR/BCR Repository includes detailed metadata (e.g. participant demographics, disease-associated parameters and tissue type). We introduce the Type 1 Diabetes AIRR Consortium goals and outline methods to use and deposit data to this comprehensive repository. Our ultimate goal is to facilitate research community access to rich, carefully annotated immune AIRR datasets to enable new scientific inquiry and insight into the natural history and pathogenesis of type 1 diabetes.

sc-SPLASH provides ultra-efficient reference-free discovery in barcoded single-cell sequencing

Typical high-throughput single-cell RNA-sequencing (scRNA-seq) analyses are primarily conducted by (pseudo)alignment, through the lens of annotated gene models, and aimed at detecting differential gene expression. This misses diversity generated by other mechanisms that diversify the transcriptome such as splicing and V(D)J recombination, and is blind to sequences missing from imperfect reference genomes. Here, we present sc-SPLASH, a highly efficient pipeline that extends our SPLASH framework for statistics-first, reference-free discovery to barcoded scRNA-seq (10x Chromium) and spatial transcriptomics (10x Visium); we also provide its optimized module for preprocessing and k-mer counting in barcoded data, BKC, as a standalone tool. sc-SPLASH rediscovers known biology including V(D)J recombination and cell-type-specific alternative splicing in human and trans-splicing in tunicate (Ciona) and when applied to spatial datasets, detects sequence variation including tumor-specific somatic mutation. In sponge (Spongilla) and tunicate (Ciona), we uncover secreted repeat proteins expressed in immune-type cells and regulated during development; the sponge genes were absent from the reference assembly. sc-SPLASH provides a powerful alternative tool for exploring transcriptomes that is applicable to the breadth of life's diversity.

Ultrasensitive allele inference from immune repertoire sequencing data with MiXCR

Allelic variability in the adaptive immune receptor loci, which harbor the gene segments that encode B cell and T cell receptors (BCR/TCR), is of critical importance for immune responses to pathogens and vaccines. Adaptive immune receptor repertoire sequencing (AIRR-seq) has become widespread in immunology research making it the most readily available source of information about allelic diversity in immunoglobulin (IG) and T cell receptor (TR) loci. Here, we present a novel algorithm for extrasensitive and specific variable (V) and joining (J) gene allele inference, allowing the reconstruction of individual high-quality gene segment libraries. The approach can be applied for inferring allelic variants from peripheral blood lymphocyte BCR and TCR repertoire sequencing data, including hypermutated isotype-switched BCR sequences, thus allowing high-throughput novel allele discovery from a wide variety of existing data sets. The developed algorithm is a part of the MiXCR software. We demonstrate the accuracy of this approach using AIRR-seq paired with long-read genomic sequencing data, comparing it to a widely used algorithm, TIgGER. We applied the algorithm to a large set of IG heavy chain (IGH) AIRR-seq data from 450 donors of ancestrally diverse population groups, and to the largest reported full-length TCR alpha and beta chain (TRA and TRB) AIRR-seq data set, representing 134 individuals. This allowed us to assess the genetic diversity within the IGH, TRA, and TRB loci in different populations and to establish a database of alleles of V and J genes inferred from AIRR-seq data and their population frequencies with free public access through VDJ.online database.

B Cell Receptor Repertoire Analysis of the CD21lo B Cell Compartment in Healthy Individuals, Patients With Sjögren's Disease, and Patients With Radiographic Axial Spondyloarthritis

B cells with low or absent expression of CD21 (CD21lo B cells) gained attention due to their expansion in the peripheral blood of patients with immune-mediated, rheumatic diseases. This is not only observed in typical autoimmune diseases like systemic lupus erythematosus and Sjögren's disease (SjD) but also in radiographic axial spondyloarthritis (r-axSpA), which is considered an autoinflammatory disease. To gain more insight into the origins of the heterogeneous CD21lo B-cell population, and its relation to the plasmablast (PB) compartment, we profiled the B-cell-receptor (BCR) repertoire in CD27- and CD27+ fractions of CD21lo B cells and early PBs using next-generation sequencing. Populations were sorted from peripheral blood of healthy individuals, SjD patients, and r-axSpA patients (n = 10 for each group). In healthy individuals and both patient groups, our findings indicate that CD27-CD21lo B cells, which exhibit few mutations in their BCR, may develop into CD27+CD21lo B cells and PBs, both marked by considerably more mutations. Given the known expansion of circulating CD27-CD21lo B cells in SjD and r-axSpA patients and clonal relationships with both CD27+CD21lo B cells and early PBs, these cells might actively contribute to (pathological) immune responses in rheumatic diseases with autoimmune and/or autoinflammatory characteristics.

Pathobiont-driven antibody sialylation through IL-10 undermines vaccination

The pathobiont Staphylococcus aureus (Sa) induces nonprotective antibody imprints that underlie ineffective staphylococcal vaccination. However, the mechanism by which Sa modifies antibody activity is not clear. Herein, we demonstrate that IL-10 is the decisive factor that abrogates antibody protection in mice. Sa-induced B10 cells drive antigen-specific vaccine suppression that affects both recalled and de novo developed B cells. Released IL-10 promotes STAT3 binding upstream of the gene encoding sialyltransferase ST3gal4 and increases its expression by B cells, leading to hyper-α2,3sialylation of antibodies and loss of protective activity. IL-10 enhances α2,3sialylation on cell-wall-associated IsdB, IsdA, and MntC antibodies along with suppression of the respective Sa vaccines. Consistent with mouse findings, human anti-Sa antibodies as well as anti-pseudomonal antibodies from cystic fibrosis subjects (high IL-10) are hypersialylated, compared with anti-Streptococcus pyogenes and pseudomonal antibodies from normal individuals. Overall, we demonstrate a pathobiont-centric mechanism that modulates antibody glycosylation through IL-10, leading to loss of staphylococcal vaccine efficacy.

The carbonyl nucleobase adduct M3Ade is a potent antigen for adaptive polyclonal MR1-restricted T cells

The major histocompatibility complex (MHC) class I-related molecule MHC-class-I-related protein 1 (MR1) presents metabolites to distinct MR1-restricted T cell subsets, including mucosal-associated invariant T (MAIT) and MR1T cells. However, self-reactive MR1T cells and the nature of recognized antigens remain underexplored. Here, we report a cell endogenous carbonyl adduct of adenine (8-(9H-purin-6-yl)-2-oxa-8-azabicyclo[3.3.1]nona-3,6-diene-4,6-dicarbaldehyde [M3Ade]) sequestered in the A' pocket of MR1. M3Ade induced in vitro MR1-mediated stimulation of MR1T cell clones that bound MR1-M3Ade tetramers. MR1-M3Ade tetramers identified heterogeneous MR1-reactive T cells ex vivo in healthy donors, individuals with acute myeloid leukemia, and tumor-infiltrating lymphocytes from non-small cell lung adenocarcinoma and hepatocarcinoma. These cells displayed phenotypic, transcriptional, and functional diversity at distinct differentiation stages, indicating their adaptive nature. They were also polyclonal, with some preferential T cell receptor (TCRαβ) pair usage. Thus, M3Ade is an MR1-presented self-metabolite that enables stimulation and tracking of human-MR1T cells from blood and tissue, aiding our understanding of their roles in health and disease.

Impaired development of memory B cells and antibody responses in humans and mice deficient in PD-1 signaling

T follicular helper (Tfh) cells abundantly express the immunoreceptor programmed cell death protein 1 (PD-1), and the impact of PD-1 deficiency on antibody (Ab)-mediated immunity in mice is associated with compromised Tfh cell functions. Here, we revisited the role of the PD-1-PD-L1 axis on Ab-mediated immunity. Individuals with inherited PD-1 or PD-L1 deficiency had fewer memory B cells and impaired Ab responses, similar to Pdcd1-/- and Cd274-/-Pdcd1lg2-/- mice. PD-1, PD-L1, or both could be detected on the surface of human naive B cells following in vitro activation. PD-1- or PD-L1-deficient B cells had reduced expression of the transcriptional regulator c-Myc and c-Myc-target genes in vivo, and PD-1 deficiency or neutralization of PD-1 or PD-L1 impeded c-Myc expression and Ab production in human B cells isolated in vitro. Furthermore, B cell-specific deletion of Pdcd1 prevented the physiological accumulation of memory B cells in mice. Thus, PD-1 shapes optimal B cell memory and Ab-mediated immunity through B cell-intrinsic and B cell-extrinsic mechanisms, suggesting that B cell dysregulation contributes to infectious and autoimmune complications following anti-PD-1-PD-L1 immunotherapy.

Defining a highly conserved B cell epitope in the receptor binding motif of SARS-CoV-2 spike glycoprotein

SARS-CoV-2 mRNA vaccines induce robust and persistent germinal centre (GC) B cell responses in humans. It remains unclear how the continuous evolution of the virus impacts the breadth of the induced GC B cell response. Using ultrasound-guided fine needle aspiration, we examined draining lymph nodes of nine healthy adults following bivalent booster immunization. We show that 77.8% of the B cell clones in the GC expressed as representative monoclonal antibodies recognized the spike protein, with a third (37.8%) of these targeting the receptor binding domain (RBD). Strikingly, only one RBD-targeting mAb, mAb-52, neutralized all tested SARS-CoV-2 strains, including the recent KP.2 variant. mAb-52 utilizes the IGHV3-66 public clonotype, protects hamsters challenged against the EG.5.1 variant and targets the class I/II RBD epitope, closely mimicking the binding footprint of ACE2. Finally, we show that the remarkable breadth of mAb-52 is due to the somatic hypermutations accumulated within vaccine-induced GC reaction.

Group 2 innate lymphoid cells are a non-redundant source of interleukin-5 required for development and function of murine B1 cells

Tissue-resident immune cells, such as innate lymphoid cells, mediate protective or detrimental immune responses at barrier surfaces. Upon activation by stromal or epithelial cell-derived alarmins, group 2 innate lymphoid cells (ILC2s) are a rapid source of type 2 cytokines, such as IL-5. However, due to the overlap in effector functions, it remains unresolved whether ILC2s are an essential component of the type 2 response or whether their function can be compensated by other cells, such as T cells. Here we show a non-redundant role of ILC2s in supporting the development and function of B1 cells. We demonstrate that B1 cells fail to develop properly in the absence of ILC2s and identify the IL-33 receptor on ILC2s as an essential cell-intrinsic regulator of IL-5 production. Further, conditional deletion of Il5 in ILC2s results in defective B1 cell development and immunoglobulin production. Consequently, B1 cells with phosphatidylcholine specific B cell receptor rearrangements are diminished in ILC2-deficient mice. Thus, our data establish an essential function of ILC2s in supporting B1 cells and antibody production at barrier surfaces.

Immunomolecular and reactivity landscapes of gut IgA subclasses in homeostasis and inflammatory bowel disease

The human gut includes plasma cells (PCs) expressing immunoglobulin A1 (IgA1) or IgA2, two structurally distinct IgA subclasses with elusive regulation, function, and reactivity. We show here that intestinal IgA1+ and IgA2+ PCs co-emerged early in life, comparably accumulated somatic mutations, and were enriched within short-lived CD19+ and long-lived CD19- PC subsets, respectively. IgA2+ PCs were extensively clonally related to IgA1+ PCs and a subset of them presumably emerged from IgA1+ precursors. Of note, secretory IgA1 (SIgA1) and SIgA2 dually coated a large fraction of mucus-embedded bacteria, including Akkermansia muciniphila. Disruption of homeostasis by inflammatory bowel disease (IBD) was associated with an increase in actively proliferating IgA1+ plasmablasts, a depletion in long-lived IgA2+ PCs, and increased SIgA1+SIgA2+ gut microbiota. Such increase featured enhanced IgA1 reactivity to pathobionts, including Escherichia coli, combined with depletion of beneficial A. muciniphila. Thus, gut IgA1 and IgA2 emerge from clonally related PCs and show unique changes in both frequency and reactivity in IBD.

Applying phylogenetic methods for species delimitation to distinguish B-cell clonal families

The adaptive immune system generates a diverse array of B-cell receptors through the processes of V(D)J recombination and somatic hypermutation. B-cell receptors that bind to an antigen will undergo clonal expansion, creating a Darwinian evolutionary dynamic within individuals. A key step in studying these dynamics is to identify sequences derived from the same ancestral V(D)J recombination event (i.e. a clonal family). There are a number of widely used methods for accomplishing this task but a major limitation of all of them is that they rely, at least in part, on the ability to map sequences to a germline reference set. This requirement is particularly problematic in non-model systems where we often know little about the germline allelic diversity in the study population. Recognizing that delimiting B-cell clonal families is analogous to delimiting species from single locus data, we propose a novel strategy of reconstructing the phylogenetic tree of all B-cell sequences in a sample and using a popular species delimitation method, multi-rate Poisson Tree Processes (mPTP), to delimit clonal families. Using extensive simulations, we show that not only does this phylogenetically explicit approach perform well for the purpose of delimiting clonal families when no reference allele set is available, it performs similarly to state-of-the-art techniques developed specifically for B-cell data even when we have a complete reference allele set. Additionally, our analysis of an empirical dataset shows that mPTP performs similarly to leading methods in the field. These findings demonstrate the utility of using off-the-shelf phylogenetic techniques for analyzing B-cell clonal dynamics in non-model systems, and suggests that phylogenetic inference techniques may be potentially combined with mapping based approaches for even more robust inferences, even in model systems.

Immune responses in checkpoint myocarditis across heart, blood and tumour

Immune checkpoint inhibitors are widely used anticancer therapies1 that can cause morbid and potentially fatal immune-related adverse events such as immune-related myocarditis (irMyocarditis)2-5. The pathogenesis of irMyocarditis and its relationship to antitumour immunity remain poorly understood. Here we sought to define immune responses in heart, tumour and blood in patients with irMyocarditis by leveraging single-cell RNA sequencing coupled with T cell receptor (TCR) sequencing, microscopy and proteomics analyses of samples from 28 patients with irMyocarditis and 41 unaffected individuals. Analyses of 84,576 cardiac cells by single-cell RNA sequencing combined with multiplexed microscopy demonstrated increased frequencies and co-localization of cytotoxic T cells, conventional dendritic cells and inflammatory fibroblasts in irMyocarditis heart tissue. Analyses of 366,066 blood cells revealed decreased frequencies of plasmacytoid dendritic cells, conventional dendritic cells and B lineage cells but an increased frequency of other mononuclear phagocytes in irMyocarditis. Fifty-two heart-expanded TCR clones from eight patients did not recognize the putative cardiac autoantigens α-myosin, troponin I or troponin T. Additionally, TCRs enriched in heart tissue were largely nonoverlapping with those enriched in paired tumour tissue. The presence of heart-expanded TCRs in a cycling blood CD8 T cell population was associated with fatal irMyocarditis case status. Collectively, these findings highlight crucial biology driving irMyocarditis and identify putative biomarkers.

Clonal landscape of autoantibody-secreting plasmablasts in COVID-19 patients

Whereas severe COVID-19 is often associated with elevated autoantibody titers, the underlying mechanism behind their generation has remained unclear. Here we report clonal composition and diversity of autoantibodies in humoral response to SARS-CoV-2. Immunoglobulin repertoire analysis and characterization of plasmablast-derived monoclonal antibodies uncovered clonal expansion of plasmablasts producing cardiolipin (CL)-reactive autoantibodies. Half of the expanded CL-reactive clones exhibited strong binding to SARS-CoV-2 antigens. One such clone, CoV1804, was reactive to both CL and viral nucleocapsid (N), and further showed anti-nucleolar activity in human cells. Notably, antibodies sharing genetic features with CoV1804 were identified in COVID-19 patient-derived immunoglobulins, thereby constituting a novel public antibody. These public autoantibodies had numerous mutations that unambiguously enhanced anti-N reactivity, when causing fluctuations in anti-CL reactivity along with the acquisition of additional self-reactivities, such as anti-nucleolar activity, in the progeny. Thus, potentially CL-reactive precursors may have developed multiple self-reactivities through clonal selection, expansion, and somatic hypermutation driven by viral antigens. Our results revealed the nature of autoantibody production during COVID-19 and provided novel insights into the origin of virus-induced autoantibodies.

Omicron-specific ultra-potent SARS-CoV-2 neutralizing antibodies targeting the N1/N2 loop of Spike N-terminal domain

A multitude of functional mutations continue to emerge on the N-terminal domain (NTD) of the spike protein in SARS-CoV-2 Omicron subvariants. Understanding the immunogenicity of Omicron NTD and the properties of antibodies elicited by it is crucial for comprehending the impact of NTD mutations on viral fitness and guiding vaccine design. In this study, we find that most of NTD-targeting antibodies isolated from individuals with BA.5/BF.7 breakthrough infection (BTI) are ancestral (wild-type or WT)-reactive and non-neutralizing. Surprisingly, we identified five ultra-potent neutralizing antibodies (NAbs) that can only bind to Omicron but not WT NTD. Structural analysis revealed that they bind to a unique epitope on the N1/N2 loop of NTD and interact with the receptor-binding domain (RBD) via the light chain. These Omicron-specific NAbs achieve neutralization through ACE2 competition and blockage of ACE2-mediated S1 shedding. However, BA.2.86 and BA.2.87.1, which carry insertions or deletions on the N1/N2 loop, can evade these antibodies. Together, we provided a detailed map of the NTD-targeting antibody repertoire in the post-Omicron era, demonstrating their vulnerability to NTD mutations enabled by its evolutionary flexibility, despite their potent neutralization. These results revealed the function of the indels in the NTD of BA.2.86/JN.1 sublineage in evading neutralizing antibodies and highlighted the importance of considering the immunogenicity of NTD in vaccine design.

Utilizing machine learning and hemagglutinin sequences to identify likely hosts of influenza H3Nx viruses

Influenza is a disease that represents both a public health and agricultural risk with pandemic potential. Among the subtypes of influenza A virus, H3 influenza virus can infect many avian and mammalian species and is therefore a virus of interest to human and veterinary public health. The primary goal of this study was to train and validate classifiers for the identification of the most likely host species using the hemagglutinin gene segment of H3 viruses. A five-step process was implemented, which included training four machine learning classifiers, testing the classifiers on the validation dataset, and further exploration of the best-performing model on three additional datasets. The gradient boosting machine classifier showed the highest host-classification accuracy with a 98.0 % (95 % CI [97.01, 98.73]) correct classification rate on an independent validation dataset. The classifications were further analyzed using the predicted probability score which highlighted sequences of particular interest. These sequences were both correctly and incorrectly classified sequences that showed considerable predicted probability for multiple hosts. This showed the potential of using these classifiers for rapid sequence classification and highlighting sequences of interest. Additionally, the classifiers were tested on a separate swine dataset composed of H3N2 sequences from 1998 to 2003 from the United States of America, and a separate canine dataset composed of canine H3N2 sequences of avian origin. These two datasets were utilized to look at the applications of predicted probability and host convergence over time. Lastly, the classifiers were used on an independent dataset of environmental sequences to explore the host identification of environmental sequences. The results of these classifiers show the potential for machine learning to be used as a host identification technique for viruses of unknown origin on a species-specific level.

Multiplexed Antibody Sequencing and Profiling of the Human Hemagglutinin-specific Memory B Cell Response following Influenza Vaccination

Influenza virus is a highly contagious respiratory pathogen causing between 9.4 and 41 million infections per year in the United States in the last decade. Annual vaccination is recommended by the World Health Organization, with the goal to reduce influenza severity and transmission. Ag-specific single B cell sequencing methodologies have opened up new avenues into the dissection of the Ab response to influenza virus. The improvement of these methodologies is pivotal to reduce the associated costs and optimize the operational workflow and throughput, especially in the context of multiple samples. In this study, PBMCs and serum samples were collected longitudinally from eight influenza vaccinees either vaccinated yearly for four consecutive influenza seasons or once for one season. Following the serological and B cell profiling of their polyclonal Ab response to a panel of historical, recent, and next-generation influenza vaccine hemagglutinin (HA) and virus strains, a single multiplexed Ag-specific single B cell sequencing run allowed to capture HA-specific memory B cells that were analyzed for preferential Ig H chain/L chain pairing, isotype/subclass usage, and the presence of public BCR clonotypes across participants. Binding and functional profiles of representative private and public clonotypes confirmed their HA specificity, and their overall binding and functional activity were consistent with those observed at the polyclonal level. Collectively, this high-resolution and multiplexed Ab repertoire analysis demonstrated the validity of this optimized methodology in capturing Ag-specific BCR clonotypes, even in the context of a rare B cell population, such as in the case of the peripheral Ag-specific memory B cells.

Bulk T cell repertoire sequencing (TCR-Seq) is a powerful technology for understanding inflammation-mediated diseases

Multiple alterations in the T cell repertoire were identified in many chronic inflammatory diseases such as inflammatory bowel disease, multiple sclerosis, and rheumatoid arthritis, suggesting that T cells might, directly or indirectly, be implicated in these pathologies. This has sparked a deep interest in characterizing disease-associated T cell clonotypes as well as in identifying and quantifying their contribution to the pathophysiology of different autoimmune and inflammation-mediated diseases. Bulk T cell repertoire sequencing (TCR-Seq) has emerged as a powerful method to profile the T cell repertoire of a sample in a high throughput fashion. Given the increasing utilization of TCR-Seq, we aimed here to provide a comprehensive, up-to-date review of the method, its extensions, and its ability to investigate chronic and autoimmune diseases. Specifically, we started by introducing the immunological basis of TCR repertoire generation and features, followed by discussing different experimental approach to perform TCR-Seq, then we describe different methods and frameworks for analyzing the generated datasets. Subsequently, different experimental techniques for investigating the antigenicity of T cell clonotypes are described. Lastly, we discuss recent studies that utilized TCR-Seq to understand different inflammation-mediated diseases, discuss fallbacks of the technology and potential future directions to overcome current limitations.

Aberrant zonal recycling of germinal center B cells impairs appropriate selection in lupus

Autoimmune diseases such as lupus are characterized by polyclonal B cell activation, leading to the production of autoantibodies. The mechanism leading to B cell dysregulation is unclear; however, the defect may lie in selection within germinal centers (GCs). GC B cells cycle between proliferation and mutation in the dark zone and selection in the light zone (LZ). Temporal assessment of GCs from mice with either persistent infection or lupus showed an accumulation of LZ B cells. Yet, only in lupus, GC B cells exhibited reduced proliferation and progressive loss of MYC and FOXO1, which regulate zonal recycling and differentiation. As lupus progressed, decreased mutational frequency and repertoire diversity were associated with reduced responsiveness to CD40 signaling, despite accumulation of plasma cells. Collectively, these findings suggest that lupus disease progression coincides with an intrinsic defect in LZ B cell signaling, altering the zonal recycling, selection, and differentiation of autoreactive B cells.

SARS-CoV-2 infection elucidates features of pregnancy-specific immunity

Pregnancy is a risk factor for increased severity of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and other respiratory infections, but the mechanisms underlying this risk are poorly understood. To gain insight into the role of pregnancy in modulating immune responses at baseline and upon SARS-CoV-2 infection, we collected peripheral blood mononuclear cells and plasma from 226 women, including 152 pregnant individuals and 74 non-pregnant women. We find that SARS-CoV-2 infection is associated with altered T cell responses in pregnant women, including a clonal expansion of CD4-expressing CD8+ T cells, diminished interferon responses, and profound suppression of monocyte function. We also identify shifts in cytokine and chemokine levels in the sera of pregnant individuals, including a robust increase of interleukin-27, known to drive T cell exhaustion. Our findings reveal nuanced pregnancy-associated immune responses, which may contribute to the increased susceptibility of pregnant individuals to viral respiratory infection.

Diverse priming outcomes under conditions of very rare precursor B cells

Rare B cells can have special pathogen-recognition features giving them the potential to make outsized contributions to protective immunity. However, rare naive B cells infrequently participate in immune responses. We investigated how germline-targeting vaccine antigen delivery and adjuvant selection affect priming of exceptionally rare BG18-like HIV broadly neutralizing antibody-precursor B cells (~1 in 50 million) in non-human primates. Only escalating dose (ED) priming immunization using the saponin adjuvant SMNP elicited detectable BG18-like cells in germinal centers (GCs). All groups had strong GC responses, but only ED+SMNP and bolus+SMNP induced BG18-like memory B cells in >50% of animals. One group had vaccine-specific GC responses equivalent to ED+SMNP, but BG18-like memory B cells were rarely detected. Following homologous boosting, BG18-like memory B cells were more frequent in a bolus priming group, but had lower somatic hypermutation and affinities. This outcome was inversely associated with post-prime antibody titers, suggesting antibody feedback can significantly influence rare precursor B cell responses.

Integrated single-cell transcriptome and TCR profiles of hepatocellular carcinoma highlight the convergence on interferon signaling during immunotherapy

BACKGROUND

Despite the success of immune checkpoint inhibitor (ICI)-based combination therapies in hepatocellular carcinoma (HCC), its effectiveness remains confined to a subset of patients. The development of reliable, predictive markers is important for accurate patient stratification and further mechanistic understanding of therapy response.

METHODS

We comprehensively analyzed paired single-cell RNA transcriptome and T-cell repertoire profiles from 14 HCC ascites samples, collected from 7 patients before and after treatment with the combination of sintilimab (anti-PD-1) and bevacizumab (anti-VEGF).

RESULTS

We identify a widespread convergence on interferon (IFN) signaling across various immune cell lineages in treatment-responsive patients with HCC, indicating a common transcriptional state transition in the immune microenvironment linked to immunotherapy response in HCC. Strong IFN signaling marks CD8+ T cells with larger clonal expansion and enhanced cytotoxicity, macrophages toward M1-like polarization and strong T-cell recruitment ability, dendritic cells with increased antigen presentation capacity, as well as highly cytotoxic natural killer cells and activated B cells. By translating our finding to cohorts of patients with HCC, we demonstrate the specificity of IFN-signaling in the prognosis of patients with HCC and its ability to predict immunotherapy response.

CONCLUSIONS

This study provides a unique single-cell resource with clonal and longitudinal resolution during ICI therapy and reveals IFN signaling as a biomarker of immunotherapy response in HCC, suggesting a beneficial effect by combining IFN inducers with ICIs for patients with HCC.

Defective germinal center selection results in persistence of self-reactive B cells from the primary to the secondary repertoire in Primary Antiphospholipid Syndrome

Primary antiphospholipid syndrome (PAPS) is a life-threatening clotting disorder mediated by pathogenic autoantibodies. Here we dissect the origin of self-reactive B cells in human PAPS using peripheral blood and bone marrow of patients with triple-positive PAPS via combined single-cell RNA sequencing, B cell receptors (BCR) repertoire profiling, CITEseq analysis and single cell immortalization. We find that antiphospholipid (aPL)-specific B cells are present in the naive compartment, polyreactive, and derived from the natural repertoire. Furthermore, B cells with aPL specificities are not eliminated in patients with PAPS, persist until the memory and long-lived plasma cell stages, likely after defective germinal center selection, while becoming less polyreactive. Lastly, compared with the non-PAPS cells, PAPS B cells exhibit distinct IFN and APRIL signature as well as dysregulated mTORC1 and MYC pathways. Our findings may thus elucidate the survival mechanisms of these autoreactive B cells and suggest potential therapeutic targets for the treatment of PAPS.

A B cell screen against endogenous retroviruses identifies glycan-reactive IgM that recognizes a broad array of enveloped viruses

Endogenous retroviruses (ERVs), comprising a substantial portion of the vertebrate genome, are remnants of ancient genetic invaders. ERVs with near-intact coding potential reactivate in B cell-deficient mice. To study how B cells contribute to host anti-ERV immunity, we used an antigen-baiting strategy to enrich B cells reactive to ERV surface antigens. We identified ERV-reactive B-1 cells expressing germline-encoded natural IgM antibodies in naïve mice, the level of which further increases upon innate immune sensor stimulation. B cell receptor repertoire profiling of ERV-reactive B-1 cells revealed increased usage of the Igh VH gene that gives rise to glycan-specific antibodies targeting terminal N-acetylglucosamine moieties on ERV glycoproteins, which further engage the complement pathway to mediate anti-ERV responses. These same antibodies also recognize glycoproteins of other enveloped viruses but not self-proteins. These results reveal an innate antiviral mechanism of germline-encoded antibodies with broad reactivity to enveloped viruses, which constitutes a natural antibody repertoire capable of preventing the emergence of infectious ERVs.

Discovery and characterization of a pan-betacoronavirus S2-binding antibody

The continued emergence of deadly human coronaviruses from animal reservoirs highlights the need for pan-coronavirus interventions for effective pandemic preparedness. Here, using linking B cell receptor to antigen specificity through sequencing (LIBRA-seq), we report a panel of 50 coronavirus antibodies isolated from human B cells. Of these, 54043-5 was shown to bind the S2 subunit of spike proteins from alpha-, beta-, and deltacoronaviruses. A cryoelectron microscopy (cryo-EM) structure of 54043-5 bound to the prefusion S2 subunit of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike defined an epitope at the apex of S2 that is highly conserved among betacoronaviruses. Although non-neutralizing, 54043-5 induced Fc-dependent antiviral responses in vitro, including antibody-dependent cellular cytotoxicity (ADCC) and antibody-dependent cellular phagocytosis (ADCP). In murine SARS-CoV-2 challenge studies, protection against disease was observed after introduction of Leu234Ala, Leu235Ala, and Pro329Gly (LALA-PG) substitutions in the Fc region of 54043-5. Together, these data provide new insights into the protective mechanisms of non-neutralizing antibodies and define a broadly conserved epitope within the S2 subunit.

Antigen specificity of clonally-enriched CD8+ T cells in multiple sclerosis

CD8+ T cells are the dominant lymphocyte population in multiple sclerosis (MS) lesions where they are highly clonally expanded. The clonal identity, function, and antigen specificity of CD8+ T cells in MS are not well understood. Here we report a comprehensive single-cell RNA-seq and T cell receptor (TCR)-seq analysis of the cerebrospinal fluid (CSF) and blood from a cohort of treatment-naïve MS patients and control participants. A small subset of highly expanded and activated CD8+ T cells were enriched in the CSF in MS that displayed high activation, cytotoxicity and tissue-homing transcriptional profiles. Using a combination of unbiased and targeted antigen discovery approaches, MS-derived CD8+ T cell clonotypes recognizing Epstein-Barr virus (EBV) antigens and multiple novel mimotopes were identified. These findings shed vital insight into the role of CD8+ T cells in MS and pave the way towards disease biomarkers and therapeutic targets.

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.

mRNA-based influenza vaccine expands breadth of B cell response in humans

Eliciting broad and durable antibody responses against rapidly evolving pathogens like influenza viruses remains a formidable challenge1,2. The germinal center (GC) reaction enables the immune system to generate broad, high-affinity, and durable antibody responses to vaccination3-5. mRNA-based severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) vaccines induce persistent GC B cell responses in humans6-9. Whether an mRNA-based influenza vaccine could induce a superior GC response in humans compared to the conventional inactivated influenza virus vaccine remains unclear. We assessed B cell responses in peripheral blood and draining lymph nodes in cohorts receiving the inactivated or mRNA-based quadrivalent seasonal influenza vaccine. Participants receiving the mRNA-based vaccine produced more robust plasmablast responses and higher antibody titers to H1N1 and H3N2 influenza A viruses and comparable antibody titers against influenza B virus strains. Importantly, mRNA-based vaccination stimulated robust recall B cell responses characterized by sustained GC reactions that lasted at least 26 weeks post-vaccination in three of six participants analyzed. In addition to promoting the maturation of responding B cell clones, these sustained GC reactions resulted in enhanced engagement of low-frequency pre-existing memory B cells, expanding the landscape of vaccine-elicited B cell clones. This translated to expansion of the serological repertoire and increased breadth of serum antibody responses. These findings reveal an important role for the induction of persistent GC responses to influenza vaccination in humans to broaden the repertoire of vaccine-induced antibodies.

Intestinal newborn regulatory B cell antibodies modulate microbiota communities

The role of immunoglobulins produced by IL-10-producing regulatory B cells remains unknown. We found that a particular newborn regulatory B cell population (nBreg) negatively regulates the production of immunoglobulin M (IgM) via IL-10 in an autocrine manner, limiting the intensity of the polyreactive antibody response following innate activation. Based on nBreg scRNA-seq signature, we identify these cells and their repertoire in fetal and neonatal intestinal tissues. By characterizing 205 monoclonal antibodies cloned from intestinal nBreg, we show that newborn germline-encoded antibodies display reactivity against bacteria representing six different phyla of the early microbiota. nBreg-derived antibodies can influence the diversity and the cooperation between members of early microbial communities, at least in part by modulating energy metabolism. These results collectively suggest that nBreg populations help facilitate early-life microbiome establishment and shed light on the paradoxical activities of regulatory B cells in early life.

Use of 3M-052-AF with Alum adjuvant in HIV trimer vaccine induces human autologous neutralizing antibodies

Stabilized trimers preserving the native-like HIV envelope structure may be key components of a preventive HIV vaccine regimen to induce broadly neutralizing antibodies (bnAbs). We evaluated trimeric BG505 SOSIP.664 gp140 formulated with a novel TLR7/8 signaling adjuvant, 3M-052-AF/Alum, for safety, adjuvant dose-finding, and immunogenicity in a first-in-healthy adult (n = 17), randomized, and placebo-controlled trial (HVTN 137A). The vaccine regimen appeared safe. Robust, trimer-specific antibody, and B cell and CD4+ T cell responses emerged after vaccination. Five vaccinees developed serum autologous tier 2 nAbs (ID50 titer, 1:28-1:8647) after two to three doses targeting C3/V5 and/or V1/V2/V3 Env regions by electron microscopy and mutated pseudovirus-based neutralization analyses. Trimer-specific, B cell-derived monoclonal antibody activities confirmed these results and showed weak heterologous neutralization in the strongest responder. Our findings demonstrate the clinical utility of the 3M-052-AF/Alum adjuvant and support further improvements of trimer-based Env immunogens to focus responses on multiple broad nAb epitopes.

Discovery, recognized antigenic structures, and evolution of cross-serotype broadly neutralizing antibodies from porcine B-cell repertoires against foot-and-mouth disease virus

It is a great challenge to isolate the broadly neutralizing antibodies (bnAbs) against foot-and-mouth disease virus (FMDV) due to its existence as seven distinct serotypes without cross-protection. Here, by vaccination of pig with FMDV serotypes O and A whole virus antigens, we obtained 10 bnAbs against serotypes O, A and/or Asia1 by dissecting 216 common clonotypes of two serotypes O and A specific porcine B-cell receptor (BCR) gene repertoires containing total 12720 B cell clones, indicating the induction of cross-serotype bnAbs after sequential vaccination with serotypes O and A antigens. The majority of porcine bnAbs (9/10) were derived from terminally differentiated B cells of different clonal lineages, which convergently targeted the conserved "RGDL" motif on structural protein VP1 of FMDV by mimicking receptor recognition to inhibit viral attachment to cells. Cryo-EM complex structures revealed that the other bnAb pOA-2 specifically targets a novel inter-pentamer antigen structure surrounding the viral three-fold axis, with a highly conserved determinant at residue 68 on VP2. This unique binding pattern enabled cross-serotype neutralization by destabilizing the viral particle. The evolutionary analysis of pOA-2 demonstrated its origin from an intermediate B-cell, emphasizing the crucial role of somatic hypermutations (SHMs) in balancing the breadth and potency of neutralization. However, excessive SHMs may deviate from the trajectory of broad neutralization. This study provides a strategy to uncover bnAbs against highly mutable pathogens and the cross-serotype antigenic structures to explore broadly protective FMDV vaccine.

Full-length single-cell BCR sequencing paired with RNA sequencing reveals convergent responses to pneumococcal vaccination

Single-cell RNA sequencing (scRNA-seq) can resolve transcriptional features from individual cells, but scRNA-seq techniques capable of resolving the variable regions of B cell receptors (BCRs) remain limited, especially from widely-used 3'-barcoded libraries. Here, we report a method that can recover paired, full-length variable region sequences of BCRs from 3'-barcoded scRNA-seq libraries. We first verify this method (B3E-seq) can produce accurate, full-length BCR sequences. We then apply this method to profile B cell responses elicited against the capsular polysaccharide of Streptococcus pneumoniae serotype 3 (ST3) by glycoconjugate vaccines in five infant rhesus macaques. We identify BCR features associated with specificity for the ST3 antigen which are present in multiple vaccinated monkeys, indicating a convergent response to vaccination. These results demonstrate the utility of our method to resolve key features of the B cell repertoire and profile antigen-specific responses elicited by vaccination.

An unbiased comparison of immunoglobulin sequence aligners

Adaptive Immune Receptor Repertoire sequencing (AIRR-seq) is critical for our understanding of the adaptive immune system's dynamics in health and disease. Reliable analysis of AIRR-seq data depends on accurate rearranged immunoglobulin (Ig) sequence alignment. Various Ig sequence aligners exist, but there is no unified benchmarking standard representing the complexities of AIRR-seq data, obscuring objective comparisons of aligners across tasks. Here, we introduce GenAIRR, a modular simulation framework for generating Ig sequences alongside their ground truths. GenAIRR realistically simulates the intricacies of V(D)J recombination, somatic hypermutation, and an array of sequence corruptions. We comprehensively assessed prominent Ig sequence aligners across various metrics, unveiling unique performance characteristics for each aligner. The GenAIRR-produced datasets, combined with the proposed rigorous evaluation criteria, establish a solid basis for unbiased benchmarking of immunogenetics computational tools. It sets up the ground for further improving the crucial task of Ig sequence alignment, ultimately enhancing our understanding of adaptive immunity.

Expansion of tumor-reactive CD8+ T cell clonotypes occurs in the spleen in response to immune checkpoint blockade

Immune checkpoint blockade (ICB) enhances T cell responses against cancer, leading to long-term survival in a fraction of patients. CD8+ T cell differentiation in response to chronic antigen stimulation is highly complex, and it remains unclear precisely which T cell differentiation states at which anatomic sites are critical for the response to ICB. We identified an intermediate-exhausted population in the white pulp of the spleen that underwent substantial expansion in response to ICB and gave rise to tumor-infiltrating clonotypes. Increased systemic antigen redirected differentiation of this population toward a more circulatory exhausted KLR state, whereas a lack of cross-presented tumor antigen reduced its differentiation in the spleen. An analogous population of exhausted KLR CD8+ T cells in human blood samples exhibited diminished tumor-trafficking ability. Collectively, our data demonstrate the critical role of antigen density within the spleen for the differentiation and expansion of T cell clonotypes in response to ICB.

Rocuronium-specific antibodies drive perioperative anaphylaxis but can also function as reversal agents in preclinical models

Neuromuscular blocking agents (NMBAs) relax skeletal muscles to facilitate surgeries and ease intubation but can lead to adverse reactions, including complications because of postoperative residual neuromuscular blockade (rNMB) and, in rare cases, anaphylaxis. Both adverse reactions vary between types of NMBAs, with rocuronium, a widely used nondepolarizing NMBA, inducing one of the longest rNMB durations and highest anaphylaxis incidences. rNMB induced by rocuronium can be reversed by the synthetic γ-cyclodextrin sugammadex. However, in rare cases, sugammadex can provoke anaphylaxis. Thus, additional therapeutic options are needed. Rocuronium-induced anaphylaxis is proposed to rely on preexisting rocuronium-binding antibodies. To understand the pathogenesis of rocuronium-induced anaphylaxis and to identify potential therapeutics, we investigated the memory B cell antibody repertoire of patients with suspected hypersensitivity to rocuronium. We identified polyclonal antibody repertoires with a high diversity among V(D)J genes without evidence of clonal groups. When recombinantly expressed, these antibodies demonstrated specificity and low affinity for rocuronium without cross-reactivity for other NMBAs. Moreover, when these antibodies were expressed as human immunoglobulin E (IgE), they triggered human mast cell activation and passive systemic anaphylaxis in transgenic mice, although their affinities were insufficient to serve as reversal agents. Rocuronium-specific, high-affinity antibodies were thus isolated from rocuronium-immunized mice. The highest-affinity antibody was able to reverse rocuronium-induced neuromuscular blockade in nonhuman primates with kinetics comparable to that of sugammadex. Together, these data support the hypothesis that antibodies cause anaphylactic reactions to rocuronium and pave the way for improved diagnostics and neuromuscular blockade reversal agents.

Interleukin-2 immunotherapy reveals human regulatory T cell subsets with distinct functional and tissue-homing characteristics

Due to its stimulatory potential for immunomodulatory CD4+ regulatory T (Treg) cells, low-dose interleukin-2 (IL-2) immunotherapy has gained considerable attention for the treatment of autoimmune diseases. In this investigator-initiated single-arm non-placebo-controlled phase-2 clinical trial of low-dose IL-2 immunotherapy in systemic lupus erythematosus (SLE) patients, we generated a comprehensive atlas of in vivo human immune responses to low-dose IL-2. We performed an in-depth study of circulating and cutaneous immune cells by imaging mass cytometry, high-parameter flow cytometry, transcriptomics, and targeted serum proteomics. Low-dose IL-2 stimulated various circulating immune cells, including Treg cells with a skin-homing phenotype that appeared in the skin of SLE patients in close interaction with endothelial cells. Analysis of surface proteins and transcriptomes revealed different IL-2-driven Treg cell activation programs, including gut-homing CD38+, skin-homing HLA-DR+, and highly proliferative inflammation-homing CD38+ HLA-DR+ Treg cells. Collectively, these data define the distinct human Treg cell subsets that are responsive to IL-2 immunotherapy.

Natural malaria infection elicits rare but potent neutralizing antibodies to the blood-stage antigen RH5

Plasmodium falciparum reticulocyte-binding protein homolog 5 (RH5) is the most advanced blood-stage malaria vaccine candidate and is being evaluated for efficacy in endemic regions, emphasizing the need to study the underlying antibody response to RH5 during natural infection, which could augment or counteract responses to vaccination. Here, we found that RH5-reactive B cells were rare, and circulating immunoglobulin G (IgG) responses to RH5 were short-lived in malaria-exposed Malian individuals, despite repeated infections over multiple years. RH5-specific monoclonal antibodies isolated from eight malaria-exposed individuals mostly targeted non-neutralizing epitopes, in contrast to antibodies isolated from five RH5-vaccinated, malaria-naive UK individuals. However, MAD8-151 and MAD8-502, isolated from two malaria-exposed Malian individuals, were among the most potent neutralizers out of 186 antibodies from both cohorts and targeted the same epitopes as the most potent vaccine-induced antibodies. These results suggest that natural malaria infection may boost RH5-vaccine-induced responses and provide a clear strategy for the development of next-generation RH5 vaccines.