Problems in using statistical analysis of replacement and silent mutations in antibody genes for determining antigen-driven affinity selection

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

Introduction

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

Methods

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

Results

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

Discussion

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

IGHV mutational status of nodal marginal zone lymphoma by NGS reveals distinct pathogenic pathways with different prognostic implications

The precise B cell of origin and molecular pathogenesis of nodal marginal zone lymphoma (NMZL) remain poorly defined. To date, due to the rarity of NMZL, the vast majority of already-published studies have been conducted on a limited number of samples and the technical approach to analyze the immunoglobulin genes was of amplifying rearranged variable region genes with the classical direct sequencing of the PCR products followed by cloning. Here, we studied the B cell Ig heavy-chain repertoires by next-generation sequencing (NGS) in 30 NMZL cases. Most of the cases were mutated (20/28; 71.5%) with homologies to the respective germ line genes ranging from 85 to 97, 83%, whereas 8/28 (28.5%) were unmutated. In addition, our results show that NMZL cases have a biased usage of specific immunoglobulin heavy-chain variable (IGHV) region genes. Moreover, we documented intraclonal diversity in all (100%) of the mutated cases and ongoing somatic hypermutations (SHM) have been confirmed by hundreds of reads. We analyzed the mutational pattern to detect and quantify antigen selection pressure and we found a positive selection in 4 cases, whereas in the remaining cases there was an unspecific stimulation. Finally, the disease-specific survival and the progression-free survival were significantly different between cases with mutated and unmutated IGHV genes, pointing out mutational status as a possible new biomarker in NMZL.

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.

On the origin of rheumatoid factors: Insights from analyses of variable region sequences

OBJECTIVES

Rheumatoid factors (RFs) are thought to play an important role in rheumatoid arthritis (RA), but are also found in healthy donors (HDs). Previous studies examined variable region sequences of these autoantibodies at a time when knowledge of the human germline repertoire was incomplete. Here we collected and analyzed RF sequence data from the literature to elucidate how RFs develop and whether their characteristics differ between RA patients and HDs.

METHODS

A database was built containing nucleotide sequences of RF heavy and light chain variable domains and characteristics including affinity, isotype and specificity, all collected from published papers. Gene usage and mutation frequencies were analyzed using IMGT/HiV-QUEST. Selection strength was assessed with the BASELINe tool.

RESULTS

Sequences were retrieved for 183 RF clones (87 RA; 67 HDs; 29 other). No biased gene usage was observed for RA and HDs. However, there does appear to be skewed gene usage in RFs from patients with mixed cryoglobulinemia. Mutation frequency varies considerably between RFs, and isotype-switched clones have significantly more mutations. Monospecific RFs carry more mutations than polyspecific RFs; no difference was found for RA- versus HD-derived RFs. Overall, reported affinity is low (median 1 µM), with a non-significant trend toward higher affinity of RA-derived RFs. Mutation frequency and affinity did not appear to be correlated. BASELINe analysis suggests an overall lack of positive selection and less negative selection strength in RA-derived RFs.

CONCLUSIONS

RFs derived from RA patients have similar properties as those derived from HDs. The RF response can be characterized as a moderately matured autoantibody response, with variable levels of somatic hypermutation, but low affinity.

Rare, high-affinity mouse anti-PD-1 antibodies that function in checkpoint blockade, discovered using microfluidics and molecular genomics

Conventionally, mouse hybridomas or well-plate screening are used to identify therapeutic monoclonal antibody candidates. In this study, we present an alternative to hybridoma-based discovery that combines microfluidics, yeast single-chain variable fragment (scFv) display, and deep sequencing to rapidly interrogate and screen mouse antibody repertoires. We used our approach on six wild-type mice to identify 269 molecules that bind to programmed cell death protein 1 (PD-1), which were present at an average of 1 in 2,000 in the pre-sort scFv libraries. Two rounds of fluorescence-activated cell sorting (FACS) produced populations of PD-1-binding scFv with a mean enrichment of 800-fold, whereas most scFv present in the pre-sort mouse repertoires were de-enriched. Therefore, our work suggests that most of the antibodies present in the repertoires of immunized mice are not strong binders to PD-1. We observed clusters of related antibody sequences in each mouse following FACS, suggesting evolution of clonal lineages. In the pre-sort repertoires, these putative clonal lineages varied in both the complementary-determining region (CDR)3K and CDR3H, while the FACS-selected PD-1-binding subsets varied primarily in the CDR3H. PD-1 binders were generally not highly diverged from germline, showing 98% identity on average with germline V-genes. Some CDR3 sequences were discovered in more than one animal, even across different mouse strains, suggesting convergent evolution. We synthesized 17 of the anti-PD-1 binders as full-length monoclonal antibodies. All 17 full-length antibodies bound recombinant PD-1 with K_D < 500 nM (average = 62 nM). Fifteen of the 17 full-length antibodies specifically bound surface-expressed PD-1 in a FACS assay, and nine of the antibodies functioned as checkpoint inhibitors in a cellular assay. We conclude that our method is a viable alternative to hybridomas, with key advantages in comprehensiveness and turnaround time.

B cell repertoire expansion occurs in meningeal ectopic lymphoid tissue

Ectopic lymphoid tissues (ELT) can be found in multiple sclerosis (MS) and other organ-specific inflammatory conditions. Whether ELT in the meninges of central nervous system (CNS) autoimmune disease exhibit local germinal center (GC) activity remains unknown. In an experimental autoimmune encephalomyelitis model of CNS autoimmunity, we found activation-induced cytidine deaminase, a GC-defining enzyme, in meningeal ELT (mELT) densely populated by B and T cells. To determine GC activity in mELT, we excised meningeal lymphoid aggregates using laser capture microscopy and evaluated B cell repertoires in mELT and secondary lymphoid organs by next-generation immune repertoire sequencing. We found immunoglobulin heavy chain variable region sequences that were unique to mELT and had accumulated functionally relevant somatic mutations, together indicating localized antigen-driven affinity maturation. Our results suggest that B cells in mELT actively participate in CNS autoimmunity, which may be relevant to mELT in MS and ELT in other chronic inflammatory conditions.

Chronic signs of memory B cell activation in patients with Behçet's disease are partially restored by anti-tumour necrosis factor treatment

OBJECTIVES

Behçet's disease (BD), an auto-inflammatory vasculitis with oro-genital ulcerations, skin lesions and uveitis, is regarded as T cell mediated. A successful trial with rituximab suggests an additive role for B cells in the pathogenesis. Therefore, we studied B cell abnormalities in BD patients and the effect of TNF-blocking therapy.

METHODS

B cells in blood (n = 36) and tissue (n = 6) of BD patients were analysed with flow cytometry and/or immunohistochemistry and compared with healthy controls (n = 22). BD current activity form (BDCAF) in relation to B cell somatic hypermutations (SHMs) and immunoglobulin class-switching were studied.

RESULTS

Thirty-six patients (17 males) were included, mean age 44 years, average disease duration 10 years and mean BDCAF 2.7. Blood B cell numbers were significantly lower in patients than in controls (P = 0.0061), mostly due to decreased CD27⁺ memory B cells expressing IgM (P = 0.0001), IgG (P = 0.0002) and IgA (P = 0.0038) B cell subsets. CD27⁺ IgA⁺ B cells showed the highest magnitude of decrease in active disease, measured with BDCAF (P = 0.02). CD27⁺ IgM⁺ IgD⁺ B cells were impaired in replication history (P = 0.0133) and selection of SHM, whereas IgA⁺ B cells carried elevated SHM levels (P = 0.04) and lower IgA2 subclass usage (P = 0.0004) than controls. Immunohistochemistry revealed B cells in tissue of active mucosal ulcers. In adalimumab-treated patients, blood B cells were similar to controls.

CONCLUSION

We show significant deviations in the memory B cell compartment, related to disease activity and therapeutic efficacy. Pronounced molecular impairments were seen in the fast-responding IgM⁺-memory and the mucosal IgA⁺-memory B cells. Because of the demonstrated abundance of B cells in affected tissue, we hypothesize relocation of memory B cells to the site of inflammation could account for the deviations found in blood of BD patients. These peripheral B cells are easily accessible as a marker to monitor therapeutic efficacy.

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.

The mutation patterns in B-cell immunoglobulin receptors reflect the influence of selection acting at multiple time-scales

During the several-week course of an immune response, B cells undergo a process of clonal expansion, somatic hypermutation of the immunoglobulin (Ig) genes and affinity-dependent selection. Over a lifetime, each B cell may participate in multiple rounds of affinity maturation as part of different immune responses. These two time-scales for selection are apparent in the structure of B-cell lineage trees, which often contain a 'trunk' consisting of mutations that are shared across all members of a clone, and several branches that form a 'canopy' consisting of mutations that are shared by a subset of clone members. The influence of affinity maturation on the B-cell population can be inferred by analysing the pattern of somatic mutations in the Ig. While global analysis of mutation patterns has shown evidence of strong selection pressures shaping the B-cell population, the effect of different time-scales of selection and diversification has not yet been studied. Analysis of B cells from blood samples of three healthy individuals identifies a range of clone sizes with lineage trees that can contain long trunks and canopies indicating the significant diversity introduced by the affinity maturation process. We here show that observed mutation patterns in the framework regions (FWRs) are determined by an almost purely purifying selection on both short and long time-scales. By contrast, complementarity determining regions (CDRs) are affected by a combination of purifying and antigen-driven positive selection on the short term, which leads to a net positive selection in the long term. In both the FWRs and CDRs, long-term selection is strongly dependent on the heavy chain variable gene family.

A repertoire of high-affinity monoclonal antibodies specific to S. typhi: as potential candidate for improved typhoid diagnostic

Typhoid fever is a significant global health problem with highest burden on the developing world. The severity of typhoid is often underestimated, and currently available serological diagnostic assays are inadequate due to lack in requisite sensitivity and specificity. This underlines an absolute need to develop a reliable and accurate diagnostics that would benefit long-term disease control and treatment and to understand the real disease burden. Here, we have utilized flagellin protein of S. typhi that is surface accessible, abundantly expressed, and highly immunogenic, for developing immunodiagnostic tests. Flagellin monomers are composed of conserved amino-terminal and carboxy-terminal, and serovar-specific middle region. We have generated a panel of murine monoclonal antibodies (mAbs) against the middle region of flagellin, purified from large culture of S. typhi to ensure its native conformation. These mAbs showed unique specificity and very high affinity toward S. typhi flagellin without showing any cross-reactivity with other serovars. Genetic analysis of mAbs also revealed high frequency of somatic mutation due to antigenic selection process across variable region to achieve high binding affinity. These antibodies also displayed stable binding in stringent reaction conditions for antigen-antibody interactions, like DMSO, urea, KSCN, guanidinium HCl, and extremes of pH. One of the mAbs potentially reversed the TLR5-mediated immune response, in vitro by inhibiting TLR5-flagellin interaction. In our study, binding of these mAbs to flagellin, with high affinity, present on bacterial surface, as well as in soluble form, validates their potential use in developing improved diagnostics with significantly higher sensitivity and specificity.

The splenic autoimmune response to ADAMTS13 in thrombotic thrombocytopenic purpura contains recurrent antigen-binding CDR3 motifs

The spleen harbors ADAMTS13-specific memory B cells following acute acquired TTP. The splenic anti-ADAMTS13 antibody repertoire is characterized by a set of unique and novel CDR3 motifs, 4 shared by 2 patients.

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.

Models of somatic hypermutation targeting and substitution based on synonymous mutations from high-throughput immunoglobulin sequencing data

Analyses of somatic hypermutation (SHM) patterns in B cell immunoglobulin (Ig) sequences contribute to our basic understanding of adaptive immunity, and have broad applications not only for understanding the immune response to pathogens, but also to determining the role of SHM in autoimmunity and B cell cancers. Although stochastic, SHM displays intrinsic biases that can confound statistical analysis, especially when combined with the particular codon usage and base composition in Ig sequences. Analysis of B cell clonal expansion, diversification, and selection processes thus critically depends on an accurate background model for SHM micro-sequence targeting (i.e., hot/cold-spots) and nucleotide substitution. Existing models are based on small numbers of sequences/mutations, in part because they depend on data from non-coding regions or non-functional sequences to remove the confounding influences of selection. Here, we combine high-throughput Ig sequencing with new computational analysis methods to produce improved models of SHM targeting and substitution that are based only on synonymous mutations, and are thus independent of selection. The resulting “S5F” models are based on 806,860 Synonymous mutations in 5-mer motifs from 1,145,182 Functional sequences and account for dependencies on the adjacent four nucleotides (two bases upstream and downstream of the mutation). The estimated profiles can explain almost half of the variance in observed mutation patterns, and clearly show that both mutation targeting and substitution are significantly influenced by neighboring bases. While mutability and substitution profiles were highly conserved across individuals, the variability across motifs was found to be much larger than previously estimated. The model and method source code are made available at http://clip.med.yale.edu/SHM

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.

Constant domain-regulated antibody catalysis

Some antibodies contain variable (V) domain catalytic sites. We report the superior amide and peptide bond-hydrolyzing activity of the same heavy and light chain V domains expressed in the IgM constant domain scaffold compared with the IgG scaffold. The superior catalytic activity of recombinant IgM was evident using two substrates, a small model peptide that is hydrolyzed without involvement of high affinity epitope binding, and HIV gp120, which is recognized specifically by noncovalent means prior to the hydrolytic reaction. The catalytic activity was inhibited by an electrophilic phosphonate diester, consistent with a nucleophilic catalytic mechanism. All 13 monoclonal IgMs tested displayed robust hydrolytic activities varying over a 91-fold range, consistent with expression of the catalytic functions at distinct levels by different V domains. The catalytic activity of polyclonal IgM was superior to polyclonal IgG from the same sera, indicating that on average IgMs express the catalytic function at levels greater than IgGs. The findings indicate a favorable effect of the remote IgM constant domain scaffold on the integrity of the V-domain catalytic site and provide a structural basis for conceiving antibody catalysis as a first line immune function expressed at high levels prior to development of mature IgG class antibodies.

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.

Somatic hypermutation of TCR γ V genes in the sandbar shark

In a recent publication we demonstrated that somatic hypermutation occurs in the V region of the TCR γ gene of the sandbar shark (Carcharhinus plumbeus). We hypothesize that similar mechanisms are used to generate somatic mutations in both immunoglobulin and TCR γ genes of the sharks. Two distinct patterns of mutation occur, single nucleotide mutations (point mutations) and mutations comprising 2-5 consecutive bases (tandem mutations). Our data indicates that point mutations occur by a mechanism similar to that of somatic hypermutation in immunoglobulin genes of mammals, whereas tandem mutations may be generated by an error-prone DNA polymerase with terminal deoxynucleotidyl transferase (TdT)-like activity. Shark hotspot motifs identical to those of higher vertebrates were identified. We confirm that, as in immunoglobulin of sharks and higher vertebrates, highly significant targeting of AID activity to the classical DGYW/WRCH motif occurs in somatic hypermutation of sandbar shark TCR γ V genes. Our analysis suggests that the purpose of somatic mutations in shark TCR γ V-regions is to generate a more diverse repertoire in γ/δ receptors, rather than receptors with higher affinity.

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.

Human memory B cells originate from three distinct germinal center-dependent and -independent maturation pathways

Multiple distinct memory B-cell subsets have been identified in humans, but it remains unclear how their phenotypic diversity corresponds to the type of responses from which they originate. Especially, the contribution of germinal center-independent responses in humans remains controversial. We defined 6 memory B-cell subsets based on their antigen-experienced phenotype and differential expression of CD27 and IgH isotypes. Molecular characterization of their replication history, Ig somatic hypermutation, and class-switch profiles demonstrated their origin from 3 different pathways. CD27⁻IgG⁺ and CD27⁺IgM⁺ B cells are derived from primary germinal center reactions, and CD27⁺IgA⁺ and CD27⁺IgG⁺ B cells are from consecutive germinal center responses (pathway 1). In contrast, natural effector and CD27⁻IgA⁺ memory B cells have limited proliferation and are also present in CD40L-deficient patients, reflecting a germinal center-independent origin. Natural effector cells at least in part originate from systemic responses in the splenic marginal zone (pathway 2). CD27⁻IgA⁺ cells share low replication history and dominant Igλ and IgA2 use with gut lamina propria IgA+ B cells, suggesting their common origin from local germinal center-independent responses (pathway 3). Our findings shed light on human germinal center-dependent and -independent B-cell memory formation and provide new opportunities to study these processes in immunologic diseases.

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.

In situ B cell-mediated immune responses and tubulointerstitial inflammation in human lupus nephritis

The most prevalent severe manifestation of systemic lupus erythematosus is nephritis, which is characterized by immune complex deposition, inflammation, and scarring in glomeruli and the tubulointerstitium. Numerous studies indicated that glomerulonephritis results from a systemic break in B cell tolerance, resulting in the local deposition of immune complexes containing Abs reactive with ubiquitous self-Ags. However, the pathogenesis of systemic lupus erythematosus tubulointerstitial disease is not known. In this article, we demonstrate that in more than half of a cohort of 68 lupus nephritis biopsies, the tubulointerstitial infiltrate was organized into well-circumscribed T:B cell aggregates or germinal centers (GCs) containing follicular dendritic cells. Sampling of the in situ-expressed Ig repertoire revealed that both histological patterns were associated with intrarenal B cell clonal expansion and ongoing somatic hypermutation. However, in the GC histology, the proliferating cells were CD138(-)CD20(+) centroblasts, whereas they were CD138(+)CD20(low/-) plasmablasts in T:B aggregates. The presence of GCs or T:B aggregates was strongly associated with tubular basement membrane immune complexes. These data implicate tertiary lymphoid neogenesis in the pathogenesis of lupus tubulointerstitial inflammation.

Ig gene diversification and selection in follicular lymphoma, diffuse large B cell lymphoma and primary central nervous system lymphoma revealed by lineage tree and mutation analyses

Follicular lymphoma (FL), diffuse large B cell lymphoma (DLBCL) and primary central nervous system lymphoma are B cell malignancies. FL and DLBCL have a germinal center origin. We have applied mutational analyses and a novel algorithm for quantifying shape properties of mutational lineage trees to investigate the nature of the diversification, somatic hypermutation and selection processes that affect B cell clones in these malignancies and reveal whether they differ from normal responses. Lineage tree analysis demonstrated higher diversification and mutations per cell in the lymphoma clones. This was caused solely by the longer diversification times of the malignant clones, as their recent diversification processes were similar to those of normal responses, implying similar mutation frequencies. Since previous analyses of antigen-driven selection were shown to yield false positives, we performed a corrected analysis of replacement and silent mutation patterns, which revealed selection against replacement mutations in the framework regions, responsible for the structural integrity of the B cell receptor, but not for positive selection for replacements in the complementary determining regions. Most replacements, however, were neutral or conservative, suggesting that if at all selection operates in these malignancies it is for structural B cell receptor integrity but not for antigen binding.

E2A-positive gastric MALT lymphoma has weaker plasmacytoid infiltrates and stronger expression of the memory B-cell-associated miR-223: possible correlation with stage and treatment response

Extranodal marginal-zone lymphoma of mucosa-associated lymphoid tissue of the stomach (gastric MALT lymphoma) is derived from memory B cells of the marginal zone. Normal memory B cells do not express markers of germinal-center B cells, such as E2A (immunoglobulin enhancer-binding factor E12/E47), B-cell chronic lymphocytic leukemia/lymphoma 6 (BCL6), or activation-induced cytidine deaminase (AID). E2A is a transcription factor that induces somatic hypermutations and blocks plasma cell differentiation. In 50 stage-I(E)/II(E1) gastric MALT lymphomas, we confirmed that all cases were BCL6(-)/AID(-), but a subset (50%, 25/50) was E2A(+). As E2A(-) and E2A(+) gastric MALT lymphomas had similar numbers of somatic hypermutations without intraclonal variations, which implied an origin from memory B cells, the expression of E2A was best regarded as a marker of aberrant follicular differentiation. Although the status of somatic hypermutation was not affected by E2A, E2A(+) gastric MALT lymphoma showed less plasmacytoid infiltrates and higher expressions of miRNA-223, a microRNA associated with memory B cells. Clinically, E2A(+) gastric MALT lymphomas were more likely to spread to perigastric lymph nodes and were less responsive to Helicobacter eradication therapy than were E2A(-) gastric MALT lymphomas. Taken together, aberrant E2A expression is a diagnostic feature of a subtype of gastric MALT lymphoma with weaker plasmacytoid infiltrates and stronger miR-223 expression. A prospective study would be necessary to verify the association between E2A expression and a poor response to Helicobacter eradication therapy.

Somatic hypermutation and antigen-driven selection of B cells are altered in autoimmune diseases

B cells have been found to play a critical role in the pathogenesis of several autoimmune (AI) diseases. A common feature amongst many AI diseases is the formation of ectopic germinal centers (GC) within the afflicted tissue or organ, in which activated B cells expand and undergo somatic hypermutation (SHM) and antigen-driven selection on their immunoglobulin variable region (IgV) genes. However, it is not yet clear whether these processes occurring in ectopic GCs are identical to those in normal GCs. The analysis of IgV mutations has aided in revealing many aspects concerning B cell expansion, mutation and selection in GC reactions. We have applied several mutation analysis methods, based on lineage tree construction, to a large set of data, containing IgV productive and non-productive heavy and light chain sequences from several different tissues, to examine three of the most profoundly studied AI diseases - Rheumatoid Arthritis (RA), Multiple Sclerosis (MS) and Sjögren's Syndrome (SS). We have found that RA and MS sequences exhibited normal mutation spectra and targeting motifs, but a stricter selection compared to normal controls, which was more apparent in RA. SS sequence analysis results deviated from normal controls in both mutation spectra and indications of selection, also showing differences between light and heavy chain IgV and between different tissues. The differences revealed between AI diseases and normal control mutation patterns may result from the different microenvironmental influences to which ectopic GCs are exposed, relative to those in normal secondary lymphoid tissues.

Toward effective HIV vaccination: induction of binary epitope reactive antibodies with broad HIV neutralizing activity

We describe murine monoclonal antibodies (mAbs) raised by immunization with an electrophilic gp120 analog (E-gp120) expressing the rare ability to neutralize genetically heterologous human immunodeficiency virus (HIV) strains. Unlike gp120, E-gp120 formed covalent oligomers. The reactivity of gp120 and E-gp120 with mAbs to reference neutralizing epitopes was markedly different, indicating their divergent structures. Epitope mapping with synthetic peptides and electrophilic peptide analogs indicated binary recognition of two distinct gp120 regions by anti-E-gp120 mAbs, the 421-433 and 288-306 peptide regions. Univalent Fab and single chain Fv fragments expressed the ability to recognize both peptides. X-ray crystallography of an anti-E-gp120 Fab fragment revealed two neighboring cavities, the typical antigen-binding cavity formed by the complementarity determining regions (CDRs) and another cavity dominated by antibody heavy chain variable (V(H)) domain framework (FR) residues. Substitution of the FR cavity V(H) Lys-19 residue by an Ala residue resulted in attenuated binding of the 421-433 region peptide probe. The CDRs and V(H) FR replacement/silent mutation ratios exceeded the ratio for a random mutation process, suggesting adaptive development of both putative binding sites. All mAbs studied were derived from V(H)1 family genes, suggesting biased recruitment of the V gene germ line repertoire by E-gp120. The conserved 421-433 region of gp120 is essential for HIV binding to host CD4 receptors. This region is recognized weakly by the FR of antibodies produced without exposure to HIV, but it usually fails to induce adaptive synthesis of neutralizing antibodies. We present models accounting for improved CD4-binding site recognition and broad HIV neutralizing activity of the mAbs, long sought goals in HIV vaccine development.

Antigen-driven selection in germinal centers as reflected by the shape characteristics of immunoglobulin gene lineage trees: a large-scale simulation study

During the immune response, the generation of memory B lymphocytes in germinal centers involves affinity maturation of the cells' antigen receptors, based on somatic hypermutation of receptor genes and antigen-driven selection of the resulting mutants. Affinity maturation is vital for immune protection, and is the basis of humoral immune learning and memory. Lineage trees of somatically hypermutated immunoglobulin genes often serve to qualitatively illustrate claims concerning the dynamics of affinity maturation in germinal centers. Here, we derive the quantitative relationships between parameters characterizing affinity maturation dynamics (proliferation, differentiation and mutation rates, initial affinity of the Ig to the antigen, and selection thresholds) and the mathematical properties of lineage trees, using a computer simulation which combines mathematical models for all mature B cell populations, stochastic models of hypermutation and selection, lineage tree generation and measurement of graphical tree characteristics. We identified seven key lineage tree properties, and found correlations of these with initial clone affinity and with the selection threshold. These two parameters were found to be the main factors affecting lineage tree shapes in both primary and secondary response trees. The results also confirm that recycling from centrocytes back to centroblasts is highly likely.

Improved methods for detecting selection by mutation analysis of Ig V region sequences

Statistical methods based on the relative frequency of replacement mutations in B lymphocyte Ig V region sequences have been widely used to detect the forces of selection that shape the B cell repertoire. However, current methods produce an unexpectedly high frequency of false positives and are sensitive to intrinsic biases of somatic hypermutation that can give the appearance of selection. The new statistical test proposed here provides a better trade-off between sensitivity and specificity compared with previous approaches. The low specificity of existing methods was shown in silico to result from an interaction between the effects of positive and negative selection. False detection of positive selection was confirmed in vivo through a re-analysis of published sequence data from diffuse large B cell lymphomas, highlighting the need for re-analysis of some existing studies. The sensitivity of the proposed method to detect selection was validated using new Ig transgenic mouse models in which positive selection was expected to be a significant force, as well as with a simulation-based approach. Previous concerns that intrinsic biases of somatic hypermutation could give the appearance of selection were addressed by extending the current mutation models to more fully account for the impact of microsequence on relative mutability and to include transition bias. High specificity was confirmed using a large set of non-productively rearranged Ig sequences. These results show that selection can be detected in vivo with high specificity using the new method proposed here, allowing greater insight into the existence and direction of antigen-driven selection.

Influence of EBV on the Peripheral Blood Memory B Cell Compartment

Peripheral blood memory B cells latently infected with EBV bear somatic mutations and are typically isotype switched consistent with being classical Ag-selected memory B cells. In this work, we performed a comparative analysis of the expressed Ig genes between large sets of EBV-infected and uninfected peripheral blood B cells, isolated from the same infectious mononucleosis patients, to determine whether differences exist that could reveal the influence of EBV on the production and maintenance of these cells. We observed that EBV(+) cells on average accumulated more somatic hypermutations than EBV(-) cells. In addition, they had more replacement mutations and a higher replacement-silent ratio of mutations in their CDRs. We also found that EBV occupies a skewed niche within the memory compartment, due to its exclusion from the CD27(+)IgD(+)IgM(+) subset, but this skewing does not affect the overall structure of the compartment. These results indicate that EBV impacts the mutation and selection process of infected cells but that once they enter memory they cannot be distinguished from uninfected cells by host homeostasis mechanisms.

Lineage tree analysis of immunoglobulin variable-region gene mutations in autoimmune diseases: chronic activation, normal selection

Autoimmune diseases show high diversity in the affected organs, clinical manifestations and disease dynamics. Yet they all share common features, such as the ectopic germinal centers found in many affected tissues. Lineage trees depict the diversification, via somatic hypermutation (SHM), of immunoglobulin variable-region (IGV) genes. We previously developed an algorithm for quantifying the graphical properties of IGV gene lineage trees, allowing evaluation of the dynamical interplay between SHM and antigen-driven selection in different lymphoid tissues, species, and disease situations. Here, we apply this method to ectopic GC B cell clones from patients with Myasthenia Gravis, Rheumatoid Arthritis, and Sjögren's Syndrome, using data scaling to minimize the effects of the large variability due to methodological differences between groups. Autoimmune trees were found to be significantly larger relative to normal controls. In contrast, comparison of the measurements for tree branching indicated that similar selection pressure operates on autoimmune and normal control clones.

A new memory CD27-IgG+ B cell population in peripheral blood expressing VH genes with low frequency of somatic mutation

In humans, up to 40% of peripheral B cells express CD27 and have hypermutated variable regions in their Ig genes. The CD27+ B cells are considered to be derived from germinal center following specific antigenic stimulation. Actually, somatic hypermutation in Ig genes and CD27 expression are hallmarks of memory B cells. However, the blood IgM+ IgD+ CD27+ B cells were recently associated to splenic marginal zone B cells and proposed to be a subset distinct from germinal center-derived memory B cells showing premutated Igs. The results presented herein further weaken this bona fide association because B cells expressing surface IgG, but not CD27, were found in human blood. Representing 1-4% of all peripheral B cells and approximately 25% of the IgG+ blood B cells, this population expressed mutated IgG genes showing antigenic selection characteristics but with lower mutation frequencies than that of CD27+ IgG+ B cells. However, their morphology and phenotype were similar to that of CD27+ IgG+ cells. Interestingly, the proportion of IgG2 over IgG3 transcripts was opposite in CD27- IgG+ and CD27+ IgG+ cells, suggesting distinct functions or origins. Overall, these findings extend the memory B cell reservoir beyond the CD27+ compartment and could provide further insights into B cell disorders of unknown etiology.

Correlation of molecular characteristics, isotype, and in vitro functional activity of human antipneumococcal monoclonal antibodies

Structure-function correlations of pneumococcal antibodies are important in predicting how changes in the pneumococcus (Pnc)-specific B-cell repertoire will influence humoral immunity against invasive Pnc disease. Using a unique panel of human hybridomas derived from memory B cells after pneumococcal conjugate vaccination, we analyzed the structure-function relationship of nine monoclonal antibodies (MAbs) reactive to Pnc polysaccharides. The avidities of the antibodies correlated with the avidity of donor immune serum (R, 0.7; P < 0.025), and this relationship was particularly strong for immunoglobulin A clones (R, 1; P < 0.0005), suggesting that the MAbs may represent important clones contributing to serological memory. Common heavy-light chain combinations and amino acid replacement mutations were seen for clones with the same serospecificity from different individuals. The two highest-avidity MAbs used Vh3-48, and two MAbs with the same serospecificity, using the same V gene pairings (Vh3-7 and Vk2A17), had similar avidities, suggesting that canonical V gene use makes an important contribution to avidity. Although all clones had mutation levels consistent with their being derived from memory B cells, low levels of replacement mutation were associated with high avidities. This relationship was strongest for Vh genes (R, 0.8; P < 0.01). Opsonophagocytosis was demonstrated for all clones, and there was a trend toward clones using canonical genes with low levels of mutation having high opsonophagocytic activities (R, 0.5). These data suggest that the use of canonical genes in the Pnc antibody response is associated with highly functional antibodies and that most somatic mutations seen in these genes are not antigen selected.