sciCSR infers B cell state transition and predicts class-switch recombination dynamics using single-cell transcriptomic data.

Class-switch recombination (CSR) is an integral part of B cell maturation. Here we present sciCSR (pronounced 'scissor', single-cell inference of class-switch recombination), a computational pipeline that analyzes CSR events and dynamics of B cells from single-cell RNA sequencing (scRNA-seq) experiments. Validated on both simulated and real data, sciCSR re-analyzes scRNA-seq alignments to differentiate productive heavy-chain immunoglobulin transcripts from germline 'sterile' transcripts. From a snapshot of B cell scRNA-seq data, a Markov state model is built to infer the dynamics and direction of CSR. Applying sciCSR on severe acute respiratory syndrome coronavirus 2 vaccination time-course scRNA-seq data, we observe that sciCSR predicts, using data from an earlier time point in the collected time-course, the isotype distribution of B cell receptor repertoires of subsequent time points with high accuracy (cosine similarity ~0.9). Using processes specific to B cells, sciCSR identifies transitions that are often missed by conventional RNA velocity analyses and can reveal insights into the dynamics of B cell CSR during immune response.

BRepertoire: a user-friendly web server for analysing antibody repertoire data.

Antibody repertoire analysis by high throughput sequencing is now widely used, but a persisting challenge is enabling immunologists to explore their data to discover discriminating repertoire features for their own particular investigations. Computational methods are necessary for large-scale evaluation of antibody properties. We have developed BRepertoire, a suite of user-friendly web-based software tools for large-scale statistical analyses of repertoire data. The software is able to use data preprocessed by IMGT, and performs statistical and comparative analyses with versatile plotting options. BRepertoire has been designed to operate in various modes, for example analysing sequence-specific V(D)J gene usage, discerning physico-chemical properties of the CDR regions and clustering of clonotypes. Those analyses are performed on the fly by a number of R packages and are deployed by a shiny web platform. The user can download the analysed data in different table formats and save the generated plots as image files ready for publication. We believe BRepertoire to be a versatile analytical tool that complements experimental studies of immune repertoires. To illustrate the server's functionality, we show use cases including differential gene usage in a vaccination dataset and analysis of CDR3H properties in old and young individuals. The server is accessible under http://mabra.biomed.kcl.ac.uk/BRepertoire.

A Reassessment of IgM Memory Subsets in Humans.

From paired blood and spleen samples from three adult donors, we performed high-throughput VH sequencing of human B cell subsets defined by IgD and CD27 expression: IgD(+)CD27(+) ("marginal zone [MZ]"), IgD(-)CD27(+) ("memory," including IgM ["IgM-only"], IgG and IgA) and IgD(-)CD27(-) cells ("double-negative," including IgM, IgG, and IgA). A total of 91,294 unique sequences clustered in 42,670 clones, revealing major clonal expansions in each of these subsets. Among these clones, we further analyzed those shared sequences from different subsets or tissues for VH gene mutation, H-CDR3-length, and VH/JH usage, comparing these different characteristics with all sequences from their subset of origin for which these parameters constitute a distinct signature. The IgM-only repertoire profile differed notably from that of MZ B cells by a higher mutation frequency and lower VH4 and higher JH6 gene usage. Strikingly, IgM sequences from clones shared between the MZ and the memory IgG/IgA compartments showed a mutation and repertoire profile of IgM-only and not of MZ B cells. Similarly, all IgM clonal relationships (among MZ, IgM-only, and double-negative compartments) involved sequences with the characteristics of IgM-only B cells. Finally, clonal relationships between tissues suggested distinct recirculation characteristics between MZ and switched B cells. The "IgM-only" subset (including cells with its repertoire signature but higher IgD or lower CD27 expression levels) thus appear as the only subset showing precursor-product relationships with CD27(+) switched memory B cells, indicating that they represent germinal center-derived IgM memory B cells and that IgM memory and MZ B cells constitute two distinct entities.

Redemption of autoantibodies on anergic B cells by variable-region glycosylation and mutation away from self-reactivity.

The best-understood mechanisms for achieving antibody self/non-self discrimination discard self-reactive antibodies before they can be tested for binding microbial antigens, potentially creating holes in the repertoire. Here we provide evidence for a complementary mechanism: retaining autoantibodies in the repertoire displayed as low levels of IgM and high IgD on anergic B cells, masking a varying proportion of autoantibody-binding sites with carbohydrates, and removing their self-reactivity by somatic hypermutation and selection in germinal centers (GCs). Analysis of human antibody sequences by deep sequencing of isotype-switched memory B cells or in IgG antibodies elicited against allogeneic RhD+ erythrocytes, vaccinia virus, rotavirus, or tetanus toxoid provides evidence for reactivation of anergic IgM(low) IgD+ IGHV4-34+ B cells and removal of cold agglutinin self-reactivity by hypermutation, often accompanied by mutations that inactivated an N-linked glycosylation sequon in complementarity-determining region 2 (CDR2). In a Hy10 antibody transgenic model where anergic B cells respond to a biophysically defined lysozyme epitope displayed on both foreign and self-antigens, cell transfers revealed that anergic IgM(low) IgD+ B cells form twice as many GC progeny as naïve IgM(hi) IgD+ counterparts. Their GC progeny were rapidly selected for CDR2 mutations that blocked 72% of antigen-binding sites with N-linked glycan, decreased affinity 100-fold, and then cleared the binding sites of blocking glycan. These results provide evidence for a mechanism to acquire self/non-self discrimination by somatic mutation away from self-reactivity, and reveal how varying the efficiency of N-glycosylation provides a mechanism to modulate antibody avidity.

Immunoglobulin light chain allelic inclusion in systemic lupus erythematosus.

The principles of allelic exclusion state that each B cell expresses a single light and heavy chain pair. Here, we show that B cells with both kappa and lambda light chains (Igκ and Igλ) are enriched in some patients with the systemic autoimmune disease systemic lupus erythematosus (SLE), but not in the systemic autoimmune disease control granulomatosis with polyangiitis. Detection of dual Igκ and Igλ expression by flow cytometry could not be abolished by acid washing or by DNAse treatment to remove any bound polyclonal antibody or complexes, and was retained after two days in culture. Both surface and intracytoplasmic dual light chain expression was evident by flow cytometry and confocal microscopy. We observed reduced frequency of rearrangements of the kappa-deleting element (KDE) in SLE and an inverse correlation between the frequency of KDE rearrangement and the frequency of dual light chain expressing B cells. We propose that dual expression of Igκ and Igλ by a single B cell may occur in some patients with SLE when this may be a consequence of reduced activity of the KDE.

Granulocyte colony stimulating factor exacerbates antineutrophil cytoplasmic antibody vasculitis.

OBJECTIVES: Granulocyte colony stimulating factor (GCSF) is important in mobilising neutrophils from the bone marrow but also has a range of proinflammatory effects. We therefore decided to investigate the role of GCSF in antineutrophil cytoplasmic antibody (ANCA) vasculitis. METHODS: We measured GCSF levels in the serum of 38 patients with active ANCA vasculitis compared with 31 age-matched controls, and assessed the effect of GCSF priming on the response of human neutrophils to ANCA. We also examined the effect of exogenous GCSF administration in a murine model of antimyeloperoxidase (anti-MPO) vasculitis, and the effect of GCSF on murine neutrophil activation. RESULTS: The serum levels of GCSF in patients with active ANCA vasculitis were significantly higher than those of age matched healthy controls (mean 38.04 vs 18.35 pg/ml, p<0.001). Furthermore, we demonstrated that GCSF primed human neutrophils in vitro for a respiratory burst in response to anti-MPO ANCA. In an anti-MPO antibody transfer model, mice given GCSF had more crescents (mean 29.1% vs 5.8% per glomerular cross section, p<0.05), more macrophages (mean 3.2 vs 1.2 per glomerular cross-section, p<0.01), higher serum creatines (mean 13.6 vs 8.3 µmol/l, p<0.05) and more haematuria (p<0.05) compared with controls. In vivo administration of GCSF with lipopolysaccharide (LPS), but not LPS alone, led to upregulation of CD11c on murine neutrophils. CONCLUSIONS: These data suggest that GCSF, which is raised in patient serum, may play an important role in exacerbating disease in ANCA vasculitis. In addition, GCSF therapy for neutropenia should be used with caution in these patients.

PCR amplification and high throughput sequencing of immunoglobulin heavy chain genes from formalin-fixed paraffin-embedded human biopsies.

The use of high throughput sequencing (HTS) technologies in biomedicine is expanding in a variety of fields in recent years. The 454 system is an HTS platform that is ideally suited to characterize B cell receptor (BCR) repertoires by sequencing of immunoglobulin (Ig) genes, as it is able to sequence stretches of several hundred nucleotides. Most studies that used this platform for antibody repertoire analyses have started from fresh or frozen tissues or peripheral blood samples, and rely on starting with optimal quality DNA. In this paper we demonstrate that BCR repertoire analysis can be done using DNA from formalin-fixed paraffin-embedded (FFPE) human tissue samples. The heterogeneity of BCR repertoires we obtained confirms the plausibility of HTS of DNA from FFPE specimens. The establishment of experimental protocols and computational tools that enable sequence data analysis from the low quality DNA of FFPE tissues is important for enabling research, as it would enable the use of the rich source of preserved samples in clinical biobanks and biopsy archives.

B cell profiles, antibody repertoire and reactivity reveal dysregulated responses with autoimmune features in melanoma.

B cells are known to contribute to the anti-tumor immune response, especially in immunogenic tumors such as melanoma, yet humoral immunity has not been characterized in these cancers to detail. Here we show comprehensive phenotyping in samples of circulating and tumor-resident B cells as well as serum antibodies in melanoma patients. Memory B cells are enriched in tumors compared to blood in paired samples and feature distinct antibody repertoires, linked to specific isotypes. Tumor-associated B cells undergo clonal expansion, class switch recombination, somatic hypermutation and receptor revision. Compared with blood, tumor-associated B cells produce antibodies with proportionally higher levels of unproductive sequences and distinct complementarity determining region 3 properties. The observed features are signs of affinity maturation and polyreactivity and suggest an active and aberrant autoimmune-like reaction in the tumor microenvironment. Consistent with this, tumor-derived antibodies are polyreactive and characterized by autoantigen recognition. Serum antibodies show reactivity to antigens attributed to autoimmune diseases and cancer, and their levels are higher in patients with active disease compared to post-resection state. Our findings thus reveal B cell lineage dysregulation with distinct antibody repertoire and specificity, alongside clonally-expanded tumor-infiltrating B cells with autoimmune-like features, shaping the humoral immune response in melanoma.

Disparate roles of ATR and ATM in immunoglobulin class switch recombination and somatic hypermutation.

Class switch recombination (CSR) and somatic hypermutation (SHM) are mechanistically related processes initiated by activation-induced cytidine deaminase. Here, we have studied the role of ataxia telangiectasia and Rad3-related protein (ATR) in CSR by analyzing the recombinational junctions, resulting from in vivo switching, in cells from patients with mutations in the ATR gene. The proportion of cells that have switched to immunoglobulin (Ig)A and IgG in the peripheral blood seems to be normal in ATR-deficient (ATRD) patients and the recombined S regions show a normal "blunt end-joining," but impaired end joining with partially complementary (1-3 bp) DNA ends. There was also an increased usage of microhomology at the mu-alpha switch junctions, but only up to 9 bp, suggesting that the end-joining pathway requiring longer microhomologies (> or =10 bp) may be ATR dependent. The SHM pattern in the Ig variable heavy chain genes is altered, with fewer mutations occurring at A and more mutations at T residues and thus a loss of strand bias in targeting A/T pairs within certain hotspots. These data suggest that the role of ATR is partially overlapping with that of ataxia telangiectasia-mutated protein, but that the former is also endowed with unique functional properties in the repair processes during CSR and SHM.

IgA-producing plasma cells originate from germinal centers that are induced by B-cell receptor engagement in humans.

BACKGROUND & AIMS: IgA contributes to homeostatic balance between host and intestinal microbiota. Mechanisms that initiate the IgA response are unclear and likely to differ between humans and animal models. We used multiple experimental approaches to investigate the origin of human intestinal plasma cells that produce IgA in the gastrointestinal tract. METHODS: Complexity of IgA-producing plasma cell populations in human gastrointestinal mucosa and bone marrow and the specific response to oral cholera vaccine were compared by analysis of immunoglobulin genes. Flow cytometry, gene expression analysis, and immunohistochemistry were used to analyze signaling pathways induced by B-cell receptor engagement in human gut-associated lymphoid tissue (GALT) and involvement of innate immunity in B-cell activation in GALT compared with nonintestinal sites. RESULTS: Human intestinal IgA-producing plasma cells appeared to be of germinal center origin; there was no evidence for the population complexity that accompanies multiple pathways of derivation observed in bone marrow. In germinal center B cells of human GALT, Btk and Erk are phosphorylated, CD22 is down-regulated, Lyn is translocated to the cell membrane, and Fos and Jun are up-regulated; these features indicate B-cell receptor ligation during germinal center evolution. No differences in innate activation of B cells were observed in GALT, compared with peripheral immune compartments. CONCLUSIONS: IgA-producing plasma cells appear to be derived from GALT germinal centers in humans. B-cell receptor engagement promotes formation of germinal centers of GALT, with no more evidence for innate immune receptor activation in the mucosa than nonintestinal immune compartments. Germinal centers in GALT should be targets of mucosal vaccinations because they are the source of human intestinal IgA response.

High-throughput immunoglobulin repertoire analysis distinguishes between human IgM memory and switched memory B-cell populations.

B-cell receptor (BCR) diversity is achieved centrally by rearrangement of Variable, Diversity, and Joining genes, and peripherally by somatic hypermutation and class-switching of the rearranged genes. Peripheral B-cell populations are subject to both negative and positive selection events in the course of their development that have the potential to shape the BCR repertoire. The origin of IgM(+)IgD(+)CD27(+) (IgM memory) cells is controversial. It has been suggested that they may be a prediversified, antigen-independent, population of cells or that they are a population of cells that develop in response to T-independent antigens. Most recently, it was suggested that the majority of IgM memory cells are directly related to switched memory cells and are early emigrants from the germinal center reaction. Advances in sequencing technology have enabled us to undertake large scale IGH repertoire analysis of transitional, naive, IgM memory and switched memory B-cell populations. We find that the memory B-cell repertoires differ from the transitional and naive repertoires, and that the IgM memory repertoire is distinct from that of class-switched memory. Thus we conclude that a large proportion of IgM memory cells develop in response to different stimuli than for class-switched memory cell development.

B cells and aging: molecules and mechanisms.

Recent advances allow aging-associated changes in B-cell function to be approached at a mechanistic level. Reduced expression of genes crucial to lineage commitment and differentiation yield diminished B-cell production. Moreover, intrinsic differences in the repertoire generated by B-cell precursors in aged individuals, coupled with falling B-cell generation rates and life-long homeostatic competition, result in narrowed clonotypic diversity. Similarly, reductions in gene products crucial for immunoglobulin class switch recombination and somatic hypermutation impact the efficacy of humoral immune responses. Together, these findings set the stage for integrated analyses of how age-related changes at the molecular, cellular and population levels interact to yield the overall aging phenotype.

Pandemic, Epidemic, Endemic: B Cell Repertoire Analysis Reveals Unique Anti-Viral Responses to SARS-CoV-2, Ebola and Respiratory Syncytial Virus.

Immunoglobulin gene heterogeneity reflects the diversity and focus of the humoral immune response towards different infections, enabling inference of B cell development processes. Detailed compositional and lineage analysis of long read IGH repertoire sequencing, combining examples of pandemic, epidemic and endemic viral infections with control and vaccination samples, demonstrates general responses including increased use of IGHV4-39 in both Zaire Ebolavirus (EBOV) and COVID-19 patient cohorts. We also show unique characteristics absent in Respiratory Syncytial Virus or yellow fever vaccine samples: EBOV survivors show unprecedented high levels of class switching events while COVID-19 repertoires from acute disease appear underdeveloped. Despite the high levels of clonal expansion in COVID-19 IgG1 repertoires there is a striking lack of evidence of germinal centre mutation and selection. Given the differences in COVID-19 morbidity and mortality with age, it is also pertinent that we find significant differences in repertoire characteristics between young and old patients. Our data supports the hypothesis that a primary viral challenge can result in a strong but immature humoral response where failures in selection of the repertoire risk off-target effects.

B cell immunosenescence: different features of naive and memory B cells in elderly.

Elderly people show a reduced protection against new infections and a decreased response to vaccines as a consequence of impairment of both cellular and humoral immunity. In this paper we have studied memory/naïve B cells in the elderly, evaluating surface immunoglobulin expression, production of the pro- and anti-inflammatory cytokines, tumor necrosis factor (TNF)-α and interleukin (IL)-10, and presence of somatic hypermutation, focusing on the IgG(+)IgD(-)CD27(-) double negative (DN) B cells that are expanded in the elderly. Our results show that naïve B cells from young donors need a sufficiently strong stimulus to be activated "in vitro", while naïve B cells from old subjects are able to produce IL-10 and TNF-α when stimulated "physiologically" (α-CD40/IL-4), suggesting that these cells might play a role in the control of the immuno-inflammatory environment in the elderly. In addition, in the elderly there is an accumulation of DN B cells with a reduced rate of somatic hypermutation. Thus, DN B lymphocytes may be exhausted cells that are expanded and accumulate as a by-product of persistent stimulation or impaired germinal center formation.

Single-Cell Transcriptomic Analyses Define Distinct Peripheral B Cell Subsets and Discrete Development Pathways.

Separation of B cells into different subsets has been useful to understand their different functions in various immune scenarios. In some instances, the subsets defined by phenotypic FACS separation are relatively homogeneous and so establishing the functions associated with them is straightforward. Other subsets, such as the "Double negative" (DN, CD19+CD27-IgD-) population, are more complex with reports of differing functionality which could indicate a heterogeneous population. Recent advances in single-cell techniques enable an alternative route to characterize cells based on their transcriptome. To maximize immunological insight, we need to match prior data from phenotype-based studies with the finer granularity of the single-cell transcriptomic signatures. We also need to be able to define meaningful B cell subsets from single cell analyses performed on PBMCs, where the relative paucity of a B cell signature means that defining B cell subsets within the whole is challenging. Here we provide a reference single-cell dataset based on phenotypically sorted B cells and an unbiased procedure to better classify functional B cell subsets in the peripheral blood, particularly useful in establishing a baseline cellular landscape and in extracting significant changes with respect to this baseline from single-cell datasets. We find 10 different clusters of B cells and applied a novel, geometry-inspired, method to RNA velocity estimates in order to evaluate the dynamic transitions between B cell clusters. This indicated the presence of two main developmental branches of memory B cells. A T-independent branch that involves IgM memory cells and two DN subpopulations, culminating in a population thought to be associated with Age related B cells and the extrafollicular response. The other, T-dependent, branch involves a third DN cluster which appears to be a precursor of classical memory cells. In addition, we identify a novel DN4 population, which is IgE rich and closely linked to the classical/precursor memory branch suggesting an IgE specific T-dependent cell population.

Tumor-Infiltrating B Lymphocyte Profiling Identifies IgG-Biased, Clonally Expanded Prognostic Phenotypes in Triple-Negative Breast Cancer.

In breast cancer, humoral immune responses may contribute to clinical outcomes, especially in more immunogenic subtypes. Here, we investigated B lymphocyte subsets, immunoglobulin expression, and clonal features in breast tumors, focusing on aggressive triple-negative breast cancers (TNBC). In samples from patients with TNBC and healthy volunteers, circulating and tumor-infiltrating B lymphocytes (TIL-B) were evaluated. CD20+CD27+IgD- isotype-switched B lymphocytes were increased in tumors, compared with matched blood. TIL-B frequently formed stromal clusters with T lymphocytes and engaged in bidirectional functional cross-talk, consistent with gene signatures associated with lymphoid assembly, costimulation, cytokine-cytokine receptor interactions, cytotoxic T-cell activation, and T-cell-dependent B-cell activation. TIL-B-upregulated B-cell receptor (BCR) pathway molecules FOS and JUN, germinal center chemokine regulator RGS1, activation marker CD69, and TNFα signal transduction via NFκB, suggesting BCR-immune complex formation. Expression of genes associated with B lymphocyte recruitment and lymphoid assembly, including CXCL13, CXCR4, and DC-LAMP, was elevated in TNBC compared with other subtypes and normal breast. TIL-B-rich tumors showed expansion of IgG but not IgA isotypes, and IgG isotype switching positively associated with survival outcomes in TNBC. Clonal expansion was biased toward IgG, showing expansive clonal families with specific variable region gene combinations and narrow repertoires. Stronger positive selection pressure was present in the complementarity determining regions of IgG compared with their clonally related IgA in tumor samples. Overall, class-switched B lymphocyte lineage traits were conspicuous in TNBC, associated with improved clinical outcomes, and conferred IgG-biased, clonally expanded, and likely antigen-driven humoral responses. SIGNIFICANCE: Tumor-infiltrating B lymphocytes assemble in clusters, undergoing B-cell receptor-driven activation, proliferation, and isotype switching. Clonally expanded, IgG isotype-biased humoral immunity associates with favorable prognosis primarily in triple-negative breast cancers.

Lambda light chain revision in the human intestinal IgA response.

Revision of Ab L chains by secondary rearrangement in mature B cells has the potential to change the specific target of the immune response. In this study, we show for the first time that L chain revision is normal and widespread in the largest Ab producing population in man: intestinal IgA plasma cells (PC). Biases in the productive and non-productive repertoire of lambda L chains, identification of the circular products of rearrangement that have the characteristic biases of revision, and identification of RAG genes and protein all reflect revision during normal intestinal IgA PC development. We saw no evidence of IgH revision, probably due to inappropriately orientated recombination signal sequences, and little evidence of kappa-chain revision, probably due to locus inactivation by the kappa-deleting element. We propose that the lambda L chain locus is available and a principal modifier and diversifier of Ab specificity in intestinal IgA PCs.

Age-Related Changes in B Cells Relevant to Vaccine Responses.

Older people have reduced immune responses to infection and vaccination. B cell activation is key for the efficacy of the vaccine response, but there are several age-related changes in B cells which may contribute to the loss of vaccine efficacy. Different subpopulations of B cells have different functions and phenotypes. These populations can change as we age; older people have been shown to have fewer "IgM memory" cells, regulatory B cells and plasma cells and more IgD-CD27- "double-negative" and "age-related B cells." While the overall quantity of antibody in the blood does not change, the quality of the B cell response changes; producing less specific antibodies upon challenge and more autoreactive antibodies. This could be due to changes in selection pressures, as has been demonstrated by repertoire sequencing of different subsets of B cells at different ages. Other changes in antibody repertoire are seen, including greater levels of IgG2 in older people and altered IgG1 IGHV gene usage. Since B cells rely on their environment for efficient responses, some of these changes may be due to age-related changes in accessory cells/signals. Other changes appear to be intrinsic to older/aged B cells themselves, such as their tendency to produce greater levels of inflammatory cytokines.