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

Single-cell RNA sequencing reveals diverse B cell phenotypes in patients with anti-NMDAR encephalitis

BACKGROUNDS

Anti-N-methyl-D-aspartate receptor encephalitis (NMDAR-E) is a severe autoimmune disorder characterized by prominent psychiatric symptoms. Although the role of NMDAR antibodies in the disease has been extensively studied, the phenotype of B cell subsets is still not fully understood.

METHODS

We utilized single-cell RNA sequencing, single-cell B cell receptor sequencing (scBCR-seq), bulk BCR sequencing, flow cytometry, and enzyme-linked immunosorbent assay to analyze samples from both NMDAR-E patients and control individuals.

RESULTS

The cerebrospinal fluid (CSF) of NMDAR-E patients showed significantly increased B cell counts, predominantly memory B (Bm) cells. CSF Bm cells in NMDAR-E patients exhibited upregulated expression of differential expression genes (DEGs) associated with immune regulatory function (TNFRSF13B and ITGB1), whereas peripheral B cells upregulated DEGs related to antigen presentation. Additionally, NMDAR-E patients displayed higher levels of IgD- CD27- double negative (DN) cells and DN3 cells in peripheral blood (PB). In vitro, DN1 cell subsets from NMDAR-E patients differentiated into DN2 and DN3 cells, while CD27+ and/or IgD+ B cells (non-DN) differentiated into antibody-secreting cells (ASCs) and DN cells. NR1-IgG antibodies were found in B cell culture supernatants from patients. Differential expression of B cell IGHV genes in CSF and PB of NMDAR-E patients suggests potential antigen class switching.

CONCLUSION

B cell subpopulations in the CSF and PB of NMDAR-E patients exhibit distinct compositions and transcriptomic features. In vitro, non-DN cells from NMDAR-E can differentiate into DN cells and ASCs, potentially producing NR1-IgG antibodies. Further research is necessary to investigate the potential contribution of DN cell subpopulations to NR1-IgG antibody production.

Paired single-B-cell transcriptomics and receptor sequencing reveal activation states and clonal signatures that characterize B cells in acute myeloid leukemia

BACKGROUND

Acute myeloid leukemia (AML) is associated with a dismal prognosis. Immune checkpoint blockade (ICB) to induce antitumor activity in AML patients has yielded mixed results. Despite the pivotal role of B cells in antitumor immunity, a comprehensive assessment of B lymphocytes within AML's immunological microenvironment along with their interaction with ICB remains rather constrained.

METHODS

We performed an extensive analysis that involved paired single-cell RNA and B-cell receptor (BCR) sequencing on 52 bone marrow aspirate samples. These samples included 6 from healthy bone marrow donors (normal), 24 from newly diagnosed AML patients (NewlyDx), and 22 from 8 relapsed or refractory AML patients (RelRef), who underwent assessment both before and after azacitidine/nivolumab treatment.

RESULTS

We delineated nine distinct subtypes of B cell lineage in the bone marrow. AML patients exhibited reduced nascent B cell subgroups but increased differentiated B cells compared with healthy controls. The limited diversity of BCR profiles and extensive somatic hypermutation indicated antigen-driven affinity maturation within the tumor microenvironment of RelRef patients. We established a strong connection between the activation or stress status of naïve and memory B cells, as indicated by AP-1 activity, and their differentiation state. Remarkably, atypical memory B cells functioned as specialized antigen-presenting cells closely interacting with AML malignant cells, correlating with AML stemness and worse clinical outcomes. In the AML microenvironment, plasma cells demonstrated advanced differentiation and heightened activity. Notably, the clinical response to ICB was associated with B cell clonal expansion and plasma cell function.

CONCLUSIONS

Our findings establish a comprehensive framework for profiling the phenotypic diversity of the B cell lineage in AML patients, while also assessing the implications of immunotherapy. This will serve as a valuable guide for future inquiries into AML treatment strategies.

B cell repertoire sequencing of HIV-1 pediatric elite-neutralizers identifies multiple broadly neutralizing antibody clonotypes

Introduction

A limited subset of HIV-1 infected adult individuals typically after at least 2-3 years of chronic infection, develop broadly neutralizing antibodies (bnAbs), suggesting that highly conserved neutralizing epitopes on the HIV-1 envelope glycoprotein are difficult for B cell receptors to effectively target, during natural infection. Recent studies have shown the evolution of bnAbs in HIV-1 infected infants.

Methods

We used bulk BCR sequencing (BCR-seq) to profile the B cell receptors from longitudinal samples (3 time points) collected from a rare pair of antiretroviralnaïve, HIV-1 infected pediatric monozygotic twins (AIIMS_329 and AIIMS_330) who displayed elite plasma neutralizing activity against HIV-1.

Results

BCR-seq of both twins revealed convergent antibody characteristics including V-gene use, CDRH3 lengths and somatic hypermutation (SHM). Further, antibody clonotypes with genetic features similar to highly potent bnAbs isolated from adults showed ongoing development in donor AIIMS_330 but not in AIIMS_329, corroborating our earlier findings based on plasma bnAbs responses. An increase in SHM was observed in sequences of the IgA isotype from AIIMS_330.

Discussion

This study suggests that children living with chronic HIV-1 can develop clonotypes of HIV-1 bnAbs against multiple envelope epitopes similar to those isolated from adults, highlighting that such B cells could be steered to elicit bnAbs responses through vaccines aimed to induce bnAbs against HIV-1 in a broad range of people including children.

Lack of affinity signature for germinal center cells that have initiated plasma cell differentiation

Long-lived plasma cells (PCs) secrete antibodies that can provide sustained immunity against infection. High-affinity cells are proposed to preferentially select into this compartment, potentiating the immune response. We used single-cell RNA-seq to track the germinal center (GC) development of Ighg2A10 B cells, specific for the Plasmodium falciparum circumsporozoite protein (PfCSP). Following immunization with Plasmodium sporozoites, we identified 3 populations of cells in the GC light zone (LZ). One LZ population expressed a gene signature associated with the initiation of PC differentiation and readily formed PCs in vitro. The estimated affinity of these pre-PC B cells was indistinguishable from that of LZ cells that remained in the GC. This remained true when high- or low-avidity recombinant PfCSP proteins were used as immunogens. These findings suggest that the initiation of PC development occurs via an affinity-independent process.

Ultrahigh frequencies of peripherally matured LGI1- and CASPR2-reactive B cells characterize the cerebrospinal fluid in autoimmune encephalitis

Intrathecal synthesis of central nervous system (CNS)-reactive autoantibodies is observed across patients with autoimmune encephalitis (AE), who show multiple residual neurobehavioral deficits and relapses despite immunotherapies. We leveraged two common forms of AE, mediated by leucine-rich glioma inactivated-1 (LGI1) and contactin-associated protein-like 2 (CASPR2) antibodies, as human models to comprehensively reconstruct and profile cerebrospinal fluid (CSF) B cell receptor (BCR) characteristics. We hypothesized that the resultant observations would both inform the observed therapeutic gap and determine the contribution of intrathecal maturation to pathogenic B cell lineages. From the CSF of three patients, 381 cognate-paired IgG BCRs were isolated by cell sorting and scRNA-seq, and 166 expressed as monoclonal antibodies (mAbs). Sixty-two percent of mAbs from singleton BCRs reacted with either LGI1 or CASPR2 and, strikingly, this rose to 100% of cells in clonal groups with ≥4 members. These autoantigen-reactivities were more concentrated within antibody-secreting cells (ASCs) versus B cells (P < 0.0001), and both these cell types were more differentiated than LGI1- and CASPR2-unreactive counterparts. Despite greater differentiation, autoantigen-reactive cells had acquired few mutations intrathecally and showed minimal variation in autoantigen affinities within clonal expansions. Also, limited CSF T cell receptor clonality was observed. In contrast, a comparison of germline-encoded BCRs versus the founder intrathecal clone revealed marked gains in both affinity and mutational distances (P = 0.004 and P < 0.0001, respectively). Taken together, in patients with LGI1 and CASPR2 antibody encephalitis, our results identify CSF as a compartment with a remarkably high frequency of clonally expanded autoantigen-reactive ASCs whose BCR maturity appears dominantly acquired outside the CNS.

Systematic evaluation of B-cell clonal family inference approaches

The reconstruction of clonal families (CFs) in B-cell receptor (BCR) repertoire analysis is a crucial step to understand the adaptive immune system and how it responds to antigens. The BCR repertoire of an individual is formed throughout life and is diverse due to several factors such as gene recombination and somatic hypermutation. The use of Adaptive Immune Receptor Repertoire sequencing (AIRR-seq) using next generation sequencing enabled the generation of full BCR repertoires that also include rare CFs. The reconstruction of CFs from AIRR-seq data is challenging and several approaches have been developed to solve this problem. Currently, most methods use the heavy chain (HC) only, as it is more variable than the light chain (LC). CF reconstruction options include the definition of appropriate sequence similarity measures, the use of shared mutations among sequences, and the possibility of reconstruction without preliminary clustering based on V- and J-gene annotation. In this study, we aimed to systematically evaluate different approaches for CF reconstruction and to determine their impact on various outcome measures such as the number of CFs derived, the size of the CFs, and the accuracy of the reconstruction. The methods were compared to each other and to a method that groups sequences based on identical junction sequences and another method that only determines subclones. We found that after accounting for data set variability, in particular sequencing depth and mutation load, the reconstruction approach has an impact on part of the outcome measures, including the number of CFs. Simulations indicate that unique junctions and subclones should not be used as substitutes for CF and that more complex methods do not outperform simpler methods. Also, we conclude that different approaches differ in their ability to correctly reconstruct CFs when not considering the LC and to identify shared CFs. The results showed the effect of different approaches on the reconstruction of CFs and highlighted the importance of choosing an appropriate method.

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

Pregnancy is a risk factor for increased severity of SARS-CoV-2 and other respiratory infections. The mechanisms underlying this risk have not been well-established, partly due to a limited understanding of how pregnancy shapes immune responses. To gain insight into the role of pregnancy in modulating immune responses at steady state and upon perturbation, we collected peripheral blood mononuclear cells (PBMC), plasma, and stool from 226 women, including 152 pregnant individuals (n = 96 with SARS-CoV-2 infection and n = 56 healthy controls) and 74 non-pregnant women (n = 55 with SARS-CoV-2 and n = 19 healthy controls). We found that SARS-CoV-2 infection was associated with altered T cell responses in pregnant compared to non-pregnant women. Differences included a lower percentage of memory T cells, a distinct clonal expansion of CD4-expressing CD8 + T cells, and the enhanced expression of T cell exhaustion markers, such as programmed cell death-1 (PD-1) and T cell immunoglobulin and mucin domain-3 (Tim-3), in pregnant women. We identified additional evidence of immune dysfunction in severely and critically ill pregnant women, including a lack of expected elevation in regulatory T cell (Treg) levels, diminished interferon responses, and profound suppression of monocyte function. Consistent with earlier data, we found maternal obesity was also associated with altered immune responses to SARS-CoV-2 infection, including enhanced production of inflammatory cytokines by T cells. Certain gut bacterial species were altered in pregnancy and upon SARS-CoV-2 infection in pregnant individuals compared to non-pregnant women. Shifts in cytokine and chemokine levels were also identified in the sera of pregnant individuals, most notably a robust increase of interleukin-27 (IL-27), a cytokine known to drive T cell exhaustion, in the pregnant uninfected control group compared to all non-pregnant groups. IL-27 levels were also significantly higher in uninfected pregnant controls compared to pregnant SARS-CoV-2-infected individuals. Using two different preclinical mouse models of inflammation-induced fetal demise and respiratory influenza viral infection, we found that enhanced IL-27 protects developing fetuses from maternal inflammation but renders adult female mice vulnerable to viral infection. These combined findings from human and murine studies reveal nuanced pregnancy-associated immune responses, suggesting mechanisms underlying the increased susceptibility of pregnant individuals to viral respiratory infections.

Computational Analysis of B-Cell Receptor (BCR) Immune Repertoires with Abalign

The widespread application of high-throughput sequencing technology has generated massive sequences of B-cell receptor (BCR) immune repertoires. Computational analysis of these data has gained significant attention due to the increasing importance of immunotherapy and precision medicine. It not only reveals the diversity and dynamic changes in immune responses, contributing to the study of associated diseases, but also provides valuable information for immunodiagnostics and drug development. Recently, we introduced a BCR-specific multiple sequence alignment (MSA) method along with a comprehensive platform software called Abalign, which stands out as an excellent choice for analyzing BCR immune repertoires due to its unique high-throughput processing capability. It offers ultra-fast MSA functionality and a wide range of analytical features, including BCR/antibody extraction, clonal grouping, lineage tree construction, mutation profiling, diversity statistics, VJ gene assignment, antibody humanization, and more. Importantly, users can perform these analyses using the graphical user interface without any programming skills or scripts. In this article, we present a series of protocols that integrate Abalign's analysis modules into a cohesive workflow. This step-by-step workflow provides detailed instructions for software installation, data preparation, and comprehensive analysis of BCR immune repertoires. This workflow facilitates the efficient acquisition of comprehensive results in profiling BCR immune repertoires, offering insights into the impacts of infectious diseases, allergies, autoimmune disorders, tumor immunology, and antibody drugs. Abalign is freely available at http://cao.labshare.cn/abalign/. © 2024 Wiley Periodicals LLC. Basic Protocol 1: Resource preparation Basic Protocol 2: Analyzing BCR immune repertoires Support Protocol 1: Aiding antibody humanization Support Protocol 2: Constructing B-cell lineage trees Alternate Protocol: Running with Linux command line Basic Protocol 3: Comparing BCR immune repertoires.

Multimodal Single-Cell Sequencing of B Cells in Primary Sjögren's Syndrome

OBJECTIVE

B cells are important in the pathogenesis of primary Sjögren's syndrome (pSS). Patients positive for Sjögren's syndrome antigen A/Sjögren syndrome antigen B (SSA/SSB) autoantibodies are more prone to systemic disease manifestations and adverse outcomes. We aimed to determine the role of B cell composition, gene expression, and B cell receptor usage in pSS subgroups stratified for SSA/SSB antibodies.

METHODS

Over 230,000 B cells were isolated from peripheral blood of patients with pSS (n = 6 SSA-, n = 8 SSA+ single positive and n = 10 SSA/SSB+ double positive) and four healthy controls and processed for single-cell RNA sequencing (scRNA-seq) and single-cell variable, diversity, and joining (VDJ) gene sequencing (scVDJ-seq).

RESULTS

We show that SSA/SSB+ patients present the highest and lowest proportion of naïve and memory B cells, respectively, and the highest up-regulation of interferon-induced genes across all B cell subtypes. Differential usage of IGHV showed that IGHV1-69 and IGHV4-30-4 were more often used in all pSS subgroups compared with controls. Memory B cells from SSA/SSB+ patients displayed a higher proportion of cells with unmutated VDJ transcripts compared with other pSS patient groups and controls, indicating altered somatic hypermutation processes. Comparison with previous studies revealed heterogeneous clonotype pools, with little overlap in CDR3 sequences. Joint analysis using scRNA-seq and scVDJ-seq data allowed unsupervised stratification of patients with pSS and identified novel parameters that correlated to disease manifestations and antibody status.

CONCLUSION

We describe heterogeneity and molecular characteristics in B cells from patients with pSS, providing clues to intrinsic differences in B cells that affect the phenotype and outcome and allowing stratification of patients with pSS at improved resolution.

Influence of pre-B cell receptor deficiency on the immunoglobulin repertoires in peripheral blood B cells before and after immunization

During B cell development, pre-B cell receptor (pre-BCR), comprising the immunoglobulin heavy chain (HC) and surrogate light chain (SLC), plays a crucial role. The expression of pre-BCR serves as a certification of HC quality, confirming its ability to associate with the SLC and light chain (LC). In mice lacking SLC, the absence of this quality control mechanism leads to a distorted repertoire of HCs in the spleen and bone marrow. In this study, we conducted a comparative analysis of the immunoglobulin gene repertoire in peripheral blood cells of both wild-type mice and pre-BCR-deficient mice. Our findings reveal differences not only in the μ HC repertoire but also in the α HC and κ LC repertoires of the pre-BCR-deficient mice. These results suggest that the pre-BCR-mediated quality check of HC influences the selection of class-switched HC and LC repertoires. To further explore the impact of pre-BCR deficiency, we immunized these mice with thymus-dependent antigens and compared the antigen-responding repertoires. Our observations indicate that the affinity maturation pathways remain consistent between wild-type mice and pre-BCR-deficient mice, albeit with variations in the degree of maturation.

B cell receptor repertoire abnormalities in autoimmune disease

B cells play a crucial role in the immune response and contribute to various autoimmune diseases. Recent studies have revealed abnormalities in the B cell receptor (BCR) repertoire of patients with autoimmune diseases, with distinct features observed among different diseases and B cell subsets. Classically, BCR repertoire was used as an identifier of distinct antigen-specific clonotypes, but the recent advancement of analyzing large-scale repertoire has enabled us to use it as a tool for characterizing cellular biology. In this review, we provide an overview of the BCR repertoire in autoimmune diseases incorporating insights from our latest research findings. In systemic lupus erythematosus (SLE), we observed a significant skew in the usage of VDJ genes, particularly in CD27+IgD+ unswitched memory B cells and plasmablasts. Notably, autoreactive clones within unswitched memory B cells were found to be increased and strongly associated with disease activity, underscoring the clinical significance of this subset. Similarly, various abnormalities in the BCR repertoire have been reported in other autoimmune diseases such as rheumatoid arthritis. Thus, BCR repertoire analysis holds potential for enhancing our understanding of the underlying mechanisms involved in autoimmune diseases. Moreover, it has the potential to predict treatment effects and identify therapeutic targets in autoimmune diseases.

nf-core/airrflow: an adaptive immune receptor repertoire analysis workflow employing the Immcantation framework

Adaptive Immune Receptor Repertoire sequencing (AIRR-seq) is a valuable experimental tool to study the immune state in health and following immune challenges such as infectious diseases, (auto)immune diseases, and cancer. Several tools have been developed to reconstruct B cell and T cell receptor sequences from AIRR-seq data and infer B and T cell clonal relationships. However, currently available tools offer limited parallelization across samples, scalability or portability to high-performance computing infrastructures. To address this need, we developed nf-core/airrflow, an end-to-end bulk and single-cell AIRR-seq processing workflow which integrates the Immcantation Framework following BCR and TCR sequencing data analysis best practices. The Immcantation Framework is a comprehensive toolset, which allows the processing of bulk and single-cell AIRR-seq data from raw read processing to clonal inference. nf-core/airrflow is written in Nextflow and is part of the nf-core project, which collects community contributed and curated Nextflow workflows for a wide variety of analysis tasks. We assessed the performance of nf-core/airrflow on simulated sequencing data with sequencing errors and show example results with real datasets. To demonstrate the applicability of nf-core/airrflow to the high-throughput processing of large AIRR-seq datasets, we validated and extended previously reported findings of convergent antibody responses to SARS-CoV-2 by analyzing 97 COVID-19 infected individuals and 99 healthy controls, including a mixture of bulk and single-cell sequencing datasets. Using this dataset, we extended the convergence findings to 20 additional subjects, highlighting the applicability of nf-core/airrflow to validate findings in small in-house cohorts with reanalysis of large publicly available AIRR datasets. nf-core/airrflow is available free of charge, under the MIT license on GitHub (https://github.com/nf-core/airrflow). Detailed documentation and example results are available on the nf-core website at (https://nf-co.re/airrflow).

The multilevel extensive diversity across the cynomolgus macaque captured by ultra-deep adaptive immune receptor repertoire sequencing

The extent to which AIRRs differ among and within individuals remains elusive. Via ultra-deep repertoire sequencing of 22 and 25 tissues in three cynomolgus macaques, respectively, we identified 84 and 114 novel IGHV and TRBV alleles, confirming 72 (85.71%) and 100 (87.72%) of them. The heterogeneous V gene usage patterns were influenced, in turn, by genetics, isotype (for BCRs only), tissue group, and tissue. A higher proportion of intragroup shared clones in the intestinal tissues than those in other tissues suggests a close intra-intestinal adaptive immunity network. Significantly higher mutation burdens in the public clones and the inter-tissue shared IgM and IgD clones indicate that they might target the shared antigens. This study reveals the extensive heterogeneity of the AIRRs at various levels and has broad fundamental and clinical implications. The data generated here will serve as an invaluable resource for future studies on adaptive immunity in health and diseases.

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

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

Anchor Clustering for million-scale immune repertoire sequencing data

BACKGROUND

The clustering of immune repertoire data is challenging due to the computational cost associated with a very large number of pairwise sequence comparisons. To overcome this limitation, we developed Anchor Clustering, an unsupervised clustering method designed to identify similar sequences from millions of antigen receptor gene sequences. First, a Point Packing algorithm is used to identify a set of maximally spaced anchor sequences. Then, the genetic distance of the remaining sequences to all anchor sequences is calculated and transformed into distance vectors. Finally, distance vectors are clustered using unsupervised clustering. This process is repeated iteratively until the resulting clusters are small enough so that pairwise distance comparisons can be performed.

RESULTS

Our results demonstrate that Anchor Clustering is faster than existing pairwise comparison clustering methods while providing similar clustering quality. With its flexible, memory-saving strategy, Anchor Clustering is capable of clustering millions of antigen receptor gene sequences in just a few minutes.

CONCLUSIONS

This method enables the meta-analysis of immune-repertoire data from different studies and could contribute to a more comprehensive understanding of the immune repertoire data space.

Discovery and Characterization of a Pan-betacoronavirus S2-binding antibody

Three coronaviruses have spilled over from animal reservoirs into the human population and caused deadly epidemics or pandemics. The continued emergence of coronaviruses highlights the need for pan-coronavirus interventions for effective pandemic preparedness. Here, using LIBRA-seq, we report a panel of 50 coronavirus antibodies isolated from human B cells. Of these antibodies, 54043-5 was shown to bind the S2 subunit of spike proteins from alpha-, beta-, and deltacoronaviruses. A cryo-EM structure of 54043-5 bound to the pre-fusion S2 subunit of the 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, including ADCC and 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.

Lupus-associated innate receptors drive extrafollicular evolution of autoreactive B cells

In systemic lupus erythematosus, recent findings highlight the extrafollicular (EF) pathway as prominent origin of autoantibody-secreting cells (ASCs). CD21loCD11c+ B cells, associated with aging, infection, and autoimmunity, are contributors to autoreactive EF ASCs but have an obscure developmental trajectory. To study EF kinetics of autoreactive B cell in tissue, we adoptively transferred WT and gene knockout B cell populations into the 564Igi mice - an autoreactive host enriched with autoantigens and T cell help. Time-stamped analyses revealed TLR7 dependence in early escape of peripheral B cell tolerance and establishment of a pre-ASC division program. We propose CD21lo cells as precursors to EF ASCs due to their elevated TLR7 sensitivity and proliferative nature. Blocking receptor function reversed CD21 loss and reduced effector cell generation, portraying CD21 as a differentiation initiator and a possible target for autoreactive B cell suppression. Repertoire analysis further delineated proto-autoreactive B cell selection and receptor evolution toward self-reactivity. This work elucidates receptor and clonal dynamics in EF development of autoreactive B cells, and establishes modular, native systems to probe mechanisms of autoreactivity.

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.

Dandelion uses the single-cell adaptive immune receptor repertoire to explore lymphocyte developmental origins

Assessment of single-cell gene expression (single-cell RNA sequencing) and adaptive immune receptor (AIR) sequencing (scVDJ-seq) has been invaluable in studying lymphocyte biology. Here we introduce Dandelion, a computational pipeline for scVDJ-seq analysis. It enables the application of standard V(D)J analysis workflows to single-cell datasets, delivering improved V(D)J contig annotation and the identification of nonproductive and partially spliced contigs. We devised a strategy to create an AIR feature space that can be used for both differential V(D)J usage analysis and pseudotime trajectory inference. The application of Dandelion improved the alignment of human thymic development trajectories of double-positive T cells to mature single-positive CD4/CD8 T cells, generating predictions of factors regulating lineage commitment. Dandelion analysis of other cell compartments provided insights into the origins of human B1 cells and ILC/NK cell development, illustrating the power of our approach. Dandelion is available at https://www.github.com/zktuong/dandelion .

Integrative single-cell analysis reveals distinct adaptive immune signatures in the cutaneous lesions of pemphigus

Pemphigus, an autoimmune bullous disease affecting the skin and mucosal membranes, is primarily driven by anti-desmoglein (Dsg) autoantibodies. However, the underlying immune mechanisms of this disease remain largely elusive. Here, we compile an unbiased atlas of immune cells in pemphigus cutaneous lesions at single-cell resolution. We reveal clonally expanded antibody-secreting cells (ASCs) that exhibit variable hypermutation and accumulation of IgG4 class-switching in their immunoglobulin genes. Importantly, pathogenic Dsg-specific ASCs are localized within pemphigus lesions and can evolve from both Dsg-autoreactive and non-binding precursors. We observe an altered distribution of CD4+ T cell subsets within pemphigus lesions, including an imbalance of Th17/Th2 cells. Significantly, we identify a distinct subpopulation of Th17 cells expressing CXCL13 and IL-21 within pemphigus lesions, implying its pivotal role in B cell recruitment and local production of autoantibodies. Furthermore, we characterize multiple clonally expanded CD8+ subpopulations, including effector GMZB+ and GMZK+ subsets with augmented cytotoxic activities, within pemphigus lesions. Chemokine-receptor mapping uncovers cell-type-specific signaling programs involved in the recruitment of T/B cells within pemphigus lesions. Our findings significantly contribute to advancing the understanding of the heterogeneous immune microenvironment and the pathogenesis of pemphigus cutaneous lesions, thereby providing valuable insights for potential therapeutic interventions in this disease.

IGH repertoire analysis at scale: deciphering the complexity of B cell infiltration and migration in esophageal squamous cell carcinoma

Tumor-infiltrating B-lineage cells have become predictors of prognosis and immunotherapy responses in various cancers. However, limited knowledge about their infiltration and migration patterns has hindered the understanding of their anti-tumor functions. Here, we examined the immunoglobulin heavy chain (IGH) repertoires in 496 multi-regional tumor, 107 normal tissue, and 48 metastatic lymph node samples obtained from 107 patients with esophageal squamous cell carcinoma (ESCC). Our study revealed higher IgG-type B-lineage cells infiltration in tumors than in healthy tissue, which was associated with improved patient outcomes. Genes such as ACTN1, COL6A5, and pathways like focal adhesion, which shapes the physical structure of tumors, could affect B-lineage cell infiltration. Notably, the IGH sequence was used as an identity-tag to monitor B cell migration, and their infiltration schema within the tumor were depicted based on our multi-regional tumor specimens. This analysis revealed an escalation in B cell clones overlapped between metastatic lymph nodes and tumors. Therefore, the Lymph Node Activation Index was defined, which could predict the outcomes of patients with lymph node metastasis. This research introduces a novel framework for probing B cell infiltration and migration within the tumor microenvironment using large-scale transcriptome data, while simultaneously providing fresh perspectives on B cell immunology within ESCC.

Follicular lymphoma evolves with a surmountable dependency on acquired glycosylation motifs in the B-cell receptor

ABSTRACT

An early event in the genesis of follicular lymphoma (FL) is the acquisition of new glycosylation motifs in the B-cell receptor (BCR) due to gene rearrangement and/or somatic hypermutation. These N-linked glycosylation motifs (N-motifs) contain mannose-terminated glycans and can interact with lectins in the tumor microenvironment, activating the tumor BCR pathway. N-motifs are stable during FL evolution, suggesting that FL tumor cells are dependent on them for their survival. Here, we investigated the dynamics and potential impact of N-motif prevalence in FL at the single-cell level across distinct tumor sites and over time in 17 patients. Although most patients had acquired at least 1 N-motif as an early event, we also found (1) cases without N-motifs in the heavy or light chains at any tumor site or time point and (2) cases with discordant N-motif patterns across different tumor sites. Inferring phylogenetic trees of the patients with discordant patterns, we observed that both N-motif-positive and N-motif-negative tumor subclones could be selected and expanded during tumor evolution. Comparing N-motif-positive with N-motif-negative tumor cells within a patient revealed higher expression of genes involved in the BCR pathway and inflammatory response, whereas tumor cells without N-motifs had higher activity of pathways involved in energy metabolism. In conclusion, although acquired N-motifs likely support FL pathogenesis through antigen-independent BCR signaling in most patients with FL, N-motif-negative tumor cells can also be selected and expanded and may depend more heavily on altered metabolism for competitive survival.

Atypical B cells consist of subsets with distinct functional profiles

Atypical B cells are a population of activated B cells that are commonly enriched in individuals with chronic immune activation but are also part of a normal immune response to infection or vaccination. To better define the role of atypical B cells in the human adaptive immune response, we performed single-cell sequencing of transcriptomes, cell surface markers, and B cell receptors in individuals with chronic exposure to the malaria parasite Plasmodium falciparum, a condition known to lead to accumulation of circulating atypical B cells. We identified three previously uncharacterized populations of atypical B cells with distinct transcriptional and functional profiles and observed marked differences among these three subsets in their ability to produce immunoglobulin G upon T-cell-dependent activation. Our findings help explain the conflicting observations in prior studies regarding the function of atypical B cells and highlight their different roles in the adaptive immune response in chronic inflammatory conditions.

Single-cell atlas of healthy human blood unveils age-related loss of NKG2C+GZMB-CD8+ memory T cells and accumulation of type 2 memory T cells

Extensive, large-scale single-cell profiling of healthy human blood at different ages is one of the critical pending tasks required to establish a framework for the systematic understanding of human aging. Here, using single-cell RNA/T cell receptor (TCR)/BCR-seq with protein feature barcoding, we profiled 317 samples from 166 healthy individuals aged 25-85 years old. From this, we generated a dataset from ∼2 million cells that described 55 subpopulations of blood immune cells. Twelve subpopulations changed with age, including the accumulation of GZMK+CD8+ T cells and HLA-DR+CD4+ T cells. In contrast to other T cell memory subsets, transcriptionally distinct NKG2C+GZMB-CD8+ T cells counterintuitively decreased with age. Furthermore, we found a concerted age-associated increase in type 2/interleukin (IL)4-expressing memory subpopulations across CD4+ and CD8+ T cell compartments (CCR4+CD8+ Tcm and Th2 CD4+ Tmem), suggesting a systematic functional shift in immune homeostasis with age. Our work provides novel insights into healthy human aging and a comprehensive annotated resource.

Spatial transcriptomics of B cell and T cell receptors reveals lymphocyte clonal dynamics

The spatial distribution of lymphocyte clones within tissues is critical to their development, selection, and expansion. We have developed spatial transcriptomics of variable, diversity, and joining (VDJ) sequences (Spatial VDJ), a method that maps B cell and T cell receptor sequences in human tissue sections. Spatial VDJ captures lymphocyte clones that match canonical B and T cell distributions and amplifies clonal sequences confirmed by orthogonal methods. We found spatial congruency between paired receptor chains, developed a computational framework to predict receptor pairs, and linked the expansion of distinct B cell clones to different tumor-associated gene expression programs. Spatial VDJ delineates B cell clonal diversity and lineage trajectories within their anatomical niche. Thus, Spatial VDJ captures lymphocyte spatial clonal architecture across tissues, providing a platform to harness clonal sequences for therapy.

Prior viral infection primes cross-reactive CD8+ T cells that respond to mouse heart allografts

Introduction

Significant evidence suggests a connection between transplant rejection and the presence of high levels of pre-existing memory T cells. Viral infection can elicit viral-specific memory T cells that cross-react with allo-MHC capable of driving allograft rejection in mice. Despite these advances, and despite their critical role in transplant rejection, a systematic study of allo-reactive memory T cells, their specificities, and the role of cross-reactivity with viral antigens has not been performed.

Methods

Here, we established a model to identify, isolate, and characterize cross-reactive T cells using Nur77 reporter mice (C57BL/6 background), which transiently express GFP exclusively upon TCR engagement. We infected Nur77 mice with lymphocytic choriomeningitis virus (LCMV-Armstrong) to generate a robust memory compartment, where quiescent LCMV-specific memory CD8+ T cells could be readily tracked with MHC tetramer staining. Then, we transplanted LCMV immune mice with allogeneic hearts and monitored expression of GFP within MHC-tetramer defined viral-specific T cells as an indicator of their ability to cross-react with alloantigens.

Results

Strikingly, prior LCMV infection significantly increased the kinetics and magnitude of rejection as well as CD8+ T cell recruitment into allogeneic, but not syngeneic, transplanted hearts, relative to non-infected controls. Interestingly, as early as day 1 after allogeneic heart transplant an average of ~8% of MHC-tetramer+ CD8+ T cells expressed GFP, in contrast to syngeneic heart transplants, where the frequency of viral-specific CD8+ T cells that were GFP+ was <1%. These data show that a significant percentage of viral-specific memory CD8+ T cells expressed T cell receptors that also recognized alloantigens in vivo. Notably, the frequency of cross-reactive CD8+ T cells differed depending upon the viral epitope. Further, TCR sequences derived from cross-reactive T cells harbored distinctive motifs that may provide insight into cross-reactivity and allo-specificity.

Discussion

In sum, we have established a mouse model to track viral-specific, allo-specific, and cross-reactive T cells; revealing that prior infection elicits substantial numbers of viral-specific T cells that cross-react to alloantigen, respond very early after transplant, and may promote rapid rejection.

Sustained Liver HBsAg Loss and Clonal T- and B-Cell Expansion upon Therapeutic DNA Vaccination Require Low HBsAg Levels

BACKGROUND

Suppression of HBV DNA, inhibition of HBV surface (HBsAg) production and therapeutic vaccination to reverse HBV-specific T-cell exhaustion in chronic HBV patients are likely required to achieve a functional cure. In the AAV-HBV mouse model, therapeutic vaccination can be effective in clearing HBV when HBsAg levels are low. Using a single-cell approach, we investigated the liver immune environment with different levels of HBsAg and sustained HBsAg loss through treatment with a GalNAc-HBV-siRNA followed by therapeutic vaccination.

METHODS

AAV-HBV-transduced C57BL/6 mice were treated with GalNAc-HBV-siRNA to lower HBsAg levels and then vaccinated using a DNA vaccine. We used single-cell RNA and V(D)J sequencing to understand liver immune microenvironment changes.

RESULTS

GalNAc-HBV-siRNA, followed by therapeutic vaccination, achieved sustained HBsAg loss in all mice. This was accompanied by CD4 follicular helper T-cell induction, polyclonal activation of CD8 T cells and clonal expansion of plasma cells that were responsible for antibody production.

CONCLUSIONS

This study provides novel insights into liver immune changes at the single-cell level, highlighting the correlation between induced reduction of HBsAg levels and clonal expansion of CD4, CD8 T cells and plasma cells in the liver upon HBV siRNA and subsequent therapeutic vaccination.

Yclon: Ultrafast clustering of B cell clones from high-throughput immunoglobulin repertoire sequencing data

MOTIVATION

The next-generation sequencing technologies have transformed our understanding of immunoglobulin (Ig) profiles in various immune states. Clonotyping, which groups Ig sequences into B cell clones, is crucial in investigating the diversity of repertoires and changes in antigen exposure. Despite its importance, there is no widely accepted method for clonotyping, and existing methods are computationally intensive for large sequencing datasets.

RESULTS

To address this challenge, we introduce YClon, a fast and efficient approach for clonotyping Ig repertoire data. YClon uses a hierarchical clustering approach, similar to other methods, to group Ig sequences into B cell clones in a highly sensitive and specific manner. Notably, our approach outperforms other methods by being more than 30 to 5000 times faster in processing the repertoires analyzed. Astonishingly, YClon can effortlessly handle up to 2 million Ig sequences on a standard laptop computer. This enables in-depth analysis of large and numerous antibody repertoires.

AVAILABILITY AND IMPLEMENTATION

YClon was implemented in Python3 and is freely available on GitHub.

Immune Responses in Checkpoint Myocarditis Across Heart, Blood, and Tumor

Immune checkpoint inhibitors (ICIs) are widely used anti-cancer therapies that can cause morbid and potentially fatal immune-related adverse events (irAEs). ICI-related myocarditis (irMyocarditis) is uncommon but has the highest mortality of any irAE. The pathogenesis of irMyocarditis and its relationship to anti-tumor immunity remain poorly understood. We sought to define immune responses in heart, tumor, and blood during irMyocarditis and identify biomarkers of clinical severity by leveraging single-cell (sc)RNA-seq coupled with T cell receptor (TCR) sequencing, microscopy, and proteomics analysis of 28 irMyocarditis patients and 23 controls. Our analysis of 284,360 cells from heart and blood specimens identified cytotoxic T cells, inflammatory macrophages, conventional dendritic cells (cDCs), and fibroblasts enriched in irMyocarditis heart tissue. Additionally, potentially targetable, pro-inflammatory transcriptional programs were upregulated across multiple cell types. TCR clones enriched in heart and paired tumor tissue were largely non-overlapping, suggesting distinct T cell responses within these tissues. We also identify the presence of cardiac-expanded TCRs in a circulating, cycling CD8 T cell population as a novel peripheral biomarker of fatality. Collectively, these findings highlight critical biology driving irMyocarditis and putative biomarkers for therapeutic intervention.

SARS-CoV-2 antibodies from children exhibit broad neutralization and belong to adult public clonotypes

From the beginning of the COVID-19 pandemic, children have exhibited different susceptibility to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, reinfection, and disease compared with adults. Motivated by the established significance of SARS-CoV-2-neutralizing antibodies in adults, here we characterize SARS-CoV-2-specific antibody repertoires in a young cohort of individuals aged from 5 months to 18 years old. Our results show that neutralizing antibodies in children possess similar genetic features compared to antibodies identified in adults, with multiple antibodies from children belonging to previously established public antibody clonotypes in adults. Notably, antibodies from children show potent neutralization of circulating SARS-CoV-2 variants that have cumulatively resulted in resistance to virtually all approved monoclonal antibody therapeutics. Our results show that children can rely on similar SARS-CoV-2 antibody neutralization mechanisms compared to adults and are an underutilized source for the discovery of effective antibody therapeutics to counteract the ever-evolving pandemic.

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), has been shown to be of critical importance for immune responses to pathogens and vaccines. In recent years, B cell and T cell receptor repertoire sequencing (Rep-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 in different populations. Here we present a novel algorithm for extra-sensitive and specific variable (V) and joining (J) gene allele inference and genotyping allowing 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 genotyping and novel allele discovery from a wide variety of existing datasets. The developed algorithm is a part of the MiXCR software ( https://mixcr.com ) and can be incorporated into any pipeline utilizing upstream processing with MiXCR. We demonstrate the accuracy of this approach using Rep-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) Rep-Seq data from 450 donors of ancestrally diverse population groups, and to the largest reported full-length TCR alpha and beta chain (TRA; TRB) Rep-Seq dataset, representing 134 individuals. This allowed us to assess the genetic diversity of genes within the IGH, TRA and TRB loci in different populations and demonstrate the connection between antibody repertoire gene usage and the number of allelic variants present in the population. Finally we established a database of allelic variants of V and J genes inferred from Rep-Seq data and their population frequencies with free public access at https://vdj.online .

Plasma cells in human pancreatic ductal adenocarcinoma secrete antibodies against self-antigens

Intratumoral B cell responses are associated with more favorable clinical outcomes in human pancreatic ductal adenocarcinoma (PDAC). However, the antigens driving these B cell responses are largely unknown. We sought to discover these antigens by using single-cell RNA sequencing (scRNA-Seq) and immunoglobulin (Ig) sequencing of tumor-infiltrating immune cells from 7 primary PDAC samples. We identified activated T and B cell responses and evidence of germinal center reactions. Ig sequencing identified plasma cell (PC) clones expressing isotype-switched and hypermutated Igs, suggesting the occurrence of T cell-dependent B cell responses. We assessed the reactivity of 41 recombinant antibodies that represented the products of 235 PCs and 12 B cells toward multiple cell lines and PDAC tissues and observed frequent staining of intracellular self-antigens. Three of these antigens were identified: the filamentous actin (F-actin), the nucleic protein RuvB like AAA ATPase 2 (RUVBL2), and the mitochondrial protein heat shock protein family D (Hsp60) member 1 (HSPD1). Antibody titers against F-actin and HSPD1 were substantially elevated in the plasma of patients with PDAC compared with healthy donors. Thus, PCs in PDAC produce autoantibodies reacting with intracellular self-antigens, which may result from promotion of preexisting, autoreactive B cell responses. These observations indicate the chronic inflammatory microenvironment of PDAC can support the adaptive immune response.

A combined adjuvant approach primes robust germinal center responses and humoral immunity in non-human primates

Adjuvants and antigen delivery kinetics can profoundly influence B cell responses and should be critically considered in rational vaccine design, particularly for difficult neutralizing antibody targets such as human immunodeficiency virus (HIV). Antigen kinetics can change depending on the delivery method. To promote extended immunogen bioavailability and to present antigen in a multivalent form, native-HIV Env trimers are modified with short phosphoserine peptide linkers that promote tight binding to aluminum hydroxide (pSer:alum). Here we explore the use of a combined adjuvant approach that incorporates pSer:alum-mediated antigen delivery with potent adjuvants (SMNP, 3M-052) in an extensive head-to-head comparison study with conventional alum to assess germinal center (GC) and humoral immune responses. Priming with pSer:alum plus SMNP induces additive effects that enhance the magnitude and persistence of GCs, which correlate with better GC-TFH cell help. Autologous HIV-neutralizing antibody titers are improved in SMNP-immunized animals after two immunizations. Over 9 months after priming immunization of pSer:alum with either SMNP or 3M-052, robust Env-specific bone marrow plasma cells (BM BPC) are observed. Furthermore, pSer-modification of Env trimer reduce targeting towards immunodominant non-neutralizing epitopes. The study shows that a combined adjuvant approach can augment humoral immunity by modulating immunodominance and shows promise for clinical translation.

Tissue-resident B cells orchestrate macrophage polarisation and function

B cells play a central role in humoral immunity but also have antibody-independent functions. Studies to date have focused on B cells in blood and secondary lymphoid organs but whether B cells reside in non-lymphoid organs (NLO) in homeostasis is unknown. Here we identify, using intravenous labeling and parabiosis, a bona-fide tissue-resident B cell population in lung, liver, kidney and urinary bladder, a substantial proportion of which are B-1a cells. Tissue-resident B cells are present in neonatal tissues and also in germ-free mice NLOs, albeit in lower numbers than in specific pathogen-free mice and following co-housing with 'pet-store' mice. They spatially co-localise with macrophages and regulate their polarization and function, promoting an anti-inflammatory phenotype, in-part via interleukin-10 production, with effects on bacterial clearance during urinary tract infection. Thus, our data reveal a critical role for tissue-resident B cells in determining the homeostatic 'inflammatory set-point' of myeloid cells, with important consequences for tissue immunity.

Mpox vaccine and infection-driven human immune signatures: an immunological analysis of an observational study

BACKGROUND

Monkeypox virus has recently infected more than 88 000 people, raising concerns about our preparedness against this emerging viral pathogen. Licensed and approved for mpox, the JYNNEOS vaccine has fewer side-effects than previous smallpox vaccines and has shown immunogenicity against monkeypox in animal models. This study aims to elucidate human immune responses to JYNNEOS vaccination compared with mpox-induced immunity.

METHODS

Peripheral blood mononuclear cells and sera were obtained from ten individuals vaccinated with one or two doses of JYNNEOS and six individuals diagnosed with monkeypox virus infection. Samples were obtained from seven individuals before vaccination to serve as a baseline. We examined the polyclonal serum (ELISA) and single B-cell (heavy chain gene and transcriptome data) antibody repertoires and T-cell responses (activation-induced marker and intracellular cytokine staining assays) induced by the JYNNEOS vaccine versus monkeypox virus infection.

FINDINGS

All participants were men between the ages of 21 and 60 years, except for one woman in the group of mpox-convalescent individuals, and none had previous orthopoxvirus exposure. All mpox cases were mild. Vaccinee samples were collected 6-33 days after the first dose and 5-40 days after the second dose. Mpox-convalescent samples were collected 20-102 days after infection. In vaccine recipients, gene-level plasmablast and antibody responses were negligible and sera displayed moderate binding to recombinant orthopoxviral proteins (A29L, A35R, E8L, A30L, A27L, A33R, B18R, and L1R) and native proteins from the 2022 monkeypox outbreak strain. By contrast, recent monkeypox virus infection (within 20-102 days) induced robust serum antibody responses to monkeypox virus proteins and to native monkeypox virus proteins from a viral isolate obtained during the 2022 outbreak. JYNNEOS vaccine recipients presented robust orthopoxviral CD4+ and CD8+ T-cell responses.

INTERPRETATION

Infection with monkeypox virus resulted in robust B-cell and T-cell responses, whereas immunisation with JYNNEOS elicited more robust T-cell responses. These data can help to inform vaccine design and policies for preventing mpox in humans.

FUNDING

National Cancer Institute (National Institutes of Health), National Institute of Allergy and Infectious Diseases (National Institutes of Health), and Icahn School of Medicine.

Cross-talk between Myeloid and B Cells Shapes the Distinct Microenvironments of Primary and Secondary Liver Cancer

The tumor microenvironment is distinctive in primary and secondary liver cancer. B cells represent an important component of immune infiltrates. Here, we demonstrated that B cells are an important regulator in hepatocellular carcinoma (HCC) and colorectal cancer liver metastasis (CRLM) microenvironments. B cells displayed distinct developmental trajectories in HCC and CRLM. Single-cell analysis revealed that IgG+ plasma cells preferentially accumulated in HCC, whereas IgA+ plasma cells were preferentially enriched in CRLM. Mechanistically, IgG+ plasma cells in HCC were recruited by tumor-associated macrophages via the CXCR3-CXCL10 axis, whereas IgA+ plasma cells in CRLM were recruited by metastatic tumor cells via CCR10-CCL28 signaling. Functionally, IgG+ plasma cells preferentially promoted protumorigenic macrophages formation in HCC, and IgA+ plasma cells preferentially induced granulocytic myeloid-derived suppressor cells activation in CRLM. Clinically, increased infiltration of IgG+ plasma cells and macrophages in HCC was correlated to worse survival, whereas increased intratumoral IgA+ plasma cells and neutrophils in CRLM indicated poor prognosis. Taken together, this study demonstrated plasma and myeloid cell-mediated immunosuppression in HCC and CRLM, suggesting that selectively modulating primary or secondary tumor-related immunosuppressive regulatory networks might reprogram the microenvironment and provide an immunotherapeutic strategy for treating liver cancer.

SIGNIFICANCE

The immunomodulatory patterns of tumor-infiltrating B cells are distinct in primary and secondary liver cancer, with plasma cells mediating important physiologic processes that drive cancer progression.

VH2+ Antigen-Experienced B Cells in the Cerebrospinal Fluid Are Expanded and Enriched in Pediatric Anti-NMDA Receptor Encephalitis

Pediatric and adult autoimmune encephalitis (AE) are often associated with Abs to the NR1 subunit of the N-methyl-d-aspartate (NMDA) receptor (NMDAR). Very little is known regarding the cerebrospinal fluid humoral immune profile and Ab genetics associated with pediatric anti-NMDAR-AE. Using a combination of cellular, molecular, and immunogenetics tools, we collected cerebrospinal fluid from pediatric subjects and generated 1) flow cytometry data to calculate the frequency of B cell subtypes in the cerebrospinal fluid of pediatric subjects with anti-NMDAR-AE and controls, 2) a panel of recombinant human Abs from a pediatric case of anti-NMDAR-AE that was refractory to treatment, and 3) a detailed analysis of the Ab genes that bound the NR1 subunit of the NMDAR. Ag-experienced B cells including memory cells, plasmablasts, and Ab-secreting cells were expanded in the pediatric anti-NMDAR-AE cohort, but not in the controls. These Ag-experienced B cells in the cerebrospinal fluid of a pediatric case of NMDAR-AE that was refractory to treatment had expanded use of variable H chain family 2 (VH2) genes with high somatic hypermutation that all bound to the NR1 subunit of the NMDAR. A CDR3 motif was identified in this refractory case that likely drove early stage activation and expansion of naive B cells to Ab-secreting cells, facilitating autoimmunity associated with pediatric anti-NMDAR-AE through the production of Abs that bind NR1. These features of humoral immune responses in the cerebrospinal fluid of pediatric anti-NMDAR-AE patients may be relevant for clinical diagnosis and treatment.

Immune landscape of tertiary lymphoid structures in hepatocellular carcinoma (HCC) treated with neoadjuvant immune checkpoint blockade

Neoadjuvant immunotherapy is thought to produce long-term remissions through induction of antitumor immune responses before removal of the primary tumor. Tertiary lymphoid structures (TLS), germinal center-like structures that can arise within tumors, may contribute to the establishment of immunological memory in this setting, but understanding of their role remains limited. Here, we investigated the contribution of TLS to antitumor immunity in hepatocellular carcinoma (HCC) treated with neoadjuvant immunotherapy. We found that neoadjuvant immunotherapy induced the formation of TLS, which were associated with superior pathologic response, improved relapse free survival, and expansion of the intratumoral T and B cell repertoire. While TLS in viable tumor displayed a highly active mature morphology, in areas of tumor regression we identified an involuted TLS morphology, which was characterized by dispersion of the B cell follicle and persistence of a T cell zone enriched for ongoing antigen presentation and T cell-mature dendritic cell interactions. Involuted TLS showed increased expression of T cell memory markers and expansion of CD8+ cytotoxic and tissue resident memory clonotypes. Collectively, these data reveal the circumstances of TLS dissolution and suggest a functional role for late-stage TLS as sites of T cell memory formation after elimination of viable tumor.

Multi-omics analysis of mucosal and systemic immunity to SARS-CoV-2 after birth

The dynamics of immunity to infection in infants remain obscure. Here, we used a multi-omics approach to perform a longitudinal analysis of immunity to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in infants and young children by analyzing blood samples and weekly nasal swabs collected before, during, and after infection with Omicron and non-Omicron variants. Infection stimulated robust antibody titers that, unlike in adults, showed no sign of decay for up to 300 days. Infants mounted a robust mucosal immune response characterized by inflammatory cytokines, interferon (IFN) α, and T helper (Th) 17 and neutrophil markers (interleukin [IL]-17, IL-8, and CXCL1). The immune response in blood was characterized by upregulation of activation markers on innate cells, no inflammatory cytokines, but several chemokines and IFNα. The latter correlated with viral load and expression of interferon-stimulated genes (ISGs) in myeloid cells measured by single-cell multi-omics. Together, these data provide a snapshot of immunity to infection during the initial weeks and months of life.

Distinct features of B cell receptors in neuromyelitis optica spectrum disorder among CNS inflammatory demyelinating diseases

BACKGROUND

Neuromyelitis optica spectrum disorder (NMOSD) stands out among CNS inflammatory demyelinating diseases (CIDDs) due to its unique disease characteristics, including severe clinical attacks with extensive lesions and its association with systemic autoimmune diseases. We aimed to investigate whether characteristics of B cell receptors (BCRs) differ between NMOSD and other CIDDs using high-throughput sequencing.

METHODS

From a prospective cohort, we recruited patients with CIDDs and categorized them based on the presence and type of autoantibodies: NMOSD with anti-aquaporin-4 antibodies, myelin oligodendrocyte glycoprotein antibody-associated disease (MOGAD) with anti-myelin oligodendrocyte glycoprotein antibodies, double-seronegative demyelinating disease (DSN), and healthy controls (HCs). The BCR features, including isotype class, clonality, somatic hypermutation (SHM), and the third complementarity-determining region (CDR3) length, were analyzed and compared among the different disease groups.

RESULTS

Blood samples from 33 patients with CIDDs (13 NMOSD, 12 MOGAD, and 8 DSN) and 34 HCs were investigated for BCR sequencing. Patients with NMOSD tended to have more activated BCR features compare to the other disease groups. They showed a lower proportion of unswitched isotypes (IgM and IgD) and a higher proportion of switched isotypes (IgG), increased clonality of BCRs, higher rates of SHM, and shorter lengths of CDR3. Notably, advanced age was identified as a clinical factor associated with these activated BCR features, including increased levels of clonality and SHM rates in the NMOSD group. Conversely, no such clinical factors were found to be associated with activated BCR features in the other CIDD groups.

CONCLUSIONS

NMOSD patients, among those with CIDDs, displayed the most pronounced B cell activation, characterized by higher levels of isotype class switching, clonality, SHM rates, and shorter CDR3 lengths. These findings suggest that B cell-mediated humoral immune responses and characteristics in NMOSD patients are distinct from those observed in the other CIDDs, including MOGAD. Age was identified as a clinical factor associated with BCR activation specifically in NMOSD, implying the significance of persistent B cell activation attributed to anti-aquaporin-4 antibodies, even in the absence of clinical relapses throughout an individual's lifetime.

IL-4 downregulates BCL6 to promote memory B cell selection in germinal centers

Germinal center (GC)-derived memory B cells (MBCs) are critical for humoral immunity as they differentiate into protective antibody-secreting cells during re-infection. GC formation and cellular interactions within the GC have been studied in detail, yet the exact signals that allow for the selection and exit of MBCs are not understood. Here, we show that IL-4 signaling in GC B cells directly downregulates BCL6 via negative autoregulation to release cells from the GC program and promote MBC formation. This selection event requires additional survival cues and can therefore result in either GC exit or death. We demonstrate that both increasing IL-4 bioavailability or limiting IL-4 signaling disrupt MBC selection stringency. In this way, IL-4 control of BCL6 expression serves as a tunable switch within the GC to tightly regulate MBC selection and affinity maturation.

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.

Neuronal binding by antibodies can be influenced by low pH stress during the isolation procedure

Low pH stress and its influence on antibody binding is a common consideration among chemists, but is only recently emerging as a consideration in Immunological studies. Antibody characterizations in Multiple Sclerosis (MS), an autoimmune disease of the Central Nervous System (CNS) has revealed that antibodies in the cerebrospinal fluid (CSF) of patients with Multiple Sclerosis bind to myelin-related and non-myelin antigen targets. Many laboratories have used molecular biology techniques to generate recombinant human antibodies (rhAbs) expressed by individual B cells from healthy donors and patients with systemic autoimmune disease to identify antigen targets. This approach has been adapted within the Neuroimmunology research community to investigate antigen targets of individual B cells in the CSF of MS patients. Our laboratory determines which antibodies to clone based on their immunogenetics and this method enriches for cloning of rhAbs that bind to neurons. However, newer technologies to assist in purification of these rhAbs from culture supernatants use an acidic elution buffer which may enhance low pH stress on the antibody structure. Our laboratory routinely uses a basic elution buffer to purify rhAbs from culture supernatants to avoid low pH stress to the antibody structure. Our goal was to investigate whether acidic elution of our rhAbs using Next Generation Chromatography would impact the rhAbs' ability to bind neurons. The limited data presented here for two neuron-binding rhAbs tested indicated that acidic elution buffers used during rhAb purification impacted the ability of rhAbs with low CDR3 charge to maintain binding to neuronal targets. Reproducibility in a larger panel of rhAbs and factors underlying these observations remain untested.

Assessing immunogenicity barriers of the HIV-1 envelope trimer

Understanding the balance between epitope shielding and accessibility on HIV-1 envelope (Env) trimers is essential to guide immunogen selection for broadly neutralizing antibody (bnAb) based vaccines. To investigate the antigenic space of Env immunogens, we created a strategy based on synthetic, high diversity, Designed Ankyrin Repeat Protein (DARPin) libraries. We show that DARPin Antigenicity Analysis (DANA), a purely in vitro screening tool, has the capability to extrapolate relevant information of antigenic properties of Env immunogens. DANA screens of stabilized, soluble Env trimers revealed that stronger trimer stabilization led to the selection of highly mutated DARPins with length variations and framework mutations mirroring observations made for bnAbs. By mimicking heterotypic prime-boost immunization regimens, DANA may be used to select immunogen combinations that favor the selection of trimer-reactive binders. This positions DANA as a versatile strategy for distilling fundamental antigenic features of immunogens, complementary to preclinical immunogenicity testing.

Recognition determinants of improved HIV-1 neutralization by a heavy chain matured pediatric antibody

The structural and characteristic features of HIV-1 broadly neutralizing antibodies (bnAbs) from chronically infected pediatric donors are currently unknown. Herein, we characterized a heavy chain matured HIV-1 bnAb 44m, identified from a pediatric elite-neutralizer. Interestingly, in comparison to its wild-type AIIMS-P01 bnAb, 44m exhibited moderately higher level of somatic hypermutations of 15.2%. The 44m neutralized 79% of HIV-1 heterologous viruses (n = 58) tested, with a geometric mean IC50 titer of 0.36 μg/mL. The cryo-EM structure of 44m Fab in complex with fully cleaved glycosylated native-like BG505.SOSIP.664.T332N gp140 envelope trimer at 4.4 Å resolution revealed that 44m targets the V3-glycan N332-supersite and GDIR motif to neutralize HIV-1 with improved potency and breadth, plausibly attributed by a matured heavy chain as compared to that of wild-type AIIMS-P01. This study further improves our understanding on pediatric HIV-1 bnAbs and structural basis of broad HIV-1 neutralization by 44m may be useful blueprint for vaccine design in future.

SARS-CoV-2 Omicron BA.1 breakthrough infection drives late remodeling of the memory B cell repertoire in vaccinated individuals

How infection by a viral variant showing antigenic drift impacts a preformed mature human memory B cell (MBC) repertoire remains an open question. Here, we studied the MBC response up to 6 months after SARS-CoV-2 Omicron BA.1 breakthrough infection in individuals previously vaccinated with three doses of the COVID-19 mRNA vaccine. Longitudinal analysis, using single-cell multi-omics and functional analysis of monoclonal antibodies from RBD-specific MBCs, revealed that a BA.1 breakthrough infection mostly recruited pre-existing cross-reactive MBCs with limited de novo response against BA.1-restricted epitopes. Reorganization of clonal hierarchy and new rounds of germinal center reactions, however, combined to maintain diversity and induce progressive maturation of the MBC repertoire against common Hu-1 and BA.1, but not BA.5-restricted, SARS-CoV-2 Spike RBD epitopes. Such remodeling was further associated with a marked improvement in overall neutralizing breadth and potency. These findings have fundamental implications for the design of future vaccination booster strategies.

Memory B cell development elicited by mRNA booster vaccinations in the elderly

Despite mRNA vaccination, elderly individuals remain especially vulnerable to severe consequences of SARS-CoV-2 infection. Here, we compare the memory B cell responses in a cohort of elderly and younger individuals who received mRNA booster vaccinations. Plasma neutralizing potency and breadth were similar between the two groups. By contrast, the absolute number of SARS-CoV-2-specific memory B cells was lower in the elderly. Antibody sequencing revealed that the SARS-CoV-2-specific elderly memory compartments were more clonal and less diverse. Notably, memory antibodies from the elderly preferentially targeted the ACE2-binding site on the RBD, while those from younger individuals targeted less accessible but more conserved epitopes. Nevertheless, individual memory antibodies elicited by booster vaccines in the elderly and younger individuals showed similar levels of neutralizing activity and breadth against SARS-CoV-2 variants. Thus, the relatively diminished protective effects of vaccination against serious disease in the elderly are associated with a smaller number of antigen-specific memory B cells that express altered antibody repertoires.

Polyclonal lymphoid expansion drives paraneoplastic autoimmunity in neuroblastoma

Neuroblastoma is a lethal childhood solid tumor of developing peripheral nerves. Two percent of children with neuroblastoma develop opsoclonus myoclonus ataxia syndrome (OMAS), a paraneoplastic disease characterized by cerebellar and brainstem-directed autoimmunity but typically with outstanding cancer-related outcomes. We compared tumor transcriptomes and tumor-infiltrating T and B cell repertoires from 38 OMAS subjects with neuroblastoma to 26 non-OMAS-associated neuroblastomas. We found greater B and T cell infiltration in OMAS-associated tumors compared to controls and showed that both were polyclonal expansions. Tertiary lymphoid structures (TLSs) were enriched in OMAS-associated tumors. We identified significant enrichment of the major histocompatibility complex (MHC) class II allele HLA-DOB∗01:01 in OMAS patients. OMAS severity scores were associated with the expression of several candidate autoimmune genes. We propose a model in which polyclonal auto-reactive B lymphocytes act as antigen-presenting cells and drive TLS formation, thereby supporting both sustained polyclonal T cell-mediated anti-tumor immunity and paraneoplastic OMAS neuropathology.

Double-negative-2 B cells are the major synovial plasma cell precursor in rheumatoid arthritis

B cells are key pathogenic drivers of chronic inflammation in rheumatoid arthritis (RA). There is limited understanding of the relationship between synovial B cell subsets and pathogenic antibody secreting cells (ASCs). This knowledge is crucial for the development of more targeted B-cell depleting therapies. While CD11c+ double-negative 2 (DN2) B cells have been suggested as an ASC precursor in lupus, to date there is no proven link between the two subsets in RA. We have used both single-cell gene expression and BCR sequencing to study synovial B cells from patients with established RA, in addition to flow cytometry of circulating B cells. To better understand the differentiation patterns within the diseased tissue, a combination of RNA-based trajectory inference and clonal lineage analysis of BCR relationships were used. Both forms of analysis indicated that DN2 B cells serve as a major precursors to synovial ASCs. This study advances our understanding of B cells in RA and reveals the origin of pathogenic ASCs in the RA synovium. Given the significant role of DN2 B cells as a progenitor to pathogenic B cells in RA, it is important to conduct additional research to investigate the origins of DN2 B cells in RA and explore their potential as therapeutic targets in place of the less specific pan-B cells depletion therapies currently in use.