IgTree: creating Immunoglobulin variable region gene lineage trees

Systematic evaluation of intratumoral and peripheral BCR repertoires in three cancers

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

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

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

New generalized metric based on branch length distance to compare B cell lineage trees

The B cell lineage tree encapsulates the successive phases of B cell differentiation and maturation, transitioning from hematopoietic stem cells to mature, antibody-secreting cells within the immune system. Mathematically, this lineage can be conceptualized as an evolutionary tree, where each node represents a distinct stage in B cell development, and the edges reflect the differentiation pathways. To compare these lineage trees, a rigorous mathematical metric is essential. Analyzing B cell lineage trees mathematically and quantifying changes in lineage attributes over time necessitates a comparison methodology capable of accurately assessing and measuring these changes. Addressing the intricacies of multiple B cell lineage tree comparisons, this study introduces a novel metric that enhances the precision of comparative analysis. This metric is formulated on principles of metric theory and evolutionary biology, quantifying the dissimilarities between lineage trees by measuring branch length distance and weight. By providing a framework for systematically classifying lineage trees, this metric facilitates the development of predictive models that are crucial for the creation of targeted immunotherapy and vaccines. To validate the effectiveness of this new metric, synthetic datasets that mimic the complexity and variability of real B cell lineage structures are employed. We demonstrated the ability of the new metric method to accurately capture the evolutionary nuances of B cell lineages.

Protocol for fast clonal family inference and analysis from large-scale B cell receptor repertoire sequencing data

The expeditious identification and comprehensive analysis of clonal families from extensive B cell receptor (BCR) repertoire sequencing data are imperative for elucidating the intricacies of B cell immune responses. Here, we introduce a computational pipeline designed to swiftly deduce clonal families from bulk BCR heavy-chain sequencing data, accompanied by a suite of functional modules tailored to streamline post-clustering analysis. The outlined methodology encompasses guidelines for software installation, meticulous data preparation, and the systematic inference and analysis of clonal families. For complete details on the use and execution of this protocol, please refer to Wang et al.1.

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

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

Anti-Citrullinated Protein Antibodies With Multiple Specificities Ameliorate Collagen Antibody-Induced Arthritis in a Time-Dependent Manner

OBJECTIVE

Anti-citrullinated protein antibodies (ACPAs) are highly specific for rheumatoid arthritis (RA) and have long been regarded as pathogenic. Despite substantial in vitro evidence supporting this claim, reports investigating the proinflammatory effects of ACPAs in animal models of arthritis are rare and include mixed results. Here, we sequenced the plasmablast antibody repertoire of a patient with RA and functionally characterized the encoded ACPAs.

METHODS

We expressed ACPAs from the antibody repertoire of a patient with RA and characterized their autoantigen specificities on antigen arrays and enzyme-linked immunosorbent assays. Binding affinities were estimated by bio-layer interferometry. Select ACPAs (n = 9) were tested in the collagen antibody-induced arthritis (CAIA) mouse model to evaluate their effects on joint inflammation.

RESULTS

Recombinant ACPAs bound preferentially and with high affinity (nanomolar range) to citrullinated (cit) autoantigens (primarily histones and fibrinogen) and to auto-cit peptidylarginine deiminase 4 (PAD4). ACPAs were grouped for in vivo testing based on their predominant cit-antigen specificities. Unexpectedly, injections of recombinant ACPAs significantly reduced paw thickness and arthritis severity in CAIA mice as compared with isotype-matched control antibodies (P ≤ 0.001). Bone erosion, synovitis, and cartilage damage were also significantly reduced (P ≤ 0.01). This amelioration of CAIA was observed for all the ACPAs tested and was independent of cit-PAD4 and cit-fibrinogen specificities. Furthermore, disease amelioration was more prominent when ACPAs were injected at earlier stages of CAIA than at later phases of the model.

CONCLUSION

Recombinant patient-derived ACPAs ameliorated CAIA. Their antiinflammatory effects were more preventive than therapeutic. This study highlights a potential protective role for ACPAs in arthritis.

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.

B cell phylogenetics in the single cell era

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

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

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

Abalign: a comprehensive multiple sequence alignment platform for B-cell receptor immune repertoires

The utilization of high-throughput sequencing (HTS) for B-cell receptor (BCR) immune repertoire analysis has become widespread in the fields of adaptive immunity and antibody drug development. However, the sheer volume of sequences generated by these experiments presents a challenge in data processing. Specifically, multiple sequence alignment (MSA), a critical aspect of BCR analysis, remains inadequate for handling massive BCR sequencing data and lacks the ability to provide immunoglobulin-specific information. To address this gap, we introduce Abalign, a standalone program specifically designed for ultrafast MSA of BCR/antibody sequences. Benchmark tests demonstrate that Abalign achieves comparable or even better accuracy than state-of-the-art MSA tools, and shows remarkable advantages in terms of speed and memory consumption, reducing the time required for high-throughput analysis from weeks to hours. In addition to its alignment capabilities, Abalign offers a broad range of BCR analysis features, including extracting BCRs, constructing lineage trees, assigning VJ genes, analyzing clonotypes, profiling mutations, and comparing BCR immune repertoires. With its user-friendly graphic interface, Abalign can be easily run on personal computers instead of computing clusters. Overall, Abalign is an easy-to-use and effective tool that enables researchers to analyze massive BCR/antibody sequences, leading to new discoveries in the field of immunoinformatics. The software is freely available at http://cao.labshare.cn/abalign/.

B cell M-CLL clones retain selection against replacement mutations in their immunoglobulin gene framework regions

Introduction

Chronic lymphocytic leukemia (CLL) is the most common adult leukemia, accounting for 30–40% of all adult leukemias. The dynamics of B-lymphocyte CLL clones with mutated immunoglobulin heavy chain variable region (IgHV) genes in their tumor (M-CLL) can be studied using mutational lineage trees.

Methods

Here, we used lineage tree-based analyses of somatic hypermutation (SHM) and selection in M-CLL clones, comparing the dominant (presumably malignant) clones of 15 CLL patients to their non-dominant (presumably normal) B cell clones, and to those of healthy control repertoires. This type of analysis, which was never previously published in CLL, yielded the following novel insights.

Results

CLL dominant clones undergo – or retain – more replacement mutations that alter amino acid properties such as charge or hydropathy. Although, as expected, CLL dominant clones undergo weaker selection for replacement mutations in the complementarity determining regions (CDRs) and against replacement mutations in the framework regions (FWRs) than non-dominant clones in the same patients or normal B cell clones in healthy controls, they surprisingly retain some of the latter selection in their FWRs. Finally, using machine learning, we show that even the non-dominant clones in CLL patients differ from healthy control clones in various features, most notably their expression of higher fractions of transition mutations.

Discussion

Overall, CLL seems to be characterized by significant loosening – but not a complete loss – of the selection forces operating on B cell clones, and possibly also by changes in SHM mechanisms.

Reconstructing B cell lineage trees with minimum spanning tree and genotype abundances

B cell receptor (BCR) genes exposed to an antigen undergo somatic hypermutations and Darwinian antigen selection, generating a large BCR-antibody diversity. This process, known as B cell affinity maturation, increases antibody affinity, forming a specific B cell lineage that includes the unmutated ancestor and mutated variants. In a B cell lineage, cells with a higher antigen affinity will undergo clonal expansion, while those with a lower affinity will not proliferate and probably be eliminated. Therefore, cellular (genotype) abundance provides a valuable perspective on the ongoing evolutionary process. Phylogenetic tree inference is often used to reconstruct B cell lineage trees and represents the evolutionary dynamic of BCR affinity maturation. However, such methods should process B-cell population data derived from experimental sampling that might contain different cellular abundances. There are a few phylogenetic methods for tracing the evolutionary events occurring in B cell lineages; best-performing solutions are time-demanding and restricted to analysing a reduced number of sequences, while time-efficient methods do not consider cellular abundances. We propose ClonalTree, a low-complexity and accurate approach to construct B-cell lineage trees that incorporates genotype abundances into minimum spanning tree (MST) algorithms. Using both simulated and experimental data, we demonstrate that ClonalTree outperforms MST-based algorithms and achieves a comparable performance to a method that explores tree-generating space exhaustively. Furthermore, ClonalTree has a lower running time, being more convenient for building B-cell lineage trees from high-throughput BCR sequencing data, mainly in biomedical applications, where a lower computational time is appreciable. It is hundreds to thousands of times faster than exhaustive approaches, enabling the analysis of a large set of sequences within minutes or seconds and without loss of accuracy. The source code is freely available at github.com/julibinho/ClonalTree.

Lineage Reconstruction of In Vitro Identified Antigen-Specific Autoreactive B Cells from Adaptive Immune Receptor Repertoires

The emergence, survival, growth and maintenance of autoreactive (AR) B-cell clones, the hallmark of humoral autoimmunity, leave their footprints in B-cell receptor repertoires. Collecting IgH sequences related to polyreactive (PR) ones from adaptive immune receptor repertoire (AIRR) datasets make the reconstruction and analysis of PR/AR B-cell lineages possible. We developed a computational approach, named ImmChainTracer, to extract members and to visualize clonal relationships of such B-cell lineages. Our approach was successfully applied on the IgH repertoires of patients suffering from monogenic hypomorphic RAG1 and 2 deficiency (pRD) or polygenic systemic lupus erythematosus (SLE) autoimmune diseases to identify relatives of AR IgH sequences and to track their fate in AIRRs. Signs of clonal expansion, affinity maturation and class-switching events in PR/AR and non-PR/AR B-cell lineages were revealed. An extension of our method towards B-cell expansion caused by any trigger (e.g., infection, vaccination or antibody development) may provide deeper insight into antigen specific B-lymphogenesis.

A scalable model for simulating multi-round antibody evolution and benchmarking of clonal tree reconstruction methods

Affinity maturation (AM) of B cells through somatic hypermutations (SHMs) enables the immune system to evolve to recognize diverse pathogens. The accumulation of SHMs leads to the formation of clonal lineages of antibody-secreting b cells that have evolved from a common naïve B cell. Advances in high-throughput sequencing have enabled deep scans of B cell receptor repertoires, paving the way for reconstructing clonal trees. However, it is not clear if clonal trees, which capture microevolutionary time scales, can be reconstructed using traditional phylogenetic reconstruction methods with adequate accuracy. In fact, several clonal tree reconstruction methods have been developed to fix supposed shortcomings of phylogenetic methods. Nevertheless, no consensus has been reached regarding the relative accuracy of these methods, partially because evaluation is challenging. Benchmarking the performance of existing methods and developing better methods would both benefit from realistic models of clonal lineage evolution specifically designed for emulating B cell evolution. In this paper, we propose a model for modeling B cell clonal lineage evolution and use this model to benchmark several existing clonal tree reconstruction methods. Our model, designed to be extensible, has several features: by evolving the clonal tree and sequences simultaneously, it allows modeling selective pressure due to changes in affinity binding; it enables scalable simulations of large numbers of cells; it enables several rounds of infection by an evolving pathogen; and, it models building of memory. In addition, we also suggest a set of metrics for comparing clonal trees and measuring their properties. Our results show that while maximum likelihood phylogenetic reconstruction methods can fail to capture key features of clonal tree expansion if applied naively, a simple post-processing of their results, where short branches are contracted, leads to inferences that are better than alternative methods.

Mutational patterns along different evolution paths of follicular lymphoma

Follicular lymphoma (FL) is an indolent disease, characterized by a median life expectancy of 18-20 years and by intermittent periods of relapse and remission. FL frequently transforms into the more aggressive diffuse large B cell lymphoma (t-FL). In previous studies, the analysis of immunoglobulin heavy chain variable region (IgHV) genes in sequential biopsies from the same patient revealed two different patterns of tumor clonal evolution: direct evolution, through acquisition of additional IgHV mutations over time, or divergent evolution, in which lymphoma clones from serial biopsies independently develop from a less-mutated common progenitor cell (CPC). Our goal in this study was to characterize the somatic hypermutation (SHM) patterns of IgHV genes in sequential FL samples from the same patients, and address the question of whether the mutation mechanisms (SHM targeting, DNA repair or both), or selection forces acting on the tumor clones, were different in FL samples compared to healthy control samples, or in late relapsed/transformed FL samples compared to earlier ones. Our analysis revealed differences in the distribution of mutations from each of the nucleotides when tumor and non-tumor clones were compared, while FL and transformed FL (t-FL) tumor clones displayed similar mutation distributions. Lineage tree measurements suggested that either initial clone affinity or selection thresholds were lower in FL samples compared to controls, but similar between FL and t-FL samples. Finally, we observed that both FL and t-FL tumor clones tend to accumulate larger numbers of potential N-glycosylation sites due to the introduction of new SHM. Taken together, these results suggest that transformation into t-FL, in contrast to initial FL development, is not associated with any major changes in DNA targeting or repair, or the selection threshold of the tumor clone.

IgTreeZ, A Toolkit for Immunoglobulin Gene Lineage Tree-Based Analysis, Reveals CDR3s Are Crucial for Selection Analysis

Somatic hypermutation (SHM) is an important diversification mechanism that plays a part in the creation of immune memory. Immunoglobulin (Ig) variable region gene lineage trees were used over the last four decades to model SHM and the selection mechanisms operating on B cell clones. We hereby present IgTreeZ (Immunoglobulin Tree analyZer), a python-based tool that analyses many aspects of Ig gene lineage trees and their repertoires. Using simulations, we show that IgTreeZ can be reliably used for mutation and selection analyses. We used IgTreeZ on empirical data, found evidence for different mutation patterns in different B cell subpopulations, and gained insights into antigen-driven selection in corona virus disease 19 (COVID-19) patients. Most importantly, we show that including the CDR3 regions in selection analyses - which is only possible if these analyses are lineage tree-based - is crucial for obtaining correct results. Overall, we present a comprehensive lineage tree analysis tool that can reveal new biological insights into B cell repertoire dynamics.

Gastric DLBCL clonal evolution as function of patient age

Diffuse large B cell lymphoma (DLBCL) is the most common type of NHL, accounting for about 40% of NHL cases, and is one of the most aggressive lymphomas. DLBCL is widespread in individuals aged more than 50 years old, with a maximum incidence in the seventh decade, but it may also occur in younger patients. DLBCL may occur in any immune system tissue, including those around the gastrointestinal tract, and even in the stomach, though gastric DLBCL has yet to be sufficiently investigated. This study aimed to understand changes in gastric Diffuse Large B cell lymphoma (gastric DLBCL) development with age. Immunoglobulin (Ig) heavy chain variable region genes were amplified from sections of nine preserved biopsies, from patients whose age varied between 25 and 89 years, sequenced and analyzed. We show first that identification of the malignant clone based on the biopsies is much less certain than was previously assumed; and second that, contrary to expectations, the repertoire of gastric B cell clones is more diverse among the elderly DLBCL patients than among the young.

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

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

Next Generation Sequencing of Cerebrospinal Fluid B Cell Repertoires in Multiple Sclerosis and Other Neuro-Inflammatory Diseases-A Comprehensive Review

During the last few decades, the role of B cells has been well established and redefined in neuro-inflammatory diseases, including multiple sclerosis and autoantibody-associated diseases. In particular, B cell maturation and trafficking across the blood–brain barrier (BBB) has recently been deciphered with the development of next-generation sequencing (NGS) approaches, which allow the assessment of representative cerebrospinal fluid (CSF) and peripheral blood B cell repertoires. In this review, we perform literature research focusing on NGS studies that allow further insights into B cell pathophysiology during neuro-inflammation. Besides the analysis of CSF B cells, the paralleled assessment of peripheral blood B cell repertoire provides deep insights into not only the CSF compartment, but also in B cell trafficking patterns across the BBB. In multiple sclerosis, CSF-specific B cell maturation, in combination with a bidirectional exchange of B cells across the BBB, is consistently detectable. These data suggest that B cells most likely encounter antigen(s) within the CSF and migrate across the BBB, with further maturation also taking place in the periphery. Autoantibody-mediated diseases, such as neuromyelitis optica spectrum disorder and LGI1 / NMDAR encephalitis, also show features of a CSF-specific B cell maturation and clonal connectivity with peripheral blood. In conclusion, these data suggest an intense exchange of B cells across the BBB, possibly feeding autoimmune circuits. Further developments in sequencing technologies will help to dissect the exact pathophysiologic mechanisms of B cells during neuro-inflammation.

Systematic memory B cell archiving and random display shape the human splenic marginal zone throughout life

Human memory B cells (MBCs) are generated and diversified in secondary lymphoid tissues throughout the organism. A paired immunoglobulin (Ig)-gene repertoire analysis of peripheral blood (PB) and splenic MBCs from infant, adult, and elderly humans revealed that throughout life, circulating MBCs are comprehensively archived in the spleen. Archive MBC clones are systematically preserved and uncoupled from class-switching. Clonality in the spleen increases steadily, but boosts at midlife, thereby outcompeting small clones. The splenic marginal zone (sMZ) represents a primed MBC compartment, generated from a stochastic exchange within the archive memory pool. This is supported by functional assays, showing that PB and splenic CD21+ MBCs acquire transient CD21high expression upon NOTCH2-stimulation. Our study provides insight that the human MBC system in PB and spleen is composed of three interwoven compartments: the dynamic relationship of circulating, archive, and its subset of primed (sMZ) memory changes with age, thereby contributing to immune aging.

Differential Effects of Fingolimod and Natalizumab on B Cell Repertoires in Multiple Sclerosis Patients

Natalizumab and fingolimod are effective multiple sclerosis (MS) therapies that disrupt lymphocyte migration but have differential effects on B cell maturation and trafficking. We investigated their effects on peripheral blood (PB) and cerebrospinal fluid (CSF) B cell repertoires using next-generation deep sequencing. Paired CSF and PB B cell subsets (naïve, CD27+ memory, and CD27-IgD- double-negative B cells and plasmablasts) were collected by applying flow cytometry at baseline and after 6 months of treatment and their respective heavy-chain variable region repertoires assessed by Illumina MiSeq. Treatment with fingolimod contracted, whereas natalizumab expanded circulating PB B cells. CSF B cell numbers remained stable following fingolimod treatment but decreased with natalizumab therapy. Clonal overlap between CSF and PB B cells was reduced with natalizumab treatment but remained stable with fingolimod therapy. Lineage analyses of pre- and posttreatment CSF B cell repertoires revealed large, clonally expanded B cell clusters in natalizumab-treated MS patients but no intrathecal clonal expansion following fingolimod therapy. Our findings suggest that natalizumab diminishes the exchange of peripheral and intrathecal B cells without impacting intrathecal clonal expansion. In contrast, fingolimod treatment fails to alter blood-brain barrier B cell exchange but diminishes intrathecal clonal expansion. Sphingosine-1 phosphate receptor inhibition may alter intrathecal B cell biology in MS.

Immune Literacy: Reading, Writing, and Editing Adaptive Immunity

Advances in reading, writing, and editing DNA are providing unprecedented insights into the complexity of immunological systems. This combination of systems and synthetic biology methods is enabling the quantitative and precise understanding of molecular recognition in adaptive immunity, thus providing a framework for reprogramming immune responses for translational medicine. In this review, we will highlight state-of-the-art methods such as immune repertoire sequencing, immunoinformatics, and immunogenomic engineering and their application toward adaptive immunity. We showcase novel and interdisciplinary approaches that have the promise of transforming the design and breadth of molecular and cellular immunotherapies.

GLaMST: grow lineages along minimum spanning tree for b cell receptor sequencing data

Background

B cell affinity maturation enables B cells to generate high-affinity antibodies. This process involves somatic hypermutation of B cell immunoglobulin receptor (BCR) genes and selection by their ability to bind antigens. Lineage trees are used to describe this microevolution of B cell immunoglobulin genes. In a lineage tree, each node is one BCR sequence that mutated from the germinal center and each directed edge represents a single base mutation, insertion or deletion. In BCR sequencing data, the observed data only contains a subset of BCR sequences in this microevolution process. Therefore, reconstructing the lineage tree from experimental data requires algorithms to build the tree based on partially observed tree nodes.

Results

We developed a new algorithm named Grow Lineages along Minimum Spanning Tree (GLaMST), which efficiently reconstruct the lineage tree given observed BCR sequences that correspond to a subset of the tree nodes. Through comparison using simulated and real data, GLaMST outperforms existing algorithms in simulations with high rates of mutation, insertion and deletion, and generates lineage trees with smaller size and closer to ground truth according to tree features that highly correlated with selection pressure.

Conclusions

GLaMST outperforms state-of-art in reconstruction of the BCR lineage tree in both efficiency and accuracy. Integrating it into existing BCR sequencing analysis frameworks can significant improve lineage tree reconstruction aspect of the analysis.

Fc receptor-like 4 and 5 define human atypical memory B cells

Atypical memory B cells accumulate in chronic infections and autoimmune conditions, and commonly express FCRL4 and FCRL5, respective IgA and IgG receptors. We characterized memory cells from tonsils on the basis of both FCRL4 and FCRL5 expression, defining three subsets with distinct surface proteins and gene expression. Atypical FCRL4+FCRL5+ memory cells had the most discrete surface protein expression and were enriched in cell adhesion pathways, consistent with functioning as tissue-resident cells. Atypical FCRL4−FCRL5+ memory cells expressed transcription factors and immunoglobulin genes that suggest poised differentiation into plasma cells. Accordingly, the FCRL4−FCRL5+ memory subset was enriched in pathways responding to endoplasmic reticulum stress and IFN-γ. We reconstructed ongoing B-cell responses as lineage trees, providing crucial in vivo developmental context. Each memory subset typically maintained its lineage, denoting mechanisms enforcing their phenotypes. Classical FCRL4−FCRL5− memory cells were infrequently detected in lineage trees, suggesting the majority were in a quiescent state. FCRL4−FCRL5+ cells were the most represented memory subset in lineage trees, indicating robust participation in ongoing responses. Together, these differences suggest FCRL4 and FCRL5 are unlikely to be passive markers but rather active drivers of human memory B-cell development and function.

AncesTree: An interactive immunoglobulin lineage tree visualizer

High-throughput sequencing of human immunoglobulin genes allows analysis of antibody repertoires and the reconstruction of clonal lineage evolution. The study of antibodies (Abs) affinity maturation is of specific interest to understand the generation of Abs with high affinity or broadly neutralizing activities. Moreover, phylogenic analysis enables the identification of the key somatic mutations required to achieve optimal antigen binding. The Immcantation framework provides a start-to-finish set of analytical methods for high-throughput adaptive immune receptor repertoire sequencing (AIRR-Seq; Rep-Seq) data. Furthermore, Immcantation's Change-O package has developed IgPhyML, an algorithm designed to build specifically immunoglobulin (Ig) phylogenic trees. Meanwhile Phylip, an algorithm that has been originally developed for applications in ecology and macroevolution, can also be used for the phylogenic reconstruction of antibodies maturation pathway. To complement Ig lineages made by IgPhyML or Dnaml (Phylip), we developed AncesTree, a graphic user interface (GUI) that aims to give researchers the opportunity to interactively explore antibodies clonal evolution. AncesTree displays interactive immunoglobulins phylogenic tree, Ig related mutations and sequence alignments using additional information coming from specialized antibody tools. The GUI is a Java standalone application allowing interaction with Ig tree that can run under Windows, Linux and Mac OS.

IGHV sequencing reveals acquired N-glycosylation sites as a clonal and stable event during follicular lymphoma evolution

Follicular lymphoma B cells undergo continuous somatic hypermutation (SHM) of their immunoglobulin variable region genes, generating a heterogeneous tumor population. SHM introduces DNA sequences encoding N-glycosylation sites asparagine-X-serine/threonine (N-gly sites) within the V-region that are rarely found in normal B-cell counterparts. Unique attached oligomannoses activate B-cell receptor signaling pathways after engagement with calcium-dependent lectins expressed by tissue macrophages. This novel interaction appears critical for tumor growth and survival. To elucidate the significance of N-gly site presence and loss during ongoing SHM, we tracked site behavior during tumor evolution and progression in a diverse group of patients through next-generation sequencing. A hierarchy of subclones was visualized through lineage trees based on SHM semblance between subclones and their discordance from the germline sequence. We observed conservation of N-gly sites in more than 96% of subclone populations within and across diagnostic, progression, and transformation events. Rare N-gly-negative subclones were lost or negligible from successive events, in contrast to N-gly-positive subclones, which could additionally migrate between anatomical sites. Ongoing SHM of the N-gly sites resulted in subclones with different amino acid compositions across disease events, yet the vast majority of resulting DNA sequences still encoded for an N-gly site. The selection and expansion of only N-gly-positive subclones is evidence of the tumor cells' dependence on sites, despite the changing genomic complexity as the disease progresses. N-gly sites were gained in the earliest identified lymphoma cells, indicating they are an early and stable event of pathogenesis. Targeting the inferred mannose-lectin interaction holds therapeutic promise.

Intrathecal B-cell activation in LGI1 antibody encephalitis

Objective

To study intrathecal B-cell activity in leucine-rich, glioma-inactivated 1 (LGI1) antibody encephalitis. In patients with LGI1 antibodies, the lack of CSF lymphocytosis or oligoclonal bands and serum-predominant LGI1 antibodies suggests a peripherally initiated immune response. However, it is unknown whether B cells within the CNS contribute to the ongoing pathogenesis of LGI1 antibody encephalitis.

Methods

Paired CSF and peripheral blood (PB) mononuclear cells were collected from 6 patients with LGI1 antibody encephalitis and 2 patients with other neurologic diseases. Deep B-cell immune repertoire sequencing was performed on immunoglobulin heavy chain transcripts from CSF B cells and sorted PB B-cell subsets. In addition, LGI1 antibody levels were determined in CSF and PB.

Results

Serum LGI1 antibody titers were on average 127-fold higher than CSF LGI1 antibody titers. Yet, deep B-cell repertoire analysis demonstrated a restricted CSF repertoire with frequent extensive clusters of clonally related B cells connected to mature PB B cells. These clusters showed intensive mutational activity of CSF B cells, providing strong evidence for an independent CNS-based antigen-driven response in patients with LGI1 antibody encephalitis but not in controls.

Conclusions

Our results demonstrate that intrathecal immunoglobulin repertoire expansion is a feature of LGI1 antibody encephalitis and suggests a need for CNS-penetrant therapies.

B cells in rheumatoid arthritis synovial tissues encode focused antibody repertoires that include antibodies that stimulate macrophage TNF-α production

Rheumatoid arthritis (RA) is characterized by the production of anti-citrullinated protein antibodies (ACPAs). To gain insights into the relationship between ACPA-expressing B cells in peripheral blood (PB) and synovial tissue (ST), we sequenced the B cell repertoire in paired PB and ST samples from five individuals with established, ACPA+ RA. Bioinformatics analysis of paired heavy- and light-chain sequences revealed clonally-related family members shared between PB and ST. ST-derived antibody repertoires exhibited reduced diversity and increased normalized clonal family size compared to PB-derived repertoires. Functional characterization showed that seven recombinant antibodies (rAbs) expressed from subject-derived sequences from both compartments bound citrullinated antigens and immune complexes (ICs) formed using one ST-derived rAb stimulated macrophage TNF-α production. Our findings demonstrate B cell trafficking between PB and ST in subjects with RA and ST repertoires include B cells that encode ACPA capable of forming ICs that stimulate cellular responses implicated in RA pathogenesis.

Repertoire-wide phylogenetic models of B cell molecular evolution reveal evolutionary signatures of aging and vaccination

In order to produce effective antibodies, B cells undergo rapid somatic hypermutation (SHM) and selection for binding affinity to antigen via a process called affinity maturation. The similarities between this process and evolution by natural selection have led many groups to use phylogenetic methods to characterize the development of immunological memory, vaccination, and other processes that depend on affinity maturation. However, these applications are limited by the fact that most phylogenetic models are designed to be applied to individual lineages comprising genetically diverse sequences, while B cell repertoires often consist of hundreds to thousands of separate low-diversity lineages. Further, several features of affinity maturation violate important assumptions in standard phylogenetic models. Here, we introduce a hierarchical phylogenetic framework that integrates information from all lineages in a repertoire to more precisely estimate model parameters while simultaneously incorporating the unique features of SHM. We demonstrate the power of this repertoire-wide approach by characterizing previously undescribed phenomena in affinity maturation. First, we find evidence consistent with age-related changes in SHM hot-spot targeting. Second, we identify a consistent relationship between increased tree length and signs of increased negative selection, apparent in the repertoires of recently vaccinated subjects and those without any known recent infections or vaccinations. This suggests that B cell lineages shift toward negative selection over time as a general feature of affinity maturation. Our study provides a framework for undertaking repertoire-wide phylogenetic testing of SHM hypotheses and provides a means of characterizing dynamics of mutation and selection during affinity maturation.

Deep sequencing of B cell receptor repertoire

Immune repertoire is a collection of enormously diverse adaptive immune cells within an individual. As the repertoire shapes and represents immunological conditions, identification of clones and characterization of diversity are critical for understanding how to protect ourselves against various illness such as infectious diseases and cancers. Over the past several years, fast growing technologies for high throughput sequencing have facilitated rapid advancement of repertoire research, enabling us to observe the diversity of repertoire at an unprecedented level. Here, we focus on B cell receptor (BCR) repertoire and review approaches to B cell isolation and sequencing library construction. These experiments should be carefully designed according to BCR regions to be interrogated, such as heavy chain full length, complementarity determining regions, and isotypes. We also highlight preprocessing steps to remove sequencing and PCR errors with unique molecular index and bioinformatics techniques. Due to the nature of massive sequence variation in BCR, caution is warranted when interpreting repertoire diversity from error-prone sequencing data. Furthermore, we provide a summary of statistical frameworks and bioinformatics tools for clonal evolution and diversity. Finally, we discuss limitations of current BCR-seq technologies and future perspectives on advances in repertoire sequencing.

Persistence of intrathecal oligoclonal B cells and IgG in multiple sclerosis

In multiple sclerosis (MS), B cells are trafficking across the blood-brain barrier, but it is not known how this relates to the synthesis of oligoclonal IgG. We used quantitative mass spectrometry of oligoclonal bands and high-throughput sequencing of immunoglobulin heavy-chain variable transcripts to study the longitudinal B cell response in the cerebrospinal fluid (CSF) and blood of two MS patients. Twenty of 22 (91%) and 25 of 28 (89%) of oligoclonal band peptides persisted in samples collected 18 months apart, in spite of a dynamic exchange across the blood-CSF barrier of B lineage cells connecting to oligoclonal IgG.

Differential transcriptome and development of human peripheral plasma cell subsets

Human antibody-secreting cells (ASCs) triggered by immunization are globally recognized as CD19loCD38hiCD27hi. Yet, different vaccines give rise to antibody responses of different longevity, suggesting ASC populations are heterogeneous. We define circulating-ASC heterogeneity in vaccine responses using multicolor flow cytometry, morphology, VH repertoire, and RNA transcriptome analysis. We also tested differential survival using a human cell-free system that mimics the bone marrow (BM) microniche. In peripheral blood, we identified 3 CD19+ and 2 CD19- ASC subsets. All subsets contributed to the vaccine-specific responses and were characterized by in vivo proliferation and activation. The VH repertoire demonstrated strong oligoclonality with extensive interconnectivity among the 5 subsets and switched memory B cells. Transcriptome analysis showed separation of CD19+ and CD19- subsets that included pathways such as cell cycle, hypoxia, TNF-α, and unfolded protein response. They also demonstrated similar long-term in vitro survival after 48 days. In summary, vaccine-induced ASCs with different surface markers (CD19 and CD138) are derived from shared proliferative precursors yet express distinctive transcriptomes. Equal survival indicates that all ASC compartments are endowed with long-lived potential. Accordingly, in vivo survival of peripheral long-lived plasma cells may be determined in part by their homing and residence in the BM microniche.

Investigation of the estuarine stonefish (Synanceia horrida) venom composition

The Estuarine stonefish (Synanceia horrida) is recognised as one of the most venomous fish species in the world but the overall venom composition has yet to be investigated using in-depth transcriptomic and proteomic methods. To date, known venom components are restricted to a hyaluronidase and a large, pore-forming toxin known as Stonustoxin (SNTX). Transcriptomic sequencing of the venom gland resulted in over 170,000 contigs with only 0.4% that were homologous to putative venom proteins. Integration of the transcriptomic data with proteomic data from the S. horrida venom confirmed the hyaluronidase and SNTX to be present, together with several other protein families including major contributions from C-type lectins. Other protein families observed included peroxiredoxin and several minor protein families such as Golgi-associated plant pathogenesis related proteins, tissue pathway factor inhibitors, and Kazal-type serine protease inhibitors that, although not putative venom proteins, may contribute to the venom's adverse effects. BIOLOGICAL SIGNIFICANCE: Proteomic analysis of milked Synanceia horrida venom, paired with transcriptomic analysis of the venom gland tissue revealed for the first time the composition of one of the world's most dangerous fish venoms. The results demonstrate that the venom is relatively less complex compared to other well-studied venomous animals with a number of unique proteins not previously found in animal venoms.

Using Genotype Abundance to Improve Phylogenetic Inference

Modern biological techniques enable very dense genetic sampling of unfolding evolutionary histories, and thus frequently sample some genotypes multiple times. This motivates strategies to incorporate genotype abundance information in phylogenetic inference. In this article, we synthesize a stochastic process model with standard sequence-based phylogenetic optimality, and show that tree estimation is substantially improved by doing so. Our method is validated with extensive simulations and an experimental single-cell lineage tracing study of germinal center B cell receptor affinity maturation.

Affinity Maturation of the Anti-Citrullinated Protein Antibody Paratope Drives Epitope Spreading and Polyreactivity in Rheumatoid Arthritis

OBJECTIVE

Anti-citrullinated protein antibodies (ACPAs) are a hallmark of rheumatoid arthritis (RA). While epitope spreading of the serum ACPA response is believed to contribute to RA pathogenesis, little is understood regarding how this phenomenon occurs. This study was undertaken to analyze the antibody repertoires of individuals with RA to gain insight into the mechanisms leading to epitope spreading of the serum ACPA response in RA.

METHODS

Plasmablasts from the blood of 6 RA patients were stained with citrullinated peptide tetramers to identify ACPA-producing B cells by flow cytometry. Plasmablasts were single-cell sorted and sequenced to obtain antibody repertoires. Sixty-nine antibodies were recombinantly expressed, and their anticitrulline reactivities were characterized using a cyclic citrullinated peptide enzyme-linked immuosorbent assay and synovial antigen arrays. Thirty-six mutated antibodies designed either to represent ancestral antibodies or to test paratope residues critical for binding, as determined from molecular modeling studies, were also tested for anticitrulline reactivities.

RESULTS

Clonally related monoclonal ACPAs and their shared ancestral antibodies each exhibited differential reactivity against citrullinated antigens. Molecular modeling identified residues within the complementarity-determining region loops and framework regions predicted to be important for citrullinated antigen binding. Affinity maturation resulted in mutations of these key residues, which conferred binding to different citrullinated epitopes and/or increased polyreactivity to citrullinated epitopes.

CONCLUSION

These results demonstrate that the different somatic hypermutations accumulated by clonally related B cells during affinity maturation alter the antibody paratope to mediate epitope spreading and polyreactivity of the ACPA response in RA, suggesting that these may be key properties that likely contribute to the pathogenicity of ACPAs.

Affinity Maturation Drives Epitope Spreading and Generation of Proinflammatory Anti-Citrullinated Protein Antibodies in Rheumatoid Arthritis

OBJECTIVE

Rheumatoid arthritis (RA) is characterized by the presence of anti-citrullinated protein antibodies (ACPAs); nevertheless, the origin, specificity, and functional properties of ACPAs remain poorly understood. The aim of this study was to characterize the evolution of ACPAs by sequencing the plasmablast antibody repertoire at serial time points in patients with established RA.

METHODS

Blood samples were obtained at up to 4 serial time points from 8 individuals with established RA who were positive for ACPAs by the anti-cyclic citrullinated peptide test. CD19+CD3-IgD-CD14-CD20-CD27+CD38++ plasmablasts were isolated by single-cell sorting and costained with citrullinated peptide tetramers to identify ACPA-expressing plasmablasts. Cell-specific oligonucleotide barcodes were utilized, followed by large-scale sequencing and bioinformatics analysis, to obtain error-corrected, paired heavy- and light-chain antibody gene sequences for each B cell.

RESULTS

Bioinformatics analysis revealed 170 persistent plasmablast lineages in the RA blood, of which 19% included multiple isotypes. Among IgG- and IgA-expressing plasmablasts, significantly more IgA-expressing than IgG-expressing persistent lineages were observed (P < 0.01). Shared complementarity-determining region 3 sequence motifs were identified across subjects. A subset of the plasmablast lineages included members derived from later time points with divergent somatic hypermutations that encoded antibodies that bind an expanded set of citrullinated antigens. Furthermore, these recombinant, differentially mutated plasmablast antibodies formed immune complexes that stimulated higher macrophage production of tumor necrosis factor (TNF) compared to antibodies representing earlier time point-derived lineage members that were less mutated.

CONCLUSION

These findings demonstrate that established RA is characterized by a persistent IgA ACPA response that exhibits ongoing affinity maturation. This observation suggests the presence of a persistent mucosal antigen that continually promotes the production of IgA plasmablasts and their affinity maturation and epitope spreading, thus leading to the generation of ACPAs that bind additional citrullinated antigens and more potently stimulate macrophage production of TNF.

Benchmarking Tree and Ancestral Sequence Inference for B Cell Receptor Sequences

B cell receptor sequences evolve during affinity maturation according to a Darwinian process of mutation and selection. Phylogenetic tools are used extensively to reconstruct ancestral sequences and phylogenetic trees from affinity-matured sequences. In addition to using general-purpose phylogenetic methods, researchers have developed new tools to accommodate the special features of B cell sequence evolution. However, the performance of classical phylogenetic techniques in the presence of B cell-specific features is not well understood, nor how much the newer generation of B cell specific tools represent an improvement over classical methods. In this paper we benchmark the performance of classical phylogenetic and new B cell-specific tools when applied to B cell receptor sequences simulated from a forward-time model of B cell receptor affinity maturation toward a mature receptor. We show that the currently used tools vary substantially in terms of tree structure and ancestral sequence inference accuracy. Furthermore, we show that there are still large performance gains to be achieved by modeling the special mutation process of B cell receptors. These conclusions are further strengthened with real data using the rules of isotype switching to count possible violations within each inferred phylogeny.

Comparative Analysis of B-Cell Receptor Repertoires Induced by Live Yellow Fever Vaccine in Young and Middle-Age Donors

Age-related changes can significantly alter the state of adaptive immune system and often lead to attenuated response to novel pathogens and vaccination. In present study we employed 5′RACE UMI-based full length and nearly error-free immunoglobulin profiling to compare plasma cell antibody repertoires in young (19–26 years) and middle-age (45–58 years) individuals vaccinated with a live yellow fever vaccine, modeling a newly encountered pathogen. Our analysis has revealed age-related differences in the responding antibody repertoire ranging from distinct IGH CDR3 repertoire properties to differences in somatic hypermutation intensity and efficiency and antibody lineage tree structure. Overall, our findings suggest that younger individuals respond with a more diverse antibody repertoire and employ a more efficient somatic hypermutation process than elder individuals in response to a newly encountered pathogen.

Sequencing HIV-neutralizing antibody exons and introns reveals detailed aspects of lineage maturation

The developmental pathways of broadly neutralizing antibodies (bNAbs) against HIV are of great importance for the design of immunogens that can elicit protective responses. Here we show the maturation features of the HIV-neutralizing anti-V1V2 VRC26 lineage by simultaneously sequencing the exon together with the downstream intron of VRC26 members. Using the mutational landscapes of both segments and the selection-free nature of the intron region, we identify multiple events of amino acid mutational convergence in the complementarity-determining region 3 (CDR3) of VRC26 members, and determine potential intermediates with diverse CDR3s to a late stage bNAb from 2 years prior to its isolation. Moreover, we functionally characterize the earliest neutralizing intermediates with critical CDR3 mutations, with some emerging only 14 weeks after initial lineage detection and containing only ~6% V gene mutations. Our results thus underscore the utility of analyzing exons and introns simultaneously for studying antibody maturation and repertoire selection.

Knowledge on how antibody responses have evolved is critical for the induction of protective immunity. Here the authors analyse, using high-throughput sequencing of both exon and intron regions, the mutation and lineage development of an HIV-neutralizing antibody to find an unexpected early emergence of broadly neutralizing species.

Tracing Antibody Repertoire Evolution by Systems Phylogeny

Antibody evolution studies have been traditionally limited to either tracing a single clonal lineage (B cells derived from a single V-(D)-J recombination) over time or examining bulk functionality changes (e.g., tracing serum polyclonal antibody proteins). Studying a single B cell disregards the majority of the humoral immune response, whereas bulk functional studies lack the necessary resolution to analyze the co-existing clonal diversity. Recent advances in high-throughput sequencing (HTS) technologies and bioinformatics have made it possible to examine multiple co-evolving antibody monoclonal lineages within the context of a single repertoire. A plethora of accompanying methods and tools have been introduced in hopes of better understanding how pathogen presence dictates the global evolution of the antibody repertoire. Here, we provide a comprehensive summary of the tremendous progress of this newly emerging field of systems phylogeny of antibody responses. We present an overview encompassing the historical developments of repertoire phylogenetics, state-of-the-art tools, and an outlook on the future directions of this fast-advancing and promising field.

Comparison of methods for phylogenetic B-cell lineage inference using time-resolved antibody repertoire simulations (AbSim)

Motivation

The evolution of antibody repertoires represents a hallmark feature of adaptive B-cell immunity. Recent advancements in high-throughput sequencing have dramatically increased the resolution to which we can measure the molecular diversity of antibody repertoires, thereby offering for the first time the possibility to capture the antigen-driven evolution of B cells. However, there does not exist a repertoire simulation framework yet that enables the comparison of commonly utilized phylogenetic methods with regard to their accuracy in inferring antibody evolution.

Results

Here, we developed AbSim, a time-resolved antibody repertoire simulation framework, which we exploited for testing the accuracy of methods for the phylogenetic reconstruction of B-cell lineages and antibody molecular evolution. AbSim enables the (i) simulation of intermediate stages of antibody sequence evolution and (ii) the modeling of immunologically relevant parameters such as duration of repertoire evolution, and the method and frequency of mutations. First, we validated that our repertoire simulation framework recreates replicates topological similarities observed in experimental sequencing data. Second, we leveraged Absim to show that current methods fail to a certain extent to predict the true phylogenetic tree correctly. Finally, we formulated simulation-validated guidelines for antibody evolution, which in the future will enable the development of accurate phylogenetic methods.

Availability and implementation

https://cran.r-project.org/web/packages/AbSim/index.html.

Contact

sai.reddy@ethz.ch.

Supplementary information

Supplementary data are available at Bioinformatics online.

The Bayesian optimist's guide to adaptive immune receptor repertoire analysis

Probabilistic modeling is fundamental to the statistical analysis of complex data. In addition to forming a coherent description of the data-generating process, probabilistic models enable parameter inference about given datasets. This procedure is well developed in the Bayesian perspective, in which one infers probability distributions describing to what extent various possible parameters agree with the data. In this paper, we motivate and review probabilistic modeling for adaptive immune receptor repertoire data then describe progress and prospects for future work, from germline haplotyping to adaptive immune system deployment across tissues. The relevant quantities in immune sequence analysis include not only continuous parameters such as gene use frequency but also discrete objects such as B-cell clusters and lineages. Throughout this review, we unravel the many opportunities for probabilistic modeling in adaptive immune receptor analysis, including settings for which the Bayesian approach holds substantial promise (especially if one is optimistic about new computational methods). From our perspective, the greatest prospects for progress in probabilistic modeling for repertoires concern ancestral sequence estimation for B-cell receptor lineages, including uncertainty from germline genotype, rearrangement, and lineage development.

Computational Evaluation of B-Cell Clone Sizes in Bulk Populations

B cell clones expand and contract during adaptive immune responses and can persist or grow uncontrollably in lymphoproliferative disorders. One way to monitor and track B cell clones is to perform large-scale sampling of bulk cell populations, amplifying, and sequencing antibody gene rearrangements by next-generation sequencing (NGS). Here, we describe a series of computational approaches for estimating B cell clone size in NGS immune repertoire profiling data of antibody heavy chain gene rearrangements. We define three different measures of B cell clone size-copy numbers, instances, and unique sequences-and show how these measures can be used to rank clones, analyze their diversity, and study their distribution within and between individuals. We provide a detailed, step-by-step procedure for performing these analyses using two different data sets of spleen samples from human organ donors. In the first data set, 19 independently generated biological replicates from a single individual are analyzed for B cell clone size, diversity and sampling sufficiency for clonal overlap analysis. In the second data set, B cell clones are compared in eight different organ donors. We comment upon frequently encountered pitfalls and offer practical advice with alternative approaches. Overall, we provide a series of pragmatic analytical approaches and show how different clone size measures can be used to study the clonal landscape in bulk B cell immune repertoire profiling data.

Acquisition of N-Glycosylation Sites in Immunoglobulin Heavy Chain Genes During Local Expansion in Parotid Salivary Glands of Primary Sjögren Patients

Previous studies revealed high incidence of acquired N-glycosylation sites acquired N-glycosylation sites in RNA transcripts encoding immunoglobulin heavy variable region (IGHV) 3 genes from parotid glands of primary Sjögren’s syndrome (pSS) patients. In this study, next generation sequencing was used to study the extent of ac-Nglycs among clonally expanded cells from all IGVH families in the salivary glands of pSS patients. RNA was isolated from parotid gland biopsies of five pSS patients and five non-pSS sicca controls. IGHV sequences covering all functional IGHV genes were amplified, sequenced, and analyzed. Each biopsy recovered 1,800–4,000 unique IGHV sequences. No difference in IGHV gene usage was observed between pSS and non-pSS sequences. Clonally related sequences with more than 0.3% of the total number of sequences per patient were referred to as dominant clone. Overall, 70 dominant clones were found in pSS biopsies, compared to 15 in non-pSS. No difference in percentage mutation in dominant clone-derived IGHV sequences was seen between pSS and non-pSS. In pSS, no evidence for antigen-driven selection in dominant clones was found. We observed a significantly higher amount of ac-Nglycs among pSS dominant clone-derived sequences compared to non-pSS. Ac-Nglycs were, however, not restricted to dominant clones or IGHV gene. Most ac-Nglycs were detected in the framework 3 region. No stereotypic rheumatoid factor rearrangements were found in dominant clones. Lineage tree analysis showed in four pSS patients, but not in non-pSS, the presence of the germline sequence from a dominant clone. Presence of germline sequence and mutated IGHV sequences in the same dominant clone provide evidence that this clone originated from a naïve B-cell recruited into the parotid gland to expand and differentiate locally into plasma cells. The increased presence of ac-Nglycs in IGHV sequences, due to somatic hypermutation, might provide B-cells an escape mechanism to survive during immune response. We speculate that glycosylation of the B-cell receptor makes the cell sensitive to environmental lectin signals to contribute to aberrant B-cell selection in pSS parotid glands.

Computational Strategies for Dissecting the High-Dimensional Complexity of Adaptive Immune Repertoires

The adaptive immune system recognizes antigens via an immense array of antigen-binding antibodies and T-cell receptors, the immune repertoire. The interrogation of immune repertoires is of high relevance for understanding the adaptive immune response in disease and infection (e.g., autoimmunity, cancer, HIV). Adaptive immune receptor repertoire sequencing (AIRR-seq) has driven the quantitative and molecular-level profiling of immune repertoires, thereby revealing the high-dimensional complexity of the immune receptor sequence landscape. Several methods for the computational and statistical analysis of large-scale AIRR-seq data have been developed to resolve immune repertoire complexity and to understand the dynamics of adaptive immunity. Here, we review the current research on (i) diversity, (ii) clustering and network, (iii) phylogenetic, and (iv) machine learning methods applied to dissect, quantify, and compare the architecture, evolution, and specificity of immune repertoires. We summarize outstanding questions in computational immunology and propose future directions for systems immunology toward coupling AIRR-seq with the computational discovery of immunotherapeutics, vaccines, and immunodiagnostics.

Studying Antibody Repertoires with Next-Generation Sequencing

Next-generation sequencing is making it possible to study the antibody repertoire of an organism in unprecedented detail, and, by so doing, to characterize its behavior in the response to infection and in pathological conditions such as autoimmunity and cancer. The polymorphic nature of the repertoire poses unique challenges that rule out the use of many commonly used NGS methods and require tradeoffs to be made when considering experimental design.We outline the main contexts in which antibody repertoire analysis has been used, and summarize the key tools that are available. The humoral immune response to vaccination has been a particular focus of repertoire analyses, and we review the key conclusions and methods used in these studies.

Human Secretory IgM Emerges from Plasma Cells Clonally Related to Gut Memory B Cells and Targets Highly Diverse Commensals

Secretory immunoglobulin A (SIgA) enhances host-microbiota symbiosis, whereas SIgM remains poorly understood. We found that gut IgM⁺ plasma cells (PCs) were more abundant in humans than mice and clonally related to a large repertoire of memory IgM⁺ B cells disseminated throughout the intestine but rare in systemic lymphoid organs. In addition to sharing a gut-specific gene signature with memory IgA⁺ B cells, memory IgM⁺ B cells were related to some IgA⁺ clonotypes and switched to IgA in response to T cell-independent or T cell-dependent signals. These signals induced abundant IgM which, together with SIgM from clonally affiliated PCs, recognized mucus-embedded commensals. Bacteria recognized by human SIgM were dually coated by SIgA and showed increased richness and diversity compared to IgA-only-coated or uncoated bacteria. Thus, SIgM may emerge from pre-existing memory rather than newly activated naive IgM⁺ B cells and could help SIgA to anchor highly diverse commensal communities to mucus.

•

IgM⁺ PCs generating SIgM are relatively abundant in human but not mouse gut

•

IgM⁺ PCs clonally relate to a large gut repertoire of memory IgM⁺ B cells

•

Gut memory IgM⁺ B cells express a tissue-specific signature and can switch to IgA

•

Human but not mouse SIgM binds a highly diverse microbiota dually coated by SIgA

Genetic background and evolution of relapses in aggressive B-cell lymphomas

Relapses of aggressive B-cell lymphomas pose a higher risk to affected patients because of potential treatment resistance and usually rapid tumor growth. Recent advances, such as targeting Bruton tyrosine kinase, have provided promising results in small numbers of cases, but treatment for the majority of patients remains challenging and outcomes are generally poor. A number of recent studies have utilized state-of-the-art genomic technologies in an attempt to better understand tumor genome evolution during relapse and to identify relapse-specific genetic alterations. It has been found that in some settings (e.g. diffuse large B-cell lymphomas in immunocompromised patients, secondary diffuse large B-cell lymphomas as Richter transformations) a significant part of the recurrences are clonally-unrelated de novo neoplasms, which might have distinct genomic and drug sensitivity profiles as well as different prognoses. Similar to earlier findings in indolent lymphomas, genetic tumor evolution of clonally-related relapsing aggressive B-cell lymphomas is predominantly characterized by two patterns: early divergence from a common progenitor and late divergence/linear evolution from a primary tumor. The clinical implications of these distinct patterns are not yet clear and will require additional investigation; however, it is plausible that these two patterns of recurrence are linked to different treatment-resistance mechanisms. Attempts to identify drivers of relapses result in a very heterogeneous list of affected genes and pathways as well as epigenetic mechanisms and suggest many ways of how recurrent tumors can adapt to treatment and expand their malignant properties.

CNS Aquaporin-4-specific B cells connect with multiple B-cell compartments in neuromyelitis optica spectrum disorder

Objectives

Neuromyelitis optica spectrum disorder (NMOSD) is a severe inflammatory disorder of the central nervous system (CNS) targeted against aquaporin‐4 (AQP4). The origin and trafficking of AQP4‐specific B cells in NMOSD remains unknown.

Methods

Peripheral (n = 7) and splenic B cells (n = 1) recovered from seven NMOSD patients were sorted into plasmablasts, naïve, memory, and CD27‐IgD‐ double negative (DN) B cells, and variable heavy chain (VH) transcriptome sequences were generated by deep sequencing. Peripheral blood (PB) VH repertoires were compared to the same patient's single‐cell cerebrospinal fluid (CSF) plasmablast (PB) VH transcriptome, CSF immunoglobulin (Ig) proteome, and serum Ig proteome. Recombinant antibodies were generated from paired CSF heavy‐ and light chains and tested for AQP4 reactivity.

Results

Approximately 9% of the CSF VH sequences aligned with PB memory B cells, DN B cells, and plasmablast VH sequences. AQP4‐specific VH sequences were observed in each peripheral B‐cell compartment. Lineage analysis of clonally related VH sequences indicates that CSF AQP4‐specific B cells are closely related to an expanded population of DN B cells that may undergo antigen‐specific B‐cell maturation within the CNS. CSF and serum Ig proteomes overlapped with the VH sequences from each B‐cell compartment; the majority of matches occurring between the PB VH sequences and serum Ig proteome.

Interpretation

During an acute NMOSD relapse, a dynamic exchange of B cells occurs between the periphery and CNS with AQP4‐specific CSF B cells emerging from postgerminal center memory B cells and plasmablasts. Expansion of the PB DN B‐cell compartment may be a potential biomarker of NMOSD activity.