The CDR1 sequences of a major proportion of human germline Ig VH genes are inherently susceptible to amino acid replacement

A humanized mouse that mounts mature class-switched, hypermutated and neutralizing antibody responses

Humanized mice are limited in terms of modeling human immunity, particularly with regards to antibody responses. Here we constructed a humanized (THX) mouse by grafting non-γ-irradiated, genetically myeloablated KitW-41J mutant immunodeficient pups with human cord blood CD34+ cells, followed by 17β-estradiol conditioning to promote immune cell differentiation. THX mice reconstitute a human lymphoid and myeloid immune system, including marginal zone B cells, germinal center B cells, follicular helper T cells and neutrophils, and develop well-formed lymph nodes and intestinal lymphoid tissue, including Peyer's patches, and human thymic epithelial cells. These mice have diverse human B cell and T cell antigen receptor repertoires and can mount mature T cell-dependent and T cell-independent antibody responses, entailing somatic hypermutation, class-switch recombination, and plasma cell and memory B cell differentiation. Upon flagellin or a Pfizer-BioNTech coronavirus disease 2019 (COVID-19) mRNA vaccination, THX mice mount neutralizing antibody responses to Salmonella or severe acute respiratory syndrome coronavirus 2 Spike S1 receptor-binding domain, with blood incretion of human cytokines, including APRIL, BAFF, TGF-β, IL-4 and IFN-γ, all at physiological levels. These mice can also develop lupus autoimmunity after pristane injection. By leveraging estrogen activity to support human immune cell differentiation and maturation of antibody responses, THX mice provide a platform to study the human immune system and to develop human vaccines and therapeutics.

The evolutionary and functional significance of germline immunoglobulin gene variation

The recombination between immunoglobulin (IG) gene segments determines an individual's naïve antibody repertoire and, consequently, (auto)antigen recognition. Emerging evidence suggests that mammalian IG germline variation impacts humoral immune responses associated with vaccination, infection, and autoimmunity - from the molecular level of epitope specificity, up to profound changes in the architecture of antibody repertoires. These links between IG germline variants and immunophenotype raise the question on the evolutionary causes and consequences of diversity within IG loci. We discuss why the extreme diversity in IG loci remains a mystery, why resolving this is important for the design of more effective vaccines and therapeutics, and how recent evidence from multiple lines of inquiry may help us do so.

The signs of negative selection in IGHV framework regions are associated with worse overall survival of chronic lymphocytic leukemia patients

The mutational status of the variable region of the immunoglobulin heavy chain (IGHV) genes remains the most significant prognostic factor in chronic lymphocytic leukemia (CLL) patients. However, the groups of mutated (M) and unmutated (UM) patients are also heterogeneous, and additional markers are used for a more accurate prognosis. The aim of our work was to determine the prognostic value of the signs of antigen selection determined by BASELINe statistics in M IGHV sequences of CLL patients. Clinical data, IGHV gene configuration, TP53, NOTCH1, SF3B1 mutations were analyzed in 127 CLL patients with M IGHV sequences. The median OS of patients with negative selection in the framework regions (FWRs) of IGHV genes was 120 months compared to 202 month in other CLL patients (P = 0.016). In multivariate Cox regression analysis Binet stage C vs A + B (P < 0.0001), SF3B1 mutations (P < 0.0001), negative selection in the FWRs (HR P = 0.007), and age ≥65 years (P = 0.034) were powerful adverse prognostic factors for OS in CLL patients with M IGHV genes. These preliminary data suggest that the signs of antigen-driven selection may be used as a prognostic factor in CLL patients with M IGHV genes in combination with other markers.

LUBAC Suppresses IL-21-Induced Apoptosis in CD40-Activated Murine B Cells and Promotes Germinal Center B Cell Survival and the T-Dependent Antibody Response

B cell activation by Tfh cells, i.e., through CD154 engagement of CD40 and IL-21, and survival within GCs are crucial for the T-dependent Ab response. LUBAC, composed of HOIP, SHARPIN, and HOIL-1, catalyzes linear ubiquitination (Linear M1-Ub) to mediate NF-κB activation and cell survival induced by TNF receptor superfamily members, which include CD40. As shown in this study, B cells expressing the Sharpin null mutation cpdm (Sharpin^cpdm) could undergo proliferation, CSR, and SHM in response to immunization by a T-dependent Ag, but were defective in survival within GCs, enrichment of a mutation enhancing the BCR affinity, and production of specific Abs. Sharpin^cpdm B cells stimulated in vitro with CD154 displayed normal proliferation and differentiation, marginally impaired NF-κB activation and survival, but markedly exacerbated death triggered by IL-21. While activating the mitochondria-dependent apoptosis pathway in both Sharpin^+/+ and Sharpin^cpdm B cells, IL-21 induced Sharpin^cpdm B cells to undergo sustained activation of caspase 9 and caspase 8 of the mitochondria-dependent and independent pathway, respectively, and ultimately caspase 3 in effecting apoptosis. These were associated with loss of the caspase 8 inhibitor cFLIP and reduction in cFLIP Linear M1-Ub, which interferes with cFLIP poly-ubiquitination at Lys48 and degradation. Finally, the viability of Sharpin^cpdm B cells was rescued by caspase inhibitors but virtually abrogated – together with Linear M1-Ub and cFLIP levels – by a small molecule HOIP inhibitor. Thus, LUBAC controls the cFLIP expression and inhibits the effects of caspase 8 and IL-21-activated caspase 9, thereby suppressing apoptosis of CD40 and IL-21-activated B cells and promoting GC B cell survival.

Lost structural and functional inter-relationships between Ig and TCR loci in mammals revealed in sharks

Immunoglobulins and T cell receptors (TCR) have obvious structural similarities as well as similar immunogenetic diversification and selection mechanisms. Nevertheless, the two receptor systems and the loci that encode them are distinct in humans and classical murine models, and the gene segments comprising each repertoire are mutually exclusive. Additionally, while both B and T cells employ recombination-activating genes (RAG) for primary diversification, immunoglobulins are afforded a supplementary set of activation-induced cytidine deaminase (AID)-mediated diversification tools. As the oldest-emerging vertebrates sharing the same adaptive B and T cell receptor systems as humans, extant cartilaginous fishes allow a potential view of the ancestral immune system. In this review, we discuss breakthroughs we have made in studies of nurse shark (Ginglymostoma cirratum) T cell receptors demonstrating substantial integration of loci and diversification mechanisms in primordial B and T cell repertoires. We survey these findings in this shark model where they were first described, while noting corroborating examples in other vertebrate groups. We also consider other examples where the gnathostome common ancestry of the B and T cell receptor systems have allowed dovetailing of genomic elements and AID-based diversification approaches for the TCR. The cartilaginous fish seem to have retained this T/B cell plasticity to a greater extent than more derived vertebrate groups, but representatives in all vertebrate taxa except bony fish and placental mammals show such plasticity.

Human C. difficile toxin-specific memory B cell repertoires encode poorly neutralizing antibodies

Clostridioides difficile is a leading cause of nosocomial infection responsible for significant morbidity and mortality with limited options for therapy. Secreted C. difficile toxin B (TcdB) is a major contributor to disease pathology, and select TcdB-specific Abs may protect against disease recurrence. However, the high frequency of recurrence suggests that the memory B cell response, essential for new Ab production following C. difficile reexposure, is insufficient. We therefore isolated TcdB-specific memory B cells from individuals with a history of C. difficile infection and performed single-cell deep sequencing of their Ab genes. Herein, we report that TcdB-specific memory B cell–encoded antibodies showed somatic hypermutation but displayed limited isotype class switch. Memory B cell–encoded mAb generated from the gene sequences revealed low to moderate affinity for TcdB and a limited ability to neutralize TcdB. These findings indicate that memory B cells are an important factor in C. difficile disease recurrence.

The C. difficile toxin-specific human memory B cell repertoire encodes low-affinity, non-neutralizing antibodies.

Grade I, II and III Follicular Lymphomas Express Ig VH Genes with Different Patterns of Somatic Mutation

Follicular lymphoma (FL) is an indolent, B-cell, non-Hodgkin’s lymphoma with varying cytological appearance and clinical behavior. The genetic hallmark of FL is the t(14;18) translocation, and as a germinal center derived entity it is also characterized by somatic hypermutation of the immunoglobulin heavy chain (IgH) gene. In an attempt to correlate this molecular signature with the cytological grading of FL, we have analyzed the IgH variable (IgV_H), regions in all cytological grades of FL. Four FL cases showing t(14;18) translocation were classified into grade I-III categories according to the current WHO guidelines. The IgV_H gene segments were PCR-amplified, sequenced, and compared to their respective germline IgV_H sequences. The neoplastic cells of grade I and II FLs revealed clonally related, but highly divergent IgV_H gene sequences indicating the ongoing nature of somatic hypermutation. Grade III FL also showed extensive presence of somatic hypermutation, but these mutations were not associated with intraclonal divergence. Thus, these results suggest that grade I-II and grade III FL may represent different biological entities. The presence of ongoing somatic hypermutation of IgV_H sequences in grade I and II FLs is compatible with direct follicular origin of these tumor cells, contrasting the homogenous, stable clones of grade III FL resembling a post-follicular stage of B-cell development. Our findings demonstrate that contrary to the three tiered cytological grading, molecular features of IgH genes classify FL into two distinct subcategories. These studies also suggest that with progression FL gains post-follicular–like molecular features and becomes independent of the germinal center microenvironment.

AID Overlapping and Polη Hotspots Are Key Features of Evolutionary Variation Within the Human Antibody Heavy Chain (IGHV) Genes

Somatic hypermutation (SHM) of the immunoglobulin variable (IgV) loci is a key process in antibody affinity maturation. The enzyme activation-induced deaminase (AID), initiates SHM by creating C → U mismatches on single-stranded DNA (ssDNA). AID has preferential hotspot motif targets in the context of WRC/GYW (W = A/T, R = A/G, Y = C/T) and particularly at WGCW overlapping hotspots where hotspots appear opposite each other on both strands. Subsequent recruitment of the low-fidelity DNA repair enzyme, Polymerase eta (Polη), during mismatch repair, creates additional mutations at WA/TW sites. Although there are more than 50 functional immunoglobulin heavy chain variable (IGHV) segments in humans, the fundamental differences between these genes and their ability to respond to all possible foreign antigens is still poorly understood. To better understand this, we generated profiles of WGCW hotspots in each of the human IGHV genes and found the expected high frequency in complementarity determining regions (CDRs) that encode the antigen binding sites but also an unexpectedly high frequency of WGCW in certain framework (FW) sub-regions. Principal Components Analysis (PCA) of these overlapping AID hotspot profiles revealed that one major difference between IGHV families is the presence or absence of WGCW in a sub-region of FW3 sometimes referred to as “CDR4.” Further differences between members of each family (e.g., IGHV1) are primarily determined by their WGCW densities in CDR1. We previously suggested that the co-localization of AID overlapping and Polη hotspots was associated with high mutability of certain IGHV sub-regions, such as the CDRs. To evaluate the importance of this feature, we extended the WGCW profiles, combining them with local densities of Polη (WA) hotspots, thus describing the co-localization of both types of hotspots across all IGHV genes. We also verified that co-localization is associated with higher mutability. PCA of the co-localization profiles showed CDR1 and CDR2 as being the main contributors to variance among IGHV genes, consistent with the importance of these sub-regions in antigen binding. Our results suggest that AID overlapping (WGCW) hotspots alone or in conjunction with Polη (WA/TW) hotspots are key features of evolutionary variation between IGHV genes.

Origin of biomolecular games: deception and molecular evolution

Biological macromolecules encode information: some of it to endow the molecule with structural flexibility, some of it to enable molecular actions as a catalyst or a substrate, but a residual part can be used to communicate with other macromolecules. Thus, macromolecules do not need to possess information only to survive in an environment, but also to strategically interact with others by sending signals to a receiving macromolecule that can properly interpret the signal and act suitably. These sender-receiver signalling games are sustained by the information asymmetry that exists among the macromolecules. In both biochemistry and molecular evolution, the important role of information asymmetry remains largely unaddressed. Here, we provide a new unifying perspective on the impact of information symmetry between macromolecules on molecular evolutionary processes, while focusing on molecular deception. Biomolecular games arise from the ability of biological macromolecules to exert precise recognition, and their role as units of selection, meaning that they are subject to competition and cooperation with other macromolecules. Thus, signalling game theory can be used to better understand fundamental features of living systems such as molecular recognition, molecular mimicry, selfish elements and 'junk' DNA. We show how deceptive behaviour at the molecular level indicates a conflict of interest, and so provides evidence of genetic conflict. This model proposes that molecular deception is diagnostic of selfish behaviour, helping to explain the evasive behaviour of transposable elements in 'junk' DNA, for example. Additionally, in this broad review, a range of major evolutionary transitions are shown to be associated with the establishment of signalling conventions, many of which are susceptible to molecular deception. These perspectives allow us to assign rudimentary behaviour to macromolecules, and show how participation in signalling games differentiates biochemistry from abiotic chemistry.

Phenotypic and Ig Repertoire Analyses Indicate a Common Origin of IgD-CD27- Double Negative B Cells in Healthy Individuals and Multiple Sclerosis Patients

IgD^-CD27^- double negative (DN) B cells with proinflammatory characteristics are abnormally elevated in a proportion of multiple sclerosis (MS) patients. In this study, the origin and selection characteristics of DN B cells were studied in MS patients and healthy controls (HC). Expression of developmental markers on peripheral blood DN, IgD^-CD27⁺ class-switched memory (CSM) and IgD⁺CD27^- naive B cells of HC (n = 48) and MS patients (n = 96) was determined by flow cytometry. High-throughput adaptive immune receptor repertoire sequencing was performed on peripheral blood DN and CSM B cells of HC and MS patients (n = 3 each). DN B cells from HC and MS patients showed similar phenotypic and Ig repertoire characteristics. Phenotypic analysis indicated a mature state of DN B cells by low CD5, CD10, and CD38 expression. However, the frequency of CD95⁺ and IgA⁺ cells was lower in DN versus CSM B cells. DN B cells are Ag experienced, as shown by somatic hypermutation of their Ig genes in adaptive immune receptor repertoire sequencing, although they showed a lower mutation load than CSM B cells. Shared clones were found between DN and CSM B cells, although >95% of the clones were unique to each population, and differences in V(D)J usage and CDR3 physicochemical properties were found. Thus, DN B cells arise in HC and MS patients via a common developmental pathway that is probably linked to immune aging. However, DN and CSM B cells develop through unique differentiation pathways, with most DN B cells representing an earlier maturation state.

Inferred Allelic Variants of Immunoglobulin Receptor Genes: A System for Their Evaluation, Documentation, and Naming

Immunoglobulins or antibodies are the main effector molecules of the B-cell lineage and are encoded by hundreds of variable (V), diversity (D), and joining (J) germline genes, which recombine to generate enormous IG diversity. Recently, high-throughput adaptive immune receptor repertoire sequencing (AIRR-seq) of recombined V-(D)-J genes has offered unprecedented insights into the dynamics of IG repertoires in health and disease. Faithful biological interpretation of AIRR-seq studies depends upon the annotation of raw AIRR-seq data, using reference germline gene databases to identify the germline genes within each rearrangement. Existing reference databases are incomplete, as shown by recent AIRR-seq studies that have inferred the existence of many previously unreported polymorphisms. Completing the documentation of genetic variation in germline gene databases is therefore of crucial importance. Lymphocyte receptor genes and alleles are currently assigned by the Immunoglobulins, T cell Receptors and Major Histocompatibility Nomenclature Subcommittee of the International Union of Immunological Societies (IUIS) and managed in IMGT®, the international ImMunoGeneTics information system® (IMGT). In 2017, the IMGT Group reached agreement with a group of AIRR-seq researchers on the principles of a streamlined process for identifying and naming inferred allelic sequences, for their incorporation into IMGT®. These researchers represented the AIRR Community, a network of over 300 researchers whose objective is to promote all aspects of immunoglobulin and T-cell receptor repertoire studies, including the standardization of experimental and computational aspects of AIRR-seq data generation and analysis. The Inferred Allele Review Committee (IARC) was established by the AIRR Community to devise policies, criteria, and procedures to perform this function. Formalized evaluations of novel inferred sequences have now begun and submissions are invited via a new dedicated portal (https://ogrdb.airr-community.org). Here, we summarize recommendations developed by the IARC-focusing, to begin with, on human IGHV genes-with the goal of facilitating the acceptance of inferred allelic variants of germline IGHV genes. We believe that this initiative will improve the quality of AIRR-seq studies by facilitating the description of human IG germline gene variation, and that in time, it will expand to the documentation of TR and IG genes in many vertebrate species.

Beyond Hot Spots: Biases in Antibody Somatic Hypermutation and Implications for Vaccine Design

The evolution of antibodies in an individual during an immune response by somatic hypermutation (SHM) is essential for the ability of the immune system to recognize and remove the diverse spectrum of antigens that may be encountered. These mutations are not produced at random; nucleotide motifs that result in increased or decreased rates of mutation were first reported in 1992. Newer models that estimate the propensity for mutation for every possible 5- or 7-nucleotide motif have emphasized the complexity of SHM targeting and suggested possible new hot spot motifs. Even with these fine-grained approaches, however, non-local context matters, and the mutations observed at a specific nucleotide motif varies between species and even by locus, gene segment, and position along the gene segment within a single species. An alternative method has been provided to further abstract away the molecular mechanisms underpinning SHM, prompted by evidence that certain stereotypical amino acid substitutions are favored at each position of a particular V gene. These "substitution profiles," whether obtained from a single B cell lineage or an entire repertoire, offer a simplified approach to predict which substitutions will be well-tolerated and which will be disfavored, without the need to consider path-dependent effects from neighboring positions. However, this comes at the cost of merging the effects of two distinct biological processes, the generation of mutations, and the selection acting on those mutations. Since selection is contingent on the particular antigens an individual has been exposed to, this suggests that SHM may have evolved to prefer mutations that are most likely to be useful against pathogens that have co-evolved with us. Alternatively, the ability to select favorable mutations may be strongly limited by the biases of SHM targeting. In either scenario, the sequence space explored by SHM is significantly limited and this consequently has profound implications for the rational design of vaccine strategies.

Immunoglobulin gene analysis as a tool for investigating human immune responses

The human immunoglobulin repertoire is a hugely diverse set of sequences that are formed by processes of gene rearrangement, heavy and light chain gene assortment, class switching and somatic hypermutation. Early B cell development produces diverse IgM and IgD B cell receptors on the B cell surface, resulting in a repertoire that can bind many foreign antigens but which has had self-reactive B cells removed. Later antigen-dependent development processes adjust the antigen affinity of the receptor by somatic hypermutation. The effector mechanism of the antibody is also adjusted, by switching the class of the antibody from IgM to one of seven other classes depending on the required function. There are many instances in human biology where positive and negative selection forces can act to shape the immunoglobulin repertoire and therefore repertoire analysis can provide useful information on infection control, vaccination efficacy, autoimmune diseases, and cancer. It can also be used to identify antigen-specific sequences that may be of use in therapeutics. The juxtaposition of lymphocyte development and numerical evaluation of immune repertoires has resulted in the growth of a new sub-speciality in immunology where immunologists and computer scientists/physicists collaborate to assess immune repertoires and develop models of immune action.

In Vitro Evolution of Antibodies Inspired by In Vivo Evolution

In vitro generation of antibodies often requires variable domain sequence evolution to adapt the protein in terms of affinity, specificity, or developability. Such antibodies, including those that are of interest for clinical development, may have their origins in a diversity of immunoglobulin germline genes. Others and we have previously shown that antibodies of different origins tend to evolve along different, preferred trajectories. Apart from substitutions within the complementary determining regions, evolution may also, in a germline gene-origin-defined manner, be focused to residues in the framework regions, and even to residues within the protein core, in many instances at a substantial distance from the antibody’s antigen-binding site. Examples of such germline origin-defined patterns of evolution are described. We propose that germline gene-preferred substitution patterns offer attractive alternatives that should be considered in efforts to evolve antibodies intended for therapeutic use with respect to appropriate affinity, specificity, and product developability. We also hypothesize that such germline gene-origin-defined in vitro evolution hold potential to result in products with limited immunogenicity, as similarly evolved antibodies will be parts of conventional, in vivo-generated antibody responses and thus are likely to have been seen by the immune system in the past.

Accurate immune repertoire sequencing reveals malaria infection driven antibody lineage diversification in young children

Accurately measuring antibody repertoire sequence composition in a small amount of blood is challenging yet important for understanding repertoire responses to infection and vaccination. We develop molecular identifier clustering-based immune repertoire sequencing (MIDCIRS) and use it to study age-related antibody repertoire development and diversification before and during acute malaria in infants (< 12 months old) and toddlers (12-47 months old) with 4-8 ml of blood. Here, we show this accurate and high-coverage repertoire-sequencing method can use as few as 1000 naive B cells. Unexpectedly, we discover high levels of somatic hypermutation in infants as young as 3 months old. Antibody clonal lineage analysis reveals that somatic hypermutation levels are increased in both infants and toddlers upon infection, and memory B cells isolated from individuals who previously experienced malaria continue to induce somatic hypermutations upon malaria rechallenge. These results highlight the potential of antibody repertoire diversification in infants and toddlers.Somatic hypermutation of antibodies can occur in infants but are difficult to track. Here the authors present a new method called MIDCIRS for deep quantitative repertoire sequencing with few cells, and show infants as young as 3 months can expand antibody lineage complexity in response to malaria infection.

Different Somatic Hypermutation Levels among Antibody Subclasses Disclosed by a New Next-Generation Sequencing-Based Antibody Repertoire Analysis

A diverse antibody repertoire is primarily generated by the rearrangement of V, D, and J genes and subsequent somatic hypermutation (SHM). Class-switch recombination (CSR) produces various isotypes and subclasses with different functional properties. Although antibody isotypes and subclasses are considered to be produced by both direct and sequential CSR, it is still not fully understood how SHMs accumulate during the process in which antibody subclasses are generated. Here, we developed a new next-generation sequencing (NGS)-based antibody repertoire analysis capable of identifying all antibody isotype and subclass genes and used it to examine the peripheral blood mononuclear cells of 12 healthy individuals. Using a total of 5,480,040 sequences, we compared percentage frequency of variable (V), junctional (J) sequence, and a combination of V and J, diversity, length, and amino acid compositions of CDR3, SHM, and shared clones in the IgM, IgD, IgG3, IgG1, IgG2, IgG4, IgA1, IgE, and IgA2 genes. The usage and diversity were similar among the immunoglobulin (Ig) subclasses. Clonally related sequences sharing identical V, D, J, and CDR3 amino acid sequences were frequently found within multiple Ig subclasses, especially between IgG1 and IgG2 or IgA1 and IgA2. SHM occurred most frequently in IgG4, while IgG3 genes were the least mutated among all IgG subclasses. The shared clones had almost the same SHM levels among Ig subclasses, while subclass-specific clones had different levels of SHM dependent on the genomic location. Given the sequential CSR, these results suggest that CSR occurs sequentially over multiple subclasses in the order corresponding to the genomic location of IGHCs, but CSR is likely to occur more quickly than SHMs accumulate within Ig genes under physiological conditions. NGS-based antibody repertoire analysis should provide critical information on how various antibodies are generated in the immune system.

Human B-cell memory is shaped by age- and tissue-specific T-independent and GC-dependent events

Switched and IgM memory B cells execute different and noninterchangeable functions. We studied memory B cells in children of different ages, in peripheral blood and spleen and compared them with those of children born asplenic or unable to build germinal centers. We show that, whereas switched memory B cells are mostly generated in the germinal centers at all ages, IgM memory B cells can be distinct in three types with different developmental history. Innate IgM memory B cells, the largest pool in infants, are generated in the spleen by a germinal center-independent mechanism. With age, if the spleen is present and germinal centers are functional, innate IgM memory B cells are remodelled and accumulate somatic mutations. The third type of IgM memory B cell is a by-product of the germinal center reaction. Our data suggest that the B-cell memory developmental program is implemented during the first 5-6 years of life.

Functional Versatility of AGY Serine Codons in Immunoglobulin Variable Region Genes

In systemic autoimmunity, autoantibodies directed against nuclear antigens (Ags) often arise by somatic hypermutation (SHM) that converts AGT and AGC (AGY) Ser codons into Arg codons. This can occur by three different single-base changes. Curiously, AGY Ser codons are far more abundant in complementarity-determining regions (CDRs) of IgV-region genes than expected for random codon use or from species-specific codon frequency data. CDR AGY codons are also more abundant than TCN Ser codons. We show that these trends hold even in cartilaginous fishes. Because AGC is a preferred target for SHM by activation-induced cytidine deaminase, we asked whether the AGY abundance was solely due to a selection pressure to conserve high mutability in CDRs regardless of codon context but found that this was not the case. Instead, AGY triplets were selectively enriched in the Ser codon reading frame. Motivated by reports implicating a functional role for poly/autoreactive specificities in antiviral antibodies, we also analyzed mutations at AGY in antibodies directed against a number of different viruses and found that mutations producing Arg codons in antiviral antibodies were indeed frequent. Unexpectedly, however, we also found that AGY codons mutated often to encode nearly all of the amino acids that are reported to provide the most frequent contacts with Ag. In many cases, mutations producing codons for these alternative amino acids in antiviral antibodies were more frequent than those producing Arg codons. Mutations producing each of these key amino acids required only single-base changes in AGY. AGY is the only codon group in which two-thirds of random mutations generate codons for these key residues. Finally, by directly analyzing X-ray structures of immune complexes from the RCSB protein database, we found that Ag-contact residues generated via SHM occurred more often at AGY than at any other codon group. Thus, preservation of AGY codons in antibody genes appears to have been driven by their exceptional functional versatility, despite potential autoreactive consequences.

Unlike in Children with Allergic Asthma, IgE Transcripts from Preschool Children with Atopic Dermatitis Display Signs of Superantigen-Driven Activation

The IgE repertoire in children with asthma reflects an adaptive B cell response, indicative of Ag-driven selection. However, the same might not apply to atopic dermatitis, which is often the first manifestation of atopy. The objective of our present study was to characterize the IgE repertoire of preschool children with atopic dermatitis with regard to signs of superantigen-like activation, clonal relationship, and indications of Ag selection. Total RNA was isolated from PBMCs of five children with atopic dermatitis. IgE transcripts were amplified, cloned, and sequenced using RT-PCR. We obtained 200 functional IgE sequences, which were compared with 1140 sequences from 11 children with asthma. Whereas variable gene segment of the H Ig chain (VH) gene usage in asthma reflected germline distribution, IgE transcripts from children with atopic dermatitis displayed a dominance of the otherwise scarcely expressed VH2 and VH4 family. Whereas IgE transcripts from children with asthma were highly mutated (7.2%), somatic mutation rate in atopic dermatitis was less than half as high (3.4%). Moreover, the proportion of transcripts that were indicative of Ag selection was reduced to 11% in atopic dermatitis (24% in asthma). In summary, IgE repertoires vary significantly between children with different atopic diseases. Compared with children with asthma, IgE transcripts from preschool children with atopic dermatitis are significantly less mutated, clonally less focused, and less indicative of Ag selection. We consider our data reconcilable with the hypothesis that a superantigen-like activation contributes to the maturation and selection of the IgE repertoire in atopic dermatitis.

Quantifying evolutionary constraints on B-cell affinity maturation

The antibody repertoire of each individual is continuously updated by the evolutionary process of B-cell receptor (BCR) mutation and selection. It has recently become possible to gain detailed information concerning this process through high-throughput sequencing. Here, we develop modern statistical molecular evolution methods for the analysis of B-cell sequence data, and then apply them to a very deep short-read dataset of BCRs. We find that the substitution process is conserved across individuals but varies significantly across gene segments. We investigate selection on BCRs using a novel method that side-steps the difficulties encountered by previous work in differentiating between selection and motif-driven mutation; this is done through stochastic mapping and empirical Bayes estimators that compare the evolution of in-frame and out-of-frame rearrangements. We use this new method to derive a per-residue map of selection, which provides a more nuanced view of the constraints on framework and variable regions.

Genomic variation in the porcine immunoglobulin lambda variable region

Production of a vast antibody repertoire is essential for the protection against pathogens. Variable region germline complexity contributes to repertoire diversity and is a standard feature of mammalian immunoglobulin loci, but functional V region genes are limited in swine. For example, the porcine lambda light chain locus is composed of 23 variable (V) genes and 4 joining (J) genes, but only 10 or 11 V and 2 J genes are functional. Allelic variation in V and J may increase overall diversity within a population, yet lead to repertoire holes in individuals lacking key alleles. Previous studies focused on heavy chain genetic variation, thus light chain allelic diversity is not known. We characterized allelic variation of the porcine immunoglobulin lambda variable (IGLV) region genes. All intact IGLV genes in 81 pigs were amplified, sequenced, and analyzed to determine their allelic variation and functionality. We observed mutational variation across the entire length of the IGLV genes, in both framework and complementarity determining regions (CDRs). Three recombination hotspot motifs were also identified suggesting that non-allelic homologous recombination is an evolutionarily alternative mechanism for generating germline antibody diversity. Functional alleles were greatest in the most highly expressed families, IGLV3 and IGLV8. At the population level, allelic variation appears to help maintain the potential for broad antibody repertoire diversity in spite of reduced gene segment choices and limited germline sequence modification. The trade-off may be a reduction in repertoire diversity within individuals that could result in an increased variation in immunity to infectious disease and response to vaccination.

B cell variable genes have evolved their codon usage to focus the targeted patterns of somatic mutation on the complementarity determining regions

The exceptional ability of B cells to diversify through somatic mutation and improve affinity of the repertoire toward the antigens is the cornerstone of adaptive immunity. Somatic mutation is not evenly distributed and exhibits certain micro-sequence specificities. We show here that the combination of somatic mutation targeting and the codon usage in human B cell receptor (BCR) Variable (V) genes create expected patterns of mutation and post mutation changes that are focused on their complementarity determining regions (CDR). T cell V genes are also skewed in targeting mutations but to a lesser extent and are lacking the codon usage bias observed in BCRs. This suggests that the observed skew in T cell receptors is due to their amino acid usage, which is similar to that of BCRs. The mutation targeting and the codon bias allow B cell CDRs to diversify by specifically accumulating nonconservative changes. We counted the distribution of mutations to CDR in 4 different human datasets. In all four cases we found that the number of actual mutations in the CDR correlated significantly with the V gene mutation biases to the CDR predicted by our models. Finally, it appears that the mutation bias in V genes indeed relates to their long-term survival in actual human repertoires. We observed that resting repertoires of B cells overexpressed V genes that were especially biased toward focused mutation and change in the CDR. This bias in V gene usage was somewhat relaxed at the height of the immune response to a vaccine, presumably because of the need for a wider diversity in a primary response. However, older patients did not retain this flexibility and were biased toward using only highly skewed V genes at all stages of their response.

Identification of tumor-binding scFv derived from clonally related B cells in tumor and lymph node of a patient with breast cancer

The purpose of this study was to determine the clonal relationship between B cells within a breast cancer and the B cells in the tumor-draining lymph node (TDLN). We also determined the binding capacity of antibodies derived from these sources to autologous cancer and autologous noncancer breast tissue. Antibody clonality of B cells derived from tumor and lymph node was determined by analyzing heavy and light chain immunoglobulin sequences. The number of shared clonal groups observed between tumor and lymph node antibodies was significant for both heavy (p = 0.004) and light chain (p = 0.012) populations. Panning with phage-displayed single-chain variable fragment libraries derived from the tumor and lymph node B cells resulted in multiple antibodies that bound autologous tumor. Sequence analysis of enriched antibodies recovered after the third round of panning the tumor and TDLN libraries against autologous tumor lysates had a genetic relationship. These results indicate that B cells infiltrating a patient's breast cancer and B cells present in the tumor-draining lymph node are clonally and functionally related.

Accumulation of self-reactive naïve and memory B cell reveals sequential defects in B cell tolerance checkpoints in Sjögren's syndrome

Sjögren's syndrome (SS) is an autoimmune disease characterised by breach of self-tolerance towards nuclear antigens resulting in high affinity circulating autoantibodies. Although peripheral B cell disturbances have been described in SS, with predominance of naïve and reduction of memory B cells, the stage at which errors in B cell tolerance checkpoints accumulate in SS is unknown. Here we determined the frequency of self- and poly-reactive B cells in the circulating naïve and memory compartment of SS patients. Single CD27-IgD+ naïve, CD27+IgD+ memory unswitched and CD27+IgD- memory switched B cells were sorted by FACS from the peripheral blood of 7 SS patients. To detect the frequency of polyreactive and autoreactive clones, paired Ig VH and VL genes were amplified, cloned and expressed as recombinant monoclonal antibodies (rmAbs) displaying identical specificity of the original B cells. IgVH and VL gene usage and immunoreactivity of SS rmAbs were compared with those obtained from healthy donors (HD). From a total of 353 VH and 293 VL individual sequences, we obtained 114 rmAbs from circulating naïve (n = 66) and memory (n = 48) B cells of SS patients. Analysis of the Ig V gene repertoire did not show significant differences in SS vs. HD B cells. In SS patients, circulating naïve B cells (with germline VH and VL genes) displayed a significant accumulation of clones autoreactive against Hep-2 cells compared to HD (43.1% vs. 25%). Moreover, we demonstrated a progressive increase in the frequency of circulating anti-nuclear naïve (9.3%), memory unswitched (22.2%) and memory switched (27.3%) B cells in SS patients. Overall, these data provide novel evidence supporting the existence of both early and late defects in B cell tolerance checkpoints in patients with SS resulting in the accumulation of autoreactive naïve and memory B cells.

Ig gene analysis reveals altered selective pressures on Ig-producing cells in parotid glands of primary Sjögren's syndrome patients

In this study, we sought to understand the selective pressures shaping the Ig-producing cell repertoire in the parotid glands of primary Sjögren's syndrome (pSS) patients before and after rituximab treatment (RTX). In particular, we evaluated the role of potential N-glycosylation motifs acquired by somatic hypermutation (ac-Nglycs) within Ig H chain V region (IGHV) genes as alternative selective pressures for B cells in pSS. Five pSS patients received RTX. Sequential parotid salivary gland biopsies were taken before RTX, at 12 wk and at 36-52 wk after treatment. Parotid biopsies from four non-pSS patients served as controls. Sequence analysis was carried out on the IgA and IgG RNA transcripts expressing IGHV3 genes in all parotid biopsies. Both IgG and IgA sequences from pSS patients exhibited no evidence for positive Ag-driven selection pressure in their CDRs in contrast to non-pSS controls. The prevalence of IgG sequences with ac-Nglycs was significantly higher in pSS patients than in non-pSS controls. Selection pressures shaping the IgG and IgA repertoire within pSS patients' parotid glands are distinct from those in non-pSS controls, with very little evidence for positive (auto)antigen selection. The higher prevalence of ac-Nglycs on pSS-IgG compared with non-pSS IgG indicates that ac-Nglycs could be an alternative form of selection pressure. We speculate that B cell hyperproliferation within parotid glands of pSS patients may result from Ag-independent interactions such as that between glycosylated B cell receptors and lectins within the microenvironment rather than (auto)antigen-specific stimulation. Our study brings a new perspective into research on pSS.

IgG4 and IgE transcripts in childhood allergic asthma reflect divergent antigen-driven selection

The physiologic function of the "odd" Ab IgG4 remains enigmatic. IgG4 mediates immunotolerance, as, for example, during specific immunotherapy of allergies, but it mediates tissue damage in autoimmune pemphigus vulgaris and "IgG4-related disease." Approximately half of the circulating IgG4 molecules are bispecific owing to their unique ability to exchange half-molecules. Better understanding of the interrelation between IgG4 and IgE repertoires may yield insight into the pathogenesis of allergies and into potential novel therapies that modulate IgG4 responses. We aimed to compare the selective forces that forge the IgG4 and IgE repertoires in allergic asthma. Using an IgG4-specific RT-PCR, we amplified, cloned, and sequenced IgG4 H chain transcripts of PBMCs from 10 children with allergic asthma. We obtained 558 functional IgG4 sequences, of which 286 were unique. Compared with previously published unique IgE transcripts from the same blood samples, the somatic mutation rate was significantly enhanced in IgG4 transcripts (62 versus 83%; p < 0.001), whereas fewer IgG4 sequences displayed statistical evidence of Ag-driven selection (p < 0.001). On average, the hypervariable CDRH3 region was four nucleotides shorter in IgG4 than in IgE transcripts (p < 0.001). IgG4 transcripts in the circulation of children with allergic asthma reflect some characteristics of classical Ag-driven B2 immune responses but display less indication of Ag selection than do IgE transcripts. Although allergen-specific IgG4 can block IgE-mediated allergen presentation and degranulation of mast cells, key factors that influence the Ag-binding properties of the Ab differ between the overall repertoires of circulating IgG4- and IgE-expressing cells.

Defining the recognition elements of Lewis Y-reactive antibodies

Antibody response to carbohydrate antigens is often independent of T cells and the process of affinity/specificity improvement is considered strictly dependent on the germinal centers. Antibodies induced during a T cell-independent type 2 (TI-2) response are less variable and less functionally versatile than those induced with T cell help. The antigen specificity consequences of accumulation of somatic mutations in antibodies during TI-2 responses of Marginal Zone (MZ) B cells is a fact that still needs explanation. Germline genes that define carbohydrate-reactive antibodies are known to sculpt antibody-combining sites containing innate, key side-chain contacts that define the antigen recognition step. However, substitutions associated with MZ B cell derived antibodies might affect the mobility and polyspecificity of the antibody. To examine this hypothesis, we analyzed antibodies reactive with the neolactoseries antigen Lewis Y (LeY) to define the residue subset required for the reactive repertoire for the LeY antigen. Our molecular simulation studies of crystallographically determined and modeled antibody-LeY complexes suggests that the heavy-chain germline gene VH7183.a13.20 and the light-chain Vκ cr1 germline gene are sufficient to account for the recognition of the trisaccharide-H determinant Types 1–4, while the specificity for LeY is driven by the CDR3 backbone conformation of the heavy chain and not the side chain interactions. These results confirm that these monoclonals use germline-encoded amino acids to recognize simple carbohydrate determinants like trisaccharide-H but relies on somatic mutations in the periphery of the combining site to modify affinity for LeY through electrostatic interactions that leads to their optimized binding. These observations bring further attention to the role of mutations in T-cell independent antibodies to distinguish self from non-self carbohydrate antigens.

Postnatal and adult immunoglobulin repertoires of innate-like CD19(+)CD45R(lo) B Cells

The diversity in antibody repertoire relies on different B cell populations working efficiently to fulfil distinct specific functions. We recently described an innate-like CD19(+)CD45R(-/lo) (19(+)45R(lo)) cell population in postnatal unstimulated adult mice, a heterogeneous population containing cells expressing immunoglobulin M (IgM) and others behaving as differentiated mature B lymphocytes (intracytoplasmic IgG1, AID(+), Blimp-1(+)RAG2(-)). In the present study, we characterized the Ig repertoire expressed by splenic 19(+)45R(lo) cells, assuming that they would bear a restricted repertoire biased for germline rearrangements and low mutation rates similar to other innate-like cells. Sequences from 19(+)45R(lo) cells displayed a variety of V, D and J regions, and the analysis of the CDR-H3 region revealed an intermediate overall CDR-H3 length and moderate hydrophobicity. Both IgM and switched sequences of PD15 19(+)45R(lo) cells had shorter CDR-H3 region and fewer non-template N nucleotides than adult sequences, as expected for profiles that correspond to an immature phenotype. Regarding the mutation rate in the VH regions, IgG1 sequences already carried a high rate of replacement mutations at PD15, which increased further in the sequences obtained from adult mice. Moreover, statistical models suggest that a proportion of the switched sequences in adult 19(+)45R(lo) cells had experienced antigen selection, unlike other innate-like B cell compartments.

Human CD19 and CD40L deficiencies impair antibody selection and differentially affect somatic hypermutation

BACKGROUND

Individuals with genetic defects in CD40 ligand (CD40L) or B-cell antigen receptor coreceptor molecules CD19 and CD81 suffer from an antibody deficiency. Still, these patients carry low levels of memory B cells and serum antibodies.

OBJECTIVE

We sought to assess why the remaining memory B cells and antibodies in the blood of these patients do not provide functional immunity.

METHODS

We included CD19-deficient patients (n = 8), CD40L-deficient patients (n = 8), and healthy controls (n = 50) to perform detailed flow cytometry on blood B cells, molecular analysis of IgA and IgG transcripts, as well as functional analysis of B-cell activation.

RESULTS

CD19-deficient and CD40L-deficient patients carried reduced numbers of all memory B-cell subsets except CD27(-)IgA(+) B cells. Their immunoglobulin heavy chain class-switched transcripts contained less somatic mutations and reduced usage of IgM-distal IgG2 and IgA2 subclasses. The selection strength of mutations for antigen binding was significantly lower than in controls, whereas selection to maintain superantigen binding was normal. Furthermore, the patients showed impaired selection against inherently autoreactive properties of their immunoglobulins. Somatic hypermutation analysis revealed decreased activation-induced cytidine deaminase and uracil-DNA glycosylase 2 activity in CD40L deficiency and increased uracil-DNA glycosylase 2 but decreased mismatch repair in CD19 deficiency. B-cell activation studies revealed that this was at least in part due to transcriptional regulation of DNA repair genes.

CONCLUSIONS

This study on CD19 and CD40L deficiencies illustrates that both the B-cell antigen receptor and CD40 signaling pathways are required for the selection of immunoglobulin reactivity. Still, they differentially mediate DNA repair pathways during somatic hypermutation, thereby together shaping the human in vivo antigen-experienced B-cell repertoire.

Integrating B cell lineage information into statistical tests for detecting selection in Ig sequences

Detecting selection in B cell Ig sequences is critical to understanding affinity maturation and can provide insights into Ag-driven selection in normal and pathologic immune responses. The most common sequence-based methods for detecting selection analyze the ratio of replacement and silent mutations using a binomial statistical analysis. However, these approaches have been criticized for low sensitivity. An alternative method is based on the analysis of lineage trees constructed from sets of clonally related Ig sequences. Several tree shape measures have been proposed as indicators of selection that can be statistically compared across cohorts. However, we show that tree shape analysis is confounded by underlying experimental factors that are difficult to control for in practice, including the sequencing depth and number of generations in each clone. Thus, although lineage tree shapes may reflect selection, their analysis alone is an unreliable measure of in vivo selection. To usefully capture the information provided by lineage trees, we propose a new method that applies the binomial statistical framework to mutations identified based on lineage tree structure. This hybrid method is able to detect selection with increased sensitivity in both simulated and experimental data sets. We anticipate that this approach will be especially useful in the analysis of large-scale Ig sequencing data sets generated by high-throughput sequencing technologies.

Combinatorial H3K9acS10ph histone modification in IgH locus S regions targets 14-3-3 adaptors and AID to specify antibody class-switch DNA recombination

Class-switch DNA recombination (CSR) is central to the antibody response, in that it changes the immunoglobulin heavy chain (IgH) constant region, thereby diversifying biological effector functions of antibodies. The activation-induced cytidine deaminase (AID)-centered CSR machinery excises and rejoins DNA between an upstream (donor) and a downstream (acceptor) S region, which precede the respective constant region DNA. AID is stabilized on S regions by 14-3-3 adaptors. These adaptors display a high affinity for 5'-AGCT-3' repeats, which recur in all S regions. However, how 14-3-3, AID, and the CSR machinery target exclusively the donor and acceptor S regions is poorly understood. Here, we show that histone methyltransferases and acetyltransferases are induced by CD40 or Toll-like receptor signaling and catalyze H3K4me3 and H3K9ac/K14ac histone modifications, which are enriched in S regions but do not specify the S region targets of CSR. By contrast, the combinatorial H3K9acS10ph modification specifically marks the S regions set to recombine and directly recruits 14-3-3 adaptors for AID stabilization there. Inhibition of the enzymatic activity of GCN5 and PCAF histone acetyltransferases reduces H3K9acS10ph in S regions, 14-3-3 and AID stabilization, and CSR. Thus, H3K9acS10ph is a histone code that is "written" specifically in S regions and is "read" by 14-3-3 adaptors to target AID for CSR as an important biological outcome.

A Temporal Model of Human IgE and IgG Antibody Function

The diversity of the human antibody repertoire that is generated by V(D)J gene rearrangement is extended by nine constant region genes that give antibodies their complex array of effector functions. The application of high throughput sequencing to the study of V(D)J gene rearrangements has led to significant recent advances in our understanding of the antigen-binding repertoire. In contrast, our understanding of antibody function has changed little, and mystery still surrounds the existence of four distinctive IgG subclasses. Recent observations from murine models and from human studies of VDJ somatic point mutations suggest that the timing of emergence of cells from the germinal center may vary as a consequence of class switching. This should lead to predictable differences in affinity between isotypes. These differences, and varying abilities of the isotypes to fix complement and bind FcRs, could help coordinate the humoral defenses over the time course of a response. We therefore propose a Temporal Model of human IgE and IgG function in which early emergence of IgE sensitizes sentinel mast cells while switching to IgG3 recruits FcγR-mediated functions to the early response. IgG1 then emerges as the major effector of antigen clearance, and subsequently IgG2 competes with IgG1 to produce immune complexes that slow the inflammatory drive. Persisting antigen may finally stimulate high affinity IgG4 that outcompetes other isotypes and can terminate IgG1/FcγR-mediated activation via the inhibitory FcγRIIB. In this way, IgG antibodies of different subclasses, at different concentrations and with sometimes opposing functions deliver cohesive, protective immune function.

Somatic hypermutation analysis in follicular lymphoma provides evidence suggesting bidirectional cell migration between lymph node and bone marrow during disease progression and relapse

In follicular lymphoma, somatic hypermutation of the immunoglobulin heavy chain genes facilitates the identification of different lymphoma cell clones, and the construction of genealogical trees. To investigate the dissemination of lymphoma cells, and the role of bone marrow in disease progression, we simultaneously analyzed the somatic hypermutation patterns of lymph node and bone marrow specimens taken from three patients at onset and relapse of their disease. Immunoglobulin heavy chain genes were amplified by polymerase chain reaction, cloned and sequenced. Mutational pedigrees were constructed in a hierarchical order. When direct transition of one mutation pattern into that of a successor clones was not feasible, hypothetical predecessor clones were created, and a probability measurement calculation was introduced. Eighty-five sequenced clones were generated. The average mutation rates were 13.45% for the lymph node specimens, and 9.78% for the bone marrow ones. Forty-two hypothetical predecessor clones were introduced into inter-compartment pedigrees. The genealogical trees showed that early lymphoma clones with a low mutational load quickly migrate from lymph nodes into the bone marrow. Bi-directional lymphoma cell migration was detectable between the two compartments. In one case of follicular lymphoma, a clone identical to the initial lymph node clone was detected 2 years later in the bone marrow. The newly introduced algorithm allows the evaluation of both time and direction of follicular lymphoma cell migration. We found evidence that follicular lymphoma originates in the lymph node, and infiltrates the bone marrow early in the course of the disease. Moreover, inter-compartment migration between lymph nodes and bone marrow occurs in both directions.

Inconclusive evidence for or against positive antigen selection in the shaping of human immunoglobulin E repertoires: a call for new approaches

BACKGROUND

The mechanisms driving the development of immunoglobulin E (IgE) antibody repertoires are a matter of debate. Alternatives to the classical view on antibody development, involving somatic mutation and antigen-driven selection of high-affinity variants in germinal centers, have been proposed.

METHODS

We have re-analyzed the pattern of mutations in previously isolated and characterized human clonally unrelated IgE-encoding transcripts using the validated focused binomial methodology to find evidence in such genes of antigen-specific selection.

RESULTS

As expected there is a selection against replacement mutations in IgE framework regions. In contrast, in all examined cases but one (assessing IgE repertoires of parasite-infected individuals) there was no evidence in favor of either positive or negative selection in complementarity determining regions. Importantly, however, the validated method also failed to detect selection for replacement mutations in two, non-IgE, hypermutated antibody populations targeting tetanus toxoid and vaccinia virus, respectively.

CONCLUSIONS

Current methodology is unable to define with certainty, using commonly assessed IgE repertoire sizes, whether antigen selection is or is not a major driving force in the establishment of human IgE. New approaches are needed to address this matter.

IgA response in preterm neonates shows little evidence of antigen-driven selection

After birth, contact to environmental Ags induces the production of IgA, which represents a first line of defense for the neonate. We sought to characterize the maturation of the repertoire of IgA H chain transcripts in circulating blood B cells during human ontogeny. We found that IgA H chain transcripts were present in cord blood as early as 27 wk of gestation and that the restrictions of the primary Ab repertoire (IgM) persisted in the IgA repertoire. Thus, B cells harboring more "mature" V(H) regions were not preferred for class switch to IgA. Preterm and term neonates expressed a unique IgA repertoire, which was characterized by short CDR-H3 regions, preference of the J(H) proximal D(H)7-27 gene segment, and very few somatic mutations. During the first postnatal months, these restrictions were slowly released. Preterm birth did not measurably accelerate the maturation of the IgA repertoire. At a postconceptional age of 60 wk, somatic mutation frequency of IgA H chain transcripts reached 25% of the adult values but still showed little evidence of Ag-driven selection. These results indicate that similar to IgG, the IgA repertoire expands in a controlled manner after birth. Thus, the IgA repertoire of the newborn has distinctive characteristics that differ from the adult IgA repertoire. These observations might explain the lower affinity and specificity of neonatal IgA Abs, which could contribute to a higher susceptibility to infections and altered responses to vaccinations, but might also prevent the development of autoimmune and allergic diseases.

Antibodies to a superantigenic glycoprotein 120 epitope as the basis for developing an HIV vaccine

Failure to induce synthesis of neutralizing Abs to the CD4 binding determinant (CD4BD) of gp120, a central objective in HIV vaccine research, has been alternately ascribed to insufficient immunogen binding to Abs in their germline V region configuration expressed as BCRs, insufficient adaptive mutations in Ab V regions, and conformational instability of gp120. We employed peptide analogs of gp120 residues 421-433 within the CD4BD (CD4BD(core)) to identify Abs produced without prior exposure to HIV (constitutive Abs). The CD4BD(core) peptide was recognized by single-chain Fv fragments from noninfected humans with lupus that neutralized genetically diverse strains belonging to various HIV subtypes. Replacing the framework region (FR) of a V(H)4-family single-chain Fv with the corresponding V(H)3-family FRs from single-chain Fv JL427 improved the CD4BD(core) peptide-binding activity, suggesting a CD4BD(core) binding site outside the pocket formed by the CDRs. Replacement mutations in the FR site vicinity suggested the potential for adaptive improvement. A very small subset of serum CD4BD(core)-specific serum IgAs from noninfected humans without autoimmune disease isolated by epitope-specific chromatography neutralized the virus potently. A CD4BD(core)-specific, HIV neutralizing murine IgM with H and L chain V regions (V(H) and V(L) regions) free of immunogen-driven somatic mutations was induced by immunization with a CD4BD(core) peptide analog containing an electrophilic group that binds B cells covalently. The studies indicate broad and potent HIV neutralization by constitutive Abs as an innate, germline-encoded activity directed to the superantigenic CD4BD(core) epitope that is available for amplification for vaccination against HIV.

Molecular characterization of immunoglobulin gene rearrangements in diffuse large B-cell lymphoma: antigen-driven origin and IGHV4-34 as a particular subgroup of the non-GCB subtype

The pathogenesis of diffuse large B-cell lymphoma (DLBCL) remains partially unknown. The analysis of the B-cell receptor of the malignant cells could contribute to a better understanding of the DLBCL biology. We studied the molecular features of the immunoglobulin heavy chain (IGH) rearrangements in 165 patients diagnosed with DLBCL not otherwise specified. Clonal IGH rearrangements were amplified according to the BIOMED-2 protocol and PCR products were sequenced directly. We also analyzed the criteria for stereotyped patterns in all complete IGHV-IGHD-IGHJ (V-D-J) sequences. Complete V-D-J rearrangements were identified in 130 of 165 patients. Most cases (89%) were highly mutated, but 12 sequences were truly unmutated or minimally mutated. Three genes, IGHV4-34, IGHV3-23, and IGHV4-39, accounted for one third of the whole cohort, including an overrepresentation of IGHV4-34 (15.5% overall). Interestingly, all IGHV4-34 rearrangements and all unmutated sequences belonged to the nongerminal center B-cell-like (non-GCB) subtype. Overall, we found three cases following the current criteria for stereotyped heavy chain VH CDR3 sequences, two of them belonging to subsets previously described in CLL. IGHV gene repertoire is remarkably biased, implying an antigen-driven origin in DLBCL. The particular features in the sequence of the immunoglobulins suggest the existence of particular subgroups within the non-GCB subtype.

Immunoglobulin analysis tool: a novel tool for the analysis of human and mouse heavy and light chain transcripts

Sequence analysis of immunoglobulin (Ig) heavy and light chain transcripts can refine categorization of B cell subpopulations and can shed light on the selective forces that act during immune responses or immune dysregulation, such as autoimmunity, allergy, and B cell malignancy. High-throughput sequencing yields Ig transcript collections of unprecedented size. The authoritative web-based IMGT/HighV-QUEST program is capable of analyzing large collections of transcripts and provides annotated output files to describe many key properties of Ig transcripts. However, additional processing of these flat files is required to create figures, or to facilitate analysis of additional features and comparisons between sequence sets. We present an easy-to-use Microsoft^® Excel^® based software, named Immunoglobulin Analysis Tool (IgAT), for the summary, interrogation, and further processing of IMGT/HighV-QUEST output files. IgAT generates descriptive statistics and high-quality figures for collections of murine or human Ig heavy or light chain transcripts ranging from 1 to 150,000 sequences. In addition to traditionally studied properties of Ig transcripts – such as the usage of germline gene segments, or the length and composition of the CDR-3 region – IgAT also uses published algorithms to calculate the probability of antigen selection based on somatic mutational patterns, the average hydrophobicity of the antigen-binding sites, and predictable structural properties of the CDR-H3 loop according to Shirai’s H3-rules. These refined analyses provide in-depth information about the selective forces acting upon Ig repertoires and allow the statistical and graphical comparison of two or more sequence sets. IgAT is easy to use on any computer running Excel^® 2003 or higher. Thus, IgAT is a useful tool to gain insights into the selective forces and functional properties of small to extremely large collections of Ig transcripts, thereby assisting a researcher to mine a data set to its fullest.

Quantifying selection in high-throughput Immunoglobulin sequencing data sets

High-throughput immunoglobulin sequencing promises new insights into the somatic hypermutation and antigen-driven selection processes that underlie B-cell affinity maturation and adaptive immunity. The ability to estimate positive and negative selection from these sequence data has broad applications not only for understanding the immune response to pathogens, but is also critical to determining the role of somatic hypermutation in autoimmunity and B-cell cancers. Here, we develop a statistical framework for Bayesian estimation of Antigen-driven SELectIoN (BASELINe) based on the analysis of somatic mutation patterns. Our approach represents a fundamental advance over previous methods by shifting the problem from one of simply detecting selection to one of quantifying selection. Along with providing a more intuitive means to assess and visualize selection, our approach allows, for the first time, comparative analysis between groups of sequences derived from different germline V(D)J segments. Application of this approach to next-generation sequencing data demonstrates different selection pressures for memory cells of different isotypes. This framework can easily be adapted to analyze other types of DNA mutation patterns resulting from a mutator that displays hot/cold-spots, substitution preference or other intrinsic biases.

Constitutive production of catalytic antibodies to a Staphylococcus aureus virulence factor and effect of infection

Antibodies that recognize microbial B lymphocyte superantigenic epitopes are produced constitutively with no requirement for adaptive immune maturation. We report cleavage of the Staphylococcus aureus virulence factor extracellular fibrinogen-binding protein (Efb) by catalytic antibodies produced with no exposure to the bacterium and reduction of the catalytic antibody activity following infection. IgG catalytic antibodies that specifically hydrolyzed Efb via a nucleophilic catalytic mechanism were found in the blood of healthy humans and aseptic mice free of S. aureus infection. IgG hydrolyzed peptide bonds on the C-terminal side of basic amino acids, including a bond located within the C3b-binding domain of Efb. Efb digested with the IgG lost its ability to bind C3b and inhibit complement-dependent antibody-mediated red blood cell lysis. In addition to catalysis, the IgG expressed saturable Efb binding activity. IgG from S. aureus-infected mice displayed reduced Efb cleaving activity and increased Efb binding activity compared with uninfected controls, suggesting differing effects of the infection on the antibody subsets responsible for the two activities. IgG from children hospitalized for S. aureus infection also displayed reduced Efb cleavage compared with healthy children. These data suggest a potential defense function for constitutively produced catalytic antibodies to a putative superantigenic site of Efb, but an adaptive catalytic response appears to be proscribed.

The evolution of adaptive immunity

The concept of adaptive immunity suggests de novo generation in each individual of extremely large repertoires of diversified receptors and selective expansion of receptors that match the antigen/pathogen. Accordingly, adaptive immune system is also called "anticipatory". It allows each individual to have a unique repertoire of immune receptors corresponding to its life history. The memory of an antigen gets encoded in the clonal composition of the organism's immune cells instead of being encoded in the genome. Consequently, the immune response to repeated encounter with the same antigen becomes stronger, a phenomenon called immunological memory. Elements of adaptive immunity are found at all taxonomical levels, whereas in vertebrates, adaptive mechanisms have become the cornerstone of the immune system. In jaw vertebrates, adaptive immune receptors of T and B lymphoid cells belong to immunoglobulin superfamily and are created by rearrangement of gene segments. In jawless vertebrates lamprey and hagfish, recombination of leucine-rich repeat modules is used to form variable lymphocyte receptors. Striking functional similarity of the cellular and humoral branches of these systems suggests similar driving forces underlying their development.

IgE sequences in individuals living in an area of endemic parasitism show little mutational evidence of antigen selection

Patterns of somatic mutation in IgE genes from allergic individuals have been a focus of study for many years, but IgE sequences have never been reported from parasitized individuals. To study the role of antigen selection in the evolution of the anti-parasite response, we therefore generated 118 IgE sequences from donors living in Papua New Guinea (PNG), an area of endemic parasitism. For comparison, we also generated IgG1, IgG2, IgG3 and IgG4 sequences from these donors, as well as IgG1 sequences from Australian donors. IgE sequences had, on average, 23.0 mutations. PNG IgG sequences had average mutation levels that varied from 17.7 (IgG3) to 27.1 (IgG4). Mean mutation levels correlated significantly with the position of their genes in the constant region gene locus (IgG3 < IgG1 < IgG2 < IgG4). Interestingly, given the heavy, life-long antigen burden experienced by PNG villagers, average mutation levels in IgG sequences were little different to that seen in Australian IgG1 sequences (19.2). Patterns of mutation provide clear evidence of antigen selection in many IgG sequences. The percentage of IgG sequences that showed significant accumulations of replacement mutations in the complementarity determining regions ranged from 22% of IgG3 sequences to 39% of IgG2 sequences. By contrast, only 12% of IgE sequences had such evidence of antigen selection, and this was significantly less than in PNG IgG1, IgG2 and IgG4 subclass sequences (P < 0.01). The anti-parasite IgE response therefore has the reduced evidence of antigen selection that has previously been reported in studies of IgE sequences from allergic individuals.

Detecting selection in immunoglobulin sequences

The ability to detect selection by analyzing mutation patterns in experimentally derived immunoglobulin (Ig) sequences is a critical part of many studies. Such techniques are useful not only for understanding the response to pathogens, but also to determine the role of antigen-driven selection in autoimmunity, B cell cancers and the diversification of pre-immune repertoires in certain species. Despite its importance, quantifying selection in experimentally derived sequences is fraught with difficulties. The necessary parameters for statistical tests (such as the expected frequency of replacement mutations in the absence of selection) are non-trivial to calculate, and results are not easily interpretable when analyzing more than a handful of sequences. We have developed a web server that implements our previously proposed Focused binomial test for detecting selection. Several features are integrated into the web site in order to facilitate analysis, including V(D)J germline segment identification with IMGT alignment, batch submission of sequences and integration of additional test statistics proposed by other groups. We also implement a Z-score-based statistic that increases the power of detecting selection while maintaining specificity, and further allows for the combined analysis of sequences from different germlines. The tool is freely available at http://clip.med.yale.edu/selection.

Evidence of canonical somatic hypermutation in hairy cell leukemia

To compare hairy cell leukemia (HCL) with chronic lymphocytic leukemia (CLL) and normal B cells with respect to their B-cell receptors, somatic hypermutation (SHM) features in HCL were examined in a series of 130 immunoglobulin gene heavy chain rearrangements, including 102 from 100 classic (HCLc) and 28 from 26 variant (HCLv) patients. The frequency of unmutated rearrangements in HCLc was much lower than that in HCLv (17% vs 54%, P < .001) or historically in CLL (17% vs 46%, P < .001), but HCLv and CLL were similar (P = .45). As previously reported for CLL, evidence of canonical SHM was observed in HCLc rearrangements, including: (1) a higher ratio of replacement to silent mutations in the complementarity determining regions than in the framework regions (2.83 vs 1.41, P < .001), (2) higher transition to transversion ratio than would be expected if mutations were random (1.49 vs 0.5, P < .001), and (3) higher than expected concentration of mutations within RGYW hot spots (13.92% vs 3.33%, P < .001). HCLv met these 3 criteria of canonical SHM to a lesser extent. These data suggest that, whereas HCLc cells may recognize antigen-like CLL and normal B cells before malignant transformation, HCLv cells from some patients may originate differently, possibly without undergoing antigen recognition.

Antigen selection in IgE antibodies from individuals with chronic rhinosinusitis with nasal polyps

BACKGROUND

It is generally assumed that all immunoglobulin isotopes develop under antigen selection pressure, leading to dramatic increases in antigen-binding affinity. As activated B cells proliferate, somatic mutations accumulate in the regions of the immunoglobulin gene associated with antigen binding. Emerging evidence from studies investigating mutations in variable region sequences of IgE antibodies suggest that IgE may develop under less selection pressure than other isotypes. Recent studies have focused on IgE mutation patterns in sequences from the blood of allergic individuals. There is, however, little evidence of these patterns in IgE sequences isolated from tissue.

METHODS

Semi-nested reverse-transcription polymerase chain reaction was used to amplify the V region of IgE sequences from nasal tissue of individuals with chronic rhinosinusitis with nasal polyps (CRSwNPs). IgE sequences were analyzed for evidence of antigen selection and compared with previously reported IgE sequences from other inflammatory and allergic disorders and nonallergic individuals.

RESULTS

IgE sequences were successfully amplified from four individuals with CRSwNPs. Of 217 sequences amplified, 38 were unique, 31 of which were from tissue. Identification of the IGHV, IGHD, and IGHJ genes making up each sequence showed overrepresentation of an unusual gene in one individual, but otherwise normal gene usage. Mutation analysis revealed that only 5 of the 31 unique sequences from tissue show clear evidence of antigen selection.

CONCLUSION

With little influence from antigen selection, IgE antibodies are unlikely to be highly specific for antigens. Consequently, these findings have significant implications for the relevance of specific IgE, e.g., Staphylococcus aureus enterotoxin or fungal-specific IgE, in CRSwNP pathogenesis. Whether specific IgE expression is tightly related to pathogenesis or is merely a byproduct of B-cell interactions in local mucosa with colonizing organisms remains unresolved.

Analysis of API2-MALT1 fusion, trisomies, and immunoglobulin VH genes in pulmonary mucosa-associated lymphoid tissue lymphoma

Pulmonary mucosa-associated lymphoid tissue lymphoma is unique in that chronic inflammation is rare and that API2-MALT1 fusion, resulting from t(11;18)(q21;q21), occurs frequently. In this study, we examined 20 cases for API2-MALT1 fusion using the multiplex reverse-transcription polymerase chain reaction and looked for trisomy 3, trisomy 18, and abnormalities of MALT1 and IGH genes using fluorescence in situ hybridization. In addition, we analyzed VH genes by subcloning of the monoclonal polymerase chain reaction products. Of 20 cases studied, we detected gene abnormalities in 16: API2-MALT1 fusion in 9, trisomy 3 in 5, trisomy 18 in 4, MALT1 abnormality in 13, and IGH abnormality in 1. MALT1 gene abnormalities were concordant with API2-MALT1 fusion or trisomy 18. One case showed API2-MALT1 fusion and trisomy 3. On detection of API2-MALT1 fusion and trisomies, we were able to divide our cases into 3 groups, API2-MALT1 positive, trisomy positive, and no detectable gene abnormality, suggesting that tumor development had processed along different genetic pathways. All 20 cases were analyzed for VH genes. Most of the VH genes selected by the lymphomas belonged to the VH3 family, but there was no restriction to any particular VH fragment. Of interest, VH genes were unmutated in 7 cases, suggesting that T-cell-independent extrafollicular B-cell maturation may be important in the development of this lymphoma. In addition, both mutated and unmutated tumor cases were found to carry the API2-MALT1 fusion and trisomy 3. This observation suggests that these gene abnormalities may occur in microenvironments found before or outside of follicular germinal centers.