Trials and Tribulations with VH Replacement

Systematic evaluation of B-cell clonal family inference approaches

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

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.

Understanding the human antibody repertoire

The origins of the various elements in the human antibody repertoire have been and still are subject to considerable uncertainty. Uncertainty in respect of whether the various elements have always served a specific defense function or whether they were co-opted from other organismal roles to form a crude naïve repertoire that then became more complex as combinatorial mechanisms were added. Estimates of the current size of the human antibody naïve repertoire are also widely debated with numbers anywhere from 10 million members, based on experimentally derived numbers, to in excess of one thousand trillion members or more, based on the different sequences derived from theoretical combinatorial calculations. There are questions that are relevant at both ends of this number spectrum. At the lower bound it could be questioned whether this is an insufficient repertoire size to counter all the potential antigen-bearing pathogens. At the upper bound the question is rather simpler: How can any individual interrogate such an astronomical number of antibody-bearing B cells in a timeframe that is meaningful? This review evaluates the evolutionary aspects of the adaptive immune system, the calculations that lead to the large repertoire estimates, some of the experimental evidence pointing to a more restricted repertoire whose variation appears to derive from convergent 'structure and specificity features', and includes a theoretical model that seems to support it. Finally, a solution that may reconcile the size difference anomaly, which is still a hot subject of debate, is suggested.

In situ detection of PR3-ANCA+ B cells and alterations in the variable region of immunoglobulin genes support a role of inflamed tissue in the emergence of auto-reactivity in granulomatosis with polyangiitis

Circulating anti-neutrophilic cytoplasmic autoantibodies targeting proteinase 3 (PR3-ANCA) are a diagnostic and pathogenic hallmark of granulomatosis with polyangiitis (GPA). It is, however, incompletely understood if inflamed tissue supports presence or emergence of PR3-ANCA⁺ B cells. In search of such cells in inflamed tissue of GPA, immunofluorescence staining for IgG and a common PR3-ANCA idiotype (5/7 Id) was undertaken. Few 5/7 Id⁺/IgG⁺ B cells were detected in respiratory and kidney tissue of GPA. To gain more insight into surrogate markers possibly indicative of an anti-PR3-response, a meta-analysis comprising IGVH and IGVL genes derived from respiratory tract tissue of GPA (231 clones) was performed. Next generation sequencing-based IGHV genes derived from peripheral blood of healthy donors (244.353 clones) and previously published IGLV genes (148 clones) served as controls. Additionally, Ig genes of three murine and five known human monoclonal anti-PR3 antibodies were analyzed. Primary and probably secondary rearrangements led to altered VDJ usage and an extended complementarity determining region 3 (CDR3) of IGHV clones from GPA tissue. Selection against amino acid exchanges was prominent in the framework region of IGHV clones from GPA tissue. The comparison of V(D)J rearrangements and deduced amino acid sequences of the CDR3 yielded no identities and few similarities between clones derived from respiratory tissue of GPA and anti-PR3 antibodies, arguing against a presence of B cells that carry PR3-ANCA-prone Ig genes among the clones. In line with the scarcity of 5/7 Id⁺ B lymphocytes in GPA tissue, the results suggest that with respect to a local anti-PR3 response, methods detecting rare clones are required.

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.

Circulating Plasmablasts from Chronically Human Immunodeficiency Virus-Infected Individuals Predominantly Produce Polyreactive/Autoreactive Antibodies

Understanding the B-cell response during chronic human immunodeficiency virus (HIV) infection is essential for eliciting broad and potent neutralizing antibodies (Abs). In this study, we analyzed the plasmablast repertoire of chronically HIV-infected individuals in combination with antiretroviral therapy (ART). Among the obtained 72 recombinant monoclonal antibodies (mAbs), 27.8% weakly bound to HIV gp140 and were non-neutralizing. Remarkably, 56.9% were polyreactive and 55.6% were autoreactive. The prominent feature of being polyreactive/autoreactive is not limited to anti-gp140 Abs. Furthermore, these polyreactive/autoreactive Abs displayed striking cross-reactivity with DWEYS in the N-methyl-d-aspartate receptor (NMDAR), and this binding induced SH-SY5Y cell apoptosis. We also found higher frequencies of VH4-34 utilization and VH replacement in the plasmablast repertoire of chronically HIV-infected individuals, which may contribute to the generation of poly/autoreactive Abs. Taken together, these data demonstrate that circulating plasmablasts in chronically HIV-infected individuals experienced with ART predominantly produce poly/autoreactive Abs with minimal anti-HIV neutralizing capacity and potential cross-reactivity with autoantigens. This may represent another dysfunction of B cells during chronic HIV infection.

An atlas of B-cell clonal distribution in the human body

B-cell responses result in clonal expansion, and can occur in a variety of tissues. To define how B-cell clones are distributed in the body, we sequenced 933,427 B-cell clonal lineages and mapped them to eight different anatomic compartments in six human organ donors. We show that large B-cell clones partition into two broad networks-one spans the blood, bone marrow, spleen and lung, while the other is restricted to tissues within the gastrointestinal (GI) tract (jejunum, ileum and colon). Notably, GI tract clones display extensive sharing of sequence variants among different portions of the tract and have higher frequencies of somatic hypermutation, suggesting extensive and serial rounds of clonal expansion and selection. Our findings provide an anatomic atlas of B-cell clonal lineages, their properties and tissue connections. This resource serves as a foundation for studies of tissue-based immunity, including vaccine responses, infections, autoimmunity and cancer.

Mechanisms That Shape Human Antibody Repertoire Development in Mice Transgenic for Human Ig H and L Chain Loci

To determine the impact of the milieu on the development of the human B cell repertoire, we carried out a comprehensive analysis of productive and nonproductive Ig gene rearrangements from transgenic mice engineered to express single copies of the unrearranged human H chain and L chain Ig gene loci. By examining the nonproductive repertoire as an indication of the immediate product of the rearrangement machinery without an impact of selection, we discovered that the distribution of human rearrangements arising in the mouse was generally comparable to that seen in humans. However, differences between the distribution of nonproductive and productive rearrangements that reflect the impact of selection suggested species-specific selection played a role in shaping the respective repertoires. Although expression of some VH genes was similar in mouse and human (IGHV3-23, IGHV3-30, and IGHV4-59), other genes behaved differently (IGHV3-33, IGHV3-48, IGHV4-31, IGHV4-34, and IGHV1-18). Gene selection differences were also noted in L chains. Notably, nonproductive human VH rearrangements in the transgenic mice expressed shorter CDRH3 with less N addition. Even the CDRH3s in the productive rearrangements were shorter in length than those of the normal human productive repertoire. Amino acids in the CDRH3s in both species showed positive selection of tyrosines and glycines, and negative selection of leucines. The data indicate that the environment in which B cells develop can affect the expressed Ig repertoire by exerting influences on the distribution of expressed VH and VL genes and by influencing the amino acid composition of the Ag binding site.

The analysis of clonal expansions in normal and autoimmune B cell repertoires

Clones are the fundamental building blocks of immune repertoires. The number of different clones relates to the diversity of the repertoire, whereas their size and sequence diversity are linked to selective pressures. Selective pressures act both between clones and within different sequence variants of a clone. Understanding how clonal selection shapes the immune repertoire is one of the most basic questions in all of immunology. But how are individual clones defined? Here we discuss different approaches for defining clones, starting with how antibodies are diversified during different stages of B cell development. Next, we discuss how clones are defined using different experimental methods. We focus on high-throughput sequencing datasets, and the computational challenges and opportunities that these data have for mining the antibody repertoire landscape. We discuss methods that visualize sequence variants within the same clone and allow us to consider collections of shared mutations to determine which sequences share a common ancestry. Finally, we comment on features of frequently encountered expanded B cell clones that may be of particular interest in the setting of autoimmunity and other chronic conditions.

Discrimination of germline V genes at different sequencing lengths and mutational burdens: A new tool for identifying and evaluating the reliability of V gene assignment

Immune repertoires are collections of lymphocytes that express diverse antigen receptor gene rearrangements consisting of Variable (V), (Diversity (D) in the case of heavy chains) and Joining (J) gene segments. Clonally related cells typically share the same germline gene segments and have highly similar junctional sequences within their third complementarity determining regions. Identifying clonal relatedness of sequences is a key step in the analysis of immune repertoires. The V gene is the most important for clone identification because it has the longest sequence and the greatest number of sequence variants. However, accurate identification of a clone's germline V gene source is challenging because there is a high degree of similarity between different germline V genes. This difficulty is compounded in antibodies, which can undergo somatic hypermutation. Furthermore, high-throughput sequencing experiments often generate partial sequences and have significant error rates. To address these issues, we describe a novel method to estimate which germline V genes (or alleles) cannot be discriminated under different conditions (read lengths, sequencing errors or somatic hypermutation frequencies). Starting with any set of germline V genes, this method measures their similarity using different sequencing lengths and calculates their likelihood of unambiguous assignment under different levels of mutation. Hence, one can identify, under different experimental and biological conditions, the germline V genes (or alleles) that cannot be uniquely identified and bundle them together into groups of specific V genes with highly similar sequences.

Virus-induced preferential antibody gene-usage and its importance in humoral autoimmunity

It is known that even the adaptive components of the immune system are based on genetic traits common to all individuals, and that diversity is shaped by the lifelong contacts with different non-self antigens, including those found on infectious pathogens. Besides the individual differences, some of these common traits may be more prone to react against a given antigen, and this may be exploited by the infectious pathogens. Indeed, viral infections can deregulate immune response by subverting antibody (Ab) gene usage, leading to the overexpression of specific Ab subfamilies. This overexpression often results in a protective antiviral response but, in some cases, also correlates with a higher likelihood of developing humoral autoimmune disorders. These aspects of virus-induced autoimmunity have never been thoroughly reviewed, and this is the main purpose of this review. An accurate examination of virus specific Ab subfamilies elicited during infections may help further characterize the complex interplay between viruses and the humoral immune response, and be useful in the design of future monoclonal antibody (mAb)-based anti-infective prophylactic and therapeutic strategies.

Antibody repertoire diversification through VH gene replacement in mice cloned from an IgA plasma cell

In mammals, VDJ recombination is responsible for the establishment of a highly diversified preimmune antibody repertoire. Acquisition of a functional Ig heavy (H) chain variable (V) gene rearrangement is thought to prevent further recombination at the IgH locus. Here, we describe VHQ52(NT); Vκgr32(NT) Ig monoclonal mice reprogrammed from the nucleus of an intestinal IgA(+) plasma cell. In VHQ52(NT) mice, IgA replaced IgM to drive early B-cell development and peripheral B-cell maturation. In VHQ52(NT) animals, over 20% of mature B cells disrupted the single productive, nonautoimmune IgH rearrangement through VH replacement and exchanged it with a highly diversified pool of IgH specificities. VH replacement occurred in early pro-B cells, was independent of pre-B-cell receptor signaling, and involved predominantly one adjacent VH germ-line gene. VH replacement was also identified in 5% of peripheral B cells of mice inheriting a different productive VH rearrangement expressed in the form of an IgM H chain. In summary, editing of a productive IgH rearrangement through VH replacement can account for up to 20% of the IgH repertoire expressed by mature B cells.

Immunologic Basis for Long HCDR3s in Broadly Neutralizing Antibodies Against HIV-1

A large number of potent broadly neutralizing antibodies (bnAbs) against HIV-1 have been reported in recent years, raising hope for the possibility of an effective vaccine based on epitopes recognized by these protective antibodies. However, many of these bnAbs contain the long heavy chain complementarity-determining region 3 (HCDR3), which is viewed as an obstacle to the development of an HIV-1 vaccine targeting the bnAb responses. This mini-review summarizes the current literature and discusses the different potential immunologic mechanisms for generating long HCDR3, including D–D fusion, VH replacement, long N region addition, and skewed D–J gene usage, among which potential VH replacement products appear to be significant contributors. VH replacement occurs through recombinase activated gene-mediated secondary recombination and contributes to the diversified naïve B cell repertoire. During VH replacement, a short stretch of nucleotides from previously rearranged VH genes remains within the newly formed HCDR3, thus elongating its length. Accumulating evidence suggests that long HCDR3s are present in significant numbers in the human mature naïve B cell repertoire and are primarily generated by recombination during B cell development. These new observations indicate that long HCDR3s, though low in frequency, are a normal feature of the human antibody naïve repertoire and they appear to be selected to target conserved epitopes located in deep, partially obscured regions of the HIV-1 envelope trimer. Therefore, the presence of long HCDR3 sequences should not necessarily be viewed as an obstacle to the development of an HIV-1 vaccine based upon bnAb responses.