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.

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.

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.

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.

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.

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.

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.

Transitional B Cells in Early Human B Cell Development - Time to Revisit the Paradigm?

The B cell repertoire is generated in the adult bone marrow by an ordered series of gene rearrangement processes that result in massive diversity of immunoglobulin (Ig) genes and consequently an equally large number of potential specificities for antigen. As the process is essentially random, the cells exhibiting excess reactivity with self-antigens are generated and need to be removed from the repertoire before the cells are fully mature. Some of the cells are deleted, and some will undergo receptor editing to see if changing the light chain can rescue an autoreactive antibody. As a consequence, the binding properties of the B cell receptor are changed as development progresses through pre-B ≫ immature ≫ transitional ≫ naïve phenotypes. Using long-read, high-throughput, sequencing we have produced a unique set of sequences from these four cell types in human bone marrow and matched peripheral blood, and our results describe the effects of tolerance selection on the B cell repertoire at the Ig gene level. Most strong effects of selection are seen within the heavy chain repertoire and can be seen both in gene usage and in CDRH3 characteristics. Age-related changes are small, and only the size of the CDRH3 shows constant and significant change in these data. The paucity of significant changes in either kappa or lambda light chain repertoires implies that either the heavy chain has more influence over autoreactivity than light chain and/or that switching between kappa and lambda light chains, as opposed to switching within the light chain loci, may effect a more successful autoreactive rescue by receptor editing. Our results show that the transitional cell population contains cells other than those that are part of the pre-B ≫ immature ≫ transitional ≫ naïve development pathway, since the population often shows a repertoire that is outside the trajectory of gene loss/gain between pre-B and naïve stages.

Significant Differences in Physicochemical Properties of Human Immunoglobulin Kappa and Lambda CDR3 Regions.

Antibody variable regions are composed of a heavy and a light chain, and in humans, there are two light chain isotypes: kappa and lambda. Despite their importance in receptor editing, the light chain is often overlooked in the antibody literature, with the focus being on the heavy chain complementarity-determining region (CDR)-H3 region. In this paper, we set out to investigate the physicochemical and structural differences between human kappa and lambda light chain CDR regions. We constructed a dataset containing over 29,000 light chain variable region sequences from IgM-transcribing, newly formed B cells isolated from human bone marrow and peripheral blood. We also used a published human naïve dataset to investigate the CDR-H3 properties of heavy chains paired with kappa and lambda light chains and probed the Protein Data Bank to investigate the structural differences between kappa and lambda antibody CDR regions. We found that kappa and lambda light chains have very different CDR physicochemical and structural properties, whereas the heavy chains with which they are paired do not differ significantly. We also observed that the mean CDR3 N nucleotide addition in the kappa, lambda, and heavy chain gene rearrangements are correlated within donors but can differ between donors. This indicates that terminal deoxynucleotidyl transferase may work with differing efficiencies between different people but the same efficiency in the different classes of immunoglobulin chain within one person. We have observed large differences in the physicochemical and structural properties of kappa and lambda light chain CDR regions. This may reflect different roles in the humoral immune response.

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.

Age-related aspects of human IgM(+) B cell heterogeneity.

The CD27(+) IgD(+) B cell population, known as IgM memory, reduces with age. It is thought that this population is responsible for pneumococcal polysaccharide T-independent responses, and that the age-related reduction might be partially responsible for the increased susceptibility of older people to bacterial pathogens. There are other IgM(+) B cell populations that do not express IgD. We compared the different IgM populations using high-throughput sequencing of the immunoglobulin (Ig) gene repertoire and multidimensional cell phenotyping and found that the different populations of IgM cells, defined by CD27 and IgD expression, have repertoire differences. Some of these differences are likely indicative of different selection pressures in an immune response, although the older individuals were found to have a changed repertoire in naive B cells, which may contribute to some of the changes seen in memory cells. In addition, even within the CD27(+) IgD(+) IgM memory population there are multiple cell types. We show that the level of IgM expression varies substantially and hypothesize that this distinguishes between T-dependent and T-independent types of IgM memory cells. Significant age-related changes in the relative proportions of these populations may exacerbate the reduction in T-independent responders in old age.

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.

Monoclonal paraprotein influences baseline B-cell repertoire diversity and perturbates influenza vaccination-induced B-cell response.

Monoclonal gammopathy of undetermined significance (MGUS) arises from a clonal expansion of plasma cells in the bone marrow, secreting monoclonal (M) paraprotein. It is associated with increased susceptibility to infections, which may reflect altered B-cell repertoire. To investigate this, we examined the immunoglobulin (Ig) M, IgG, and IgA B-cell repertoire diversity in MGUS at baseline and after influenza vaccination (n = 16) in comparison with healthy controls (HCs; n = 16). The Complementary Determining Region 3 region of the immunoglobulin heavy chain variable region gene was amplified and B-cell spectratypes analyzed by high-resolution electrophoresis. Spectratype Gaussian distribution, kurtosis, and skewness were quantified to measure repertoire shifts. Both HC and MGUS baseline spectratypes show interindividual variability that is more pronounced in the IGHG and IGHA repertoires. Overall, baseline B-cell repertoire is more altered in MGUS, with oligoclonality observed in 50% (p = 0.01). Postvaccination, significant differences emerged in MGUS in relation to M-protein levels. High M-protein concentration is associated with a more oligoclonal IgG and IgA response at day 7 postvaccination, and, in contrast to HCs, vaccination also induced significant perturbations in the MGUS IgM repertoire at day 7 (p = 0.005). Monoclonal expansion in MGUS thus has an effect on the baseline B-cell repertoire and influences the recruited repertoire upon vaccination.

Immunoglobulin kappa variable region gene selection during early human B cell development in health and systemic lupus erythematosus.

The unique specificity of the B cell receptor is generated by an ordered sequence of gene rearrangement events. Once IGH genes have rearranged, rearrangement at the IGK locus is initiated followed by the IGL locus if functional IGK rearrangement is not achieved. Receptor specificity can subsequently be altered by secondary light chain editing based on the features of the heavy and light chain combination. The final profile of expressed genes is not random and biases in this profile are associated with several autoimmune diseases. However, how and when biases are created is not known. To increase our understanding of the processes of selection and editing of IGK rearrangements, we compared four groups of rearrangements of IGK acquired by next generation sequencing. First, expressed rearrangements of IGK from cDNA of IGK expressing B cells. Second, productive rearrangements of IGK from DNA of the same kappa expressing B cells. Third, non-productive rearrangements of IGK from DNA of IGK and IGL expressing B cells, and fourth productively rearranged IGK from DNA of IGL expressing B cells. The latter group would have been rejected during B cell development in favour of rearrangement at the IGL locus and are therefore selected against. We saw evidence that rearranged IGK segments can be selected at a checkpoint where the decision to rearrange the IGL locus is made. In addition, our data suggest that mechanisms regulating the expression or not of IGK rearrangements may also contribute to repertoire development and also that this latter component of the selection process is defective in SLE.

Immunoglobulin gene repertoire diversification and selection in the stomach - from gastritis to gastric lymphomas.

Chronic gastritis is characterized by gastric mucosal inflammation due to autoimmune responses or infection, frequently with Helicobacter pylori. Gastritis with H. pylori background can cause gastric mucosa-associated lymphoid tissue lymphoma (MALT-L), which sometimes further transforms into diffuse large B-cell lymphoma (DLBCL). However, gastric DLBCL can also be initiated de novo. The mechanisms underlying transformation into DLBCL are not completely understood. We analyzed immunoglobulin repertoires and clonal trees to investigate whether and how immunoglobulin gene repertoires, clonal diversification, and selection in gastritis, gastric MALT-L, and DLBCL differ from each other and from normal responses. The two gastritis types (positive or negative for H. pylori) had similarly diverse repertoires. MALT-L dominant clones (defined as the largest clones in each sample) presented higher diversification and longer mutational histories compared with all other conditions. DLBCL dominant clones displayed lower clonal diversification, suggesting the transforming events are triggered by similar responses in different patients. These results are surprising, as we expected to find similarities between the dominant clones of gastritis and MALT-L and between those of MALT-L and DLBCL.

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.

Human lymphocyte repertoires in ageing.

Deterioration of adaptive immunity with ageing may reflect changes in the repertoire of T cells and B cells available to respond to antigenic challenges, due to altered proportions and absolute numbers of lymphocyte subpopulations as well as changes in the repertoire of antigen receptor genes expressed by these cells. High-throughput DNA sequencing (HTS) now facilitates examination of immunoglobulin and T cell receptor gene rearrangements, and initial studies using these methods to study immune system ageing in humans have demonstrated age-related alterations in the receptor populations within lymphocyte subsets, as well as in repertoires responding to vaccination. Accurate measurement of repertoire diversity remains an experimental challenge. Studies of larger numbers of human subjects, analysis of defined lymphocyte subpopulations including antigen-specific populations, and controlling for factors such as chronic viral infections, will be important for gaining additional understanding of the impact of ageing on human lymphocyte populations.